Human genomics. Effect of predicted protein-truncating genetic variants on the human transcriptome

Genomics of periodontal disease and tooth morbidity

In this review we critically summarize the evidence base and the progress to date regarding the genomic basis of periodontal disease and tooth morbidity (ie, dental caries and tooth loss), and discuss future applications and research directions in the context of precision oral health and care. Evidence for these oral/dental traits from genome-wide association studies first emerged less than a decade ago. Basic and translational research activities in this domain are now under way by multiple groups around the world. Key departure points in the oral health genomics discourse are: (a) some heritable variation exists for periodontal and dental diseases; (b) the environmental component (eg, social determinants of health and behavioral risk factors) has a major influence on the population distribution but probably interacts with factors of innate susceptibility at the person-level; (c) sizeable, multi-ethnic, well-characterized samples or cohorts with high-quality measures on oral health outcomes and genomics information are required to make decisive discoveries; (d) challenges remain in the measurement of oral health and disease, with current periodontitis and dental caries traits capturing only a part of the health-disease continuum, and are little or not informed by the underlying biology; (e) the substantial individual heterogeneity that exists in the clinical presentation and lifetime trajectory of oral disease can be identified and leveraged in a precision medicine framework or, if unappreciated, can hamper translational efforts. In this review we discuss how composite or biologically informed traits may offer improvements over clinically defined ones for the genomic interrogation of oral diseases. We demonstrate the utility of the results of genome-wide association studies for the development and testing of a genetic risk score for severe periodontitis. We conclude that exciting opportunities lie ahead for improvements in the oral health of individual patients and populations via advances in our understanding of the genomic basis of oral health and disease. The pace of new discoveries and their equitable translation to practice will largely depend on investments in the education and training of the oral health care workforce, basic and population research, and sustained collaborative efforts..

Cancer biology as revealed by the research autopsy

A research autopsy is a post-mortem medical procedure performed on a deceased individual with the primary goal of collecting tissue to support basic and translational research. This approach has increasingly been used to investigate the pathophysiological mechanisms of cancer evolution, metastasis and treatment resistance. In this Review, we discuss the rationale for the use of research autopsies in cancer research and provide an evidence-based discussion of the quality of post-mortem tissues compared with other types of biospecimens. We also discuss the advantages of using post-mortem tissues over other types of biospecimens, including the large amounts of tissue that can be obtained and the extent of multiregion sampling that is achievable, which is not otherwise possible in living patients. We highlight how the research autopsy has supported the identification of the clonal origins and modes of spread among metastases, the extent that selective pressures imposed by treatments cause bottlenecks leading to parallel and convergent tumour evolution, and the creation of rare tissue banks and patient-derived model systems. Finally, we comment on the future of the research autopsy as an integral component of precision medicine strategies.

The impact of nonsense-mediated mRNA decay on genetic disease, gene editing and cancer immunotherapy

Premature termination codons (PTCs) can result in the production of truncated proteins or the degradation of messenger RNAs by nonsense-mediated mRNA decay (NMD). Which of these outcomes occurs can alter the effect of a mutation, with the engagement of NMD being dependent on a series of rules. Here, by applying these rules genome-wide to obtain a resource called NMDetective, we explore the impact of NMD on genetic disease and approaches to therapy. First, human genetic diseases differ in whether NMD typically aggravates or alleviates the effects of PTCs. Second, failure to trigger NMD is a cause of ineffective gene inactivation by CRISPR-Cas9 gene editing. Finally, NMD is a determinant of the efficacy of cancer immunotherapy, with only frameshifted transcripts that escape NMD predicting a response. These results demonstrate the importance of incorporating the rules of NMD into clinical decision-making. Moreover, they suggest that inhibiting NMD may be effective in enhancing cancer immunotherapy.

DNA Variant in the RPGRIP1L Gene Influences Alternative Splicing

The retinitis pigmentosa GTPase regulator interacting protein 1-like (RPGRIP1L) gene encodes a ciliary protein that is critical for processes related to brain development, including development of left-right asymmetry, sonic hedgehog signaling, and neural tube formation. RPGRIP1L is a risk factor for retinal degeneration, and rare, deleterious variants in the RPGRIP1L gene cause Joubert syndrome and Meckel syndrome, both autosomal recessive disorders. These syndromes are characterized by dysfunctional primary cilia that result in abnormal development - and even lethality in the case of Meckel syndrome. Genetic studies have also implicated RPGRIP1L in psychiatric disorders by suggestive findings from genome-wide association studies and findings from rare-variant exome analyses for bipolar disorder and de novo mutations in autism. In this study we identify a common variant in RPGRIP1L, rs7203525, that influences alternative splicing, increasing the inclusion of exon 20 of RPGRIP1L. We detected this alternative splicing association in human postmortem brain tissue samples and, using a minigene assay combined with in vitro mutagenesis, confirmed that the alternative splicing is attributable to the alleles of this variant. The predominate RPGRIP1L isoform expressed in adult brains does not contain exon 20; thus, a shift to include this exon may impact brain function.

Components of genetic associations across 2,138 phenotypes in the UK Biobank highlight adipocyte biology

Population-based biobanks with genomic and dense phenotype data provide opportunities for generating effective therapeutic hypotheses and understanding the genomic role in disease predisposition. To characterize latent components of genetic associations, we apply truncated singular value decomposition (DeGAs) to matrices of summary statistics derived from genome-wide association analyses across 2,138 phenotypes measured in 337,199 White British individuals in the UK Biobank study. We systematically identify key components of genetic associations and the contributions of variants, genes, and phenotypes to each component. As an illustration of the utility of the approach to inform downstream experiments, we report putative loss of function variants, rs114285050 (GPR151) and rs150090666 (PDE3B), that substantially contribute to obesity-related traits and experimentally demonstrate the role of these genes in adipocyte biology. Our approach to dissect components of genetic associations across the human phenome will accelerate biomedical hypothesis generation by providing insights on previously unexplored latent structures.

Systematic Analysis of Splice-Site-Creating Mutations in Cancer

