A validated heart-specific model for splice-disrupting variants in childhood heart disease.

BACKGROUND: Congenital heart disease (CHD) is the most common congenital anomaly. Almost 90% of isolated cases have an unexplained genetic etiology after clinical testing. Non-canonical splice variants that disrupt mRNA splicing through the loss or creation of exon boundaries are not routinely captured and/or evaluated by standard clinical genetic tests. Recent computational algorithms such as SpliceAI have shown an ability to predict such variants, but are not specific to cardiac-expressed genes and transcriptional isoforms. METHODS: We used genome sequencing (GS) (n = 1101 CHD probands) and myocardial RNA-Sequencing (RNA-Seq) (n = 154 CHD and n = 43 cardiomyopathy probands) to identify and validate splice disrupting variants, and to develop a heart-specific model for canonical and non-canonical splice variants that can be applied to patients with CHD and cardiomyopathy. Two thousand five hundred seventy GS samples from the Medical Genome Reference Bank were analyzed as healthy controls. RESULTS: Of 8583 rare DNA splice-disrupting variants initially identified using SpliceAI, 100 were associated with altered splice junctions in the corresponding patient myocardium affecting 95 genes. Using strength of myocardial gene expression and genome-wide DNA variant features that were confirmed to affect splicing in myocardial RNA, we trained a machine learning model for predicting cardiac-specific splice-disrupting variants (AUC 0.86 on internal validation). In a validation set of 48 CHD probands, the cardiac-specific model outperformed a SpliceAI model alone (AUC 0.94 vs 0.67 respectively). Application of this model to an additional 947 CHD probands with only GS data identified 1% patients with canonical and 11% patients with non-canonical splice-disrupting variants in CHD genes. Forty-nine percent of predicted splice-disrupting variants were intronic and > 10 bp from existing splice junctions. The burden of high-confidence splice-disrupting variants in CHD genes was 1.28-fold higher in CHD cases compared with healthy controls. CONCLUSIONS: A new cardiac-specific in silico model was developed using complementary GS and RNA-Seq data that improved genetic yield by identifying a significant burden of non-canonical splice variants associated with CHD that would not be detectable through panel or exome sequencing.

Genome-wide enhancer-associated tandem repeats are expanded in cardiomyopathy.

BACKGROUND: Cardiomyopathy is a clinically and genetically heterogeneous heart condition that can lead to heart failure and sudden cardiac death in childhood. While it has a strong genetic basis, the genetic aetiology for over 50% of cardiomyopathy cases remains unknown. METHODS: In this study, we analyse the characteristics of tandem repeats from genome sequence data of unrelated individuals diagnosed with cardiomyopathy from Canada and the United Kingdom (n = 1216) and compare them to those found in the general population. We perform burden analysis to identify genomic and epigenomic features that are impacted by rare tandem repeat expansions (TREs), and enrichment analysis to identify functional pathways that are involved in the TRE-associated genes in cardiomyopathy. We use Oxford Nanopore targeted long-read sequencing to validate repeat size and methylation status of one of the most recurrent TREs. We also compare the TRE-associated genes to those that are dysregulated in the heart tissues of individuals with cardiomyopathy. FINDINGS: We demonstrate that tandem repeats that are rarely expanded in the general population are predominantly expanded in cardiomyopathy. We find that rare TREs are disproportionately present in constrained genes near transcriptional start sites, have high GC content, and frequently overlap active enhancer H3K27ac marks, where expansion-related DNA methylation may reduce gene expression. We demonstrate the gene silencing effect of expanded CGG tandem repeats in DIP2B through promoter hypermethylation. We show that the enhancer-associated loci are found in genes that are highly expressed in human cardiomyocytes and are differentially expressed in the left ventricle of the heart in individuals with cardiomyopathy. INTERPRETATION: Our findings highlight the underrecognized contribution of rare tandem repeat expansions to the risk of cardiomyopathy and suggest that rare TREs contribute to ∼4% of cardiomyopathy risk. FUNDING: Government of Ontario (RKCY), The Canadian Institutes of Health Research PJT 175329 (RKCY), The Azrieli Foundation (RKCY), SickKids Catalyst Scholar in Genetics (RKCY), The University of Toronto McLaughlin Centre (RKCY, SM), Ted Rogers Centre for Heart Research (SM), Data Sciences Institute at the University of Toronto (SM), The Canadian Institutes of Health Research PJT 175034 (SM), The Canadian Institutes of Health Research ENP 161429 under the frame of ERA PerMed (SM, RL), Heart and Stroke Foundation of Ontario & Robert M Freedom Chair in Cardiovascular Science (SM), Bitove Family Professorship of Adult Congenital Heart Disease (EO), Canada Foundation for Innovation (SWS, JR), Canada Research Chair (PS), Genome Canada (PS, JR), The Canadian Institutes of Health Research (PS).

