Identifying modifier genes of monogenic disease: strategies and difficulties

How do stochastic processes and genetic threshold effects explain incomplete penetrance and inform causal disease mechanisms?

Incomplete penetrance is the rule rather than the exception in Mendelian disease. In syndromic monogenic disorders, phenotypic variability can be viewed as the combination of incomplete penetrance for each of multiple independent clinical features. Within genetically identical individuals, such as isogenic model organisms, stochastic variation at molecular and cellular levels is the primary cause of incomplete penetrance according to a genetic threshold model. By defining specific probability distributions of causal biological readouts and genetic liability values, stochasticity and incomplete penetrance provide information about threshold values in biological systems. Ascertainment of threshold values has been achieved by simultaneous scoring of relatively simple phenotypes and quantitation of molecular readouts at the level of single cells. However, this is much more challenging for complex morphological phenotypes using experimental and reductionist approaches alone, where cause and effect are separated temporally and across multiple biological modes and scales. Here I consider how causal inference, which integrates observational data with high confidence causal models, might be used to quantify the relative contribution of different sources of stochastic variation to phenotypic diversity. Collectively, these approaches could inform disease mechanisms, improve predictions of clinical outcomes and prioritize gene therapy targets across modes and scales of gene function. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Sterol O-Acyltransferase 1 (SOAT1): A Genetic Modifier of Niemann-Pick Disease, Type C1

Niemann-Pick disease type C1 (NPC1) is a lysosomal disorder due to impaired intracellular cholesterol transport out of the endolysosomal compartment.. Marked heterogeneity has been observed in individuals with the same NPC1 genotype, thus suggesting a significant effect of modifier genes. Prior work demonstrated that decreased SOAT1 activity decreased disease severity in an NPC1 mouse model. Thus, we hypothesized that a polymorphism associated with decreased SOAT1 expression might influence the NPC1 phenotype. Phenotyping and genomic sequencing of 117 individuals with NPC1 was performed as part of a Natural History trial. Phenotyping included determination of disease severity and disease burden. Significant clinical heterogeneity is present in individuals homozygous for the NPC1I1061T variant and in siblings. Analysis of the SOAT1 polymorphism, rs1044925 (A>C), showed a significant association of the C-allele with earlier age of neurological onset. The C-allele may be associated with a higher Annualized Severity Index Score as well as increased frequency of liver disease and seizures. A polymorphism associated with decreased expression of SOAT1 appears to be a genetic modifier of the NPC1 phenotype. This finding is consistent with prior data showing decreased phenotypic severity in Npc1-/-:Soat1-/- mice and supports efforts to investigate the potential of SOAT1 inhibitors as a potential therapy for NPC1.

Stochastic gene expression and environmental stressors trigger variable somite segmentation phenotypes

Mutations of several genes cause incomplete penetrance and variable expressivity of phenotypes, which are usually attributed to modifier genes or gene-environment interactions. Here, we show stochastic gene expression underlies the variability of somite segmentation defects in embryos mutant for segmentation clock genes her1 or her7. Phenotypic strength is further augmented by low temperature and hypoxia. By performing live imaging of the segmentation clock reporters, we further show that groups of cells with higher oscillation amplitudes successfully form somites while those with lower amplitudes fail to do so. In unfavorable environments, the number of cycles with high amplitude oscillations and the number of successful segmentations proportionally decrease. These results suggest that individual oscillation cycles stochastically fail to pass a threshold amplitude, resulting in segmentation defects in mutants. Our quantitative methodology is adaptable to investigate variable phenotypes of mutant genes in different tissues.

Global analysis of suppressor mutations that rescue human genetic defects

BACKGROUND

Genetic suppression occurs when the deleterious effects of a primary "query" mutation, such as a disease-causing mutation, are rescued by a suppressor mutation elsewhere in the genome.

METHODS

To capture existing knowledge on suppression relationships between human genes, we examined 2,400 published papers for potential interactions identified through either genetic modification of cultured human cells or through association studies in patients.

RESULTS

The resulting network encompassed 476 unique suppression interactions covering a wide spectrum of diseases and biological functions. The interactions frequently linked genes that operate in the same biological process. Suppressors were strongly enriched for genes with a role in stress response or signaling, suggesting that deleterious mutations can often be buffered by modulating signaling cascades or immune responses. Suppressor mutations tended to be deleterious when they occurred in absence of the query mutation, in apparent contrast with their protective role in the presence of the query. We formulated and quantified mechanisms of genetic suppression that could explain 71% of interactions and provided mechanistic insight into disease pathology. Finally, we used these observations to predict suppressor genes in the human genome.

CONCLUSIONS

The global suppression network allowed us to define principles of genetic suppression that were conserved across diseases, model systems, and species. The emerging frequency of suppression interactions among human genes and range of underlying mechanisms, together with the prevalence of suppression in model organisms, suggest that compensatory mutations may exist for most genetic diseases.

Hypertrophic Cardiomyopathy versus Storage Diseases with Myocardial Involvement

One of the main causes of heart failure is cardiomyopathies. Among them, the most common is hypertrophic cardiomyopathy (HCM), characterized by thickening of the left ventricular muscle. This article focuses on HCM and other cardiomyopathies with myocardial hypertrophy, including Fabry disease, Pompe disease, and Danon disease. The genetics and pathogenesis of these diseases are described, as well as current and experimental treatment options, such as pharmacological intervention and the potential of gene therapies. Although genetic approaches are promising and have the potential to become the best treatments for these diseases, further research is needed to evaluate their efficacy and safety. This article describes current knowledge and advances in the treatment of the aforementioned cardiomyopathies.

A Role for Genetic Modifiers in Tubulointerstitial Kidney Diseases

With the increased availability of genomic sequencing technologies, the molecular bases for kidney diseases such as nephronophthisis and mitochondrially inherited and autosomal-dominant tubulointerstitial kidney diseases (ADTKD) has become increasingly apparent. These tubulointerstitial kidney diseases (TKD) are monogenic diseases of the tubulointerstitium and result in interstitial fibrosis and tubular atrophy (IF/TA). However, monogenic inheritance alone does not adequately explain the highly variable onset of kidney failure and extra-renal manifestations. Phenotypes vary considerably between individuals harbouring the same pathogenic variant in the same putative monogenic gene, even within families sharing common environmental factors. While the extreme end of the disease spectrum may have dramatic syndromic manifestations typically diagnosed in childhood, many patients present a more subtle phenotype with little to differentiate them from many other common forms of non-proteinuric chronic kidney disease (CKD). This review summarises the expanding repertoire of genes underpinning TKD and their known phenotypic manifestations. Furthermore, we collate the growing evidence for a role of modifier genes and discuss the extent to which these data bridge the historical gap between apparently rare monogenic TKD and polygenic non-proteinuric CKD (excluding polycystic kidney disease).

Genetic modifiers of synucleinopathies-lessons from experimental models

α-Synuclein is a pleiotropic protein underlying a group of progressive neurodegenerative diseases, including Parkinson's disease and dementia with Lewy bodies. Together, these are known as synucleinopathies. Like all neurological diseases, understanding of disease mechanisms is hampered by the lack of access to biopsy tissues, precluding a real-time view of disease progression in the human body. This has driven researchers to devise various experimental models ranging from yeast to flies to human brain organoids, aiming to recapitulate aspects of synucleinopathies. Studies of these models have uncovered numerous genetic modifiers of α-synuclein, most of which are evolutionarily conserved. This review discusses what we have learned about disease mechanisms from these modifiers, and ways in which the study of modifiers have supported ongoing efforts to engineer disease-modifying interventions for synucleinopathies.

Assessing the digenic model in rare disorders using population sequencing data

