Digenic inheritance and Mendelian disease

A knowledge graph approach to predict and interpret disease-causing gene interactions

BACKGROUND

Understanding the impact of gene interactions on disease phenotypes is increasingly recognised as a crucial aspect of genetic disease research. This trend is reflected by the growing amount of clinical research on oligogenic diseases, where disease manifestations are influenced by combinations of variants on a few specific genes. Although statistical machine-learning methods have been developed to identify relevant genetic variant or gene combinations associated with oligogenic diseases, they rely on abstract features and black-box models, posing challenges to interpretability for medical experts and impeding their ability to comprehend and validate predictions. In this work, we present a novel, interpretable predictive approach based on a knowledge graph that not only provides accurate predictions of disease-causing gene interactions but also offers explanations for these results.

RESULTS

We introduce BOCK, a knowledge graph constructed to explore disease-causing genetic interactions, integrating curated information on oligogenic diseases from clinical cases with relevant biomedical networks and ontologies. Using this graph, we developed a novel predictive framework based on heterogenous paths connecting gene pairs. This method trains an interpretable decision set model that not only accurately predicts pathogenic gene interactions, but also unveils the patterns associated with these diseases. A unique aspect of our approach is its ability to offer, along with each positive prediction, explanations in the form of subgraphs, revealing the specific entities and relationships that led to each pathogenic prediction.

CONCLUSION

Our method, built with interpretability in mind, leverages heterogenous path information in knowledge graphs to predict pathogenic gene interactions and generate meaningful explanations. This not only broadens our understanding of the molecular mechanisms underlying oligogenic diseases, but also presents a novel application of knowledge graphs in creating more transparent and insightful predictors for genetic research.

Identifying digenic disease genes via machine learning in the Undiagnosed Diseases Network

Rare diseases affect millions of people worldwide, and discovering their genetic causes is challenging. More than half of the individuals analyzed by the Undiagnosed Diseases Network (UDN) remain undiagnosed. The central hypothesis of this work is that many of these rare genetic disorders are caused by multiple variants in more than one gene. However, given the large number of variants in each individual genome, experimentally evaluating combinations of variants for potential to cause disease is currently infeasible. To address this challenge, we developed the digenic predictor (DiGePred), a random forest classifier for identifying candidate digenic disease gene pairs by features derived from biological networks, genomics, evolutionary history, and functional annotations. We trained the DiGePred classifier by using DIDA, the largest available database of known digenic-disease-causing gene pairs, and several sets of non-digenic gene pairs, including variant pairs derived from unaffected relatives of UDN individuals. DiGePred achieved high precision and recall in cross-validation and on a held-out test set (PR area under the curve > 77%), and we further demonstrate its utility by using digenic pairs from the recent literature. In contrast to other approaches, DiGePred also appropriately controls the number of false positives when applied in realistic clinical settings. Finally, to enable the rapid screening of variant gene pairs for digenic disease potential, we freely provide the predictions of DiGePred on all human gene pairs. Our work enables the discovery of genetic causes for rare non-monogenic diseases by providing a means to rapidly evaluate variant gene pairs for the potential to cause digenic disease.

Variant filtering, digenic variants, and other challenges in clinical sequencing: a lesson from fibrillinopathies

Genome-scale high-throughput sequencing enables the detection of unprecedented numbers of sequence variants. Variant filtering and interpretation are facilitated by mutation databases, in silico tools, and population-based reference datasets such as ExAC/gnomAD, while variants are classified using the ACMG/AMP guidelines. These methods, however, pose clinically relevant challenges. We queried the gnomAD dataset for (likely) pathogenic variants in genes causing autosomal-dominant disorders. Furthermore, focusing on the fibrillinopathies Marfan syndrome (MFS) and congenital contractural arachnodactyly (CCA), we screened 500 genomes of our patients for co-occurring variants in FBN1 and FBN2. In gnomAD, we detected 2653 (likely) pathogenic variants in 253 genes associated with autosomal-dominant disorders, enabling the estimation of variant-filtering thresholds and disease predisposition/prevalence rates. In our database, we discovered two families with hitherto unreported co-occurrence of FBN1/FBN2 variants causing phenotypes with mixed or modified MFS/CCA clinical features. We show that (likely) pathogenic gnomAD variants may be more frequent than expected and are challenging to classify according to the ACMG/AMP guidelines as well as that fibrillinopathies are likely underdiagnosed and may co-occur. Consequently, selection of appropriate frequency cutoffs, recognition of digenic variants, and variant classification represent considerable challenges in variant interpretation. Neglecting these challenges may lead to incomplete or missed diagnoses.

Digenic inheritance of human primary microcephaly delineates centrosomal and non-centrosomal pathways

