Thrombosis risk in single- and double-heterozygous carriers of factor V Leiden and prothrombin G20210A in FinnGen and the UK Biobank

ABSTRACT

The factor V Leiden (FVL; rs6025) and prothrombin G20210A (PTGM; rs1799963) polymorphisms are 2 of the most well-studied genetic risk factors for venous thromboembolism (VTE). However, double heterozygosity (DH) for FVL and PTGM remains poorly understood, with previous studies showing marked disagreement regarding thrombosis risk conferred by the DH genotype. Using multidimensional data from the UK Biobank (UKB) and FinnGen biorepositories, we evaluated the clinical impact of DH carrier status across 937 939 individuals. We found that 662 participants (0.07%) were DH carriers. After adjustment for age, sex, and ancestry, DH individuals experienced a markedly elevated risk of VTE compared with wild-type individuals (odds ratio [OR] = 5.24; 95% confidence interval [CI], 4.01-6.84; P = 4.8 × 10-34), which approximated the risk conferred by FVL homozygosity. A secondary analysis restricted to UKB participants (N = 445 144) found that effect size estimates for the DH genotype remained largely unchanged (OR = 4.53; 95% CI, 3.42-5.90; P < 1 × 10-16) after adjustment for commonly cited VTE risk factors, such as body mass index, blood type, and markers of inflammation. In contrast, the DH genotype was not associated with a significantly higher risk of any arterial thrombosis phenotype, including stroke, myocardial infarction, and peripheral artery disease. In summary, we leveraged population-scale genomic data sets to conduct, to our knowledge, the largest study to date on the DH genotype and were able to establish far more precise effect size estimates than previously possible. Our findings indicate that the DH genotype may occur as frequently as FVL homozygosity and may confer a similarly increased risk of VTE.

Risk of meningomyelocele mediated by the common 22q11.2 deletion

Meningomyelocele is one of the most severe forms of neural tube defects (NTDs) and the most frequent structural birth defect of the central nervous system. We assembled the Spina Bifida Sequencing Consortium to identify causes. Exome and genome sequencing of 715 parent-offspring trios identified six patients with chromosomal 22q11.2 deletions, suggesting a 23-fold increased risk compared with the general population. Furthermore, analysis of a separate 22q11.2 deletion cohort suggested a 12- to 15-fold increased NTD risk of meningomyelocele. The loss of Crkl, one of several neural tube-expressed genes within the minimal deletion interval, was sufficient to replicate NTDs in mice, where both penetrance and expressivity were exacerbated by maternal folate deficiency. Thus, the common 22q11.2 deletion confers substantial meningomyelocele risk, which is partially alleviated by folate supplementation.

The expanded spectrum of human disease associated with GREB1L likely includes complex congenital heart disease

OBJECTIVE

GREB1L has been linked prenatally to Potter's sequence, as well as less severe anomalies of the kidney, uterus, inner ear, and heart. The full phenotypic spectrum is unknown. The purpose of this study was to characterize known and novel pre- and postnatal phenotypes associated with GREB1L.

METHODS

We solicited cases from the Fetal Sequencing Consortium, screened a population-based genomic database, and conducted a comprehensive literature search to identify disease cases associated with GREB1L. We present a detailed phenotypic spectrum and molecular changes.

RESULTS

One hundred twenty-seven individuals with 51 unique pathogenic or likely pathogenic GREB1L variants were identified. 24 (47%) variants were associated with isolated kidney anomalies, 19 (37%) with anomalies of multiple systems, including one case of hypoplastic left heart syndrome, five (10%) with isolated sensorineural hearing loss, two (4%) with isolated uterine agenesis; and one (2%) with isolated tetralogy of Fallot.

CONCLUSION

GREB1L may cause complex congenital heart disease (CHD) in humans. Clinicians should consider GREB1L testing in the setting of CHD, and cardiac screening in the setting of GREB1L variants.

Copy number variation analysis in 138 families with steroid-resistant nephrotic syndrome identifies causal homozygous deletions in PLCE1 and NPHS2 in two families

BACKGROUND

Steroid-resistant nephrotic syndrome (SRNS) is the second most common cause of kidney failure in children and adults under the age of 20 years. Previously, we were able to detect by exome sequencing (ES) a known monogenic cause of SRNS in 25-30% of affected families. However, ES falls short of detecting copy number variants (CNV). Therefore, we hypothesized that causal CNVs could be detected in a large SRNS cohort.