An important fraction of patients with rare disorders remains with no clear genetic diagnostic, even after whole-exome or whole-genome sequencing, posing a difficulty in giving adequate treatment and genetic counseling. The analysis of genomic data in rare disorders mostly considers the presence of single gene variants in coding regions that follow a concrete monogenic mode of inheritance. A digenic inheritance, with variants in two functionally-related genes in the same individual, is a plausible alternative that might explain the genetic basis of the disease in some cases. In this case, digenic disease combinations should be absent or underrepresented in healthy individuals. We develop a framework to evaluate the significance of digenic combinations and test its statistical power in different scenarios. We suggest that this approach will be relevant with the advent of new sequencing efforts including hundreds of thousands of samples.

Investigation on the Effects of Modifying Genes on the Spinal Muscular Atrophy Phenotype

Introduction  Spinal muscular atrophy (SMA) is an autosomal recessive neuromuscular disorder caused by the degeneration of motor neurons, muscle weakness, and atrophy that leads to infant's death. The duplication of exon 7/8 in the SMN2 gene reduces the clinical severity of disease, and it is defined as modifying effect. In this study, we aim to investigate the expression of modifying genes related to the prognosis of SMA like PLS3 , PFN2 , ZPR1 , CORO1C , GTF2H2 , NRN1 , SERF1A , NCALD , NAIP , and TIA1.

Methods  Seventeen patients, who came to Trakya University, Faculty of Medicine, Medical Genetics Department, with a preliminary diagnosis of SMA disease, and eight healthy controls were included in this study after multiplex ligation-dependent probe amplification analysis. Gene expression levels were determined by real-time reverse transcription polymerase chain reaction and delta–delta CT method by the isolation of RNA from peripheral blood of patients and controls.

Results   SERF1A and NAIP genes compared between A group and B + C + D groups, and A group of healthy controls, showed statistically significant differences ( p  = 0.037, p  = 0.001).

Discussion   PLS3, NAIP , and NRN1 gene expressions related to SMA disease have been reported before in the literature. In our study, the expression levels of SERF1A , GTF2H2 , NCALD , ZPR1 , TIA1 , PFN2 , and CORO1C genes have been studied for the first time in SMA patients.

Prediction of Regulatory SNPs in Putative Minor Genes of the Neuro-Cardiovascular Variant in Fabry Reveals Insights into Autophagy/Apoptosis and Fibrosis

Even though a mutation in monogenic diseases leads to a “classic” manifestation, many disorders exhibit great clinical variability that could be due to modifying genes also called minor genes. Fabry disease (FD) is an X-linked inborn error resulting from the deficient or absent activity of alpha-galactosidase A (α-GAL) enzyme, that leads to deposits of globotriaosylceramide. With our proprietary software SNPclinic v.1.0, we analyzed 110 single nucleotide polymorphisms (SNPs) in the proximal promoter of 14 genes that could modify the FD phenotype FD. We found seven regulatory-SNP (rSNPs) in three genes (IL10, TGFB1 and EDN1) in five cell lines relevant to FD (Cardiac myocytes and fibroblasts, Astrocytes-cerebellar, endothelial cells and T helper cells 1-TH1). Each SNP was confirmed as a true rSNP in public eQTL databases, and additional software suggested the prediction of variants. The two proposed rSNPs in IL10, could explain components for the regulation of active B cells that influence the fibrosis process. The three predicted rSNPs in TGFB1, could act in apoptosis-autophagy regulation. The two putative rSNPs in EDN1, putatively regulate chronic inflammation. The seven rSNPs described here could act to modulate Fabry’s clinical phenotype so we propose that IL10, TGFB1 and EDN1 be considered minor genes in FD.

The False Dawn of Polygenic Risk Scores for Human Disease Prediction

Polygenic risk scores (PRSs) are being constructed for many diseases and are presented today as a promising avenue in the field of human genetics. These scores aim at predicting the risk of developing a disease by leveraging the many genome-wide association studies (GWAS) conducted during the two last decades. Important investments are being made to improve score estimates by increasing GWAS sample sizes, by developing more sophisticated methods, and by proposing different corrections for potential biases. PRSs have entered the market with direct-to-consumer companies proposing to compute them from saliva samples and even recently to help parents select the healthiest embryos. In this paper, we recall how PRSs arose and question the credit they are given by revisiting underlying assumptions in light of the history of human genetics and by comparing them with estimated breeding values (EBVs) used for selection in livestock.

Exome Sequencing Identifies Genetic Variants Associated with Extreme Manifestations of the Cardiovascular Phenotype in Marfan Syndrome

Marfan Syndrome (MFS) is an autosomal dominant condition caused by variants in the fibrillin-1 (FBN1) gene. Cardinal features of MFS include ectopia lentis (EL), musculoskeletal features and aortic root aneurysm and dissection. Although dissection of the ascending aorta is the main cause of mortality in MFS, the clinical course differs considerably in age of onset and severity, even among individuals who share the same causative variant, suggesting the existence of additional genetic variants that modify the severity of the cardiovascular phenotype in MFS. We recruited MFS patients and classified them into severe (n = 8) or mild aortic phenotype (n = 14) according to age of presentation of the first aorta-related incident. We used Exome Sequencing to identify the genetic variants associated with the severity of aortic manifestations and we performed linkage analysis where suitable. We found five genes associated with severe aortic phenotype and three genes that could be protective for this phenotype in MFS. These genes regulate components of the extracellular matrix, TGFβ pathway and other signaling pathways that are involved in the maintenance of the ECM or angiogenesis. Further studies will be required to understand the functional effect of these variants and explore novel, personalized risk management and, potentially, therapies for these patients.

Pyridoxine-dependent epilepsy (PDE-ALDH7A1) in adulthood: A Dutch pilot study exploring clinical and patient-reported outcomes

Background

Little is known about pyridoxine-dependent epilepsy due to α-aminoadipic semialdehyde dehydrogenase deficiency (PDE-ALDH7A1) in adulthood, as the genetic basis of the disorder has only been elucidated 15 years ago. This creates a knowledge gap for physicians, pediatric patients and their parents, which was aimed to address in this study using clinical data as well as patient-reported outcome measures (PROMs) for the patient's perspective.

Methods

Dutch, genetically confirmed PDE-ALDH7A1 patients ≥18 years were eligible for inclusion. Clinical data were collected as well as PROMs (PROMIS item banks Anxiety, Depression, Anger, Physical Functioning, Cognitive Functioning, Cognitive Abilities, Ability to Participate and Satisfaction with Social Roles).

Results

Ten out of 11 patients agreed to participate (91% response rate). Seizure control at last follow up (median age 25.2 years, range 17.8–29.8 years) was achieved with pyridoxine monotherapy in 70%, 20% with adjunct common-anti epileptic drugs and 10% did not obtain complete seizure control. Neurologic symptoms were present in all but one patient (90%) and included tremors, noted in 40%. Neuro-imaging abnormalities were present in 80%. Intellectual disability was present in 70%. One patient (10%) attended university, three maintained a job without assistance, five maintained a job with assistance or attended social daycare, and one patient never followed regular education. The cohort scored significantly lower on the PROMIS Cognitive Functioning compared to the general (age-related) population. Distribution of scores was wide on all PROMIS item banks.

Discussion & conclusion

Outcomes of this young adult cohort are heterogeneous and individualized approaches are therefore needed. Long-term seizure control with pyridoxine was achieved for almost all patients. Neurologic symptoms were noted in the majority, including tremors, as well as neuro-imaging abnormalities and intellectual disability, additionally reflected by the PROMIS Cognitive Functioning. PDE-ALDH7A1 patients scored comparable to the general population on all other PROMs, especially regarding Ability to Participate and Satisfaction with Social Roles this may indicate a positive interpretation of their functioning. The aim is to expand this pilot study to larger populations to obtain more solid data, and to advance the use of PROMs to engage patients in research and provide the opportunity for personalized care.

