Population data improves variant interpretation in autosomal dominant polycystic kidney disease

Exome Sequencing of a Clinical Population for Autosomal Dominant Polycystic Kidney Disease

IMPORTANCE

Most studies of autosomal dominant polycystic kidney disease (ADPKD) genetics have used kidney specialty cohorts, focusing on PKD1 and PKD2. These can lead to biased estimates of population prevalence of ADPKD-associated gene variants and their phenotypic expression.

OBJECTIVE

To determine the prevalence of ADPKD and contributions of PKD1, PKD2, and other genes related to cystic kidney disease in a large, unselected cohort.

DESIGN, SETTING, AND PARTICIPANTS

This retrospective observational study used an unselected health system-based cohort in central and northeast Pennsylvania with exome sequencing (enrolled from 2004 to 2020) and electronic health record data (up to October 2021). The genotype-first approach included the entire cohort and the phenotype-first approach focused on patients with ADPKD diagnosis codes, confirmed by chart and imaging review.

EXPOSURES

Loss-of-function (LOF) variants in PKD1, PKD2, and other genes associated with cystic kidney disease (ie, ALG8, ALG9, DNAJB11, GANAB, HNF1B, IFT140, SEC61B, PKHD1, PRKCSH, SEC63); likely pathogenic missense variants in PKD1 and PKD2.

MAIN OUTCOMES AND MEASURES

Genotype-first analysis: ADPKD diagnosis code (Q61.2, Q61.3, 753.13, 753.12); phenotype-first analysis: presence of a rare variant in PKD1, PKD2, or other genes associated with cystic kidney disease.

RESULTS

Of 174 172 patients (median age, 60 years; 60.6% female; 93% of European ancestry), 303 patients had ADPKD diagnosis codes, including 235 with sufficient chart review data for confirmation. In addition to PKD1 and PKD2, LOF variants in IFT140, GANAB, and HNF1B were associated with ADPKD diagnosis after correction for multiple comparisons. Among patients with LOF variants in PKD1, 66 of 68 (97%) had ADPKD; 43 of 43 patients (100%) with LOF variants in PKD2 had ADPKD. In contrast, only 24 of 77 patients (31.2%) with a PKD1 missense variant previously classified as "likely pathogenic" had ADPKD, suggesting misclassification or variable penetrance. Among patients with ADPKD diagnosis confirmed by chart review, 180 of 235 (76.6%) had a potential genetic cause, with the majority being rare variants in PKD1 (127 patients) or PKD2 (34 patients); 19 of 235 (8.1%) had variants in other genes associated with cystic kidney disease. Of these 235 patients with confirmed ADPKD, 150 (63.8%) had a family history of ADPKD. The yield for a genetic determinant of ADPKD was higher for those with a family history of ADPKD compared with those without family history (91.3% [137/150] vs 50.6% [43/85]; difference, 40.7% [95% CI, 29.2%-52.3%]; P < .001). Previously unreported PKD1, PKD2, and GANAB variants were identified with pedigree data suggesting pathogenicity, and several PKD1 missense variants previously reported as likely pathogenic appeared to be benign.

CONCLUSIONS AND RELEVANCE

This study demonstrates substantial genetic and phenotypic variability in ADPKD among patients within a regional health system in the US.

An intermediate-effect size variant in UMOD confers risk for chronic kidney disease