For the past decade, cancer genomic studies have focused on mutations leading to splice-site disruption, overlooking those having splice-creating potential. Here, we applied a bioinformatic tool, MiSplice, for the large-scale discovery of splice-site-creating mutations (SCMs) across 8,656 TCGA tumors. We report 1,964 originally mis-annotated mutations having clear evidence of creating alternative splice junctions. TP53 and GATA3 have 26 and 18 SCMs, respectively, and ATRX has 5 from lower-grade gliomas. Mutations in 11 genes, including PARP1, BRCA1, and BAP1, were experimentally validated for splice-site-creating function. Notably, we found that neoantigens induced by SCMs are likely several folds more immunogenic compared to missense mutations, exemplified by the recurrent GATA3 SCM. Further, high expression of PD-1 and PD-L1 was observed in tumors with SCMs, suggesting candidates for immune blockade therapy. Our work highlights the importance of integrating DNA and RNA data for understanding the functional and the clinical implications of mutations in human diseases.

Hotspot SF3B1 mutations induce metabolic reprogramming and vulnerability to serine deprivation

Cancer-associated mutations in the spliceosome gene SF3B1 create a neomorphic protein that produces aberrant mRNA splicing in hundreds of genes, but the ensuing biologic and therapeutic consequences of this missplicing are not well understood. Here we have provided evidence that aberrant splicing by mutant SF3B1 altered the transcriptome, proteome, and metabolome of human cells, leading to missplicing-associated downregulation of metabolic genes, decreased mitochondrial respiration, and suppression of the serine synthesis pathway. We also found that mutant SF3B1 induces vulnerability to deprivation of the nonessential amino acid serine, which was mediated by missplicing-associated downregulation of the serine synthesis pathway enzyme PHGDH. This vulnerability was manifest both in vitro and in vivo, as dietary restriction of serine and glycine in mice was able to inhibit the growth of SF3B1MUT xenografts. These findings describe a role for SF3B1 mutations in altered energy metabolism, and they offer a new therapeutic strategy against SF3B1MUT cancers.

Associations of CDH1 germline variant location and cancer phenotype in families with hereditary diffuse gastric cancer (HDGC)

INTRODUCTION

Hereditary diffuse gastric cancer (HDGC) is a cancer syndrome associated with variants in E-cadherin (CDH1), diffuse gastric cancer and lobular breast cancer. There is considerable heterogeneity in its clinical manifestations. This study aimed to determine associations between CDH1 germline variant status and clinical phenotypes of HDGC.

METHODS

One hundred and fifty-two HDGC families, including six previously unreported families, were identified. CDH1 gene-specific guidelines released by the Clinical Genome Resource (ClinGen) CDH1 Variant Curation Expert Panel were applied for pathogenicity classification of truncating, missense and splice site CDH1 germline variants. We evaluated ORs between location of truncating variants of CDH1 and incidence of colorectal cancer, breast cancer and cancer at young age (gastric cancer at <40 or breast cancer <50 years of age).

RESULTS

Frequency of truncating germline CDH1 variants varied across functional domains of the E-cadherin receptor gene and was highest in linker (0.05785 counts/base pair; p=0.0111) and PRE regions (0.10000; p=0.0059). Families with truncating CDH1 germline variants located in the PRE-PRO region were six times more likely to have family members affected by colorectal cancer (OR 6.20, 95% CI 1.79 to 21.48; p=0.004) compared with germline variants in other regions. Variants in the intracellular E-cadherin region were protective for cancer at young age (OR 0.2, 95% CI 0.06 to 0.64; p=0.0071) and in the linker regions for breast cancer (OR 0.35, 95% CI 0.12 to 0.99; p=0.0493). Different CDH1 genotypes were associated with different intracellular signalling activation levels including different p-ERK, p-mTOR and β-catenin levels in early submucosal T1a lesions of HDGC families with different CDH1 variants.

CONCLUSION

Type and location of CDH1 germline variants may help to identify families at increased risk for concomitant cancers that might benefit from individualised surveillance and intervention strategies.

A patient with homozygous nonsense variants in two Leigh syndrome disease genes: Distinguishing a dual diagnosis from a hypomorphic protein-truncating variant

Leigh syndrome is a mitochondrial disease caused by pathogenic variants in over 85 genes. Whole exome sequencing of a patient with Leigh-like syndrome identified homozygous protein-truncating variants in two genes associated with Leigh syndrome; a reported pathogenic variant in PDHX (NP_003468.2:p.(Arg446*)), and an uncharacterized variant in complex I (CI) assembly factor TIMMDC1 (NP_057673.2:p.(Arg225*)). The TIMMDC1 variant was predicted to truncate 61 amino acids at the C-terminus and functional studies demonstrated a hypomorphic impact of the variant on CI assembly. However, the mutant protein could still rescue CI assembly in TIMMDC1 knockout cells and the patient's clinical phenotype was not clearly distinct from that of other patients with the same PDHX defect. Our data suggest that the hypomorphic effect of the TIMMDC1 protein-truncating variant does not constitute a dual diagnosis in this individual. We recommend cautious assessment of variants in the C-terminus of TIMMDC1 and emphasize the need to consider the caveats detailed within the American College of Medical Genetics and Genomics (ACMG) criteria when assessing variants.

Adjusting for Principal Components of Molecular Phenotypes Induces Replicating False Positives

High-throughput measurements of molecular phenotypes provide an unprecedented opportunity to model cellular processes and their impact on disease. These highly structured datasets are usually strongly confounded, creating false positives and reducing power. This has motivated many approaches based on principal components analysis (PCA) to estimate and correct for confounders, which have become indispensable elements of association tests between molecular phenotypes and both genetic and nongenetic factors. Here, we show that these correction approaches induce a bias, and that it persists for large sample sizes and replicates out-of-sample. We prove this theoretically for PCA by deriving an analytic, deterministic, and intuitive bias approximation. We assess other methods with realistic simulations, which show that perturbing any of several basic parameters can cause false positive rate (FPR) inflation. Our experiments show the bias depends on covariate and confounder sparsity, effect sizes, and their correlation. Surprisingly, when the covariate and confounder have [Formula: see text], standard two-step methods all have [Formula: see text]-fold FPR inflation. Our analysis informs best practices for confounder correction in genomic studies, and suggests many false discoveries have been made and replicated in some differential expression analyses.

