Inhibition of YTHDF2 triggers proteotoxic cell death in MYC-driven breast cancer.

RNA-binding proteins (RBPs) are critical regulators of post-transcriptional gene expression, and aberrant RBP-RNA interactions can promote cancer progression. Here, we interrogate the function of RBPs in cancer using pooled CRISPR-Cas9 screening and identify 57 RBP candidates with distinct roles in supporting MYC-driven oncogenic pathways. We find that disrupting YTHDF2-dependent mRNA degradation triggers apoptosis in triple-negative breast cancer (TNBC) cells and tumors. eCLIP and m6A sequencing reveal that YTHDF2 interacts with mRNAs encoding proteins in the MAPK pathway that, when stabilized, induce epithelial-to-mesenchymal transition and increase global translation rates. scRibo-STAMP profiling of translating mRNAs reveals unique alterations in the translatome of single cells within YTHDF2-depleted solid tumors, which selectively contribute to endoplasmic reticulum stress-induced apoptosis in TNBC cells. Thus, our work highlights the therapeutic potential of RBPs by uncovering a critical role for YTHDF2 in counteracting the global increase of mRNA synthesis in MYC-driven breast cancers.

Spliceosome-targeted therapies trigger an antiviral immune response in triple-negative breast cancer.

Many oncogenic insults deregulate RNA splicing, often leading to hypersensitivity of tumors to spliceosome-targeted therapies (STTs). However, the mechanisms by which STTs selectively kill cancers remain largely unknown. Herein, we discover that mis-spliced RNA itself is a molecular trigger for tumor killing through viral mimicry. In MYC-driven triple-negative breast cancer, STTs cause widespread cytoplasmic accumulation of mis-spliced mRNAs, many of which form double-stranded structures. Double-stranded RNA (dsRNA)-binding proteins recognize these endogenous dsRNAs, triggering antiviral signaling and extrinsic apoptosis. In immune-competent models of breast cancer, STTs cause tumor cell-intrinsic antiviral signaling, downstream adaptive immune signaling, and tumor cell death. Furthermore, RNA mis-splicing in human breast cancers correlates with innate and adaptive immune signatures, especially in MYC-amplified tumors that are typically immune cold. These findings indicate that dsRNA-sensing pathways respond to global aberrations of RNA splicing in cancer and provoke the hypothesis that STTs may provide unexplored strategies to activate anti-tumor immune pathways.

Engineering and Functional Characterization of Fusion Genes Identifies Novel Oncogenic Drivers of Cancer.

Oncogenic gene fusions drive many human cancers, but tools to more quickly unravel their functional contributions are needed. Here we describe methodology permitting fusion gene construction for functional evaluation. Using this strategy, we engineered the known fusion oncogenes, BCR-ABL1, EML4-ALK, and ETV6-NTRK3, as well as 20 previously uncharacterized fusion genes identified in The Cancer Genome Atlas datasets. In addition to confirming oncogenic activity of the known fusion oncogenes engineered by our construction strategy, we validated five novel fusion genes involving MET, NTRK2, and BRAF kinases that exhibited potent transforming activity and conferred sensitivity to FDA-approved kinase inhibitors. Our fusion construction strategy also enabled domain-function studies of BRAF fusion genes. Our results confirmed other reports that the transforming activity of BRAF fusions results from truncation-mediated loss of inhibitory domains within the N-terminus of the BRAF protein. BRAF mutations residing within this inhibitory region may provide a means for BRAF activation in cancer, therefore we leveraged the modular design of our fusion gene construction methodology to screen N-terminal domain mutations discovered in tumors that are wild-type at the BRAF mutation hotspot, V600. We identified an oncogenic mutation, F247L, whose expression robustly activated the MAPK pathway and sensitized cells to BRAF and MEK inhibitors. When applied broadly, these tools will facilitate rapid fusion gene construction for subsequent functional characterization and translation into personalized treatment strategies. Cancer Res; 77(13); 3502-12. ©2017 AACR.

The spliceosome is a therapeutic vulnerability in MYC-driven cancer.

