Alternative RNA splicing defects in pediatric cancers: new insights in tumorigenesis and potential therapeutic vulnerabilities.

BACKGROUND: Compared with adult cancers, pediatric cancers are uniquely characterized by a genomically stable landscape and lower tumor mutational burden. Alternative splicing, however, a global cellular process that produces different messenger RNA/protein isoforms from a single messenger RNA transcript, has been increasingly implicated in the development of pediatric cancers. DESIGN: We review the current literature on the role of alternative splicing in adult cancer, cancer predisposition syndromes, and pediatric cancers. We also describe multiple splice variants identified in adult cancers and confirmed through comprehensive genomic profiling in our institutional cohort of rare, refractory, and relapsed pediatric and adolescent young adult cancer patients. Finally, we summarize the contributions of alternative splicing events to neoantigens and chemoresistance and prospects for splicing-based therapies. RESULTS: Published dysregulated splicing events can be categorized as exon inclusion, exon exclusion, splicing factor up-regulation, or splice site alterations. We observe these phenomena in cancer predisposition syndromes (Lynch syndrome, Li-Fraumeni syndrome, CHEK2) and pediatric leukemia (B-cell acute lymphoblastic leukemia), sarcomas (Ewing sarcoma, rhabdomyosarcoma, osteosarcoma), retinoblastoma, Wilms' tumor, and neuroblastoma. Within our institutional cohort, we demonstrate splice variants in key regulatory genes (CHEK2, TP53, PIK3R1, MDM2, KDM6A, NF1) that resulted in exon exclusion or splice site alterations, which were predicted to impact functional protein expression and promote tumorigenesis. Differentially spliced isoforms and splicing proteins also impact neoantigen creation and treatment resistance, such as imatinib or glucocorticoid regimens. Additionally, splice-altering strategies with the potential to change the therapeutic landscape of pediatric cancers include antisense oligonucleotides, adeno-associated virus gene transfers, and small molecule inhibitors. CONCLUSIONS: Alternative splicing plays a critical role in the formation and growth of pediatric cancers, and our institutional cohort confirms and highlights the broad spectrum of affected genes in a variety of cancers. Further studies that elucidate the mechanisms of disease-inducing splicing events will contribute toward the development of novel therapeutics.

Genomic characterization of HER2-positive breast cancer and response to neoadjuvant trastuzumab and chemotherapy-results from the ACOSOG Z1041 (Alliance) trial.

Background: HER2 (ERBB2) gene amplification and its corresponding overexpression are present in 15-30% of invasive breast cancers. While HER2-targeted agents are effective treatments, resistance remains a major cause of death. The American College of Surgeons Oncology Group Z1041 trial (NCT00513292) was designed to compare the pathologic complete response (pCR) rate of distinct regimens of neoadjuvant chemotherapy and trastuzumab, but ultimately identified no difference. Patients and methods: In supplement to tissues from 37 Z1041 cases, 11 similarly treated cases were obtained from a single institution study (NCT00353483). We have extracted genomic DNA from both pre-treatment tumor biopsies and blood of these 48 cases, and performed whole genome (WGS) and exome sequencing. Coincident with these efforts, we have generated RNA-seq profiles from 42 of the tumor biopsies. Among patients in this cohort, 24 (50%) achieved a pCR. Results: We have characterized the genomic landscape of HER2-positive breast cancer and investigated associations between genomic features and pCR. Cases assigned to the HER2-enriched subtype by RNA-seq analysis were more likely to achieve a pCR compared to the luminal, basal-like, or normal-like subtypes (19/27 versus 3/15; P = 0.0032). Mutational events led to the generation of putatively active neoantigens, but were overall not associated with pCR. ERBB2 and GRB7 were the genes most commonly observed in fusion events, and genomic copy number analysis of the ERBB2 locus indicated that cases with either no observable or low-level ERBB2 amplification were less likely to achieve a pCR (7/8 versus 17/40; P = 0.048). Moreover, among cases that achieved a pCR, tumors consistently expressed immune signatures that may contribute to therapeutic response. Conclusion: The identification of these features suggests that it may be possible to predict, at the time of diagnosis, those HER2-positive breast cancer patients who will not respond to treatment with chemotherapy and trastuzumab. ClinicalTrials.gov identifiers: NCT00513292, NCT00353483.

