The continuum of causality in human genetic disorders

The co-occurrence of genetic variants in the TYR and OCA2 genes confers susceptibility to albinism

Although rare genetic conditions are mostly caused by DNA sequence alterations that functionally disrupt individual genes, large-scale studies using genome sequencing have started to unmask additional complexity. Understanding how combinations of variants in different genes shape human phenotypes is expected to provide important insights into the clinical and genetic heterogeneity of rare disorders. Here, we use albinism, an archetypal rare condition associated with hypopigmentation, as an exemplar for the study of genetic interactions. We analyse data from the Genomics England 100,000 Genomes Project alongside a cohort of 1120 individuals with albinism, and investigate the effect of dual heterozygosity for the combination of two established albinism-related variants: TYR:c.1205 G > A (p.Arg402Gln) [rs1126809] and OCA2:c.1327 G > A (p.Val443Ile) [rs74653330]. As each of these changes alone is insufficient to cause disease when present in the heterozygous state, we sought evidence of synergistic effects. We show that, when both variants are present, the probability of receiving a diagnosis of albinism is significantly increased (odds ratio 12.8; 95% confidence interval 6.0 - 24.7; p-value 2.1 ×10-8). Further analyses in an independent cohort, the UK Biobank, support this finding and highlight that heterozygosity for the TYR:c.1205 G > A and OCA2:c.1327 G > A variant combination is associated with statistically significant alterations in visual acuity and central retinal thickness (traits that are considered albinism endophenotypes). The approach discussed in this report opens up new avenues for the investigation of oligogenic patterns in apparently Mendelian disorders.

CRISPR screening uncovers a long-range enhancer for ONECUT1 in pancreatic differentiation and links a diabetes risk variant

Functional enhancer annotation is critical for understanding tissue-specific transcriptional regulation and prioritizing disease-associated non-coding variants. However, unbiased enhancer discovery in disease-relevant contexts remains challenging. To identify enhancers pertinent to diabetes, we conducted a CRISPR interference (CRISPRi) screen in the human pluripotent stem cell (hPSC) pancreatic differentiation system. Among the enhancers identified, we focused on an enhancer we named ONECUT1e-664kb, ∼664 kb from the ONECUT1 promoter. Previous studies have linked ONECUT1 coding mutations to pancreatic hypoplasia and neonatal diabetes. We found that homozygous deletion of ONECUT1e-664kb in hPSCs leads to a near-complete loss of ONECUT1 expression and impaired pancreatic differentiation. ONECUT1e-664kb contains a type 2 diabetes-associated variant (rs528350911) disrupting a GATA motif. Introducing the risk variant into hPSCs reduced binding of key pancreatic transcription factors (GATA4, GATA6, and FOXA2), supporting its causal role in diabetes. This work highlights the utility of unbiased enhancer discovery in disease-relevant settings for understanding monogenic and complex disease.

Exploring the Genotype-Phenotype Correlations in a Child with Inherited Seizure and Thrombocytopenia by Digenic Network Analysis

Understanding the correlation between genotype and phenotype remains challenging for modern genetics. Digenic network analysis may provide useful models for understanding complex phenotypes that traditional Mendelian monogenic models cannot explain. Clinical data, whole exome sequencing data, in silico, and machine learning analysis were combined to construct a digenic network that may help unveil the complex genotype-phenotype correlations in a child presenting with inherited seizures and thrombocytopenia. The proband inherited a maternal heterozygous missense variant in SCN1A (NM_001165963.4:c.2722G>A) and a paternal heterozygous missense variant in MYH9 (NM_002473.6:c.3323A>C). In silico analysis showed that these two variants may be pathogenic for inherited seizures and thrombocytopenia in the proband. Moreover, focusing on 230 epilepsy-associated genes and 35 thrombopoiesis genes, variant call format data of the proband were analyzed using machine learning tools (VarCoPP 2.0) and Digenic Effect predictor. A digenic network was constructed, and SCN1A and MYH9 were found to be core genes in the network. Further analysis showed that MYH9 might be a modifier of SCN1A, and the variant in MYH9 might not only influence the severity of SCN1A-related seizure but also lead to thrombocytopenia in the bone marrow. In addition, another eight variants might also be co-factors that account for the proband's complex phenotypes. Our data show that as a supplement to the traditional Mendelian monogenic model, digenic network analysis may provide reasonable models for the explanation of complex genotype-phenotype correlations.

Polygenic Risk Scores in Predicting Coronary Artery Disease in Symptomatic Patients. A Validation Study

AIM

Clinical risk scores for coronary artery disease (CAD) are used in clinical practice to select patients for diagnostic testing and therapy. Several studies have proposed that polygenic risk scores (PRSs) can improve the prediction of CAD, but the scores need to be validated in clinical populations with accurately characterized phenotypes. We assessed the predictive power of the three most promising PRSs for the prediction of coronary atherosclerosis and obstructive CAD.

METHODS

This study was conducted on 943 symptomatic patients with suspected CAD for whom the phenotype was accurately characterized using anatomic and functional imaging. Previously published genome-wide polygenic scores were generated to compare a genetic model based on PRSs with a model based on clinical data. The test and PRS cohorts were predominantly Caucasian of northern European ancestry.

RESULTS

All three PRSs predicted coronary atherosclerosis and obstructive CAD statistically significantly. The predictive accuracy of the models combining clinical data and different PRSs varied between 0.778 and 0.805 in terms of the area under the receiver operating characteristic (AUROC), being close to the model including only clinical variables (AUROC 0.769). The difference between the clinical model and combined clinical + PRS model was not significant for PRS1 (p=0.627) and PRS3 (p=0.061). Only PRS2 slightly improved the predictive power of the model (p=0.04). The likelihood ratios showed the very weak diagnostic power of all PRSs.

CONCLUSION

The addition of PRSs to conventional risk factors did not clinically significantly improve the predictive accuracy for either coronary atherosclerosis or obstructive CAD, showing that current PRSs are not justified for routine clinical use in CAD.

Lethal phenotypes in Mendelian disorders

PURPOSE

Existing resources that characterize the essentiality status of genes are based on either proliferation assessment in human cell lines, viability evaluation in mouse knockouts, or constraint metrics derived from human population sequencing studies. Several repositories document phenotypic annotations for rare disorders; however, there is a lack of comprehensive reporting on lethal phenotypes.

METHODS

We queried Online Mendelian Inheritance in Man for terms related to lethality and classified all Mendelian genes according to the earliest age of death recorded for the associated disorders, from prenatal death to no reports of premature death. We characterized the genes across these lethality categories, examined the evidence on viability from mouse models and explored how this information could be used for novel gene discovery.

RESULTS

We developed the Lethal Phenotypes Portal to showcase this curated catalog of human essential genes. Differences in the mode of inheritance, physiological systems affected, and disease class were found for genes in different lethality categories, as well as discrepancies between the lethal phenotypes observed in mouse and human.

CONCLUSION

We anticipate that this resource will aid clinicians in the diagnosis of early lethal conditions and assist researchers in investigating the properties that make these genes essential for human development.

Variant Effect Prediction in the Age of Machine Learning