Primary microcephaly (PM) is characterized by a small head since birth and is vastly heterogeneous both genetically and phenotypically. While most cases are monogenic, genetic interactions between Aspm and Wdr62 have recently been described in a mouse model of PM. Here, we used two complementary, holistic in vivo approaches: high throughput DNA sequencing of multiple PM genes in human patients with PM, and genome‐edited zebrafish modeling for the digenic inheritance of PM. Exomes of patients with PM showed a significant burden of variants in 75 PM genes, that persisted after removing monogenic causes of PM (e.g., biallelic pathogenic variants in CEP152). This observation was replicated in an independent cohort of patients with PM, where a PM gene panel showed in addition that the burden was carried by six centrosomal genes. Allelic frequencies were consistent with digenic inheritance. In zebrafish, non‐centrosomal gene casc5 −/− produced a severe PM phenotype, that was not modified by centrosomal genes aspm or wdr62 invalidation. A digenic, quadriallelic PM phenotype was produced by aspm and wdr62. Our observations provide strong evidence for digenic inheritance of human PM, involving centrosomal genes. Absence of genetic interaction between casc5 and aspm or wdr62 further delineates centrosomal and non‐centrosomal pathways in PM.

Recent Advances in Craniosynostosis

Craniosynostosis is a pathologic craniofacial disorder and is defined as the premature fusion of one or more cranial (calvarial) sutures. Cranial sutures are fibrous joints consisting of nonossified mesenchymal cells that play an important role in the development of healthy craniofacial skeletons. Early fusion of these sutures results in incomplete brain development that may lead to complications of several severe medical conditions including seizures, brain damage, mental delay, complex deformities, strabismus, and visual and breathing problems. As a congenital disease, craniosynostosis has a heterogeneous origin that can be affected by genetic and epigenetic alterations, teratogens, and environmental factors and make the syndrome highly complex. To date, approximately 200 syndromes have been linked to craniosynostosis. In addition to being part of a syndrome, craniosynostosis can be nonsyndromic, formed without any additional anomalies. More than 50 nuclear genes that relate to craniosynostosis have been identified. Besides genetic factors, epigenetic factors like microRNAs and mechanical forces also play important roles in suture fusion. As craniosynostosis is a multifactorial disorder, evaluating the craniosynostosis syndrome requires and depends on all the information obtained from clinical findings, genetic analysis, epigenetic or environmental factors, or gene modulators. In this review, we will focus on embryologic and genetic studies, as well as epigenetic and environmental studies. We will discuss published studies and correlate the findings with unknown aspects of craniofacial disorders.

Genome sequencing and implications for rare disorders

The practice of genomic medicine stands to revolutionize our approach to medical care, and to realize this goal will require discovery of the relationship between rare variation at each of the ~ 20,000 protein-coding genes and their consequent impact on individual health and expression of Mendelian disease. The step-wise evolution of broad-based, genome-wide cytogenetic and molecular genomic testing approaches (karyotyping, chromosomal microarray [CMA], exome sequencing [ES]) has driven much of the rare disease discovery to this point, with genome sequencing representing the newest member of this team. Each step has brought increased sensitivity to interrogate individual genomic variation in an unbiased method that does not require clinical prediction of the locus or loci involved. Notably, each step has also brought unique limitations in variant detection, for example, the low sensitivity of ES for detection of triploidy, and of CMA for detection of copy neutral structural variants. The utility of genome sequencing (GS) as a clinical molecular diagnostic test, and the increased sensitivity afforded by addition of long-read sequencing or other -omics technologies such as RNAseq or metabolomics, are not yet fully explored, though recent work supports improved sensitivity of variant detection, at least in a subset of cases. The utility of GS will also rely upon further elucidation of the complexities of genetic and allelic heterogeneity, multilocus rare variation, and the impact of rare and common variation at a locus, as well as advances in functional annotation of identified variants. Much discovery remains to be done before the potential utility of GS is fully appreciated.

Insights into genetics, human biology and disease gleaned from family based genomic studies

Identifying genes and variants contributing to rare disease phenotypes and Mendelian conditions informs biology and medicine, yet potential phenotypic consequences for variation of >75% of the ~20,000 annotated genes in the human genome are lacking. Technical advances to assess rare variation genome-wide, particularly exome sequencing (ES), enabled establishment in the United States of the National Institutes of Health (NIH)-supported Centers for Mendelian Genomics (CMGs) and have facilitated collaborative studies resulting in novel "disease gene" discoveries. Pedigree-based genomic studies and rare variant analyses in families with suspected Mendelian conditions have led to the elucidation of hundreds of novel disease genes and highlighted the impact of de novo mutational events, somatic variation underlying nononcologic traits, incompletely penetrant alleles, phenotypes with high locus heterogeneity, and multilocus pathogenic variation. Herein, we highlight CMG collaborative discoveries that have contributed to understanding both rare and common diseases and discuss opportunities for future discovery in single-locus Mendelian disorder genomics. Phenotypic annotation of all human genes; development of bioinformatic tools and analytic methods; exploration of non-Mendelian modes of inheritance including reduced penetrance, multilocus variation, and oligogenic inheritance; construction of allelic series at a locus; enhanced data sharing worldwide; and integration with clinical genomics are explored. Realizing the full contribution of rare disease research to functional annotation of the human genome, and further illuminating human biology and health, will lay the foundation for the Precision Medicine Initiative.

