Integrative Proteogenomics for Differential Expression and Splicing Variation in a DM1 Mouse Model.

Dysregulated mRNA splicing is involved in the pathogenesis of many diseases including cancer, neurodegenerative diseases, and muscular dystrophies such as myotonic dystrophy type 1 (DM1). Comprehensive assessment of dysregulated splicing on the transcriptome and proteome level has been methodologically challenging, and thus investigations have often been targeting only few genes. Here, we performed a large-scale coordinated transcriptomic and proteomic analysis to characterize a DM1 mouse model (HSALR) in comparison to wild type. Our integrative proteogenomics approach comprised gene- and splicing-level assessments for mRNAs and proteins. It recapitulated many known instances of aberrant mRNA splicing in DM1 and identified new ones. It enabled the design and targeting of splicing-specific peptides and confirmed the translation of known instances of aberrantly spliced disease-related genes (e.g., Atp2a1, Bin1, Ryr1), complemented by novel findings (Flnc and Ywhae). Comparative analysis of large-scale mRNA and protein expression data showed quantitative agreement of differentially expressed genes and splicing patterns between disease and wild type. We hence propose this work as a suitable blueprint for a robust and scalable integrative proteogenomic strategy geared toward advancing our understanding of splicing-based disorders. With such a strategy, splicing-based biomarker candidates emerge as an attractive and accessible option, as they can be efficiently asserted on the mRNA and protein level in coordinated fashion.

The ZFP-1(AF10)/DOT-1 complex opposes H2B ubiquitination to reduce Pol II transcription.

The inhibition of transcriptional elongation plays an important role in gene regulation in metazoans, including C. elegans. Here, we combine genomic and biochemical approaches to dissect a role of ZFP-1, the C. elegans AF10 homolog, in transcriptional control. We show that ZFP-1 and its interacting partner DOT-1.1 have a global role in negatively modulating the level of polymerase II (Pol II) transcription on essential widely expressed genes. Moreover, the ZFP-1/DOT-1.1 complex contributes to progressive Pol II pausing on essential genes during development and to rapid Pol II pausing during stress response. The slowing down of Pol II transcription by ZFP-1/DOT-1.1 is associated with an increase in H3K79 methylation and a decrease in H2B monoubiquitination, which promotes transcription. We propose a model wherein the recruitment of ZFP-1/DOT-1.1 and deposition of H3K79 methylation at highly expressed genes initiates a negative feedback mechanism for the modulation of their expression.

Nkx2-1 represses a latent gastric differentiation program in lung adenocarcinoma.

Tissue-specific differentiation programs become dysregulated during cancer evolution. The transcription factor Nkx2-1 is a master regulator of pulmonary differentiation that is downregulated in poorly differentiated lung adenocarcinoma. Here we use conditional murine genetics to determine how the identity of lung epithelial cells changes upon loss of their master cell-fate regulator. Nkx2-1 deletion in normal and neoplastic lungs causes not only loss of pulmonary identity but also conversion to a gastric lineage. Nkx2-1 is likely to maintain pulmonary identity by recruiting transcription factors Foxa1 and Foxa2 to lung-specific loci, thus preventing them from binding gastrointestinal targets. Nkx2-1-negative murine lung tumors mimic mucinous human lung adenocarcinomas, which express gastric markers. Loss of the gastrointestinal transcription factor Hnf4α leads to derepression of the embryonal proto-oncogene Hmga2 in Nkx2-1-negative tumors. These observations suggest that loss of both active and latent differentiation programs is required for tumors to reach a primitive, poorly differentiated state.

Differential Tks5 isoform expression contributes to metastatic invasion of lung adenocarcinoma.

