Mendelian Disease Genes Research Articles

A library of induced pluripotent stem cells from clinically well-characterized, diverse healthy human individuals.

SummaryA library of well-characterized human induced pluripotent stem cell (hiPSC) lines from clinically healthy human subjects could serve as a useful resource of normal controls for in vitro human development, disease modeling, genotype-phenotype association studies, and drug response evaluation. We report generation and extensive characterization of a gender-balanced, racially/ethnically diverse library of hiPSC lines from 40 clinically healthy human individuals who range in age from 22 to 61 years. The hiPSCs match the karyotype and short tandem repeat identities of their parental fibroblasts, and have a transcription profile characteristic of pluripotent stem cells. We provide whole-genome sequencing data for one hiPSC clone from each individual, genomic ancestry determination, and analysis of mendelian disease genes and risks. We document similar transcriptomic profiles, single-cell RNA-sequencing-derived cell clusters, and physiology of cardiomyocytes differentiated from multiple independent hiPSC lines. This extensive characterization makes this hiPSC library a valuable resource for many studies on human biology.

Decreased recent adaptation at human mendelian disease genes as a possible consequence of interference between advantageous and deleterious variants

Advances in genome sequencing have improved our understanding of the genetic basis of human diseases, and thousands of human genes have been associated with different diseases. Recent genomic adaptation at disease genes has not been well characterized. Here, we compare the rate of strong recent adaptation in the form of selective sweeps between mendelian, non-infectious disease genes and non-disease genes across distinct human populations from the 1000 Genomes Project. We find that mendelian disease genes have experienced far less selective sweeps compared to non-disease genes especially in Africa. Investigating further the possible causes of the sweep deficit at disease genes, we find that this deficit is very strong at disease genes with both low recombination rates and with high numbers of associated disease variants, but is almost non-existent at disease genes with higher recombination rates or lower numbers of associated disease variants. Because segregating recessive deleterious variants have the ability to interfere with adaptive ones, these observations strongly suggest that adaptation has been slowed down by the presence of interfering recessive deleterious variants at disease genes. These results suggest that disease genes suffer from a transient inability to adapt as fast as the rest of the genome.

The Path to Conserved Extended Haplotypes: Megabase-Length Haplotypes at High Population Frequency.

This minireview describes the history of the conceptual development of conserved extended haplotypes (CEHs): megabase-length haplotypes that exist at high (≥0.5%) population frequency. My career began in internal medicine, shifted to pediatrics, and clinical practice changed to research. My research interest was initially in hematology: on plasma proteins, their metabolism, synthesis, and function. This narrowed to a focus on proteins of the human complement system, their role in immunity and their genetics, beginning with polymorphism and deficiency of C3. My group identified genetic polymorphisms and/or inherited deficiencies of C2, C4, C6, and C8. After defining glycine-rich beta glycoprotein as factor B (Bf) in the properdin system, we found that the genes for Bf (CFB), C2, C4A, and C4B were inherited as a single haplotypic unit which we named the “complotype.” Complotypes are located within the major histocompatibility complex (MHC) between HLA-B and HLA-DRB1 and are designated (in arbitrary order) by their CFB, C2, C4A, and C4B types. Pedigree analysis revealed long stretches (several megabases) of apparently fixed DNA within the MHC that we referred to as “extended haplotypes” (later as “CEHs”). About 10 to 12 common CEHs constitute at least 25 – 30% of MHC haplotypes among European Caucasian populations. These CEHs contain virtually all the most common markers of MHC-associated diseases. In the case of type 1 diabetes, we have proposed a purely genetic and epigenetic model (with a small number of Mendelian recessive disease genes) that explains all the puzzling features of the disease, including its rising incidence.

Data-driven modelling of mutational hotspots and in silico predictors in hypertrophic cardiomyopathy

BackgroundAlthough rare missense variants in Mendelian disease genes often cluster in specific regions of proteins, it is unclear how to consider this when evaluating the pathogenicity of a gene or...

Mutation spectrum of NOD2 reveals recessive inheritance as a main driver of Early Onset Crohn\u2019s Disease

