Genome-Wide Study Updates in the International Genetics and Translational Research in Transplantation Network (iGeneTRAiN)
- 1Division of Transplantation Department of Surgery, University of Pennsylvania, Philadelphia, PA, United States
- 2Department of Cardiology, University Medical Center Utrecht, University of Utrecht, Utrecht, Netherlands
- 3Department of Medicine, Division of Nephrology, Vagelos College of Physicians & Surgeons, Columbia University, New York, NY, United States
- 4Immunogenetics and Biology of Transplantation, Città della Salute e della Scienza, University Hospital of Turin, Turin, Italy
- 5Medical Genetics, Department of Medical Sciences, University Turin, Turin, Italy
- 6Department of Hospital Pharmacy, University Medical Center Rotterdam, Rotterdam, Netherlands
- 7Department of Nephrology, Khung Hee University, Seoul, South Korea
- 8Division of Nephrology, St. Joseph's Healthcare Hamilton, McMaster University, Hamilton, ON, Canada
The prevalence of end-stage renal disease (ESRD) and the number of kidney transplants performed continues to rise every year, straining the procurement of deceased and living kidney allografts and health systems. Genome-wide genotyping and sequencing of diseased populations have uncovered genetic contributors in substantial proportions of ESRD patients. A number of these discoveries are beginning to be utilized in risk stratification and clinical management of patients. Specifically, genetics can provide insight into the primary cause of chronic kidney disease (CKD), the risk of progression to ESRD, and post-transplant outcomes, including various forms of allograft rejection. The International Genetics & Translational Research in Transplantation Network (iGeneTRAiN), is a multi-site consortium that encompasses >45 genetic studies with genome-wide genotyping from over 51,000 transplant samples, including genome-wide data from >30 kidney transplant cohorts (n = 28,015). iGeneTRAiN is statistically powered to capture both rare and common genetic contributions to ESRD and post-transplant outcomes. The primary cause of ESRD is often difficult to ascertain, especially where formal biopsy diagnosis is not performed, and is unavailable in ∼2% to >20% of kidney transplant recipients in iGeneTRAiN studies. We overview our current copy number variant (CNV) screening approaches from genome-wide genotyping datasets in iGeneTRAiN, in attempts to discover and validate genetic contributors to CKD and ESRD. Greater aggregation and analyses of well phenotyped patients with genome-wide datasets will undoubtedly yield insights into the underlying pathophysiological mechanisms of CKD, leading the way to improved diagnostic precision in nephrology.
The global prevalence of end-stage renal disease (ESRD) continues to climb. In 2016, 19,301 kidney transplants were performed in the United States, and approximately five times as many were performed worldwide1,2. Due to improvements in surgical techniques, immunosuppression protocols, and clinical management of post-transplant complications, the five-year graft survival rates for kidneys obtained from deceased and living donors reached highs of 75.3% and 85.3%, respectively (Cohen et al.,2006;Serur et al., 2011; Vignolini et al., 2019). However, the prevalence of ESRD cases in the US has continued to rise by ∼20,000 cases per year over the past three decades, creating an increased need for kidney allografts1. This increase is believed to be due primarily to worsening diets and other modifiable factors associated with Western lifestyle but also to an increase in the longevity of pre-transplant ESRD cases.
It is well established that genetic factors contribute to the development and progression of specific types of chronic kidney disease (CKD), yet many previous studies have been limited in scope due to small sample sizes and genotyping strategies (Azarpira et al., 2014; Misra et al., 2014; Phelan et al 2014; Parsa et al., 2017; Stapleton et al., 2018). Studies of families with severe phenotypes of diseases, such as Alport’s Syndrome and Fabry Disease, have significantly contributed to the understanding of the genetic characteristics of these conditions (Gillion et al., 2018; Kashiwagi et al., 2018; McCloskey et al., 2018). However, milder forms of these diseases and their role in the development of ESRD have yet to be explored in great depth.
Genome-Wide Genotyping Arrays
Array based genome-wide genotyping from diverse patient populations facilitates very precise ancestry determination using methods such as principal component analysis (Cai et al., 2013; Li et al., 2015). Genome-wide association studies (GWAS) among patients with CKD have detected both rare and common genetic variants significantly associated with estimated glomerular filtration rate (eGFR) decline and microalbuminuria, some of the strongest predictors of CKD outcomes, despite >80% of GWAS participants having eGFRs in the normal range (Boger et al., 2011a;Boger et al., 2011b;Reznichenko et al., 2012; Gorski et al., 2015; Parsa et al., 2017; Limou et al., 2018).
The findings of genome-wide studies may also provide new therapeutic targets to slow the progression of CKD to ESRD, which may delay or impact the need for transplantation in some patient populations (Wuttke & Kottgen, 2016; Kalatharan et al., 2018). For example, nephropathic cystinosis, a rare autosomal recessive disease, is caused by a 57-kb deletion in the CTNS gene in ∼75% of patients of European ancestry and progresses to ESRD if left untreated (Brodin-Sartorius et al., 2012). However, treatment with oral cysteamine by five years of age has been found to significantly decrease the prevalence and delay the onset of ESRD (Brodin-Sartorius et al., 2012). Additionally, at least 38 genes have been associated with the development of genetic focal segmental glomerulosclerosis (FSGS), some of which have been shown to be responsive to glucocorticoid treatment (Rosenberg & Kopp, 2017). GWAS findings can also provide insight into the biology of ESRD, helping to remove diagnostic heterogeneity. The two APOL1 risk alleles (G1 and G2) found in high frequency in sub-Saharan African populations and strongly associated with FSGS and HIV nephropathy were found to activate protein kinase R, thus inducing glomerular injury and proteinuria (Kopp et al., 2011; Limou et al., 2014; Okamoto et al., 2018). Overall, results from genome-wide screening can enable physicians to provide accurate genetic diagnoses for the primary cause of ESRD, enabling timely and effective therapeutic managemenvwt and aiding in the evaluation of family members as living donors (Snoek et al., 2018).
Whole-Exome and Whole-Genome Sequencing
In the last decade, whole-exome sequencing (WES) and whole-genome sequencing (WGS) approaches have been used very successfully to discover and diagnose genetic disorders in a clinical context (Mallawaarachchi et al., 2016; Lata et al., 2018; Warejko et al., 2018;Groopman et al., 2019). WES typically yields sufficient depth of sequencing coverage across ∼95% of nucleotides in coding regions captured and has been used to diagnose rare high penetrant, Mendelian disorders, discover common variants, and identify causal mutations in cancer (Huang et al., 2018; Zhang et al., 2018). WES has recently been implemented as a first-line diagnostic tool in clinical medicine. In a study on fetuses with congenital anomalies of the kidney and urinary tract (CAKUT), pathogenic variants were discovered in 13% of cases (Lei et al., 2017). WES has also been applied to adult-onset CKD and ESRD, in which ∼10% of cases are caused by Mendelian mutations (Wuhl et al., 2014; Lata et al., 2018; Groopman et al., 2019). In a cohort of >3,000 patients with advanced CKD and ESRD ascertained for a clinical trial, WES identified diagnostic variants in 9.3% of patients encompassing 66 monogenic disorders (Groopman et al., 2019). Of the 343 detected variants, 141 (41%) had not been previously reported as pathogenic. Additionally, diagnostic variants were identified in 17.1% of individuals with nephropathy of unknown origin, altering medical management by initiating multidisciplinary care, prompting referral to clinical trials, and guiding donor selection for transplantation (Groopman et al., 2019). However, it should be noted that many CKD studies using WES have struggled to obtain adequate control populations. iGeneTRAiN has a large pool of healthy donors (in kidney and in other organs), which represents a strong advantage for our study designs.
WGS is the most comprehensive approach for the detection of inherited variants due to more complete genome-wide coverage, although there are additional challenges compared to WES. WGS can capture single nucleotide genetic variants, small Insertions and Deletions (Indels), and Copy-Number Variants (Cnvs) throughout the human genome. Although it has a higher cost per sample and can be more difficult to analyze than wes, greater diagnostic yields are evident in patients with negative or inconclusive wes results (Alfares et al., 2018; Lionel et al., 2018). WGS has been shown to identify a diagnostic genetic variant in ∼10–50% of individuals with a suspected genetic disorder, depending on the clinical study population(S-) being screened (van Der Ven et al., 2018; Groopman et al., 2019; Mann et al., 2019).
International Genetics and Translational Research in Transplantation Network
Despite technological advances that enable research to be carried out on a genome-wide scale, many studies have been hindered by small sample sizes in single transplant sites, as well as the vast number of complex donor and recipient clinical covariates and disease-related phenotypes observed in transplantation. The International Genetics & Translational Research in Transplantation Network (iGeneTRAiN) is a multi-site consortium that encompasses >45 genetic studies with ∼51,210 solid-organ transplant subjects (International Genetics and Translational Research in Transplantation Network (iGeneTRAiN), 2015). The iGeneTRAiN consortium aims to discover and validate solid organ transplant related genetic factors and post-transplant complications, including primary disease, disease recurrence, drug- and cardio-metabolic related adverse events, and different forms of allograft rejection (International Genetics and Translational Research in Transplantation Network (iGeneTRAiN), 2015). Of the iGeneTRAiN samples, 54% (n = 28,015) are from kidney transplant cohorts and include 17,742 (63.3%) recipients and 10,273 (36.7%) donors. The genotyped donor DNA provides control samples for all iGeneTRAiN studies, a large advantage over previously published genetic studies.
The iGeneTRAiN consortium designed and developed a genome-wide genotyping array, the “TxArray,” which was enriched with content relating to known or putative transplant-specific genetic associations (Li et al., 2015). The TxArray version 1 contains ∼782,000 genetic markers, with tailored transplant-specific content to capture variants across HLA, KIR, loss-of-function, pharmacogenomic, and cardio-metabolic loci. The array also contains extensive overlap with the UK Biobank Axiom® Array and the Axiom Biobank Genotyping Array, enabling future joint studies or meta-analyses using conventional, hypothesis-free GWAS approaches (Li et al., 2015).
The first wave of iGeneTRAiN kidney cohorts had a wide geographic representation with participants from various sites in the United States, Canada, Australia, The Netherlands, United Kingdom, and Ireland, including both adult and pediatric sites (Figure 1). Over the past few years, many genetic discoveries have been made within the iGeneTRAiN cohorts related to kidney, heart, liver, and lung transplants (Oetting et al., 2016; Greenland et al., 2017; Shaked et al., 2017; Hernandez-Fuentes et al., 2018; Oetting et al., 2018; Snoek et al., 2018; Oetting et al., 2019; Reindl-Schwaighofer et al., 2019). The Wellcome Trust Case Control Consortium (WTCCC) carried out the first large-scale GWAS with both kidney transplant donor and recipient DNA with the goal of identifying genetic variants, in addition to the HLA regions, that significantly contribute to long- and/or short-term renal allograft survival (Hernandez-Fuentes et al., 2018). No non-HLA signals were observed at genome-wide significance in this initial study, illustrating the need for harmonization of larger, well-phenotyped kidney transplant cohorts. In addition to the previously discovered common loss-of-function variant CYP3A5*3 allele (rs776746), the Deterioration of Kidney Allograft Function (DeKAF) Trial identified two CYP3A5 variants, rs10264272 and rs41303343, and one CYP3A4 variant, rs35599367, that explain additional portions of variance observed in dose-adjusted tacrolimus (TAC) through blood concentrations for African American (AA) and European ancestry (EA) kidney transplant recipients, respectively (Dai et al., 2006; Jacobson et al., 2011; Oetting et al., 2016; Oetting et al., 2018). These findings illustrate the utility of genome-wide studies when determining immunosuppression therapy regimens post-transplant, potentially contributing to improvements in renal allograft survival. Another iGeneTRAiN study showed that GWAS performed in nontransplant settings can predict post-transplant complications. Polygenic risk scores calculated from non-melanoma skin cancer (NMSC) GWAS in the general population predicted risk of and time to post-transplant NMSC and added additional predictive value beyond that explained by clinical variables (Stapleton et al., 2019).
Ongoing iGeneTRAiN Kidney Genome-Wide Studies
Recently, additional kidney transplant cohorts from Austria, Belgium, Germany, France, Italy, The Netherlands, Saudi Arabia, South Korea, Switzerland, and additional United States sites have joined iGeneTRAiN. This greatly increases ancestral diversity of recipients and donors, as well as statistical power to detect transplant related genetic variants that impact primary disease and transplant outcomes (Figure 1). Our large sample sizes are enabling us to investigate both donor and recipient characteristics that effect ESRD cause, treatment, and transplant-related outcomes.
Where available, we obtained formal clinical diagnoses of primary cause of ESRD, organized into disease categories of diabetic, arteriopathic, glomerular, acute kidney injury, infective and obstructive nephropathy, congenital, familial, toxic nephropathy, and malignancies, for all iGeneTRAiN kidney cohort subjects (Table 1). With these datasets, we are working to increase our understanding of the genetic underpinnings of ESRD and primary disease through single nucleotide polymorphism (SNP) based GWAS, copy number variant (CNV) screening, donor-recipient properties, allogenicity, and transplant outcomes.
Copy Number Variant Screening in iGeneTRAiN Cohorts
Genome-wide genotyping arrays are well established as an effective means for identification of known and novel CNVs (Sallustio et al., 2015; Ai et al., 2016; Verbitsky et al., 2019). CNV screening within iGeneTRAiN subjects is of major interest for both assessing the genetic architecture of primary disease and for allogenicity studies. iGeneTRAiN has developed an extensive loss-of-function (LoF) pipeline which includes haplotype phasing of over 10 million directly genotyped and imputed variants. We are particularly interested in two copy LoF (by single-nucleotide variants and/or CNVs) and integration of one or two copy LoF variants for donor-recipient interaction analyses, for association with time-to-rejection and graft loss events (International Genetics and Translational Research in Transplantation Network (iGeneTRAiN), 2015).
CNV screening in a priori regions for primary disease has been performed in iGeneTRAiN cohorts. For example, we performed CNV screening in patients with nephronophthisis (NPH), the most common genetic cause of ESRD in children and often caused by homozygous NPHP1 full gene deletions (Levy and Feingold, 2000; Hildebrandt, 2010; Wolf & Hildebrandt, 2011). In iGeneTRAiN, we previously examined this region in a subset of iGeneTRAiN studies for adult-onset ESRD (n = 5,606 patients). Of the subjects analyzed, 26 patients showed homozygous NPHP1 CNV deletions. Interestingly only 12% of these patients were previously diagnosed as having NPH and many presented with ESRD later in adulthood (Snoek et al., 2018). Thus, using the two copy gene loss of NPHP1 from GWG arrays to ascertain NPH status and examine NPH-related information in iGeneTRAiN studies, including accuracy of case-ascertainment and age-of-onset, shows a strong proof-of-principle for use in other high penetrant autosomal recessive/dominant cases, and the need for further sequencing for rare single-nucleotide variants in adult-onset ESRD patients. Furthermore, in a recent genome-wide analysis of CNVs in almost 3,000 cases of CAKUT, 45 distinct, known genomic disorders at 37 independent genomic loci were identified in 4% of CAKUT cases, and novel genomic disorders were found in an additional ∼2% of cases (Verbitsky et al., 2019). Genome-wide genotyping and imputation using large whole-genome sequencing (WGS) datasets, such as the 1000 genomes project (1KGP), typically cannot identify variants in the most common ancestral populations to a minor allele frequency (MAF) of <0.005, yet it is often possible to identify rare CNVs using monomorphic or SNP based probes across loci.
Our previous analyses of the Axiom TxArray genome-wide genotyping data was primarily limited to approximately 2,000 a priori CNV regions of interest that had specific probes designed onto the TxArray. Initial analyses used an adaption of the BRLMM-P CNV algorithm adapted from algorithms previously used to cluster genotypes across many samples (Yeung et al., 2008). However, BRLMM-P could only identify up to three clusters and thus was only able to detect 0, 1, and 2 copy deletions. The newer Axiom Analysis Suite 4.0 software allows streamlined, targeted, and de-novo whole-genome CNV region analysis3. A major advantage of the newer software is the ability to detect duplications as well as 0, 1, and 2 copy deletions (Figure 2).
Figure 2 (A) Copy Number Variant Region Analyses of the NPHP1 locus using the Axiom Analysis Suite 4.0 to call diploid state (CN=2), one copy deletion state (CN=1), two copy deletion state (CN=0), and duplication state (CN=3). The x-axis shows the chromosomal location on Chromosome 2 while the y-axis shows the standard Log2Ratio intensity. (B) In depth illustration of the duplication (CN=3) state.
Genome-wide genotyping studies have become very affordable and streamlined. However, large sample sizes, on the order of 10,000–100,000, are needed in order to detect both rare variants with large contributions and common variants with minor contributions to a specific phenotype(s) (Korte and Farlow, 2013). While it is very important to bolster statistical power to detect genetic underpinnings of transplant-related phenotypes by aggregating similar cohorts, great caution must be exercised when combining genotyping and phenotyping datasets, especially as transplant study covariates are very complex and can vary greatly by era and geographical region. iGeneTRAiN does have a unified quality control/quality assurance GWAS pipeline, including adjustment for population-based stratification (International Genetics and Translational Research in Transplantation Network (iGeneTRAiN), 2015). Association study analyses do adjust for all known/available study covariates, including patient demographics and clinical characteristics, and we adjust for each transplant site alone to look for confounders. Genome-wide genotyping arrays are generally poor at detecting rarer frequency pathogenic variants, with the exception of medium to large CNVs. Significant advances in genomic technologies and the decreasing cost of WES/WGS efforts over the past several years have made it increasingly feasible to carry out better designed genome-wide studies in a clinical environment (Gumpinger et al., 2018). However, there are still significant advantages to having genome-wide genotyping array datasets, as rigorous quality control and quality assurance measures are generally performed on the original DNA, and gender, ancestry, and HLA (amino acid imputation) concordance checks can be performed before progressing to WES or WGS pipelines for deeper genetic characterization. GWAS are able to provide insight into genetic risk scores and pathogenic CNVs, as genome-wide variants are covered in conventional genome-wide genotyping arrays (Sampson and Juppner, 2013; Li et al., 2015 Marigorta et al., 2018; Snoek et al., 2018; Stapleton et al., 2019). For example, a meta-analysis across 36 articles identified three genetic variants that are significantly associated with new onset diabetes after transplantation (NODAT), all of which are also known risk factor variants for Type 2 diabetes. The integration and analysis of large and complex multi-omic datasets has been demonstrated in a number of recent high impact publications, which in general increase, by approximately 10-fold, the statistical power to detect and illustrate functional variants (Chen et al., 2012; Piening et al., 2018; Schüssler-Fiorenza Rose et al., 2019; Zhou et al., 2019). iGeneTRAiN genomic data can be integrated with results from proteome-, metabolome-, and transcriptome-wide transplant studies to further characterize clinical risks and allow for personalized treatments, as a number of iGeneTRAiN studies have multi-omic datasets/samples (International Genetics and Translational Research in Transplantation Network (iGeneTRAiN), 2015).
The advent of single-cell RNA sequencing (scRNASeq) has yielded major insights into the biology of CKD. Expression quantitative trait loci (eQTL) atlases have been generated for glomerular and tubular compartments from human kidney cells. Integrating results from genome-wide studies of CKD with eQTL from scRNAseq as well as known regulatory region maps has been shown to identify novel CKD genes (Qiu et al., 2018). The Human Cell Atlas project is a major international initiative which aims to create comprehensive reference maps of all human cells to gain fundamental insight into the understanding of human health and will undoubtedly aid in the diagnoses and surveillance of a range of diseases (Regev et al., 2017).
Future of iGeneTRAiN Kidney Cohorts Analyses
As the population of kidney transplant recipients and donors continues to grow in the iGeneTRAiN consortium and as post-transplant outcomes accrue, we will be able to further increase our knowledge of the genetic underpinnings of ESRD, primary disease, and post-transplant outcomes, such as acute rejection and graft loss. These sequencing approaches may provide additional insight into donor-recipient (D-R) interactions that influence graft outcomes. Although it is well established that allelic matches across HLA loci impact clinical outcomes post-transplant, there is a paucity of genome-wide research conducted to identify donor-recipient interactions independent of HLA (Thorsby, 2009; Chan-On and Sarwal, 2016;Stapleton et al., 2018). One recent iGeneTRAiN kidney D-R study showed decreased allograft survival of recipients with increased D-R kidney transmembrane non-synonymous SNPs (nsSNPs). We further demonstrated that we could detect alloantibodies against customized amino-acid peptides designed with a number of these kidney transmembrane nsSNPs using sera from these patients (Reindl-Schwaighofer, et al., 2019). Finally, data from all solid-organ transplant studies in the iGeneTRAiN consortium will be utilized in cross-organ studies in order to gain additional insight into the genetics of acute rejection, allograft/patient survival, and pharmacogenomic outcomes.
All data used in this publication was collected in accordance with local IRB stipulations.
CF, MM, JS, FZ, LC, CW, SD, AA, TG, SL KK, ML and BK all provided data relating to their respective cohorts and all read and provided feedback on the manuscript. BK, MM, and CF performed CNV analyses.
Support was received from the Philadelphia Gift-of-Life Organ Procurement Organization.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank the Gift of Life Organ Procurement Organization, Philadelphia for funding which enabled this research.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.01084/full#supplementary-material
Ai, Z., Li, M., Liu, W., Foo, J. N., Mansouri, O., Yin, P., et al. (2016). Low alpha-defensin gene copy number increases the risk for IgA nephropathy and renal dysfunction. Sci. Transl. Med. 8 (345), 345ra388. doi: 10.1126/scitranslmed.aaf2106
Alfares, A., Aloraini, T., Subaie, L. A., Alissa, A., Qudsi, A. A., Alahmad, A., et al. (2018). Whole-genome sequencing offers additional but limited clinical utility compared with reanalysis of whole-exome sequencing. Genet. Med. 20 (11), 1328–1333. doi: 10.1038/gim.2018.41
Boger, C. A., Gorski, M., Li, M., Hoffmann, M. M., Huang, C., Yang, Q., et al. (2011b). Association of eGFR-Related Loci Identified by GWAS with Incident CKD and ESRD. PloS Genet. 7 (9), e1002292. doi: 10.1371/journal.pgen.1002292
Brodin-Sartorius, A., Tete, M. J., Niaudet, P., Antignac, C., Guest, G., Ottolenghi, C., et al. (2012). Cysteamine therapy delays the progression of nephropathic cystinosis in late adolescents and adults. Kidney Int. 81 (2), 179–189. doi: 10.1038/ki.2011.277
Calabrese, D. R., Florez, R., Dewey, K., Hui, C., Torgerson, D., Chong, T., et al. (2018). Genotypes associated with tacrolimus pharmacokinetics impact clinical outcomes in lung transplant recipients. Clin. Transplant. 32 (8), e13332. doi: 10.1111/ctr.13332
Chen, R., Mias, G. I., Pook-Than, J., Jiang, L., Lam, H. Y., Chen, R., et al. (2012). Personal omics profiling reveals dynamic molecular and medical phenotypes. Cell 148 (6), 1293–1307. doi: 10.1016/j.cell.2012.02.009
Cohen, D. J., St Martin, L., Christensen, L. L., Bloom, R. D., Sung, R. S. (2006). Kidney and pancreas transplantation in the United States, 1995-2004. Am. J. Transplant. 6 (5 Pt 2), 1153–1169. doi: 10.1111/j.1600-6143.2006.01272.x
Dai, Y., Hebert, M. F., Isoherranen, N., Davis, C. L., Marsh, C., Shen, D. D., et al. (2006). Effect of CYP3A5 polymorphism on tacrolimus metabolic clearance in vitro. Drug Metab. Dispos 34 (5), 836–847. doi: 10.1124/dmd.105.008680
Gillion, V., Dahan, K., Cosyns, J. P., Hilbert, P., Jadoul, M., Goffin, E., et al. (2018). Genotype and outcome after kidney transplantation in alport syndrome. Kidney Int. Rep. 3 (3), 652–660. doi: 10.1016/j.ekir.2018.01.008
Gorski, M., Tin, A., Garnaas, M., McMahon, G. M., Chu, A. Y., Tayo, B. O., et al. (2015). Genome-wide association study of kidney function decline in individuals of European descent. Kidney Int. 87 (5), 1017–1029. doi: 10.1038/ki.2014.361
Greenland, J. R., Sun, H., Calabrese, D., Chong, T., Singer, J. P., Kukreja, J., et al. (2017). HLA mismatching favoring host-versus-graft NK cell activity via KIR3DL1 is associated with improved outcomes following lung transplantation. Am. J. Transplant. 17 (8), 2192–2199. doi: 10.1111/ajt.14295
Groopman, E. E., Marasa, M., Cameron-Christie, S., Petrovski, S., Aggarwal, V. S., Milo-Rasouly, H., et al. (2019). Diagnostic utility of exome sequencing for kidney disease. N Engl. J. Med. 380 (2), 142–151. doi: 10.1056/NEJMoa1806891
Hernandez-Fuentes, M. P., Franklin, C., Rebollo-Mesa, I., Mollon, J., Delaney, F., Perucha, E., et al. (2018). Long- and short-term outcomes in renal allografts with deceased donors: a large recipient and donor genome-wide association study. Am. J. Transplant. 18 (6), 1370–1379. doi: 10.1111/ajt.14594
International Genetics and Translational Research in Transplantation Network (iGeneTRAiN). (2015). Design and Implementation of the International Genetics and Translational Research in Transplantation Network. Transplantation 99 (11), 2401–2412. doi: 10.1097/TP.0000000000000913
Jacobson, P. A., Oetting, W. S., Brearley, A. M., Leduc, R., Guan, W., Schladt, D., et al. (2011). Novel polymorphisms associated with tacrolimus trough concentrations: results from a multicenter kidney transplant consortium. Transplantation 91 (3), 300–308. doi: 10.1097/TP.0b013e318200e991
Kalatharan, V., Lemaire, M., Lanktree, M. B. (2018). Opportunities and challenges for genetic studies of end-stage renal disease in Canada. Can. J. Kidney Health Dis. 5, 2054358118789368. doi: 10.1177/2054358118789368
Kashiwagi, Y., Suzuki, S., Agata, K., Morishima, Y., Inagaki, N., Numabe, H., et al. (2018). A family case of X-linked Alport syndrome patients with a novel variant in COL4A5. CEN Case Rep. 8 (2), 75–78. doi: 10.1007/s13730-018-0368-4
Kopp, J. B., Nelson, G. W., Sampath, K., Johnson, R. C., Genovese, G., An, P., et al. (2011). APOL1 genetic variants in focal segmental glomerulosclerosis and HIV-associated nephropathy. J. Am. Soc. Nephrol. 22 (11), 2129–2137. doi: 10.1681/ASN.2011040388
Lata, S., Marasa, M., Li, Y., Fasel, D. A., Groopman, E., Jobanputra, V., et al. (2018). Whole-Exome Sequencing in Adults With Chronic Kidney Disease: A Pilot Study. Ann. Intern Med. 168 (2), 100–109. doi: 10.7326/M17-1319
Lei, T. Y., Fu, F., Li, R., Wang, D., Wang, R. Y., Jing, X. Y., et al. (2017). Whole-exome sequencing for prenatal diagnosis of fetuses with congenital anomalies of the kidney and urinary tract. Nephrol. Dial Transplant. 32 (10), 1665–1675. doi: 10.1093/ndt/gfx031
Li, Y. R., van Setten, J., Verma, S. S., Lu, Y., Holmes, M. V., Gao, H., et al. (2015). Concept and design of a genome-wide association genotyping array tailored for transplantation-specific studies. Genome Med. 7, 90. doi: 10.1186/s13073-015-0211-x
Limou, S., Nelson, G. W., Kopp, J. B., Winkler, C. A. (2014). APOL1 kidney risk alleles: population genetics and disease associations. Adv. Chronic Kidney Dis. 21 (5), 426–433. doi: 10.1053/j.ackd.2014.06.005
Lionel, A. C., Costain, G., Monfared, N., Walker, S., Reuter, M. S., Hosseini, S. M., et al. (2018). Improved diagnostic yield compared with targeted gene sequencing panels suggests a role for whole-genome sequencing as a first-tier genetic test. Genet. Med. 20 (4), 435–443. doi: 10.1038/gim.2017.119
Mallawaarachchi, A. C., Hort, Y., Cowley, M. J., McCabe, M. J., Minoche, A., Dinger, M. E., et al. (2016). Whole-genome sequencing overcomes pseudogene homology to diagnose autosomal dominant polycystic kidney disease. Eur. J. Hum. Genet. 24 (11), 1584–1590. doi: 10.1038/ejhg.2016.48
Mann, N., Braun, D. A., Amann, K., Tan, W., Shril, S., Connaughton, D. M., et al. (2019). whole-exome sequencing enables a precision medicine approach for kidney transplant recipients. J. Am. Soc. Nephrol. 30 (2), 201–215. doi: 10.1681/ASN.2018060575
Misra, M. K., Kapoor, R., Pandey, S. K., Sharma, R. K., Agrawal, S. (2014). Association of CTLA-4 gene polymorphism with end-stage renal disease and renal allograft outcome. J. Interferon Cytokine Res. 34 (3), 148–161. doi: 10.1089/jir.2013.0069
Oetting, W. S., Schladt, D. P., Guan, W., Miller, M. B., Remmel, R. P., Dorr, C., et al. (2016). Genomewide association study of tacrolimus concentrations in african american kidney transplant recipients identifies multiple CYP3A5 alleles. Am. J. Transplant. 16 (2), 574–582. doi: 10.1111/ajt.13495
Oetting, W. S., Wu, B., Schladt, D. P., Guan, W., Remmel, R. P., Mannon, R. B., et al. (2018). Genome-wide association study identifies the common variants in CYP3A4 and CYP3A5 responsible for variation in tacrolimus trough concentration in Caucasian kidney transplant recipients. Pharmacogenomics J. 18 (3), 501–505. doi: 10.1038/tpj.2017.49
Oetting, W. S., Wu, B., Schladt, D. P., Guan, W., van Setten, J., Keating, B. J., et al. (2019). Genetic variants associated with immunosuppressant pharmacokinetics and adverse effects in the DeKAF genomics genome wide association studies. Transplantation 103 (6), 1131–1139. doi: 10.1097/TP.0000000000002625
Okamoto, K., Rausch, J. W., Wakashin, H., Fu, Y., Chung, J. Y., Dummer, P. D., et al. (2018). APOL1 risk allele RNA contributes to renal toxicity by activating protein kinase R. Commun. Biol. 1, 188. doi: 10.1038/s42003-018-0188-2
Parsa, A., Kanetsky, P. A., Xiao, R., Gupta, J., Mitra, N., Limou, S., et al. (2017). Genome-wide association of CKD progression: the chronic renal insufficiency cohort study. J. Am. Soc. Nephrol. 28 (3), 923–934. doi: 10.1681/ASN.2015101152
Piening, B. D., Zhou, W., Contrepois, K., Röst, H., Gu Urban, G. J., Mishra, T., et al. (2018). Integrative personal omics profiles during periods of weight gain and loss. Cell Syst. 6 (2), 157–170. doi: 10.1016/j.cels.2017.12.013
Qiu, C., Huang, S., Park, J., Park, Y., Ko, Y. A., Seasock, M. J., et al. (2018). Renal compartment-specific genetic variation analyses identify new pathways in chronic kidney disease. Nat. Med. 24 (11), 1721–1731. doi: 10.1038/s41591-018-0194-4
Reindl-Schwaighofer, R., Heinzel, A., Kainz, A., van Setten, J., Jelencsics, K., Hu, K., et al. (2019). Contribution of non-HLA incompatibility between donor and recipient to kidney allograft survival: genome-wide analysis in a prospective cohort. Lancet 393 (10174), 910–917. doi: 10.1016/S0140-6736(18)32473-5
Reznichenko, A., Boger, C. A., Snieder, H., van den Born, J., de Borst, M. H., Damman, J., et al. (2012). UMOD as a susceptibility gene for end-stage renal disease. BMC Med. Genet. 13, 78. doi: 10.1186/1471-2350-13-78
Sallustio, F., Cox, S. N., Serino, G., Curci, C., Pesce, F., De Palma, G., et al. (2015). Genome-wide scan identifies a copy number variable region at 3p21.1 that influences the TLR9 expression levels in IgA nephropathy patients. Eur. J. Hum. Genet. 23 (7), 940–948. doi: 10.1038/ejhg.2014.208
Schüssler-Fiorenza Rose, S. M., Contrepois, K., Moneghetti, K. J., Zhou, K., Mishra, T., Mataraso, S., et al. (2019). A longitudinal big data approach for precision health. Nat. Med. 25 (5), 792–804. doi: 10.1038/s41591-019-0414-6
Serur, D., Saal, S., Wang, J., Sullivan, J., Bologa, R., Hartono, C., et al. (2011). Deceased-donor kidney transplantation: improvement in long-term survival. Nephrol. Dial Transplant. 26 (1), 317–324. doi: 10.1093/ndt/gfq415
Shaked, A., Chang, B. L., Barnes, M. R., Sayre, P., Li, Y. R., Asare, S., et al. (2017). An ectopically expressed serum miRNA signature is prognostic, diagnostic, and biologically related to liver allograft rejection. Hepatology 65 (1), 269–280. doi: 10.1002/hep.28786
Snoek, R., van Setten, J., Keating, B. J., Israni, A. K., Jacobson, P. A., Oetting, W. S., et al. (2018). NPHP1 (Nephrocystin-1) gene deletions cause adult-onset ESRD. J. Am. Soc. Nephrol. 29 (6), 1772–1779. doi: 10.1681/ASN.2017111200
Stapleton, C. P., Birdwell, K. A., McKnight, A. J., Maxwell, A. P., Mark, P. B., Sanders, M. L., et al. (2019). Polygenic risk score as a determinant of risk of non-melanoma skin cancer in a European-descent renal transplant cohort. Am. J. Transplant. 19 (3), 801–810. doi: 10.1111/ajt.15057
van der Ven, A. T., Connaughton, D. M., Ityel, H., Mann, N., Nakayama, M., Chen, J., et al. (2018). Whole-exome sequencing identifies causative mutations in families with congenital anomalies of the kidney and urinary tract. J. Am. Soc. Nephrol. 29 (9), 2348–2361. doi: 10.1681/ASN.2017121265
Verbitsky, M., Westland, R., Perez, A., Kiryluk, K., Liu, Q., Krithivasan, P., et al. (2019). The copy number variation landscape of congenital anomalies of the kidney and urinary tract. Nat. Genet. 51 (1), 117–127. doi: 10.1038/s41588-018-0281-y
Vignolini, G., Campi, R., Sessa, F., Greco, I., Larti, A., Giancane, S., et al. (2019). Development of a robot-assisted kidney transplantation programme from deceased donors in a referral academic centre: technical nuances and preliminary results. BJU Int. 123 (3), 474–484. doi: 10.1111/bju.14588
Warejko, J. K., Tan, W., Daga, A., Schapiro, D., Lawson, J. A., Shril, S., et al. (2018). Whole exome sequencing of patients with steroid-resistant nephrotic syndrome. Clin. J. Am. Soc. Nephrol. 13 (1), 53–62. doi: 10.2215/CJN.04120417
Wuhl, E., van Stralen, K. J., Wanner, C., Ariceta, G., Heaf, J. G., Bjerre, A. K., et al. (2014). Renal replacement therapy for rare diseases affecting the kidney: an analysis of the ERA-EDTA Registry. Nephrol. Dial Transplant. 29 Suppl 4, iv1–iv8. doi: 10.1093/ndt/gfu030
Yeung, J. M., Sham, P. C., Chan, A. S., Cherny, S. S. (2008). OpenADAM: an open source genome-wide association data management system for Affymetrix SNP arrays. BMC Genomics 9, 636. doi: 10.1186/1471-2164-9-636
Zhang, J., Qiu, W., Liu, H., Qian, C., Liu, D., Wang, H., et al. (2018). Genomic alterations in gastric cancers discovered via whole-exome sequencing. BMC Cancer 18 (1), 1270. doi: 10.1186/s12885-018-5097-8
Keywords: genomics, kidney disease, GWAS, whole exome sequencing analyses, whole genome sequencing
Citation: Fishman CE, Mohebnasab M, van Setten J, Zanoni F, Wang C, Deaglio S, Amoroso A, Callans L, van Gelder T, Lee S, Kiryluk K, Lanktree MB and Keating BJ (2019) Genome-Wide Study Updates in the International Genetics and Translational Research in Transplantation Network (iGeneTRAiN). Front. Genet. 10:1084. doi: 10.3389/fgene.2019.01084
Received: 26 March 2019; Accepted: 09 October 2019;
Published: 15 November 2019.
Edited by:Harvest F. Gu, China Pharmaceutical University, China
Reviewed by:Muhammad Jawad Hassan, National University of Medical Sciences (NUMS), Pakistan
Theodora Katsila, National Hellenic Research Foundation, Greece
Copyright © 2019 Fishman, Mohebnasab, van Setten, Zanoni, Wang, Deaglio, Amoroso, Callans, van Gelder, Lee, Kiryluk, Lanktree and Keating. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Brendan J. Keating, email@example.com