Metastasis accounts for the vast majority of cancer-related deaths, yet the molecular mechanisms that drive metastatic spread remain poorly understood. Here we report that Tks5, which has been linked to the formation of proteolytic cellular protrusions known as invadopodia, undergoes an isoform switch during metastatic progression in a genetically engineered mouse model of lung adenocarcinoma. Nonmetastatic primary tumor-derived cells predominantly expressed a short isoform, Tks5short, while metastatic primary tumor- and metastasis-derived cells acquired increased expression of the full-length isoform Tks5long. This elevation of Tks5long to Tks5short ratio correlated with a commensurate increase in invadopodia activity in metastatic cells compared with nonmetastatic cells. Further characterization of these isoforms by knockdown and overexpression experiments demonstrated that Tks5long promoted invadopodia in vitro and increased metastasis in transplant models and an autochthonous model of lung adenocarcinoma. Conversely, Tks5short decreased invadopodia stability and proteolysis, acting as a natural dominant-negative inhibitor to Tks5long. Importantly, high Tks5long and low Tks5short expressions in human lung adenocarcinomas correlated with metastatic disease and predicted worse survival of early stage patients. These data indicate that tipping the Tks5 isoform balance to a high Tks5long to Tks5short ratio promotes invadopodia-mediated invasion and metastasis.

Suppression of lung adenocarcinoma progression by Nkx2-1.

Despite the high prevalence and poor outcome of patients with metastatic lung cancer the mechanisms of tumour progression and metastasis remain largely uncharacterized. Here we modelled human lung adenocarcinoma, which frequently harbours activating point mutations in KRAS and inactivation of the p53 pathway, using conditional alleles in mice. Lentiviral-mediated somatic activation of oncogenic Kras and deletion of p53 in the lung epithelial cells of Kras(LSL-G12D/+);p53(flox/flox) mice initiates lung adenocarcinoma development. Although tumours are initiated synchronously by defined genetic alterations, only a subset becomes malignant, indicating that disease progression requires additional alterations. Identification of the lentiviral integration sites allowed us to distinguish metastatic from non-metastatic tumours and determine the gene expression alterations that distinguish these tumour types. Cross-species analysis identified the NK2-related homeobox transcription factor Nkx2-1 (also called Ttf-1 or Titf1) as a candidate suppressor of malignant progression. In this mouse model, Nkx2-1 negativity is pathognomonic of high-grade poorly differentiated tumours. Gain- and loss-of-function experiments in cells derived from metastatic and non-metastatic tumours demonstrated that Nkx2-1 controls tumour differentiation and limits metastatic potential in vivo. Interrogation of Nkx2-1-regulated genes, analysis of tumours at defined developmental stages, and functional complementation experiments indicate that Nkx2-1 constrains tumours in part by repressing the embryonically restricted chromatin regulator Hmga2. Whereas focal amplification of NKX2-1 in a fraction of human lung adenocarcinomas has focused attention on its oncogenic function, our data specifically link Nkx2-1 downregulation to loss of differentiation, enhanced tumour seeding ability and increased metastatic proclivity. Thus, the oncogenic and suppressive functions of Nkx2-1 in the same tumour type substantiate its role as a dual function lineage factor.

Comparative oncogenomic analysis of copy number alterations in human and zebrafish tumors enables cancer driver discovery.

The identification of cancer drivers is a major goal of current cancer research. Finding driver genes within large chromosomal events is especially challenging because such alterations encompass many genes. Previously, we demonstrated that zebrafish malignant peripheral nerve sheath tumors (MPNSTs) are highly aneuploid, much like human tumors. In this study, we examined 147 zebrafish MPNSTs by massively parallel sequencing and identified both large and focal copy number alterations (CNAs). Given the low degree of conserved synteny between fish and mammals, we reasoned that comparative analyses of CNAs from fish versus human MPNSTs would enable elimination of a large proportion of passenger mutations, especially on large CNAs. We established a list of orthologous genes between human and zebrafish, which includes approximately two-thirds of human protein-coding genes. For the subset of these genes found in human MPNST CNAs, only one quarter of their orthologues were co-gained or co-lost in zebrafish, dramatically narrowing the list of candidate cancer drivers for both focal and large CNAs. We conclude that zebrafish-human comparative analysis represents a powerful, and broadly applicable, tool to enrich for evolutionarily conserved cancer drivers.

Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1.