Over the years, many computational methods have been created for the analysis of the impact of single amino acid substitutions resulting from single-nucleotide variants in genome coding regions. Historically, all methods have been supervised and thus limited by the inadequate sizes of experimentally curated data sets and by the lack of a standardized definition of variant effect. The emergence of unsupervised, deep learning (DL)-based methods raised an important question: Can machines learn the language of life from the unannotated protein sequence data well enough to identify significant errors in the protein "sentences"? Our analysis suggests that some unsupervised methods perform as well or better than existing supervised methods. Unsupervised methods are also faster and can, thus, be useful in large-scale variant evaluations. For all other methods, however, their performance varies by both evaluation metrics and by the type of variant effect being predicted. We also note that the evaluation of method performance is still lacking on less-studied, nonhuman proteins where unsupervised methods hold the most promise.

CRISPR Screening Uncovers a Long-Range Enhancer for ONECUT1 in Pancreatic Differentiation and Links a Diabetes Risk Variant

Functional enhancer annotation is a valuable first step for understanding tissue-specific transcriptional regulation and prioritizing disease-associated non-coding variants for investigation. However, unbiased enhancer discovery in physiologically relevant contexts remains a major challenge. To discover regulatory elements pertinent to diabetes, we conducted a CRISPR interference screen in the human pluripotent stem cell (hPSC) pancreatic differentiation system. Among the enhancers uncovered, we focused on a long-range enhancer ∼664 kb from the ONECUT1 promoter, since coding mutations in ONECUT1 cause pancreatic hypoplasia and neonatal diabetes. Homozygous enhancer deletion in hPSCs was associated with a near-complete loss of ONECUT1 gene expression and compromised pancreatic differentiation. This enhancer contains a confidently fine-mapped type 2 diabetes associated variant (rs528350911) which disrupts a GATA motif. Introduction of the risk variant into hPSCs revealed substantially reduced binding of key pancreatic transcription factors (GATA4, GATA6 and FOXA2) on the edited allele, accompanied by a slight reduction of ONECUT1 transcription, supporting a causal role for this risk variant in metabolic disease. This work expands our knowledge about transcriptional regulation in pancreatic development through the characterization of a long-range enhancer and highlights the utility of enhancer discovery in disease-relevant settings for understanding monogenic and complex disease.

Lethal phenotypes in Mendelian disorders

Essential genes are those whose function is required for cell proliferation and/or organism survival. A gene's intolerance to loss-of-function can be allocated within a spectrum, as opposed to being considered a binary feature, since this function might be essential at different stages of development, genetic backgrounds or other contexts. Existing resources that collect and characterise the essentiality status of genes are based on either proliferation assessment in human cell lines, embryonic and postnatal viability evaluation in different model organisms, and gene metrics such as intolerance to variation scores derived from human population sequencing studies. There are also several repositories available that document phenotypic annotations for rare disorders in humans such as the Online Mendelian Inheritance in Man (OMIM) and the Human Phenotype Ontology (HPO) knowledgebases. This raises the prospect of being able to use clinical data, including lethality as the most severe phenotypic manifestation, to further our characterisation of gene essentiality. Here we queried OMIM for terms related to lethality and classified all Mendelian genes into categories, according to the earliest age of death recorded for the associated disorders, from prenatal death to no reports of premature death. To showcase this curated catalogue of human essential genes, we developed the Lethal Phenotypes Portal (https://lethalphenotypes.research.its.qmul.ac.uk), where we also explore the relationships between these lethality categories, constraint metrics and viability in cell lines and mouse. Further analysis of the genes in these categories reveals differences in the mode of inheritance of the associated disorders, physiological systems affected and disease class. We highlight how the phenotypic similarity between genes in the same lethality category combined with gene family/group information can be used for novel disease gene discovery. Finally, we explore the overlaps and discrepancies between the lethal phenotypes observed in mouse and human and discuss potential explanations that include differences in transcriptional regulation, functional compensation and molecular disease mechanisms. We anticipate that this resource will aid clinicians in the diagnosis of early lethal conditions and assist researchers in investigating the properties that make these genes essential for human development.

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.

The emergence of genotypic divergence and future precision medicine applications

Genotypic divergence is a term adapted from population genetics and intimately linked to evolution. We use divergence here to emphasize the differences that set individuals apart in any cohort. The history of genetics is filled with descriptions of genotypic differences, but causal inference of interindividual biological variation has been scarce. We suggest that the practice of precision medicine requires a divergent approach, an approach dependent on the causal interpretation of previous convergent (and preliminary) knowledge in the field. This knowledge has relied on convergent descriptive syndromology (lumping), which has overemphasized a reductionistic gene determinism on the quest of seeking associations without causal understanding. Regulatory variants with small effect and somatic mutations are some of the modifying factors that lead to incomplete penetrance and intrafamilial variable expressivity often observed in apparently monogenic clinical disorders. A truly divergent approach to precision medicine requires splitting, that is, the consideration of different layers of genetic phenomena that interact causally in a nonlinear fashion. This chapter reviews convergences and divergences in genetics and genomics, aiming to discuss what can be causally understood to approximate the as-yet utopian lands of Precision Medicine for patients with neurodegenerative disorders.

Genetically transitional disease: a new concept in genomic medicine

Traditional classification of genetic diseases as monogenic and polygenic has lagged far behind scientific progress. In this opinion article, we propose and define a new terminology, genetically transitional disease (GTD), referring to cases where a large-effect mutation is necessary, but not sufficient, to cause disease. This leads to a working disease nosology based on gradients of four types of genetic architecture: monogenic, polygenic, GTD, and mixed. We present four scenarios under which GTD may occur; namely, subsets of traditionally Mendelian disease, modifiable Tier 1 monogenic conditions, variable penetrance, and situations where a genetic mutational spectrum produces qualitatively divergent pathologies. The implications of the new nosology in precision medicine are discussed, in which therapeutic options may target the molecular cause or the disease phenotype.

A genome-wide association study identifies distinct variants associated with pulmonary function among European and African ancestries from the UK Biobank

Pulmonary function is an indicator of well-being, and pulmonary pathologies are the third major cause of death worldwide. We analysed the UK Biobank genome-wide association summary statistics of pulmonary function for Europeans and individuals of recent African descent to identify variants associated with the trait in the two ancestries. Here, we show 627 variants in Europeans and 3 in Africans associated with three pulmonary function parameters. In addition to the 110 variants in Europeans previously reported to be associated with phenotypes related to pulmonary function, we identify 279 novel loci, including an ISX intergenic variant rs369476290 on chromosome 22 in Africans. Remarkably, we find no shared variants among Africans and Europeans. Furthermore, enrichment analyses of variants separately for each ancestry background reveal significant enrichment for terms related to pulmonary phenotypes in Europeans but not Africans. Further analysis of studies of pulmonary phenotypes reveals that individuals of European background are disproportionally overrepresented in datasets compared to Africans, with the gap widening over the past five years. Our findings extend our understanding of the different variants that modify the pulmonary function in Africans and Europeans, a promising finding for future GWASs and medical studies.

A Patient with Corticobasal Syndrome and Progressive Non-Fluent Aphasia (CBS-PNFA), with Variants in ATP7B, SETX, SORL1, and FOXP1 Genes

