RNA-Seq reveals differentially expressed isoforms and novel splice variants in buccal mucosal cancer

Splicing Machinery Is Impaired in Oral Squamous Cell Carcinomas and Linked to Key Pathophysiological Features

Alternative splicing dysregulation is an emerging cancer hallmark, potentially serving as a source of novel diagnostic, prognostic, or therapeutic tools. Inhibitors of the activity of the splicing machinery can exert antitumoral effects in cancer cells. We aimed to characterize the splicing machinery (SM) components in oral squamous cell carcinoma (OSCC) and to evaluate the direct impact of the inhibition of SM-activity on OSCC-cells. The expression of 59 SM-components was assessed using a prospective case-control study of tumor and healthy samples from 37 OSCC patients, and the relationship with clinical and histopathological features was assessed. The direct effect of pladienolide-B (SM-inhibitor) on the proliferation rate of primary OSCC cell cultures was evaluated. A significant dysregulation in several SM components was found in OSCC vs. adjacent-healthy tissues [i.e., 12 out of 59 (20%)], and their expression was associated with clinical and histopathological features of less aggressiveness and overall survival. Pladienolide-B treatment significantly decreased OSCC-cell proliferation. Our data reveal a significantly altered expression of several SM-components and link it to pathophysiological features, reinforcing a potential clinical and pathophysiological relevance of the SM dysregulation in OSCC. The inhibition of SM-activity might be a therapeutic avenue in OSCC, offering a clinically relevant opportunity to be explored.

Identification of a Novel Splice Variant of Neural Cell Adhesion Molecule in Glioblastoma Through Proteogenomics Analysis

Splice variants are known to be important in the pathophysiology of tumors, including the brain cancers. We applied a proteogenomics pipeline to identify splice variants in glioblastoma (GBM, grade IV glioma), a highly malignant brain tumor, using in-house generated mass spectrometric proteomic data and public domain RNASeq dataset. Our analysis led to the identification of a novel exon that maps to the long isoform of Neural cell adhesion molecule 1 (NCAM1), expressed on the surface of glial cells and neurons, important for cell adhesion and cell signaling. The presence of the novel exon is supported with the identification of five peptides spanning it. Additional peptides were also detected in sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) gel separated proteins from GBM patient tissue, underscoring the presence of the novel peptides in the intact brain protein. The novel exon was detected in the RNASeq dataset in 18 of 25 GBM samples and separately validated in additional 10 GBM tumor tissues using quantitative real-time-polymerase chain reaction (qRT-PCR). Both transcriptomic and proteomic data indicate downregulation of NCAM1, including the novel variant, in GBM. Domain analysis of the novel NCAM1 sequence indicates that the insertion of the novel exon contributes extra low-complexity region in the protein that may be important for protein-protein interactions and hence for cell signaling associated with tumor development. Taken together, the novel NCAM1 variant reported in this study exemplifies the importance of future multiomics research and systems biology applications in GBM.

Significantly increased amino acid accumulation in a novel albino branch of the tea plant (Camellia sinensis)

A normal tea plant with one albino branch was discovered. RNA sequencing, albinism phenotype and ultrastructural observations provided a valuable understanding of the albino mechanism in tea plants. Tea plants with a specific color (white or yellow) have been studied extensively. A normal tea plant (Camellia sinensis cv. quntizhong) with one albino branch was discovered in a local tea plantation in Huangshan, Anhui, China. The pure albino leaves on this special branch had accumulated a fairly high content of amino acids, especially theanine (45.31 mg/g DW), and had a low concentration of polyphenols and an extremely low chlorophyll (Chl) content compared with control leaves. Ultrastructural observation of an albino leaf revealed no chloroplasts, whereas it was viable in the control leaf. RNA sequencing and differentially expressed gene (DEG) analysis were performed on the albino leaves and on control leaves from a normal green branch. The related genes involved in theanine and polyphenol biosynthesis were also investigated in this study. DEG expression patterns in Chl biosynthesis or degradation, carotenoid biosynthesis or degradation, chloroplast development, and biosynthesis were influenced in the albino leaves. Chloroplast deletion in albino leaves had probably destroyed the balance of carbon and nitrogen metabolism, leading to a high accumulation of free amino acids and a low concentration of polyphenols in the albino leaves. The obtained results can provide insight into the mechanism underlying this special albino branch phenotype, and are a valuable contribution toward understanding the albino mechanism in tea plants.