METHODS

We performed genome-wide single nucleotide polymorphism (SNP)-based CNV analysis on a cohort of 138 SRNS families, in whom we previously did not identify a genetic cause through ES. We evaluated ES and CNV data for variants in 60 known SRNS genes and in 13 genes in which variants are known to cause a phenocopy of SRNS. We applied previously published, predefined criteria for CNV evaluation.

RESULTS

We detected a novel CNV in two genes in 2 out of 138 families (1.5%). The 9,673 bp homozygous deletion in PLCE1 and the 6,790 bp homozygous deletion in NPHS2 were confirmed across the breakpoints by PCR and Sanger sequencing.

CONCLUSIONS

We confirmed that CNV analysis can identify the genetic cause in SRNS families that remained unsolved after ES. Though the rate of detected CNVs is minor, CNV analysis can be used when there are no other genetic causes identified. Causative CNVs are less common in SRNS than in other monogenic kidney diseases, such as congenital anomalies of the kidneys and urinary tract, where the detection rate was 5.3%. A higher resolution version of the Graphical abstract is available as Supplementary information.

Mouse and human studies support DSTYK loss of function as a low-penetrance and variable expressivity risk factor for congenital urinary tract anomalies

PURPOSE

Previous work identified rare variants in DSTYK associated with human congenital anomalies of the kidney and urinary tract (CAKUT). Here, we present a series of mouse and human studies to clarify the association, penetrance, and expressivity of DSTYK variants.

METHODS

We phenotypically characterized Dstyk knockout mice of 3 separate inbred backgrounds and re-analyzed the original family segregating the DSTYK c.654+1G>A splice-site variant (referred to as "SSV" below). DSTYK loss of function (LOF) and SSVs were annotated in individuals with CAKUT, epilepsy, or amyotrophic lateral sclerosis vs controls. A phenome-wide association study analysis was also performed using United Kingdom Biobank (UKBB) data.

RESULTS

Results demonstrate ∼20% to 25% penetrance of obstructive uropathy, at least, in C57BL/6J and FVB/NJ Dstyk-/- mice. Phenotypic penetrance increased to ∼40% in C3H/HeJ mutants, with mild-to-moderate severity. Re-analysis of the original family segregating the rare SSV showed low penetrance (43.8%) and no alternative genetic causes for CAKUT. LOF DSTYK variants burden showed significant excess for CAKUT and epilepsy vs controls and an exploratory phenome-wide association study supported association with neurological disorders.

CONCLUSION

These data support causality for DSTYK LOF variants and highlights the need for large-scale sequencing studies (here >200,000 cases) to accurately assess causality for genes and variants to lowly penetrant traits with common population prevalence.

Strong protective effect of the APOL1 p.N264K variant against G2-associated focal segmental glomerulosclerosis and kidney disease

African Americans have a significantly higher risk of developing chronic kidney disease, especially focal segmental glomerulosclerosis -, than European Americans. Two coding variants (G1 and G2) in the APOL1 gene play a major role in this disparity. While 13% of African Americans carry the high-risk recessive genotypes, only a fraction of these individuals develops FSGS or kidney failure, indicating the involvement of additional disease modifiers. Here, we show that the presence of the APOL1 p.N264K missense variant, when co-inherited with the G2 APOL1 risk allele, substantially reduces the penetrance of the G1G2 and G2G2 high-risk genotypes by rendering these genotypes low-risk. These results align with prior functional evidence showing that the p.N264K variant reduces the toxicity of the APOL1 high-risk alleles. These findings have important implications for our understanding of the mechanisms of APOL1-associated nephropathy, as well as for the clinical management of individuals with high-risk genotypes that include the G2 allele.

Strong protective effect of the APOL1 p.N264K variant against G2-associated focal segmental glomerulosclerosis and kidney disease

Black Americans have a significantly higher risk of developing chronic kidney disease (CKD), especially focal segmental glomerulosclerosis (FSGS), than European Americans. Two coding variants (G1 and G2) in the APOL1 gene play a major role in this disparity. While 13% of Black Americans carry the high-risk recessive genotypes, only a fraction of these individuals develops FSGS or kidney failure, indicating the involvement of additional disease modifiers. Here, we show that the presence of the APOL1 p.N264K missense variant, when co-inherited with the G2 APOL1 risk allele, substantially reduces the penetrance of the G1G2 and G2G2 high-risk genotypes by rendering these genotypes low-risk. These results align with prior functional evidence showing that the p.N264K variant reduces the toxicity of the APOL1 high-risk alleles. These findings have important implications for our understanding of the mechanisms of APOL1 -associated nephropathy, as well as for the clinical management of individuals with high-risk genotypes that include the G2 allele.