Our aim was to analyze the phenotypic-genetic correlations in a patient diagnosed with early onset corticobasal syndrome with progressive non-fluent aphasia (CBS-PNFA), characterized by predominant apraxia of speech, accompanied by prominent right-sided upper-limb limb-kinetic apraxia, alien limb phenomenon, synkinesis, myoclonus, mild cortical sensory loss, and right-sided hemispatial neglect. Whole-exome sequencing (WES) identified rare single heterozygous variants in ATP7B (c.3207C>A), SORL1 (c.352G>A), SETX (c.2385_2387delAAA), and FOXP1 (c.1762G>A) genes. The functional analysis revealed that the deletion in the SETX gene changed the splicing pattern, which was accompanied by lower SETX mRNA levels in the patient's fibroblasts, suggesting loss-of-function as the underlying mechanism. In addition, the patient's fibroblasts demonstrated altered mitochondrial architecture with decreased connectivity, compared to the control individuals. This is the first association of the CBS-PNFA phenotype with the most common ATP7B pathogenic variant p.H1069Q, previously linked to Wilson's disease, and early onset Parkinson's disease. This study expands the complex clinical spectrum related to variants in well-known disease genes, such as ATP7B, SORL1, SETX, and FOXP1, corroborating the hypothesis of oligogenic inheritance. To date, the FOXP1 gene has been linked exclusively to neurodevelopmental speech disorders, while our study highlights its possible relevance for adult-onset progressive apraxia of speech, which guarantees further study.

Achievements, prospects and challenges in precision care for monogenic insulin-deficient and insulin-resistant diabetes

Integration of genomic and other data has begun to stratify type 2 diabetes in prognostically meaningful ways, but this has yet to impact on mainstream diabetes practice. The subgroup of diabetes caused by single gene defects thus provides the best example to date of the vision of 'precision diabetes'. Monogenic diabetes may be divided into primary pancreatic beta cell failure, and primary insulin resistance. In both groups, clear examples of genotype-selective responses to therapy have been advanced. The benign trajectory of diabetes due to pathogenic GCK mutations, and the sulfonylurea-hyperresponsiveness conferred by activating KCNJ11 or ABCC8 mutations, or loss-of-function HNF1A or HNF4A mutations, often decisively guide clinical management. In monogenic insulin-resistant diabetes, subcutaneous leptin therapy is beneficial in some severe lipodystrophy. Increasing evidence also supports use of 'obesity therapies' in lipodystrophic people even without obesity. In beta cell diabetes the main challenge is now implementation of the precision diabetes vision at scale. In monogenic insulin-resistant diabetes genotype-specific benefits are proven in far fewer patients to date, although further genotype-targeted therapies are being evaluated. The conceptual paradigm established by the insulin-resistant subgroup with 'adipose failure' may have a wider influence on precision therapy for common type 2 diabetes, however. For all forms of monogenic diabetes, population-wide genome sequencing is currently forcing reappraisal of the importance assigned to pathogenic mutations when gene sequencing is uncoupled from prior suspicion of monogenic diabetes.

The Role of Genome Sequencing in Neonatal Intensive Care Units

Genetic diseases disrupt the functionality of an infant's genome during fetal-neonatal adaptation and represent a leading cause of neonatal and infant mortality in the United States. Due to disease acuity, gene locus and allelic heterogeneity, and overlapping and diverse clinical phenotypes, diagnostic genome sequencing in neonatal intensive care units has required the development of methods to shorten turnaround times and improve genomic interpretation. From 2012 to 2021, 31 clinical studies documented the diagnostic and clinical utility of first-tier rapid or ultrarapid whole-genome sequencing through cost-effective identification of pathogenic genomic variants that change medical management, suggest new therapeutic strategies, and refine prognoses. Genomic diagnosis also permits prediction of reproductive recurrence risk for parents and surviving probands. Using implementation science and quality improvement, deployment of a genomic learning healthcare system will contribute to a reduction of neonatal and infant mortality through the integration of genome sequencing into best-practice neonatal intensive care.

Mendel paved the path toward understanding genetic diseases

Background

July 20th, 2022, marks the 200th anniversary of the “Father of Genetics,” Gregor Mendel’s birth. His experiments with pea plants established the fundamental principles of genetic inheritance.

Main text

In this article, a succinct review of literature is hereby done to answer two key questions: (1) How Mendel’s principles of genetic inheritance helped us understand Mendelian disorders? and (2) How the study of Mendelian disorders can help us understand complex diseases?

Conclusion

This literature review concludes that continued effort to understand the genetic basis of Mendelian disorders will improve our understanding and treatment strategies for complex diseases.

TREM2 Gene Compound Heterozygosity in Neurodegenerative Disorders

BACKGROUND

Homozygous variants of the TREM2 and TYROBP genes have been shown to be causative for multiple bone cysts and neurodegeneration leading to progressive dementia (NHD, Nasu-Hakola disease).

OBJECTIVE

To determine if biallelic variants of these genes and/or oligogenic inheritance could be responsible for a wider spectrum of neurodegenerative conditions.

METHODS

