Rare variants in the sodium-dependent phosphate transporter gene SLC34A3 explain missing heritability of urinary stone disease.

Urinary stone disease (USD) is a major health burden affecting over 10% of the United Kingdom population. While stone disease is associated with lifestyle, genetic factors also strongly contribute. Common genetic variants at multiple loci from genome-wide association studies account for 5% of the estimated 45% heritability of the disorder. Here, we investigated the extent to which rare genetic variation contributes to the unexplained heritability of USD. Among participants of the United Kingdom 100,000-genome project, 374 unrelated individuals were identified and assigned diagnostic codes indicative of USD. Whole genome gene-based rare variant testing and polygenic risk scoring against a control population of 24,930 ancestry-matched controls was performed. We observed (and replicated in an independent dataset) exome-wide significant enrichment of monoallelic rare, predicted damaging variants in the SLC34A3 gene for a sodium-dependent phosphate transporter that were present in 5% cases compared with 1.6% of controls. This gene was previously associated with autosomal recessive disease. The effect on USD risk of having a qualifying SLC34A3 variant was greater than that of a standard deviation increase in polygenic risk derived from GWAS. Addition of the rare qualifying variants in SLC34A3 to a linear model including polygenic score increased the liability-adjusted heritability from 5.1% to 14.2% in the discovery cohort. We conclude that rare variants in SLC34A3 represent an important genetic risk factor for USD, with effect size intermediate between the fully penetrant rare variants linked with Mendelian disorders and common variants associated with USD. Thus, our findings explain some of the heritability unexplained by prior common variant genome-wide association studies.

Shared genetic risk across different presentations of gene test-negative idiopathic nephrotic syndrome.

BACKGROUND: Idiop athic nephrotic syndrome (INS) is classified in children according to response to initial corticosteroid therapy into steroid-sensitive (SSNS) and steroid-resistant nephrotic syndrome (SRNS), and in adults according to histology into minimal change disease (MCD) and focal segmental glomerulosclerosis (FSGS). However, there is well-recognised phenotypic overlap between these entities. Genome-wide association studies (GWAS) have shown a strong association between SSNS and variation at HLA, suggesting an underlying immunological basis. We sought to determine whether a risk score generated from genetic variants associated with SSNS could be used to gain insight into the pathophysiology of INS presenting in other ways. METHODS: We developed an SSNS genetic risk score (SSNS-GRS) from the five variants independently associated with childhood SSNS in a previous European GWAS. We quantified SSNS-GRS in independent cohorts of European individuals with childhood SSNS, non-monogenic SRNS, MCD, and FSGS, and contrasted them with SSNS-GRS quantified in individuals with monogenic SRNS, membranous nephropathy (a different immune-mediated disease-causing nephrotic syndrome), and healthy controls. RESULTS: The SSNS-GRS was significantly elevated in cohorts with SSNS, non-monogenic SRNS, MCD, and FSGS compared to healthy participants and those with membranous nephropathy. The SSNS-GRS in all cohorts with non-monogenic INS were also significantly elevated compared to those with monogenic SRNS. CONCLUSIONS: The shared genetic risk factors among patients with different presentations of INS strongly suggests a shared autoimmune pathogenesis when monogenic causes are excluded. Use of the SSNS-GRS, in addition to testing for monogenic causes, may help to classify patients presenting with INS. A higher resolution version of the Graphical abstract is available as Supplementary information.

Diverse ancestry whole-genome sequencing association study identifies TBX5 and PTK7 as susceptibility genes for posterior urethral valves.