Genomic Analysis in the Age of Human Genome Sequencing

Affordable genome sequencing technologies promise to revolutionize the field of human genetics by enabling comprehensive studies that interrogate all classes of genome variation, genome-wide, across the entire allele frequency spectrum. Ongoing projects worldwide are sequencing many thousands-and soon millions-of human genomes as part of various gene mapping studies, biobanking efforts, and clinical programs. However, while genome sequencing data production has become routine, genome analysis and interpretation remain challenging endeavors with many limitations and caveats. Here, we review the current state of technologies for genetic variant discovery, genotyping, and functional interpretation and discuss the prospects for future advances. We focus on germline variants discovered by whole-genome sequencing, genome-wide functional genomic approaches for predicting and measuring variant functional effects, and implications for studies of common and rare human disease.

Genome, transcriptome and proteome: the rise of omics data and their integration in biomedical sciences

Advances in the technologies and informatics used to generate and process large biological data sets (omics data) are promoting a critical shift in the study of biomedical sciences. While genomics, transcriptomics and proteinomics, coupled with bioinformatics and biostatistics, are gaining momentum, they are still, for the most part, assessed individually with distinct approaches generating monothematic rather than integrated knowledge. As other areas of biomedical sciences, including metabolomics, epigenomics and pharmacogenomics, are moving towards the omics scale, we are witnessing the rise of inter-disciplinary data integration strategies to support a better understanding of biological systems and eventually the development of successful precision medicine. This review cuts across the boundaries between genomics, transcriptomics and proteomics, summarizing how omics data are generated, analysed and shared, and provides an overview of the current strengths and weaknesses of this global approach. This work intends to target students and researchers seeking knowledge outside of their field of expertise and fosters a leap from the reductionist to the global-integrative analytical approach in research.

The Continuing Evolution of Precision Health in Type 2 Diabetes: Achievements and Challenges

PURPOSE OF REVIEW

The purpose of this review was to summarize recent advances in the genomics of type 2 diabetes (T2D) and to highlight current initiatives to advance precision health.

RECENT FINDINGS

Generation of multi-omic data to measure each of the "biologic layers," developments in describing genomic function and annotation in T2D relevant tissue, along with the increasing recognition that T2D is a heterogeneous disease, and large-scale collaborations have all contributed to advancing our understanding of the molecular basis of T2D. Substantial advances have been made in understanding the molecular basis of T2D pathogenesis, such that precision health diabetes is increasingly becoming a reality. For precision diabetes to become a routine in clinical and public health, additional large-scale multi-omic initiatives are needed along with better assessment of our environment to delineate an individual's diabetes subtype for improved detection and management.

Protocols, Methods, and Tools for Genome-Wide Association Studies (GWAS) of Dental Traits

Oral health and disease are known to be influenced by complex interactions between environmental (e.g., social and behavioral) factors and innate susceptibility. Although the exact contribution of genomics and other layers of "omics" to oral health is an area of active research, it is well established that the susceptibility to dental caries, periodontal disease, and other oral and craniofacial traits is substantially influenced by the human genome. A comprehensive understanding of these genomic factors is necessary for the realization of precision medicine in the oral health domain. To aid in this direction, the advent and increasing affordability of high-throughput genotyping has enabled the simultaneous interrogation of millions of genetic polymorphisms for association with oral and craniofacial traits. Specifically, genome-wide association studies (GWAS) of dental caries and periodontal disease have provided initial insights into novel loci and biological processes plausibly implicated in these two common, complex, biofilm-mediated diseases. This paper presents a summary of protocols, methods, tools, and pipelines for the conduct of GWAS of dental caries, periodontal disease, and related traits. The protocol begins with the consideration of different traits for both diseases and outlines procedures for genotyping, quality control, adjustment for population stratification, heritability and association analyses, annotation, reporting, and interpretation. Methods and tools available for GWAS are being constantly updated and improved; with this in mind, the presented approaches have been successfully applied in numerous GWAS and meta-analyses among tens of thousands of individuals, including dental traits such as dental caries and periodontal disease. As such, they can serve as a guide or template for future genomic investigations of these and other traits.

Breast cancer risk in neurofibromatosis type 1 is a function of the type of NF1 gene mutation: a new genotype-phenotype correlation

Background

Neurofibromatosis type 1 (NF1) predisposes to breast cancer (BC), but no genotype-phenotype correlations have been described.

Methods

Constitutional NF1 mutations in 78 patients with NF1 with BC (NF1-BC) were compared with the NF1 Leiden Open Variation Database (n=3432).

Results

No cases were observed with whole or partial gene deletions (HR 0.10; 95% CI 0.006 to 1.63; p=0.014, Fisher’s exact test). There were no gross relationships with mutation position. Forty-five (64.3%; HR 6.4–83) of the 70 different mutations were more frequent than expected (p<0.05), while 52 (74.3%; HR 5.3–83) were significant when adjusted for multiple comparisons (adjusted p≤0.125; Benjamini-Hochberg). Higher proportions of both nonsense and missense mutations were also observed (adjusted p=0.254; Benjamini-Hochberg). Ten of the 11 missense cases with known age of BC occurred at <50 years (p=0.041). Eighteen cases had BRCA1/2 testing, revealing one BRCA2 mutation.

Discussion

These data strongly support the hypothesis that certain constitutional mutation types, and indeed certain specific variants in NF1 confer different risks of BC. The lack of large deletions and excess of nonsenses and missenses is consistent with gain of function mutations conferring risk of BC, and also that neurofibromin may function as a dimer. The observation that somatic NF1 amplification can occur independently of ERBB2 amplification in sporadic BC supports this concept. A prospective clinical-molecular study of NF1-BC needs to be established to confirm and build on these findings, but regardless of NF1 mutation status patients with NF1-BC warrant testing of other BC-predisposing genes.

Allelic imbalance and haploinsufficiency in MYBPC3-linked hypertrophic cardiomyopathy