Inflammatory bowel disease (IBD), clinically defined as Crohn’s disease (CD), ulcerative colitis (UC), or IBD-unclassified, results in chronic inflammation of the gastrointestinal tract in genetically susceptible hosts. Pediatric onset IBD represents ≥ 25% of all IBD diagnoses and often presents with intestinal stricturing, perianal disease, and failed response to conventional treatments. NOD2 was the first and is the most replicated locus associated with adult IBD, to date. However, its role in pediatric onset IBD is not well understood. We performed whole-exome sequencing on a cohort of 1,183 patients with pediatric onset IBD (ages 0–18.5 years). We identified 92 probands with biallelic rare and low frequency NOD2 variants accounting for approximately 8% of our cohort, suggesting a Mendelian inheritance pattern of disease. Additionally, we investigated the contribution of recessive inheritance of NOD2 alleles in adult IBD patients from a large clinical population cohort. We found that recessive inheritance of NOD2 variants explains ~ 7% of cases in this adult IBD cohort, including ~ 10% of CD cases, confirming the observations from our pediatric IBD cohort. Exploration of EHR data showed that several of these adult IBD patients obtained their initial IBD diagnosis before 18 years of age, consistent with early onset disease. While it has been previously reported that carriers of more than one NOD2 risk alleles have increased susceptibility to Crohn’s Disease (CD), our data formally demonstrate that recessive inheritance of NOD2 alleles is a mechanistic driver of early onset IBD, specifically CD, likely due to loss of NOD2 protein function. Collectively, our findings show that recessive inheritance of rare and low frequency deleterious NOD2 variants account for 7–10% of CD cases and implicate NOD2 as a Mendelian disease gene for early onset Crohn’s Disease.

Prioritizing non-coding regions based on human genomic constraint and sequence context with deep learning

Elucidating functionality in non-coding regions is a key challenge in human genomics. It has been shown that intolerance to variation of coding and proximal non-coding sequence is a strong predictor of human disease relevance. Here, we integrate intolerance to variation, functional genomic annotations and primary genomic sequence to build JARVIS: a comprehensive deep learning model to prioritize non-coding regions, outperforming other human lineage-specific scores. Despite being agnostic to evolutionary conservation, JARVIS performs comparably or outperforms conservation-based scores in classifying pathogenic single-nucleotide and structural variants. In constructing JARVIS, we introduce the genome-wide residual variation intolerance score (gwRVIS), applying a sliding-window approach to whole genome sequencing data from 62,784 individuals. gwRVIS distinguishes Mendelian disease genes from more tolerant CCDS regions and highlights ultra-conserved non-coding elements as the most intolerant regions in the human genome. Both JARVIS and gwRVIS capture previously inaccessible human-lineage constraint information and will enhance our understanding of the non-coding genome.

MendelVar: gene prioritization at GWAS loci using phenotypic enrichment of Mendelian disease genes.

MotivationGene prioritization at human GWAS loci is challenging due to linkage-disequilibrium and long-range gene regulatory mechanisms. However, identifying the causal gene is crucial to enable identification of potential drug targets and better understanding of molecular mechanisms. Mapping GWAS traits to known phenotypically relevant Mendelian disease genes near a locus is a promising approach to gene prioritization.ResultsWe present MendelVar, a comprehensive tool that integrates knowledge from four databases on Mendelian disease genes with enrichment testing for a range of associated functional annotations such as Human Phenotype Ontology, Disease Ontology and variants from ClinVar. This open web-based platform enables users to strengthen the case for causal importance of phenotypically matched candidate genes at GWAS loci. We demonstrate the use of MendelVar in post-GWAS gene annotation for type 1 diabetes, type 2 diabetes, blood lipids and atopic dermatitis.Availability and implementationMendelVar is freely available at https://mendelvar.mrcieu.ac.ukSupplementary information Supplementary data are available at Bioinformatics online.

Whole-Exome Sequencing Uncovers Novel Causative Variants and Additional Findings in Three Patients Affected by Glycogen Storage Disease Type VI and Fanconi-Bickel Syndrome.