Genome-wide isoform-level analysis reveals tumor-specific isoforms for lung adenocarcinoma diagnosis and prognosis

Last decades have witnessed the great progress in exploration of tumor transcriptome. However, most researches were restricted in gene-level expression. mRNA isoforms, especially tumor-specific isoforms have not been fully explored in tumor. Here, by analyzing RNA-seq data derived from hundreds of samples in TCGA projects, we comprehensively characterized the expression variations of mRNA isoforms in adenocarcinoma of lung (LUAD), which is one of leading causes of cancer-related death. Our analysis found that a variety of mRNA isoforms showed differential expression in LUAD tumor samples. Some of them even showed distinct variations compared to their host genes. Further analysis of functional enrichment revealed that up- and down-regulated mRNA isoforms took part in different types of biological process. In addition, we also identified hundreds of isoforms that expressed exclusively in LUAD tumor samples. Furthermore, the expression level of several isoforms, such as uc001kuk.3 and uc003yls.2, could separate tumor patients by overall survival periods. Our study provided new candidates for the diagnosis and prognosis of lung cancer.

Identification of the key genes implicated in the transformation of OLP to OSCC using RNA-sequencing

Oral lichen planus (OLP) is a chronic inflammatory disease that may transform to oral squamous cell carcinoma (OSCC), while its carcinogenesis mechanisms are not entirely clear. This study was designed to identify the important genes involved in the malignant transformation of OLP to OSCC. After RNA-sequencing, the differently expressed genes (DEGs) in OLP vs. normal and OSCC vs. normal groups, respectively, were identified by limma package in R language, and then clustering analysis were conducted by Pheatmap package in R language. Weighed gene co-expression network analysis (WGCNA) was performed for the DEGs to screen disease-associated modules. Using Cytoscape software, co-expression networks were constructed for the genes involved in the modules. Enrichment analysis was conducted for the genes involved in the co-expression networks using GOstat package in R language. Finally, quantitative real-time PCR (qRT-PCR) experiments were conducted to validate the key genes. There were, respectively, 223 and 548 DEGs in OLP vs. normal and OSCC vs. normal groups. WGCNA identified the blue modules for the DEGs in the two groups as disease-associated modules. Moreover, 19 common DEGs (including upregulated BCL9L, PER2 and TSPAN33, and downregulated GMPS and HES1) associated with both OLP and OSCC were identified. In the co-expression networks, BCL9L, HES1, PER2 and TSPAN33 might function in OLP via interactions (such as BCL9L-TSPAN33 and HES1-PER2). qRT-PCR analysis showed that BCL9L, PER2 and TSPAN33 were significantly upregulated, and GMP and HES1 were downregulated. These findings indicated that BCL9L, GMPS, HES1, PER2 and TSPAN33 affected the transformation of OLP to OSCC.

Identification and quantification of novel RNA isoforms in horn cancer of Bos indicus by comprehensive RNA-Seq

Horn cancer (HC) is a squamous cell carcinoma of horn, commonly observed in Bos indicus of the Asian countries. To elucidate the complexity of alternative splicing present in the HC, high-throughput sequencing and analysis of HC and matching horn normal (HN) tissue were carried out. A total of 535,067 and 849,077 reads were analysed after stringent quality filtering for HN and HC, respectively. Cufflinks pipeline for transcriptome analysis revealed 4786 novel splice isoforms comprising 2432 exclusively in HC, 2055 exclusively in HN and 298 in both the conditions. Based on pathway clustering and in silico verification, 102 novel splice isoforms were selected and further analysed with respect to change in protein sequence using Blastp. Finally, fourteen novel splicing events supported both by Cufflinks and UCSC genome browser were selected and confirmed expression by RT-qPCR. Future studies targeted at in-depth characterization of these potential candidate splice isoforms might be helpful in the development of relevant biomarkers for early diagnosis of HC. The results reported in this study refine the available information on transcriptome repertoire of bovine species and boost the research in the line of development of relevant biomarkers for early diagnosis of HC.

De novo Sequencing of the Leaf Transcriptome Reveals Complex Light-Responsive Regulatory Networks in Camellia sinensis cv. Baijiguan