PIEZO1 Hypomorphic Variants in Congenital Lymphatic Dysplasia Cause Shape and Hydration Alterations of Red Blood Cells

PIEZO1 is a cation channel activated by mechanical force. It plays an important physiological role in several biological processes such as cardiovascular, renal, endothelial and hematopoietic systems. Two different diseases are associated with alteration in the DNA sequence of PIEZO1: (i) dehydrated hereditary stomatocytosis (DHS1, #194380), an autosomal dominant hemolytic anemia caused by gain-of-function mutations; (ii) lymphatic dysplasia with non-immune fetal hydrops (LMPH3, #616843), an autosomal recessive condition caused by biallelic loss-of-function mutations. We analyzed a 14-year-old boy affected by severe lymphatic dysplasia already present prenatally, with peripheral edema, hydrocele, and chylothoraces. By whole exome sequencing, we identified compound heterozygosity for PIEZO1, with one splicing and one deletion mutation, the latter causing the formation of a premature stop codon that leads to mRNA decay. The functional analysis of the erythrocytes of the patient highlighted altered hydration with the intracellular loss of the potassium content and structural abnormalities with anisopoikolocytosis and presence of both spherocytes and stomatocytes. This novel erythrocyte trait, sharing features with both hereditary spherocytosis and overhydrated hereditary stomatocytosis, complements the clinical features associated with loss-of-function mutations of PIEZO1 in the context of the generalized lymphatic dysplasia of LMPH3 type.

De novo mutations in Caudal Type Homeo Box transcription Factor 2 (CDX2) in patients with persistent cloaca

The cloaca is an embryonic cavity that is divided into the urogenital sinus and rectum upon differentiation of the cloacal epithelium triggered by tissue-specific transcription factors including CDX2. Defective differentiation leads to persistent cloaca in humans (PC), a phenotype recapitulated in Cdx2 mutant mice. PC is linked to hypo/hyper-vitaminosis A. Although no gene has ever been identified, there is a strong evidence for a genetic contribution to PC. We applied whole-exome sequencing and copy-number-variants analyses to 21 PC patients and their unaffected parents. The damaging p.Cys132* and p.Arg237His de novo CDX2 variants were identified in two patients. These variants altered the expression of CYP26A1, a direct CDX2 target encoding the major retinoic acid (RA)-degrading enzyme. Other RA genes, including the RA-receptor alpha, were also mutated. Genes governing the development of cloaca-derived structures were recurrently mutated and over-represented in the basement-membrane components set (q-value < 1.65 × 10-6). Joint analysis of the patients' profile highlighted the extracellular matrix-receptor interaction pathway (MsigDBID: M7098, FDR: q-value < 7.16 × 10-9). This is the first evidence that PC is genetic, with genes involved in the RA metabolism at the lead. Given the CDX2 de novo variants and the role of RA, our observations could potentiate preventive measures. For the first time, a gene recapitulating PC in mouse models is found mutated in humans.

The classification, genetic diagnosis and modelling of monogenic autoinflammatory disorders

Monogenic autoinflammatory disorders are an increasingly heterogeneous group of conditions characterised by innate immune dysregulation. Improved genetic sequencing in recent years has led not only to the discovery of a plethora of conditions considered to be 'autoinflammatory', but also the broadening of the clinical and immunological phenotypic spectra seen in these disorders. This review outlines the classification strategies that have been employed for monogenic autoinflammatory disorders to date, including the primary innate immune pathway or the dominant cytokine implicated in disease pathogenesis, and highlights some of the advantages of these models. Furthermore, the use of the term 'autoinflammatory' is discussed in relation to disorders that cross the innate and adaptive immune divide. The utilisation of next-generation sequencing (NGS) in this population is examined, as are potential in vivo and in vitro methods of modelling to determine pathogenicity of novel genetic findings. Finally, areas where our understanding can be improved are highlighted, such as phenotypic variability and genotype-phenotype correlations, with the aim of identifying areas of future research.

Early-onset epileptic encephalopathy and severe developmental delay in an association with de novo double mutations in NF1 and MAGEL2

Advance in the exome-wide sequencing analysis contributes to identifying hundreds of genes that are associated with early-onset epileptic encephalopathy and neurodevelopmental disorders. On the basis of massive sequencing data, functional interactions among different genes are suggested to explain the common molecular pathway underlying the pathogenic process of these disorders. However, the relevance of such interactions with the phenotypic severity or variety in an affected individual remains elusive. In this report, we present a 45-year-old woman with neurofibromatosis type 1 (NF1), infantile-onset epileptic encephalopathy, and severe developmental delay. Whole-exome sequencing identified de novo pathogenic mutations in NF1 and the Schaaf-Yang syndrome-associated gene, MAGEL2. Literature-curated interaction data predicted that NF1 and MAGEL2 proteins were closely connected in this network via their common interacting proteins. Direct conversion of fibroblasts into neurons in vitro showed that neuronal cells from 9 patients with NF1 expressed significantly lower levels of MAGEL2 (54%, p = 0.0047) than those from healthy individuals. These data provide the first evidence that pathogenic mutations of NF1 deregulate the expression of other neurodevelopmental disease-associated genes. De novo mutations in multiple genes may lead to severe developmental phenotypes through their cumulative effects or synergistic interactions.

Digenic inheritance and genetic modifiers

Digenic inheritance (DI) concerns pathologies with the simplest form of multigenic etiology, implicating more than 1 gene (and perhaps the environment). True DI is when biallelic or even triallelic mutations in 2 distinct genes, in cis or in trans, are necessary and sufficient to cause pathology with a defined diagnosis. In true DI, a heterozygous mutation in each of 2 genes alone is not associated with a recognizable phenotype. Well-documented diseases with true DI are so far rare and follow non-Mendelian inheritance. DI is also encountered when by serendipity, pathogenic mutations responsible for 2 distinct disease entities are co-inherited, leading to a mixed phenotype. Also, we can consider many true monogenic Mendelian conditions, which show impressively broad spectrum of phenotypes due to pseudo-DI, as a result of co-inheriting genetic modifiers (GMs). I am herewith reviewing examples of GM and embark on presenting some recent notable examples of true DI, with wider discussion of the literature. Undeniably, the advent of high throughput sequencing is bound to unravel more patients suffering with true DI conditions and elucidate many important GM, thus impacting precision medicine.

Hereditary stomatocytosis: An underdiagnosed condition

Hereditary stomatocytoses are a wide class of hemolytic anemias characterized by alterations of ionic flux with increased cation permeability that results in inappropriate shrinkage or swelling of the erythrocytes, and water lost or gained osmotically. The last few years have been crucial for new acquisitions in this field in terms of identifying new causative genes and of studying their pathogenetic mechanisms. This review summarizes the main features of erythrocyte membrane transport diseases, dividing them into forms with either isolated erythroid phenotype (nonsyndromic) or extra-hematological manifestations (syndromic), and focusing particularly on the most recent advances regarding dehydrated forms of hereditary stomatocytosis and familial pseudohyperkalemia.

The genetics of congenitally small brains

Primary microcephaly (PM) refers to a congenitally small brain, resulting from insufficient prenatal production of neurons, and serves as a model disease for brain volumic development. Known PM genes delineate several cellular pathways, among which the centriole duplication pathway, which provide interesting clues about the cellular mechanisms involved. The general interest of the genetic dissection of PM is illustrated by the convergence of Zika virus infection and PM gene mutations on congenital microcephaly, with CENPJ/CPAP emerging as a key target. Physical (protein-protein) and genetic (digenic inheritance) interactions of Wdr62 and Aspm have been demonstrated in mice, and should now be sought in humans using high throughput parallel sequencing of multiple PM genes in PM patients and control subjects, in order to categorize mutually interacting genes, hence delineating functional pathways in vivo in humans.

Understanding mutational effects in digenic diseases

To further our understanding of the complexity and genetic heterogeneity of rare diseases, it has become essential to shed light on how combinations of variants in different genes are responsible for a disease phenotype. With the appearance of a resource on digenic diseases, it has become possible to evaluate how digenic combinations differ in terms of the phenotypes they produce. All instances in this resource were assigned to two classes of digenic effects, annotated as true digenic and composite classes. Whereas in the true digenic class variants in both genes are required for developing the disease, in the composite class, a variant in one gene is sufficient to produce the phenotype, but an additional variant in a second gene impacts the disease phenotype or alters the age of onset. We show that a combination of variant, gene and higher-level features can differentiate between these two classes with high accuracy. Moreover, we show via the analysis of three digenic disorders that a digenic effect decision profile, extracted from the predictive model, motivates why an instance was assigned to either of the two classes. Together, our results show that digenic disease data generates novel insights, providing a glimpse into the oligogenic realm.

A systematic review of genetic syndromes with obesity

Syndromic monogenic obesity typically follows Mendelian patterns of inheritance and involves the co-presentation of other characteristics, such as mental retardation, dysmorphic features and organ-specific abnormalities. Previous reviews on obesity have reported 20 to 30 syndromes but no systematic review has yet been conducted on syndromic obesity. We searched seven databases using terms such as 'obesity', 'syndrome' and 'gene' to conduct a systematic review of literature on syndromic obesity. Our literature search identified 13,719 references. After abstract and full-text review, 119 relevant papers were eligible, and 42 papers were identified through additional searches. Our analysis of these 161 papers found that 79 obesity syndromes have been reported in literature. Of the 79 syndromes, 19 have been fully genetically elucidated, 11 have been partially elucidated, 27 have been mapped to a chromosomal region and for the remaining 22, neither the gene(s) nor the chromosomal location(s) have yet been identified. Interestingly, 54.4% of the syndromes have not been assigned a name, whereas 13.9% have more than one name. We report on organizational inconsistencies (e.g. naming discrepancies and syndrome classification) and provide suggestions for improvements. Overall, this review illustrates the need for increased clinical and genetic research on syndromes with obesity.

Genetic advances in craniosynostosis

Craniosynostosis, the premature ossification of one or more skull sutures, is a clinically and genetically heterogeneous congenital anomaly affecting approximately one in 2,500 live births. In most cases, it occurs as an isolated congenital anomaly, that is, nonsyndromic craniosynostosis (NCS), the genetic, and environmental causes of which remain largely unknown. Recent data suggest that, at least some of the midline NCS cases may be explained by two loci inheritance. In approximately 25-30% of patients, craniosynostosis presents as a feature of a genetic syndrome due to chromosomal defects or mutations in genes within interconnected signaling pathways. The aim of this review is to provide a detailed and comprehensive update on the genetic and environmental factors associated with NCS, integrating the scientific findings achieved during the last decade. Focus on the neurodevelopmental, imaging, and treatment aspects of NCS is also provided.

Genetic study of congenital bile-duct dilatation identifies de novo and inherited variants in functionally related genes

Background

Congenital dilatation of the bile-duct (CDD) is a rare, mostly sporadic, disorder that results in bile retention with severe associated complications. CDD affects mainly Asians. To our knowledge, no genetic study has ever been conducted.

Methods

We aim to identify genetic risk factors by a “trio-based” exome-sequencing approach, whereby 31 CDD probands and their unaffected parents were exome-sequenced. Seven-hundred controls from the local population were used to detect gene-sets significantly enriched with rare variants in CDD patients.

Results

Twenty-one predicted damaging de novo variants (DNVs; 4 protein truncating and 17 missense) were identified in several evolutionarily constrained genes (p < 0.01). Six genes carrying DNVs were associated with human developmental disorders involving epithelial, connective or bone morphologies (PXDN, RTEL1, ANKRD11, MAP2K1, CYLD, ACAN) and four linked with cholangio- and hepatocellular carcinomas (PIK3CA, TLN1 CYLD, MAP2K1). Importantly, CDD patients have an excess of DNVs in cancer-related genes (p < 0.025). Thirteen genes were recurrently mutated at different sites, forming compound heterozygotes or functionally related complexes within patients.

Conclusions

Our data supports a strong genetic basis for CDD and show that CDD is not only genetically heterogeneous but also non-monogenic, requiring mutations in more than one genes for the disease to develop. The data is consistent with the rarity and sporadic presentation of CDD.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-016-0236-z) contains supplementary material, which is available to authorized users.