MYC (also known as c-MYC) overexpression or hyperactivation is one of the most common drivers of human cancer. Despite intensive study, the MYC oncogene remains recalcitrant to therapeutic inhibition. MYC is a transcription factor, and many of its pro-tumorigenic functions have been attributed to its ability to regulate gene expression programs. Notably, oncogenic MYC activation has also been shown to increase total RNA and protein production in many tissue and disease contexts. While such increases in RNA and protein production may endow cancer cells with pro-tumour hallmarks, this increase in synthesis may also generate new or heightened burden on MYC-driven cancer cells to process these macromolecules properly. Here we discover that the spliceosome is a new target of oncogenic stress in MYC-driven cancers. We identify BUD31 as a MYC-synthetic lethal gene in human mammary epithelial cells, and demonstrate that BUD31 is a component of the core spliceosome required for its assembly and catalytic activity. Core spliceosomal factors (such as SF3B1 and U2AF1) associated with BUD31 are also required to tolerate oncogenic MYC. Notably, MYC hyperactivation induces an increase in total precursor messenger RNA synthesis, suggesting an increased burden on the core spliceosome to process pre-mRNA. In contrast to normal cells, partial inhibition of the spliceosome in MYC-hyperactivated cells leads to global intron retention, widespread defects in pre-mRNA maturation, and deregulation of many essential cell processes. Notably, genetic or pharmacological inhibition of the spliceosome in vivo impairs survival, tumorigenicity and metastatic proclivity of MYC-dependent breast cancers. Collectively, these data suggest that oncogenic MYC confers a collateral stress on splicing, and that components of the spliceosome may be therapeutic entry points for aggressive MYC-driven cancers.

Clinical features associated with copy number variations of the 14q32 imprinted gene cluster.

Uniparental disomy (UPD) for imprinted chromosomes can cause abnormal phenotypes due to absent or overexpression of imprinted genes. UPD(14)pat causes a unique constellation of features including thoracic skeletal anomalies, polyhydramnios, placentomegaly, and limited survival; its hypothesized cause is overexpression of paternally expressed RTL1, due to absent regulatory effects of maternally expressed RTL1as. UPD(14)mat causes a milder condition with hypotonia, growth failure, and precocious puberty; its hypothesized cause is absence of paternally expressed DLK1. To more clearly establish how gains and losses of imprinted genes can cause disease, we report six individuals with copy number variations of the imprinted 14q32 region identified through clinical microarray-based comparative genomic hybridization. Three individuals presented with UPD(14)mat-like phenotypes (Temple syndrome) and had apparently de novo deletions spanning the imprinted region, including DLK1. One of these deletions was shown to be on the paternal chromosome. Two individuals with UPD(14)pat-like phenotypes had 122-154kb deletions on their maternal chromosomes that included RTL1as but not the differentially methylated regions that regulate imprinted gene expression, providing further support for RTL1 overexpression as a cause for the UPD(14)pat phenotype. The sixth individual is tetrasomic for a 1.7Mb segment, including the imprinted region, and presents with intellectual disability and seizures but lacks significant phenotypic overlap with either UPD(14) syndrome. Therefore, the 14q32 imprinted region is dosage sensitive, with deletions of different critical regions causing UPD(14)mat- and UPD(14)pat-like phenotypes, while copy gains are likely insufficient to recapitulate these phenotypes.

Structural variation and missense mutation in SBDS associated with Shwachman-Diamond syndrome.

BACKGROUND: Shwachman-Diamond syndrome (SDS) is an autosomal recessive ribosomopathy caused mainly by compound heterozygous mutations in SBDS. Structural variation (SV) involving the SBDS locus has been rarely reported in association with the disease. We aimed to determine whether an SV contributed to the pathogenesis of a case lacking biallelic SBDS point mutations. CASE PRESENTATION: Whole exome sequencing was performed in a patient with SDS lacking biallelic SBDS point mutations. Array comparative genomic hybridization and Southern blotting were used to seek SVs across the SBDS locus. Locus-specific polymerase chain reaction (PCR) encompassing flanking intronic sequence was also performed to investigate mutation within the locus. RNA expression and Western blotting were performed to analyze allele and protein expression. We found the child harbored a single missense mutation in SBDS (c.98A > C; p.K33T), inherited from the mother, and an SV in the SBDS locus, inherited from the father. The missense allele and SV segregated in accordance with Mendelian expectations for autosomal recessive SDS. Complementary DNA and western blotting analysis and locus specific PCR support the contention that the SV perturbed SBDS protein expression in the father and child. CONCLUSION: Our findings implicate genomic rearrangements in the pathogenesis of some cases of SDS and support patients lacking biallelic SBDS point mutations be tested for SV within the SBDS locus.