Tea plants (Camellia sinensis L.) possess high genetic diversity that is important for breeding. One cultivar, Baijiguan, exhibits a yellow leaf phenotype, reduced chlorophyll (Chl) content, and aberrant chloroplast structures under high light intensity. In contrast, under low light intensity, the flush shoot from Baijiguan becomes green, the Chl content increases significantly, and the chloroplasts exhibit normal structures. To understand the underlying molecular mechanisms for these observations, we performed de novo transcriptome sequencing and digital gene expression (DGE) profiling using Illumina sequencing technology. De novo transcriptome assembly identified 88,788 unigenes, including 1652 transcription factors from 25 families. In total, 1993 and 2576 differentially expressed genes (DEGs) were identified in Baijiguan plants treated with 3 and 6 days of shade, respectively. Gene Ontology (GO) and pathway enrichment analyses indicated that the DEGs are predominantly involved in the ROS scavenging system, chloroplast development, photosynthetic pigment synthesis, secondary metabolism, and circadian systems. The light-responsive gene POR (protochlorophyllide oxidoreductase) and transcription factor HY5 were identified. Quantitative real-time PCR (qRT-PCR) analysis of 20 selected DEGs confirmed the RNA-sequencing (RNA-Seq) results. Overall, these findings suggest that high light intensity inhibits the expression of photosystem II 10-kDa protein (PsbR) in Baijiguan, thus affecting PSII stability, chloroplast development and chlorophyll biosynthesis.

Comprehensive analysis of the whole transcriptomes from two different pig breeds using RNA-Seq method

Next-generation sequencing RNA-Seq technology is a powerful tool that creates new possibilities for whole-transcriptome analysis. In our study, the RNA-Seq method was applied to analyze global changes in transcriptome from muscle tissue (m. semimembranosus) in two pig breeds (Pietrain and Polish Landrace, PL). The breeds differ in terms of muscularity, growth rate and reproduction traits. Using three different approaches (deseq, cufflinks and edger) and taking into account the most restrictive criteria, 35 genes differentially expressed between Pietrain and PL pigs were identified. In both breeds, the most abundant were transcripts encoding ribosomal and cytoskeletal proteins (TPM3, TCAP, TMOD4, TPM2, TNNC1) and calcium-binding proteins involved in muscle contraction, calcium-mediated signaling or cation transport (CASQ1, MLC2V, SLC25A4, MYL3). In PL pigs, we identified up-regulation of several genes that play crucial roles in reproduction: female gamete generation (BDP1, PTPN21, USP9X), fertilization (EGFR) and embryonic development (CPEB4). In the Pietrain breed, only seven genes were over-expressed (CISH, SPP1, TUBA8, ATP6V1C2, IGKC, predicted LOC100510960 and LOC100626400), and they play important roles in, for example, negative regulation of apoptosis, immune response, cell-cell signaling, cell growth and migration as well as the metabolic process. The functions of the majority of selected genes were consistent with phenotypic variation in investigated breeds; thus, we proposed a new panel of candidate genes that can be associated with economically important pig traits.

Analysis of the transcriptome in molecular epidemiology studies

The human transcriptome is complex, comprising multiple transcript types, mostly in the form of non-coding RNA (ncRNA). The majority of ncRNA is of the long form (lncRNA, ≥ 200 bp), which plays an important role in gene regulation through multiple mechanisms including epigenetics, chromatin modification, control of transcription factor binding, and regulation of alternative splicing. Both mRNA and ncRNA exhibit additional variability in the form of alternative splicing and RNA editing. All aspects of the human transcriptome can potentially be dysregulated by environmental exposures. Next-generation RNA sequencing (RNA-Seq) is the best available methodology to measure this although it has limitations, including experimental bias. The third phase of the MicroArray Quality Control Consortium project (MAQC-III), also called Sequencing Quality Control (SeQC), aims to address these limitations through standardization of experimental and bioinformatic methodologies. A limited number of toxicogenomic studies have been conducted to date using RNA-Seq. This review describes the complexity of the human transcriptome, the application of transcriptomics by RNA-Seq or microarray in molecular epidemiology studies, and limitations of these approaches including the type of cell or tissue analyzed, experimental variation, and confounding. By using good study designs with precise, individual exposure measurements, sufficient power and incorporation of phenotypic anchors, studies in human populations can identify biomarkers of exposure and/or early effect and elucidate mechanisms of action underlying associated diseases, even at low doses. Analysis of datasets at the pathway level can compensate for some of the limitations of RNA-Seq and, as more datasets become available, will increasingly elucidate the exposure-disease continuum.