Two locus inheritance of non-syndromic midline craniosynostosis via rare SMAD6 and common BMP2 alleles

Premature fusion of the cranial sutures (craniosynostosis), affecting 1 in 2000 newborns, is treated surgically in infancy to prevent adverse neurologic outcomes. To identify mutations contributing to common non-syndromic midline (sagittal and metopic) craniosynostosis, we performed exome sequencing of 132 parent-offspring trios and 59 additional probands. Thirteen probands (7%) had damaging de novo or rare transmitted mutations in SMAD6, an inhibitor of BMP – induced osteoblast differentiation (p<10⁻²⁰). SMAD6 mutations nonetheless showed striking incomplete penetrance (<60%). Genotypes of a common variant near BMP2 that is strongly associated with midline craniosynostosis explained nearly all the phenotypic variation in these kindreds, with highly significant evidence of genetic interaction between these loci via both association and analysis of linkage. This epistatic interaction of rare and common variants defines the most frequent cause of midline craniosynostosis and has implications for the genetic basis of other diseases.

DOI:http://dx.doi.org/10.7554/eLife.20125.001

The bones in the front, back and sides of the human skull are not fused to one another at birth in order to allow the brain to double in size during the first year of life and continue growing into adulthood. However, one in 2,000 infants is born with a condition called craniosynostosis in which some of these bones have already fused. This fusion prevents the skull from growing properly, and can lead to the brain becoming compressed. As such, surgeons routinely undo the fusion in these infants to allow the brain and skull to grow normally.