Implementation and Feasibility of Clinical Genome Sequencing Embedded Into the Outpatient Nephrology Care for Patients With Proteinuric Kidney Disease

Introduction

The diagnosis and management of proteinuric kidney diseases such as focal segmental glomerulosclerosis (FSGS) are challenging. Genetics holds the promise to improve clinical decision making for these diseases; however, it is often performed too late to enable timely clinical action and it is not implemented within routine outpatient nephrology visits.

Methods

We sought to test the implementation and feasibility of clinical rapid genome sequencing (GS) in guiding decision making in patients with proteinuric kidney disease in real-time and embedded in the outpatient nephrology setting.

Results

We enrolled 10 children or young adults with biopsy-proven FSGS (9 cases) or minimal change disease (1 case). The mean age at enrollment was 16.2 years (range 2-30). The workflow did not require referral to external genetics clinics but was conducted entirely during the nephrology standard-of-care appointments. The total turn-around-time from enrollment to return-of-results and clinical decision averaged 21.8 days (12.4 for GS), which is well within a time frame that allows clinically relevant treatment decisions. A monogenic or APOL1-related form of kidney disease was diagnosed in 5 of 10 patients. The genetic findings resulted in a rectified diagnosis in 6 patients. Both positive and negative GS findings determined a change in pharmacological treatment. In 3 patients, the results were instrumental for transplant evaluation, donor selection, and the immunosuppressive treatment. All patients and families received genetic counseling.

Conclusion

Clinical GS is feasible and can be implemented in real-time in the outpatient care to help guiding clinical management. Additional studies are needed to confirm the cost-effectiveness and broader utility of clinical GS across the phenotypic and demographic spectrum of kidney diseases.

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.

CERT1 mutations perturb human development by disrupting sphingolipid homeostasis

Neural differentiation, synaptic transmission, and action potential propagation depend on membrane sphingolipids, whose metabolism is tightly regulated. Mutations in the ceramide transporter CERT (CERT1), which is involved in sphingolipid biosynthesis, are associated with intellectual disability, but the pathogenic mechanism remains obscure. Here, we characterize 31 individuals with de novo missense variants in CERT1. Several variants fall into a previously uncharacterized dimeric helical domain that enables CERT homeostatic inactivation, without which sphingolipid production goes unchecked. The clinical severity reflects the degree to which CERT autoregulation is disrupted, and inhibiting CERT pharmacologically corrects morphological and motor abnormalities in a Drosophila model of the disease, which we call ceramide transporter (CerTra) syndrome. These findings uncover a central role for CERT autoregulation in the control of sphingolipid biosynthetic flux, provide unexpected insight into the structural organization of CERT, and suggest a possible therapeutic approach for patients with CerTra syndrome.

Multi-population genome-wide association study implicates immune and non-immune factors in pediatric steroid-sensitive nephrotic syndrome

Pediatric steroid-sensitive nephrotic syndrome (pSSNS) is the most common childhood glomerular disease. Previous genome-wide association studies (GWAS) identified a risk locus in the HLA Class II region and three additional independent risk loci. But the genetic architecture of pSSNS, and its genetically driven pathobiology, is largely unknown. Here, we conduct a multi-population GWAS meta-analysis in 38,463 participants (2440 cases). We then conduct conditional analyses and population specific GWAS. We discover twelve significant associations-eight from the multi-population meta-analysis (four novel), two from the multi-population conditional analysis (one novel), and two additional novel loci from the European meta-analysis. Fine-mapping implicates specific amino acid haplotypes in HLA-DQA1 and HLA-DQB1 driving the HLA Class II risk locus. Non-HLA loci colocalize with eQTLs of monocytes and numerous T-cell subsets in independent datasets. Colocalization with kidney eQTLs is lacking but overlap with kidney cell open chromatin suggests an uncharacterized disease mechanism in kidney cells. A polygenic risk score (PRS) associates with earlier disease onset. Altogether, these discoveries expand our knowledge of pSSNS genetic architecture across populations and provide cell-specific insights into its molecular drivers. Evaluating these associations in additional cohorts will refine our understanding of population specificity, heterogeneity, and clinical and molecular associations.