The proto-oncogene KRAS is mutated in a wide array of human cancers, most of which are aggressive and respond poorly to standard therapies. Although the identification of specific oncogenes has led to the development of clinically effective, molecularly targeted therapies in some cases, KRAS has remained refractory to this approach. A complementary strategy for targeting KRAS is to identify gene products that, when inhibited, result in cell death only in the presence of an oncogenic allele. Here we have used systematic RNA interference to detect synthetic lethal partners of oncogenic KRAS and found that the non-canonical IkappaB kinase TBK1 was selectively essential in cells that contain mutant KRAS. Suppression of TBK1 induced apoptosis specifically in human cancer cell lines that depend on oncogenic KRAS expression. In these cells, TBK1 activated NF-kappaB anti-apoptotic signals involving c-Rel and BCL-XL (also known as BCL2L1) that were essential for survival, providing mechanistic insights into this synthetic lethal interaction. These observations indicate that TBK1 and NF-kappaB signalling are essential in KRAS mutant tumours, and establish a general approach for the rational identification of co-dependent pathways in cancer.

The matrisome: in silico definition and in vivo characterization by proteomics of normal and tumor extracellular matrices.

The extracellular matrix (ECM) is a complex meshwork of cross-linked proteins providing both biophysical and biochemical cues that are important regulators of cell proliferation, survival, differentiation, and migration. We present here a proteomic strategy developed to characterize the in vivo ECM composition of normal tissues and tumors using enrichment of protein extracts for ECM components and subsequent analysis by mass spectrometry. In parallel, we have developed a bioinformatic approach to predict the in silico "matrisome" defined as the ensemble of ECM proteins and associated factors. We report the characterization of the extracellular matrices of murine lung and colon, each comprising more than 100 ECM proteins and each presenting a characteristic signature. Moreover, using human tumor xenografts in mice, we show that both tumor cells and stromal cells contribute to the production of the tumor matrix and that tumors of differing metastatic potential differ in both the tumor- and the stroma-derived ECM components. The strategy we describe and illustrate here can be broadly applied and, to facilitate application of these methods by others, we provide resources including laboratory protocols, inventories of ECM domains and proteins, and instructions for bioinformatically deriving the human and mouse matrisome.

A conserved PHD finger protein and endogenous RNAi modulate insulin signaling in Caenorhabditis elegans.

Insulin signaling has a profound effect on longevity and the oxidative stress resistance of animals. Inhibition of insulin signaling results in the activation of DAF-16/FOXO and SKN-1/Nrf transcription factors and increased animal fitness. By studying the biological functions of the endogenous RNA interference factor RDE-4 and conserved PHD zinc finger protein ZFP-1 (AF10), which regulate overlapping sets of genes in Caenorhabditis elegans, we identified an important role for these factors in the negative modulation of transcription of the insulin/PI3 signaling-dependent kinase PDK-1. Consistently, increased expression of pdk-1 in zfp-1 and rde-4 mutants contributed to their reduced lifespan and sensitivity to oxidative stress and pathogens due to the reduction in the expression of DAF-16 and SKN-1 targets. We found that the function of ZFP-1 in modulating pdk-1 transcription was important for the extended lifespan of the age-1(hx546) reduction-of-function PI3 kinase mutant, since the lifespan of the age-1; zfp-1 double mutant strain was significantly shorter compared to age-1(hx546). We further demonstrate that overexpression of ZFP-1 caused an increased resistance to oxidative stress in a DAF-16-dependent manner. Our findings suggest that epigenetic regulation of key upstream signaling components in signal transduction pathways through chromatin and RNAi may have a large impact on the outcome of signaling and expression of numerous downstream genes.