Recurrent HERV-H-mediated 3q13.2-q13.31 deletions cause a syndrome of hypotonia and motor, language, and cognitive delays.

We describe the molecular and clinical characterization of nine individuals with recurrent, 3.4-Mb, de novo deletions of 3q13.2-q13.31 detected by chromosomal microarray analysis. All individuals have hypotonia and language and motor delays; they variably express mild to moderate cognitive delays (8/9), abnormal behavior (7/9), and autism spectrum disorders (3/9). Common facial features include downslanting palpebral fissures with epicanthal folds, a slightly bulbous nose, and relative macrocephaly. Twenty-eight genes map to the deleted region, including four strong candidate genes, DRD3, ZBTB20, GAP43, and BOC, with important roles in neural and/or muscular development. Analysis of the breakpoint regions based on array data revealed directly oriented human endogenous retrovirus (HERV-H) elements of ~5 kb in size and of >95% DNA sequence identity flanking the deletion. Subsequent DNA sequencing revealed different deletion breakpoints and suggested nonallelic homologous recombination (NAHR) between HERV-H elements as a mechanism of deletion formation, analogous to HERV-I-flanked and NAHR-mediated AZFa deletions. We propose that similar HERV elements may also mediate other recurrent deletion and duplication events on a genome-wide scale. Observation of rare recurrent chromosomal events such as these deletions helps to further the understanding of mechanisms behind naturally occurring variation in the human genome and its contribution to genetic disease.

Investigation of NRXN1 deletions: clinical and molecular characterization.

Deletions at 2p16.3 involving exons of NRXN1 are associated with susceptibility for autism and schizophrenia, and similar deletions have been identified in individuals with developmental delay and dysmorphic features. We have identified 34 probands with exonic NRXN1 deletions following referral for clinical microarray-based comparative genomic hybridization. To more firmly establish the full phenotypic spectrum associated with exonic NRXN1 deletions, we report the clinical features of 27 individuals with NRXN1 deletions, who represent 23 of these 34 families. The frequency of exonic NRXN1 deletions among our postnatally diagnosed patients (0.11%) is significantly higher than the frequency among reported controls (0.02%; P = 6.08 × 10(-7) ), supporting a role for these deletions in the development of abnormal phenotypes. Generally, most individuals with NRXN1 exonic deletions have developmental delay (particularly speech), abnormal behaviors, and mild dysmorphic features. In our cohort, autism spectrum disorders were diagnosed in 43% (10/23), and 16% (4/25) had epilepsy. The presence of NRXN1 deletions in normal parents and siblings suggests reduced penetrance and/or variable expressivity, which may be influenced by genetic, environmental, and/or stochastic factors. The pathogenicity of these deletions may also be affected by the location of the deletion within the gene. Counseling should appropriately represent this spectrum of possibilities when discussing recurrence risks or expectations for a child found to have a deletion in NRXN1.

Exonic deletions in AUTS2 cause a syndromic form of intellectual disability and suggest a critical role for the C terminus.

Genomic rearrangements involving AUTS2 (7q11.22) are associated with autism and intellectual disability (ID), although evidence for causality is limited. By combining the results of diagnostic testing of 49,684 individuals, we identified 24 microdeletions that affect at least one exon of AUTS2, as well as one translocation and one inversion each with a breakpoint within the AUTS2 locus. Comparison of 17 well-characterized individuals enabled identification of a variable syndromic phenotype including ID, autism, short stature, microcephaly, cerebral palsy, and facial dysmorphisms. The dysmorphic features were more pronounced in persons with 3'AUTS2 deletions. This part of the gene is shown to encode a C-terminal isoform (with an alternative transcription start site) expressed in the human brain. Consistent with our genetic data, suppression of auts2 in zebrafish embryos caused microcephaly that could be rescued by either the full-length or the C-terminal isoform of AUTS2. Our observations demonstrate a causal role of AUTS2 in neurocognitive disorders, establish a hitherto unappreciated syndromic phenotype at this locus, and show how transcriptional complexity can underpin human pathology. The zebrafish model provides a valuable tool for investigating the etiology of AUTS2 syndrome and facilitating gene-function analysis in the future.