Glycogen storage diseases (GSDs) are the heterogeneous group of disorders caused by mutations in at least 30 different genes. Different types of GSDs, especially liver GSDs, take overlapping symptoms and can be clinically indistinguishable. This survey evaluated the use of whole-exome sequencing (WES) for the genetic analysis of the liver GSD-suspected patients in three unrelated families. An in-house filtering pipeline was used to assess rare pathogenic variants in GSD-associated genes, autosomal recessive/mendelian disorder genes (carrier status for genetic counseling subjects), and the ACMG’s list of 59 actionable genes. For the interpretation of the causative variants and the incidental/secondary findings, ACMG guidelines were applied. Additionally, we have explored PharmGKB class IA/IB pharmacogenetic variants. The segregation analysis was performed using Sanger sequencing for the novel causative variants. Bioinformatics analysis of the exome data in three individuals revealed three novel homozygous causative variants in the GSD-associated genes. The first variant, c.298_307delATGATCAACC in PYGL gene has related to HERS disease (GSD VI). Both variants of c.1043dupT and c.613-1G > C in SLC2A2 gene have been associated with Fanconi-Bickel syndrome (GSDXI). Eight pathogenic/likely pathogenic medical actionable findings in Mendelian disease genes and 10 pharmacogenetic variants with underlying drug response phenotypes have been identified. No known/expected pathogenic variants were detected in the ACMG’s list of 59 actionable genes. The logical filtering steps can help in finding other medical actionable secondary/incidental findings as well as effectively identifying the causative variants in heterogeneous conditions such as GSDs. Three novel variants related to GSD genes recognized in liver GSD-suspected patients with early infantile and childhood-age onset.

Evolutionary Properties of Biological Pathways in Health and Disease

Evolutionary properties of genes are related to their function in health and disease. Previous work has studied these properties separately, focusing on specific gene categories. This prevents a systems-level understanding of their impact on biological processes and involvement in disease. Here, we develop a random forest classifier integrating nine properties into a cumulative score. This score divides 25 functional pathways encompassing 10,334 genes into three groups. High-scoring pathways which accomplish basic cell functions, are enriched in tumour suppressors and core essential genes with little tolerance towards loss-of-function alterations. Intermediate-scoring pathways which contribute to metabolism, development and immune system, include context essential genes, oncogenes and dominant Mendelian disease genes. Low-scoring pathways perform organ-specific functions and include drivers of somatic clonal growth and recessive Mendelian disease genes. Our integrated analysis rationalises the evolution of human genes within broader functional groups and enables to predict their role in germline and somatic diseases.

Long-read trio sequencing of individuals with unsolved intellectual disability

Long-read sequencing (LRS) has the potential to comprehensively identify all medically relevant genome variation, including variation commonly missed by short-read sequencing (SRS) approaches. To determine this potential, we performed LRS around 15×–40× genome coverage using the Pacific Biosciences Sequel I System for five trios. The respective probands were diagnosed with intellectual disability (ID) whose etiology remained unresolved after SRS exomes and genomes. Systematic assessment of LRS coverage showed that ~35 Mb of the human reference genome was only accessible by LRS and not SRS. Genome-wide structural variant (SV) calling yielded on average 28,292 SV calls per individual, totaling 12.9 Mb of sequence. Trio-based analyses which allowed to study segregation, showed concordance for up to 95% of these SV calls across the genome, and 80% of the LRS SV calls were not identified by SRS. De novo mutation analysis did not identify any de novo SVs, confirming that these are rare events. Because of high sequence coverage, we were also able to call single nucleotide substitutions. On average, we identified 3 million substitutions per genome, with a Mendelian inheritance concordance of up to 97%. Of these, ~100,000 were located in the ~35 Mb of the genome that was only captured by LRS. Moreover, these variants affected the coding sequence of 64 genes, including 32 known Mendelian disease genes. Our data show the potential added value of LRS compared to SRS for identifying medically relevant genome variation.

The Impact of Next-Generation Sequencing on the Diagnosis, Treatment, and Prevention of Hereditary Neuromuscular Disorders

The impact of high-throughput sequencing in genetic neuromuscular disorders cannot be overstated. The ability to rapidly and affordably sequence multiple genes simultaneously has enabled a second golden age of Mendelian disease gene discovery, with flow-on impacts for rapid genetic diagnosis, evidence-based treatment, tailored therapy development, carrier-screening, and prevention of disease recurrence in families. However, there are likely many more neuromuscular disease genes and mechanisms to be discovered. Many patients and families remain without a molecular diagnosis following targeted panel sequencing, clinical exome sequencing, or even genome sequencing. Here we review how massively parallel, or next-generation, sequencing has changed the field of genetic neuromuscular disorders, and anticipate future benefits of recent technological innovations such as RNA-seq implementation and detection of tandem repeat expansions from short-read sequencing.

