Diagnostic Outcomes of Concurrent DNA and RNA Sequencing in Individuals Undergoing Hereditary Cancer Testing.

Importance: Personalized surveillance, prophylaxis, and cancer treatment options for individuals with hereditary cancer predisposition are informed by results of germline genetic testing. Improvements to genomic technology, such as the availability of RNA sequencing, may increase identification of individuals eligible for personalized interventions by improving the accuracy and yield of germline testing. Objective: To assess the cumulative association of paired DNA and RNA testing with detection of disease-causing germline genetic variants and resolution of variants of uncertain significance (VUS). Design, Setting, and Participants: Paired DNA and RNA sequencing was performed on individuals undergoing germline testing for hereditary cancer indication at a single diagnostic laboratory from March 2019 through April 2020. Demographic characteristics, clinical data, and test results were curated as samples were received, and changes to variant classification were assessed over time. Data analysis was performed from May 2020 to June 2023. Main Outcomes and Measures: Main outcomes were increase in diagnostic yield, decrease in VUS rate, the overall results by variant type, the association of RNA evidence with variant classification, and the corresponding predicted effect on cancer risk management. Results: A total of 43 524 individuals were included (median [range] age at testing, 54 [2-101] years; 37 373 female individuals [85.7%], 6224 male individuals [14.3%], and 2 individuals of unknown sex [<0.1%]), with 43 599 tests. A total of 2197 (5.0%) were Ashkenazi Jewish, 1539 (3.5%) were Asian, 3077 (7.1%) were Black, 2437 (5.6%) were Hispanic, 27 793 (63.7%) were White, and 2049 (4.7%) were other race, and for 4507 individuals (10.3%), race and ethnicity were unknown. Variant classification was impacted in 549 individuals (1.3%). Medically significant upgrades were made in 97 individuals, including 70 individuals who had a variant reclassified from VUS to pathogenic/likely pathogenic (P/LP) and 27 individuals who had a novel deep intronic P/LP variant that would not have been detected using DNA sequencing alone. A total of 93 of 545 P/LP splicing variants (17.1%) were dependent on RNA evidence for classification, and 312 of 439 existing splicing VUS (71.1%) were resolved by RNA evidence. Notably, the increase in positive rate (3.1%) and decrease in VUS rate (-3.9%) was higher in Asian, Black, and Hispanic individuals combined compared to White individuals (1.6%; P = .02; and -2.5%; P < .001). Conclusions and Relevance: Findings of this diagnostic study demonstrate that the ability to perform RNA sequencing concurrently with DNA sequencing represents an important advancement in germline genetic testing by improving detection of novel variants and classification of existing variants. This expands the identification of individuals with hereditary cancer predisposition and increases opportunities for personalization of therapeutics and surveillance.

Mutational and splicing landscape in a cohort of 43,000 patients tested for hereditary cancer.

DNA germline genetic testing can identify individuals with cancer susceptibility. However, DNA sequencing alone is limited in its detection and classification of mRNA splicing variants, particularly those located far from coding sequences. Here we address the limitations of splicing variant identification and interpretation by pairing DNA and RNA sequencing and describe the mutational and splicing landscape in a clinical cohort of 43,524 individuals undergoing genetic testing for hereditary cancer predisposition.

Splicing profile by capture RNA-seq identifies pathogenic germline variants in tumor suppressor genes.

Germline variants in tumor suppressor genes (TSGs) can result in RNA mis-splicing and predisposition to cancer. However, identification of variants that impact splicing remains a challenge, contributing to a substantial proportion of patients with suspected hereditary cancer syndromes remaining without a molecular diagnosis. To address this, we used capture RNA-sequencing (RNA-seq) to generate a splicing profile of 18 TSGs (APC, ATM, BRCA1, BRCA2, BRIP1, CDH1, CHEK2, MLH1, MSH2, MSH6, MUTYH, NF1, PALB2, PMS2, PTEN, RAD51C, RAD51D, and TP53) in 345 whole-blood samples from healthy donors. We subsequently demonstrated that this approach can detect mis-splicing by comparing splicing profiles from the control dataset to profiles generated from whole blood of individuals previously identified with pathogenic germline splicing variants in these genes. To assess the utility of our TSG splicing profile to prospectively identify pathogenic splicing variants, we performed concurrent capture DNA and RNA-seq in a cohort of 1000 patients with suspected hereditary cancer syndromes. This approach improved the diagnostic yield in this cohort, resulting in a 9.1% relative increase in the detection of pathogenic variants, demonstrating the utility of performing simultaneous DNA and RNA genetic testing in a clinical context.

Transposable element expression in tumors is associated with immune infiltration and increased antigenicity.