Consensus on precision medicine for metastatic cancers: a report from the MAP conference.

Recent advances in biotechnologies have led to the development of multiplex genomic and proteomic analyses for clinical use. Nevertheless, guidelines are currently lacking to determine which molecular assays should be implemented in metastatic cancers. The first MAP conference was dedicated to exploring the use of genomics to better select therapies in the treatment of metastatic cancers. Sixteen consensus items were covered. There was a consensus that new technologies like next-generation sequencing of tumors and ddPCR on circulating free DNA have convincing analytical validity. Further work needs to be undertaken to establish the clinical utility of liquid biopsies and the added clinical value of expanding from individual gene tests into large gene panels. Experts agreed that standardized bioinformatics methods for biological interpretation of genomic data are needed and that precision medicine trials should be stratified based on the level of evidence available for the genomic alterations identified.

A genomic case study of mixed fibrolamellar hepatocellular carcinoma.

BACKGROUND: Mixed fibrolamellar hepatocellular carcinoma (mFL-HCC) is a rare liver tumor defined by the presence of both pure FL-HCC and conventional HCC components, represents up to 25% of cases of FL-HCC, and has been associated with worse prognosis. Recent genomic characterization of pure FL-HCC identified a highly recurrent transcript fusion (DNAJB1:PRKACA) not found in conventional HCC. PATIENTS AND METHODS: We performed exome and transcriptome sequencing of a case of mFL-HCC. A novel BAC-capture approach was developed to identify a 400 kb deletion as the underlying genomic mechanism for a DNAJB1:PRKACA fusion in this case. A sensitive Nanostring Elements assay was used to screen for this transcript fusion in a second case of mFL-HCC, 112 additional HCC samples and 44 adjacent non-tumor liver samples. RESULTS: We report the first comprehensive genomic analysis of a case of mFL-HCC. No common HCC-associated mutations were identified. The very low mutation rate of this case, large number of mostly single-copy, long-range copy number variants, and high expression of ERBB2 were more consistent with previous reports of pure FL-HCC than conventional HCC. In particular, the DNAJB1:PRKACA fusion transcript specifically associated with pure FL-HCC was detected at very high expression levels. Subsequent analysis revealed the presence of this fusion in all primary and metastatic samples, including those with mixed or conventional HCC pathology. A second case of mFL-HCC confirmed our finding that the fusion was detectable in conventional components. An expanded screen identified a third case of fusion-positive HCC, which upon review, also had both conventional and fibrolamellar features. This screen confirmed the absence of the fusion in all conventional HCC and adjacent non-tumor liver samples. CONCLUSION: These results indicate that mFL-HCC is similar to pure FL-HCC at the genomic level and the DNAJB1:PRKACA fusion can be used as a diagnostic tool for both pure and mFL-HCC.

Predicting response and survival in chemotherapy-treated triple-negative breast cancer.

BACKGROUND: In this study, we evaluated the ability of gene expression profiles to predict chemotherapy response and survival in triple-negative breast cancer (TNBC). METHODS: Gene expression and clinical-pathological data were evaluated in five independent cohorts, including three randomised clinical trials for a total of 1055 patients with TNBC, basal-like disease (BLBC) or both. Previously defined intrinsic molecular subtype and a proliferation signature were determined and tested. Each signature was tested using multivariable logistic regression models (for pCR (pathological complete response)) and Cox models (for survival). Within TNBC, interactions between each signature and the basal-like subtype (vs other subtypes) for predicting either pCR or survival were investigated. RESULTS: Within TNBC, all intrinsic subtypes were identified but BLBC predominated (55-81%). Significant associations between genomic signatures and response and survival after chemotherapy were only identified within BLBC and not within TNBC as a whole. In particular, high expression of a previously identified proliferation signature, or low expression of the luminal A signature, was found independently associated with pCR and improved survival following chemotherapy across different cohorts. Significant interaction tests were only obtained between each signature and the BLBC subtype for prediction of chemotherapy response or survival. CONCLUSIONS: The proliferation signature predicts response and improved survival after chemotherapy, but only within BLBC. This highlights the clinical implications of TNBC heterogeneity, and suggests that future clinical trials focused on this phenotypic subtype should consider stratifying patients as having BLBC or not.

