Alternative splicing and differential gene expression in colon cancer detected by a whole genome exon array. BMC Genomics. 2006;7:325. [PMID: 17192196 PMCID: PMC1769375 DOI: 10.1186/1471-2164-7-325]

A reliable and quick method for screening alternative splicing variants for low-abundance genes

Alternative splicing (AS) is a universal phenomenon in eukaryotes, and it is still challenging to identify AS events. Several methods have been developed to identify AS events, such as expressed sequence tags (EST), microarrays and RNA-seq. However, EST has limitations in identifying low-abundance genes, while microarray and RNA-seq are high-throughput technologies, and PCR-based technology is needed for validation. To overcome the limitations of EST and shortcomings of high-throughput technologies, we established a method to identify AS events, especially for low-abundance genes, by reverse transcription (RT) PCR with gene-specific primers (GSPs) followed by nested PCR. This process includes two major steps: 1) the use of GSPs to amplify as long as the specific gene segment and 2) multiple rounds of nested PCR to screen the AS and confirm the unknown splicing variants. With this method, we successfully identified three new splicing variants, namely, GenBank Accession No. HM623886 for the bdnf gene (GenBank GeneID: 12064), GenBank Accession No. JF417977 for the trkc gene (GenBank GeneID: 18213) and GenBank Accession No. HM623888 for the glb-18 gene (GenBank GeneID: 172485). In addition to its reliability and simplicity, the method is also cost-effective and labor-intensive. In conclusion, we developed an RT-nested PCR method using gene-specific primers to efficiently identify known and novel AS variants. This approach overcomes the limitations of existing methods for detecting rare transcripts. By enabling the discovery of new isoforms, especially for low-abundance genes, this technique can aid research into aberrant splicing in disease. Future studies can apply this method to uncover AS variants involved in cancer, neurodegeneration, and other splicing-related disorders.

Alternative splicing in prostate cancer progression and therapeutic resistance

Prostate cancer (CaP) remains the second leading cause of cancer deaths in western men. CaP mortality results from diverse molecular mechanisms that mediate resistance to the standard of care treatments for metastatic disease. Recently, alternative splicing has been recognized as a hallmark of CaP aggressiveness. Alternative splicing events cause treatment resistance and aggressive CaP behavior and are determinants of the emergence of the two major types of late-stage treatment-resistant CaP, namely castration-resistant CaP (CRPC) and neuroendocrine CaP (NEPC). Here, we review recent multi-omics data that are uncovering the complicated landscape of alternative splicing events during CaP progression and the impact that different gene transcript isoforms can have on CaP cell biology and behavior. We discuss renewed insights in the molecular machinery by which alternative splicing occurs and contributes to the failure of systemic CaP therapies. The potential for alternative splicing events to serve as diagnostic markers and/or therapeutic targets is explored. We conclude by considering current challenges and promises associated with splicing-modulating therapies, and their potential for clinical translation into CaP patient care.

Alternatively Spliced Isoforms of MUC4 and ADAM12 as Biomarkers for Colorectal Cancer Metastasis

There is a pertinent need to develop prognostic biomarkers for practicing predictive, preventive and personalized medicine (PPPM) in colorectal cancer metastasis. The analysis of isoform expression data governed by alternative splicing provides a high-resolution picture of mRNAs in a defined condition. This information would not be available by studying gene expression changes alone. Hence, we utilized our prior data from an exon microarray and found ADAM12 and MUC4 to be strong biomarker candidates based on their alternative splicing scores and pattern. In this study, we characterized their isoform expression in a cell line model of metastatic colorectal cancer (SW480 & SW620). These two genes were found to be good prognostic indicators in two cohorts from The Cancer Genome Atlas database. We studied their exon structure using sequence information in the NCBI and ENSEMBL genome databases to amplify and validate six isoforms each for the ADAM12 and MUC4 genes. The differential expression of these isoforms was observed between normal, primary and metastatic colorectal cancer cell lines. RNA-Seq analysis further proved the differential expression of the gene isoforms. The isoforms of MUC4 and ADAM12 were found to change expression levels in response to 5-Fluorouracil (5-FU) treatment in a dose-, time- and cell line-dependent manner. Furthermore, we successfully detected the protein isoforms of ADAM12 and MUC4 in cell lysates, reflecting the differential expression at the protein level. The change in the mRNA and protein expression of MUC4 and ADAM12 in primary and metastatic cells and in response to 5-FU qualifies them to be studied as potential biomarkers. This comprehensive study underscores the importance of studying alternatively spliced isoforms and their potential use as prognostic and/or predictive biomarkers in the PPPM approach towards cancer.

Knowledge mapping of alternative splicing of cancer from 2012 to 2021: A bibliometric analysis

Background

As a processing method of RNA precursors, alternative splicing (AS) is critical to normal cellular activities. Aberrant AS events are associated with cancer development and can be promising targets to treat cancer. However, no detailed and unbiased study describes the current state of AS of cancer research. We aim to measure and recognize the current state and trends of AS cancer research in this study.

Methods

The Web of Science Core Collection was used to acquire the articles. Utilizing three bibliometric tools (CiteSpace, VOSviewer, R-bibliometrix), we were able to measure and recognize the influence and collaboration data of individual articles, journals, and co-citations. Analysis of co-occurrence and burst information helped us identify the trending research areas related to AS of cancer.

Results

From 2012 to 2021, the total number of papers on AS of cancer published in 766 academic journals was 3,507, authored by 20,406 researchers in 405 institutions from 80 countries/regions. Research involving AS of cancer genes was primarily conducted in the United States and China; simultaneously, the Chinese Academy of Sciences, Fudan University, and National Cancer Institute were the institutions with strong research capabilities. Scorilas Andreas is the scholar with the most publications, while the most co-citations were generated by Wang, Eric T. Plos One published the most papers on AS of cancer, while J Biol Chem was the most co-cited academic journal in this field. The results of keyword co-occurrence analysis can be divided into three types: molecular (P53, CD44, androgen receptor, srsf3, esrp1), pathological process (apoptosis, EMT, metastasis, angiogenesis, proliferation), and disease (breast cancer, colorectal cancer, prostate cancer, hepatocellular carcinoma, gastric cancer).

Conclusion

Research on AS of cancer has been increasing in intensity over the past decade. Current AS of cancer studies focused on the hallmarks of AS in cancer and AS signatures including diagnostic and therapeutic targets. Among them, the current trends are splicing factors regulating epithelial-mesenchymal transition and other hallmarks, aberrant splicing events in tumors, and further mechanisms. These might give researchers interested in this field a forward-looking perspective and inform further research.

Aerocyte specification and lung adaptation to breathing is dependent on alternative splicing changes

Adaptation to breathing is a critical step in lung function and it is crucial for organismal survival. Alveoli are the lung gas exchange units and their development, from late embryonic to early postnatal stages, requires feedbacks between multiple cell types. However, how the crosstalk between the alveolar cell types is modulated to anticipate lung adaptation to breathing is still unclear. Here, we uncovered a synchronous alternative splicing switch in multiple genes in the developing mouse lungs at the transition to birth, and we identified hnRNP A1, Cpeb4, and Elavl2/HuB as putative splicing regulators of this transition. Notably, we found that Vegfa switches from the Vegfa 164 isoform to the longer Vegfa 188 isoform exclusively in lung alveolar epithelial AT1 cells. Functional analysis revealed that VEGFA 188 (and not VEGFA 164) drives the specification of Car4-positive aerocytes, a subtype of alveolar endothelial cells specialized in gas exchanges. Our results reveal that the cell type-specific regulation of Vegfa alternative splicing just before birth modulates the epithelial-endothelial crosstalk in the developing alveoli to promote lung adaptation to breathing.

DJExpress: An Integrated Application for Differential Splicing Analysis and Visualization

RNA-seq analysis of alternative pre-mRNA splicing has facilitated an unprecedented understanding of transcriptome complexity in health and disease. However, despite the availability of countless bioinformatic pipelines for transcriptome-wide splicing analysis, the use of these tools is often limited to expert bioinformaticians. The need for high computational power, combined with computational outputs that are complicated to visualize and interpret present obstacles to the broader research community. Here we introduce DJExpress, an R package for differential expression analysis of transcriptomic features and expression-trait associations. To determine gene-level differential junction usage as well as associations between junction expression and molecular/clinical features, DJExpress uses raw splice junction counts as input data. Importantly, DJExpress runs on an average laptop computer and provides a set of interactive and intuitive visualization formats. In contrast to most existing pipelines, DJExpress can handle both annotated and de novo identified splice junctions, thereby allowing the quantification of novel splice events. Moreover, DJExpress offers a web-compatible graphical interface allowing the analysis of user-provided data as well as the visualization of splice events within our custom database of differential junction expression in cancer (DJEC DB). DJEC DB includes not only healthy and tumor tissue junction expression data from TCGA and GTEx repositories but also cancer cell line data from the DepMap project. The integration of DepMap functional genomics data sets allows association of junction expression with molecular features such as gene dependencies and drug response profiles. This facilitates identification of cancer cell models for specific splicing alterations that can then be used for functional characterization in the lab. Thus, DJExpress represents a powerful and user-friendly tool for exploration of alternative splicing alterations in RNA-seq data, including multi-level data integration of alternative splicing signatures in healthy tissue, tumors and cancer cell lines.