Eighty-five percent of craniosynostosis cases occur in infants with no other abnormalities, (called non-syndromic cases) and most have no other affected family member. It has therefore been unclear whether these infants have craniosynostosis due to a genetic or non-genetic cause. If the cause is genetic, it is also not clear whether a mutation in a single gene, the combined effects of many genes, or something in between is responsible.

Now, by focusing on a group of 191 infants with premature fusion of bones joined at the midline of the skull, Timberlake et al. asked if any of the approximately 20,000 genes in the human genome were altered more frequently in these infants than would be expected by chance. This search revealed that rare mutations that disable one copy of a gene called SMAD6 in combination with a common DNA variant near another gene called BMP2 account for about 7% of infants with midline forms of craniosynostosis. These genes are both known to regulate how bones form, which explains how the mutation of these genes could lead to craniosynostosis.

In all cases, the parents of these children were unaffected. This was typically because one parent had only the SMAD6 mutation while the other had only the common BMP2 variant; the transmission of both to their offspring resulted in craniosynostosis. The finding that a rare mutation’s effect is strongly modified by a common variant from another site in the genome is unprecedented. These findings will allow doctors to counsel families about the risk of having additional children with craniosynostosis.

Timberlake et al. next plan to study more patients with craniosynostosis to identify additional genes that contribute to this disease. They will also look at other diseases to see whether the combination of rare mutation and common DNA variant could be behind other unexplained disorders.