Genomic Disorders in CKD across the Lifespan

SIGNIFICANCE STATEMENT

Pathogenic structural genetic variants, also known as genomic disorders, have been associated with pediatric CKD. This study extends those results across the lifespan, with genomic disorders enriched in both pediatric and adult patients compared with controls. In the Chronic Renal Insufficiency Cohort study, genomic disorders were also associated with lower serum Mg, lower educational performance, and a higher risk of death. A phenome-wide association study confirmed the link between kidney disease and genomic disorders in an unbiased way. Systematic detection of genomic disorders can provide a molecular diagnosis and refine prediction of risk and prognosis.

BACKGROUND

Genomic disorders (GDs) are associated with many comorbid outcomes, including CKD. Identification of GDs has diagnostic utility.

METHODS

We examined the prevalence of GDs among participants in the Chronic Kidney Disease in Children (CKiD) cohort II ( n =248), Chronic Renal Insufficiency Cohort (CRIC) study ( n =3375), Columbia University CKD Biobank (CU-CKD; n =1986), and the Family Investigation of Nephropathy and Diabetes (FIND; n =1318) compared with 30,746 controls. We also performed a phenome-wide association analysis (PheWAS) of GDs in the electronic MEdical Records and GEnomics (eMERGE; n =11,146) cohort.

RESULTS

We found nine out of 248 (3.6%) CKiD II participants carried a GD, replicating prior findings in pediatric CKD. We also identified GDs in 72 out of 6679 (1.1%) adult patients with CKD in the CRIC, CU-CKD, and FIND cohorts, compared with 199 out of 30,746 (0.65%) GDs in controls (OR, 1.7; 95% CI, 1.3 to 2.2). Among adults with CKD, we found recurrent GDs at the 1q21.1, 16p11.2, 17q12, and 22q11.2 loci. The 17q12 GD (diagnostic of renal cyst and diabetes syndrome) was most frequent, present in 1:252 patients with CKD and diabetes. In the PheWAS, dialysis and neuropsychiatric phenotypes were the top associations with GDs. In CRIC participants, GDs were associated with lower serum magnesium, lower educational achievement, and higher mortality risk.

CONCLUSION

Undiagnosed GDs are detected both in children and adults with CKD. Identification of GDs in these patients can enable a precise genetic diagnosis, inform prognosis, and help stratify risk in clinical studies. GDs could also provide a molecular explanation for nephropathy and comorbidities, such as poorer neurocognition for a subset of patients.

Copy Number Variation Analysis Facilitates Identification of Genetic Causation in Patients with Congenital Anomalies of the Kidney and Urinary Tract

Background

Congenital anomalies of the kidneys and urinary tract (CAKUT) are the most common cause of chronic kidney disease among children and adults younger than 30 yr. In our previous study, whole-exome sequencing (WES) identified a known monogenic cause of isolated or syndromic CAKUT in 13% of families with CAKUT. However, WES has limitations and detection of copy number variations (CNV) is technically challenging, and CNVs causative of CAKUT have previously been detected in up to 16% of cases.

Objective

To detect CNVs causing CAKUT in this WES cohort and increase the diagnostic yield.

Design setting and participants

We performed a genome-wide single nucleotide polymorphism (SNP)-based CNV analysis on the same CAKUT cohort for whom WES was previously conducted.

Outcome measurements and statistical analysis

We evaluated and classified the CNVs using previously published predefined criteria.

Results and limitations

In a cohort of 170 CAKUT families, we detected a pathogenic CNV known to cause CAKUT in nine families (5.29%, 9/170). There were no competing variants on genome-wide CNV analysis or WES analysis. In addition, we identified novel likely pathogenic CNVs that may cause a CAKUT phenotype in three of the 170 families (1.76%).

Conclusions

CNV analysis in this cohort of 170 CAKUT families previously examined via WES increased the rate of diagnosis of genetic causes of CAKUT from 13% on WES to 18% on WES + CNV analysis combined. We also identified three candidate loci that may potentially cause CAKUT.

Patient summary