NEW GENES AND DISEASES / NGS & RELATED TECHNIQUES: P.333 Progressive proximal muscle weakness, myopathic muscle biopsy and a novel candidate gene

More than 200 new Mendelian disease gene discoveries have been published every year since the introduction of next-generation sequencing technologies. To date, almost 600 genes are associated with neuromuscular disorders. We present the phenotype of an Estonian girl harboring variants in a novel muscle disease candidate gene. The female patient was first noted to have exercise intolerance at age 3 years. Two years later, she presented with abnormal gait, toe walking, frequent falls, lumbar hyperlordosis, hyperextention of the knees, and valgus feet.

A polyclonal allelic expression assay for detecting regulatory effects of transcript variants

We present an assay to experimentally test the regulatory effects of genetic variants within transcripts using CRISPR/Cas9 followed by targeted sequencing. We applied the assay to 32 premature stop-gained variants across the genome and in two Mendelian disease genes, 33 putative causal variants of eQTLs, and 62 control variants in HEK293T cells, replicating a subset of variants in HeLa cells. We detected significant effects in the expected direction (in 60% of variants), demonstrating the ability of the assay to capture regulatory effects of eQTL variants and nonsense-mediated decay triggered by premature stop-gained variants. The results suggest a utility for validating transcript-level effects of genetic variants.

Beyond medically actionable results: an analytical pipeline for decreasing the burden of returning all clinically significant secondary findings.

Genomic sequencing advances have increased the potential to identify secondary findings (SFs). Current guidelines recommend the analysis of 59 medically actionable genes; however, patient preferences indicate interest in learning a broader group of SFs. We aimed to develop an analytical pipeline for the efficient analysis and return of all clinically significant SFs. We developed a pipeline consisting of comprehensive gene lists for five categories of SFs and filtration parameters for prioritization of variants in each category. We applied the pipeline to 42 exomes to assess feasibility and efficiency. Comprehensive lists of clinically significant SF genes were curated for each category: (1) 90 medically actionable genes and 28 pharmacogenomic variants; (2) 17 common disease risk variants; (3) 3166 Mendelian disease genes, (4) 7 early onset neurodegenerative disorder genes; (5) 688 carrier status results. Analysis of 42 exomes using our pipeline resulted in a significant decrease (> 98%) in variants compared to the raw analysis (13,036.56 ± 59.72 raw variants/exome vs 161.32 ± 7.68 filtered variants/exome), and aided in time and costs savings for the overall analysis process. Our pipeline represents a critical step in overcoming the analytic challenge associated with returning all clinically relevant SFs to allow for its routine implementation in clinical practice.

Fitting a naturally scaled point system to the ACMG/AMP variant classification guidelines.

Recently, we demonstrated that the qualitative American College of Medical Genetics and Genomics/Association for Medical Pathology (ACMG/AMP) guidelines for evaluation of Mendelian disease gene variants are fundamentally compatible with a quantitative Bayesian formulation. Here, we show that the underlying ACMG/AMP "strength of evidence categories" can be abstracted into a point system. These points are proportional to Log(odds), are additive, and produce a system that recapitulates the Bayesian formulation of the ACMG/AMP guidelines. The strengths of this system are its simplicity and that the connection between point values and odds of pathogenicity allows empirical calibration of the strength of evidence for individual data types. Weaknesses include that a narrow range of prior probabilities is locked in and that the Bayesian nature of the system is inapparent. We conclude that a points-based system has the practical attribute of user-friendliness and can be useful so long as the underlying Bayesian principles are acknowledged.

Gene-wise Variant Burden and Genomic Characterization of Nearly Every Gene

