A Neanderthal haplotype introgressed into the human genome confers protection against membranous nephropathy.

Class 2 HLA and PLA2R1 alleles are exceptionally strong genetic risk factors for membranous nephropathy (MN), leading, through an unknown mechanism, to a targeted autoimmune response. Introgressed archaic haplotypes (introduced from an archaic human genome into the modern human genome) might influence phenotypes through gene dysregulation. Here, we investigated the genomic region surrounding the PLA2R1 gene. We reconstructed the phylogeny of Neanderthal and modern haplotypes in this region and calculated the probability of the observed clustering being the result of introgression or common descent. We imputed variants for the participants in our previous genome-wide association study and we compared the distribution of Neanderthal variants between MN cases and controls. The region associated with the lead MN risk locus in the PLA2R1 gene was confirmed and showed that, within a 507 kb region enriched in introgressed sequence, a stringently defined 105 kb haplotype, intersecting the coding regions for PLA2R1 and ITGB6, is inherited from Neanderthals. Thus, introgressed Neanderthal haplotypes overlapping PLA2R1 are differentially represented in MN cases and controls, with enrichment In controls suggesting a protective effect.

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.

A Founder Mutation in EHD1 Presents with Tubular Proteinuria and Deafness.

BACKGROUND: The endocytic reabsorption of proteins in the proximal tubule requires a complex machinery and defects can lead to tubular proteinuria. The precise mechanisms of endocytosis and processing of receptors and cargo are incompletely understood. EHD1 belongs to a family of proteins presumably involved in the scission of intracellular vesicles and in ciliogenesis. However, the relevance of EHD1 in human tissues, in particular in the kidney, was unknown. METHODS: Genetic techniques were used in patients with tubular proteinuria and deafness to identify the disease-causing gene. Diagnostic and functional studies were performed in patients and disease models to investigate the pathophysiology. RESULTS: We identified six individuals (5-33 years) with proteinuria and a high-frequency hearing deficit associated with the homozygous missense variant c.1192C>T (p.R398W) in EHD1. Proteinuria (0.7-2.1 g/d) consisted predominantly of low molecular weight proteins, reflecting impaired renal proximal tubular endocytosis of filtered proteins. Ehd1 knockout and Ehd1R398W/R398W knockin mice also showed a high-frequency hearing deficit and impaired receptor-mediated endocytosis in proximal tubules, and a zebrafish model showed impaired ability to reabsorb low molecular weight dextran. Interestingly, ciliogenesis appeared unaffected in patients and mouse models. In silico structural analysis predicted a destabilizing effect of the R398W variant and possible inference with nucleotide binding leading to impaired EHD1 oligomerization and membrane remodeling ability. CONCLUSIONS: A homozygous missense variant of EHD1 causes a previously unrecognized autosomal recessive disorder characterized by sensorineural deafness and tubular proteinuria. Recessive EHD1 variants should be considered in individuals with hearing impairment, especially if tubular proteinuria is noted.

IgA Nephropathy Genetic Risk Score to Estimate the Prevalence of IgA Nephropathy in UK Biobank.