Profound global loss of DNA methylation is a hallmark of many cancers. One potential consequence of this is the reactivation of transposable elements (TEs) which could stimulate the immune system via cell-intrinsic antiviral responses. Here, we develop REdiscoverTE, a computational method for quantifying genome-wide TE expression in RNA sequencing data. Using The Cancer Genome Atlas database, we observe increased expression of over 400 TE subfamilies, of which 262 appear to result from a proximal loss of DNA methylation. The most recurrent TEs are among the evolutionarily youngest in the genome, predominantly expressed from intergenic loci, and associated with antiviral or DNA damage responses. Treatment of glioblastoma cells with a demethylation agent results in both increased TE expression and de novo presentation of TE-derived peptides on MHC class I molecules. Therapeutic reactivation of tumor-specific TEs may synergize with immunotherapy by inducing inflammation and the display of potentially immunogenic neoantigens.

Wnt11 regulates cardiac chamber development and disease during perinatal maturation.

Ventricular chamber growth and development during perinatal circulatory transition is critical for functional adaptation of the heart. However, the chamber-specific programs of neonatal heart growth are poorly understood. We used integrated systems genomic and functional biology analyses of the perinatal chamber specific transcriptome and we identified Wnt11 as a prominent regulator of chamber-specific proliferation. Importantly, downregulation of Wnt11 expression was associated with cyanotic congenital heart defect (CHD) phenotypes and correlated with O2 saturation levels in hypoxemic infants with Tetralogy of Fallot (TOF). Perinatal hypoxia treatment in mice suppressed Wnt11 expression and induced myocyte proliferation more robustly in the right ventricle, modulating Rb1 protein activity. Wnt11 inactivation was sufficient to induce myocyte proliferation in perinatal mouse hearts and reduced Rb1 protein and phosphorylation in neonatal cardiomyocytes. Finally, downregulated Wnt11 in hypoxemic TOF infantile hearts was associated with Rb1 suppression and induction of proliferation markers. This study revealed a previously uncharacterized function of Wnt11-mediated signaling as an important player in programming the chamber-specific growth of the neonatal heart. This function influences the chamber-specific development and pathogenesis in response to hypoxia and cyanotic CHDs. Defining the underlying regulatory mechanism may yield chamber-specific therapies for infants born with CHDs.

Recurrent patterns of DNA copy number alterations in tumors reflect metabolic selection pressures.

Copy number alteration (CNA) profiling of human tumors has revealed recurrent patterns of DNA amplifications and deletions across diverse cancer types. These patterns are suggestive of conserved selection pressures during tumor evolution but cannot be fully explained by known oncogenes and tumor suppressor genes. Using a pan-cancer analysis of CNA data from patient tumors and experimental systems, here we show that principal component analysis-defined CNA signatures are predictive of glycolytic phenotypes, including 18F-fluorodeoxy-glucose (FDG) avidity of patient tumors, and increased proliferation. The primary CNA signature is enriched for p53 mutations and is associated with glycolysis through coordinate amplification of glycolytic genes and other cancer-linked metabolic enzymes. A pan-cancer and cross-species comparison of CNAs highlighted 26 consistently altered DNA regions, containing 11 enzymes in the glycolysis pathway in addition to known cancer-driving genes. Furthermore, exogenous expression of hexokinase and enolase enzymes in an experimental immortalization system altered the subsequent copy number status of the corresponding endogenous loci, supporting the hypothesis that these metabolic genes act as drivers within the conserved CNA amplification regions. Taken together, these results demonstrate that metabolic stress acts as a selective pressure underlying the recurrent CNAs observed in human tumors, and further cast genomic instability as an enabling event in tumorigenesis and metabolic evolution.

Decoding the Long Noncoding RNA During Cardiac Maturation: A Roadmap for Functional Discovery.

BACKGROUND: Cardiac maturation during perinatal transition of heart is critical for functional adaptation to hemodynamic load and nutrient environment. Perturbation in this process has major implications in congenital heart defects. Transcriptome programming during perinatal stages is an important information but incomplete in current literature, particularly, the expression profiles of the long noncoding RNAs (lncRNAs) are not fully elucidated. METHODS AND RESULTS: From comprehensive analysis of transcriptomes derived from neonatal mouse heart left and right ventricles, a total of 45 167 unique transcripts were identified, including 21 916 known and 2033 novel lncRNAs. Among these lncRNAs, 196 exhibited significant dynamic regulation along maturation process. By implementing parallel weighted gene co-expression network analysis of mRNA and lncRNA data sets, several lncRNA modules coordinately expressed in a developmental manner similar to protein coding genes, while few lncRNAs revealed chamber-specific patterns. Out of 2262 lncRNAs located within 50 kb of protein coding genes, 5% significantly correlate with the expression of their neighboring genes. The impact of Ppp1r1b-lncRNA on the corresponding partner gene Tcap was validated in cultured myoblasts. This concordant regulation was also conserved in human infantile hearts. Furthermore, the Ppp1r1b-lncRNA/Tcap expression ratio was identified as a molecular signature that differentiated congenital heart defect phenotypes. CONCLUSIONS: The study provides the first high-resolution landscape on neonatal cardiac lncRNAs and reveals their potential interaction with mRNA transcriptome during cardiac maturation. Ppp1r1b-lncRNA was identified as a regulator of Tcap expression, with dynamic interaction in postnatal cardiac development and congenital heart defects.