Posterior urethral valves (PUV) are the commonest cause of end-stage renal disease in children, but the genetic architecture of this rare disorder remains unknown. We performed a sequencing-based genome-wide association study (seqGWAS) in 132 unrelated male PUV cases and 23,727 controls of diverse ancestry, identifying statistically significant associations with common variants at 12q24.21 (p=7.8 × 10-12; OR 0.4) and rare variants at 6p21.1 (p=2.0 × 10-8; OR 7.2), that were replicated in an independent European cohort of 395 cases and 4151 controls. Fine mapping and functional genomic data mapped these loci to the transcription factor TBX5 and planar cell polarity gene PTK7, respectively, the encoded proteins of which were detected in the developing urinary tract of human embryos. We also observed enrichment of rare structural variation intersecting with candidate cis-regulatory elements, particularly inversions predicted to affect chromatin looping (p=3.1 × 10-5). These findings represent the first robust genetic associations of PUV, providing novel insights into the underlying biology of this poorly understood disorder and demonstrate how a diverse ancestry seqGWAS can be used for disease locus discovery in a rare disease.

GWAS meta-analysis of intrahepatic cholestasis of pregnancy implicates multiple hepatic genes and regulatory elements.

Intrahepatic cholestasis of pregnancy (ICP) is a pregnancy-specific liver disorder affecting 0.5-2% of pregnancies. The majority of cases present in the third trimester with pruritus, elevated serum bile acids and abnormal serum liver tests. ICP is associated with an increased risk of adverse outcomes, including spontaneous preterm birth and stillbirth. Whilst rare mutations affecting hepatobiliary transporters contribute to the aetiology of ICP, the role of common genetic variation in ICP has not been systematically characterised to date. Here, we perform genome-wide association studies (GWAS) and meta-analyses for ICP across three studies including 1138 cases and 153,642 controls. Eleven loci achieve genome-wide significance and have been further investigated and fine-mapped using functional genomics approaches. Our results pinpoint common sequence variation in liver-enriched genes and liver-specific cis-regulatory elements as contributing mechanisms to ICP susceptibility.

Genotype-phenotype correlations for COL4A3-COL4A5 variants resulting in Gly substitutions in Alport syndrome.

Alport syndrome is the commonest inherited kidney disease and nearly half the pathogenic variants in the COL4A3-COL4A5 genes that cause Alport syndrome result in Gly substitutions. This study examined the molecular characteristics of Gly substitutions that determine the severity of clinical features. Pathogenic COL4A5 variants affecting Gly in the Leiden Open Variation Database in males with X-linked Alport syndrome were correlated with age at kidney failure (n = 157) and hearing loss diagnosis (n = 80). Heterozygous pathogenic COL4A3 and COL4A4 variants affecting Gly (n = 304) in autosomal dominant Alport syndrome were correlated with the risk of haematuria in the UK 100,000 Genomes Project. Gly substitutions were stratified by exon location (1 to 20 or 21 to carboxyl terminus), being adjacent to a non-collagenous region (interruption or terminus), and the degree of instability caused by the replacement residue. Pathogenic COL4A5 variants that resulted in a Gly substitution with a highly destabilising residue reduced the median age at kidney failure by 7 years (p = 0.002), and age at hearing loss diagnosis by 21 years (p = 0.004). Substitutions adjacent to a non-collagenous region delayed kidney failure by 19 years (p = 0.014). Heterozygous pathogenic COL4A3 and COL4A4 variants that resulted in a Gly substitution with a highly destabilising residue (Arg, Val, Glu, Asp, Trp) were associated with an increased risk of haematuria (p = 0.018), and those adjacent to a non-collagenous region were associated with a reduced risk (p = 0.046). Exon location had no effect. In addition, COL4A5 variants adjacent to non-collagenous regions were over-represented in the normal population in gnomAD (p < 0.001). The nature of the substitution and of nearby residues determine the risk of haematuria, early onset kidney failure and hearing loss for Gly substitutions in X-linked and autosomal dominant Alport syndrome.

Gitelman-Like Syndrome Caused by Pathogenic Variants in mtDNA.