Beyond Mendelian Inheritance: Genetic Buffering and Phenotype Variability

Understanding the way genes work amongst individuals and across generations to shape form and function is a common theme for many genetic studies. The recent advances in genetics, genome engineering and DNA sequencing reinforced the notion that genes are not the only players that determine a phenotype. Due to physiological or pathological fluctuations in gene expression, even genetically identical cells can behave and manifest different phenotypes under the same conditions. Here, we discuss mechanisms that can influence or even disrupt the axis between genotype and phenotype; the role of modifier genes, the general concept of genetic redundancy, genetic compensation, the recently described transcriptional adaptation, environmental stressors, and phenotypic plasticity. We furthermore highlight the usage of induced pluripotent stem cells (iPSCs), the generation of isogenic lines through genome engineering, and sequencing technologies can help extract new genetic and epigenetic mechanisms from what is hitherto considered 'noise'.

OUP accepted manuscript

The mutations in modifier genes may contribute to some inherited diseases including Wilson disease (WD). This study was designed to identify potential modifier genes that contribute to WD. A total of 10 WD patients with single or no heterozygous ATP7B mutations were recruited for whole-exome sequencing (WES). Five hundred and thirteen candidate genes, of which the genetic variants present in at least two patients, were identified. In order to clarify which proteins might be involved in copper transfer or metabolism processes, the isobaric tags for relative and absolute quantitation (iTRAQ) was performed to identify the differentially expressed proteins between normal and CuSO₄-treated cell lines. Thirteen genes/proteins were identified by both WES and iTRAQ, indicating that disease-causing variants of these genes may actually contribute to the aberrant copper ion accumulation. Additionally, the c.86C > T (p.S29L) mutation in the SLC31A2 gene (coding CTR2) has a relative higher frequency in our cohort of WD patients (6/191) than reported (0.0024 in gnomAD database) in our healthy donors (0/109), and CTR2^S29L leads to increased intracellular Cu concentration and Cu-induced apoptosis in cultured cell lines. In conclusion, the WES and iTRAQ approaches successfully identified several disease-causing variants in potential modifier genes that may be involved in the WD phenotype.

Novel Machado-Joseph disease-modifying genes and pathways identified by whole-exome sequencing

Machado-Joseph disease (MJD/SCA3) is a neurodegenerative polyglutamine disorder exhibiting a wide spectrum of phenotypes. The abnormal size of the (CAG)n at ATXN3 explains ~55% of the age at onset variance, suggesting the involvement of other factors, namely genetic modifiers, whose identification remains limited. Our aim was to find novel genetic modifiers, analyse their epistatic effects and identify disease-modifying pathways contributing to MJD variable expressivity. We performed whole-exome sequencing in a discovery sample of four age at onset concordant and four discordant first-degree relative pairs of Azorean patients, to identify candidate variants which genotypes differed for each discordant pair but were shared in each concordant pair. Variants identified by this approach were then tested in an independent multi-origin cohort of 282 MJD patients. Whole-exome sequencing identified 233 candidate variants, from which 82 variants in 53 genes were prioritized for downstream analysis. Eighteen disease-modifying pathways were identified; two of the most enriched pathways were relevant for the nervous system, namely the neuregulin signaling and the agrin interactions at neuromuscular junction. Variants at PARD3, NFKB1, CHD5, ACTG1, CFAP57, DLGAP2, ITGB1, DIDO1 and CERS4 modulate age at onset in MJD, with those identified in CFAP57, ACTG1 and DIDO1 showing consistent effects across cohorts of different geographical origins. Network analyses of the nine novel MJD modifiers highlighted several important molecular interactions, including genes/proteins previously related with MJD pathogenesis, namely between ACTG1/APOE and VCP/ITGB1. We describe novel pathways, modifiers, and their interaction partners, providing a broad molecular portrait of age at onset modulation to be further exploited as new disease-modifying targets for MJD and related diseases.

Systems Genetics Analysis Defines Importance Of TMEM43/LUMA For Cardiac And Metabolic Related Pathways

Broad cellular functions and diseases including muscular dystrophy, arrhythmogenic right ventricular cardiomyopathy (ARVC5) and cancer are associated with transmembrane protein43 (TMEM43/LUMA). The study aimed to investigate biological roles of TMEM43 through genetic regulation, gene pathways and gene networks, candidate interacting genes, and up- or downstream regulators. Cardiac transcriptomes from 40 strains of recombinant inbred BXD mice and two parental strains representing murine genetic reference population (GRP) were applied for genetic correlation, functional enrichment, and coexpression network analysis using systems genetics approach. The results were validated in a newly created knock-in Tmem43-S358L mutation mouse model (Tmem43^S358L) that displayed signs of cardiac dysfunction, resembling ARVC5 phenotype seen in humans. We found high Tmem43 levels among BXDs with broad variability in expression. Expression of Tmem43 highly negatively correlated with heart mass and heart rate among BXDs, whereas levels of Tmem43 highly positively correlated with plasma high-density lipoproteins (HDL). Through finding differentially expressed genes (DEGs) between Tmem43^S358L mutant and wild-type (Tmem43^WT) lines, 18 pathways (out of 42 found in BXDs GRP) that are involved in ARVC, hypertrophic cardiomyopathy, dilated cardiomyopathy, nonalcoholic fatty liver disease, Alzheimer’s disease, Parkinson’s disease, and Huntington’s disease were verified. We further constructed Tmem43-mediated gene network, in which Ctnna1, Adcy6, Gnas, Ndufs6, and Uqcrc2 were significantly altered in Tmem43^S358L mice versus Tmem43^WT controls. Our study defined the importance of Tmem43 for cardiac- and metabolism-related pathways, suggesting that cardiovascular disease-relevant risk factors may also increase risk of metabolic and neurodegenerative diseases via TMEM43-mediated pathways.

Targeted Gene Delivery: Where to Land

Genome-editing technologies have the potential to correct most genetic defects involved in blood disorders. In contrast to mutation-specific editing, targeted gene insertion can correct most of the mutations affecting the same gene with a single therapeutic strategy (gene replacement) or provide novel functions to edited cells (gene addition). Targeting a selected genomic harbor can reduce insertional mutagenesis risk, while enabling the exploitation of endogenous promoters, or selected chromatin contexts, to achieve specific transgene expression levels/patterns and the modulation of disease-modifier genes. In this review, we will discuss targeted gene insertion and the advantages and limitations of different genomic harbors currently under investigation for various gene therapy applications.

Mimicking of splicing-related retinitis pigmentosa mutations in C. elegans allow drug screens and identification of disease modifiers

CRISPR/Cas and the high conservation of the spliceosome components facilitate the mimicking of human pathological mutations in splicing factors of model organisms. The degenerative retinal disease retinitis pigmentosa (RP) is caused by mutations in distinct types of genes, including missense mutations in splicing factors that provoke RP in an autosomal dominant form (s-adRP). Using CRISPR in Caenorhabditis elegans, we generated mutant strains to mimic s-adRP mutations reported in PRPF8 and SNRNP200. Whereas these inherited mutations are present in heterozygosis in patients, C. elegans allows the maintenance of these mutations as homozygotes, which is advantageous for genetic and drug screens. We found that snrp-200(cer23[V676L]) and prp-8(cer14[H2302del]) display pleiotropic phenotypes, including reduced fertility. However, snrp-200(cer24[S1080L]) and prp-8(cer22[R2303G]) are weak alleles suitable for RNAi screens for identifying genetic interactions, which could uncover potential disease modifiers. We screened a collection of RNAi clones for splicing-related genes and identified three splicing factors: isy-1/ISY1, cyn-15/PPWD1 and mog-2/SNRPA1, whose partial inactivation may modify the course of the disease. Interestingly, these three genes act as modifiers of prp-8(cer22) but not of snrp-200(cer24). Finally, a screen of the strong allele prp-8(cer14) with FDA-approved drugs did not identify molecules capable of alleviating the temperature-sensitive sterility. Instead, we detected drugs, such as dequalinium chloride, which exacerbated the phenotype, and therefore, are potentially harmful to s-adRP patients since they may accelerate the progression of the disease.

