Developing a genetic testing panel for evaluation of morbidities in kidney transplant recipients.

Cardiovascular disease, infection, malignancy, and thromboembolism are major causes of morbidity and mortality in kidney transplant recipients (KTR). Prospectively identifying monogenic conditions associated with post-transplant complications may enable personalized management. Therefore, we developed a transplant morbidity panel (355 genes) associated with major post-transplant complications including cardiometabolic disorders, immunodeficiency, malignancy, and thrombophilia. This gene panel was then evaluated using exome sequencing data from 1590 KTR. Additionally, genes associated with monogenic kidney and genitourinary disorders along with American College of Medical Genetics (ACMG) secondary findings v3.2 were annotated. Altogether, diagnostic variants in 37 genes associated with Mendelian kidney and genitourinary disorders were detected in 9.9% (158/1590) of KTR; 25.9% (41/158) had not been clinically diagnosed. Moreover, the transplant morbidity gene panel detected diagnostic variants for 56 monogenic disorders in 9.1% KTRs (144/1590). Cardiovascular disease, malignancy, immunodeficiency, and thrombophilia variants were detected in 5.1% (81), 2.1% (34), 1.8% (29) and 0.2% (3) among 1590 KTRs, respectively. Concordant phenotypes were present in half of these cases. Reviewing implications for transplant care, these genetic findings would have allowed physicians to set specific risk factor targets in 6.3% (9/144), arrange intensive surveillance in 97.2% (140/144), utilize preventive measures in 13.2% (19/144), guide disease-specific therapy in 63.9% (92/144), initiate specialty referral in 90.3% (130/144) and alter immunosuppression in 56.9% (82/144). Thus, beyond diagnostic testing for kidney disorders, sequence annotation identified monogenic disorders associated with common post-transplant complications in 9.1% of KTR, with important clinical implications. Incorporating genetic diagnostics for transplant morbidities would enable personalized management in pre- and post-transplant care.

Genetic versus self-reported African ancestry of the recipient and neighborhood predictors of kidney transplantation outcomes in 2 multiethnic urban cohorts.

African American (AA) kidney recipients have a higher risk of allograft rejection and failure compared to non-AAs, but to what extent these outcomes are due to genetic versus environmental effects is currently unknown. Herein, we tested the effects of recipient self-reported race versus genetic proportion of African ancestry (pAFR), and neighborhood socioeconomic status (SES) on kidney allograft outcomes in multiethnic kidney transplant recipients from Columbia University (N = 1083) and the University of Pennsylvania (N = 738). All participants were genotyped with SNP arrays to estimate genetic admixture proportions. US census tract variables were used to analyze the effect of neighborhood factors. In both cohorts, self-reported recipient AA race and pAFR were individually associated with increased risk of rejection and failure after adjustment for known clinical risk factors and neighborhood SES factors. Joint analysis confirmed that self-reported recipient AA race and pAFR were both associated with a higher risk of allograft rejection (AA: HR 1.61 (1.31-1.96), P = 4.05E-06; pAFR: HR 1.90 (1.46-2.48), P = 2.40E-06) and allograft failure (AA: HR 1.52 (1.18-1.97), P = .001; pAFR: HR 1.70 (1.22-2.35), P = .002). Further research is needed to disentangle the role of genetics versus environmental, social, and structural factors contributing to poor transplantation outcomes in kidney recipients of African ancestry.

Genome-wide association analyses define pathogenic signaling pathways and prioritize drug targets for IgA nephropathy.

IgA nephropathy (IgAN) is a progressive form of kidney disease defined by glomerular deposition of IgA. Here we performed a genome-wide association study of 10,146 kidney-biopsy-diagnosed IgAN cases and 28,751 controls across 17 international cohorts. We defined 30 genome-wide significant risk loci explaining 11% of disease risk. A total of 16 loci were new, including TNFSF4/TNFSF18, REL, CD28, PF4V1, LY86, LYN, ANXA3, TNFSF8/TNFSF15, REEP3, ZMIZ1, OVOL1/RELA, ETS1, IGH, IRF8, TNFRSF13B and FCAR. The risk loci were enriched in gene orthologs causing abnormal IgA levels when genetically manipulated in mice. We also observed a positive genetic correlation between IgAN and serum IgA levels. High polygenic score for IgAN was associated with earlier onset of kidney failure. In a comprehensive functional annotation analysis of candidate causal genes, we observed convergence of biological candidates on a common set of inflammatory signaling pathways and cytokine ligand-receptor pairs, prioritizing potential new drug targets.