DOI:http://dx.doi.org/10.7554/eLife.20125.002

Clinical genomics: from a truly personal genome viewpoint

The path to Clinical Genomics is punctuated by our understanding of what types of DNA structural and sequence variation contribute to disease, the many technical challenges to detect such variation genome-wide, and the initial struggles to interpret personal genome variation in the context of disease. This review describes one perspective of the development of clinical genomics; whereas the experimental challenges, and hurdles to overcoming them, might be deemed readily apparent, the non-technical issues for clinical implementation may be less obvious. Some of these latter challenges, including: (1) informed consent, (2) privacy, (3) what constitutes potentially pathogenic variation contributing to disease, (4) disease penetrance in populations, and (5) the genetic architecture of disease, and the struggles sometimes faced for solutions, are highlighted using illustrative examples.

DIDA: A curated and annotated digenic diseases database

DIDA (DIgenic diseases DAtabase) is a novel database that provides for the first time detailed information on genes and associated genetic variants involved in digenic diseases, the simplest form of oligogenic inheritance. The database is accessible via http://dida.ibsquare.be and currently includes 213 digenic combinations involved in 44 different digenic diseases. These combinations are composed of 364 distinct variants, which are distributed over 136 distinct genes. The web interface provides browsing and search functionalities, as well as documentation and help pages, general database statistics and references to the original publications from which the data have been collected. The possibility to submit novel digenic data to DIDA is also provided. Creating this new repository was essential as current databases do not allow one to retrieve detailed records regarding digenic combinations. Genes, variants, diseases and digenic combinations in DIDA are annotated with manually curated information and information mined from other online resources. Next to providing a unique resource for the development of new analysis methods, DIDA gives clinical and molecular geneticists a tool to find the most comprehensive information on the digenic nature of their diseases of interest.

Heterozygous RTEL1 mutations are associated with familial pulmonary fibrosis

Pulmonary fibrosis is a fatal disease with progressive loss of respiratory function. Defective telomere maintenance leading to telomere shortening is a cause of pulmonary fibrosis, as mutations in the telomerase component genes TERT (reverse transcriptase) and TERC (RNA component) are found in 15% of familial pulmonary fibrosis (FPF) cases. However, so far, about 85% of FPF remain genetically uncharacterised.Here, in order to identify new genetic causes of FPF, we performed whole-exome sequencing, with a candidate-gene approach, of 47 affected subjects from 35 families with FPF without TERT and TERC mutations.We identified heterozygous mutations in regulator of telomere elongation helicase 1 (RTEL1) in four families. RTEL1 is a DNA helicase with roles in DNA replication, genome stability, DNA repair and telomere maintenance. The heterozygous RTEL1 mutations segregated as an autosomal dominant trait in FPF, and were predicted by structural analyses to severely affect the function and/or stability of RTEL1. In agreement with this, RTEL1-mutated patients exhibited short telomeres in comparison with age-matched controls.Our results provide evidence that heterozygous RTEL1 mutations are responsible for FPF and, thereby, extend the clinical spectrum of RTEL1 deficiency. Thus, RTEL1 enlarges the number of telomere-associated genes implicated in FPF.

Excess of rare, inherited truncating mutations in autism

To assess the relative impact of inherited and de novo variants on autism risk, we generated a comprehensive set of exonic single nucleotide variants (SNVs) and copy number variants (CNVs) from 2,377 autism families. We find that private, inherited truncating SNVs in conserved genes are enriched in probands (odds ratio=1.14, p=0.0002) compared to unaffected siblings, an effect with significant maternal transmission bias to sons. We also observe a bias for inherited CNVs, specifically for small (<100 kbp), maternally inherited events (p=0.01) that are enriched in CHD8 target genes (p=7.4×10⁻³). Using a logistic regression model, we show that private truncating SNVs and rare, inherited CNVs are statistically independent autism risk factors, with odds ratios of 1.11 (p=0.0002) and 1.23 (p=0.01), respectively. This analysis identifies a second class of candidate genes (e.g., RIMS1, CUL7, and LZTR1) where transmitted mutations may create a sensitized background but are unlikely to be completely penetrant.

Somatic mosaicism: implications for disease and transmission genetics

Nearly all of the genetic material among cells within an organism is identical. However, single-nucleotide variants (SNVs), small insertions/deletions (indels), copy-number variants (CNVs), and other structural variants (SVs) continually accumulate as cells divide during development. This process results in an organism composed of countless cells, each with its own unique personal genome. Thus, every human is undoubtedly mosaic. Mosaic mutations can go unnoticed, underlie genetic disease or normal human variation, and may be transmitted to the next generation as constitutional variants. We review the influence of the developmental timing of mutations, the mechanisms by which they arise, methods for detecting mosaic variants, and the risk of passing these mutations on to the next generation.