An ensemble machine learning model based on multiple filtering and supervised attribute clustering algorithm for classifying cancer samples

BACKGROUND

Machine learning is one kind of machine intelligence technique that learns from data and detects inherent patterns from large, complex datasets. Due to this capability, machine learning techniques are widely used in medical applications, especially where large-scale genomic and proteomic data are used. Cancer classification based on bio-molecular profiling data is a very important topic for medical applications since it improves the diagnostic accuracy of cancer and enables a successful culmination of cancer treatments. Hence, machine learning techniques are widely used in cancer detection and prognosis.

METHODS

In this article, a new ensemble machine learning classification model named Multiple Filtering and Supervised Attribute Clustering algorithm based Ensemble Classification model (MFSAC-EC) is proposed which can handle class imbalance problem and high dimensionality of microarray datasets. This model first generates a number of bootstrapped datasets from the original training data where the oversampling procedure is applied to handle the class imbalance problem. The proposed MFSAC method is then applied to each of these bootstrapped datasets to generate sub-datasets, each of which contains a subset of the most relevant/informative attributes of the original dataset. The MFSAC method is a feature selection technique combining multiple filters with a new supervised attribute clustering algorithm. Then for every sub-dataset, a base classifier is constructed separately, and finally, the predictive accuracy of these base classifiers is combined using the majority voting technique forming the MFSAC-based ensemble classifier. Also, a number of most informative attributes are selected as important features based on their frequency of occurrence in these sub-datasets.

RESULTS

To assess the performance of the proposed MFSAC-EC model, it is applied on different high-dimensional microarray gene expression datasets for cancer sample classification. The proposed model is compared with well-known existing models to establish its effectiveness with respect to other models. From the experimental results, it has been found that the generalization performance/testing accuracy of the proposed classifier is significantly better compared to other well-known existing models. Apart from that, it has been also found that the proposed model can identify many important attributes/biomarker genes.

Circular RNA hsa_circ_0006401 promotes proliferation and metastasis in colorectal carcinoma

Dysregulation of circular RNA (circRNA) expression is involved in the progression of cancer. Here, we aimed to study the potential function of hsa_circ_0006401 in colorectal cancer (CRC). CircRNA hsa_circ_0006401 expression levels in CRC and adjacent nontumor tissues were analyzed by real-time quantitative PCR (qRT-PCR) and circRNA in situ hybridization (RNA-ISH). Then, CRC cell proliferation was assessed by cell counting. Wound-healing and transwell assays were utilized to detect the effect of hsa_circ_0006401 on CRC migration. A circRNA-ORF construct was created, and a specific antibody against the splice junction of hsa_circ_0006401 was prepared. Finally, the proteins directly binding to hsa_circ_0006401 peptides were identified by immunoprecipitation combined with mass spectrometry. In our study, we found hsa_circ_0006401 was closely related to CRC metastasis and exhibited upregulated expression in metastatic CRC tissue samples. Proliferation and migration were inhibited in vitro when hsa_circ_0006401 expression was silenced. Downregulation of hsa_circ_0006401 expression decreased CRC proliferation and liver metastasis in vivo. A 198-aa peptide was encoded by sequences of the splice junction absent from col6a3. Hsa_circ_0006401 promoted CRC proliferation and migration by encoding the hsa_circ_0006401 peptide. Hsa_circ_0006401 peptides decreased the mRNA and protein level of the host gene col6a3 by promoting col6a3 mRNA stabilation. In conclusion, our study revealed that circRNAs generated from col6a3 that contain an open-reading frame (ORF) encode a novel 198-aa functional peptide and hsa_circ_0006401 peptides promote stability of the host gene col6a3 mRNA to promote CRC proliferation and metastasis.

Alternative splicing in normal and pathological human placentas is correlated to genetic variants

Two major obstetric diseases, preeclampsia (PE), a pregnancy-induced endothelial dysfunction leading to hypertension and proteinuria, and intra-uterine growth-restriction (IUGR), a failure of the fetus to acquire its normal growth, are generally triggered by placental dysfunction. Many studies have evaluated gene expression deregulations in these diseases, but none has tackled systematically the role of alternative splicing. In the present study, we show that alternative splicing is an essential feature of placental diseases, affecting 1060 and 1409 genes in PE vs controls and IUGR vs controls, respectively, many of those involved in placental function. While in IUGR placentas, alternative splicing affects genes specifically related to pregnancy, in preeclamptic placentas, it impacts a mix of genes related to pregnancy and brain diseases. Also, alternative splicing variations can be detected at the individual level as sharp splicing differences between different placentas. We correlate these variations with genetic variants to define splicing Quantitative Trait Loci (sQTL) in the subset of the 48 genes the most strongly alternatively spliced in placental diseases. We show that alternative splicing is at least partly piloted by genetic variants located either in cis (52 QTL identified) or in trans (52 QTL identified). In particular, we found four chromosomal regions that impact the splicing of genes in the placenta. The present work provides a new vision of placental gene expression regulation that warrants further studies.

Identification of critical molecular pathways involved in exosome-mediated improvement of cardiac function in a mouse model of muscular dystrophy

Duchenne muscular dystrophy (DMD) is a progressive disease characterized by skeletal muscle atrophy, respiratory failure, and cardiomyopathy. Our previous studies have shown that transplantation with allogeneic myogenic progenitor-derived exosomes (MPC-Exo) can improve cardiac function in X-linked muscular dystrophy (Mdx) mice. In the present study we explored the molecular mechanisms underlying this beneficial effect. We quantified gene expression in the hearts of two strains of Mdx mice (D2.B10-DmdMdx/J and Utrntm1Ked-DmdMdx/J). Two days after MPC-Exo or control treatment, we performed unbiased next-generation RNA-sequencing to identify differentially expressed genes (DEGs) in treated Mdx hearts. Venn diagrams show a set of 780 genes that were ≥2-fold upregulated, and a set of 878 genes that were ≥2-fold downregulated, in both Mdx strains following MPC-Exo treatment as compared with control. Gene ontology (GO) and protein-protein interaction (PPI) network analysis showed that these DEGs were involved in a variety of physiological processes and pathways with a complex connection. qRT-PCR was performed to verify the upregulated ATP2B4 and Bcl-2 expression, and downregulated IL-6, MAPK8 and Wnt5a expression in MPC-Exo-treated Mdx hearts. Western blot analysis verified the increased level of Bcl-2 and decreased level of IL-6 protein in MPC-Exo-treated Mdx hearts compared with control treatment, suggesting that anti-apoptotic and anti-inflammatory effects might be responsible for heart function improvement by MPC-Exo. Based on these findings, we believed that these DEGs might be therapeutic targets that can be explored to develop new strategies for treating DMD.

Construction and prognostic analysis of miRNA-mRNA regulatory network in liver metastasis from colorectal cancer

Background

Colorectal cancer (CRC) is one of the most common cancers in the world, and liver metastasis is the leading cause of colorectal cancer-related deaths. However, the mechanism of liver metastasis in CRC has not been clearly elucidated.

Methods

Three datasets from the Gene Expression Omnibus (GEO) database were analyzed to obtain differentially expressed genes (DEGs), which were subjected to functional enrichment analysis and protein-protein interaction analysis. Subsequently, mRNA-miRNA network was constructed, and the associated DEGs and DEMs were performed for prognostic analysis. Finally, we did infiltration analysis of growth arrest specific 1 (GAS1)-associated immune cells.

Results

We obtained 325 DEGs and 9 differentially expressed miRNAs (DEMs) between primary CRC and liver metastases. Enrichment analysis and protein-protein interactions (PPI) further revealed the involvement of DEGs in the formation of the inflammatory microenvironment and epithelial-mesenchymal transition (EMT) during the liver metastases process in CRC. Survival analysis demonstrated that low-expressed GAS1 as well as low-expressed hsa-miR-33b-5p was a favorable prognostic indicator of overall survival. Further exploration of GAS1 revealed that its expression was interrelated with the infiltration of immune cells in tumor tissues.

Conclusions

In summary, DEGs, DEMs, and their interactions found in liver metastasis of CRC may provide a basis for further understanding of the mechanism of CRC metastasis.

A Meta-Analysis of Proteomic Blood Markers of Colorectal Cancer

BACKGROUND

Early diagnosis will significantly improve the survival rate of colorectal cancer (CRC); however, the existing methods for CRC screening were either invasive or inefficient. There is an emergency need for novel markers in CRC's early diagnosis. Serum proteomics has gained great potential in discovering novel markers, providing markers that reflect the early stage of cancer and prognosis prediction of CRC. In this paper, the results of proteomics of CRC studies were summarized through a meta-analysis in order to obtain the diagnostic efficiency of novel markers.

METHODS

A systematic search on bibliographic databases was performed to collect the studies that explore blood-based markers for CRC applying proteomics. The detection and validation methods, as well as the specificity and sensitivity of the biomarkers in these studies, were evaluated. Newcastle- Ottawa Scale (NOS) case-control studies version was used for quality assessment of included studies.