Mutations in cardiac myosin binding protein C (MYBPC3) represent the most frequent cause of familial hypertrophic cardiomyopathy (HCM), making up approximately 50% of identified HCM mutations. MYBPC3 is distinct among other sarcomere genes associated with HCM in that truncating mutations make up the vast majority, whereas nontruncating mutations predominant in other sarcomere genes. Several studies using myocardial tissue from HCM patients have found reduced abundance of wild-type MYBPC3 compared to control hearts, suggesting haploinsufficiency of full-length MYBPC3. Further, decreased mutant versus wild-type mRNA and lack of truncated mutant MYBPC3 protein has been demonstrated, highlighting the presence of allelic imbalance. In this review, we will begin by introducing allelic imbalance and haploinsufficiency, highlighting the broad role each plays within the spectrum of human disease. We will subsequently focus on the roles allelic imbalance and haploinsufficiency play within MYBPC3-linked HCM. Finally, we will explore the implications of these findings on future directions of HCM research. An improved understanding of allelic imbalance and haploinsufficiency may help us better understand genotype-phenotype relationships in HCM and develop novel targeted therapies, providing exciting future research opportunities.

Genetic analysis of isoform usage in the human anti-viral response reveals influenza-specific regulation of ERAP2 transcripts under balancing selection

While genetic variants are known to be associated with overall gene abundance in stimulated immune cells, less is known about their effects on alternative isoform usage. By analyzing RNA-seq profiles of monocyte-derived dendritic cells from 243 individuals, we uncovered thousands of unannotated isoforms synthesized in response to influenza infection and type 1 interferon stimulation. We identified more than a thousand quantitative trait loci (QTLs) associated with alternate isoform usage (isoQTLs), many of which are independent of expression QTLs (eQTLs) for the same gene. Compared with eQTLs, isoQTLs are enriched for splice sites and untranslated regions, but depleted of sequences upstream of annotated transcription start sites. Both eQTLs and isoQTLs explain a significant proportion of the disease heritability attributed to common genetic variants. At the ERAP2 locus, we shed light on the function of the gene and how two frequent, highly differentiated haplotypes with intermediate frequencies could be maintained by balancing selection. At baseline and following type 1 interferon stimulation, the major haplotype is associated with low ERAP2 expression caused by nonsense-mediated decay, while the minor haplotype, known to increase Crohn's disease risk, is associated with high ERAP2 expression. In response to influenza infection, we found two uncharacterized isoforms expressed from the major haplotype, likely the result of multiple perfectly linked variants affecting the transcription and splicing at the locus. Thus, genetic variants at a single locus could modulate independent gene regulatory processes in innate immune responses and, in the case of ERAP2, may confer a historical fitness advantage in response to virus.

Specifications of the ACMG/AMP variant curation guidelines for the analysis of germline CDH1 sequence variants

The variant curation guidelines published in 2015 by the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) provided the genetics community with a framework to assess variant pathogenicity; however, these rules are not gene specific. Germline pathogenic variants in the CDH1 gene cause hereditary diffuse gastric cancer and lobular breast cancer, a clinically challenging cancer predisposition syndrome that often requires a multidisciplinary team of experts to be properly managed. Given this challenge, the Clinical Genome Resource (ClinGen) Hereditary Cancer Domain prioritized the development of the CDH1 variant curation expert panel (VCEP) to develop and implement rules for CDH1 variant classifications. Here, we describe the CDH1 specifications of the ACMG/AMP guidelines, which were developed and validated after a systematic evaluation of variants obtained from a cohort of clinical laboratory data encompassing ∼827,000 CDH1 sequenced alleles. Comparing previously reported germline variants that were classified using the 2015 ACMG/AMP guidelines to the CDH1 VCEP recommendations resulted in reduced variants of uncertain significance and facilitated resolution of variants with conflicted assertions in ClinVar. The ClinGen CDH1 VCEP recommends the use of these CDH1-specific guidelines for the assessment and classification of variants identified in this clinically actionable gene.

Association of angiotensin-converting enzyme 2 gene polymorphism and enzymatic activity with essential hypertension in different gender: A case-control study

Angiotensin-converting enzyme 2 (ACE2) plays an important role in the development of essential hypertension (EH). The aim of this study was to investigate the relationship of ACE2 gene polymorphisms and enzymatic activity with EH in the northeastern Chinese Han population. 34 single-nucleotide polymorphism (SNP) loci of ACE2 were detected in 1024 EH patients and 956 normotensive (NT) controls by Sequenom Mass-ARRAY RS1000. Five SNPs (rs1514283, rs4646155, rs4646176, rs2285666, and rs879922) in ACE2 gene were determined to significantly associate with EH in female participants, while no SNP locus was linked to male group. Specifically, it was the first time to report that rs4646155 was significantly associated with EH in females. Furthermore, the correlation between ACE2 activity and clinical parameters were performed by Pearson correlation analysis in EH patients. We found that the ACE2 activity level was negatively correlated with body mass index (BMI), DBP, and pulse pressure, and significantly positively with ACE2 concentration, blood glucose and estrogen level in female EH patients. These results demonstrated that the genetic variants of ACE2 played vital roles in the development of EH. And the serum ACE2 activity can predict the development of cardiac dysfunction in EH patients.

Recommendations for interpreting the loss of function PVS1 ACMG/AMP variant criterion

The 2015 ACMG/AMP sequence variant interpretation guideline provided a framework for classifying variants based on several benign and pathogenic evidence criteria, including a pathogenic criterion (PVS1) for predicted loss of function variants. However, the guideline did not elaborate on specific considerations for the different types of loss of function variants, nor did it provide decision-making pathways assimilating information about variant type, its location, or any additional evidence for the likelihood of a true null effect. Furthermore, this guideline did not take into account the relative strengths for each evidence type and the final outcome of their combinations with respect to PVS1 strength. Finally, criteria specifying the genes for which PVS1 can be applied are still missing. Here, as part of the ClinGen Sequence Variant Interpretation (SVI) Workgroup's goal of refining ACMG/AMP criteria, we provide recommendations for applying the PVS1 criterion using detailed guidance addressing the above-mentioned gaps. Evaluation of the refined criterion by seven disease-specific groups using heterogeneous types of loss of function variants (n = 56) showed 89% agreement with the new recommendation, while discrepancies in six variants (11%) were appropriately due to disease-specific refinements. Our recommendations will facilitate consistent and accurate interpretation of predicted loss of function variants.