A community-based resource for automatic exome variant-calling and annotation in Mendelian disorders

Background

Mendelian disorders are mostly caused by single mutations in the DNA sequence of a gene, leading to a phenotype with pathologic consequences. Whole Exome Sequencing of patients can be a cost-effective alternative to standard genetic screenings to find causative mutations of genetic diseases, especially when the number of cases is limited. Analyzing exome sequencing data requires specific expertise, high computational resources and a reference variant database to identify pathogenic variants.

Results

We developed a database of variations collected from patients with Mendelian disorders, which is automatically populated thanks to an associated exome-sequencing pipeline. The pipeline is able to automatically identify, annotate and store insertions, deletions and mutations in the database. The resource is freely available online http://exome.tigem.it. The exome sequencing pipeline automates the analysis workflow (quality control and read trimming, mapping on reference genome, post-alignment processing, variation calling and annotation) using state-of-the-art software tools. The exome-sequencing pipeline has been designed to run on a computing cluster in order to analyse several samples simultaneously. The detected variants are annotated by the pipeline not only with the standard variant annotations (e.g. allele frequency in the general population, the predicted effect on gene product activity, etc.) but, more importantly, with allele frequencies across samples progressively collected in the database itself, stratified by Mendelian disorder.

Conclusions

We aim at providing a resource for the genetic disease community to automatically analyse whole exome-sequencing samples with a standard and uniform analysis pipeline, thus collecting variant allele frequencies by disorder. This resource may become a valuable tool to help dissecting the genotype underlying the disease phenotype through an improved selection of putative patient-specific causative or phenotype-associated variations.

From genes to health - challenges and opportunities

In genome science, the advancement in high-throughput sequencing technologies and bioinformatics analysis is facilitating the better understanding of Mendelian and complex trait inheritance. Charting the genetic basis of complex diseases - including pediatric cancer, and interpreting huge amount of next-generation sequencing data are among the major technical challenges to be overcome in order to understand the molecular basis of various diseases and genetic disorders. In this review, we provide insights into some major challenges currently hindering a better understanding of Mendelian and complex trait inheritance, and thus impeding medical benefits to patients.

Multiplicity of hormone-secreting tumors: common themes about cause, expression, and management

CONTEXT

Multiplicity of hormone-secreting tumors occurs in a substantial portion of hormone-excess states. Multiplicity increases the difficulty of management and drives the selection of special strategies.

EVIDENCE ACQUISITION

This is a synthesis from publications about tumor development and expression, and also about types of clinical strategy for hormone-secreting tumors.

EVIDENCE SYNTHESIS

Comparisons were made between patient groups with solitary tumors vs those with multiple tumors. Major themes with clinical relevance emerged. Usually, tumor multiplicity develops from a genetic susceptibility in all cells of a tissue. This applies to hormone-secreting tumors that begin as either polyclonal (such as in the parathyroids of familial hypocalciuric hypercalcemia) or monoclonal tumors (such as in the parathyroids of multiple endocrine neoplasia type 1 [MEN1]). High penetrance of a hereditary tumor frequently results in bilaterality and in several other types of multiplicity. Managements are better for the hormone excess than for the associated cancers. Management strategies can be categorized broadly as ablation that is total, subtotal, or zero. Examples are discussed for each category, and 1 example of each category is named here: 1) total ablation of the entire tissue with effort to replace ablated functions (for example, in C-cell neoplasia of multiple endocrine neoplasia type 2); 2) subtotal ablation with increased likelihood of persistent disease or recurrent disease (for example, in the parathyroid tumors of MEN1); or 3) no ablation of tissue with or without the use of pharmacotherapy (for example, with blockers for secretion of stomach acid in gastrinomas of MEN1).

CONCLUSIONS

Tumor multiplicity usually arises from defects in all cells of the precursor tissue. Even the optimized managements involve compromises. Still, an understanding of pathophysiology and of therapeutic options should guide optimized management.

Genetic biomarkers in aortopathy

The field of aortopathy, in common with other genomic disorders, is undergoing a revolution. This is largely driven by the implementation of newer forms of genetic sequencing (massively parallel or next-generation sequencing). Advantages conferred by this technology include reduced costs, reduced sequencing time and the ability to simultaneously test multiple genes. This has a significant advantage in the identification of genes disrupted in heritable aortopathies. These advances are enabling scientists and clinicians to identify key molecular pathways; translating fundamental genetic findings into a better understanding of disease mechanisms is ultimately leading to effective treatments. In outlining contemporary knowledge of genetic biomarkers in aortopathy we seek to demonstrate that the era of genomically orientated decision-making is here.

Digenic mutational inheritance of the integrin alpha 7 and the myosin heavy chain 7B genes causes congenital myopathy with left ventricular non-compact cardiomyopathy

Background