RESULTS

Thirty-four studies were selected from 751 studies, in which markers detected by proteomics were summarized. In total, fifty-nine proteins were classified according to their biological function. The sensitivity, specificity, or AUC varied among these markers. Among them, Mammalian STE20-like protein kinase 1/ Serine threonine kinase 4 (MST1/STK4), S100 calcium-binding protein A9 (S100A9), and Tissue inhibitor of metalloproteinases 1 (TIMP1) were suitable for effect sizes merging, and their diagnostic efficiencies were recalculated after merging. MST1/STK4 obtained a sensitivity of 68% and a specificity of 78%. S100A9 achieved a sensitivity of 72%, a specificity of 83%, and an AUC of 0.88. TIMP1 obtained a sensitivity of 42%, a specificity of 88%, and an AUC of 0.71.

CONCLUSION

MST1/STK4, S100A9, and TIMP1 showed excellent performance for CRC detection. Several other markers also presented optimized diagnostic efficacy for CRC early detection, but further verification is still needed before they are suitable for clinical use. The discovering of more efficient markers will benefit CRC treatment.

Alternative Splicing: Expanding the Landscape of Cancer Biomarkers and Therapeutics

Alternative splicing (AS) is a critical post-transcriptional regulatory mechanism used by more than 95% of transcribed human genes and responsible for structural transcript variation and proteome diversity. In the past decade, genome-wide transcriptome sequencing has revealed that AS is tightly regulated in a tissue- and developmental stage-specific manner, and also frequently dysregulated in multiple human cancer types. It is currently recognized that splicing defects, including genetic alterations in the spliced gene, altered expression of both core components or regulators of the precursor messenger RNA (pre-mRNA) splicing machinery, or both, are major drivers of tumorigenesis. Hence, in this review we provide an overview of our current understanding of splicing alterations in cancer, and emphasize the need to further explore the cancer-specific splicing programs in order to obtain new insights in oncology. Furthermore, we also discuss the recent advances in the identification of dysregulated splicing signatures on a genome-wide scale and their potential use as biomarkers. Finally, we highlight the therapeutic opportunities arising from dysregulated splicing and summarize the current approaches to therapeutically target AS in cancer.

Genome-Wide Transcriptional Analysis Reveals Alternative Splicing Event Profiles in Hepatocellular Carcinoma and Their Prognostic Significance

Accumulating evidence indicates an unexpected role of aberrant splicing in hepatocellular carcinoma (HCC) that has been seriously neglected in previous studies. There is a need for a detailed analysis of alternative splicing (AS) and its underlying biological and clinical relevance in HCC. In this study, clinical information and corresponding RNA sequencing data of HCC patients were obtained from The Cancer Genome Atlas. Percent spliced in (PSI) values and transcriptional splicing patterns of genes were determined from the original RNA sequencing data using SpliceSeq. Then, based on the PSI values of AS events in different patients, a series of bioinformatics methods was used to identify differentially expressed AS events (DEAS), determine potential regulatory relationships, and investigate the correlation between DEAS and the patients' clinicopathological features. Finally, 25,934 AS events originating from 8,795 genes were screened with high reliability; 263 of these AS events were identified as DEAS. The parent genes of these DEAS formed an intricate network with roles in the regulation of cancer-related pathway and liver metabolism. In HCC, 36 splicing factors were involved in the dysregulation of part DEAS, 100 DEAS events were correlated with overall survival, and 71 DEAS events were correlated with disease-free survival. Stratifying HCC patients according to DEAS resulted in four clusters with different survival patterns. Significant variations in AS occurred during HCC initiation and maintenance; these are likely to be vital both for biological processes and in prognosis. The HCC-related AS events identified here and the splicing networks constructed will be valuable in deciphering the underlying role of AS in HCC.

The Landscape of Tumor-Specific Antigens in Colorectal Cancer

Over the last few decades, major efforts in cancer research and treatment have intensified. Apart from standard chemotherapy approaches, immunotherapy has gained substantial traction. Personalized immunotherapy has become an important tool for cancer therapy with the discovery of immune checkpoint inhibitors. Traditionally, tumor-associated antigens are used in immunotherapy-based treatments. Nevertheless, these antigens lack specificity and may have increased toxicity. With the advent of next-generation technologies, the identification of new tumor-specific antigens is becoming more important. In colorectal cancer, several tumor-specific antigens were identified and functionally validated. Multiple clinical trials from vaccine-based and adoptive cell therapy utilizing tumor-specific antigens have commenced. Herein, we will summarize the current landscape of tumor-specific antigens particularly in colorectal cancer.

Analysis of relapse-associated alternative mRNA splicing and construction of a prognostic signature predicting relapse in I-III colon cancer

The literature comprehensively analyzed alternative splicing (AS) events in colon cancer is little and corresponding prognostic signature is still a lack. Based on data of TCGA, the relapse-associated ASs were comprehensively analyzed and a signature was further constructed to predict the relapse in I-III colon cancer. In total 1912 ASs of 1384 mRNA were identified as relapse-associated ASs, protein-protein interactions (PPI) and ASs-splicing factors (SF) interactions network were identified. We finally built a robust signature to predict the relapse of I-III colon cancer with a considerable AUC value in both the training group and the test group. The AUC in the entire set at 1, 3 and 5 year was 0.85, 0.83 and 0.836. Our study provided a profile of relapse-associated ASs in I-III colon cancer and built a robust signature to predict the relapse of I-III colon cancer.

The Biological Role of the Collagen Alpha-3 (VI) Chain and Its Cleaved C5 Domain Fragment Endotrophin in Cancer

The collagen alpha-3 (VI) chain encoded by the gene COL6A3 is one of the 3 subunits of collagen VI which is a microfibrillar component of the extracellular matrix and is essential for the stable assembly process of collagen VI. The collagen alpha-3 (VI) chain and the cleaved C5 domain fragment, called endotrophin, are highly expressed in a variety of cancers and play a crucial role in cancer progression. The biological functions of endotrophin in tumors can be driven by adipocytes. Studies have demonstrated that endotrophin can directly affect the malignancy of cancer cells through TGF-β-dependent mechanisms, inducing epithelial–mesenchymal transition and fibrosis of the tumor microenvironment. In addition, endotrophin can also recruit macrophages and endothelial cells through chemotaxis to regulate the tumor microenvironment and ultimately promote tumor inflammation and angiogenesis. Furthermore, COL6A3 and endotrophin serve as novel diagnostic and prognostic biomarkers in cancer and contribute to clinical therapeutic applications in the future. In summary, in this review, we discuss the importance of the collagen alpha-3 (VI) chain and endotrophin in cancer progression, the future clinical applications of endotrophin and the remaining challenges in this field.

A signature predicting relapse based on integrated analysis on relapse-associated alternative mRNA splicing in I-III rectal cancer

Researches focusing on the effects of alternative splicing (AS) on relapse of rectal cancer is little and signature based on the AS is blank. In this study, bioinformatic analysis was performed to identify and analyze the relapse-associated ASs, a signature was also constructed. In conclusion, 829 relapse-associated ASs of 676 mRNA were identified. 603 proteins with 2119 interactions were involved in the PPI (protein-protein interactions) network. 43 relapse-associated ASs and 64 SFs (splicing factors) with 160 interactions were indicated. Finally, we built a robust signature to predict the relapse of I-III rectal cancer with a high AUC (0.98) of ROC at 1 year. Based on the ASs involved in the signature, 4 molecular subgroups that could distinguish the relapse rate in diverse groups were identified. Our research provided an overview of relapse-associated ASs in I-III rectal cancer.

Translational activation of maternally derived mRNAs in oocytes and early embryos and the role of embryonic poly(A) binding protein (EPAB)

Transcription ceases upon stimulation of oocyte maturation and gene expression during oocyte maturation, fertilization, and early cleavage relies on translational activation of maternally derived mRNAs. Two key mechanisms that mediate translation of mRNAs in oocytes have been described in detail: cytoplasmic polyadenylation-dependent and -independent. Both of these mechanisms utilize specific protein complexes that interact with cis-acting sequences located on 3'-untranslated region (3'-UTR), and both involve embryonic poly(A) binding protein (EPAB), the predominant poly(A) binding protein during early development. While mechanistic details of these pathways have primarily been elucidated using the Xenopus model, their roles are conserved in mammals and targeted disruption of key regulators in mouse results in female infertility. Here, we provide a detailed account of the molecular mechanisms involved in translational activation during oocyte and early embryo development, and the role of EPAB in this process.

Aberrant expression of CPSF1 promotes head and neck squamous cell carcinoma via regulating alternative splicing