Identifying Genes Whose Mutant Transcripts Cause Dominant Disease Traits by Potential Gain-of-Function Alleles

Premature termination codon (PTC)-bearing transcripts are often degraded by nonsense-mediated decay (NMD) resulting in loss-of-function (LoF) alleles. However, not all PTCs result in LoF mutations, i.e., some such transcripts escape NMD and are translated to truncated peptide products that result in disease due to gain-of-function (GoF) effects. Since the location of the PTC is a major factor determining transcript fate, we hypothesized that depletion of protein-truncating variants (PTVs) within the gene region predicted to escape NMD in control databases could provide a rank for genic susceptibility for disease through GoF versus LoF. We developed an NMD escape intolerance score to rank genes based on the depletion of PTVs that would render them able to escape NMD using the Atherosclerosis Risk in Communities Study (ARIC) and the Exome Aggregation Consortium (ExAC) control databases, which was further used to screen the Baylor-Center for Mendelian Genomics disease database. This analysis revealed 1,996 genes significantly depleted for PTVs that are predicted to escape from NMD, i.e., PTVesc; further studies provided evidence that revealed a subset as candidate genes underlying Mendelian phenotypes. Importantly, these genes have characteristically low pLI scores, which can cause them to be overlooked as candidates for dominant diseases. Collectively, we demonstrate that this NMD escape intolerance score is an effective and efficient tool for gene discovery in Mendelian diseases due to production of truncated or altered proteins. More importantly, we provide a complementary analytical tool to aid identification of genes associated with dominant traits through a mechanism distinct from LoF.

A novel de novo CDH1 germline variant aids in the classification of carboxy-terminal E-cadherin alterations predicted to escape nonsense-mediated mRNA decay

Most truncating cadherin 1 (CDH1) pathogenic alterations confer an elevated lifetime risk of diffuse gastric cancer (DGC) and lobular breast cancer (LBC). However, transcripts containing carboxy-terminal premature stop codons have been demonstrated to escape the nonsense-mediated mRNA decay pathway, and gastric and breast cancer risks associated with these truncations should be carefully evaluated. A female patient underwent multigene panel testing because of a personal history of invasive LBC diagnosed at age 54, which identified the germline CDH1 nonsense alteration, c.2506G>T (p.Glu836*), in the last exon of the gene. Subsequent parental testing for the alteration was negative and additional short tandem repeat analysis confirmed the familial relationships and the de novo occurrence in the proband. Based on the de novo occurrence, clinical history, and rarity in general population databases, this alteration was classified as a likely pathogenic variant. This is the most carboxy-terminal pathogenic alteration reported to date. Additionally, this alteration contributed to the classification of six other upstream CDH1 carboxy-terminal truncating variants as pathogenic or likely pathogenic. Identifying the most distal pathogenic alteration provides evidence to classify other carboxy-terminal truncating variants as either pathogenic or benign, a fundamental step to offering presymptomatic screening and prophylactic procedures to the appropriate patients.

Splicing heterogeneity: separating signal from noise

Single-cell analyses have revealed a tremendous variety among cells in the abundance and chemical composition of RNA. Much of this heterogeneity is due to alternative splicing by the spliceosome. Little is known about how many of the resulting isoforms are biologically functional or just provide noise with little to no impact. The dynamic nature of the spliceosome provides numerous opportunities for regulation but is also the source of stochastic fluctuations. We discuss possible origins of splicing stochasticity, the experimental approaches for studying heterogeneity in isoforms, and the potential biological significance of noisy splicing in development and disease.

Base-specific mutational intolerance near splice sites clarifies the role of nonessential splice nucleotides

Variation in RNA splicing (i.e., alternative splicing) plays an important role in many diseases. Variants near 5' and 3' splice sites often affect splicing, but the effects of these variants on splicing and disease have not been fully characterized beyond the two "essential" splice nucleotides flanking each exon. Here we provide quantitative measurements of tolerance to mutational disruptions by position and reference allele-alternative allele combinations. We show that certain reference alleles are particularly sensitive to mutations, regardless of the alternative alleles into which they are mutated. Using public RNA-seq data, we demonstrate that individuals carrying such variants have significantly lower levels of the correctly spliced transcript, compared to individuals without them, and confirm that these specific substitutions are highly enriched for known Mendelian mutations. Our results propose a more refined definition of the "splice region" and offer a new way to prioritize and provide functional interpretation of variants identified in diagnostic sequencing and association studies.

Cancer driver mutation prediction through Bayesian integration of multi-omic data

Identification of cancer driver mutations is critical for advancing cancer research and personalized medicine. Due to inter-tumor genetic heterogeneity, many driver mutations occur at low frequencies, which make it challenging to distinguish them from passenger mutations. Here, we show that a novel Bayesian hierarchical modeling approach, named rDriver can achieve enhanced prediction accuracy by identifying mutations that not only have high functional impact scores but also are associated with systemic variation in gene expression levels. In examining 3,080 tumor samples from 8 cancer types in The Cancer Genome Atlas, rDriver predicted 1,389 driver mutations. Compared with existing tools, rDriver identified more low frequency mutations associated with lineage specific functional properties, timing of occurrence and patient survival. Evaluation of rDriver predictions using engineered cell-line models resulted in a positive predictive value of 0.94 in PIK3CA genes. Our study highlights the importance of integrating multi-omic data in predicting cancer driver mutations and provides a statistically rigorous solution for cancer target discovery and development.

Role of duplicate genes in determining the tissue-selectivity of hereditary diseases