The kidney-specific gene UMOD encodes for uromodulin, the most abundant protein excreted in normal urine. Rare large-effect variants in UMOD cause autosomal dominant tubulointerstitial kidney disease (ADTKD), while common low-impact variants strongly associate with kidney function and the risk of chronic kidney disease (CKD) in the general population. It is unknown whether intermediate-effect variants in UMOD contribute to CKD. Here, candidate intermediate-effect UMOD variants were identified using large-population and ADTKD cohorts. Biological and phenotypical effects were investigated using cell models, in silico simulations, patient samples, and international databases and biobanks. Eight UMOD missense variants reported in ADTKD are present in the Genome Aggregation Database (gnomAD), with minor allele frequency (MAF) ranging from 10-5 to 10-3. Among them, the missense variant p.Thr62Pro is detected in ∼1/1,000 individuals of European ancestry, shows incomplete penetrance but a high genetic load in familial clusters of CKD, and is associated with kidney failure in the 100,000 Genomes Project (odds ratio [OR] = 3.99 [1.84 to 8.98]) and the UK Biobank (OR = 4.12 [1.32 to 12.85). Compared with canonical ADTKD mutations, the p.Thr62Pro carriers displayed reduced disease severity, with slower progression of CKD and an intermediate reduction of urinary uromodulin levels, in line with an intermediate trafficking defect in vitro and modest induction of endoplasmic reticulum (ER) stress. Identification of an intermediate-effect UMOD variant completes the spectrum of UMOD-associated kidney diseases and provides insights into the mechanisms of ADTKD and the genetic architecture of CKD.

Effects of Suramin on Polycystic Kidney Disease in a Mouse Model of Polycystin-1 Deficiency

The aberrant activation of the purinergic signaling pathway has been shown to promote cyst growth and fluid secretion in autosomal dominant polycystic kidney disease (ADPKD). Suramin is an anti-parasitic drug that has strong anti-purinergic properties. Whether suramin could have a therapeutic effect on ADPKD has not been fully investigated. We examined the effect of suramin on cyst progression in a Pkd1 microRNAs transgenic mouse model that presented stable Pkd1 knockdown and moderate disease progression. The Pkd1-deficient mice were treated with suramin (60 mg/kg) by intraperitoneal injection twice a week from postnatal days 35 to 90. Kidney-to-body weight ratios, cyst indices, and blood urea nitrogen (BUN) levels were measured. Cell proliferation and macrophage infiltration were determined by immunohistochemistry. The suramin-treated group had significantly lower renal cyst densities, cell proliferation, and macrophage infiltration compared with saline-treated controls. Suramin significantly inhibited ERK phosphorylation and the expression of Il1b, Il6, Nlrp3, Tgfb, Fn1, P2rx7, and P2ry2 mRNAs in the kidneys. However, BUN levels remained high despite the reduction in cyst growth. Furthermore, plasma cystatin C and neutrophil gelatinase-associated lipocalin (NGAL) levels were significantly higher in the suramin-treated group compared with the control group. Periodic acid-Schiff staining revealed degenerative changes and epithelial cell vacuolation in the non-cystic renal tubules, which indicated phospholipidosis following suramin treatment. These results suggest that suramin may reduce renal cyst growth and inflammation, but the associated tubular cell injuries could limit its therapeutic potential. Other purinergic receptor antagonists with less nephrotoxicity may deserve further investigation for the treatment of ADPKD.

Genetic analysis of a family presenting with coexisting cerebral cavernous malformations and polycystic kidney disease

Hereditary cerebral cavernous malformations (CCMs) are characterized by clustered dilated capillary-like vessels in the brain. Autosomal dominant polycystic kidney disease (PKD) is characterized by renal cysts and extra-renal abnormalities. We report a Taiwanese family in which the index case exhibited coexisting phenotypes of both CCMs and PKD. The index case was a 55-year-old woman with known PKD who developed an intracerebral hemorrhage (ICH) in the right medulla. Neuroimaging revealed numerous microbleeds in the bilateral cerebrum and cerebellum. Radiological CCMs were suspected given the absence of other imaging markers of small vessel disease. A comprehensive panel of 183 cerebral vascular malformation genes were investigated through genome sequencing. A novel CCM2 frameshift variant (c.607_608delCT, p.Leu203Valfs∗53) causing a pathogenic premature stop codon, and a known PKD2 nonsense variant (c.2407C > T, p.Arg803∗), were found. Segregation analysis revealed that four siblings were affected by either isolated aforementioned PKD2 or CCM2 variant. Notably, radiological CCMs were exclusively found in siblings who had this CCM2 variant, and bilateral internal carotid artery aneurysms were restricted to one sibling who had the PKD2 variant but not the CCM2 variant. Our study expands the genetic spectrum of CCM2 and demonstrates unambiguous cosegregation of CCM2 and PKD2 variants with their respective phenotypes.

Analysis of rare genetic variation underlying cardiometabolic diseases and traits among 200,000 individuals in the UK Biobank

Cardiometabolic diseases are the leading cause of death worldwide. Despite a known genetic component, our understanding of these diseases remains incomplete. Here, we analyzed the contribution of rare variants to 57 diseases and 26 cardiometabolic traits, using data from 200,337 UK Biobank participants with whole-exome sequencing. We identified 57 gene-based associations, with broad replication of novel signals in Geisinger MyCode. There was a striking risk associated with mutations in known Mendelian disease genes, including MYBPC3, LDLR, GCK, PKD1 and TTN. Many genes showed independent convergence of rare and common variant evidence, including an association between GIGYF1 and type 2 diabetes. We identified several large effect associations for height and 18 unique genes associated with blood lipid or glucose levels. Finally, we found that between 1.0% and 2.4% of participants carried rare potentially pathogenic variants for cardiometabolic disorders. These findings may facilitate studies aimed at therapeutics and screening of these common disorders.

Genomic diagnostics in polycystic kidney disease: an assessment of real-world use of whole-genome sequencing

Autosomal Dominant Polycystic Kidney Disease (ADPKD) is common, with a prevalence of 1/1000 and predominantly caused by disease-causing variants in PKD1 or PKD2. Clinical diagnosis is usually by age-dependent imaging criteria, which is challenging in patients with atypical clinical features, without family history, or younger age. However, there is increasing need for definitive diagnosis of ADPKD with new treatments available. Sequencing is complicated by six pseudogenes that share 97% homology to PKD1 and by recently identified phenocopy genes. Whole-genome sequencing can definitively diagnose ADPKD, but requires validation for clinical use. We initially performed a validation study, in which 42 ADPKD patients underwent sequencing of PKD1 and PKD2 by both whole-genome and Sanger sequencing, using a blinded, cross-over method. Whole-genome sequencing identified all PKD1 and PKD2 germline pathogenic variants in the validation study (sensitivity and specificity 100%). Two mosaic variants outside pipeline thresholds were not detected. We then examined the first 144 samples referred to a clinically-accredited diagnostic laboratory for clinical whole-genome sequencing, with targeted-analysis to a polycystic kidney disease gene-panel. In this unselected, diagnostic cohort (71 males :73 females), the diagnostic rate was 70%, including a diagnostic rate of 81% in patients with typical ADPKD (98% with PKD1/PKD2 variants) and 60% in those with atypical features (56% PKD1/PKD2; 44% PKHD1/HNF1B/GANAB/ DNAJB11/PRKCSH/TSC2). Most patients with atypical disease did not have clinical features that predicted likelihood of a genetic diagnosis. These results suggest clinicians should consider diagnostic genomics as part of their assessment in polycystic kidney disease, particularly in atypical disease.

A phenome-wide association study of 26 mendelian genes reveals phenotypic expressivity of common and rare variants within the general population

The clinical evaluation of a genetic syndrome relies upon recognition of a characteristic pattern of signs or symptoms to guide targeted genetic testing for confirmation of the diagnosis. However, individuals displaying a single phenotype of a complex syndrome may not meet criteria for clinical diagnosis or genetic testing. Here, we present a phenome-wide association study (PheWAS) approach to systematically explore the phenotypic expressivity of common and rare alleles in genes associated with four well-described syndromic diseases (Alagille (AS), Marfan (MS), DiGeorge (DS), and Noonan (NS) syndromes) in the general population.

Using human phenotype ontology (HPO) terms, we systematically mapped 60 phenotypes related to AS, MS, DS and NS in 337,198 unrelated white British from the UK Biobank (UKBB) based on their hospital admission records, self-administrated questionnaires, and physiological measurements. We performed logistic regression adjusting for age, sex, and the first 5 genetic principal components, for each phenotype and each variant in the target genes (JAG1, NOTCH2 FBN1, PTPN1 and RAS-opathy genes, and genes in the 22q11.2 locus) and performed a gene burden test.

Overall, we observed multiple phenotype-genotype correlations, such as the association between variation in JAG1, FBN1, PTPN11 and SOS2 with diastolic and systolic blood pressure; and pleiotropy among multiple variants in syndromic genes. For example, rs11066309 in PTPN11 was significantly associated with a lower body mass index, an increased risk of hypothyroidism and a smaller size for gestational age, all in concordance with NS-related phenotypes. Similarly, rs589668 in FBN1 was associated with an increase in body height and blood pressure, and a reduced body fat percentage as observed in Marfan syndrome.

Our findings suggest that the spectrum of associations of common and rare variants in genes involved in syndromic diseases can be extended to individual phenotypes within the general population.

Standard medical evaluation of genetic syndromes relies upon recognizing a characteristic pattern of signs or symptoms to guide targeted genetic testing for confirmation of the diagnosis. This may lead to missing diagnoses in patients with silent or a low expressed form of the syndrome. Here we take advantage of a rich electronic health record, various phenotypic measurements, and genetic information in 337,198 unrelated white British from the UKBB, to study the relation between single syndromic disease phenotypes and genes related to syndromic disease. We show multiple phenotype-genotype associations in concordance with phenotypes variations found in syndromic diseases. For example, we show that a commonly found variant in FBN1 was associated with high standing/sitting height ratio and reduced body fat percentage as observed in individuals with Marfan syndrome. Our findings suggest that common and rare alleles in syndromic disease genes are causative of individual component phenotypes present in a general population; further research is needed to characterize the pleiotropic effect of alleles in syndromic genes in persons without the syndromic disease.

Identification of a pathogenic mutation in a Chinese pedigree with polycystic kidney disease

Polycystic kidney disease (PKD) is a life-threatening inherited disease with a morbidity of 1:500–1,000 worldwide. Numerous progressively enlarging cysts are observed in the bilateral kidneys of patients with PKD, inducing structural damage and loss of kidney function. The present study analyzed one family with PKD. Whole exome sequencing of the proband was performed to detect the pathogenic gene present in the family. Candidate gene segments for lineal consanguinity in the family were amplified by nest polymerase chain reaction, followed by Sanger sequencing. One novel duplication variant (NM_001009944.2:c.9359dupA:p.Y3120_E3121delinsX) and one missense mutation (c.G9022A:p.V3008M) were detected in PKD1. Additionally, the pathogenic substitutions in PKD1 published from the dataset were analyzed. Following analysis and confirmation, the duplication variant NM_001009944.2:c.9359dupA:p.Y3120_E3121delinsX in PKD1, within the polycystin-1, lipoxygenase, α-toxin domain, was considered to be the pathogenic factor in the examined family with autosomal dominant PKD. Additionally, based on the analysis of 4,805 pathogenic substitutions in PKD1 within various regions, the presence of the missense mutation in the N-terminal domain of polycystin-1 may present high pathogenicity in ADPKD.