We analyzed 52 genes associated with neurodegenerative disorders using targeted next generation sequencing in a selected group of 29 patients (n = 14 Alzheimer's disease, n = 8 frontotemporal dementia, n = 7 amyotrophic lateral sclerosis) carrying diverse already determined rare variants in exon 2 of TREM2. Molecular modeling was used to get an insight into the potential effects of the mutation.

RESULTS

We identified a novel mutation c.401_406delinsTCTAT; p.(Asp134Valfs*55) in exon 3 of TREM2 in an Alzheimer's disease patient also carrying the p.Arg62His TREM2 variant. Molecular modeling revealed that the identified mutation prevents anchoring of the TREM2 protein in the membrane, leaving the core of the Ig-like domain intact.

CONCLUSION

Our results expand the spectrum of neurodegenerative diseases, where the carriers of biallelic mutations in TREM2 have been described for Alzheimer's disease, and highlight the impact of variant burden in other genes on phenotypic heterogeneity.

An intermediate-effect size variant in UMOD confers risk for chronic kidney disease

The kidney-specific gene UMOD encodes for uromodulin, the most abundant protein excreted in normal urine. Rare large-effect variants in UMOD cause autosomal dominant tubulointerstitial kidney disease (ADTKD), while common low-impact variants strongly associate with kidney function and the risk of chronic kidney disease (CKD) in the general population. It is unknown whether intermediate-effect variants in UMOD contribute to CKD. Here, candidate intermediate-effect UMOD variants were identified using large-population and ADTKD cohorts. Biological and phenotypical effects were investigated using cell models, in silico simulations, patient samples, and international databases and biobanks. Eight UMOD missense variants reported in ADTKD are present in the Genome Aggregation Database (gnomAD), with minor allele frequency (MAF) ranging from 10-5 to 10-3. Among them, the missense variant p.Thr62Pro is detected in ∼1/1,000 individuals of European ancestry, shows incomplete penetrance but a high genetic load in familial clusters of CKD, and is associated with kidney failure in the 100,000 Genomes Project (odds ratio [OR] = 3.99 [1.84 to 8.98]) and the UK Biobank (OR = 4.12 [1.32 to 12.85). Compared with canonical ADTKD mutations, the p.Thr62Pro carriers displayed reduced disease severity, with slower progression of CKD and an intermediate reduction of urinary uromodulin levels, in line with an intermediate trafficking defect in vitro and modest induction of endoplasmic reticulum (ER) stress. Identification of an intermediate-effect UMOD variant completes the spectrum of UMOD-associated kidney diseases and provides insights into the mechanisms of ADTKD and the genetic architecture of CKD.

UMOD and the architecture of kidney disease

The identification of genetic factors associated with the risk, onset, and progression of kidney disease has the potential to provide mechanistic insights and therapeutic perspectives. In less than two decades, technological advances yielded a trove of information on the genetic architecture of chronic kidney disease. The spectrum of genetic influence ranges from (ultra)rare variants with large effect size, involved in Mendelian diseases, to common variants, often non-coding and with small effect size, which contribute to polygenic diseases. Here, we review the paradigm of UMOD, the gene coding for uromodulin, to illustrate how a kidney-specific protein of major physiological importance is involved in a spectrum of kidney disorders. This new field of investigation illustrates the importance of genetic variation in the pathogenesis and prognosis of disease, with therapeutic implications.

Aggregated Genomic Data as Cohort-Specific Allelic Frequencies can Boost Variants and Genes Prioritization in Non-Solved Cases of Inherited Retinal Dystrophies

The introduction of NGS in genetic diagnosis has increased the repertoire of variants and genes involved and the amount of genomic information produced. We built an allelic-frequency (AF) database for a heterogeneous cohort of genetic diseases to explore the aggregated genomic information and boost diagnosis in inherited retinal dystrophies (IRD). We retrospectively selected 5683 index-cases with clinical exome sequencing tests available, 1766 with IRD and the rest with diverse genetic diseases. We calculated a subcohort’s IRD-specific AF and compared it with suitable pseudocontrols. For non-solved IRD cases, we prioritized variants with a significant increment of frequencies, with eight variants that may help to explain the phenotype, and 10/11 of uncertain significance that were reclassified as probably pathogenic according to ACMG. Moreover, we developed a method to highlight genes with more frequent pathogenic variants in IRD cases than in pseudocontrols weighted by the increment of benign variants in the same comparison. We identified 18 genes for further studies that provided new insights in five cases. This resource can also help one to calculate the carrier frequency in IRD genes. A cohort-specific AF database assists with variants and genes prioritization and operates as an engine that provides a new hypothesis in non-solved cases, augmenting the diagnosis rate.

Specificity of Psychiatric Polygenic Risk Scores and their Effects on Associated Risk Phenotypes

Background

Polygenic risk scores (PRSs) are indices of genetic liability for illness, but their clinical utility for predicting risk for a specific psychiatric disorder is limited. Genetic overlap among disorders and their effects on allied phenotypes may be a possible explanation, but this has been difficult to quantify given focus on singular disorders and/or allied phenotypes.

Methods

We constructed PRSs for 5 psychiatric disorders (schizophrenia, bipolar disorder, major depressive disorder, autism spectrum disorder, attention-deficit/hyperactivity disorder) and 3 nonpsychiatric control traits (height, type II diabetes, irritable bowel disease) in the UK Biobank (N = 31,616) and quantified associations between PRSs and phenotypes allied with mental illness: behavioral (symptoms, cognition, trauma) and brain measures from magnetic resonance imaging. We then evaluated the extent of specificity among PRSs and their effects on these allied phenotypes.

Results

Correlations among psychiatric PRSs replicated previous work, with overlap between schizophrenia and bipolar disorder, which was distinct from overlap between autism spectrum disorder and attention-deficit/hyperactivity disorder; overlap between psychiatric and control PRSs was minimal. There was, however, substantial overlap of PRS effects on allied phenotypes among psychiatric disorders and among psychiatric disorders and control traits, where the extent and pattern of overlap was phenotype specific.

Conclusions

Results show that genetic distinctions between psychiatric disorders and between psychiatric disorders and control traits exist, but this does not extend to their effects on allied phenotypes. Although overlap can be informative, work is needed to construct PRSs that will function at the level of specificity needed for clinical application.

Inherited variants in CHD3 show variable expressivity in Snijders Blok-Campeau syndrome

PURPOSE

Common diagnostic next-generation sequencing strategies are not optimized to identify inherited variants in genes associated with dominant neurodevelopmental disorders as causal when the transmitting parent is clinically unaffected, leaving a significant number of cases with neurodevelopmental disorders undiagnosed.

METHODS

We characterized 21 families with inherited heterozygous missense or protein-truncating variants in CHD3, a gene in which de novo variants cause Snijders Blok-Campeau syndrome.

RESULTS

Computational facial and Human Phenotype Ontology-based comparisons showed that the phenotype of probands with inherited CHD3 variants overlaps with the phenotype previously associated with de novo CHD3 variants, whereas heterozygote parents are mildly or not affected, suggesting variable expressivity. In addition, similarly reduced expression levels of CHD3 protein in cells of an affected proband and of healthy family members with a CHD3 protein-truncating variant suggested that compensation of expression from the wild-type allele is unlikely to be an underlying mechanism. Notably, most inherited CHD3 variants were maternally transmitted.

CONCLUSION

Our results point to a significant role of inherited variation in Snijders Blok-Campeau syndrome, a finding that is critical for correct variant interpretation and genetic counseling and warrants further investigation toward understanding the broader contributions of such variation to the landscape of human disease.

Positron Emission Tomography in the Neuroimaging of Autism Spectrum Disorder: A Review

Autism spectrum disorder (ASD) is a basket term for neurodevelopmental disorders characterized by marked impairments in social interactions, repetitive and stereotypical behaviors, and restricted interests and activities. Subtypes include (A) disorders with known genetic abnormalities including fragile X syndrome, Rett syndrome, and tuberous sclerosis and (B) idiopathic ASD, conditions with unknown etiologies. Positron emission tomography (PET) is a molecular imaging technology that can be utilized in vivo for dynamic and quantitative research, and is a valuable tool for exploring pathophysiological mechanisms, evaluating therapeutic efficacy, and accelerating drug development in ASD. Recently, several imaging studies on ASD have been published and physiological changes during ASD progression was disclosed by PET. This paper reviews the specific radioligands for PET imaging of critical biomarkers in ASD, and summarizes and discusses the similar and different discoveries in outcomes of previous studies. It is of great importance to identify general physiological changes in cerebral glucose metabolism, cerebral blood flow perfusion, abnormalities in neurotransmitter systems, and inflammation in the central nervous system in ASD, which may provide excellent points for further ASD research.

Scaling up oligogenic diseases research with OLIDA: the Oligogenic Diseases Database

Improving the understanding of the oligogenic nature of diseases requires access to high-quality, well-curated Findable, Accessible, Interoperable, Reusable (FAIR) data. Although first steps were taken with the development of the Digenic Diseases Database, leading to novel computational advancements to assist the field, these were also linked with a number of limitations, for instance, the ad hoc curation protocol and the inclusion of only digenic cases. The OLIgogenic diseases DAtabase (OLIDA) presents a novel, transparent and rigorous curation protocol, introducing a confidence scoring mechanism for the published oligogenic literature. The application of this protocol on the oligogenic literature generated a new repository containing 916 oligogenic variant combinations linked to 159 distinct diseases. Information extracted from the scientific literature is supplemented with current knowledge support obtained from public databases. Each entry is an oligogenic combination linked to a disease, labelled with a confidence score based on the level of genetic and functional evidence that supports its involvement in this disease. These scores allow users to assess the relevance and proof of pathogenicity of each oligogenic combination in the database, constituting markers for reporting improvements on disease-causing oligogenic variant combinations. OLIDA follows the FAIR principles, providing detailed documentation, easy data access through its application programming interface and website, use of unique identifiers and links to existing ontologies.

Known allosteric proteins have central roles in genetic disease

Allostery is a form of protein regulation, where ligands that bind sites located apart from the active site can modify the activity of the protein. The molecular mechanisms of allostery have been extensively studied, because allosteric sites are less conserved than active sites, and drugs targeting them are more specific than drugs binding the active sites. Here we quantify the importance of allostery in genetic disease. We show that 1) known allosteric proteins are central in disease networks, contribute to genetic disease and comorbidities much more than non-allosteric proteins, and there is an association between being allosteric and involvement in disease; 2) they are enriched in many major disease types like hematopoietic diseases, cardiovascular diseases, cancers, diabetes, or diseases of the central nervous system; 3) variants from cancer genome-wide association studies are enriched near allosteric proteins, indicating their importance to polygenic traits; and 4) the importance of allosteric proteins in disease is due, at least partly, to their central positions in protein-protein interaction networks, and less due to their dynamical properties.