We conducted a genetics study on families with congenital anomalies of the kidney and urinary tract (CAKUT). We identified gene mutations that can explain CAKUT symptoms in 5.29% of the families, which increased the percentage of genetic causes of CAKUT to 18% from a previous study, so roughly one in five of our patients with CAKUT had a genetic cause. These analyses can help patients with CAKUT and their families in identifying a possible genetic cause.

Incorporating genetics services into adult kidney disease care

Studies have shown that as many as 1 in 10 adults with chronic kidney disease has a monogenic form of disease. However, genetic services in adult nephrology are limited. An adult Kidney Genetics Clinic was established within the nephrology division at a large urban academic medical center to increase access to genetic services and testing in adults with kidney disease. Between June 2019 and December 2021, a total of 363 patients were referred to the adult Kidney Genetics Clinic. Of those who completed genetic testing, a positive diagnostic finding was identified in 27.1%, a candidate diagnostic finding was identified in 6.7% of patients, and a nondiagnostic positive finding was identified in an additional 8.6% of patients, resulting in an overall yield of 42.4% for clinically relevant genetic findings in tested patients. A genetic diagnosis had implications for medical management, family member testing, and eligibility for clinical trials. With the utilization of telemedicine, genetic services reached a diverse geographic and patient population. Genetic education efforts were integral to the clinic's success, as they increased visibility and helped providers identify appropriate referrals. Ongoing access to genomic services will remain a fundamental component of patient care in adults with kidney disease.

Multidisciplinary approaches for elucidating genetics and molecular pathogenesis of urinary tract malformations

Advances in clinical diagnostics and molecular tools have improved our understanding of the genetically heterogeneous causes underlying congenital anomalies of kidney and urinary tract (CAKUT). However, despite a sharp incline of CAKUT reports in the literature within the past 2 decades, there remains a plateau in the genetic diagnostic yield that is disproportionate to the accelerated ability to generate robust genome-wide data. Explanations for this observation include (i) diverse inheritance patterns with incomplete penetrance and variable expressivity, (ii) rarity of single-gene drivers such that large sample sizes are required to meet the burden of proof, and (iii) multigene interactions that might produce either intra- (e.g., copy number variants) or inter- (e.g., effects in trans) locus effects. These challenges present an opportunity for the community to implement innovative genetic and molecular avenues to explain the missing heritability and to better elucidate the mechanisms that underscore CAKUT. Here, we review recent multidisciplinary approaches at the intersection of genetics, genomics, in vivo modeling, and in vitro systems toward refining a blueprint for overcoming the diagnostic hurdles that are pervasive in urinary tract malformation cohorts. These approaches will not only benefit clinical management by reducing age at molecular diagnosis and prompting early evaluation for comorbid features but will also serve as a springboard for therapeutic development.

GWAS in Mice Maps Susceptibility to HIV-Associated Nephropathy to the Ssbp2 Locus

BACKGROUND

To gain insight into the pathogenesis of collapsing glomerulopathy, a rare form of FSGS that often arises in the setting of viral infections, we performed a genome-wide association study (GWAS) among inbred mouse strains using a murine model of HIV-1 associated nephropathy (HIVAN).

METHODS

We first generated F1 hybrids between HIV-1 transgenic mice on the FVB/NJ background and 20 inbred laboratory strains. Analysis of histology, BUN, and urinary NGAL demonstrated marked phenotypic variation among the transgenic F1 hybrids, providing strong evidence for host genetic factors in the predisposition to nephropathy. A GWAS in 365 transgenic F1 hybrids generated from these 20 inbred strains was performed.

RESULTS

We identified a genome-wide significant locus on chromosome 13-C3 and multiple additional suggestive loci. Crossannotation of the Chr. 13 locus, including single-cell transcriptomic analysis of wildtype and HIV-1 transgenic mouse kidneys, nominated Ssbp2 as the most likely candidate gene. Ssbp2 is highly expressed in podocytes, encodes a transcriptional cofactor that interacts with LDB1 and LMX1B, which are both previously implicated in FSGS. Consistent with these data, older Ssbp2 null mice spontaneously develop glomerulosclerosis, tubular casts, interstitial fibrosis, and inflammation, similar to the HIVAN mouse model.

CONCLUSIONS

These findings demonstrate the utility of GWAS in mice to uncover host genetic factors for rare kidney traits and suggest Ssbp2 as susceptibility gene for HIVAN, potentially acting via the LDB1-LMX1B transcriptional network.

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.