Rare Single Nucleotide and Copy Number Variants and the Etiology of Congenital Obstructive Uropathy: Implications for Genetic Diagnosis.

SIGNIFICANCE STATEMENT: Congenital obstructive uropathy (COU) is a prevalent human developmental defect with highly heterogeneous clinical presentations and outcomes. Genetics may refine diagnosis, prognosis, and treatment, but the genomic architecture of COU is largely unknown. Comprehensive genomic screening study of 733 cases with three distinct COU subphenotypes revealed disease etiology in 10.0% of them. We detected no significant differences in the overall diagnostic yield among COU subphenotypes, with characteristic variable expressivity of several mutant genes. Our findings therefore may legitimize a genetic first diagnostic approach for COU, especially when burdening clinical and imaging characterization is not complete or available. BACKGROUND: Congenital obstructive uropathy (COU) is a common cause of developmental defects of the urinary tract, with heterogeneous clinical presentation and outcome. Genetic analysis has the potential to elucidate the underlying diagnosis and help risk stratification. METHODS: We performed a comprehensive genomic screen of 733 independent COU cases, which consisted of individuals with ureteropelvic junction obstruction ( n =321), ureterovesical junction obstruction/congenital megaureter ( n =178), and COU not otherwise specified (COU-NOS; n =234). RESULTS: We identified pathogenic single nucleotide variants (SNVs) in 53 (7.2%) cases and genomic disorders (GDs) in 23 (3.1%) cases. We detected no significant differences in the overall diagnostic yield between COU sub-phenotypes, and pathogenic SNVs in several genes were associated to any of the three categories. Hence, although COU may appear phenotypically heterogeneous, COU phenotypes are likely to share common molecular bases. On the other hand, mutations in TNXB were more often identified in COU-NOS cases, demonstrating the diagnostic challenge in discriminating COU from hydronephrosis secondary to vesicoureteral reflux, particularly when diagnostic imaging is incomplete. Pathogenic SNVs in only six genes were found in more than one individual, supporting high genetic heterogeneity. Finally, convergence between data on SNVs and GDs suggest MYH11 as a dosage-sensitive gene possibly correlating with severity of COU. CONCLUSIONS: We established a genomic diagnosis in 10.0% of COU individuals. The findings underscore the urgent need to identify novel genetic susceptibility factors to COU to better define the natural history of the remaining 90% of cases without a molecular diagnosis.

Clinical and Genetic Characteristics of CKD Patients with High-Risk APOL1 Genotypes.