Age and Sex Differences in the Genetics of Cardiomyopathy.

Cardiomyopathy has variable penetrance. We analyzed age and sex-related genetic differences in 1,397 cardiomyopathy patients (Ontario, UK) with whole genome sequencing. Pediatric cases (n = 471) harbored more deleterious protein-coding variants in Tier 1 cardiomyopathy genes compared to adults (n = 926) (34.6% vs 25.9% respectively, p = 0.0015), with variant enrichment in constrained coding regions. Pediatric patients had a higher burden of sarcomere and lower burden of channelopathy gene variants compared to adults. Specifically, pediatric patients had more MYH7 and MYL3 variants in hypertrophic cardiomyopathy, and fewer TTN truncating variants in dilated cardiomyopathy. MYH7 variants clustered in the myosin head and neck domains in children. OBSCN was a top mutated gene in adults, enriched for protein-truncating variants. In dilated cardiomyopathy, female patients had a higher burden of z-disc gene variants compared to males. Genetic differences may explain age and sex-related variability in cardiomyopathy penetrance. Genotype-guided predictions of age of onset can inform pre-test genetic counseling. Pediatric cardiomyopathy patients were more likely to be genotype-positive than adults with a higher burden of variants in MYH7, MYL3, TNNT2, VCL. Adults had a higher burden of OBSCN and TTN variants. Females with dilated cardiomyopathy (DCM) had a higher burden of z-disc gene variants compared to males.

Whole genome sequencing delineates regulatory, copy number, and cryptic splice variants in early onset cardiomyopathy.

Cardiomyopathy (CMP) is a heritable disorder. Over 50% of cases are gene-elusive on clinical gene panel testing. The contribution of variants in non-coding DNA elements that result in cryptic splicing and regulate gene expression has not been explored. We analyzed whole-genome sequencing (WGS) data in a discovery cohort of 209 pediatric CMP patients and 1953 independent replication genomes and exomes. We searched for protein-coding variants, and non-coding variants predicted to affect the function or expression of genes. Thirty-nine percent of cases harbored pathogenic coding variants in known CMP genes, and 5% harbored high-risk loss-of-function (LoF) variants in additional candidate CMP genes. Fifteen percent harbored high-risk regulatory variants in promoters and enhancers of CMP genes (odds ratio 2.25, p = 6.70 × 10-7 versus controls). Genes involved in α-dystroglycan glycosylation (FKTN, DTNA) and desmosomal signaling (DSC2, DSG2) were most highly enriched for regulatory variants (odds ratio 6.7-58.1). Functional effects were confirmed in patient myocardium and reporter assays in human cardiomyocytes, and in zebrafish CRISPR knockouts. We provide strong evidence for the genomic contribution of functionally active variants in new genes and in regulatory elements of known CMP genes to early onset CMP.

The telomere length landscape of prostate cancer.

Replicative immortality is a hallmark of cancer, and can be achieved through telomere lengthening and maintenance. Although the role of telomere length in cancer has been well studied, its association to genomic features is less well known. Here, we report the telomere lengths of 392 localized prostate cancer tumours and characterize their relationship to genomic, transcriptomic and proteomic features. Shorter tumour telomere lengths are associated with elevated genomic instability, including single-nucleotide variants, indels and structural variants. Genes involved in cell proliferation and signaling are correlated with tumour telomere length at all levels of the central dogma. Telomere length is also associated with multiple clinical features of a tumour. Longer telomere lengths in non-tumour samples are associated with a lower rate of biochemical relapse. In summary, we describe the multi-level integration of telomere length, genomics, transcriptomics and proteomics in localized prostate cancer.