A longstanding puzzle in human genetics is what limits the clinical manifestation of hundreds of hereditary diseases to certain tissues, while their causal genes are expressed throughout the human body. A general conception is that tissue-selective disease phenotypes emerge when masking factors operate in unaffected tissues, but are specifically absent or insufficient in disease-manifesting tissues. Although this conception has critical impact on the understanding of disease manifestation, it was never challenged in a systematic manner across a variety of hereditary diseases and affected tissues. Here, we address this gap in our understanding via rigorous analysis of the susceptibility of over 30 tissues to 112 tissue-selective hereditary diseases. We focused on the roles of paralogs of causal genes, which are presumably capable of compensating for their aberration. We show for the first time at large-scale via quantitative analysis of omics datasets that, preferentially in the disease-manifesting tissues, paralogs are under-expressed relative to causal genes in more than half of the diseases. This was observed for several susceptible tissues and for causal genes with varying number of paralogs, suggesting that imbalanced expression of paralogs increases tissue susceptibility. While for many diseases this imbalance stemmed from up-regulation of the causal gene in the disease-manifesting tissue relative to other tissues, it was often combined with down-regulation of its paralog. Notably in roughly 20% of the cases, this imbalance stemmed only from significant down-regulation of the paralog. Thus, dosage relationships between paralogs appear as important, yet currently under-appreciated, modifiers of disease manifestation.

Medical relevance of protein-truncating variants across 337,205 individuals in the UK Biobank study

Protein-truncating variants can have profound effects on gene function and are critical for clinical genome interpretation and generating therapeutic hypotheses, but their relevance to medical phenotypes has not been systematically assessed. Here, we characterize the effect of 18,228 protein-truncating variants across 135 phenotypes from the UK Biobank and find 27 associations between medical phenotypes and protein-truncating variants in genes outside the major histocompatibility complex. We perform phenome-wide analyses and directly measure the effect in homozygous carriers, commonly referred to as “human knockouts,” across medical phenotypes for genes implicated as being protective against disease or associated with at least one phenotype in our study. We find several genes with strong pleiotropic or non-additive effects. Our results illustrate the importance of protein-truncating variants in a variety of diseases.

Protein-truncating variants (PTVs) are predicted to significantly affect a gene’s function and, thus, human traits. Here, DeBoever et al. systematically analyze PTVs in more than 300,000 individuals across 135 phenotypes and identify 27 associations between PTVs and medical conditions.

Mutations in the pancreatic secretory enzymes CPA1 and CPB1 are associated with pancreatic cancer

To evaluate whether germline variants in genes encoding pancreatic secretory enzymes contribute to pancreatic cancer susceptibility, we sequenced the coding regions of CPB1 and other genes encoding pancreatic secretory enzymes and known pancreatitis susceptibility genes (PRSS1, CPA1, CTRC, and SPINK1) in a hospital series of pancreatic cancer cases and controls. Variants in CPB1, CPA1 (encoding carboxypeptidase B1 and A1), and CTRC were evaluated in a second set of cases with familial pancreatic cancer and controls. More deleterious CPB1 variants, defined as having impaired protein secretion and induction of endoplasmic reticulum (ER) stress in transfected HEK 293T cells, were found in the hospital series of pancreatic cancer cases (5/986, 0.5%) than in controls (0/1,045, P = 0.027). Among familial pancreatic cancer cases, ER stress-inducing CPB1 variants were found in 4 of 593 (0.67%) vs. 0 of 967 additional controls (P = 0.020), with a combined prevalence in pancreatic cancer cases of 9/1,579 vs. 0/2,012 controls (P < 0.01). More ER stress-inducing CPA1 variants were also found in the combined set of hospital and familial cases with pancreatic cancer than in controls [7/1,546 vs. 1/2,012; P = 0.025; odds ratio, 9.36 (95% CI, 1.15-76.02)]. Overall, 16 (1%) of 1,579 pancreatic cancer cases had an ER stress-inducing CPA1 or CPB1 variant, compared with 1 of 2,068 controls (P < 0.00001). No other candidate genes had statistically significant differences in variant prevalence between cases and controls. Our study indicates ER stress-inducing variants in CPB1 and CPA1 are associated with pancreatic cancer susceptibility and implicate ER stress in pancreatic acinar cells in pancreatic cancer development.

Genome-wide CRISPR-Cas9 Screen Identifies Leukemia-Specific Dependence on a Pre-mRNA Metabolic Pathway Regulated by DCPS

To identify novel targets for acute myeloid leukemia (AML) therapy, we performed genome-wide CRISPR-Cas9 screening using AML cell lines, followed by a second screen in vivo. Here, we show that the mRNA decapping enzyme scavenger (DCPS) gene is essential for AML cell survival. The DCPS enzyme interacted with components of pre-mRNA metabolic pathways, including spliceosomes, as revealed by mass spectrometry. RG3039, a DCPS inhibitor originally developed to treat spinal muscular atrophy, exhibited anti-leukemic activity via inducing pre-mRNA mis-splicing. Humans harboring germline biallelic DCPS loss-of-function mutations do not exhibit aberrant hematologic phenotypes, indicating that DCPS is dispensable for human hematopoiesis. Our findings shed light on a pre-mRNA metabolic pathway and identify DCPS as a target for AML therapy.

Role of titin in cardiomyopathy: from DNA variants to patient stratification

Dilated cardiomyopathy (DCM) affects approximately 1 in 250 individuals and is the leading indication for heart transplantation. DCM is often familial, and the most common genetic predisposition is a truncating variation in the giant sarcomeric protein, titin, which occurs in up to 15% of ambulant patients with DCM and 25% of end-stage or familial cases. In this article, we review the evidence for the role of titin truncation in the pathogenesis of DCM and our understanding of the molecular mechanisms and pathophysiological consequences of variation in the gene encoding titin (TTN). Such variation is common in the general population (up to 1% of individuals), and we consider key features that discriminate variants with disease-causing potential from those that are benign. We summarize strategies for clinical interpretation of genetic variants for use in the diagnosis of patients and the evaluation of their relatives. Finally, we consider the contemporary and potential future role for genetic stratification in cardiomyopathy and in the general population, evaluating titin variation as a predictor of outcome and treatment response for precision medicine.

The contribution of alternative splicing to genetic risk for psychiatric disorders