Alternative mRNA splicing increases protein diversity, and alternative splicing events (ASEs) drive oncogenesis in multiple tumor types. However, the driving alterations that underlie the broad dysregulation of ASEs are incompletely defined. Using head and neck squamous cell carcinoma (HNSCC) as a model, we hypothesized that the genomic alteration of genes associated with the spliceosome may broadly induce ASEs across a broad range of target genes, driving an oncogenic phenotype. We identified 319 spliceosome genes and employed a discovery pipeline to identify 13 candidate spliceosome genes altered in HNSCC using The Cancer Genome Atlas (TCGA) HNSCC data. Phenotypic screens identified amplified and overexpressed CPSF1 as a target gene alteration that was validated in proliferation, colony formation, and apoptosis assays in cell line and xenograft systems as well as in primary HNSCC. We employed knockdown and overexpression assays followed by identification of ASEs regulated by CPSF1 overexpression to identify changes in ASEs, and the expression of these ASEs was validated using RNA from cell line models. Alterations in expression of spliceosome genes, including CPSF1, may contribute to HNSCC by mediating aberrant ASE expression.

MyoMiner: explore gene co-expression in normal and pathological muscle

BACKGROUND

High-throughput transcriptomics measures mRNA levels for thousands of genes in a biological sample. Most gene expression studies aim to identify genes that are differentially expressed between different biological conditions, such as between healthy and diseased states. However, these data can also be used to identify genes that are co-expressed within a biological condition. Gene co-expression is used in a guilt-by-association approach to prioritize candidate genes that could be involved in disease, and to gain insights into the functions of genes, protein relations, and signaling pathways. Most existing gene co-expression databases are generic, amalgamating data for a given organism regardless of tissue-type.

METHODS

To study muscle-specific gene co-expression in both normal and pathological states, publicly available gene expression data were acquired for 2376 mouse and 2228 human striated muscle samples, and separated into 142 categories based on species (human or mouse), tissue origin, age, gender, anatomic part, and experimental condition. Co-expression values were calculated for each category to create the MyoMiner database.

RESULTS

Within each category, users can select a gene of interest, and the MyoMiner web interface will return all correlated genes. For each co-expressed gene pair, adjusted p-value and confidence intervals are provided as measures of expression correlation strength. A standardized expression-level scatterplot is available for every gene pair r-value. MyoMiner has two extra functions: (a) a network interface for creating a 2-shell correlation network, based either on the most highly correlated genes or from a list of genes provided by the user with the option to include linked genes from the database and (b) a comparison tool from which the users can test whether any two correlation coefficients from different conditions are significantly different.

CONCLUSIONS

These co-expression analyses will help investigators to delineate the tissue-, cell-, and pathology-specific elements of muscle protein interactions, cell signaling and gene regulation. Changes in co-expression between pathologic and healthy tissue may suggest new disease mechanisms and help define novel therapeutic targets. Thus, MyoMiner is a powerful muscle-specific database for the discovery of genes that are associated with related functions based on their co-expression. MyoMiner is freely available at https://www.sys-myo.com/myominer.

Whole-exome and RNA sequencing of pulmonary carcinoid reveals chromosomal rearrangements associated with recurrence

INTRODUCTION

The majority of pulmonary carcinoid (PC) tumors can be cured by surgical resection alone, but a significant proportion of patients experience recurrence. As PC is insensitive to conventional chemotherapy, further clarification of the molecular mechanisms of metastasis is needed in order to develop targeted therapeutics.

METHODS

We performed comprehensive whole-exome sequencing (WES) of primary tumors and corresponding normal lung tissues from 14 PC patients (including 4 patients who developed postsurgical distant metastasis) and RNA sequencing of primary tumors from 6 PC patients (including 4 patients who developed postsurgical distant metastasis). Exon array-based gene expression analysis was performed in 25 cases of PC.

RESULTS

We identified a total of 139 alterations in 136 genes. MUC6 and SPTA1 were recurrently mutated at a frequency of 21% (3/14) and 14% (2/14), respectively. Mucin protein family genes including MUC2, MUC4 and MUC6 were mutated in a mutually exclusive manner in 36% (5/14). Pathway analysis of the mutated genes revealed enrichment of genes involved in mitogen-activated protein kinase (MAPK) signaling, regulation of the actin cytoskeleton and focal adhesion, and transforming growth factor (TGF)-β signaling. RNA sequencing revealed a total of 8 novel fusion transcripts including one derived from a chromosomal translocation between the TRIB2 and PRKCE genes. All of the 8 fusion genes were detected in primary PCs that had developed metastasis after surgical resection. We identified 14 genes (DENND1B, GRID1, CLMN, DENND1B, NRP1, SEL1L3, C5orf13, TNFRSF21, TES, STK39, MTHFD2, OPN3, MET, and HIST1H3C) up-regulated in 5 PCs that had relapsed after surgical resection.

CONCLUSIONS

In this study we identified novel somatic mutations and chromosomal rearrangements in PC by examining clinically aggressive cases that had developed postsurgical metastasis. It will be essential to validate the clinical significance of these genetic changes in a larger independent patient cohort.

A functional network of gastric-cancer-associated splicing events controlled by dysregulated splicing factors

Comprehensive understanding of aberrant splicing in gastric cancer is lacking. We RNA-sequenced 19 gastric tumor–normal pairs and identified 118 high-confidence tumor-associated (TA) alternative splicing events (ASEs) based on high-coverage sequencing and stringent filtering, and also identified 8 differentially expressed splicing factors (SFs). The TA ASEs occurred in genes primarily involved in cytoskeletal organization. We constructed a correlative network between TA ASE splicing ratios and SF expression, replicated it in independent gastric cancer data from The Cancer Genome Atlas and experimentally validated it by knockdown of the nodal SFs (PTBP1, ESRP2 and MBNL1). Each SF knockdown drove splicing alterations in several corresponding TA ASEs and led to alterations in cellular migration consistent with the role of TA ASEs in cytoskeletal organization. We have therefore established a robust network of dysregulated splicing associated with tumor invasion in gastric cancer. Our work is a resource for identifying oncogenic splice forms, SFs and splicing-generated tumor antigens as biomarkers and therapeutic targets.

2,3,7,8-Tetrachlorodibenzo-p-dioxin modifies alternative splicing in mouse liver

Alternative splicing is a co-transcriptional mechanism that generates protein diversity by including or excluding exons in different combinations, thereby expanding the diversity of protein isoforms of a single gene. Abnormalities in this process can result in deleterious effects to human health, and several xenobiotics are known to interfere with splicing regulation through multiple mechanisms. These changes could lead to human diseases such as cancer, neurological disorders, autoimmune diseases, and developmental disorders. 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is an environmental contaminant generated as a byproduct of various industrial activities. Exposure to this dioxin has been linked to a wide range of pathologies through the alteration of multiple cellular processes. However, the effects of TCDD exposure on alternative splicing have not yet been studied. Here, we investigated whether a single po. dose of 5 μg/kg or 500 μg/kg TCDD influence hepatic alternative splicing in adult male C57BL/6Kou mouse. We identified several genes whose alternative splicing of precursor messenger RNAs was modified following TCDD exposure. In particular, we demonstrated that alternative splicing of Cyp1a1, Ahrr, and Actn1 was significantly altered after TCDD treatment. These findings show that the exposure to TCDD has an impact on alternative-splicing, and suggest a new avenue for understanding TCDD-mediated toxicity and pathogenesis.

Preliminary In Vitro Effects of CD8+ T Lymphocyte Specific for the CD20 Alternative Splicing D393-CD20 Peptide Expressed on Burkitt Lymphoma Cells

The effective discovery of clinically relevant tumor antigens holds a fundamental role for the development of new diagnostic tools and anticancer immunotherapies. D393-CD20 mRNA is absent from normal resting B cells but present in various malignant or transformed B cells. CD8+T lymphocytes play a central role in immunity to cancer. In this study, we want use from T CD8+ against D393-CD20 for effect in RAMOS cell line. After isolation and expanding of specific TCD8 + Lymphocyte against D393-CD20 antigen, for examining the effect of specialized T lymphocyte clone of D393-CD20 antigen on RAMOS cell line, we co-cultured them together, and the rate of apoptosis were examined by flow cytometry and cytotoxicity techniques by using MTT technique. We observed that specialized TCD8+ lymphocyte of D393-CD20 antigen can induce apoptosis in malignant B-lymphocytes, and this antigen can be a proper target for immunotherapy.

Alternative splicing events during adipogenesis from hMSCs

Adipogenesis, the developmental process of progenitor-cell differentiating into adipocytes, leads to fat metabolic disorders. Alternative splicing (AS), a ubiquitous regulatory mechanism of gene expression, allows the generation of more than one unique messenger RNA (mRNA) species from a single gene. Till now, alternative splicing events during adipogenesis from human mesenchymal stem cells (hMSCs) are not yet fully elucidated. We performed RNA-Seq coupled with bioinformatics analysis to identify the differentially expressed AS genes and events during adipogenesis from hMSCs. A global survey separately identified 1262, 1181, 1167, and 1227 ASE involved in the most common types of AS including cassette exon, alt3, and alt5, especially with cassette exon the most prevalent, at 7, 14, 21, and 28 days during adipogenesis. Interestingly, 122 differentially expressed ASE referred to 118 genes, and the three genes including ACTN1 (alt3 and cassette), LRP1 (alt3 and alt5), and LTBP4 (cassette, cassette_multi, and unknown), appeared in multiple AS types of ASE during adipogenesis. Except for all the identified ASE of LRP1 occurred in the extracellular topological domain, alt3 (84) in transmembrane domain significantly differentially expressed was the potential key event during adipogenesis. Overall, we have, for the first time, conducted the global transcriptional profiling during adipogenesis of hMSCs to identify differentially expressed ASE and ASE-related genes. This finding would provide extensive ASE as the regulator of adipogenesis and the potential targets for future molecular research into adipogenesis-related metabolic disorders.