Rare variants in KDR, encoding VEGF Receptor 2, are associated with tetralogy of Fallot.

PURPOSE: Rare genetic variants in KDR, encoding the vascular endothelial growth factor receptor 2 (VEGFR2), have been reported in patients with tetralogy of Fallot (TOF). However, their role in disease causality and pathogenesis remains unclear. METHODS: We conducted exome sequencing in a familial case of TOF and large-scale genetic studies, including burden testing, in >1,500 patients with TOF. We studied gene-targeted mice and conducted cell-based assays to explore the role of KDR genetic variation in the etiology of TOF. RESULTS: Exome sequencing in a family with two siblings affected by TOF revealed biallelic missense variants in KDR. Studies in knock-in mice and in HEK 293T cells identified embryonic lethality for one variant when occurring in the homozygous state, and a significantly reduced VEGFR2 phosphorylation for both variants. Rare variant burden analysis conducted in a set of 1,569 patients of European descent with TOF identified a 46-fold enrichment of protein-truncating variants (PTVs) in TOF cases compared to controls (P = 7 × 10-11). CONCLUSION: Rare KDR variants, in particular PTVs, strongly associate with TOF, likely in the setting of different inheritance patterns. Supported by genetic and in vivo and in vitro functional analysis, we propose loss-of-function of VEGFR2 as one of the mechanisms involved in the pathogenesis of TOF.

Molecular landmarks of tumor hypoxia across cancer types.

Many primary-tumor subregions have low levels of molecular oxygen, termed hypoxia. Hypoxic tumors are at elevated risk for local failure and distant metastasis, but the molecular hallmarks of tumor hypoxia remain poorly defined. To fill this gap, we quantified hypoxia in 8,006 tumors across 19 tumor types. In ten tumor types, hypoxia was associated with elevated genomic instability. In all 19 tumor types, hypoxic tumors exhibited characteristic driver-mutation signatures. We observed widespread hypoxia-associated dysregulation of microRNAs (miRNAs) across cancers and functionally validated miR-133a-3p as a hypoxia-modulated miRNA. In localized prostate cancer, hypoxia was associated with elevated rates of chromothripsis, allelic loss of PTEN and shorter telomeres. These associations are particularly enriched in polyclonal tumors, representing a constellation of features resembling tumor nimbosus, an aggressive cellular phenotype. Overall, this work establishes that tumor hypoxia may drive aggressive molecular features across cancers and shape the clinical trajectory of individual tumors.

Melorheostosis: Exome sequencing of an associated dermatosis implicates postzygotic mosaicism of mutated KRAS.

Melorheostosis (MEL) is the rare sporadic dysostosis characterized by monostotic or polyostotic osteosclerosis and hyperostosis often distributed in a sclerotomal pattern. The prevailing hypothesis for MEL invokes postzygotic mosaicism. Sometimes scleroderma-like skin changes, considered a representation of the pathogenetic process of MEL, overlie the bony changes, and sometimes MEL becomes malignant. Osteopoikilosis (OPK) is the autosomal dominant skeletal dysplasia that features symmetrically distributed punctate osteosclerosis due to heterozygous loss-of-function mutation within LEMD3. Rarely, radiographic findings of MEL occur in OPK. However, germline mutation of LEMD3 does not explain sporadic MEL. To explore if mosaicism underlies MEL, we studied a boy with polyostotic MEL and characteristic overlying scleroderma-like skin, a few bony lesions consistent with OPK, and a large epidermal nevus known to usually harbor a HRAS, FGFR3, or PIK3CA gene mutation. Exome sequencing was performed to ~100× average read depth for his two dermatoses, two areas of normal skin, and peripheral blood leukocytes. As expected for non-malignant tissues, the patient's mutation burden in his normal skin and leukocytes was low. He, his mother, and his maternal grandfather carried a heterozygous, germline, in-frame, 24-base-pair deletion in LEMD3. Radiographs of the patient and his mother revealed bony foci consistent with OPK, but she showed no MEL. For the patient, somatic variant analysis, using four algorithms to compare all 20 possible pairwise combinations of his five DNA samples, identified only one high-confidence mutation, heterozygous KRAS Q61H (NM_033360.3:c.183A>C, NP_203524.1:p.Gln61His), in both his dermatoses but absent in his normal skin and blood. Thus, sparing our patient biopsy of his MEL bone, we identified a heterozygous somatic KRAS mutation in his scleroderma-like dermatosis considered a surrogate for MEL. This implicates postzygotic mosaicism of mutated KRAS, perhaps facilitated by germline LEMD3 haploinsufficiency, causing his MEL.