A genetic contribution to psychiatric disorders has clearly been established and genome-wide association studies now provide the location of risk genes and genetic variants associated with risk. However, the mechanism by which these genes and variants contribute to psychiatric disorders is mostly undetermined. This is in part because non-synonymous protein coding changes cannot explain the majority of variants associated with complex genetic traits. Based on this, it is predicted that these variants are causing gene expression changes, including changes to alternative splicing. Genetic changes influencing alternative splicing have been identified as risk factors in Mendelian disorders; however, currently there is a paucity of research on the role of alternative splicing in complex traits. This stems partly from the difficulty of predicting the role of genetic variation in splicing. Alterations to canonical splice site sequences, nucleotides adjacent to splice junctions, and exonic and intronic splicing regulatory sequences can influence splice site choice. Recent studies have identified global changes in alternatively spliced transcripts in brain tissues, some of which correlate with altered levels of splicing trans factors. Disease-associated variants have also been found to affect cis-acting splicing regulatory sequences and alter the ratio of alternatively spliced transcripts. These findings are reviewed here, as well as the current datasets and resources available to study alternative splicing in psychiatric disorders. Identifying and understanding risk variants that cause alternative splicing is critical to understanding the mechanisms of risk as well as to pave the way for new therapeutic options.

Enhancing GTEx by bridging the gaps between genotype, gene expression, and disease

Genetic variants have been associated with myriad molecular phenotypes that provide new insight into the range of mechanisms underlying genetic traits and diseases. Identifying any particular genetic variant's cascade of effects, from molecule to individual, requires assaying multiple layers of molecular complexity. We introduce the Enhancing GTEx (eGTEx) project that extends the GTEx project to combine gene expression with additional intermediate molecular measurements on the same tissues to provide a resource for studying how genetic differences cascade through molecular phenotypes to impact human health.

RNA expression in human retina

Recent Genome-wide Association Studies (GWASs) for eye diseases/traits have delivered a number of novel findings across a diverse range of diseases, including age-related macular degeneration (AMD), glaucoma and refractive error. However, despite this astonishing rate of success, the major challenge still remains to not only confirm that the genes implicated in these studies are truly the genes conferring protection from or risk of disease but also to define the functional roles these genes play in disease. Ongoing evidence is accumulating that the single nucleotide polymorphisms (SNPs) used in GWAS and fine mapping studies have causal effects through their influence on gene expression rather than affecting protein function. The biological interpretation of SNP regulatory effects for a tissue requires knowledge of the transcriptome for that tissue. We summarize the reasons to characterize the complete retinal transcriptome as well as the evidence to include an assessment of differences in regional retinal expression.

A Comparison of mRNA Sequencing with Random Primed and 3'-Directed Libraries