Compensatory epistasis explored by molecular dynamics simulations

A non-negligible proportion of human pathogenic variants are known to be present as wild type in at least some non-human mammalian species. The standard explanation for this finding is that molecular mechanisms of compensatory epistasis can alleviate the mutations' otherwise pathogenic effects. Examples of compensated variants have been described in the literature but the interacting residue(s) postulated to play a compensatory role have rarely been ascertained. In this study, the examination of five human X-chromosomally encoded proteins (FIX, GLA, HPRT1, NDP and OTC) allowed us to identify several candidate compensated variants. Strong evidence for a compensated/compensatory pair of amino acids in the coagulation FIXa protein (involving residues 270 and 271) was found in a variety of mammalian species. Both amino acid residues are located within the 60-loop, spatially close to the 39-loop that performs a key role in coagulation serine proteases. To understand the nature of the underlying interactions, molecular dynamics simulations were performed. The predicted conformational change in the 39-loop consequent to the Glu270Lys substitution (associated with hemophilia B) appears to impair the protein's interaction with its substrate but, importantly, such steric hindrance is largely mitigated in those proteins that carry the compensatory residue (Pro271) at the neighboring amino acid position.

Genome Editing Human Pluripotent Stem Cells to Model β-Cell Disease and Unmask Novel Genetic Modifiers

A mechanistic understanding of the genetic basis of complex diseases such as diabetes mellitus remain elusive due in large part to the activity of genetic disease modifiers that impact the penetrance and/or presentation of disease phenotypes. In the face of such complexity, rare forms of diabetes that result from single-gene mutations (monogenic diabetes) can be used to model the contribution of individual genetic factors to pancreatic β-cell dysfunction and the breakdown of glucose homeostasis. Here we review the contribution of protein coding and non-protein coding genetic disease modifiers to the pathogenesis of diabetes subtypes, as well as how recent technological advances in the generation, differentiation, and genome editing of human pluripotent stem cells (hPSC) enable the development of cell-based disease models. Finally, we describe a disease modifier discovery platform that utilizes these technologies to identify novel genetic modifiers using induced pluripotent stem cells (iPSC) derived from patients with monogenic diabetes caused by heterozygous mutations.

Embryology, Malformations, and Rare Diseases of the Cochlea

Despite the low overall prevalence of individual rare diseases, cochlear dysfunction leading to hearing loss represents a symptom in a large proportion. The aim of this work was to provide a clear overview of rare cochlear diseases, taking into account the embryonic development of the cochlea and the systematic presentation of the different disorders. Although rapid biotechnological and bioinformatic advances may facilitate the diagnosis of a rare disease, an interdisciplinary exchange is often required to raise the suspicion of a rare disease. It is important to recognize that the phenotype of rare inner ear diseases can vary greatly not only in non-syndromic but also in syndromic hearing disorders. Finally, it becomes clear that the phenotype of the individual rare diseases cannot be determined exclusively by classical genetics even in monogenetic disorders.

Genetics and Genomics Approaches for Diagnosis and Research Into Hereditary Anemias

The hereditary anemias are a relatively heterogeneous set of disorders that can show wide clinical and genetic heterogeneity, which often hampers correct clinical diagnosis. The classical diagnostic workflow for these conditions generally used to start with analysis of the family and personal histories, followed by biochemical and morphological evaluations, and ending with genetic testing. However, the diagnostic framework has changed more recently, and genetic testing is now a suitable approach for differential diagnosis of these patients. There are several approaches to this genetic testing, the choice of which depends on phenotyping, genetic heterogeneity, and gene size. For patients who show complete phenotyping, single-gene testing remains recommended. However, genetic analysis now includes next-generation sequencing, which is generally based on custom-designed targeting panels and whole-exome sequencing. The use of next-generation sequencing also allows the identification of new causative genes, and of polygenic conditions and genetic factors that modify disease severity of hereditary anemias. In the research field, whole-genome sequencing is useful for the identification of non-coding causative mutations, which might account for the disruption of transcriptional factor occupancy sites and cis-regulatory elements. Moreover, advances in high-throughput sequencing techniques have now resulted in the identification of genome-wide profiling of the chromatin structures known as the topologically associating domains. These represent a recurrent disease mechanism that exposes genes to inappropriate regulatory elements, causing errors in gene expression. This review focuses on the challenges of diagnosis and research into hereditary anemias, with indications of both the advantages and disadvantages. Finally, we consider the future perspectives for the use of next-generation sequencing technologies in this era of precision medicine.

The GENDULF algorithm: mining transcriptomics to uncover modifier genes for monogenic diseases

Modifier genes are believed to account for the clinical variability observed in many Mendelian disorders, but their identification remains challenging due to the limited availability of genomics data from large patient cohorts. Here, we present GENDULF (GENetic moDULators identiFication), one of the first methods to facilitate prediction of disease modifiers using healthy and diseased tissue gene expression data. GENDULF is designed for monogenic diseases in which the mechanism is loss of function leading to reduced expression of the mutated gene. When applied to cystic fibrosis, GENDULF successfully identifies multiple, previously established disease modifiers, including EHF, SLC6A14, and CLCA1. It is then utilized in spinal muscular atrophy (SMA) and predicts U2AF1 as a modifier whose low expression correlates with higher SMN2 pre-mRNA exon 7 retention. Indeed, knockdown of U2AF1 in SMA patient-derived cells leads to increased full-length SMN2 transcript and SMN protein expression. Taking advantage of the increasing availability of transcriptomic data, GENDULF is a novel addition to existing strategies for prediction of genetic disease modifiers, providing insights into disease pathogenesis and uncovering novel therapeutic targets.

Multi-Omic Approach to Identify Phenotypic Modifiers Underlying Cerebral Demyelination in X-Linked Adrenoleukodystrophy

X-linked adrenoleukodystrophy (ALD) is a peroxisomal metabolic disorder with a highly complex clinical presentation. ALD is caused by mutations in the ABCD1 gene, and is characterized by the accumulation of very long-chain fatty acids in plasma and tissues. Disease-causing mutations are 'loss of function' mutations, with no prognostic value with respect to the clinical outcome of an individual. All male patients with ALD develop spinal cord disease and a peripheral neuropathy in adulthood, although age of onset is highly variable. However, the lifetime prevalence to develop progressive white matter lesions, termed cerebral ALD (CALD), is only about 60%. Early identification of transition to CALD is critical since it can be halted by allogeneic hematopoietic stem cell therapy only in an early stage. The primary goal of this study is to identify molecular markers which may be prognostic of cerebral demyelination from a simple blood sample, with the hope that blood-based assays can replace the current protocols for diagnosis. We collected six well-characterized brother pairs affected by ALD and discordant for the presence of CALD and performed multi-omic profiling of blood samples including genome, epigenome, transcriptome, metabolome/lipidome, and proteome profiling. In our analysis we identify discordant genomic alleles present across all families as well as differentially abundant molecular features across the omics technologies. The analysis was focused on univariate modeling to discriminate the two phenotypic groups, but was unable to identify statistically significant candidate molecular markers. Our study highlights the issues caused by a large amount of inter-individual variation, and supports the emerging hypothesis that cerebral demyelination is a complex mix of environmental factors and/or heterogeneous genomic alleles. We confirm previous observations about the role of immune response, specifically auto-immunity and the potential role of PFN1 protein overabundance in CALD in a subset of the families. We envision our methodology as well as dataset has utility to the field for reproducing previous or enabling future modifier investigations.