BACKGROUND: Gitelman syndrome is the most frequent hereditary salt-losing tubulopathy characterized by hypokalemic alkalosis and hypomagnesemia. Gitelman syndrome is caused by biallelic pathogenic variants in SLC12A3, encoding the Na+-Cl- cotransporter (NCC) expressed in the distal convoluted tubule. Pathogenic variants of CLCNKB, HNF1B, FXYD2, or KCNJ10 may result in the same renal phenotype of Gitelman syndrome, as they can lead to reduced NCC activity. For approximately 10 percent of patients with a Gitelman syndrome phenotype, the genotype is unknown. METHODS: We identified mitochondrial DNA (mtDNA) variants in three families with Gitelman-like electrolyte abnormalities, then investigated 156 families for variants in MT-TI and MT-TF, which encode the transfer RNAs for phenylalanine and isoleucine. Mitochondrial respiratory chain function was assessed in patient fibroblasts. Mitochondrial dysfunction was induced in NCC-expressing HEK293 cells to assess the effect on thiazide-sensitive 22Na+ transport. RESULTS: Genetic investigations revealed four mtDNA variants in 13 families: m.591C>T (n=7), m.616T>C (n=1), m.643A>G (n=1) (all in MT-TF), and m.4291T>C (n=4, in MT-TI). Variants were near homoplasmic in affected individuals. All variants were classified as pathogenic, except for m.643A>G, which was classified as a variant of uncertain significance. Importantly, affected members of six families with an MT-TF variant additionally suffered from progressive chronic kidney disease. Dysfunction of oxidative phosphorylation complex IV and reduced maximal mitochondrial respiratory capacity were found in patient fibroblasts. In vitro pharmacological inhibition of complex IV, mimicking the effect of the mtDNA variants, inhibited NCC phosphorylation and NCC-mediated sodium uptake. CONCLUSION: Pathogenic mtDNA variants in MT-TF and MT-TI can cause a Gitelman-like syndrome. Genetic investigation of mtDNA should be considered in patients with unexplained Gitelman syndrome-like tubulopathies.

Prevalence Estimates of Predicted Pathogenic COL4A3-COL4A5 Variants in a Population Sequencing Database and Their Implications for Alport Syndrome.

BACKGROUND: The reported prevalence of Alport syndrome varies from one in 5000 to one in 53,000 individuals. This study estimated the frequencies of predicted pathogenic COL4A3-COL4A5 variants in sequencing databases of populations without known kidney disease. METHODS: Predicted pathogenic variants were identified using filtering steps based on the ACMG/AMP criteria, which considered collagen IV α3-α5 position 1 Gly to be critical domains. The population frequencies of predicted pathogenic COL4A3-COL4A5 variants were then determined per mean number of sequenced alleles. Population frequencies for compound heterozygous and digenic combinations were calculated from the results for heterozygous variants. RESULTS: COL4A3-COL4A5 variants resulting in position 1 Gly substitutions were confirmed to be associated with hematuria (for each, P<0.001). Predicted pathogenic COL4A5 variants were found in at least one in 2320 individuals. p.(Gly624Asp) represented nearly half (16 of 33, 48%) of the variants in Europeans. Most COL4A5 variants (54 of 59, 92%) had a biochemical feature that potentially mitigated the clinical effect. The predicted pathogenic heterozygous COL4A3 and COL4A4 variants affected one in 106 of the population, consistent with the finding of thin basement membrane nephropathy in normal donor kidney biopsy specimens. Predicted pathogenic compound heterozygous variants occurred in one in 88,866 individuals, and digenic variants in at least one in 44,793. CONCLUSIONS: The population frequencies for Alport syndrome are suggested by the frequencies of predicted pathogenic COL4A3-COL4A5 variants, but must be adjusted for the disease penetrance of individual variants and for the likelihood of already diagnosed disease and non-Gly substitutions. Disease penetrance may depend on other genetic and environmental factors.

Genetics of renovascular hypertension in children.