Molecular Features of Subtype-Specific Progression from Ductal Carcinoma In Situ to Invasive Breast Cancer.

Breast cancer consists of at least five main molecular "intrinsic" subtypes that are reflected in both pre-invasive and invasive disease. Although previous studies have suggested that many of the molecular features of invasive breast cancer are established early, it is unclear what mechanisms drive progression and whether the mechanisms of progression are dependent or independent of subtype. We have generated mRNA, miRNA, and DNA copy-number profiles from a total of 59 in situ lesions and 85 invasive tumors in order to comprehensively identify those genes, signaling pathways, processes, and cell types that are involved in breast cancer progression. Our work provides evidence that there are molecular features associated with disease progression that are unique to the intrinsic subtypes. We additionally establish subtype-specific signatures that are able to identify a small proportion of pre-invasive tumors with expression profiles that resemble invasive carcinoma, indicating a higher likelihood of future disease progression.

GenVisR: Genomic Visualizations in R.

UNLABELLED: Visualizing and summarizing data from genomic studies continues to be a challenge. Here, we introduce the GenVisR package to addresses this challenge by providing highly customizable, publication-quality graphics focused on cohort level genome analyses. GenVisR provides a rapid and easy-to-use suite of genomic visualization tools, while maintaining a high degree of flexibility by leveraging the abilities of ggplot2 and Bioconductor. AVAILABILITY AND IMPLEMENTATION: GenVisR is an R package available via Bioconductor (https://bioconductor.org/packages/GenVisR) under GPLv3. Support is available via GitHub (https://github.com/griffithlab/GenVisR/issues) and the Bioconductor support website. CONTACTS: obigriffith@wustl.edu or mgriffit@wustl.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Comprehensive genomic analysis reveals FLT3 activation and a therapeutic strategy for a patient with relapsed adult B-lymphoblastic leukemia.

The genomic events responsible for the pathogenesis of relapsed adult B-lymphoblastic leukemia (B-ALL) are not yet clear. We performed integrative analysis of whole-genome, whole-exome, custom capture, whole-transcriptome (RNA-seq), and locus-specific genomic assays across nine time points from a patient with primary de novo B-ALL. Comprehensive genome and transcriptome characterization revealed a dramatic tumor evolution during progression, yielding a tumor with complex clonal architecture at second relapse. We observed and validated point mutations in EP300 and NF1, a highly expressed EP300-ZNF384 gene fusion, a microdeletion in IKZF1, a focal deletion affecting SETD2, and large deletions affecting RB1, PAX5, NF1, and ETV6. Although the genome analysis revealed events of potential biological relevance, no clinically actionable treatment options were evident at the time of the second relapse. However, transcriptome analysis identified aberrant overexpression of the targetable protein kinase encoded by the FLT3 gene. Although the patient had refractory disease after salvage therapy for the second relapse, treatment with the FLT3 inhibitor sunitinib rapidly induced a near complete molecular response, permitting the patient to proceed to a matched-unrelated donor stem cell transplantation. The patient remains in complete remission more than 4 years later. Analysis of this patient's relapse genome revealed an unexpected, actionable therapeutic target that led to a specific therapy associated with a rapid clinical response. For some patients with relapsed or refractory cancers, this approach may indicate a novel therapeutic intervention that could alter outcome.

ORegAnno 3.0: a community-driven resource for curated regulatory annotation.

The Open Regulatory Annotation database (ORegAnno) is a resource for curated regulatory annotation. It contains information about regulatory regions, transcription factor binding sites, RNA binding sites, regulatory variants, haplotypes, and other regulatory elements. ORegAnno differentiates itself from other regulatory resources by facilitating crowd-sourced interpretation and annotation of regulatory observations from the literature and highly curated resources. It contains a comprehensive annotation scheme that aims to describe both the elements and outcomes of regulatory events. Moreover, ORegAnno assembles these disparate data sources and annotations into a single, high quality catalogue of curated regulatory information. The current release is an update of the database previously featured in the NAR Database Issue, and now contains 1 948 307 records, across 18 species, with a combined coverage of 334 215 080 bp. Complete records, annotation, and other associated data are available for browsing and download at http://www.oreganno.org/.

The prognostic ease and difficulty of invasive breast carcinoma.

Breast carcinoma (BC) has been extensively profiled by high-throughput technologies for over a decade, and broadly speaking, these studies can be grouped into those that seek to identify patient subtypes (studies of heterogeneity) or those that seek to identify gene signatures with prognostic or predictive capacity. The sheer number of reported signatures has led to speculation that everything is prognostic in BC. Here, we show that this ubiquity is an apparition caused by a poor understanding of the interrelatedness between subtype and the molecular determinants of prognosis. Our approach constructively shows how to avoid confounding due to a patient's subtype, clinicopathological profile, or treatment profile. The approach identifies patients who are predicted to have good outcome at time of diagnosis by all available clinical and molecular markers but who experience a distant metastasis within 5 years. These inherently difficult patients (~7% of BC) are prioritized for investigations of intratumoral heterogeneity.

Autocrine Activation of the Wnt/ß-Catenin Pathway by CUX1 and GLIS1 in Breast Cancers.

Autocrine activation of the Wnt/ß-catenin pathway occurs in several cancers, notably in breast tumors, and is associated with higher expression of various Wnt ligands. Using various inhibitors of the FZD/LRP receptor complex, we demonstrate that some adenosquamous carcinomas that develop in MMTV-CUX1 transgenic mice represent a model for autocrine activation of the Wnt/ß-catenin pathway. By comparing expression profiles of laser-capture microdissected mammary tumors, we identify Glis1 as a transcription factor that is highly expressed in the subset of tumors with elevated Wnt gene expression. Analysis of human cancer datasets confirms that elevated WNT gene expression is associated with high levels of CUX1 and GLIS1 and correlates with genes of the epithelial-to-mesenchymal transition (EMT) signature: VIM, SNAI1 and TWIST1 are elevated whereas CDH1 and OCLN are decreased. Co-expression experiments demonstrate that CUX1 and GLIS1 cooperate to stimulate TCF/ß-catenin transcriptional activity and to enhance cell migration and invasion. Altogether, these results provide additional evidence for the role of GLIS1 in reprogramming gene expression and suggest a hierarchical model for transcriptional regulation of the Wnt/ß-catenin pathway and the epithelial-to-mesenchymal transition.

Ets2 in tumor fibroblasts promotes angiogenesis in breast cancer.

Tumor fibroblasts are active partners in tumor progression, but the genes and pathways that mediate this collaboration are ill-defined. Previous work demonstrates that Ets2 function in stromal cells significantly contributes to breast tumor progression. Conditional mouse models were used to study the function of Ets2 in both mammary stromal fibroblasts and epithelial cells. Conditional inactivation of Ets2 in stromal fibroblasts in PyMT and ErbB2 driven tumors significantly reduced tumor growth, however deletion of Ets2 in epithelial cells in the PyMT model had no significant effect. Analysis of gene expression in fibroblasts revealed a tumor- and Ets2-dependent gene signature that was enriched in genes important for ECM remodeling, cell migration, and angiogenesis in both PyMT and ErbB2 driven-tumors. Consistent with these results, PyMT and ErbB2 tumors lacking Ets2 in fibroblasts had fewer functional blood vessels, and Ets2 in fibroblasts elicited changes in gene expression in tumor endothelial cells consistent with this phenotype. An in vivo angiogenesis assay revealed the ability of Ets2 in fibroblasts to promote blood vessel formation in the absence of tumor cells. Importantly, the Ets2-dependent gene expression signatures from both mouse models were able to distinguish human breast tumor stroma from normal stroma, and correlated with patient outcomes in two whole tumor breast cancer data sets. The data reveals a key function for Ets2 in tumor fibroblasts in signaling to endothelial cells to promote tumor angiogenesis. The results highlight the collaborative networks that orchestrate communication between stromal cells and tumor cells, and suggest that targeting tumor fibroblasts may be an effective strategy for developing novel anti-angiogenic therapies.

ß-Catenin signaling is a critical event in ErbB2-mediated mammary tumor progression.

Although ERBB2 amplification and overexpression is correlated with poor outcome in breast cancer, the molecular mechanisms underlying the aggressive nature of these tumors has not been fully elucidated. To investigate this further, we have used a transgenic mouse model of ErbB2-driven tumor progression (ErbB2(KI) model) that recapitulates clinically relevant events, including selective amplification of the core erbB2 amplicon. By comparing the transcriptional profiles of ErbB2(KI) mammary tumors and human ERBB2-positive breast cancers, we show that ErbB2(KI) tumors possess molecular features of the basal subtype of ERBB2-positive human breast cancer, including activation of canonical ß-catenin signaling. Inhibition of ß-catenin-dependent signaling in ErbB2(KI)-derived tumor cells using RNA interference impaired tumor initiation and metastasis. Furthermore, treatment of ErbB2(KI) or human ERBB2-overexpressing tumor cells with a selective ß-catenin/CBP inhibitor significantly decreased proliferation and ErbB2 expression. Collectively, our data indicate that ERBB2-mediated breast cancer progression requires ß-catenin signaling and can be therapeutically targeted by selective ß-catenin/CBP inhibitors.

Met synergizes with p53 loss to induce mammary tumors that possess features of claudin-low breast cancer.

Triple-negative breast cancer (TNBC) accounts for ∼20% of cases and contributes to basal and claudin-low molecular subclasses of the disease. TNBCs have poor prognosis, display frequent mutations in tumor suppressor gene p53 (TP53), and lack targeted therapies. The MET receptor tyrosine kinase is elevated in TNBC and transgenic Met models (Met(mt)) develop basal-like tumors. To investigate collaborating events in the genesis of TNBC, we generated Met(mt) mice with conditional loss of murine p53 (Trp53) in mammary epithelia. Somatic Trp53 loss, in combination with Met(mt), significantly increased tumor penetrance over Met(mt) or Trp53 loss alone. Unlike Met(mt) tumors, which are histologically diverse and enriched in a basal-like molecular signature, the majority of Met(mt) tumors with Trp53 loss displayed a spindloid pathology with a distinct molecular signature that resembles the human claudin-low subtype of TNBC, including diminished claudins, an epithelial-to-mesenchymal transition signature, and decreased expression of the microRNA-200 family. Moreover, although mammary specific loss of Trp53 promotes tumors with diverse pathologies, those with spindloid pathology and claudin-low signature display genomic Met amplification. In both models, MET activity is required for maintenance of the claudin-low morphological phenotype, in which MET inhibitors restore cell-cell junctions, rescue claudin 1 expression, and abrogate growth and dissemination of cells in vivo. Among human breast cancers, elevated levels of MET and stabilized TP53, indicative of mutation, correlate with highly proliferative TNBCs of poor outcome. This work shows synergy between MET and TP53 loss for claudin-low breast cancer, identifies a restricted claudin-low gene signature, and provides a rationale for anti-MET therapies in TNBC.

Met induces mammary tumors with diverse histologies and is associated with poor outcome and human basal breast cancer.

Elevated MET receptor tyrosine kinase correlates with poor outcome in breast cancer, yet the reasons for this are poorly understood. We thus generated a transgenic mouse model targeting expression of an oncogenic Met receptor (Met(mt)) to the mammary epithelium. We show that Met(mt) induces mammary tumors with multiple phenotypes. These reflect tumor subtypes with gene expression and immunostaining profiles sharing similarities to human basal and luminal breast cancers. Within the basal subtype, Met(mt) induces tumors with signatures of WNT and epithelial to mesenchymal transition (EMT). Among human breast cancers, MET is primarily elevated in basal and ERBB2-positive subtypes with poor prognosis, and we show that MET, together with EMT marker, SNAIL, are highly predictive of poor prognosis in lymph node-negative patients. By generating a unique mouse model in which the Met receptor tyrosine kinase is expressed in the mammary epithelium, along with the examination of MET expression in human breast cancer, we have established a specific link between MET and basal breast cancer. This work identifies basal breast cancers and, additionally, poor-outcome breast cancers, as those that may benefit from anti-MET receptor therapies.