Exome sequencing in genetic disease: recent advances and considerations

Over the past decade, exome sequencing (ES) has allowed significant advancements to the field of disease research. By targeting the protein-coding regions of the genome, ES combines the depth of knowledge on protein-altering variants with high-throughput data generation and ease of analysis. New discoveries continue to be made using ES, and medical science has benefitted both theoretically and clinically from its continued use. In this review, we describe recent advances and successes of ES in disease research. Through selected examples of recent publications, we explore how ES continues to be a valuable tool to find variants that might explain disease etiology or provide insight into the biology underlying the disease. We then discuss shortcomings of ES in terms of variant discoveries made by other sequencing technologies that would be missed because of the scope and techniques of ES. We conclude with a brief outlook on the future of ES, suggesting that although newer and more thorough sequencing methods will soon supplant ES, its results will continue to be useful for disease research.

The BMP-SMAD pathway mediates the impaired hepatic iron metabolism associated with the ERFE-A260S variant

The erythroferrone (ERFE) is the erythroid regulator of hepatic iron metabolism by suppressing the expression of hepcidin. Congenital dyserythropoietic anemia type II (CDAII) is an inherited hyporegenerative anemia due to biallelic mutations in the SEC23B gene. Patients with CDAII exhibit marked clinical variability, even among individuals sharing the same pathogenic variants. The ERFE expression in CDAII is increased and related to abnormal erythropoiesis. We identified a recurrent low-frequency variant, A260S, in the ERFE gene in 12.5% of CDAII patients with a severe phenotype. We demonstrated that the ERFE-A260S variant leads to increased levels of ERFE, with subsequently marked impairment of iron regulation pathways at the hepatic level. Functional characterization of ERFE-A260S in the hepatic cell system demonstrated its modifier role in iron overload by impairing the BMP/SMAD pathway. We herein described for the first time an ERFE polymorphism as a genetic modifier variant. This was with a mild effect on disease expression, under a multifactorial-like model, in a condition of iron-loading anemia due to ineffective erythropoiesis.

Coordinate regulation of ELF5 and EHF at the chr11p13 CF modifier region

E74-like factor 5 (ELF5) and ETS-homologous factor (EHF) are epithelial selective ETS family transcription factors (TFs) encoded by genes at chr11p13, a region associated with cystic fibrosis (CF) lung disease severity. EHF controls many key processes in lung epithelial function so its regulatory mechanisms are important. Using CRISPR/Cas9 technology, we removed three key cis-regulatory elements (CREs) from the chr11p13 region and also activated multiple open chromatin sites with CRISPRa in airway epithelial cells. Deletion of the CREs caused subtle changes in chromatin architecture and site-specific increases in EHF and ELF5. CRISPRa had most effect on ELF5 transcription. ELF5 levels are low in airway cells but higher in LNCaP (prostate) and T47D (breast) cancer cells. ATAC-seq in these lines revealed novel peaks of open chromatin at the 5' end of chr11p13 associated with an expressed ELF5 gene. Furthermore, 4C-seq assays identified direct interactions between the active ELF5 promoter and sites within the EHF locus, suggesting coordinate regulation between these TFs. ChIP-seq for ELF5 in T47D cells revealed ELF5 occupancy within EHF introns 1 and 6, and siRNA-mediated depletion of ELF5 enhanced EHF expression. These results define a new role for ELF5 in lung epithelial biology.

Gene modifiers of cystic fibrosis lung disease: A systematic review

BACKGROUND

Lung disease is the major source of morbidity and mortality in cystic fibrosis (CF), with large variability in severity between patients. Although accurate prediction of lung disease severity would be extremely useful, no robust methods exist. Twin and sibling studies have highlighted the importance of non-cystic fibrosis transmembrane conductance regulator (CFTR) genes in determining lung disease severity but how these impact on the severity in CF remains unclear.

METHODS

A systematic review was undertaken to answer the question "In patients with CF which non-CFTR genes modify the severity of lung disease?" The method for this systematic review was based upon the "Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA)" statement, with a narrative synthesis of results planned.

RESULTS

A total of 1168 articles were screened for inclusion, with 275 articles undergoing detailed assessment for inclusion. One hundred and forty articles were included. Early studies focused on candidate genes, whereas more recent studies utilized genome-wide approaches and also examined epigenetic mechanisms, gene expression, and therapeutic response.

DISCUSSION

A large body of evidence regarding non-CFTR gene modifiers of lung disease severity has been generated, examining a wide array of genes. Limitations to existing studies include heterogeneity in outcome measures used, limited replication, and relative lack of clinical impact. Future work examining non-CFTR gene modifiers will have to overcome these limitations if gene modifiers are to have a meaningful role in the care of patients with CF.

Phenotypic severity in a family with MEND syndrome is directly associated with the accumulation of potentially functional variants of cholesterol homeostasis genes

Background

Male EBP disorder with neurologic defects (MEND) syndrome is an X‐linked disease caused by hypomorphic mutations in the EBP (emopamil‐binding protein) gene. Modifier genes may explain the clinical variability among individuals who share a primary mutation.

Methods

We studied four males (Patient 1 to Patient 4) exhibiting a descending degree of phenotypic severity from a family with MEND syndrome. To identify candidate modifier genes that explain the phenotypic variability, variants of homeostasis cholesterol genes identified by whole‐exome sequencing (WES) were ranked according to the predicted magnitude of their effect through an in‐house scoring system.

Results

Twenty‐seven from 105 missense variants found in 45 genes of the four exomes were considered significant (−5 to −9 scores). We found a direct genotype–phenotype association based on the differential accumulation of potentially functional gene variants among males. Patient 1 exhibited 17 variants, both Patients 2 and 3 exhibited nine variants, and Patient 4 exhibited only five variants.

Conclusion

We conclude that APOA5 (rs3135506), ABCA1 (rs9282541), and APOB (rs679899 and rs12714225) are the most relevant candidate modifier genes in this family. Relative accumulation of the deficiencies associated with variants of these genes along with other lesser deficiencies in other genes appears to explain the variable expressivity in MEND syndrome.

GEMPROT: visualization of the impact on the protein of the genetic variants found on each haplotype

SUMMARY

When analyzing sequence data, genetic variants are considered one by one, taking no account of whether or not they are found in the same individual. However, variant combinations might be key players in some diseases as variants that are neutral on their own can become deleterious when associated together. GEMPROT is a new analysis tool that allows, from a phased vcf file, to visualize the consequences of the genetic variants on the protein. At the level of an individual, the program shows the variants on each of the two protein sequences and the Pfam functional protein domains. When data on several individuals are available, GEMPROT lists the haplotypes found in the sample and can compare the haplotype distributions between different sub-groups of individuals. By offering a global visualization of the gene with the genetic variants present, GEMPROT makes it possible to better understand the impact of combinations of genetic variants on the protein sequence.

AVAILABILITY AND IMPLEMENTATION

GEMPROT is freely available at https://github.com/TaniaCuppens/GEMPROT. An on-line version is also available at http://med-laennec.univ-brest.fr/GEMPROT/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

100 Years of evolving gene-disease complexities and scientific debutants