BACKGROUND: IgA nephropathy (IgAN) is the commonest glomerulonephritis worldwide. Its prevalence is difficult to estimate, as people with mild disease do not commonly receive a biopsy diagnosis. We aimed to generate an IgA nephropathy genetic risk score (IgAN-GRS) and estimate the proportion of people with hematuria who had IgAN in the UK Biobank (UKBB). METHODS: We calculated an IgAN-GRS using 14 single-nucleotide polymorphisms (SNPs) drawn from the largest European Genome-Wide Association Study (GWAS) and validated the IgAN-GRS in 464 biopsy-proven IgAN European cases from the UK Glomerulonephritis DNA Bank (UKGDB) and in 379,767 Europeans in the UKBB. We used the mean of IgAN-GRS to calculate the proportion of potential IgAN in 14,181 with hematuria and other nonspecific renal phenotypes from 379,767 Europeans in the UKBB. RESULTS: The IgAN-GRS was higher in the IgAN cohort (4.30; 95% confidence interval [95% CI: 4.23-4.38) than in controls (3.98; 3.97-3.98; P < 0.0001). The mean GRS in UKBB participants with hematuria (n = 12,858) was higher (4.04; 4.02-4.06) than UKBB controls (3.98; 3.97-3.98; P < 0.0001) and higher in those with hematuria, hypertension, and microalbuminuria (n = 1323) (4.07; 4.02-4.13) versus (3.98; 3.97-3.98; P = 0.0003). Using the difference in these means, we estimated that IgAN accounted for 19% of noncancer hematuria and 28% with hematuria, hypertension, and microalbuminuria in UKBB. CONCLUSIONS: We used an IgAN-GRS to estimate the prevalence of IgAN contributing to common phenotypes that are not always biopsied. The noninvasive use of polygenic risk in this setting may have further utility to identify likely etiology of nonspecific renal phenotypes in large population cohorts.

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.

Genetic Identification of Two Novel Loci Associated with Steroid-Sensitive Nephrotic Syndrome.

BACKGROUND: Steroid-sensitive nephrotic syndrome (SSNS), the most common form of nephrotic syndrome in childhood, is considered an autoimmune disease with an established classic HLA association. However, the precise etiology of the disease is unclear. In other autoimmune diseases, the identification of loci outside the classic HLA region by genome-wide association studies (GWAS) has provided critical insights into disease pathogenesis. Previously conducted GWAS of SSNS have not identified non-HLA loci achieving genome-wide significance. METHODS: In an attempt to identify additional loci associated with SSNS, we conducted a GWAS of a large cohort of European ancestry comprising 422 ethnically homogeneous pediatric patients and 5642 ethnically matched controls. RESULTS: The GWAS found three loci that achieved genome-wide significance, which explain approximately 14% of the genetic risk for SSNS. It confirmed the previously reported association with the HLA-DR/DQ region (lead single-nucleotide polymorphism [SNP] rs9273542, P=1.59×10-43; odds ratio [OR], 3.39; 95% confidence interval [95% CI], 2.86 to 4.03) and identified two additional loci outside the HLA region on chromosomes 4q13.3 and 6q22.1. The latter contains the calcium homeostasis modulator family member 6 gene CALHM6 (previously called FAM26F). CALHM6 is implicated in immune response modulation; the lead SNP (rs2637678, P=1.27×10-17; OR, 0.51; 95% CI, 0.44 to 0.60) exhibits strong expression quantitative trait loci effects, the risk allele being associated with lower lymphocytic expression of CALHM6. CONCLUSIONS: Because CALHM6 is implicated in regulating the immune response to infection, this may provide an explanation for the typical triggering of SSNS onset by infections. Our results suggest that a genetically conferred risk of immune dysregulation may be a key component in the pathogenesis of SSNS.

Glycine Amidinotransferase (GATM), Renal Fanconi Syndrome, and Kidney Failure.

Background For many patients with kidney failure, the cause and underlying defect remain unknown. Here, we describe a novel mechanism of a genetic order characterized by renal Fanconi syndrome and kidney failure.Methods We clinically and genetically characterized members of five families with autosomal dominant renal Fanconi syndrome and kidney failure. We performed genome-wide linkage analysis, sequencing, and expression studies in kidney biopsy specimens and renal cells along with knockout mouse studies and evaluations of mitochondrial morphology and function. Structural studies examined the effects of recognized mutations.Results The renal disease in these patients resulted from monoallelic mutations in the gene encoding glycine amidinotransferase (GATM), a renal proximal tubular enzyme in the creatine biosynthetic pathway that is otherwise associated with a recessive disorder of creatine deficiency. In silico analysis showed that the particular GATM mutations, identified in 28 members of the five families, create an additional interaction interface within the GATM protein and likely cause the linear aggregation of GATM observed in patient biopsy specimens and cultured proximal tubule cells. GATM aggregates-containing mitochondria were elongated and associated with increased ROS production, activation of the NLRP3 inflammasome, enhanced expression of the profibrotic cytokine IL-18, and increased cell death.Conclusions In this novel genetic disorder, fully penetrant heterozygous missense mutations in GATM trigger intramitochondrial fibrillary deposition of GATM and lead to elongated and abnormal mitochondria. We speculate that this renal proximal tubular mitochondrial pathology initiates a response from the inflammasome, with subsequent development of kidney fibrosis.

OVAS: an open-source variant analysis suite with inheritance modelling.

BACKGROUND: The advent of modern high-throughput genetics continually broadens the gap between the rising volume of sequencing data, and the tools required to process them. The need to pinpoint a small subset of functionally important variants has now shifted towards identifying the critical differences between normal variants and disease-causing ones. The ever-increasing reliance on cloud-based services for sequence analysis and the non-transparent methods they utilize has prompted the need for more in-situ services that can provide a safer and more accessible environment to process patient data, especially in circumstances where continuous internet usage is limited. RESULTS: To address these issues, we herein propose our standalone Open-source Variant Analysis Sequencing (OVAS) pipeline; consisting of three key stages of processing that pertain to the separate modes of annotation, filtering, and interpretation. Core annotation performs variant-mapping to gene-isoforms at the exon/intron level, append functional data pertaining the type of variant mutation, and determine hetero/homozygosity. An extensive inheritance-modelling module in conjunction with 11 other filtering components can be used in sequence ranging from single quality control to multi-file penetrance model specifics such as X-linked recessive or mosaicism. Depending on the type of interpretation required, additional annotation is performed to identify organ specificity through gene expression and protein domains. In the course of this paper we analysed an autosomal recessive case study. OVAS made effective use of the filtering modules to recapitulate the results of the study by identifying the prescribed compound-heterozygous disease pattern from exome-capture sequence input samples. CONCLUSION: OVAS is an offline open-source modular-driven analysis environment designed to annotate and extract useful variants from Variant Call Format (VCF) files, and process them under an inheritance context through a top-down filtering schema of swappable modules, run entirely off a live bootable medium and accessed locally through a web-browser.

HaploForge: a comprehensive pedigree drawing and haplotype visualization web application.

MOTIVATION: Haplotype reconstruction is an important tool for understanding the aetiology of human disease. Haplotyping infers the most likely phase of observed genotypes conditional on constraints imposed by the genotypes of other pedigree members. The results of haplotype reconstruction, when visualized appropriately, show which alleles are identical by descent despite the presence of untyped individuals. When used in concert with linkage analysis, haplotyping can help delineate a locus of interest and provide a succinct explanation for the transmission of the trait locus. Unfortunately, the design choices made by existing haplotype visualization programs do not scale to large numbers of markers. Indeed, following haplotypes from generation to generation requires excessive scrolling back and forth. In addition, the most widely used program for haplotype visualization produces inconsistent recombination artefacts for the X chromosome. RESULTS: To resolve these issues, we developed HaploForge, a novel web application for haplotype visualization and pedigree drawing. HaploForge takes advantage of HTML5 to be fast, portable and avoid the need for local installation. It can accurately visualize autosomal and X-linked haplotypes from both outbred and consanguineous pedigrees. Haplotypes are coloured based on identity by descent using a novel A* search algorithm and we provide a flexible viewing mode to aid visual inspection. HaploForge can currently process haplotype reconstruction output from Allegro, GeneHunter, Merlin and Simwalk. AVAILABILITY AND IMPLEMENTATION: HaploForge is licensed under GPLv3 and is hosted and maintained via GitHub. https://github.com/mtekman/haploforge. CONTACT: r.kleta@ucl.ac.uk. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Polycystic Kidney Disease with Hyperinsulinemic Hypoglycemia Caused by a Promoter Mutation in Phosphomannomutase 2.

Hyperinsulinemic hypoglycemia (HI) and congenital polycystic kidney disease (PKD) are rare, genetically heterogeneous disorders. The co-occurrence of these disorders (HIPKD) in 17 children from 11 unrelated families suggested an unrecognized genetic disorder. Whole-genome linkage analysis in five informative families identified a single significant locus on chromosome 16p13.2 (logarithm of odds score 6.5). Sequencing of the coding regions of all linked genes failed to identify biallelic mutations. Instead, we found in all patients a promoter mutation (c.-167G>T) in the phosphomannomutase 2 gene (PMM2), either homozygous or in trans with PMM2 coding mutations. PMM2 encodes a key enzyme in N-glycosylation. Abnormal glycosylation has been associated with PKD, and we found that deglycosylation in cultured pancreatic ß cells altered insulin secretion. Recessive coding mutations in PMM2 cause congenital disorder of glycosylation type 1a (CDG1A), a devastating multisystem disorder with prominent neurologic involvement. Yet our patients did not exhibit the typical clinical or diagnostic features of CDG1A. In vitro, the PMM2 promoter mutation associated with decreased transcriptional activity in patient kidney cells and impaired binding of the transcription factor ZNF143. In silico analysis suggested an important role of ZNF143 for the formation of a chromatin loop including PMM2 We propose that the PMM2 promoter mutation alters tissue-specific chromatin loop formation, with consequent organ-specific deficiency of PMM2 leading to the restricted phenotype of HIPKD. Our findings extend the spectrum of genetic causes for both HI and PKD and provide insights into gene regulation and PMM2 pleiotropy.

Generation and validation of a zebrafish model of EAST (epilepsy, ataxia, sensorineural deafness and tubulopathy) syndrome.

Recessive mutations in KCNJ10, which encodes an inwardly rectifying potassium channel, were recently identified as the cause of EAST syndrome, a severe and disabling multi-organ disorder consisting of epilepsy, ataxia, sensorineural deafness and tubulopathy that becomes clinically apparent with seizures in infancy. A Kcnj10 knockout mouse shows postnatal mortality and is therefore not suitable for drug discovery. Because zebrafish are ideal for in vivo screening for potential therapeutics, we tested whether kcnj10 knockdown in zebrafish would fill this need. We cloned zebrafish kcnj10 and demonstrated that its function is equivalent to that of human KCNJ10. We next injected splice- and translation-blocking kcnj10 antisense morpholino oligonucleotides and reproduced the cardinal symptoms of EAST syndrome - ataxia, epilepsy and renal tubular defects. Several of these phenotypes could be assayed in an automated manner. We could rescue the morphant phenotype with complementary RNA (cRNA) encoding human wild-type KCNJ10, but not with cRNA encoding a KCNJ10 mutation identified in individuals with EAST syndrome. Our results suggest that zebrafish will be a valuable tool to screen for compounds that are potentially therapeutic for EAST syndrome or its individual symptoms. Knockdown of kcnj10 represents the first zebrafish model of a salt-losing tubulopathy, which has relevance for blood pressure control.