Comparative RNA-Seq transcriptome analyses reveal dynamic time-dependent effects of 56Fe, 16O, and 28Si irradiation on the induction of murine hepatocellular carcinoma.

BACKGROUND: One of the health risks posed to astronauts during deep space flights is exposure to high charge, high-energy (HZE) ions (Z > 13), which can lead to the induction of hepatocellular carcinoma (HCC). However, little is known on the molecular mechanisms of HZE irradiation-induced HCC. RESULTS: We performed comparative RNA-Seq transcriptomic analyses to assess the carcinogenic effects of 600 MeV/n 56Fe (0.2 Gy), 1 GeV/n 16O (0.2 Gy), and 350 MeV/n 28Si (0.2 Gy) ions in a mouse model for irradiation-induced HCC. C3H/HeNCrl mice were subjected to total body irradiation to simulate space environment HZE-irradiation, and liver tissues were extracted at five different time points post-irradiation to investigate the time-dependent carcinogenic response at the transcriptomic level. Our data demonstrated a clear difference in the biological effects of these HZE ions, particularly immunological, such as Acute Phase Response Signaling, B Cell Receptor Signaling, IL-8 Signaling, and ROS Production in Macrophages. Also seen in this study were novel unannotated transcripts that were significantly affected by HZE. To investigate the biological functions of these novel transcripts, we used a machine learning technique known as self-organizing maps (SOMs) to characterize the transcriptome expression profiles of 60 samples (45 HZE-irradiated, 15 non-irradiated control) from liver tissues. A handful of localized modules in the maps emerged as groups of co-regulated and co-expressed transcripts. The functional context of these modules was discovered using overrepresentation analysis. We found that these spots typically contained enriched populations of transcripts related to specific immunological molecular processes (e.g., Acute Phase Response Signaling, B Cell Receptor Signaling, IL-3 Signaling), and RNA Transcription/Expression. CONCLUSIONS: A large number of transcripts were found differentially expressed post-HZE irradiation. These results provide valuable information for uncovering the differences in molecular mechanisms underlying HZE specific induced HCC carcinogenesis. Additionally, a handful of novel differentially expressed unannotated transcripts were discovered for each HZE ion. Taken together, these findings may provide a better understanding of biological mechanisms underlying risks for HCC after HZE irradiation and may also have important implications for the discovery of potential countermeasures against and identification of biomarkers for HZE-induced HCC.

A transcriptome atlas of silkworm silk glands revealed by PacBio single-molecule long-read sequencing.

The silk gland of the silkworm Bombyx mori is a specialized organ where silk proteins are efficiently synthesized under precise regulation that largely determines the properties of silk fibers. To understand the genes involved in the regulation of silk protein synthesis, considerable research has focused on the transcripts expressed in silk glands; however, the complete transcriptome profile of this organ has yet to be elucidated. Here, we report a full-length silk gland transcriptome obtained by PacBio single-molecule long-read sequencing technology. In total, 11,697 non-redundant transcripts were identified in mixed samples of silk glands dissected from larvae at five developmental stages. When compared with the published reference, the full-length transcripts optimized the structures of 3002 known genes, and a total of 9061 novel transcripts with an average length of 2171 bp were detected. Among these, 1403 (15.5%) novel transcripts were computationally revealed to be lncRNAs, 8135 (89.8%) novel transcripts were annotated to different protein and nucleotide databases, and 5655 (62.4%) novel transcripts were predicted to have complete ORFs. Furthermore, we found 1867 alternative splicing events, 2529 alternative polyadenylation events, 784 fusion events and 6596 SSRs. This study provides a comprehensive set of reference transcripts and greatly revises and expands the available silkworm transcript data. In addition, these data will be very useful for studying the regulatory mechanisms of silk protein synthesis.

A pan-transcriptome analysis shows that disease resistance genes have undergone more selection pressure during barley domestication.

BACKGROUND: It has become clear in recent years that many genes in a given species may not be found in a single genotype thus using sequences from a single genotype as reference may not be adequate for various applications. RESULTS: In this study we constructed a pan-transcriptome for barley by de novo assembling 288 sets of RNA-seq data from 32 cultivated barley genotypes and 31 wild barley genotypes. The pan-transcriptome consists of 756,632 transcripts with an average N50 length of 1240 bp. Of these, 289,697 (38.2%) were not found in the genome of the international reference genotype Morex. The novel transcripts are enriched with genes associated with responses to different stresses and stimuli. At the pan-transcriptome level, genotypes of wild barley have a higher proportion of disease resistance genes than cultivated ones. CONCLUSIONS: We demonstrate that the use of the pan-transcriptome dramatically improved the efficiency in detecting variation in barley. Analysing the pan-transcriptome also found that, compared with those in other categories, disease resistance genes have gone through stronger selective pressures during domestication.