Benchmarking tools for transcription factor prioritization.

Spatiotemporal regulation of gene expression is controlled by transcription factor (TF) binding to regulatory elements, resulting in a plethora of cell types and cell states from the same genetic information. Due to the importance of regulatory elements, various sequencing methods have been developed to localise them in genomes, for example using ChIP-seq profiling of the histone mark H3K27ac that marks active regulatory regions. Moreover, multiple tools have been developed to predict TF binding to these regulatory elements based on DNA sequence. As altered gene expression is a hallmark of disease phenotypes, identifying TFs driving such gene expression programs is critical for the identification of novel drug targets. In this study, we curated 84 chromatin profiling experiments (H3K27ac ChIP-seq) where TFs were perturbed through e.g., genetic knockout or overexpression. We ran nine published tools to prioritize TFs using these real-world datasets and evaluated the performance of the methods in identifying the perturbed TFs. This allowed the nomination of three frontrunner tools, namely RcisTarget, MEIRLOP and monaLisa. Our analyses revealed opportunities and commonalities of tools that will help to guide further improvements and developments in the field.

In vivo effects on intron retention and exon skipping by the U2AF large subunit and SF1/BBP in the nematode Caenorhabditis elegans.

The in vivo analysis of the roles of splicing factors in regulating alternative splicing in animals remains a challenge. Using a microarray-based screen, we identified a Caenorhabditis elegans gene, tos-1, that exhibited three of the four major types of alternative splicing: intron retention, exon skipping, and, in the presence of U2AF large subunit mutations, the use of alternative 3' splice sites. Mutations in the splicing factors U2AF large subunit and SF1/BBP altered the splicing of tos-1. 3' splice sites of the retained intron or before the skipped exon regulate the splicing pattern of tos-1. Our study provides in vivo evidence that intron retention and exon skipping can be regulated largely by the identities of 3' splice sites.

Highly aneuploid zebrafish malignant peripheral nerve sheath tumors have genetic alterations similar to human cancers.

Aneuploidy is a hallmark of human cancers, but most mouse cancer models lack the extensive aneuploidy seen in many human tumors. The zebrafish is becoming an increasingly popular model for studying cancer. Here we report that malignant peripheral nerve sheath tumors (MPNSTs) that arise in zebrafish as a result of mutations in either ribosomal protein (rp) genes or in p53 are highly aneuploid. Karyotyping reveals that these tumors frequently harbor near-triploid numbers of chromosomes, and they vary in chromosome number from cell to cell within a single tumor. Using array comparative genomic hybridization, we found that, as in human cancers, certain fish chromosomes are preferentially overrepresented, whereas others are underrepresented in many MPNSTs. In addition, we obtained evidence for recurrent subchromosomal amplifications and deletions that may contain genes involved in cancer initiation or progression. These focal amplifications encompassed several genes whose amplification is observed in human tumors, including met, cyclinD2, slc45a3, and cdk6. One focal amplification included fgf6a. Increasing fgf signaling via a mutation that overexpresses fgf8 accelerated the onset of MPNSTs in fish bearing a mutation in p53, suggesting that fgf6a itself may be a driver of MPNSTs. Our results suggest that the zebrafish is a useful model in which to study aneuploidy in human cancer and in which to identify candidate genes that may act as drivers in fish and potentially also in human tumors.

RNA interference and retinoblastoma-related genes are required for repression of endogenous siRNA targets in Caenorhabditis elegans.