It's been over 100 years since the word `gene' is around and progressively evolving in several scientific directions. Time-to-time technological advancements have heavily revolutionized the field of genomics, especially when it's about, e.g. triple code development, gene number proposition, genetic mapping, data banks, gene-disease maps, catalogs of human genes and genetic disorders, CRISPR/Cas9, big data and next generation sequencing, etc. In this manuscript, we present the progress of genomics from pea plant genetics to the human genome project and highlight the molecular, technical and computational developments. Studying genome and epigenome led to the fundamentals of development and progression of human diseases, which includes chromosomal, monogenic, multifactorial and mitochondrial diseases. World Health Organization has classified, standardized and maintained all human diseases, when many academic and commercial online systems are sharing information about genes and linking to associated diseases. To efficiently fathom the wealth of this biological data, there is a crucial need to generate appropriate gene annotation repositories and resources. Our focus has been how many gene-disease databases are available worldwide and which sources are authentic, timely updated and recommended for research and clinical purposes. In this manuscript, we have discussed and compared 43 such databases and bioinformatics applications, which enable users to connect, explore and, if possible, download gene-disease data.

Genetics and epigenetic factors of Wilson disease

Wilson disease (WD) is a complex condition due to copper accumulation mainly in the liver and brain. The genetic base of WD is represented by pathogenic mutations of the copper-transporting gene ATP7B with consequent lack of copper excretion through the biliary tract. ATP7B is the only gene so far identified and known to be responsible for the development of the disease. Our understanding of the disease has been evolving as functional studies have associated specific disease-causing mutations with specific copper-transporter impairments. The most frequent variant in patients of European descent is the H1069Q missense mutation and it has been associated with protein misfolding, aberrant phosphorylation of the P-domain, and altered ATP binding orientation and affinity. Conversely, there is much less understanding of the relation between the genotype and the clinical manifestations of WD. WD is characterized by a highly varied and unpredictable presentation with different combined hepatic, neurological, and psychiatric symptoms. Several studies have attempted to correlate genotype and phenotype but the most recent evidences on larger populations failed to identify a relation between genotype and clinical presentations. Given that so far also modifier genes have not shown convincing association with WD, there is growing interest to identify epigenetic mechanisms of gene expression regulation as underlying the onset and progression of WD phenotype. Evidence from animal models indicated changes in methionine metabolism regulation with possible effects on DNA methylation. Mouse models of WD have indicated transcript level changes of genes related to DNA methylation in fetal and adult livers. And finally, evidence is accumulating regarding DNA methylation changes in patients with WD. It is unexplored how ATP7B genetic mutations combine with epigenetic changes to affect the phenotype. In conclusion, WD is a genetic disease with a complex regulation of its phenotype that includes molecular genetics and epigenetic mechanisms.

Association of GNAS imprinting defects and deletions of chromosome 2 in two patients: clues explaining phenotypic heterogeneity in pseudohypoparathyroidism type 1B/iPPSD3

BACKGROUND

The term pseudohypoparathyroidism (PHP) describes disorders derived from resistance to the parathyroid hormone. Albright hereditary osteodystrophy (AHO) is a disorder with several physical features that can occur alone or in association with PHP. The subtype 1B, classically associated with resistance to PTH and TSH, derives from the epigenetic dysregulation of the GNAS locus. Patients showing features of AHO were described, but no explanation for such phenotypic heterogeneity is available. An AHO-like phenotype was associated with the loss of genetic information stored in chromosome 2q37, making this genomic region an interesting object of study as it could contain modifier genes involved in the development of AHO features in patients with GNAS imprinting defects. The present study aimed to screen a series of 65 patients affected with GNAS imprinting defects, with or without signs of AHO, for the presence of 2q37 deletions in order to find genes involved in the clinical variability.

RESULTS

The molecular investigations performed on our cohort of patients with GNAS imprinting defects identified two overlapping terminal deletions of the long arm of chromosome 2. The smaller deletion was of approximately 3 Mb and contained 38 genes, one or more of which is potentially involved in the clinical presentation. Patients with the deletions were both affected by a combination of the most pathognomic AHO-like features, brachydactyly, cognitive impairment and/or behavioural defects. Our results support the hypothesis that additional genetic factors besides GNAS methylation defects are involved in the development of a complex phenotype in the subgroup of patients showing signs of AHO.

CONCLUSIONS

For the first time, the present work describes PHP patients with hormone resistance and AHO signs simultaneously affected by GNAS imprinting defects and 2q37 deletions. Although further studies are needed to confirm the cause of these two rare molecular alterations and to identify candidate genes, this finding provides novel interesting clues for the identification of factors involved in the still unexplained clinical variability observed in PHP1B.

Whole-exome sequencing identifies novel pathogenic variants across the ATP7B gene and some modifiers of Wilson's disease phenotype

BACKGROUND & AIMS

Wilson's disease (WD) is an autosomal recessive disorder associated with disease-causing alterations across the ATP7B gene, with highly variable symptoms and age of onset. We aimed to assess whether the clinical variability of WD relates to modifier genes.

METHODS

A total of 248 WD patients were included, of whom 148 were diagnosed after age of 17. Human exome libraries were constructed using AmpliSeq technology and sequenced using the IonProton platform.

RESULTS

ATP7B p.His1069Gln mutation was present in 215 patients, with 112 homozygotes and 103 heterozygotes. Three other mutations: p.Gln1351Ter, p.Trp779Ter and c.3402delC were identified in >10 patients. Among patients, 117 had a homozygous mutation, 101 were compound heterozygotes, 27 had one heterozygous mutation, and 3 other patients had no identifiable pathogenic variant of ATP7B. Sixteen mutations were novel, found as part of a compound mutation or as a sole, homozygous mutation. For disease phenotype prediction, age at diagnosis was a deciding factor, while frameshift allelic variants of ATP7B and being male increased the odds of developing a neurological phenotype. Rare allelic variants in ESD and INO80 increased and decreased chances for the neurological phenotype, respectively, while rare variants in APOE and MBD6 decreased the chances of WD early manifestation. Compound mutations contributed to earlier age of onset.

CONCLUSIONS

In a Polish population, genetic screening for WD may help genotype for four variants (p.His1069Gln, p.Gln1351Ter, p.Trp779Ter and c.3402delC), with direct sequencing of all ATP7B amplicons as a second diagnostic step. We also identified some allelic variants that may modify a WD phenotype.

Exploring genetic modifiers of Gaucher disease: The next horizon

Gaucher disease is an autosomal recessive lysosomal storage disorder resulting from mutations in the gene GBA1 that lead to a deficiency in the enzyme glucocerebrosidase. Accumulation of the enzyme's substrates, glucosylceramide and glucosylsphingosine, results in symptoms ranging from skeletal and visceral involvement to neurological manifestations. Nonetheless, there is significant variability in clinical presentations amongst patients, with limited correlation between genotype and phenotype. Contributing to this clinical variation are genetic modifiers that influence the phenotypic outcome of the disorder. In this review, we explore the role of genetic modifiers in Mendelian disorders and describe methods to facilitate their discovery. In addition, we provide examples of candidate modifiers of Gaucher disease, explore their relevance in the development of potential therapeutics, and discuss the impact of GBA1 and modifying mutations on other more common diseases like Parkinson disease. Identifying these important modulators of Gaucher phenotype may ultimately unravel the complex relationship between genotype and phenotype and lead to improved counseling and treatments.

Lnc-EPB41-Protein Interactions Associated with Congenital Pouch Colon