Aim: Current gene-level prioritization methods aim to provide information for further prioritization of 'disease-causing' mutations. Since, they are inherently biased toward disease genes, methods specific to pharmacogenetic (PGx) genes are required. Methods: We proposed a gene-wise variant burden (GVB) method that integrates in silico deleteriousness scores of the multitude of variants of a given gene at a personal-genome level. Results: GVB in its simplest form outperformed the two state-of-the-art methods with regard to predicting pharmacogenes and complex disease genes but not for rare Mendelian disease genes. GVB* adjusted by paralog counts robustly performed well in most of the pharmacogenetic subcategories. Seven molecular genetic features well characterized the unique genomic properties of PGx, complex, and Mendelian disease genes. Conclusion: Altogether, GVB is an individual-specific genescore, especially advantageous for PGx studies.

Electronic health record phenotypes associated with genetically regulated expression of CFTR and application to cystic fibrosis

PurposeThe increasing use of electronic health records (EHRs) and biobanks offers unique opportunities to study Mendelian diseases. We described a novel approach to summarize clinical manifestations from patient EHRs into phenotypic evidence for cystic fibrosis (CF) with potential to alert unrecognized patients of the disease. MethodsWe estimated genetically predicted expression (GReX) of cystic fibrosis transmembrane conductance regulator (CFTR) and tested for association with clinical diagnoses in the Vanderbilt University biobank (N=9142 persons of European descent with 71 cases of CF). The top associated EHR phenotypes were assessed in combination as a phenotype risk score (PheRS) for discriminating CF case status in an additional 2.8 million patients from Vanderbilt University Medical Center (VUMC) and 125,305 adult patients including 25,314 CF cases from MarketScan, an independent external cohort. ResultsGReX of CFTR was associated with EHR phenotypes consistent with CF. PheRS constructed using the EHR phenotypes and weights discovered by the genetic associations improved discriminative power for CF over the initially proposed PheRS in both VUMC and MarketScan. ConclusionOur study demonstrates the power of EHRs for clinical description of CF and the benefits of using a genetics-informed weighing scheme in construction of a phenotype risk score. This research may find broad applications for phenomic studies of Mendelian disease genes.

Hemizygous mutations in L1CAM in two unrelated male probands with childhood onset psychosis

To identify genes underlying childhood onset psychosis. Patients with onset of psychosis at age 13 or younger were identified from clinics across England, and they and their parents were exome sequenced and analysed for possible highly penetrant genetic contributors. We report two male childhood onset psychosis patients of different ancestries carrying hemizygous very rare possibly damaging missense variants (p.Arg846His and p.Pro145Ser) in the L1CAM gene. L1CAM is an X-linked Mendelian disease gene in which both missense and loss of function variants are associated with syndromic forms of intellectual disability and developmental disorder. This is the first study reporting a possible extension of the phenotype of L1CAM variant carriers to childhood onset psychosis. The family history and presence of other significant rare genetic variants in the patients suggest that there may be genetic interactions modulating the presentation.

De novo variants in SIAH1, encoding an E3 ubiquitin ligase, are associated with developmental delay, hypotonia and dysmorphic features

BackgroundUbiquitination has a central role in numerous biological processes, including cell development, stress responses and ageing. Perturbed ubiquitination has been implicated in human diseases ranging from cancer to neurodegenerative diseases....

Genomic Data Sharing for Novel Mendelian Disease Gene Discovery: The Matchmaker Exchange.

In the last decade, exome and/or genome sequencing has become a common test in the diagnosis of individuals with features of a rare Mendelian disorder. Despite its success, this test leaves the majority of tested individuals undiagnosed. This review describes the Matchmaker Exchange (MME), a federated network established to facilitate the solving of undiagnosed rare-disease cases through data sharing. MME supports genomic matchmaking, the act of connecting two or more parties looking for cases with similar phenotypes and variants in the same candidate genes. An application programming interface currently connects six matchmaker nodes-the Database of Chromosomal Imbalance and Phenotype in Humans Using Ensembl Resources (DECIPHER), GeneMatcher, PhenomeCentral, seqr, MyGene2, and the Initiative on Rare and Undiagnosed Diseases (IRUD) Exchange-resulting in a collective data set spanning more than 150,000 cases from more than 11,000 contributors in 88 countries. Here, we describe the successes and challenges of MME, its individual matchmaking nodes, plans for growing the network, and considerations for future directions.