In Caenorhabditis elegans, a vast number of endogenous short RNAs corresponding to thousands of genes have been discovered recently. This finding suggests that these short interfering RNAs (siRNAs) may contribute to regulation of many developmental and other signaling pathways in addition to silencing viruses and transposons. Here, we present a microarray analysis of gene expression in RNA interference (RNAi)-related mutants rde-4, zfp-1, and alg-1 and the retinoblastoma (Rb) mutant lin-35. We found that a component of Dicer complex RDE-4 and a chromatin-related zinc finger protein ZFP-1, not implicated in endogenous RNAi, regulate overlapping sets of genes. Notably, genes a) up-regulated in the rde-4 and zfp-1 mutants and b) up-regulated in the lin-35(Rb) mutant, but not the down-regulated genes are highly represented in the set of genes with corresponding endogenous siRNAs (endo-siRNAs). Our study suggests that endogenous siRNAs cooperate with chromatin factors, either C. elegans ortholog of acute lymphoblastic leukemia-1 (ALL-1)-fused gene from chromosome 10 (AF10), ZFP-1, or tumor suppressor Rb, to regulate overlapping sets of genes and predicts a large role for RNAi-based chromatin silencing in control of gene expression in C. elegans.

New insights into molecular changes in skeletal muscle aging and disease: Differential alternative splicing and senescence.

Progressive loss of muscle mass and function due to muscle fiber atrophy and loss in the elderly and chronically ill is now defined as sarcopenia. It is a major contributor to loss of independence, disability, need of long-term care as well as overall mortality. Sarcopenia is a heterogenous disease and underlying mechanisms are not completely understood. Here, we newly identified and used Tmem158, alongside Cdkn1a, as relevant senescence and denervation markers (SDMs), associated with muscle fiber atrophy. Subsequent application of laser capture microdissection (LCM) and RNA analyses revealed age- and disease-associated differences in gene expression and alternative splicing patterns in a rodent sarcopenia model. Of note, genes exhibiting such differential alternative splicing (DAS) are mainly involved in the contractile function of the muscle. Many of these splicing events are also found in a mouse model for myotonic dystrophy type 1 (DM1), underscoring the premature aging phenotype of this disease. We propose to add differential alternative splicing to the hallmarks of aging.

Periostin shows increased evolutionary plasticity in its alternatively spliced region.

BACKGROUND: Periostin (POSTN) is a secreted extracellular matrix protein of poorly defined function that has been related to bone and heart development as well as to cancer. In human and mouse, it is known to undergo alternative splicing in its C-terminal region, which is devoid of known protein domains. Differential expression of periostin, sometimes of specific splicing isoforms, is observed in a broad range of human cancers, including breast, pancreatic, and colon cancer. Here, we combine genomic and transcriptomic sequence data from vertebrate organisms to study the evolution of periostin and particularly of its C-terminal region. RESULTS: We found that the C-terminal part of periostin is markedly more variable among vertebrates than the rest of periostin in terms of exon count, length, and splicing pattern, which we interpret as a consequence of neofunctionalization after the split between periostin and its paralog transforming growth factor, beta-induced (TGFBI). We also defined periostin's sequential 13-amino acid repeat units--well conserved in teleost fish, but more obscure in higher vertebrates--whose secondary structure is predicted to be consecutive beta strands. We suggest that these beta strands may mediate binding interactions with other proteins through an extended beta-zipper in a manner similar to the way repeat units in bacterial cell wall proteins have been reported to bind human fibronectin. CONCLUSIONS: Our results, obtained with the help of the increasingly large collection of complete vertebrate genomes, document the evolutionary plasticity of periostin's C-terminal region, and for the first time suggest a basis for its functional role.

Gene-signature-derived IC50s/EC50s reflect the potency of causative upstream targets and downstream phenotypes.

Multiplexed gene-signature-based phenotypic assays are increasingly used for the identification and profiling of small molecule-tool compounds and drugs. Here we introduce a method (provided as R-package) for the quantification of the dose-response potency of a gene-signature as EC50 and IC50 values. Two signaling pathways were used as models to validate our methods: beta-adrenergic agonistic activity on cAMP generation (dedicated dataset generated for this study) and EGFR inhibitory effect on cancer cell viability. In both cases, potencies derived from multi-gene expression data were highly correlated with orthogonal potencies derived from cAMP and cell growth readouts, and superior to potencies derived from single individual genes. Based on our results we propose gene-signature potencies as a novel valid alternative for the quantitative prioritization, optimization and development of novel drugs.