Clinical utility of chromosomal microarray analysis.

OBJECTIVE: To test the hypothesis that chromosomal microarray analysis frequently diagnoses conditions that require specific medical follow-up and that referring physicians respond appropriately to abnormal test results. METHODS: A total of 46,298 postnatal patients were tested by chromosomal microarray analysis for a variety of indications, most commonly intellectual disability/developmental delay, congenital anomalies, dysmorphic features, and neurobehavioral problems. The frequency of detection of abnormalities associated with actionable clinical features was tallied, and the rate of physician response to a subset of abnormal tests results was monitored. RESULTS: A total of 2088 diagnoses were made of more than 100 different disorders that have specific clinical features that warrant follow-up. The detection rate for these conditions using high-resolution whole-genome microarrays was 5.4%, which translates to 35% of all clinically significant abnormal test results identified in our laboratory. In a subset of cases monitored for physician response, appropriate clinical action was taken more than 90% of the time as a direct result of the microarray finding. CONCLUSIONS: The disorders diagnosed by chromosomal microarray analysis frequently have clinical features that need medical attention, and physicians respond to the diagnoses with specific clinical actions, thus arguing that microarray testing provides clinical utility for a significant number of patients tested.

Referral patterns for microarray testing in prenatal diagnosis.

OBJECTIVE: To understand the prenatal referral patterns from the United States, Canada, and Israel for two whole-genome microarray platforms, each with a different resolution. METHOD: Physicians selected one of the two array designs to be performed on 1483 prenatal specimens for a 1-year period. We retrospectively examined detection rates, indications for study, and physician array selection. RESULTS: The lower resolution array (55 K) showed an ~32% decrease in the detection of results of unclear clinical significance while retaining the ability to detect all but one significant abnormality identified by the higher resolution array (135 K). A majority of samples were referred for abnormal ultrasound findings. Whereas the United States and Canada utilized the higher resolution array more often for this indication, Israel preferred the 55 K array. Referral patterns for parental anxiety were similar for the United States and Israel, with most cases being tested on the 55 K array. Few cases were referred for advanced maternal age or family history of a genetic condition from either Canada or Israel. CONCLUSION: Referral patterns varied between the countries and between indications for study. Understanding these differences will provide laboratories the critical information needed to develop array designs to meet the medical needs and patient desires for prenatal testing.

Haploinsufficiency of SOX5 at 12p12.1 is associated with developmental delays with prominent language delay, behavior problems, and mild dysmorphic features.

SOX5 encodes a transcription factor involved in the regulation of chondrogenesis and the development of the nervous system. Despite its important developmental roles, SOX5 disruption has yet to be associated with human disease. We report one individual with a reciprocal translocation breakpoint within SOX5, eight individuals with intragenic SOX5 deletions (four are apparently de novo and one inherited from an affected parent), and seven individuals with larger 12p12 deletions encompassing SOX5. Common features in these subjects include prominent speech delay, intellectual disability, behavior abnormalities, and dysmorphic features. The phenotypic impact of the deletions may depend on the location of the deletion and, consequently, which of the three major SOX5 protein isoforms are affected. One intragenic deletion, involving only untranslated exons, was present in a more mildly affected subject, was inherited from a healthy parent and grandparent, and is similar to a deletion found in a control cohort. Therefore, some intragenic SOX5 deletions may have minimal phenotypic effect. Based on the location of the deletions in the subjects compared to the controls, the de novo nature of most of these deletions, and the phenotypic similarities among cases, SOX5 appears to be a dosage-sensitive, developmentally important gene.

Evaluation of chronic lymphocytic leukemia by oligonucleotide-based microarray analysis uncovers novel aberrations not detected by FISH or cytogenetic analysis.