Congenital Pouch Colon (CPC) is a rare anorectal anomaly common to northwestern India, specifically Rajasthan. Despite efforts to understand the clinical genetic makeup of CPC, no attempt on identifying non-coding RNAs was done. We have earlier reported CPC’s rare variants from whole exome sequencing (WES) across 18 affected samples in a total of 64 subjects. A Smith–Waterman algorithm was used to infer a couple of lncRNAs from WES samples of CPC with predictions from the Noncode database. Further screening and quantification using polymerase chain reaction (PCR), we ascertained interactions using Micro Scale Thermophoresis (MST). We report the role of lnc-EPB41-1-1 shown to be promiscuously interacting with KIF13A substantiating their role in regulation.

Advances in the discovery of genetic risk factors for complex forms of neurodegenerative disorders: contemporary approaches, success, challenges and prospects

Neurodegenerative diseases constitute a large proportion of disorders in elderly, majority being sporadic in occurrence with ∼5-10% familial. A strong genetic component underlies the Mendelian forms but nongenetic factors together with genetic vulnerability contributes to the complex sporadic forms. Several gene discoveries in the familial forms have provided novel insights into the pathogenesis of neurodegeneration with implications for treatment. Conversely, findings from genetic dissection of the sporadic forms, despite large genomewide association studies and more recently whole exome and whole genome sequencing, have been limited. This review provides a concise account of the genetics that we know, the pathways that they implicate, the challenges that are faced and the prospects that are envisaged for the sporadic, complex forms of neurodegenerative diseases, taking four most common conditions, namely Alzheimer's disease, Parkinson's disease, amyotrophic lateral sclerosis and Huntington disease as examples. Poor replication across studies, inability to establish genotype-phenotype correlations and the overall failure to predict risk and/or prevent disease in this group poses a continuing challenge. Among others, clinical heterogeneity emerges as the most important impediment warranting newer approaches. Advanced computational and system biology tools to analyse the big data are being generated and the alternate strategy such as subgrouping of case-control cohorts based on deep phenotyping using the principles of Ayurveda to overcome current limitation of phenotype heterogeneity seem to hold promise. However, at this point, with advances in discovery genomics and functional analysis of putative determinants with translation potential for the complex forms being minimal, stem cell therapies are being attempted as potential interventions. In this context, the possibility to generate patient derived induced pluripotent stem cells, mutant/gene/genome correction through CRISPR/Cas9 technology and repopulating the specific brain regions with corrected neurons, which may fulfil the dream of personalized medicine have been mentioned briefly. Understanding disease pathways/biology using this technology, with implications for development of novel therapeutics are optimistic expectations in the near future.

Evaluating the Genetics of Common Variable Immunodeficiency: Monogenetic Model and Beyond

Common variable immunodeficiency (CVID) is the most frequent symptomatic primary immunodeficiency characterized by recurrent infections, hypogammaglobulinemia and poor response to vaccines. Its diagnosis is made based on clinical and immunological criteria, after exclusion of other diseases that can cause similar phenotypes. Currently, less than 20% of cases of CVID have a known underlying genetic cause. We have analyzed whole-exome sequencing and copy number variants data of 36 children and adolescents diagnosed with CVID and healthy relatives to estimate the proportion of monogenic cases. We have replicated an association of CVID to p.C104R in TNFRSF13B and reported the second case of homozygous patient to date. Our results also identify five causative genetic variants in LRBA, CTLA4, NFKB1, and PIK3R1, as well as other very likely causative variants in PRKCD, MAPK8, or DOCK8 among others. We experimentally validate the effect of the LRBA stop-gain mutation which abolishes protein production and downregulates the expression of CTLA4, and of the frameshift indel in CTLA4 producing expression downregulation of the protein. Our results indicate a monogenic origin of at least 15–24% of the CVID cases included in the study. The proportion of monogenic patients seems to be lower in CVID than in other PID that have also been analyzed by whole exome or targeted gene panels sequencing. Regardless of the exact proportion of CVID monogenic cases, other genetic models have to be considered for CVID. We propose that because of its prevalence and other features as intermediate penetrancies and phenotypic variation within families, CVID could fit with other more complex genetic scenarios. In particular, in this work, we explore the possibility of CVID being originated by an oligogenic model with the presence of heterozygous mutations in interacting proteins or by the accumulation of detrimental variants in particular immunological pathways, as well as perform association tests to detect association with rare genetic functional variation in the CVID cohort compared to healthy controls.

Exploring digenic inheritance in arrhythmogenic cardiomyopathy

Background

Arrhythmogenic cardiomyopathy (ACM) is an inherited genetic disorder, characterized by the substitution of heart muscle with fibro-fatty tissue and severe ventricular arrhythmias, often leading to heart failure and sudden cardiac death. ACM is considered a monogenic disorder, but the low penetrance of mutations identified in patients suggests the involvement of additional genetic or environmental factors.

Methods

We used whole exome sequencing to investigate digenic inheritance in two ACM families where previous diagnostic tests have revealed a PKP2 mutation in all affected and some healthy individuals. In family members with PKP2 mutations we determined all genes that harbor variants in affected but not in healthy carriers or vice versa. We computationally prioritized the most likely candidates, focusing on known ACM genes and genes related to PKP2 through protein interactions, functional relationships, or shared biological processes.

Results

We identified four candidate genes in family 1, namely DAG1, DAB2IP, CTBP2 and TCF25, and eleven candidate genes in family 2. The most promising gene in the second family is TTN, a gene previously associated with ACM, in which the affected individual harbors two rare deleterious-predicted missense variants, one of which is located in the protein’s only serine kinase domain.

Conclusions

In this study we report genes that might act as digenic players in ACM pathogenesis, on the basis of co-segregation with PKP2 mutations. Validation in larger cohorts is still required to prove the utility of this model.

Electronic supplementary material

The online version of this article (10.1186/s12881-017-0503-7) contains supplementary material, which is available to authorized users.

Transcriptomic profile of cystic fibrosis patients identifies type I interferon response and ribosomal stalk proteins as potential modifiers of disease severity

Cystic Fibrosis (CF) is the most common monogenic disease among people of Western European descent and caused by mutations in the CFTR gene. However, the disease severity is immensely variable even among patients with similar CFTR mutations due to the possible effect of 'modifier genes'. To identify genetic modifiers, we applied RNA-seq based transcriptomic analyses in CF patients with a mild and severe lung phenotype. Global gene expression and enrichment analyses revealed that genes of the type I interferon response and ribosomal stalk proteins are potential modifiers of CF related lung dysfunction. The results provide a new set of CF modifier genes with possible implications as new therapeutic targets for the treatment of CF.

The evolving role of ankyrin-B in cardiovascular disease

Over the past decade, ankyrin-B has been identified as a prominent player in cardiac physiology. Ankyrin-B has a multitude of functions, with roles in expression, localization, and regulation of proteins critical for cardiac excitability, cytoskeletal integrity, and signaling. Furthermore, human ANK2 variants that result in ankyrin-B loss of function are associated with "ankyrin-B syndrome," a complex cardiac phenotype that may include bradycardia and heart rate variability, conduction block, atrial fibrillation, QT interval prolongation, and potentially fatal catecholaminergic polymorphic ventricular tachycardia. However, our understanding of the molecular mechanisms underlying ankyrin-B function at baseline and in disease is still not fully developed owing to the complexity of ankyrin-B gene regulation, number of ankyrin-B-associated molecules, multiple roles of ankyrin-B in the heart and other organs that modulate cardiac function, and a host of unexpected clinical phenotypes. In this review, we summarize known roles of ankyrin-B in the heart and the impact of ankyrin-B dysfunction in animal models and in human disease as well as highlight important new findings illustrating the complexity of ankyrin-B signaling.

Rare variants of small effect size in neuronal excitability genes influence clinical outcome in Japanese cases of SCN1A truncation-positive Dravet syndrome