Resolving the problem of multiple accessions of the same transcript deposited across various public databases.

Maintaining the consistency of genomic annotations is an increasingly complex task because of the iterative and dynamic nature of assembly and annotation, growing numbers of biological databases and insufficient integration of annotations across databases. As information exchange among databases is poor, a 'novel' sequence from one reference annotation could be annotated in another. Furthermore, relationships to nearby or overlapping annotated transcripts are even more complicated when using different genome assemblies. To better understand these problems, we surveyed current and previous versions of genomic assemblies and annotations across a number of public databases containing long noncoding RNA. We identified numerous discrepancies of transcripts regarding their genomic locations, transcript lengths and identifiers. Further investigation showed that the positional differences between reference annotations of essentially the same transcript could lead to differences in its measured expression at the RNA level. To aid in resolving these problems, we present the algorithm 'Universal Genomic Accession Hash (UGAHash)' and created an open source web tool to encourage the usage of the UGAHash algorithm. The UGAHash web tool (http://ugahash.uni-frankfurt.de) can be accessed freely without registration. The web tool allows researchers to generate Universal Genomic Accessions for genomic features or to explore annotations deposited in the public databases of the past and present versions. We anticipate that the UGAHash web tool will be a valuable tool to check for the existence of transcripts before judging the newly discovered transcripts as novel.

The identification and characterization of novel transcripts from RNA-seq data.

Owing greatly to the advancement of next-generation sequencing (NGS), the amount of NGS data is increasing rapidly. Although there are many NGS applications, one of the most commonly used techniques 'RNA sequencing (RNA-seq)' is rapidly replacing microarray-based techniques in laboratories around the world. As more and more of such techniques are standardized, allowing technicians to perform these experiments with minimal hands-on time and reduced experimental/operator-dependent biases, the bottleneck of such techniques is clearly visible; that is, data analysis. Further complicating the matter, increasing evidence suggests most of the genome is transcribed into RNA; however, the majority of these RNAs are not translated into proteins. These RNAs that do not become proteins are called 'noncoding RNAs (ncRNAs)'. Although some time has passed since the discovery of ncRNAs, their annotations remain poor, making analysis of RNA-seq data challenging. Here, we examine the current limitations of RNA-seq analysis using case studies focused on the detection of novel transcripts and examination of their characteristics. Finally, we validate the presence of novel transcripts using biological experiments, showing novel transcripts can be accurately identified when a series of filters is applied. In conclusion, novel transcripts that are identified from RNA-seq must be examined carefully before proceeding to biological experiments.

Transcriptome analysis of Cd-treated switchgrass root revealed novel transcripts and the importance of HSF/HSP network in switchgrass Cd tolerance.

KEY MESSAGE: Transcriptome analysis of Cd-treated switchgrass roots not only revealed novel switchgrass transcripts and gene structures but also highlighted the indispensable role of HSF/HSP network in switchgrass Cd tolerance. Switchgrass (Panicum virgatum L.), a C4 perennial tall grass, can be used for revegetation of Cd-contaminated soil. In the present study, a comparative transcriptome analysis of Cd-treated switchgrass roots was conducted. The result revealed a total of 462 novel transcripts and refined gene structures of 2337 transcripts. KEGG pathway and Gene Ontology analyses of the differentially expressed genes (DEGs) suggested that activation of redox homeostasis and oxidation-related metabolic processes were the primary response to Cd stress in switchgrass roots. In particular, 21 out of 23 differentially expressed shock transcription factor genes (HSFs), and 22 out of 23 differentially expressed heat shock protein genes (HSPs) had increased expression levels after Cd treatment. Furthermore, over-expressing one HSP-encoding gene in Arabidopsis significantly improved plant Cd tolerance. The result highlighted the activation of the redox homeostasis and the involvement of the HSF/HSP network in re-establishing normal protein conformation and thus cellular homeostasis in switchgrass upon Cd stress. These DEGs, especially those of the HSF/HSP network, could be used as candidate genes for further functional studies toward improved plant Cd tolerance in switchgrass and related species.

In silico characterisation of novel rice transcripts differentially expressed in phosphorus dificient conditions suggests a role of these transcripts in multiple abiotic stresses.

Phosphorus deficiency adversely affects crop productivity. The mechanism of tolerance in plants is not well understood. The current study successfully annotated a set of highly significant (Log2 RPKM ≥3) nine novel sequences up-regulated in P deficient condition identified from a low P tolerant rice genotype. Sequence annotation identified two transcripts (Os01g37260 and Os02g11060) carrying known domains, F-box and WD, respectively. Multiple Expectation maximization for Motif Elicitation (MEME) revealed presence of conserved domains like D[LP][HY][CL]D[CM][DT]C[AP][DQ][IQ]C, [EH][DN]HN[HS] [ER][FY][EP]I[HN]H which might play a role in phosphorus deficiency tolerance. Analysis of the upstream regions indicated presence of stress responsive elements like E Box, ABRE, and MYBCORE suggesting regulation of the novel transcripts by DNA binding. Protein localization prediction tool suggests that these novel proteins might be targeted to nucleus, chloroplast and cell wall. Transcripts Os02g03640 and Os02g10250 revealed potential target sites for microRNA binding suggesting role of novel miRNAs in low phosphorus response. Our analysis suggests that an F-box protein, Os01g37260 (OSFBx14) might be a promising candidate gene playing a role in multiple abiotic stresses including P deficiency.

Identification of Novel mRNA Isoforms Associated with Acute Heat Stress Response Using RNA Sequencing Data in Sprague Dawley Rats.

The molecular mechanisms underlying heat stress tolerance in animals to high temperatures remain unclear. This study identified the differentially expressed mRNA isoforms which narrowed down the most reliable DEG markers and molecular pathways that underlie the mechanisms of thermoregulation. This experiment was performed on Sprague Dawley rats housed at 22 °C (control group; CT), and three acute heat-stressed groups housed at 42 °C for 30 min (H30), 60 min (H60), and 120 min (H120). Earlier, we demonstrated that acute heat stress increased the rectal temperature of rats, caused abnormal changes in the blood biochemical parameters, as well as induced dramatic changes in the expression levels of genes through epigenetics and post-transcriptional regulation. Transcriptomic analysis using RNA-Sequencing (RNA-Seq) data obtained previously from blood (CT and H120), liver (CT, H30, H60, and H120), and adrenal glands (CT, H30, H60, and H120) was performed. The differentially expressed mRNA isoforms (DEIs) were identified and annotated by the CLC Genomics Workbench. Biological process and metabolic pathway analyses were performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 225, 5764, and 4988 DEIs in the blood, liver, and adrenal glands were observed. Furthermore, the number of novel differentially expressed transcript lengths with annotated genes and novel differentially expressed transcript with non-annotated genes were 136 and 8 in blood, 3549 and 120 in the liver, as well as 3078 and 220 in adrenal glands, respectively. About 35 genes were involved in the heat stress response, out of which, Dnaja1, LOC680121, Chordc1, AABR07011951.1, Hsp90aa1, Hspa1b, Cdkn1a, Hmox1, Bag3, and Dnaja4 were commonly identified in the liver and adrenal glands, suggesting that these genes may regulate heat stress response through interactions between the liver and adrenal glands. In conclusion, this study would enhance our understanding of the complex underlying mechanisms of acute heat stress, and the identified mRNA isoforms and genes can be used as potential candidates for thermotolerance selection in mammals.

Identification of novel mRNA isoforms associated with meat tenderness using RNA sequencing data in beef cattle.

The Warner-Bratzler shear force (WBSF) and myofibrillar fragmentation index (MFI) are complementary methodologies used to measure beef tenderness. Longissimus thoracis samples from the 20 most extreme bulls (out of 80 bulls set) for WBSF (tender (n = 10) and tough (n = 10)) and MFI (high (n = 10) and low (n = 10)) traits were collected to perform transcriptomic analysis using RNA-Sequencing. All analysis were performed through CLC Genomics Workbench. A total of 39 and 27 transcripts for WBSF and MFI phenotypes were DE, respectively. The possible DE novel mRNA isoforms, for WBSF and MFI traits, are myosin encoders (e.g. MYL1 and MYL6). In addition, we identified potential mRNA isoforms related to genes affecting the speed fibers degradation during the meat aging process. The DE novel transcripts are transcripted by genes with biological functions related to oxidative process, energy production and striated muscle contraction. The results suggest that the identified mRNA isoforms could be used as potential candidate to select animals in order to improve meat tenderness.

Determination of genetic aberrations and novel transcripts involved in the pathogenesis of oligodendroglioma using array comparative genomic hybridization and next generation sequencing.

The aim of the present study was to determine the genetic aberrations and novel transcripts, particularly the fusion transcripts, involved in the pathogenesis of low-grade and anaplastic oligodendroglioma. In the present study, tissue samples were obtained from patients with oligodendroglioma and additionally from archived tissue samples from the Brain Tumor Tissue Bank of the Brain Tumor Foundation of Canada. Six samples were obtained, three of which were low-grade oligodendroglioma and the other three anaplastic oligodendroglioma. DNA and RNA were extracted from each tissue sample. The resulting genomic DNA was then hybridized using the Agilent CytoSure 4×180K oligonucleotide array. Human reference DNA and samples were labeled using Cy3 cytidine 5'-triphosphate (CTP) and Cy5 CTP, respectively, while human Cot-1 DNA was used to reduce non-specific binding. Microarray-based comparative genomic hybridization data was then analyzed for genetic aberrations using the Agilent Cytosure Interpret software v3.4.2. The total RNA isolated from each sample was mixed with oligo dT magnetic beads to enrich for poly(A) mRNA. cDNAs were then synthesized and subjected to end-repair, poly(A) addition and connected using sequencing adapters using the Illumina TruSeq RNA Sample Preparation kit. The fragments were then purified and selected as templates for polymerase chain reaction amplification. The final library was constructed with fragments between 350-450 base pairs and sequenced using deep transcriptome sequencing on an Illumina HiSeq 2500 sequencer. The array comparative genomic hybridization revealed numerous amplifications and deletions on several chromosomes in all samples. However, the most interesting result was from the next generation sequencing, where one anaplastic oligodendroglioma sample was demonstrated to have five novel fusion genes that may potentially serve a critical role in tumor pathogenesis and progression.

Rule-based integration of RNA-Seq analyses tools for identification of novel transcripts.

Recent evidences suggest that a substantial amount of genome is transcribed more than that was anticipated, giving rise to a large number of unknown or novel transcripts. Identification of novel transcripts can provide key insights into understanding important cellular functions as well as molecular mechanisms underlying complex diseases like cancer. RNA-Seq has emerged as a powerful tool to detect novel transcripts, which previous profiling techniques failed to identify. A number of tools are available for enabling identification of novel transcripts at different levels. Read mappers such as TopHat, MapSplice, and SOAPsplice predict novel junctions, which are the indicators of novel transcripts. Cufflinks assembles novel transcripts based on alignment information and Oases performs de novo construction of transcripts. A common limitation of all these tools is prediction of sizable number of spurious or false positive (FP) novel transcripts. An approach that integrates information from all above sources and simultaneously scrutinizes FPs to correctly identify authentic novel transcripts of high confidence is proposed.

In vitro transcriptome analysis of two Chinese isolates of Streptococcus suis serotype 2.

The Streptococcus suis serotype 2 (S. suis 2) isolates 05ZYH33 and 98HAH33 have caused severe human infections in China. Using a strand-specific RNA-seq analysis, we compared the in vitro transcriptomes of these two Chinese isolates with that of a reference strain (P1/7). In the 89K genomic island that is specific to these Chinese isolates, a toxin-antitoxin system showed relatively high levels of transcription among the S. suis. The known virulence factors with high transcriptional activity in these two highly-pathogenic strains are mainly involved in adhesion, biofilm formation, hemolysis and the synthesis and transport of the outer membrane protein. Furthermore, our analysis of novel transcripts identified over 50 protein-coding genes with one of them encoding a toxin protein. We also predicted over 30 small RNAs (sRNAs) in each strain, and most of them are involved in riboswitches. We found that six sRNA candidates that are related to bacterial virulence, including cspA and rli38, are specific to Chinese isolates. These results provide insight into the factors responsible for the difference in virulence among the different S. suis 2 isolates.

Thousands of Novel Transcripts Identified in Mouse Cerebrum, Testis, and ES Cells Based on ribo-minus RNA Sequencing.

The high-throughput next-generation sequencing technologies provide an excellent opportunity for the detection of less-abundance transcripts that may not be identifiable by previously available techniques. Here, we report a discovery of thousands of novel transcripts (mostly non-coding RNAs) that are expressed in mouse cerebrum, testis, and embryonic stem (ES) cells, through an in-depth analysis of rmRNA-seq data. These transcripts show significant associations with transcriptional start and elongation signals. At the upstream of these transcripts we observed significant enrichment of histone marks (histone H3 lysine 4 trimethylation, H3K4me3), RNAPII binding sites, and cap analysis of gene expression tags that mark transcriptional start sites. Along the length of these transcripts, we also observed enrichment of histone H3 lysine 36 trimethylation (H3K36me3). Moreover, these transcripts show strong purifying selection in their genomic loci, exonic sequences, and promoter regions, implying functional constraints on the evolution of these transcripts. These results define a collection of novel transcripts in the mouse genome and indicate their potential functions in the mouse tissues and cells.