BACKGROUND: Cytogenetic evaluation is a key component of the diagnosis and prognosis of chronic lymphocytic leukemia (CLL). We performed oligonucleotide-based comparative genomic hybridization microarray analysis on 34 samples with CLL and known abnormal karyotypes previously determined by cytogenetics and/or fluorescence in situ hybridization (FISH). RESULTS: Using a custom designed microarray that targets >1800 genes involved in hematologic disease and other malignancies, we identified additional cryptic aberrations and novel findings in 59% of cases. These included gains and losses of genes associated with cell cycle regulation, apoptosis and susceptibility loci on 3p21.31, 5q35.2q35.3, 10q23.31q23.33, 11q22.3, and 22q11.23. CONCLUSIONS: Our results show that microarray analysis will detect known aberrations, including microscopic and cryptic alterations. In addition, novel genomic changes will be uncovered that may become important prognostic predictors or treatment targets for CLL in the future.

Recurrence, submicroscopic complexity, and potential clinical relevance of copy gains detected by array CGH that are shown to be unbalanced insertions by FISH.

Insertions occur when a segment of one chromosome is translocated and inserted into a new region of the same chromosome or a non-homologous chromosome. We report 71 cases with unbalanced insertions identified using array CGH and FISH in 4909 cases referred to our laboratory for array CGH and found to have copy-number abnormalities. Although the majority of insertions were non-recurrent, several recurrent unbalanced insertions were detected, including three der(Y)ins(Y;18)(q?11.2;p11.32p11.32)pat inherited from parents carrying an unbalanced insertion. The clinical significance of these recurrent rearrangements is unclear, although the small size, limited gene content, and inheritance pattern of each suggests that the phenotypic consequences may be benign. Cryptic, submicroscopic duplications were observed at or near the insertion sites in two patients, further confounding the clinical interpretation of these insertions. Using FISH, linear amplification, and array CGH, we identified a 126-kb duplicated region from 19p13.3 inserted into MECP2 at Xq28 in a patient with symptoms of Rett syndrome. Our results demonstrate that although the interpretation of most non-recurrent insertions is unclear without high-resolution insertion site characterization, the potential for an otherwise benign duplication to result in a clinically relevant outcome through the disruption of a gene necessitates the use of FISH to determine whether copy-number gains detected by array CGH represent tandem duplications or unbalanced insertions. Further follow-up testing using techniques such as linear amplification or sequencing should be used to determine gene involvement at the insertion site after FISH has identified the presence of an insertion.

Comparative analysis of copy number detection by whole-genome BAC and oligonucleotide array CGH.

BACKGROUND: Microarray-based comparative genomic hybridization (aCGH) is a powerful diagnostic tool for the detection of DNA copy number gains and losses associated with chromosome abnormalities, many of which are below the resolution of conventional chromosome analysis. It has been presumed that whole-genome oligonucleotide (oligo) arrays identify more clinically significant copy-number abnormalities than whole-genome bacterial artificial chromosome (BAC) arrays, yet this has not been systematically studied in a clinical diagnostic setting. RESULTS: To determine the difference in detection rate between similarly designed BAC and oligo arrays, we developed whole-genome BAC and oligonucleotide microarrays and validated them in a side-by-side comparison of 466 consecutive clinical specimens submitted to our laboratory for aCGH. Of the 466 cases studied, 67 (14.3%) had a copy-number imbalance of potential clinical significance detectable by the whole-genome BAC array, and 73 (15.6%) had a copy-number imbalance of potential clinical significance detectable by the whole-genome oligo array. However, because both platforms identified copy number variants of unclear clinical significance, we designed a systematic method for the interpretation of copy number alterations and tested an additional 3,443 cases by BAC array and 3,096 cases by oligo array. Of those cases tested on the BAC array, 17.6% were found to have a copy-number abnormality of potential clinical significance, whereas the detection rate increased to 22.5% for the cases tested by oligo array. In addition, we validated the oligo array for detection of mosaicism and found that it could routinely detect mosaicism at levels of 30% and greater. CONCLUSIONS: Although BAC arrays have faster turnaround times, the increased detection rate of oligo arrays makes them attractive for clinical cytogenetic testing.