Dravet syndrome (DS) is a rare, devastating form of childhood epilepsy that is often associated with mutations in the voltage-gated sodium channel gene, SCN1A. There is considerable variability in expressivity within families, as well as among individuals carrying the same primary mutation, suggesting that clinical outcome is modulated by variants at other genes. To identify modifier gene variants that contribute to clinical outcome, we sequenced the exomes of 22 individuals at both ends of a phenotype distribution (i.e., mild and severe cognitive condition). We controlled for variation associated with different mutation types by limiting inclusion to individuals with a de novo truncation mutation resulting in SCN1A haploinsufficiency. We performed tests aimed at identifying 1) single common variants that are enriched in either phenotypic group, 2) sets of common or rare variants aggregated in and around genes associated with clinical outcome, and 3) rare variants in 237 candidate genes associated with neuronal excitability. While our power to identify enrichment of a common variant in either phenotypic group is limited as a result of the rarity of mild phenotypes in individuals with SCN1A truncation variants, our top candidates did not map to functional regions of genes, or in genes that are known to be associated with neurological pathways. In contrast, we found a statistically-significant excess of rare variants predicted to be damaging and of small effect size in genes associated with neuronal excitability in severely affected individuals. A KCNQ2 variant previously associated with benign neonatal seizures is present in 3 of 12 individuals in the severe category. To compare our results with the healthy population, we performed a similar analysis on whole exome sequencing data from 70 Japanese individuals in the 1000 genomes project. Interestingly, the frequency of rare damaging variants in the same set of neuronal excitability genes in healthy individuals is nearly as high as in severely affected individuals. Rather than a single common gene/variant modifying clinical outcome in SCN1A-related epilepsies, our results point to the cumulative effect of rare variants with little to no measurable phenotypic effect (i.e., typical genetic background) unless present in combination with a disease-causing truncation mutation in SCN1A.

Simultaneous Expression of ABCA4 and GPR143 Mutations: A Complex Phenotypic Manifestation

Purpose

To describe the complex, overlapping phenotype expressed in a two generation family harboring pathogenic mutations in the ABCA4 and GPR143 genes.

Methods

Clinical evaluation of a two generation family included quantitative autofluorescence imaging (qAF, 488-nm excitation) using a modified confocal scanning laser ophthalmoscope equipped with an internal fluorescent reference to account for varying laser power detector sensitivity, spectral-domain optical coherence tomography, and full-field ERG testing. Complete sequencing of the ABCA4 and GPR143 genes was carried out in each individual.

Results

Affected individuals presented with bull's eye lesions and qAF levels above the 95% confidence interval for healthy eyes; full-field ERG revealed no generalized rod dysfunction but mild implicit time delays in cone responses. Complete sequencing of the ABCA4 gene revealed two disease-causing mutations, p.L541P and p.G1961E; and mutational phase was confirmed in each unaffected parent. Further examination in the affected patients revealed a peripheral “mud-splattered” pattern of hypopigmented RPE after which sequencing of GPR143 revealed a novel missense variant, p.Y157C. The GPR143 variant segregated from the father who did not exhibit any indications of retinal disease with the exception of an abnormal near-infrared autofluorescence (NIR-AF) signal distribution in the macula.

Conclusions

An individual carrying both ABCA4 and GPR143 disease-causing mutations can express a complex, overlapping phenotype associated with both Stargardt disease and X-linked ocular albinism (OA1). The absence of OA1-related disease changes (with the exception of NIR-AF changes associated with melanin distribution) in the father may be indicative of mild expressivity or variable gene penetrance.

Mechanisms of suppression: The wiring of genetic resilience

Recent analysis of genome sequences has identified individuals that are healthy despite carrying severe disease-associated mutations. A possible explanation is that these individuals carry a second genomic perturbation that can compensate for the detrimental effects of the disease allele, a phenomenon referred to as suppression. In model organisms, suppression interactions are generally divided into two classes: genomic suppressors which are secondary mutations in the genome that bypass a mutant phenotype, and dosage suppression interactions in which overexpression of a suppressor gene rescues a mutant phenotype. Here, we describe the general properties of genomic and dosage suppression, with an emphasis on the budding yeast. We propose that suppression interactions between genetic variants are likely relevant for determining the penetrance of human traits. Consequently, an understanding of suppression mechanisms may guide the discovery of protective variants in healthy individuals that carry disease alleles, which could direct the rational design of new therapeutics.

Modeling diseases in multiple mouse strains for precision medicine studies

The genetic basis of the phenotypic variability observed in patients can be studied in mice by generating disease models through genetic or chemical interventions in many genetic backgrounds where the clinical phenotypes can be assessed and used for genome-wide association studies (GWAS). This is particularly relevant for rare disorders, where patients sharing identical mutations can present with a wide variety of symptoms, but there are not enough number of patients to ensure statistical power of GWAS. Inbred strains are homozygous for each loci, and their single nucleotide polymorphisms catalogs are known and freely available, facilitating the bioinformatics and reducing the costs of the study, since it is not required to genotype every mouse. This kind of approach can be applied to pharmacogenomics studies as well.

Transcript levels of plastin 3 and neuritin 1 modifier genes in spinal muscular atrophy siblings

BACKGROUND

In single gene disorders, patients with the same genotype may have variations in severity. One of the main factors affecting disease severity is modifier genes. Spinal muscular atrophy (SMA) is an autosomal recessive neuromuscular disorder caused by degeneration of alpha motor neurons. Plastin 3 (PLS3) is a phenotypic modifier of SMA, and neuritin 1 (NRN1) has also been suggested as a possible modifier gene. The aim of the present study was therefore to analyze PLS3 and NRN1 expression in SMA siblings in four families.

METHODS

The study group consisted of four SMA families with seven with discordant phenotype and two affected siblings. Total RNA was isolated from whole blood. PLS3 and NRN1 expression was analyzed on quantitative real-time polymerase chain reaction.

RESULTS

In family 1 only NRN1 expression was increased in the mildly affected sister. In family 2 only PLS3 had a modifier effect. Family 3, which had type III siblings with identical clinical phenotypes, had similar PLS3 expression between the siblings but no NRN1 expression. In family 4, neither PLS3 nor NRN1 had any correlation with severity.

CONCLUSION

On analysis of the expression of NRN1 in SMA patients for the first time, NRN1 could be a potential modifier gene. PLS3 expression does not always modify SMA phenotype. In patients with no modifier effect of known genes, genome sequencing and transcriptome analysis are promising for the identification of novel modifiers and understanding of SMA pathophysiology.

Genetics of Combined Pituitary Hormone Deficiency: Roadmap into the Genome Era

The genetic basis for combined pituitary hormone deficiency (CPHD) is complex, involving 30 genes in a variety of syndromic and nonsyndromic presentations. Molecular diagnosis of this disorder is valuable for predicting disease progression, avoiding unnecessary surgery, and family planning. We expect that the application of high throughput sequencing will uncover additional contributing genes and eventually become a valuable tool for molecular diagnosis. For example, in the last 3 years, six new genes have been implicated in CPHD using whole-exome sequencing. In this review, we present a historical perspective on gene discovery for CPHD and predict approaches that may facilitate future gene identification projects conducted by clinicians and basic scientists. Guidelines for systematic reporting of genetic variants and assigning causality are emerging. We apply these guidelines retrospectively to reports of the genetic basis of CPHD and summarize modes of inheritance and penetrance for each of the known genes. In recent years, there have been great improvements in databases of genetic information for diverse populations. Some issues remain that make molecular diagnosis challenging in some cases. These include the inherent genetic complexity of this disorder, technical challenges like uneven coverage, differing results from variant calling and interpretation pipelines, the number of tolerated genetic alterations, and imperfect methods for predicting pathogenicity. We discuss approaches for future research in the genetics of CPHD.