We report an Italian family in which the proband showed a severe phenotype characterized by the association of congenital fiber type disproportion (CFTD) with a left ventricular non-compaction cardiomyopathy (LVNC). This study was focused on the identification of the responsible gene/s.

Methods and results

Using the whole-exome sequencing approach, we identified the proband homozygous missense mutations in two genes, the myosin heavy chain 7B (MYH7B) and the integrin alpha 7 (ITGA7). Both genes are expressed in heart and muscle tissues, and both mutations were predicted to be deleterious and were not found in the healthy population.

The R890C mutation in the MYH7B gene segregated with the LVNC phenotype in the examined family. It was also found in one unrelated patient affected by LVNC, confirming a causative role in cardiomyopathy.

The E882K mutation in the ITGA7 gene, a key component of the basal lamina of muscle fibers, was found only in the proband, suggesting a role in CFTD.

Conclusions

This study identifies two novel disease genes. Mutation in MYH7B causes a classical LVNC phenotype, whereas mutation in ITGA7 causes CFTD. Both phenotypes represent alterations of skeletal and cardiac muscle maturation and are usually not severe. The severe phenotype of the proband is most likely due to a synergic effect of these two mutations.

This study provides new insights into the genetics underlying Mendelian traits and demonstrates a role for digenic inheritance in complex phenotypes.

Deletions of recessive disease genes: CNV contribution to carrier states and disease-causing alleles

Over 1200 recessive disease genes have been described in humans. The prevalence, allelic architecture, and per-genome load of pathogenic alleles in these genes remain to be fully elucidated, as does the contribution of DNA copy-number variants (CNVs) to carrier status and recessive disease. We mined CNV data from 21,470 individuals obtained by array-comparative genomic hybridization in a clinical diagnostic setting to identify deletions encompassing or disrupting recessive disease genes. We identified 3212 heterozygous potential carrier deletions affecting 419 unique recessive disease genes. Deletion frequency of these genes ranged from one occurrence to 1.5%. When compared with recessive disease genes never deleted in our cohort, the 419 recessive disease genes affected by at least one carrier deletion were longer and located farther from known dominant disease genes, suggesting that the formation and/or prevalence of carrier CNVs may be affected by both local and adjacent genomic features and by selection. Some subjects had multiple carrier CNVs (307 subjects) and/or carrier deletions encompassing more than one recessive disease gene (206 deletions). Heterozygous deletions spanning multiple recessive disease genes may confer carrier status for multiple single-gene disorders, for complex syndromes resulting from the combination of two or more recessive conditions, or may potentially cause clinical phenotypes due to a multiply heterozygous state. In addition to carrier mutations, we identified homozygous and hemizygous deletions potentially causative for recessive disease. We provide further evidence that CNVs contribute to the allelic architecture of both carrier and recessive disease-causing mutations. Thus, a complete recessive carrier screening method or diagnostic test should detect CNV alleles.

How genetic modifiers influence the phenotype of spinal muscular atrophy and suggest future therapeutic approaches

Both complex disorders and monogenetic diseases are often modulated in their phenotype by further genetic, epigenetic or extrinsic factors. This gives rise to extensive phenotypic variability and potentially protection from disease manifestations, known as incomplete penetrance. Approaches including whole transcriptome, exome, genome, methylome or proteome analyses of highly discordant phenotypes in a few individuals harboring mutations at the same locus can help to identify these modifiers. This review describes the complexity of modifying factors of one of the most frequent autosomal recessively inherited disorders in humans, spinal muscular atrophy (SMA). We will outline how this knowledge contributes to understanding of the regulatory networks and molecular pathology of SMA and how this knowledge will influence future approaches to therapies.

New insights into the genetic basis of TAR (thrombocytopenia-absent radii) syndrome

Thrombocytopenia with absent radii (TAR) syndrome is a rare disorder combining specific skeletal abnormalities with a reduced platelet count. Rare proximal microdeletions of 1q21.1 are found in the majority of patients but are also found in unaffected parents. Recently it was shown that TAR syndrome is caused by the compound inheritance of a low-frequency noncoding SNP and a rare null allele in RBM8A, a gene encoding the exon-junction complex subunit member Y14 located in the deleted region. This finding provides new insight into the complex inheritance pattern and new clues to the molecular mechanisms underlying TAR syndrome. We discuss TAR syndrome in the context of abnormal phenotypes associated with proximal and distal 1q21.1 microdeletion and microduplications with incomplete penetrance and variable expressivity.

Molecular consequences of animal breeding

The phenotypic diversity in domestic animals provides a unique opportunity to study genotype-phenotype relationships. The identification of causal mutations provides an insight into what types of mutations have contributed to phenotypic evolution in domestic animals. Whole genome sequencing has revealed that fixation of null alleles that inactivate genes, which are essential under natural conditions but disadvantageous on the farm, has not been a common mechanism for genetic adaptation in domestic animals. Numerous examples have been revealed where structural changes cause specific phenotypic effects by altering transcriptional regulation. An emerging feature is also the evolution of alleles by the accumulation of several consecutive mutations which affect gene function.