OBJECTIVE: In most cases of renovascular hypertension in children, the cause is unclear. The aim of this study was to investigate genetic variation as a factor in the development of renovascular hypertension in children. METHODS: In a cohort of 37 unrelated children from a single tertiary referral center, exome sequencing was performed. We assessed variants in recognized and suspected disease genes and searched for novel ones with a gene-based variant-burden analysis. RESULTS: In the majority of patients, exome sequencing could not identify causative variants. We found a pathogenic variant in a recognized associated disease gene in five patients (three pathogenic variants in NF1, one in ELN and a deletion of chromosome 7q11.23, consistent with Williams syndrome). In two other patients, (likely) pathogenic variants were found in putative renovascular hypertension genes (SMAD6 and GLA), with clinical implications for both. Ten additional patients carried variants of uncertain significance (VUS) in known (n = 4) or putative (n = 6) renovascular hypertension disease genes. Rare variant burden analysis yielded no further candidate genes. CONCLUSION: Genetic contributors, such as germline mutations in NF1, ELN, 7q11.23del were present in only 5 out of 37 (14%) children with renovascular hypertension. Twelve other children (32%) had potentially causal variants identified, including a pathogenic variant in SMAD6; a vasculopathy gene hitherto unknown to link with renovascular hypertension. Most importantly, our data show that exome sequencing can rarely identify the cause of renovascular hypertension in nonsyndromic children. We suggest that nongenetic factors or somatic genetic variation will play a more important role.

Large-Scale Whole-Genome Sequencing Reveals the Genetic Architecture of Primary Membranoproliferative GN and C3 Glomerulopathy.

BACKGROUND: Primary membranoproliferative GN, including complement 3 (C3) glomerulopathy, is a rare, untreatable kidney disease characterized by glomerular complement deposition. Complement gene mutations can cause familial C3 glomerulopathy, and studies have reported rare variants in complement genes in nonfamilial primary membranoproliferative GN. METHODS: We analyzed whole-genome sequence data from 165 primary membranoproliferative GN cases and 10,250 individuals without the condition (controls) as part of the National Institutes of Health Research BioResource-Rare Diseases Study. We examined copy number, rare, and common variants. RESULTS: Our analysis included 146 primary membranoproliferative GN cases and 6442 controls who were unrelated and of European ancestry. We observed no significant enrichment of rare variants in candidate genes (genes encoding components of the complement alternative pathway and other genes associated with the related disease atypical hemolytic uremic syndrome; 6.8% in cases versus 5.9% in controls) or exome-wide. However, a significant common variant locus was identified at 6p21.32 (rs35406322) (P=3.29×10-8; odds ratio [OR], 1.93; 95% confidence interval [95% CI], 1.53 to 2.44), overlapping the HLA locus. Imputation of HLA types mapped this signal to a haplotype incorporating DQA1*05:01, DQB1*02:01, and DRB1*03:01 (P=1.21×10-8; OR, 2.19; 95% CI, 1.66 to 2.89). This finding was replicated by analysis of HLA serotypes in 338 individuals with membranoproliferative GN and 15,614 individuals with nonimmune renal failure. CONCLUSIONS: We found that HLA type, but not rare complement gene variation, is associated with primary membranoproliferative GN. These findings challenge the paradigm of complement gene mutations typically causing primary membranoproliferative GN and implicate an underlying autoimmune mechanism in most cases.

Cascade Fumarate Hydratase mutation screening allows early detection of kidney tumour: a case report.

BACKGROUND: Fumarate hydratase (FH) deficiency is a rare autosomal recessive disorder which results in a major defect in cellular metabolism. It presents in infancy with progressive encephalopathy, hypotonia, seizures and failure to thrive and is often fatal in childhood. It is caused by mutations in the FH gene (1q42.1) that result in deficiency of the citric acid cycle enzyme fumarate hydratase, resulting in accumulation of fumaric acid. Heterozygous germline mutations in the FH gene predispose to an aggressive autosomal dominant inherited early-onset kidney cancer syndrome: hereditary leiomyomatosis and renal cell cancer (HLRCC). CASE PRESENTATION: Cascade FH mutation screening enabled the early diagnosis of a renal tumour in an asymptomatic parent of a child with fumarate hydratase deficiency, resulting in timely and possibly life-saving treatment. CONCLUSION: While the theoretical risk of kidney cancer in parents of children with recessive fumarate hydratase deficiency is well recognized, to our knowledge this is the first report of a kidney tumour being detected in a parent by screening performed for this indication. This underscores the importance of offering lifelong kidney surveillance to such parents and other heterozygous relatives of children born with fumarate hydratase deficiency.