Colorectal cancer vaccines: Tumor-associated antigens vs neoantigens

Therapeutic options for the treatment of colorectal cancer (CRC) are diverse but still not always satisfying. Recent success of immune checkpoint inhibition treatment for the subgroup of CRC patients suffering from hyper-mutated tumors suggests a permanent role of immune therapy in the clinical management of CRC. Substantial improvement in treatment outcome could be achieved by development of efficient patient-individual CRC vaccination strategies. This mini-review summarizes the current knowledge on the two general classes of targets: tumor-associated antigens (TAAs) and tumor-specific antigens. TAAs like carcinoembryonic antigen and melanoma associated antigen are present in and shared by a subgroup of patients and a variety of clinical studies examined the efficacy of different TAA-derived peptide vaccines. Combinations of several TAAs as the next step and the development of personalized TAA-based peptide vaccines are discussed. Improvements of peptide-based vaccines achievable by adjuvants and immune-stimulatory chemotherapeutics are highlighted. Finally, we sum up clinical studies using tumor-specific antigens - in CRC almost exclusively neoantigens - which revealed promising results; particularly no severe adverse events were reported so far. Critical progress for clinical outcomes can be expected by individualizing neoantigen-based peptide vaccines and combining them with immune-stimulatory chemotherapeutics and immune checkpoint inhibitors. In light of these data and latest developments, truly personalized neoantigen-based peptide vaccines can be expected to fulfill modern precision medicine’s requirements and will manifest as treatment pillar for routine clinical management of CRC.

Gene expression signatures for HOXA4, HOXA9, and HOXD10 reveal alterations in transcriptional regulatory networks in colon cancer

We previously reported that HOXA4, HOXA9, and HOXD10 are selectively expressed in colonic stem cells (SCs) and their overexpression contributes to colorectal cancer (CRC). Our goals here were to determine how these HOX genes are transcriptionally regulated and whether transcriptional dysregulation of HOX genes occurs in CRC. Accordingly, we used correlation analysis to identify genes that are expression-correlated or anticorrelated with HOXA4, HOXA9, and HOXD10. We then used Gene Ontology (GO) analysis to functionally classify these genes. The GO results for both HOXA4 and HOXD10 correlated gene sets for normal colon and CRC show functions mostly classified as developmental, transcriptional regulation, and DNA binding. This raised the question: Are these gene sets regulated by the same transcription factors (TFs)? Consequently, we used promoter analysis and interaction network toolset (PAINT) to identify commonly shared transcription response elements. The results indicated that completely different sets of TFs coregulate HOXA4 and HOXD10 (but not HOXA9) and their expression-correlated genes. And predicted TFs are altered in CRC compared with normal colon. Taken together, analysis of gene signatures correlated with expression of HOXA4 and HOXD10 indicates how these HOX genes are: (a) transcriptionally regulated in the normal colon; (b) dysregulated in CRC. This discovery provides a mechanism for targeting CRC SCs.

Comparison of RNA-seq and microarray platforms for splice event detection using a cross-platform algorithm

Background

RNA-seq is a reference technology for determining alternative splicing at genome-wide level. Exon arrays remain widely used for the analysis of gene expression, but show poor validation rate with regard to splicing events. Commercial arrays that include probes within exon junctions have been developed in order to overcome this problem.

We compare the performance of RNA-seq (Illumina HiSeq) and junction arrays (Affymetrix Human Transcriptome array) for the analysis of transcript splicing events. Three different breast cancer cell lines were treated with CX-4945, a drug that severely affects splicing. To enable a direct comparison of the two platforms, we adapted EventPointer, an algorithm that detects and labels alternative splicing events using junction arrays, to work also on RNA-seq data. Common results and discrepancies between the technologies were validated and/or resolved by over 200 PCR experiments.

Results

As might be expected, RNA-seq appears superior in cases where the technologies disagree and is able to discover novel splicing events beyond the limitations of physical probe-sets. We observe a high degree of coherence between the two technologies, however, with correlation of EventPointer results over 0.90. Through decimation, the detection power of the junction arrays is equivalent to RNA-seq with up to 60 million reads.

Conclusions

Our results suggest, therefore, that exon-junction arrays are a viable alternative to RNA-seq for detection of alternative splicing events when focusing on well-described transcriptional regions.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5082-2) contains supplementary material, which is available to authorized users.

Profiles of alternative splicing in colorectal cancer and their clinical significance: A study based on large-scale sequencing data

Background

Alternative splicing (AS), as a potent and pervasive mechanism of transcriptional regulatory, expands the genome's coding capacity and involves in the initiation and progression of cancer. Systematic analysis of alternative splicing in colorectal cancer (CRC) is lacking and greatly needed.

Methods

RNA-Seq data and corresponding clinical information of CRC cohort were downloaded from the TCGA data portal. Then, a java application, known as SpliceSeq, was used to evaluate the RNA splicing patterns and calculate the Percent Spliced In (PSI) value. Differently expressed AS events (DEAS) were identified based on PSI value between paired CRC and adjacent tissues. DEAS and its splicing networks were further analyzed by bioinformatics methods. Kaplan-Meier, Cox proportional regression and unsupervised clustering analysis were used to evaluate the association between DEAS and patients' clinical features.

Results

After strict filtering, a total of 34,334 AS events were identified, among which 421 AS events were found expressed differently. Parent genes of these DEAS play a important role in regulating CRC-related processes such as protein kinase activity (FDR<0.0001), PI3K-Akt signaling pathway (FDR = 0.0024) and p53 signaling pathway (FDR = 0.0143). 37 DEAS events were found to be associated with OS, and 68 DEAS events were found to be associated with DFS. Stratifying patients according to the PSI value of AT in CXCL12 and RI in CSTF3 formed significant Kaplan-Meier curves in both OS and DFS survival analysis. Unsupervised clustering analysis using DEAS revealed four clusters with distinct survival patterns, and associated with consensus molecular subtypes.

Conclusions

Large differences of AS events in CRC appear to exist, and these differences are likely to be important determinants of both prognosis and biological regulation. Our identified CRC-related AS events and uncovered splicing networks are valuable in deciphering the underlying mechanisms of AS in CRC, and provide clues of therapeutic targets to further validations.

Methylmercury exposure alters RNA splicing in human neuroblastoma SK-N-SH cells: Implications from proteomic and post-transcriptional responses

The neurotoxic effects of methylmercury (MeHg) have been intensively studied. However, the molecular mechanisms responsible for the neurotoxicity of MeHg are not fully understood. To decipher these mechanisms, proteomic and high-throughput mRNA sequencing (RNA-seq) technique were utilized, comprehensively evaluating the cellular responses of human neuroblastoma SK-N-SH cells to MeHg exposure. Proteomic results revealed that MeHg exposure interfered with RNA splicing via splicesome, along with the known molecular mechanisms of mercury-related neurotoxicity (e.g. oxidative stress, protein folding, immune system processes, and cytoskeletal organization). The effects of MeHg on RNA splicing were further verified using RNA-seq. Compared to control, a total of 658 aberrant RNA alternative splicing (AS) events were observed after MeHg exposure. Proteomics and RNA-seq results also demonstrated that mercury chloride (HgCl₂) influenced the expression levels of several RNA splicing related proteins and 676 AS events compared to control. These results suggested that RNA splicing could be a new molecular mechanism involved in MeHg and HgCl₂ neurotoxicity.

The Usage of Exon-Exon Splice Junctions for the Detection of Alternative Splicing using the REIDS model

Alternative gene splicing is a common phenomenon in which a single gene gives rise to multiple transcript isoforms. The process is strictly guided and involves a multitude of proteins and regulatory complexes. Unfortunately, aberrant splicing events have been linked to genetic disorders. Therefore, understanding mechanisms of alternative splicing regulation and differences in splicing events between diseased and healthy tissues is crucial in advancing personalized medicine and drug developments. We propose a linear mixed model, Random Effects for the Identification of Differential Splicing (REIDS), for the identification of alternative splicing events using Human Transcriptome Arrays (HTA). For each exon, a splicing score is calculated based on two scores, an exon score and an array score. The junction information is used to rank the identified exons from strongly confident to less confident candidates for alternative splicing. The design of junctions was also discussed to highlight the complexity of exon-exon and exon-junction interactions. Based on a list of Rt-PCR validated probe sets, REIDS outperforms AltAnalyze and iGems in the % recall rate.

The awakening of the CDK10/Cyclin M protein kinase