Prioritizing targets for precision cancer medicine.

The implementation of cancer genomic testing into the clinical setting has brought major opportunities. However, as our understanding of cancer initiation, maintenance and progression improves through detailed cancer genomic studies, the challenges associated with driver identification and target classification in the clinical setting become clearer. Here, we review recent insights into cancer genomic testing in the clinical setting, and suggest a target classification approach that considers the levels of evidence supporting the prioritization of tumour drivers for therapeutic targeting in light of complex cancer clonal and sub-clonal structures and clinical successes and failures in the field. We argue that such classification approaches, together with transparent reporting of both positive and negative clinical data and continued research to identify the sub-clonal dynamics of driver events during the disease course, will facilitate inter-trial comparisons, optimize patient informed consent and provide a critically balanced evaluation of genomic testing in clinical practice.

An encyclopedia of mouse genes.

The laboratory mouse is the premier model system for studies of mammalian development due to the powerful classical genetic analysis possible (see also the Jackson Laboratory web site, http://www.jax.org/) and the ever-expanding collection of molecular tools. To enhance the utility of the mouse system, we initiated a program to generate a large database of expressed sequence tags (ESTs) that can provide rapid access to genes. Of particular significance was the possibility that cDNA libraries could be prepared from very early stages of development, a situation unrealized in human EST projects. We report here the development of a comprehensive database of ESTs for the mouse. The project, initiated in March 1996, has focused on 5' end sequences from directionally cloned, oligo-dT primed cDNA libraries. As of 23 October 1998, 352,040 sequences had been generated, annotated and deposited in dbEST, where they comprised 93% of the total ESTs available for mouse. EST data are versatile and have been applied to gene identification, comparative sequence analysis, comparative gene mapping and candidate disease gene identification, genome sequence annotation, microarray development and the development of gene-based map resources.

CTCF genetic alterations in endometrial carcinoma are pro-tumorigenic.

CTCF is a haploinsufficient tumour suppressor gene with diverse normal functions in genome structure and gene regulation. However the mechanism by which CTCF haploinsufficiency contributes to cancer development is not well understood. CTCF is frequently mutated in endometrial cancer. Here we show that most CTCF mutations effectively result in CTCF haploinsufficiency through nonsense-mediated decay of mutant transcripts, or loss-of-function missense mutation. Conversely, we identified a recurrent CTCF mutation K365T, which alters a DNA binding residue, and acts as a gain-of-function mutation enhancing cell survival. CTCF genetic deletion occurs predominantly in poor prognosis serous subtype tumours, and this genetic deletion is associated with poor overall survival. In addition, we have shown that CTCF haploinsufficiency also occurs in poor prognosis endometrial clear cell carcinomas and has some association with endometrial cancer relapse and metastasis. Using shRNA targeting CTCF to recapitulate CTCF haploinsufficiency, we have identified a novel role for CTCF in the regulation of cellular polarity of endometrial glandular epithelium. Overall, we have identified two novel pro-tumorigenic roles (promoting cell survival and altering cell polarity) for genetic alterations of CTCF in endometrial cancer.

Automated methods for single-stranded DNA isolation and dideoxynucleotide DNA sequencing reactions on a robotic workstation.

Automated procedures have been developed for both the simultaneous isolation of 96 single-stranded M13 chimeric template DNAs in less than two hours, and for simultaneously pipetting 24 dideoxynucleotide sequencing reactions on a commercially available laboratory workstation. The DNA sequencing results obtained by either radiolabeled or fluorescent methods are consistent with the premise that automation of these portions of DNA sequencing projects will improve the reproducibility of the DNA isolation and the procedures for these normally labor-intensive steps provides an approach for rapid acquisition of large amounts of high quality, reproducible DNA sequence data.

Prep-A-Gene: a superior matrix for the purification of DNA and DNA fragments.

A new chromatographic matrix, Prep-A-Gene, is described for the isolation and purification of high quality DNA suitable for restriction analysis, ligation, transformation and sequencing protocols. This matrix selectively binds DNA greater than approximately 200 base pairs in length, while RNA, proteins, cellular components, agarose and other contaminants are washed free in minutes. This eliminates the need for time-consuming and laborious RNase treatments, gel extractions and phenol extractions. The DNA that is desorbed from the matrix is available immediately as a substrate for subsequent protocols. DNA purified in this manner exhibits no detectable shearing, even with more fragile chromosomal DNA.

Capillary electrophoresis platforms for DNA sequence analysis.

Commercial availability of multicapillary DNA sequencing instruments can significantly impact the existing paradigm in high-throughput DNA sequencing facilities, shifting the rate-limiting step away from the electrophoresis and detection step. These instruments, although novel and poorly understood, represent a significant capital investment.This review informs the reader about the state-of-the-art instruments, although capillary sequencers are in a state of almost constant flux, with improvements emerging rapidly.

Association between DQB1 and cervical cancer in patients with human papillomavirus and family controls.

OBJECTIVE: The role of human leukocyte antigen (HLA) DQB1 alleles and human papillomavirus (HPV) as contributing factors to invasive cervical cancer was investigated. To overcome problems of misleading causal inferences common in traditional case-control studies, a family-based test, the transmission/disequilibrium test, was used. METHODS: Ninety-six patients with pathologically confirmed invasive cervical cancer were ascertained. Human papillomavirus types were determined in 80 patients, of whom 81.25% were HPV-positive, and 18.75% were HPV-negative. Deoxyribonucleic acid was extracted from samples, taken from patients and their parents, and sequenced to determine DQB1 genotypes. Nuclear family data were used to test whether the DQB1 locus is associated with invasive cervical cancer while controlling for high-risk HPV-positive patients. The transmission/disequilibrium test evaluates whether the frequency of transmission of parental marker alleles to their affected offspring deviates from the expected Mendelian frequency of 50%. RESULTS: The HLA DQB1 locus showed evidence for allelic association with invasive cervical cancer in high-risk HPV-positive patients (P = .006). The transmission/disequilibrium test showed that the DQB1*0303 allele was transmitted to high-risk HPV patients more often than expected by chance, chi2(1) = 8.0, P = .005 (P = .035 when correcting for multiple tests). Tests of association were negative when applied to all 96 patients, irrespective of HPV status. No significant differences were found in the distribution of the DQB1 alleles among HPV-positive patients compared with those who were HPV-negative, indicating that HLA alleles are not associated with susceptibility to HPV infection. CONCLUSION: These results suggest that the DQB1*0303 allele increases the risk for invasive cervical cancer in women who are HPV-positive.

Clonal diversity of recurrently mutated genes in myelodysplastic syndromes.

Recent studies suggest that most cases of myelodysplastic syndrome (MDS) are clonally heterogeneous, with a founding clone and multiple subclones. It is not known whether specific gene mutations typically occur in founding clones or subclones. We screened a panel of 94 candidate genes in a cohort of 157 patients with MDS or secondary acute myeloid leukemia (sAML). This included 150 cases with samples obtained at MDS diagnosis and 15 cases with samples obtained at sAML transformation (8 were also analyzed at the MDS stage). We performed whole-genome sequencing (WGS) to define the clonal architecture in eight sAML genomes and identified the range of variant allele frequencies (VAFs) for founding clone mutations. At least one mutation or cytogenetic abnormality was detected in 83% of the 150 MDS patients and 17 genes were significantly mutated (false discovery rate ≤0.05). Individual genes and patient samples displayed a wide range of VAFs for recurrently mutated genes, indicating that no single gene is exclusively mutated in the founding clone. The VAFs of recurrently mutated genes did not fully recapitulate the clonal architecture defined by WGS, suggesting that comprehensive sequencing may be required to accurately assess the clonal status of recurrently mutated genes in MDS.

Recurrent DNMT3A mutations in patients with myelodysplastic syndromes.

Alterations in DNA methylation have been implicated in the pathogenesis of myelodysplastic syndromes (MDS), although the underlying mechanism remains largely unknown. Methylation of CpG dinucleotides is mediated by DNA methyltransferases, including DNMT1, DNMT3A and DNMT3B. DNMT3A mutations have recently been reported in patients with de novo acute myeloid leukemia (AML), providing a rationale for examining the status of DNMT3A in MDS samples. In this study, we report the frequency of DNMT3A mutations in patients with de novo MDS, and their association with secondary AML. We sequenced all coding exons of DNMT3A using DNA from bone marrow and paired normal cells from 150 patients with MDS and identified 13 heterozygous mutations with predicted translational consequences in 12/150 patients (8.0%). Amino acid R882, located in the methyltransferase domain of DNMT3A, was the most common mutation site, accounting for 4/13 mutations. DNMT3A mutations were expressed in the majority of cells in all tested mutant samples regardless of myeloblast counts, suggesting that DNMT3A mutations occur early in the course of MDS. Patients with DNMT3A mutations had worse overall survival compared with patients without DNMT3A mutations (P=0.005) and more rapid progression to AML (P=0.007), suggesting that DNMT3A mutation status may have prognostic value in de novo MDS.

End sequence determination from large insert clones using energy transfer fluorescent primers.

Genome mapping strategies depend heavily on confirmatory data of several types to establish overlaps between contiguous stretches of cloned DNA derived from genomic regions. One type of ancillary data that can contribute to establishing these overlaps is DNA sequence data derived from the ends of large (> 30 kb) inserts in genomic clones. This type of data can be difficult to obtain routinely, because large clones are often unstable and microgram quantities of highly purified DNA are required in each sequencing reaction to obtain sufficient signal for accurate base calling and maximum read length. Recently, we have been experimenting with methods to consistently obtain up to 800 bases of high-quality sequence data from the ends of large insert clones using ThermoSequenase DNA polymerase and Energy Transfer fluorescent primers. Our experimental approach and results, described in this paper, indicate that routinely obtaining high-quality sequence data from the ends of large insert genomic clones is feasible. Such data can contribute to the assessment of common regions between large insert clones, to the establishment of conservation of synteny between closely related species, and to the detection of additional contiguous clones.

Cloning and DNA sequence of the Mycobacterium fortuitum var fortuitum plasmid pAL5000.

The complete nucleotide sequence of the Mycobacterium fortuitum var fortuitum plasmid pAL5000 has been determined. Computer analysis of this 4821-bp plasmid for protein coding regions, based on mycobacterial codon usage preferences, reveals the presence of two putative protein coding regions immediately downstream from typical mycobacterial promoter and ribosome binding sites. Both open reading frames, ORF1 and ORF2, produced proteins with the predicted respective sizes previously shown in minicell expression experiments [A. Labidi et al. (1985) FEMS Microbiol. Lett. 30, 221-225]. ORF1 encodes a putative 20-kDa basic protein with characteristics of a DNA binding protein involved in plasmid DNA replication. ORF2 encodes a 67-kDa protein with an amino-terminal sequence suggestive of a transported protein and a possible transmembrane anchor near its carboxyl-terminal. The current sequence and its analysis are more consistent with the minicell expression experiments than the previously published sequence of the pAL5000 plasmid [J. Rauzier et al. (1988) Gene 71, 315-321].