Gene expression profiles predict early relapse in ovarian cancer after platinum-paclitaxel chemotherapy.

PURPOSE: Women with advanced epithelial ovarian cancer are routinely treated with platinum-paclitaxel chemotherapy following cytoreductive surgery, yet only approximately 20% achieve long-term disease-free survival. We hypothesized that differences in gene expression before treatment could distinguish patients with short versus long time to recurrence after administration of platinum-paclitaxel combination chemotherapy. EXPERIMENTAL DESIGN: To test this hypothesis, gene expression profiling of 79 primary surgically resected tumors from women with advanced-stage, high-grade epithelial ovarian cancer was done using cDNA microarrays containing 30,721 genes. Supervised learning algorithms were applied in an effort to develop a binary classifier that could discriminate women at risk for early (< or =21 months) versus late (>21 months) relapse after initial chemotherapy. RESULTS: A 14-gene predictive model was developed using a set of training samples (n = 51) and subsequently tested using an independent set of test samples (n = 28). This model correctly predicted the outcome of 24 of the 28 test samples (86% accuracy) with 95% positive predictive value for early relapse. CONCLUSIONS: Predictive markers for early recurrence can be identified for platinum-paclitaxel combination chemotherapy in primary ovarian carcinoma. The proposed 14-gene model requires further validation.

Comparison of the predictive accuracy of DNA array-based multigene classifiers across cDNA arrays and Affymetrix GeneChips.

We examined how well differentially expressed genes and multigene outcome classifiers retain their class-discriminating values when tested on data generated by different transcriptional profiling platforms. RNA from 33 stage I-III breast cancers was hybridized to both Affymetrix GeneChip and Millennium Pharmaceuticals cDNA arrays. Only 30% of all corresponding gene expression measurements on the two platforms had Pearson correlation coefficient r >or= 0.7 when UniGene was used to match probes. There was substantial variation in correlation between different Affymetrix probe sets matched to the same cDNA probe. When cDNA and Affymetrix probes were matched by basic local alignment tool (BLAST) sequence identity, the correlation increased substantially. We identified 182 genes in the Affymetrix and 45 in the cDNA data (including 17 common genes) that accurately separated 91% of cases in supervised hierarchical clustering in each data set. Cross-platform testing of these informative genes resulted in lower clustering accuracy of 45 and 79%, respectively. Several sets of accurate five-gene classifiers were developed on each platform using linear discriminant analysis. The best 100 classifiers showed average misclassification error rate of 2% on the original data that rose to 19.5% when tested on data from the other platform. Random five-gene classifiers showed misclassification error rate of 33%. We conclude that multigene predictors optimized for one platform lose accuracy when applied to data from another platform due to missing genes and sequence differences in probes that result in differing measurements for the same gene.

Global effects of the CSR-1 RNA interference pathway on the transcriptional landscape.

Argonaute proteins and their small RNA cofactors short interfering RNAs are known to inhibit gene expression at the transcriptional and post-transcriptional levels. In Caenorhabditis elegans, the Argonaute CSR-1 binds thousands of endogenous siRNAs (endo-siRNAs) that are antisense to germline transcripts. However, its role in gene expression regulation remains controversial. Here we used genome-wide profiling of nascent RNA transcripts and found that the CSR-1 RNA interference pathway promoted sense-oriented RNA polymerase II transcription. Moreover, a loss of CSR-1 function resulted in global increase in antisense transcription and ectopic transcription of silent chromatin domains, which led to reduced chromatin incorporation of centromere-specific histone H3. On the basis of these findings, we propose that the CSR-1 pathway helps maintain the directionality of active transcription, thereby propagating the distinction between transcriptionally active and silent genomic regions.