Cyclin-dependent kinases (CDKs) play important roles in the control of fundamental cellular processes. Some of the most characterized CDKs are considered to be pertinent therapeutic targets for cancers and other diseases, and first clinical successes have recently been obtained with CDK inhibitors. Although discovered in the pre-genomic era, CDK10 attracted little attention until it was identified as a major determinant of resistance to endocrine therapy for breast cancer. In some studies, CDK10 has been shown to promote cell proliferation whereas other studies have revealed a tumor suppressor function. The recent discovery of Cyclin M as a CDK10 activating partner has allowed the unveiling of a protein kinase activity against the ETS2 oncoprotein, whose degradation is activated by CDK10/Cyclin M-mediated phosphorylation. CDK10/Cyclin M has also been shown to repress ciliogenesis and to maintain actin network architecture, through the phoshorylation of the PKN2 protein kinase and the control of RhoA stability. These findings shed light on the molecular mechanisms underlying STAR syndrome, a severe human developmental genetic disorder caused by mutations in the Cyclin M coding gene. They also pave the way to a better understanding of the role of CDK10/Cyclin M in cancer.

High expression of PTBP1 promote invasion of colorectal cancer by alternative splicing of cortactin

Polypyrimidine tract-binding protein 1 (PTBP1) involving in almost all steps of mRNA regulation including alternative splicing metabolism during tumorigenesis due to its RNA-binding activity. Initially, we found that high expressed PTBP1 and poor prognosis was interrelated in colorectal cancer (CRC) patients with stages II and III CRC, which widely different in prognosis and treatment, by immunohistochemistry. PTBP1 was also upregulated in colon cancer cell lines. In our study, knockdown of PTBP1 by siRNA transfection decreased cell proliferation and invasion in vitro. Denovirus shRNA knockdown of PTBP1 inhibited colorectal cancer growth in vivo. Furthermore, PTBP1 regulates alternative splicing of many target genes involving in tumorgenesis in colon cancer cells. We confirmed that the splicing of cortactin exon 11 which was only contained in cortactin isoform-a, as a PTBP1 target. Knockdown of PTBP1 decreased the expression of cortactin isoform-a by exclusion of exon 11. Also the mRNA levels of PTBP1 and cortactin isoform-a were cooperatively expressed in colorectal cancer tissues. Knocking down cortactin isoform-a significantly decreased cell migration and invasion in colorectal cancer cells. Overexpression of cortactin isoform-a could rescue PTBP1-knockdown effect of cell motility. In summary the study revealed that PTBP1 facilitates colorectal cancer migration and invasion activities by inclusion of cortactin exon 11.

Alternative transcription of a shorter, non-anti-angiogenic thrombospondin-2 variant in cancer-associated blood vessels

Thrombospondin-2 (TSP2) is an anti-angiogenic matricellular protein that inhibits tumor growth and angiogenesis. Tumor-associated blood vascular endothelial cells (BECs) were isolated from human invasive bladder cancers and from matched normal bladder tissue by immuno-laser capture microdissection. Exon expression profiling analyses revealed a particularly high expression of a short TSP2 transcript containing only the last 9 (3′) exons of the full-length TSP2 transcript. Using 5′ and 3′ RACE (rapid amplification of cDNA ends) and Sanger sequencing, we confirmed the existence of the shorter transcript of TSP2 (sTSP2) and determined its sequence which completely lacked the anti-angiogenic thrombospondin type 1 repeats domain. The largest open reading frame predicted within the transcript comprises 209 amino acids and matches almost completely the C-terminal lectin domain of full-length TSP2. We produced recombinant sTSP2 and found that unlike the full-length TSP2, sTSP2 did not inhibit vascular endothelial growth factor-A-induced proliferation of cultured human BECs, but in contrast when combined with TSP2 blocked the inhibitory effects of TSP2 on BEC proliferation. In vivo studies with stably transfected A431 squamous cell carcinoma cells revealed that full-length TSP2, but not sTSP2, inhibited tumor growth and angiogenesis. This study reveals that the transcriptional program of tumor stromal cells can change to transcribe a new version of an endogenous angiogenesis inhibitor that has lost its anti-angiogenic activity.

Differentially expressed alternatively spliced genes in skeletal muscle from cancer patients with cachexia

BACKGROUND

Alternative splicing (AS) is a post-transcriptional gene regulatory mechanism that contributes to proteome diversity. Aberrant splicing mechanisms contribute to various cancers and muscle-related conditions such as Duchenne muscular dystrophy. However, dysregulation of AS in cancer cachexia (CC) remains unexplored. Our objectives were (i) to profile alternatively spliced genes (ASGs) on a genome-wide scale and (ii) to identify differentially expressed alternatively spliced genes (DASGs) associated with CC.

METHODS

Rectus abdominis muscle biopsies obtained from cancer patients were stratified into cachectic cases (n = 21, classified based on International consensus diagnostic framework for CC) and non-cachectic controls (n = 19, weight stable cancer patients). Human transcriptome array 2.0 was used for profiling ASGs using the total RNA isolated from muscle biopsies. Representative DASG signatures were validated using semi-quantitative RT-PCR.

RESULTS

We identified 8960 ASGs, of which 922 DASGs (772 up-regulated and 150 down-regulated) were identified at ≥1.4 fold-change and P < 0.05. Representative DASGs validated by semi-quantitative RT-PCR confirmed the primary findings from the human transcriptome arrays. Identified DASGs were associated with myogenesis, adipogenesis, protein ubiquitination, and inflammation. Up to 10% of the DASGs exhibited cassette exon (exon included or skipped) as a predominant form of AS event. We also observed other forms of AS events such as intron retention, alternate promoters.

CONCLUSIONS

Overall, we have, for the first time, conducted global profiling of muscle tissue to identify DASGs associated with CC. The mechanistic roles of the identified DASGs in CC pathophysiology using model systems is warranted, as well as replication of findings in independent cohorts.

Simultaneous genome-wide gene expression and transcript isoform profiling in the human malaria parasite

Gene expression DNA microarrays have been vital for characterizing whole-genome transcriptional profiles. Nevertheless, their effectiveness relies heavily on the accuracy of genome sequences, the annotation of gene structures, and the sequence-dependent performance of individual probes. Currently available gene expression arrays for the malaria parasite Plasmodium falciparum rely on an average of 2 probes per gene, usually positioned near the 3′ end of genes; consequently, existing designs are prone to measurement bias and cannot capture complexities such as the occurrence of transcript isoforms arising from alternative splicing or alternative start/ stop sites. Here, we describe two novel gene expression arrays with exon-focused probes designed with an average of 12 and 20 probes spanning each gene. This high probe density minimizes signal noise inherent in probe-to-probe sequence-dependent hybridization intensity. We demonstrate that these exon arrays accurately profile genome-wide expression, comparing favorably to currently available arrays and RNA-seq profiling, and can detect alternatively spliced transcript isoforms as well as non-coding RNAs (ncRNAs). Of the 964 candidate alternate splicing events from published RNA-seq studies, 162 are confirmed using the exon array. Furthermore, the exon array predicted 330 previously unidentified alternate splicing events. Gene expression microarrays continue to offer a cost-effective alternative to RNA-seq for the simultaneous monitoring of gene expression and alternative splicing events. Microarrays may even be preferred in some cases due to their affordability and the rapid turn-around of results when hundreds of samples are required for fine-scale systems biology investigations, including the monitoring of the networks of gene co-expression in the emergence of drug resistance.

Transcriptional signature of lymphoblastoid cell lines of BRCA1, BRCA2 and non-BRCA1/2 high risk breast cancer families

Approximately 25% of hereditary breast cancer cases are associated with a strong familial history which can be explained by mutations in BRCA1 or BRCA2 and other lower penetrance genes. The remaining high-risk families could be classified as BRCAX (non-BRCA1/2) families.

Gene expression involving alternative splicing represents a well-known mechanism regulating the expression of multiple transcripts, which could be involved in cancer development. Thus using RNA-seq methodology, the analysis of transcriptome was undertaken to potentially reveal transcripts implicated in breast cancer susceptibility and development.

RNA was extracted from immortalized lymphoblastoid cell lines of 117 women (affected and unaffected) coming from BRCA1, BRCA2 and BRCAX families. Anova analysis revealed a total of 95 transcripts corresponding to 85 different genes differentially expressed (Bonferroni corrected p-value <0.01) between those groups. Hierarchical clustering allowed distinctive subgrouping of BRCA1/2 subgroups from BRCAX individuals. We found 67 transcripts, which could discriminate BRCAX from BRCA1/BRCA2 individuals while 28 transcripts discriminate affected from unaffected BRCAX individuals.

To our knowledge, this represents the first study identifying transcripts differentially expressed in lymphoblastoid cell lines from major classes of mutation-related breast cancer subgroups, namely BRCA1, BRCA2 and BRCAX. Moreover, some transcripts could discriminate affected from unaffected BRCAX individuals, which could represent potential therapeutic targets for breast cancer treatment.

Splice variants as novel targets in pancreatic ductal adenocarcinoma