SIGNIFICANCE STATEMENT: APOL1 high-risk genotypes confer a significant risk of kidney disease, but variability in patient outcomes suggests the presence of modifiers of the APOL1 effect. We show that a diverse population of CKD patients with high-risk APOL1 genotypes have an increased lifetime risk of kidney failure and higher eGFR decline rates, with a graded risk among specific high-risk genotypes. CKD patients with high-risk APOL1 genotypes have a lower diagnostic yield for monogenic kidney disease. Exome sequencing revealed enrichment of rare missense variants within the inflammasome pathway modifying the effect of APOL1 risk genotypes, which may explain some clinical heterogeneity. BACKGROUND: APOL1 genotype has significant effects on kidney disease development and progression that vary among specific causes of kidney disease, suggesting the presence of effect modifiers. METHODS: We assessed the risk of kidney failure and the eGFR decline rate in patients with CKD carrying high-risk ( N =239) and genetically matched low-risk ( N =1187) APOL1 genotypes. Exome sequencing revealed monogenic kidney diseases. Exome-wide association studies and gene-based and gene set-based collapsing analyses evaluated genetic modifiers of the effect of APOL1 genotype on CKD. RESULTS: Compared with genetic ancestry-matched patients with CKD with low-risk APOL1 genotypes, those with high-risk APOL1 genotypes had a higher risk of kidney failure (Hazard Ratio [HR]=1.58), a higher decline in eGFR (6.55 versus 3.63 ml/min/1.73 m 2 /yr), and were younger at time of kidney failure (45.1 versus 53.6 years), with the G1/G1 genotype demonstrating the highest risk. The rate for monogenic kidney disorders was lower among patients with CKD with high-risk APOL1 genotypes (2.5%) compared with those with low-risk genotypes (6.7%). Gene set analysis identified an enrichment of rare missense variants in the inflammasome pathway in individuals with high-risk APOL1 genotypes and CKD (odds ratio=1.90). CONCLUSIONS: In this genetically matched cohort, high-risk APOL1 genotypes were associated with an increased risk of kidney failure and eGFR decline rate, with a graded risk between specific high-risk genotypes and a lower rate of monogenic kidney disease. Rare missense variants in the inflammasome pathway may act as genetic modifiers of APOL1 effect on kidney disease.

Copy Number Variant Analysis and Genome-wide Association Study Identify Loci with Large Effect for Vesicoureteral Reflux.