The Sm protein methyltransferase PRMT5 is not required for primordial germ cell specification in mice.

PRMT5 is a type II protein arginine methyltransferase with roles in stem cell biology, reprograming, cancer and neurogenesis. During embryogenesis in the mouse, it was hypothesized that PRMT5 functions with the master germline determinant BLIMP1 to promote primordial germ cell (PGC) specification. Using a Blimp1-Cre germline conditional knockout, we discovered that Prmt5 has no major role in murine germline specification, or the first global epigenetic reprograming event involving depletion of cytosine methylation from DNA and histone H3 lysine 9 dimethylation from chromatin. Instead, we discovered that PRMT5 functions at the conclusion of PGC reprograming I to promote proliferation, survival and expression of the gonadal germline program as marked by MVH. We show that PRMT5 regulates gene expression by promoting methylation of the Sm spliceosomal proteins and significantly altering the spliced repertoire of RNAs in mammalian embryonic cells and primordial cells.

COX-2 inhibition prevents the appearance of cutaneous squamous cell carcinomas accelerated by BRAF inhibitors.

Keratoacanthomas (KAs) and cutaneous squamous cell carcinomas (cuSCCs) develop in 15-30% of patients with BRAF(V600E) metastatic melanoma treated with BRAF inhibitors (BRAFi). These lesions resemble mouse skin tumors induced by the two-stage DMBA/TPA skin carcinogenesis protocol; in this protocol BRAFi accelerates tumor induction. Since prior studies demonstrated cyclooxygenase 2 (COX-2) is necessary for DMBA/TPA tumor induction, we hypothesized that COX-2 inhibition might prevent BRAFi-accelerated skin tumors. Celecoxib, a COX-2 inhibitor, significantly delayed tumor acceleration by the BRAFi inhibitor PLX7420 and decreased tumor number by 90%. Tumor gene expression profiling demonstrated that celecoxib partially reversed the PLX4720-induced gene signature. In PDV cuSCC cells, vemurafenib (a clinically approved BRAFi) increased ERK phosphorylation and soft agar colony formation; both responses were greatly decreased by celecoxib. In clinical trials trametinib, a MEK inhibitor (MEKi) increases BRAFi therapy efficacy in BRAF(V600E) melanomas and reduces BRAFi-induced KA and cuSCC frequency. Trametinib also reduced vemurafenib-induced PDV soft agar colonies, but less efficiently than celecoxib. The trametinb/celecoxib combination was more effective than either inhibitor alone. In conclusion, celecoxib suppressed both BRAFi-accelerated skin tumors and soft-agar colonies, warranting its testing as a chemopreventive agent for non-melanoma skin lesions in patients treated with BRAFi alone or in combination with MEKi.

Doxycycline alters metabolism and proliferation of human cell lines.

The tetracycline antibiotics are widely used in biomedical research as mediators of inducible gene expression systems. Despite many known effects of tetracyclines on mammalian cells-including inhibition of the mitochondrial ribosome-there have been few reports on potential off-target effects at concentrations commonly used in inducible systems. Here, we report that in human cell lines, commonly used concentrations of doxycycline change gene expression patterns and concomitantly shift metabolism towards a more glycolytic phenotype, evidenced by increased lactate secretion and reduced oxygen consumption. We also show that these concentrations are sufficient to slow proliferation. These findings suggest that researchers using doxycycline in inducible expression systems should design appropriate controls to account for potential confounding effects of the drug on cellular metabolism.

Host immunity contributes to the anti-melanoma activity of BRAF inhibitors.

The BRAF mutant, BRAF(V600E), is expressed in nearly half of melanomas, and oral BRAF inhibitors induce substantial tumor regression in patients with BRAF(V600E) metastatic melanoma. The inhibitors are believed to work primarily by inhibiting BRAF(V600E)-induced oncogenic MAPK signaling; however, some patients treated with BRAF inhibitors exhibit increased tumor immune infiltration, suggesting that a combination of BRAF inhibitors and immunotherapy may be beneficial. We used two relatively resistant variants of Braf(V600E)-driven mouse melanoma (SM1 and SM1WT1) and melanoma-prone mice to determine the role of host immunity in type I BRAF inhibitor PLX4720 antitumor activity. We found that PLX4720 treatment downregulated tumor Ccl2 gene expression and decreased tumor CCL2 expression in both Braf(V600E) mouse melanoma transplants and in de novo melanomas in a manner that was coincident with reduced tumor growth. While PLX4720 did not directly increase tumor immunogenicity, analysis of SM1 tumor-infiltrating leukocytes in PLX4720-treated mice demonstrated a robust increase in CD8(+) T/FoxP3(+)CD4(+) T cell ratio and NK cells. Combination therapy with PLX4720 and anti-CCL2 or agonistic anti-CD137 antibodies demonstrated significant antitumor activity in mouse transplant and de novo tumorigenesis models. These data elucidate a role for host CCR2 in the mechanism of action of type I BRAF inhibitors and support the therapeutic potential of combining BRAF inhibitors with immunotherapy.