Despite a wealth of genomic information, a comprehensive alternative splicing (AS) analysis of pancreatic ductal adenocarcinoma (PDAC) has not been performed yet. In the present study, we assessed whole exome-based transcriptome and AS profiles of 43 pancreas tissues using Affymetrix exon array. The AS analysis of PDAC indicated on average two AS probe-sets (ranging from 1-28) in 1,354 significantly identified protein-coding genes, with skipped exon and alternative first exon being the most frequently utilised. In addition to overrepresented extracellular matrix (ECM)-receptor interaction and focal adhesion that were also seen in transcriptome differential expression (DE) analysis, Fc gamma receptor-mediated phagocytosis and axon guidance AS genes were also highly represented. Of note, the highest numbers of AS probe-sets were found in collagen genes, which encode the characteristically abundant stroma seen in PDAC. We also describe a set of 37 'hypersensitive' genes which were frequently targeted by somatic mutations, copy number alterations, DE and AS, indicating their propensity for multidimensional regulation. We provide the most comprehensive overview of the AS landscape in PDAC with underlying changes in the spliceosomal machinery. We also collate a set of AS and DE genes encoding cell surface proteins, which present promising diagnostic and therapeutic targets in PDAC.

A random effects model for the identification of differential splicing (REIDS) using exon and HTA arrays

Background

Alternative gene splicing is a common phenomenon in which a single gene gives rise to multiple transcript isoforms. The process is strictly guided and involves a multitude of proteins and regulatory complexes. Unfortunately, aberrant splicing events do occur which have been linked to genetic disorders, such as several types of cancer and neurodegenerative diseases (Fan et al., Theor Biol Med Model 3:19, 2006). Therefore, understanding the mechanism of alternative splicing and identifying the difference in splicing events between diseased and healthy tissue is crucial in biomedical research with the potential of applications in personalized medicine as well as in drug development.

Results

We propose a linear mixed model, Random Effects for the Identification of Differential Splicing (REIDS), for the identification of alternative splicing events. Based on a set of scores, an exon score and an array score, a decision regarding alternative splicing can be made. The model enables the ability to distinguish a differential expressed gene from a differential spliced exon. The proposed model was applied to three case studies concerning both exon and HTA arrays.

Conclusion

The REIDS model provides a work flow for the identification of alternative splicing events relying on the established linear mixed model. The model can be applied to different types of arrays.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1687-8) contains supplementary material, which is available to authorized users.

hnRNP I regulates neonatal immune adaptation and prevents colitis and colorectal cancer

The intestinal epithelium plays a critical role in host-microbe homeostasis by sensing gut microbes and subsequently initiating proper immune responses. During the neonatal stage, the intestinal epithelium is under immune repression, allowing the transition for newborns from a relatively sterile intra-uterine environment to one that is rich in foreign antigens. The mechanism underlying such immune repression remains largely unclear, but involves downregulation of IRAK1 (interleukin-1 receptor-associated kinase), an essential component of toll-like receptor-mediated NF-κB signaling. We report here that heterogeneous nuclear ribonucleoprotein I (hnRNPI), an RNA binding protein, is essential for regulating neonatal immune adaptation. We generated a mouse model in which hnRNPI is ablated specifically in the intestinal epithelial cells, and characterized intestinal defects in the knockout mice. We found that loss of hnRNPI function in mouse intestinal epithelial cells results in early onset of spontaneous colitis followed by development of invasive colorectal cancer. Strikingly, the epithelium-specific hnRNPI knockout neonates contain aberrantly high IRAK1 protein levels in the colons and fail to develop immune tolerance to environmental microbes. Our results demonstrate that hnRNPI plays a critical role in establishing neonatal immune adaptation and preventing colitis and colorectal cancer.

Precisely controlled host-microbe interactions in the gastrointestinal tract are crucial for human overall health and well-being. Dysregulated host responses to gut microbiota are the major cause of autoimmune diseases, inflammatory disorders and cancers. The intestinal epithelium lines the gastrointestinal tract and plays a critical role in sensing gut microbes and subsequently developing a balance of immune tolerance and active immune responses. During the neonatal stage, the immune system in the gastrointestinal tract must be temporally suppressed to accommodate the large number of newly arrived microbes. This process is known as neonatal immune adaptation, and is critical for the establishment of proper host- microbe interactions. We studied the function of hnRNPI in the intestinal epithelium by genetically ablating it in the intestinal epithelial cells of mouse. We found that loss of hnRNPI in intestinal epithelial cells disrupts neonatal immune adaptation, resulting in spontaneous colitis and early onset of invasive colorectal cancer. We show that hnRNPI is required for the neonatal immune suppression through decreasing the protein level of IRAK1, an essential component of toll-like receptor-mediated NF-κB signaling. Our studies demonstrate a critical role of hnRNPI in establishing neonatal immune adaptation and preventing colitis and colorectal cancer.

Splicing imbalances in basal-like breast cancer underpin perturbation of cell surface and oncogenic pathways and are associated with patients' survival

Despite advancements in the use of transcriptional information to understand and classify breast cancers, the contribution of splicing to the establishment and progression of these tumours has only recently starting to emerge. Our work explores this lesser known landscape, with special focus on the basal-like breast cancer subtype where limited therapeutic opportunities and no prognostic biomarkers are currently available. Using ExonArray analysis of 176 breast cancers and 9 normal breast tissues we demonstrate that splicing levels significantly contribute to the diversity of breast cancer molecular subtypes and explain much of the differences compared with normal tissues. We identified pathways specifically affected by splicing imbalances whose perturbation would be hidden from a conventional gene-centric analysis of gene expression. We found that a large fraction of them involve cell-to-cell communication, extracellular matrix and transport, as well as oncogenic and immune-related pathways transduced by plasma membrane receptors. We identified 247 genes in which splicing imbalances are associated with clinical patients’ outcome, whilst no association was detectable at the gene expression level. These include the signaling gene TGFBR1, the proto-oncogene MYB as well as many immune-related genes such as CCR7 and FCRL3, reinforcing evidence for a role of immune components in influencing breast cancer patients’ prognosis.

Phenethyl isothiocyanate alters the gene expression and the levels of protein associated with cell cycle regulation in human glioblastoma GBM 8401 cells

Glioblastoma is the most common and aggressive primary brain malignancy. Phenethyl isothiocyanate (PEITC), a member of the isothiocyanate family, can induce apoptosis in many human cancer cells. Our previous study disclosed that PEITC induces apoptosis through the extrinsic pathway, dysfunction of mitochondria, reactive oxygen species (ROS)-induced endoplasmic reticulum (ER) stress, and intrinsic (mitochondrial) pathway in human brain glioblastoma multiforme (GBM) 8401 cells. To the best of our knowledge, we first investigated the effects of PEITC on the genetic levels of GBM 8401 cells in vitro. PEITC may induce G0/G1 cell-cycle arrest through affecting the proteins such as cdk2, cyclin E, and p21 in GBM 8401 cells. Many genes associated with cell-cycle regulation of GBM 8401 cells were changed after PEITC treatment: 48 genes were upregulated and 118 were downregulated. The cell-division cycle protein 20 (CDC20), Budding uninhibited by benzimidazole 1 homolog beta (BUB1B), and cyclin B1 were downregulated, and clusterin was upregulated in GBM 8401 cells treated with PEITC. These changes of gene expression can provide the effects of PEITC on the genetic levels and potential biomarkers for glioblastoma. © 2015 Wiley Periodicals, Inc. Environ Toxicol 32: 176-187, 2017.

Differential expression of alternatively spliced transcripts related to energy metabolism in colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is one of the most common malignant tumors worldwide. CRC molecular pathogenesis is heterogeneous and may be followed by mutations in oncogenes and tumor suppressor genes, chromosomal and microsatellite instability, alternative splicing alterations, hypermethylation of CpG islands, oxidative stress, impairment of different signaling pathways and energy metabolism. In the present work, we have studied the alterations of alternative splicing patterns of genes related to energy metabolism in CRC.

RESULTS

Using CrossHub software, we analyzed The Cancer Genome Atlas (TCGA) RNA-Seq datasets derived from colon tumor and matched normal tissues. The expression of 1014 alternative mRNA isoforms involved in cell energy metabolism was examined. We found 7 genes with differentially expressed alternative transcripts whereas overall expression of these genes was not significantly altered in CRC. A set of 8 differentially expressed transcripts of interest has been validated by qPCR. These eight isoforms encoded by OGDH, COL6A3, ICAM1, PHPT1, PPP2R5D, SLC29A1, and TRIB3 genes were up-regulated in colorectal tumors, and this is in concordance with the bioinformatics data. The alternative transcript NM_057167 of COL6A3 was also strongly up-regulated in breast, lung, prostate, and kidney tumors. Alternative transcript of SLC29A1 (NM_001078177) was up-regulated only in CRC samples, but not in the other tested tumor types.

CONCLUSIONS

We identified tumor-specific expression of alternative spliced transcripts of seven genes involved in energy metabolism in CRC. Our results bring new knowledge on alternative splicing in colorectal cancer and suggest a set of mRNA isoforms that could be used for cancer diagnosis and development of treatment methods.

Transcriptional Elongation Regulator 1 Affects Transcription and Splicing of Genes Associated with Cellular Morphology and Cytoskeleton Dynamics and Is Required for Neurite Outgrowth in Neuroblastoma Cells and Primary Neuronal Cultures