BACKGROUND: Vesicoureteral reflux (VUR) is a common, familial genitourinary disorder, and a major cause of pediatric urinary tract infection (UTI) and kidney failure. The genetic basis of VUR is not well understood. METHODS: A diagnostic analysis sought rare, pathogenic copy number variant (CNV) disorders among 1737 patients with VUR. A GWAS was performed in 1395 patients and 5366 controls, of European ancestry. RESULTS: Altogether, 3% of VUR patients harbored an undiagnosed rare CNV disorder, such as the 1q21.1, 16p11.2, 22q11.21, and triple X syndromes ((OR, 3.12; 95% CI, 2.10 to 4.54; P=6.35×10-8) The GWAS identified three study-wide significant and five suggestive loci with large effects (ORs, 1.41-6.9), containing canonical developmental genes expressed in the developing urinary tract (WDPCP, OTX1, BMP5, VANGL1, and WNT5A). In particular, 3.3% of VUR patients were homozygous for an intronic variant in WDPCP (rs13013890; OR, 3.65; 95% CI, 2.39 to 5.56; P=1.86×10-9). This locus was associated with multiple genitourinary phenotypes in the UK Biobank and eMERGE studies. Analysis of Wnt5a mutant mice confirmed the role of Wnt5a signaling in bladder and ureteric morphogenesis. CONCLUSIONS: These data demonstrate the genetic heterogeneity of VUR. Altogether, 6% of patients with VUR harbored a rare CNV or a common variant genotype conferring an OR >3. Identification of these genetic risk factors has multiple implications for clinical care and for analysis of outcomes in VUR.

The genetic architecture of membranous nephropathy and its potential to improve non-invasive diagnosis.

Membranous Nephropathy (MN) is a rare autoimmune cause of kidney failure. Here we report a genome-wide association study (GWAS) for primary MN in 3,782 cases and 9,038 controls of East Asian and European ancestries. We discover two previously unreported loci, NFKB1 (rs230540, OR = 1.25, P = 3.4 × 10-12) and IRF4 (rs9405192, OR = 1.29, P = 1.4 × 10-14), fine-map the PLA2R1 locus (rs17831251, OR = 2.25, P = 4.7 × 10-103) and report ancestry-specific effects of three classical HLA alleles: DRB1*1501 in East Asians (OR = 3.81, P = 2.0 × 10-49), DQA1*0501 in Europeans (OR = 2.88, P = 5.7 × 10-93), and DRB1*0301 in both ethnicities (OR = 3.50, P = 9.2 × 10-23 and OR = 3.39, P = 5.2 × 10-82, respectively). GWAS loci explain 32% of disease risk in East Asians and 25% in Europeans, and correctly re-classify 20-37% of the cases in validation cohorts that are antibody-negative by the serum anti-PLA2R ELISA diagnostic test. Our findings highlight an unusual genetic architecture of MN, with four loci and their interactions accounting for nearly one-third of the disease risk.

Genomic Mismatch at LIMS1 Locus and Kidney Allograft Rejection.

BACKGROUND: In the context of kidney transplantation, genomic incompatibilities between donor and recipient may lead to allosensitization against new antigens. We hypothesized that recessive inheritance of gene-disrupting variants may represent a risk factor for allograft rejection. METHODS: We performed a two-stage genetic association study of kidney allograft rejection. In the first stage, we performed a recessive association screen of 50 common gene-intersecting deletion polymorphisms in a cohort of kidney transplant recipients. In the second stage, we replicated our findings in three independent cohorts of donor-recipient pairs. We defined genomic collision as a specific donor-recipient genotype combination in which a recipient who was homozygous for a gene-intersecting deletion received a transplant from a nonhomozygous donor. Identification of alloantibodies was performed with the use of protein arrays, enzyme-linked immunosorbent assays, and Western blot analyses. RESULTS: In the discovery cohort, which included 705 recipients, we found a significant association with allograft rejection at the LIMS1 locus represented by rs893403 (hazard ratio with the risk genotype vs. nonrisk genotypes, 1.84; 95% confidence interval [CI], 1.35 to 2.50; P = 9.8×10-5). This effect was replicated under the genomic-collision model in three independent cohorts involving a total of 2004 donor-recipient pairs (hazard ratio, 1.55; 95% CI, 1.25 to 1.93; P = 6.5×10-5). In the combined analysis (discovery cohort plus replication cohorts), the risk genotype was associated with a higher risk of rejection than the nonrisk genotype (hazard ratio, 1.63; 95% CI, 1.37 to 1.95; P = 4.7×10-8). We identified a specific antibody response against LIMS1, a kidney-expressed protein encoded within the collision locus. The response involved predominantly IgG2 and IgG3 antibody subclasses. CONCLUSIONS: We found that the LIMS1 locus appeared to encode a minor histocompatibility antigen. Genomic collision at this locus was associated with rejection of the kidney allograft and with production of anti-LIMS1 IgG2 and IgG3. (Funded by the Columbia University Transplant Center and others.).

Author Correction: The copy number variation landscape of congenital anomalies of the kidney and urinary tract.

In the version of this article initially published, affiliation 38 incorrectly read "ICNU-Nephrology and Urology Department, Barcelona, Spain"; "Renal Division, Hospital Clinic, IDIBAPS, University of Barcelona, Barcelona, Spain" is the correct affiliation. The error has been corrected in the HTML and PDF versions of the article.

The copy number variation landscape of congenital anomalies of the kidney and urinary tract.

Congenital anomalies of the kidney and urinary tract (CAKUT) are a major cause of pediatric kidney failure. We performed a genome-wide analysis of copy number variants (CNVs) in 2,824 cases and 21,498 controls. Affected individuals carried a significant burden of rare exonic (that is, affecting coding regions) CNVs and were enriched for known genomic disorders (GD). Kidney anomaly (KA) cases were most enriched for exonic CNVs, encompassing GD-CNVs and novel deletions; obstructive uropathy (OU) had a lower CNV burden and an intermediate prevalence of GD-CNVs; and vesicoureteral reflux (VUR) had the fewest GD-CNVs but was enriched for novel exonic CNVs, particularly duplications. Six loci (1q21, 4p16.1-p16.3, 16p11.2, 16p13.11, 17q12 and 22q11.2) accounted for 65% of patients with GD-CNVs. Deletions at 17q12, 4p16.1-p16.3 and 22q11.2 were specific for KA; the 16p11.2 locus showed extensive pleiotropy. Using a multidisciplinary approach, we identified TBX6 as a driver for the CAKUT subphenotypes in the 16p11.2 microdeletion syndrome.