Creating a cDNA library for deep mRNA sequencing (mRNAseq) is generally done by random priming, creating multiple sequencing fragments along each transcript. A 3'-end-focused library approach cannot detect differential splicing, but has potentially higher throughput at a lower cost, along with the ability to improve quantification by using transcript molecule counting with unique molecular identifiers (UMI) that correct PCR bias. Here, we compare an implementation of such a 3'-digital gene expression (3'-DGE) approach with "conventional" random primed mRNAseq. Given our particular datasets on cultured human cardiomyocyte cell lines, we find that, while conventional mRNAseq detects ~15% more genes and needs ~500,000 fewer reads per sample for equivalent statistical power, the resulting differentially expressed genes, biological conclusions, and gene signatures are highly concordant between two techniques. We also find good quantitative agreement at the level of individual genes between two techniques for both read counts and fold changes between given conditions. We conclude that, for high-throughput applications, the potential cost savings associated with 3'-DGE approach are likely a reasonable tradeoff for modest reduction in sensitivity and inability to observe alternative splicing, and should enable many larger scale studies focusing on not only differential expression analysis, but also quantitative transcriptome profiling.

Integrating RNA sequencing into neuro-oncology practice

Malignant tumors of the central nervous system (CNS) cause substantial morbidity and mortality, yet efforts to optimize chemo- and radiotherapy have largely failed to improve dismal prognoses. Over the past decade, RNA sequencing (RNA-seq) has emerged as a powerful tool to comprehensively characterize the transcriptome of CNS tumor cells in one high-throughput step, leading to improved understanding of CNS tumor biology and suggesting new routes for targeted therapies. RNA-seq has been instrumental in improving the diagnostic classification of brain tumors, characterizing oncogenic fusion genes, and shedding light on intratumor heterogeneity. Currently, RNA-seq is beginning to be incorporated into regular neuro-oncology practice in the form of precision neuro-oncology programs, which use information from tumor sequencing to guide implementation of personalized targeted therapies. These programs show great promise in improving patient outcomes for tumors where single agent trials have been ineffective. As RNA-seq is a relatively new technique, many further applications yielding new advances in CNS tumor research and management are expected in the coming years.

Quantifying the regulatory effect size of cis-acting genetic variation using allelic fold change

Mapping cis-acting expression quantitative trait loci (cis-eQTL) has become a popular approach for characterizing proximal genetic regulatory variants. In this paper, we describe and characterize log allelic fold change (aFC), the magnitude of expression change associated with a given genetic variant, as a biologically interpretable unit for quantifying the effect size of cis-eQTLs and a mathematically convenient approach for systematic modeling of cis-regulation. This measure is mathematically independent from expression level and allele frequency, additive, applicable to multiallelic variants, and generalizable to multiple independent variants. We provide efficient tools and guidelines for estimating aFC from both eQTL and allelic expression data sets and apply it to Genotype Tissue Expression (GTEx) data. We show that aFC estimates independently derived from eQTL and allelic expression data are highly consistent, and identify technical and biological correlates of eQTL effect size. We generalize aFC to analyze genes with two eQTLs in GTEx and show that in nearly all cases the two eQTLs act independently in regulating gene expression. In summary, aFC is a solid measure of cis-regulatory effect size that allows quantitative interpretation of cellular regulatory events from population data, and it is a valuable approach for investigating novel aspects of eQTL data sets.

Whole exome sequencing in 342 congenital cardiac left sided lesion cases reveals extensive genetic heterogeneity and complex inheritance patterns

BACKGROUND

Left-sided lesions (LSLs) account for an important fraction of severe congenital cardiovascular malformations (CVMs). The genetic contributions to LSLs are complex, and the mutations that cause these malformations span several diverse biological signaling pathways: TGFB, NOTCH, SHH, and more. Here, we use whole exome sequence data generated in 342 LSL cases to identify likely damaging variants in putative candidate CVM genes.

METHODS

Using a series of bioinformatics filters, we focused on genes harboring population-rare, putative loss-of-function (LOF), and predicted damaging variants in 1760 CVM candidate genes constructed a priori from the literature and model organism databases. Gene variants that were not observed in a comparably sequenced control dataset of 5492 samples without severe CVM were then subjected to targeted validation in cases and parents. Whole exome sequencing data from 4593 individuals referred for clinical sequencing were used to bolster evidence for the role of candidate genes in CVMs and LSLs.

RESULTS

Our analyses revealed 28 candidate variants in 27 genes, including 17 genes not previously associated with a human CVM disorder, and revealed diverse patterns of inheritance among LOF carriers, including 9 confirmed de novo variants in both novel and newly described human CVM candidate genes (ACVR1, JARID2, NR2F2, PLRG1, SMURF1) as well as established syndromic CVM genes (KMT2D, NF1, TBX20, ZEB2). We also identified two genes (DNAH5, OFD1) with evidence of recessive and hemizygous inheritance patterns, respectively. Within our clinical cohort, we also observed heterozygous LOF variants in JARID2 and SMAD1 in individuals with cardiac phenotypes, and collectively, carriers of LOF variants in our candidate genes had a four times higher odds of having CVM (odds ratio = 4.0, 95% confidence interval 2.5-6.5).

CONCLUSIONS

Our analytical strategy highlights the utility of bioinformatic resources, including human disease records and model organism phenotyping, in novel gene discovery for rare human disease. The results underscore the extensive genetic heterogeneity underlying non-syndromic LSLs, and posit potential novel candidate genes and complex modes of inheritance in this important group of birth defects.

Autoregulation of RBM10 and cross-regulation of RBM10/RBM5 via alternative splicing-coupled nonsense-mediated decay

Mutations in the spliceosomal RNA binding protein RBM10 cause TARP syndrome and are frequently observed in lung adenocarcinoma (LUAD). We have previously shown that RBM10 enhances exon skipping of its target genes, including its paralog RBM5. Here, we report that RBM10 negatively regulates its own mRNA and protein expression and that of RBM5 by promoting alternative splicing-coupled nonsense-mediated mRNA decay (AS-NMD). Through computational analysis and experimental validation, we identified RBM10-promoted skipping of exon 6 or 12 in RBM10 and exon 6 or 16 in RBM5 as the underlying AS-NMD events. Importantly, we showed that LUAD-associated mutations affecting splice sites of RBM10 exons 6 or 12 abolished exon inclusion and correlated with reduced expression of RBM10 RNA. Together, our investigations have revealed novel molecular mechanisms underlying RBM10 autoregulation and cross-regulation of RBM5, thereby providing insights concerning the functions of RBM10 under various physiological and pathological conditions. Our combined computational and experimental approach should be useful for elucidating the role of AS-NMD in auto- and cross-regulation by other splicing regulators.

Genetic effects on gene expression across human tissues

Characterization of the molecular function of the human genome and its variation across individuals is essential for identifying the cellular mechanisms that underlie human genetic traits and diseases. The Genotype-Tissue Expression (GTEx) project aims to characterize variation in gene expression levels across individuals and diverse tissues of the human body, many of which are not easily accessible. Here we describe genetic effects on gene expression levels across 44 human tissues. We find that local genetic variation affects gene expression levels for the majority of genes, and we further identify inter-chromosomal genetic effects for 93 genes and 112 loci. On the basis of the identified genetic effects, we characterize patterns of tissue specificity, compare local and distal effects, and evaluate the functional properties of the genetic effects. We also demonstrate that multi-tissue, multi-individual data can be used to identify genes and pathways affected by human disease-associated variation, enabling a mechanistic interpretation of gene regulation and the genetic basis of disease.

Landscape of X chromosome inactivation across human tissues

X chromosome inactivation (XCI) silences transcription from one of the two X chromosomes in female mammalian cells to balance expression dosage between XX females and XY males. XCI is, however, incomplete in humans: up to one-third of X-chromosomal genes are expressed from both the active and inactive X chromosomes (Xa and Xi, respectively) in female cells, with the degree of 'escape' from inactivation varying between genes and individuals. The extent to which XCI is shared between cells and tissues remains poorly characterized, as does the degree to which incomplete XCI manifests as detectable sex differences in gene expression and phenotypic traits. Here we describe a systematic survey of XCI, integrating over 5,500 transcriptomes from 449 individuals spanning 29 tissues from GTEx (v6p release) and 940 single-cell transcriptomes, combined with genomic sequence data. We show that XCI at 683 X-chromosomal genes is generally uniform across human tissues, but identify examples of heterogeneity between tissues, individuals and cells. We show that incomplete XCI affects at least 23% of X-chromosomal genes, identify seven genes that escape XCI with support from multiple lines of evidence and demonstrate that escape from XCI results in sex biases in gene expression, establishing incomplete XCI as a mechanism that is likely to introduce phenotypic diversity. Overall, this updated catalogue of XCI across human tissues helps to increase our understanding of the extent and impact of the incompleteness in the maintenance of XCI.

Insights into beta cell regeneration for diabetes via integration of molecular landscapes in human insulinomas

Although diabetes results in part from a deficiency of normal pancreatic beta cells, inducing human beta cells to regenerate is difficult. Reasoning that insulinomas hold the “genomic recipe” for beta cell expansion, we surveyed 38 human insulinomas to obtain insights into therapeutic pathways for beta cell regeneration. An integrative analysis of whole-exome and RNA-sequencing data was employed to extensively characterize the genomic and molecular landscape of insulinomas relative to normal beta cells. Here, we show at the pathway level that the majority of the insulinomas display mutations, copy number variants and/or dysregulation of epigenetic modifying genes, most prominently in the polycomb and trithorax families. Importantly, these processes are coupled to co-expression network modules associated with cell proliferation, revealing candidates for inducing beta cell regeneration. Validation of key computational predictions supports the concept that understanding the molecular complexity of insulinoma may be a valuable approach to diabetes drug discovery.

Diabetes results in part from a deficiency of functional pancreatic beta cells. Here, the authors study the genomic and epigenetic landscapes of human insulinomas to gain insight into possible pathways for therapeutic beta cell regeneration, highlighting epigenetic genes and pathways.