TCERG1 is a highly conserved human protein implicated in interactions with the transcriptional and splicing machinery that is associated with neurodegenerative disorders. Biochemical, neuropathological, and genetic evidence suggests an important role for TCERG1 in Huntington's disease (HD) pathogenesis. At present, the molecular mechanism underlying TCERG1-mediated neuronal effects is unknown. Here, we show that TCERG1 depletion led to widespread alterations in mRNA processing that affected different types of alternative transcriptional or splicing events, indicating that TCERG1 plays a broad role in the regulation of alternative splicing. We observed considerable changes in the transcription and alternative splicing patterns of genes involved in cytoskeleton dynamics and neurite outgrowth. Accordingly, TCERG1 depletion in the neuroblastoma SH-SY5Y cell line and primary mouse neurons affected morphogenesis and resulted in reduced dendritic outgrowth, with a major effect on dendrite ramification and branching complexity. These defects could be rescued by ectopic expression of TCERG1. Our results indicate that TCERG1 affects expression of multiple mRNAs involved in neuron projection development, whose misregulation may be involved in TCERG1-linked neurological disorders.

Characterization of visceral and subcutaneous adipose tissue transcriptome in pregnant women with and without spontaneous labor at term: implication of alternative splicing in the metabolic adaptations of adipose tissue to parturition

OBJECTIVE

The aim of this study was to determine gene expression and splicing changes associated with parturition and regions (visceral vs. subcutaneous) of the adipose tissue of pregnant women.

STUDY DESIGN

The transcriptome of visceral and abdominal subcutaneous adipose tissue from pregnant women at term with (n=15) and without (n=25) spontaneous labor was profiled with the Affymetrix GeneChip Human Exon 1.0 ST array. Overall gene expression changes and the differential exon usage rate were compared between patient groups (unpaired analyses) and adipose tissue regions (paired analyses). Selected genes were tested by quantitative reverse transcription-polymerase chain reaction.

RESULTS

Four hundred and eighty-two genes were differentially expressed between visceral and subcutaneous fat of pregnant women with spontaneous labor at term (q-value <0.1; fold change >1.5). Biological processes enriched in this comparison included tissue and vasculature development as well as inflammatory and metabolic pathways. Differential splicing was found for 42 genes [q-value <0.1; differences in Finding Isoforms using Robust Multichip Analysis scores >2] between adipose tissue regions of women not in labor. Differential exon usage associated with parturition was found for three genes (LIMS1, HSPA5, and GSTK1) in subcutaneous tissues.

CONCLUSION

We show for the first time evidence of implication of mRNA splicing and processing machinery in the subcutaneous adipose tissue of women in labor compared to those without labor.

A method for identifying discriminative isoform-specific peptides for clinical proteomics application

Background

Clinical proteomics application aims at solving a specific clinical problem within the context of a clinical study. It has been growing rapidly in the field of biomarker discovery, especially in the area of cancer diagnostics. Until recently, protein isoform has not been viewed as a new class of early diagnostic biomarkers for clinical proteomics. A protein isoform is one of different forms of the same protein. Different forms of a protein may be produced from single-nucleotide polymorphisms (SNPs), alternative splicing, or post-translational modifications (PTMs). Previous studies have shown that protein isoforms play critical roles in tumorigenesis, disease diagnosis, and prognosis. Identifying and characterizing protein isoforms are essential to the study of molecular mechanisms and early detection of complex diseases such as breast cancer.

However, there are limitations with traditional methods such as EST sequencing, Microarray profiling (exon array, Exon-exon junction array), mRNA next-generation sequencing used for protein isoform determination: 1) not in the protein level, 2) no connectivity about connection of nonadjacent exons, 3) no SNPs and PTMs, and 4) low reproducibility. Moreover, there exist the computational challenges of clinical proteomics studies: 1) low sensitivity of instruments, 2) high data noise, and 3) high variability and low repeatability, although recent advances in clinical proteomics technology, LC-MS/MS proteomics, have been used to identify candidate molecular biomarkers in diverse range of samples, including cells, tissues, serum/plasma, and other types of body fluids.

Results

Therefore, in the paper, we presented a peptidomics method for identifying cancer-related and isoform-specific peptide for clinical proteomics application from LC-MS/MS. First, we built a Peptidomic Database of Human Protein Isoforms, then created a peptidomics approach to perform large-scale screen of breast cancer-associated alternative splicing isoform markers in clinical proteomics, and lastly performed four kinds of validations: biological validation (explainable index), exon array, statistical validation of independent samples, and extensive pathway analysis.

Conclusions

Our results showed that alternative splicing isoform makers can act as independent markers of breast cancer and that the method for identifying cancer-specific protein isoform biomarkers from clinical proteomics application is an effective one for increasing the number of identified alternative splicing isoform markers in clinical proteomics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2907-8) contains supplementary material, which is available to authorized users.

Systems Pharmacogenomics Finds RUNX1 Is an Aspirin-Responsive Transcription Factor Linked to Cardiovascular Disease and Colon Cancer

Aspirin prevents cardiovascular disease and colon cancer; however aspirin's inhibition of platelet COX-1 only partially explains its diverse effects. We previously identified an aspirin response signature (ARS) in blood consisting of 62 co-expressed transcripts that correlated with aspirin's effects on platelets and myocardial infarction (MI). Here we report that 60% of ARS transcripts are regulated by RUNX1 - a hematopoietic transcription factor - and 48% of ARS gene promoters contain a RUNX1 binding site. Megakaryocytic cells exposed to aspirin and its metabolite (salicylic acid, a weak COX-1 inhibitor) showed up regulation in the RUNX1 P1 isoform and MYL9, which is transcriptionally regulated by RUNX1. In human subjects, RUNX1 P1 expression in blood and RUNX1-regulated platelet proteins, including MYL9, were aspirin-responsive and associated with platelet function. In cardiovascular disease patients RUNX1 P1 expression was associated with death or MI. RUNX1 acts as a tumor suppressor gene in gastrointestinal malignancies. We show that RUNX1 P1 expression is associated with colon cancer free survival suggesting a role for RUNX1 in aspirin's protective effect in colon cancer. Our studies reveal an effect of aspirin on RUNX1 and gene expression that may additionally explain aspirin's effects in cardiovascular disease and cancer.

Genome-wide profiling of chemoradiation‑induced changes in alternative splicing in colon cancer cells

Alternative splicing is a key mechanism that regulates protein diversity and has been found to be associated with colon cancer progression and metastasis. However, the function of alternative splicing in chemoradiation‑resistant colon cancer remains elusive. In this study, we constructed a chemoradiation‑resistant colon cancer cell line. Through RNA-sequencing of normal and chemoradiation‑resistant colon cancer cells (HCT116), we found 818 genes that were highly expressed in the normal HCT116 cells, whereas 285 genes were highly expressed in the chemoradiation-resistant HCT116 (RCR-HCT116) cells. Gene ontology (GO) analysis showed that genes that were highly expressed in the HCT116 cells were enriched in GO categories related to cell cycle and cell division, whereas genes that were highly expressed in the RCR-HCT116 cells were associated with regulation of system processes and response to wounding. Analysis of alternative splicing events revealed that exon skipping was significantly increased in the chemoradiation‑resistant colon cancer cells. Moreover, we identified 323 alternative splicing events in 293 genes that were significantly different between the two different HCT116 cell types. These alternative splicing‑related genes were clustered functionally into several groups related with DNA replication, such as deoxyribonucleotide metabolic/catabolic processes, response to DNA damage stimulus and helicase activity. These findings enriched our knowledge by elucidating the function of alternative splicing in chemoradiation-resistant colon cancer.

Stroma derived COL6A3 is a potential prognosis marker of colorectal carcinoma revealed by quantitative proteomics

Colorectal cancer (CRC) represents the third most common cancer in males and second in females worldwide. Here, we performed a quantitative 8-plex iTRAQ proteomics analysis of the secreted proteins from five colonic fibroblast cultures and three colon cancer epithelial cell lines. We identified 1114 proteins at 0% FDR, including 587 potential secreted proteins. We further recognized 116 fibroblast-enriched proteins which were significantly associated with cell movement, angiogenesis, proliferation and wound healing, and 44 epithelial cell-enriched proteins. By interrogation of Oncomine database, we found that 20 and 8 fibroblast-enriched proteins were up- and downregulated in CRC, respectively. Western blots confirmed the fibroblast-specific secretion of filamin C, COL6A3, COL4A1 and spondin-2. Upregulated mRNA and stroma expression of COL6A3 in CRC, which were revealed by Oncomine analyses and tissue-microarray-immunohistochemistry, indicated poor prognosis. COL6A3 expression was significantly associated with Dukes stage, T stage, stage, recurrence and smoking status. Circulating plasma COL6A3 in CRC patients was upregulated significantly comparing with healthy peoples. Receiver operating characteristic curve analysis revealed that COL6A3 has better predictive performance for CRC with an area under the curve of 0.885 and the best sensitivity/specificity of 92.9%/81.3%. Thus we demonstrated that COL6A3 was a potential diagnosis and prognosis marker of CRC.