Patterning of the antero-ventral mammalian brain: Lessons from holoprosencephaly comparative biology in man and mouse

Adult form and function are dependent upon the activity of specialized signaling centers that act early in development at the embryonic midline. These centers instruct the surrounding cells to adopt a positional fate and to form the patterned structures of the phylotypic embryo. Abnormalities in these processes have devastating consequences for the individual, as exemplified by holoprosencephaly in which anterior midline development fails, leading to structural defects of the brain and/or face. In the 25 years since the first association between human holoprosencephaly and the sonic hedgehog gene, a combination of human and animal genetic studies have enhanced our understanding of the genetic and embryonic causation of this congenital defect. Comparative biology has extended the holoprosencephaly network via the inclusion of gene mutations from multiple signaling pathways known to be required for anterior midline formation. It has also clarified aspects of holoprosencephaly causation, showing that it arises when a deleterious variant is present within a permissive genome, and that environmental factors, as well as embryonic stochasticity, influence the phenotypic outcome of the variant. More than two decades of research can now be distilled into a framework of embryonic and genetic causation. This framework means we are poised to move beyond our current understanding of variants in signaling pathway molecules. The challenges now at the forefront of holoprosencephaly research include deciphering how the mutation of genes involved in basic cell processes can also cause holoprosencephaly, determining the important constituents of the holoprosencephaly permissive genome, and identifying environmental compounds that promote holoprosencephaly. This article is categorized under: Congenital Diseases > Stem Cells and Development Congenital Diseases > Genetics/Genomics/Epigenetics Congenital Diseases > Molecular and Cellular Physiology Congenital Diseases > Environmental Factors.

Contribution of Rare Genetic Variation to Disease Susceptibility in a Large Scandinavian Myositis Cohort

OBJECTIVE

Idiopathic inflammatory myopathies (IIMs) are a heterogeneous group of complex autoimmune conditions characterized by inflammation in skeletal muscle and extramuscular compartments, and interferon (IFN) system activation. We undertook this study to examine the contribution of genetic variation to disease susceptibility and to identify novel avenues for research in IIMs.

METHODS

Targeted DNA sequencing was used to mine coding and potentially regulatory single nucleotide variants from ~1,900 immune-related genes in a Scandinavian case-control cohort of 454 IIM patients and 1,024 healthy controls. Gene-based aggregate testing, together with rare variant- and gene-level enrichment analyses, was implemented to explore genotype-phenotype relations.

RESULTS

Gene-based aggregate tests of all variants, including rare variants, identified IFI35 as a potential genetic risk locus for IIMs, suggesting a genetic signature of type I IFN pathway activation. Functional annotation of the IFI35 locus highlighted a regulatory network linked to the skeletal muscle-specific gene PTGES3L, as a potential candidate for IIM pathogenesis. Aggregate genetic associations with AGER and PSMB8 in the major histocompatibility complex locus were detected in the antisynthetase syndrome subgroup, which also showed a less marked genetic signature of the type I IFN pathway. Enrichment analyses indicated a burden of synonymous and noncoding rare variants in IIM patients, suggesting increased disease predisposition associated with these classes of rare variants.

CONCLUSION

Our study suggests the contribution of rare genetic variation to disease susceptibility in IIM and specific patient subgroups, and pinpoints genetic associations consistent with previous findings by gene expression profiling. These features highlight genetic profiles that are potentially relevant to disease pathogenesis.

Stem Cell-Derived β Cells: A Versatile Research Platform to Interrogate the Genetic Basis of β Cell Dysfunction

Pancreatic β cell dysfunction is a central component of diabetes progression. During the last decades, the genetic basis of several monogenic forms of diabetes has been recognized. Genome-wide association studies (GWAS) have also facilitated the identification of common genetic variants associated with an increased risk of diabetes. These studies highlight the importance of impaired β cell function in all forms of diabetes. However, how most of these risk variants confer disease risk, remains unanswered. Understanding the specific contribution of genetic variants and the precise role of their molecular effectors is the next step toward developing treatments that target β cell dysfunction in the era of personalized medicine. Protocols that allow derivation of β cells from pluripotent stem cells, represent a powerful research tool that allows modeling of human development and versatile experimental designs that can be used to shed some light on diabetes pathophysiology. This article reviews different models to study the genetic basis of β cell dysfunction, focusing on the recent advances made possible by stem cell applications in the field of diabetes research.

Contribution of schizophrenia polygenic burden to longitudinal phenotypic variance in 22q11.2 deletion syndrome

While the recurrent 22q11.2 deletion is one of the strongest genetic risk factors for schizophrenia (SCZ), variability of its associated neuropsychiatric endophenotypes reflects its incomplete penetrance for psychosis development. To assess whether this phenotypic variability is linked to common variants associated with SCZ, we studied the association between SCZ polygenic risk score (PRS) and longitudinally acquired phenotypic information of the Swiss 22q11.2DS cohort (n = 97, 50% females, mean age 17.7 yr, mean visit interval 3.8 yr). The SCZ PRS with the best predictive performance was ascertained in the Estonian Biobank (n = 201,146) with LDpred. The infinitesimal SCZ PRS model showed the strongest capacity in discriminating SCZ cases from controls with one SD difference in SCZ PRS corresponding to an odds ratio (OR) of 1.73 (95% CI 1.57-1.90, P = 1.47 × 10^-29). In 22q11.2 patients, random-effects ordinal regression modelling using longitudinal data showed SCZ PRS to have the strongest effect on social anhedonia (OR = 2.09, P = 0.0002), and occupational functioning (OR = 1.82, P = 0.0003) within the negative symptoms course, and dysphoric mood (OR = 2.00, P = 0.002) and stress intolerance (OR = 1.76, P = 0.0002) within the general symptoms course. Genetic liability for SCZ was additionally associated with full scale cognitive decline (β = -0.25, P = 0.02) and with longitudinal volumetric reduction of the right and left hippocampi (β = -0.28, P = 0.005; β = -0.23, P = 0.02, respectively). Our results indicate that the polygenic contribution to SCZ acts upon the threshold-lowering first hit (i.e., the deletion). It modifies the endophenotypes of 22q11.2DS and augments the derailment of developmental trajectories of negative and general symptoms, cognition, and hippocampal volume.

Two Rare Variants in PLAU and BACE1 Genes-Do They Contribute to Semantic Dementia Clinical Phenotype?

We have performed whole-genome sequencing to identify the genetic variants potentially contributing to the early-onset semantic dementia phenotype in a patient with family history of dementia and episodic memory deficit accompanied with profound semantic loss. Only very rare variants of unknown significance (VUS) have been identified: a nonsense variant c.366C>A/p.Cys122* in plasminogen activator, urokinase (PLAU) and a missense variant c.944C>T/p.Thr315Met in β-site APP-cleaving enzyme 1 (BACE1)-along with known disease-modifying variants of moderate penetrance. Patient-derived fibroblasts showed reduced PLAU and elevated BACE1 mRNA and protein levels compared to control fibroblasts. Successful rescue of PLAU mRNA levels by nonsense-mediated mRNA decay (NMD) inhibitor (puromycin) confirmed NMD as the underlying mechanism. This is the first report of the PLAU variant with the confirmed haploinsufficiency, associated with semantic dementia phenotype. Our results suggest that rare variants in the PLAU and BACE1 genes should be considered in future studies on early-onset dementias.

[Genetic Analyses of Complex Phenotypes Through the Example of Anorexia Nervosa and Bodyweight Regulation]

Genetic Analyses of Complex Phenotypes Through the Example of Anorexia Nervosa and Bodyweight Regulation Abstract. Genetics variants are important for the regulation of bodyweight and also contribute to the genetic architecture of eating disorders. For many decades, family studies, a subentity of so-called formal genetic studies, were employed to determine the genetic share of bodyweight and eating disorders and found heritability rates exceeding 50 % with both phenotypes. Because of this significant contribution of genetics, the search for those genes and their variants related to the variance in bodyweight and the etiology of eating disorders - or both - was commenced by the early 1990s. Initially, candidate genes studies were conducted targeting those genes most plausibly related to either phenotype, especially based on pathophysiological considerations. This approach, however, implicated only a few genes in the regulation of bodyweight and did not provide significant insights into the genetics of eating disorders. Driven by considerable methodological advances in genetic research, especially related to the introduction of so-called genome-wide association studies by the beginning of the 21st century, today more than 1,000 variants/loci have been detected that affect the regulation of bodyweight. Eight such loci have been identified regarding anorexia nervosa (AN). These results as well as those from cross-disorder analyses provide insights into the complex regulation of bodyweight and demonstrated unforeseen pathomechanisms for AN.

Using Phecodes for Research with the Electronic Health Record: From PheWAS to PheRS

Electronic health records (EHRs) are a rich source of data for researchers, but extracting meaningful information out of this highly complex data source is challenging. Phecodes represent one strategy for defining phenotypes for research using EHR data. They are a high-throughput phenotyping tool based on ICD (International Classification of Diseases) codes that can be used to rapidly define the case/control status of thousands of clinically meaningful diseases and conditions. Phecodes were originally developed to conduct phenome-wide association studies to scan for phenotypic associations with common genetic variants. Since then, phecodes have been used to support a wide range of EHR-based phenotyping methods, including the phenotype risk score. This review aims to comprehensively describe the development, validation, and applications of phecodes and suggest some future directions for phecodes and high-throughput phenotyping.

Advancing clinical genomics and precision medicine with GVViZ: FAIR bioinformatics platform for variable gene-disease annotation, visualization, and expression analysis

Background

Genetic disposition is considered critical for identifying subjects at high risk for disease development. Investigating disease-causing and high and low expressed genes can support finding the root causes of uncertainties in patient care. However, independent and timely high-throughput next-generation sequencing data analysis is still a challenge for non-computational biologists and geneticists.

Results

In this manuscript, we present a findable, accessible, interactive, and reusable (FAIR) bioinformatics platform, i.e., GVViZ (visualizing genes with disease-causing variants). GVViZ is a user-friendly, cross-platform, and database application for RNA-seq-driven variable and complex gene-disease data annotation and expression analysis with a dynamic heat map visualization. GVViZ has the potential to find patterns across millions of features and extract actionable information, which can support the early detection of complex disorders and the development of new therapies for personalized patient care. The execution of GVViZ is based on a set of simple instructions that users without a computational background can follow to design and perform customized data analysis. It can assimilate patients’ transcriptomics data with the public, proprietary, and our in-house developed gene-disease databases to query, easily explore, and access information on gene annotation and classified disease phenotypes with greater visibility and customization. To test its performance and understand the clinical and scientific impact of GVViZ, we present GVViZ analysis for different chronic diseases and conditions, including Alzheimer’s disease, arthritis, asthma, diabetes mellitus, heart failure, hypertension, obesity, osteoporosis, and multiple cancer disorders. The results are visualized using GVViZ and can be exported as image (PNF/TIFF) and text (CSV) files that include gene names, Ensembl (ENSG) IDs, quantified abundances, expressed transcript lengths, and annotated oncology and non-oncology diseases.

Conclusions

We emphasize that automated and interactive visualization should be an indispensable component of modern RNA-seq analysis, which is currently not the case. However, experts in clinics and researchers in life sciences can use GVViZ to visualize and interpret the transcriptomics data, making it a powerful tool to study the dynamics of gene expression and regulation. Furthermore, with successful deployment in clinical settings, GVViZ has the potential to enable high-throughput correlations between patient diagnoses based on clinical and transcriptomics data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-021-00336-1.

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.

Pathogenic missense protein variants affect different functional pathways and proteomic features than healthy population variants

Missense variants are present amongst the healthy population, but some of them are causative of human diseases. A classification of variants associated with “healthy” or “diseased” states is therefore not always straightforward. A deeper understanding of the nature of missense variants in health and disease, the cellular processes they may affect, and the general molecular principles which underlie these differences is essential to offer mechanistic explanations of the true impact of pathogenic variants. Here, we have formalised a statistical framework which enables robust probabilistic quantification of variant enrichment across full-length proteins, their domains, and 3D structure-defined regions. Using this framework, we validate and extend previously reported trends of variant enrichment in different protein structural regions (surface/core/interface). By examining the association of variant enrichment with available functional pathways and transcriptomic and proteomic (protein half-life, thermal stability, abundance) data, we have mined a rich set of molecular features which distinguish between pathogenic and population variants: Pathogenic variants mainly affect proteins involved in cell proliferation and nucleotide processing and are enriched in more abundant proteins. Additionally, rare population variants display features closer to common than pathogenic variants. We validate the association between these molecular features and variant pathogenicity by comparing against existing in silico variant impact annotations. This study provides molecular details into how different proteins exhibit resilience and/or sensitivity towards missense variants and provides the rationale to prioritise variant-enriched proteins and protein domains for therapeutic targeting and development. The ZoomVar database, which we created for this study, is available at fraternalilab.kcl.ac.uk/ZoomVar. It allows users to programmatically annotate missense variants with protein structural information and to calculate variant enrichment in different protein structural regions.

How do can one improve the classification of genetic variants as harmful or harmless? This study uses a robust statistical analysis to exploit the interplay between protein structure, proteomic measurements and functional pathways to enable better discrimination between missense variants in health and disease.

Pathogenic Mutations and Putative Phenotype-Affecting Variants in Polish Myofibrillar Myopathy Patients

Myofibrillar myopathies (MFM) are heterogeneous hereditary muscle diseases with characteristic myopathological features of Z-disk dissolution and aggregates of its degradation products. The onset and progression of the disease are variable, with an elusive genetic background, and around half of the cases lacking molecular diagnosis. Here, we attempted to establish possible genetic foundations of MFM by performing whole exome sequencing (WES) in eleven unrelated families of 13 patients clinically diagnosed as MFM spectrum. A filtering strategy aimed at identification of variants related to the disease was used and included integrative analysis of WES data and human phenotype ontology (HPO) terms, analysis of muscle-expressed genes, and analysis of the disease-associated interactome. Genetic diagnosis was possible in eight out of eleven cases. Putative causative mutations were found in the DES (two cases), CRYAB, TPM3, and SELENON (four cases) genes, the latter typically presenting with a rigid spine syndrome. Moreover, a variety of additional, possibly phenotype-affecting variants were found. These findings indicate a markedly heterogeneous genetic background of MFM and show the usefulness of next generation sequencing in the identification of disease-associated mutations. Finally, we discuss the emerging concept of variant load as the basis of phenotypic heterogeneity.

The Use of Whole Genome and Exome Sequencing for Newborn Screening: Challenges and Opportunities for Population Health

Newborn screening (NBS) is a population-based program with a goal of reducing the burden of disease for conditions with significant clinical impact on neonates. Screening tests were originally developed and implemented one at a time, but newer methods have allowed the use of multiplex technologies to expand additions more rapidly to standard panels. Recent improvements in next-generation sequencing are also evolving rapidly from first focusing on individual genes, then panels, and finally all genes as encompassed by whole exome and genome sequencing. The intersection of these two technologies brings the revolutionary possibility of identifying all genetic disorders in newborns, allowing implementation of therapies at the optimum time regardless of symptoms. This article reviews the history of newborn screening and early studies examining the use of whole genome and exome sequencing as a screening tool. Lessons learned from these studies are discussed, along with technical, ethical, and societal challenges to broad implementation.

Improving the informativeness of Mendelian disease-derived pathogenicity scores for common disease

Despite considerable progress on pathogenicity scores prioritizing variants for Mendelian disease, little is known about the utility of these scores for common disease. Here, we assess the informativeness of Mendelian disease-derived pathogenicity scores for common disease and improve upon existing scores. We first apply stratified linkage disequilibrium (LD) score regression to evaluate published pathogenicity scores across 41 common diseases and complex traits (average N = 320K). Several of the resulting annotations are informative for common disease, even after conditioning on a broad set of functional annotations. We then improve upon published pathogenicity scores by developing AnnotBoost, a machine learning framework to impute and denoise pathogenicity scores using a broad set of functional annotations. AnnotBoost substantially increases the informativeness for common disease of both previously uninformative and previously informative pathogenicity scores, implying that Mendelian and common disease variants share similar properties. The boosted scores also produce improvements in heritability model fit and in classifying disease-associated, fine-mapped SNPs. Our boosted scores may improve fine-mapping and candidate gene discovery for common disease.

Rare genetic causes of complex kidney and urological diseases

Although often considered a single-entity, chronic kidney disease (CKD) comprises many pathophysiologically distinct disorders that result in persistently abnormal kidney structure and/or function, and encompass both monogenic and polygenic aetiologies. Rare inherited forms of CKD frequently span diverse phenotypes, reflecting genetic phenomena including pleiotropy, incomplete penetrance and variable expressivity. Use of chromosomal microarray and massively parallel sequencing technologies has revealed that genomic disorders and monogenic aetiologies contribute meaningfully to seemingly complex forms of CKD across different clinically defined subgroups and are characterized by high genetic and phenotypic heterogeneity. Investigations of prevalent genomic disorders in CKD have integrated genetic, bioinformatic and functional studies to pinpoint the genetic drivers underlying their renal and extra-renal manifestations, revealing both monogenic and polygenic mechanisms. Similarly, massively parallel sequencing-based analyses have identified gene- and allele-level variation that contribute to the clinically diverse phenotypes observed for many monogenic forms of nephropathy. Genome-wide sequencing studies suggest that dual genetic diagnoses are found in at least 5% of patients in whom a genetic cause of disease is identified, highlighting the fact that complex phenotypes can also arise from multilocus variation. A multifaceted approach that incorporates genetic and phenotypic data from large, diverse cohorts will help to elucidate the complex relationships between genotype and phenotype for different forms of CKD, supporting personalized medicine for individuals with kidney disease.

Practicing precision medicine with intelligently integrative clinical and multi-omics data analysis

Precision medicine aims to empower clinicians to predict the most appropriate course of action for patients with complex diseases like cancer, diabetes, cardiomyopathy, and COVID-19. With a progressive interpretation of the clinical, molecular, and genomic factors at play in diseases, more effective and personalized medical treatments are anticipated for many disorders. Understanding patient’s metabolomics and genetic make-up in conjunction with clinical data will significantly lead to determining predisposition, diagnostic, prognostic, and predictive biomarkers and paths ultimately providing optimal and personalized care for diverse, and targeted chronic and acute diseases. In clinical settings, we need to timely model clinical and multi-omics data to find statistical patterns across millions of features to identify underlying biologic pathways, modifiable risk factors, and actionable information that support early detection and prevention of complex disorders, and development of new therapies for better patient care. It is important to calculate quantitative phenotype measurements, evaluate variants in unique genes and interpret using ACMG guidelines, find frequency of pathogenic and likely pathogenic variants without disease indicators, and observe autosomal recessive carriers with a phenotype manifestation in metabolome. Next, ensuring security to reconcile noise, we need to build and train machine-learning prognostic models to meaningfully process multisource heterogeneous data to identify high-risk rare variants and make medically relevant predictions. The goal, today, is to facilitate implementation of mainstream precision medicine to improve the traditional symptom-driven practice of medicine, and allow earlier interventions using predictive diagnostics and tailoring better-personalized treatments. We strongly recommend automated implementation of cutting-edge technologies, utilizing machine learning (ML) and artificial intelligence (AI) approaches for the multimodal data aggregation, multifactor examination, development of knowledgebase of clinical predictors for decision support, and best strategies for dealing with relevant ethical issues.

ZMAT2 in Humans and Other Primates: A Highly Conserved and Understudied Gene

Recent advances in genetics present unique opportunities for enhancing our understanding of human physiology and disease predisposition through detailed analysis of gene structure, expression, and population variation via examination of data in publicly accessible genome and gene expression repositories. Yet, the vast majority of human genes remain understudied. Here, we show the scope of these genomic and genetic resources by evaluating ZMAT2, a member of a 5-gene family that through May 2020 had been the focus of only 4 peer-reviewed scientific publications. Using analysis of information extracted from public databases, we show that human ZMAT2 is a 6-exon gene and find that it exhibits minimal genetic variation in human populations and in disease states, including cancer. We further demonstrate that the gene and its encoded protein are highly conserved among nonhuman primates and define a cohort of ZMAT2 pseudogenes in the marmoset genome. Collectively, our investigations illustrate how complementary use of genomic, gene expression, and population genetic resources can lead to new insights about human and mammalian biology and evolution, and when coupled with data supporting key roles for ZMAT2 in keratinocyte differentiation and pre-RNA splicing argue that this gene is worthy of further study.

Genetic modifiers and phenotypic variability in neuromuscular disorders

Neuromuscular disorders are mostly rare diseases with autosomal dominant, recessive, or X-linked inheritance. Interestingly, among patients carrying the same mutations, a range of phenotypic severity is reported. This phenotypic variability in neuromuscular disorders is still not fully understood. This review will focus on genetic modifiers and will briefly describe metabolic pathways, in which they are involved. Genetic modifiers are variants in the same or other genes that modulate the phenotype. Proteins encoded by genetic modifiers in neuromuscular diseases are taking part in different metabolic processes, most commonly in inflammation, growth and regeneration, endoplasmic reticulum metabolism, and cytoskeletal activities. Recent advances in omics technologies, development of computational algorithms, and establishing large international consortia intensified discovery sped up investigation of genetic modifiers. As more individuals affected by neuromuscular disorders are tested, it is often suggested that classic models of genetic causation cannot explain phenotypic variability. There is a growing interest in their discovery and identifying shared metabolic pathways can contribute to design targeted therapies. We provide an update on variants acting as genetic modifiers in neuromuscular disorders and strategies used for their discovery.

Structure and expression of the long noncoding RNA gene MIR503 in humans and non-human primates

Recent technical and other advances in genomics provide unique opportunities to improve our understanding of human physiology and disease predisposition through a detailed analysis of gene structure and expression by examining data in public genome and gene-expression repositories. Yet, the vast majority of human genes remain understudied. This is particularly true of genes for long noncoding RNAs (lncRNAs). Here, we describe the detailed characterization of MIR503HG, a lncRNA gene found on the X chromosome in humans. Using information extracted from public databases, we show that human MIR503HG is a 5-exon gene, and that it is highly conserved among 5 non-human primates spanning over 85 million years ago of evolutionary diversification. MIR503HG is transcribed and processed into multiple distinct RNAs in each of these species through differential exon use and alternative RNA splicing, with a higher abundance of transcripts being found in reproductive tissues, especially during the early stages of ovary and testis development, indicating a possible role in reproductive biology. Furthermore, in select reproductive system cancers, MIR503HG transcripts are downregulated, with higher levels of RNA expression being associated with clinical outcomes. Collectively, these investigations show how the use of genomic, gene expression, and other genetic resources can lead to new insights about human biology and disease, and argue that MIR503HG is worthy of additional study.

Exome Sequencing Identifies Genes and Gene Sets Contributing to Severe Childhood Obesity, Linking PHIP Variants to Repressed POMC Transcription

Obesity is genetically heterogeneous with monogenic and complex polygenic forms. Using exome and targeted sequencing in 2,737 severely obese cases and 6,704 controls, we identified three genes (PHIP, DGKI, and ZMYM4) with an excess burden of very rare predicted deleterious variants in cases. In cells, we found that nuclear PHIP (pleckstrin homology domain interacting protein) directly enhances transcription of pro-opiomelanocortin (POMC), a neuropeptide that suppresses appetite. Obesity-associated PHIP variants repressed POMC transcription. Our demonstration that PHIP is involved in human energy homeostasis through transcriptional regulation of central melanocortin signaling has potential diagnostic and therapeutic implications for patients with obesity and developmental delay. Additionally, we found an excess burden of predicted deleterious variants involving genes nearest to loci from obesity genome-wide association studies. Genes and gene sets influencing obesity with variable penetrance provide compelling evidence for a continuum of causality in the genetic architecture of obesity, and explain some of its missing heritability.

A Bayesian method to estimate variant-induced disease penetrance

A major challenge emerging in genomic medicine is how to assess best disease risk from rare or novel variants found in disease-related genes. The expanding volume of data generated by very large phenotyping efforts coupled to DNA sequence data presents an opportunity to reinterpret genetic liability of disease risk. Here we propose a framework to estimate the probability of disease given the presence of a genetic variant conditioned on features of that variant. We refer to this as the penetrance, the fraction of all variant heterozygotes that will present with disease. We demonstrate this methodology using a well-established disease-gene pair, the cardiac sodium channel gene SCN5A and the heart arrhythmia Brugada syndrome. From a review of 756 publications, we developed a pattern mixture algorithm, based on a Bayesian Beta-Binomial model, to generate SCN5A penetrance probabilities for the Brugada syndrome conditioned on variant-specific attributes. These probabilities are determined from variant-specific features (e.g. function, structural context, and sequence conservation) and from observations of affected and unaffected heterozygotes. Variant functional perturbation and structural context prove most predictive of Brugada syndrome penetrance.

The Zmat2 gene in non-mammalian vertebrates: Organizational simplicity within a divergent locus in fish

ZMAT2 is among the least-studied of mammalian proteins and genes, even though it is the ortholog of Snu23, a protein involved in pre-mRNA splicing in yeast. Here we have used data from genomic and gene expression repositories to examine the Zmat2 gene and locus in 8 terrestrial vertebrates, 10 ray-finned fish, and 1 lobe-finned fish representing > 500 million years of evolutionary diversification. The analyses revealed that vertebrate Zmat2 genes are similar to their mammalian counterparts, as in 16/19 species studied they contain 6 exons, and in 18/19 encode a single conserved protein. However, unlike in mammals, no Zmat2 pseudogenes were identified in these vertebrates, although an expressed Zmat2 paralog was characterized in flycatcher that resembled a DNA copy of a processed and retro-transposed mRNA, and thus could be a proto-pseudogene captured during its evolutionary journey from active to inert. The Zmat2 locus in terrestrial vertebrates, and in spotted gar and coelacanth, also shares additional genes with its mammalian counterparts, including Histidyl-tRNA synthetase (Hars), Hars2, and others, but these are absent from the Zmat2 locus in teleost fish, in which Stem-loop-binding protein 2 (Slbp2) and Lymphocyte cytosolic protein 2a (Lcp2a) are present instead. Taken together, these observations argue that a recognizable Zmat2 was present in the earliest vertebrate ancestors, and postulate that during chromosomal tetraploidization and subsequent re-diploidization during modern teleost evolution, the duplicated Zmat2 gene was retained and the original lost. This study also highlights how information from genomic resources can be leveraged to reveal new biologically significant insights.

ORVAL: a novel platform for the prediction and exploration of disease-causing oligogenic variant combinations

A tremendous amount of DNA sequencing data is being produced around the world with the ambition to capture in more detail the mechanisms underlying human diseases. While numerous bioinformatics tools exist that allow the discovery of causal variants in Mendelian diseases, little to no support is provided to do the same for variant combinations, an essential task for the discovery of the causes of oligogenic diseases. ORVAL (the Oligogenic Resource for Variant AnaLysis), which is presented here, provides an answer to this problem by focusing on generating networks of candidate pathogenic variant combinations in gene pairs, as opposed to isolated variants in unique genes. This online platform integrates innovative machine learning methods for combinatorial variant pathogenicity prediction with visualization techniques, offering several interactive and exploratory tools, such as pathogenic gene and protein interaction networks, a ranking of pathogenic gene pairs, as well as visual mappings of the cellular location and pathway information. ORVAL is the first web-based exploration platform dedicated to identifying networks of candidate pathogenic variant combinations with the sole ambition to help in uncovering oligogenic causes for patients that cannot rely on the classical disease analysis tools. ORVAL is available at https://orval.ibsquare.be.