Analysis of the mouse transcriptome based on functional annotation of 60,770 full-length cDNAs

The long noncoding RNA lncR492 inhibits neural differentiation of murine embryonic stem cells

RNA interference (RNAi) screens have been shown to be valuable to study embryonic stem cell (ESC) self-renewal and they have been successfully applied to identify coding as well as noncoding genes required for maintaining pluripotency. Here, we used an RNAi library targeting >640 long noncoding RNAs (lncRNA) to probe for their role in early cell differentiation. Utilizing a Sox1-GFP ESC reporter cell line, we identified the lncRNA lncR492 as lineage-specific inhibitor of neuroectodermal differentiation. Molecular characterization showed that lncR492 interacts with the mRNA binding protein HuR and facilitates its inhibitory function by activation of Wnt signaling. Thus, lncRNAs modulate the fate decision of pluripotent stem cells.

Technological Developments in lncRNA Biology

It is estimated that more than 90% of the mammalian genome is transcribed as non-coding RNAs. Recent evidences have established that these non-coding transcripts are not junk or just transcriptional noise, but they do serve important biological purpose. One of the rapidly expanding fields of this class of transcripts is the regulatory lncRNAs, which had been a major challenge in terms of their molecular functions and mechanisms of action. The emergence of high-throughput technologies and the development in various conventional approaches have led to the expansion of the lncRNA world. The combination of multidisciplinary approaches has proven to be essential to unravel the complexity of their regulatory networks and helped establish the importance of their existence. Here, we review the current methodologies available for discovering and investigating functions of long non-coding RNAs (lncRNAs) and focus on the powerful technological advancement available to specifically address their functional importance.

A novel brown adipocyte-enriched long non-coding RNA that is required for brown adipocyte differentiation and sufficient to drive thermogenic gene program in white adipocytes

The thermogenic activities of brown and beige adipocytes can be exploited to reduce energy surplus and counteract obesity. Recent RNA sequencing studies have uncovered a number of long noncoding RNAs (lncRNAs) uniquely expressed in white and brown adipose tissues (WAT and BAT), but whether and how these lncRNAs function in adipogenesis remain largely unknown. Here, we report the identification of a novel brown adipocyte-enriched LncRNA (AK079912), and its nuclear localization, function and regulation. The expression of AK079912 increases during brown preadipocyte differentiation and in response to cold-stimulated browning of white adipocytes. Knockdown of AK079912 inhibits brown preadipocyte differentiation, manifested by reductions in lipid accumulation and down-regulation of adipogenic and BAT-specific genes. Conversely, ectopic expression of AK079912 in white preadipocytes up-regulates the expression of genes involved in thermogenesis. Mechanistically, inhibition of AK079912 reduces mitochondrial copy number and protein levels of mitochondria electron transport chain (ETC) complexes, whereas AK079912 overexpression increases the levels of ETC proteins. Lastly, reporter and pharmacological assays identify Pparγ as an upstream regulator of AK079912. These results provide new insights into the function of non-coding RNAs in brown adipogenesis and regulating browning of white adipocytes.

Lentivirus-mediated silencing of HOTAIR lncRNA restores gefitinib sensitivity by activating Bax/Caspase-3 and suppressing TGF-α/EGFR signaling in lung adenocarcinoma

Secondary resistance is a major limitation in the efficacy of epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor treatment of lung cancer. Previous studies have shown that expression of the long non-coding RNA HOX transcript antisense RNA (HOTAIR) is upregulated in lung cancer, which is correlated with metastasis and poor prognosis. However, the precise role of HOTAIR and its effects on gefitinib resistance in human lung adenocarcinoma are not known. To address this issue, in the present study we established a gefitinib-resistant (R)PC-9 human lung adenocarcinoma cell line and examined cell viability with the 3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium assay. We found that gefitinib concentrations <10 µM inhibited the viability of PC-9 but not RPC-9 cells in a dose-dependent manner. Lentivirus-mediated HOTAIR RNA interference induced cell apoptosis and S-phase arrest, as determined by terminal deoxynucleotidyl transferase-mediated dUTP nick end labeling and flow cytometry. Consistent with these observations, HOTAIR suppression was associated with tumor shrinkage and restoration of gefitinib sensitivity in RPC-9 xenograft mice. Immunohistochemical analyses and western blot revealed that HOTAIR silencing resulted in the upregulation of B cell lymphoma 2-associated X protein (Bax), Caspase-3 and transforming growth factor α (TGF-α) and downregulation of EGFR and B cell lymphoma 2 (Bcl-2) levels. These results indicate that HOTAIR normally prevents the activation of Bax/Caspase-3 while inducing TGF-α/EGFR signaling. Thus, targeting HOTAIR may be a novel therapeutic strategy for treating gefitinib-resistant lung adenocarcinoma.

Noncoding RNAs: Master Regulators of Inflammatory Signaling

Inflammatory signaling underlies many diseases, from arthritis to cancer. Our understanding of inflammation has thus far been limited to the world of proteins, because we are only just beginning to understand the role that noncoding RNAs (ncRNA) might play. It is now clear that ncRNA do not constitute transcriptional 'noise' but instead harbor physiological functions in controlling signaling pathways. In this review, we cover the newly discovered mechanisms and functions of ncRNAs in the regulation of inflammatory signaling. We also describe advances in experimental techniques allowing this field of research to take root. These findings have opened new avenues for putative therapeutic intervention in inflammatory diseases, which may be seen translated into clinical outcomes in the future.

Long Non-Coding RNAs: Emerging and Versatile Regulators in Host-Virus Interactions

Long non-coding RNAs (lncRNAs) are a class of non-protein-coding RNA molecules, which are involved in various biological processes, including chromatin modification, cell differentiation, pre-mRNA transcription and splicing, protein translation, etc. During the last decade, increasing evidence has suggested the involvement of lncRNAs in both immune and antiviral responses as positive or negative regulators. The immunity-associated lncRNAs modulate diverse and multilayered immune checkpoints, including activation or repression of innate immune signaling components, such as interleukin (IL)-8, IL-10, retinoic acid inducible gene I, toll-like receptors 1, 3, and 8, and interferon (IFN) regulatory factor 7, transcriptional regulation of various IFN-stimulated genes, and initiation of the cell apoptosis pathways. Additionally, some virus-encoded lncRNAs facilitate viral replication through individually or synergistically inhibiting the host antiviral responses or regulating multiple steps of the virus life cycle. Moreover, some viruses are reported to hijack host-encoded lncRNAs to establish persistent infections. Based on these amazing discoveries, lncRNAs are an emerging hotspot in host–virus interactions. In this review, we summarized the current findings of the host- or virus-encoded lncRNAs and the underlying mechanisms, discussed their impacts on immune responses and viral replication, and highlighted their critical roles in host–virus interactions.

Monitoring transcription initiation activities in rat and dog

The promoter landscape of several non-human model organisms is far from complete. As a part of FANTOM5 data collection, we generated 13 profiles of transcription initiation activities in dog and rat aortic smooth muscle cells, mesenchymal stem cells and hepatocytes by employing CAGE (Cap Analysis of Gene Expression) technology combined with single molecule sequencing. Our analyses show that the CAGE profiles recapitulate known transcription start sites (TSSs) consistently, in addition to uncover novel TSSs. Our dataset can be thus used with high confidence to support gene annotation in dog and rat species. We identified 28,497 and 23,147 CAGE peaks, or promoter regions, for rat and dog respectively, and associated them to known genes. This approach could be seen as a standard method for improvement of existing gene models, as well as discovery of novel genes. Given that the FANTOM5 data collection includes dog and rat matched cell types in human and mouse as well, this data would also be useful for cross-species studies.

Natural Antisense Transcripts at the Interface between Host Genome and Mobile Genetic Elements

Non-coding RNAs are involved in epigenetic processes, playing a role in the regulation of gene expression at the transcriptional and post-transcriptional levels. A particular group of ncRNA are natural antisense transcripts (NATs); these are transcribed in the opposite direction to protein coding transcripts and are widespread in eukaryotes. Their abundance, evidence of phylogenetic conservation and an increasing number of well-characterized examples of antisense-mediated gene regulation are indicative of essential biological roles of NATs. There is evidence to suggest that they interfere with their corresponding sense transcript to elicit concordant and discordant regulation. The main mechanisms involved include transcriptional interference as well as dsRNA formation. Sense–antisense hybrid formation can trigger RNA interference, RNA editing or protein kinase R. However, the exact molecular mechanisms elicited by NATs in the context of these regulatory roles are currently poorly understood. Several examples confirm that ectopic expression of antisense transcripts trigger epigenetic silencing of the related sense transcript. Genomic approaches suggest that the antisense transcriptome carries a broader biological significance which goes beyond the physiological regulation of the directly related sense transcripts. Because NATs show evidence of conservation we speculate that they played a role in evolution, with early eukaryotes gaining selective advantage through the regulatory effects. With the surge of genome and transcriptome sequencing projects, there is promise of a more comprehensive understanding of the biological role of NATs and the regulatory mechanisms involved.

The functions and unique features of long intergenic non-coding RNA

Long intergenic non-coding RNA (lincRNA) genes have diverse features that distinguish them from mRNA-encoding genes and exercise functions such as remodelling chromatin and genome architecture, RNA stabilization and transcription regulation, including enhancer-associated activity. Some genes currently annotated as encoding lincRNAs include small open reading frames (smORFs) and encode functional peptides and thus may be more properly classified as coding RNAs. lincRNAs may broadly serve to fine-tune the expression of neighbouring genes with remarkable tissue specificity through a diversity of mechanisms, highlighting our rapidly evolving understanding of the non-coding genome.

PARI Regulates Stalled Replication Fork Processing To Maintain Genome Stability upon Replication Stress in Mice

DNA replication is frequently perturbed by intrinsic, as well as extrinsic, genotoxic stress. At damaged forks, DNA replication and repair activities require proper coordination to maintain genome integrity. We show here that PARI antirecombinase plays an essential role in modulating the initial response to replication stress in mice. PARI is functionally dormant at replisomes during normal replication, but upon replication stress, it enhances nascent-strand shortening that is regulated by RAD51 and MRE11. PARI then promotes double-strand break induction, followed by new origin firing instead of replication restart. Such PARI function is apparently obstructive to replication but is nonetheless physiologically required for chromosome stability in vivo and ex vivo Of note, Pari-deficient embryonic stem cells exhibit spontaneous chromosome instability, which is attenuated by differentiation induction, suggesting that pluripotent stem cells have a preferential requirement for PARI that acts against endogenous replication stress. PARI is a latent modulator of stalled fork processing, which is required for stable genome inheritance under both endogenous and exogenous replication stress in mice.

Integration of lncRNA-miRNA-mRNA reveals novel insights into oviposition regulation in honey bees

BACKGROUND

The honey bee (Apis mellifera) is a highly diverse species commonly used for honey production and pollination services. The oviposition of the honey bee queen affects the development and overall performance of the colony. To investigate the ovary activation and oviposition processes on a molecular level, a genome-wide analysis of lncRNAs, miRNAs and mRNA expression in the ovaries of the queens was performed to screen for differentially expressed coding and noncoding RNAs. Further analysis identified relevant candidate genes or RNAs.

RESULTS

The analysis of the RNA profiles in different oviposition phase of the queens revealed that 740 lncRNAs, 81 miRNAs and 5,481 mRNAs were differently expressed during the ovary activation; 88 lncRNAs, 13 miRNAs and 338 mRNAs were differently expressed during the oviposition inhibition process; and finally, 100 lncRNAs, four miRNAs and 497 mRNAs were differently expressed during the oviposition recovery process. In addition, functional annotation of differentially expressed RNAs revealed several pathways that are closely related to oviposition, including hippo, MAPK, notch, Wnt, mTOR, TGF-beta and FoxO signaling pathways. Furthermore, in the QTL region for ovary size, 73 differentially expressed genes and 14 differentially expressed lncRNAs were located, which are considered as candidate genes affecting ovary size and oviposition. Moreover, a core set of genes served as bridges among different miRNAs were identified through the integrated analysis of lncRNA-miRNA-mRNA network.

CONCLUSION

The observed dramatic expression changes of coding and noncoding RNAs suggest that they may play a critical role in honey bee queens' oviposition. The identified candidate genes for oviposition activation and regulation could serve as a resource for further studies of genetic markers of oviposition in honey bees.

Reciprocal regulation of chromatin state and architecture by HOTAIRM1 contributes to temporal collinear HOXA gene activation

Thousands of long non-coding RNAs (lncRNAs) have been identified in mammals, many of which represent important regulators of gene expression. However, the mechanisms used by lncRNAs to control transcription remain largely uncharacterized. Here, we report on HOTAIRM1, a promising lncRNA biomarker in leukemia and solid tumors. We find that HOTAIRM1 contributes to three-dimensional chromatin organization changes required for the temporal collinear activation of HOXA genes. We show that distinct HOTAIRM1 variants preferentially associate with either UTX/MLL or PRC2 complexes to modulate the levels of activating and silencing marks at the bivalent domain. HOTAIRM1 contributes to physical dissociation of chromatin loops at the cluster proximal end, which delays recruitment of the histone demethylase UTX and transcription of central HOXA genes. Interestingly, we find overall proximal HOXA gene activation without chromatin conformation changes by HOTAIRM1 in a different cell type. Our results reveal a previously unappreciated relationship between chromatin structure, architecture and lncRNA function.

Long and small noncoding RNAs during oocyte-to-embryo transition in mammals

Oocyte-to-embryo transition is a process during which an oocyte ovulates, is fertilized, and becomes a developing embryo. It involves the first major genome reprogramming event in life of an organism where gene expression, which gave rise to a differentiated oocyte, is remodeled in order to establish totipotency in blastomeres of an early embryo. This remodeling involves replacement of maternal RNAs with zygotic RNAs through maternal RNA degradation and zygotic genome activation. This review is focused on expression and function of long noncoding RNAs (lncRNAs) and small RNAs during oocyte-to-embryo transition in mammals. LncRNAs are an assorted rapidly evolving collection of RNAs, which have no apparent protein-coding capacity. Their biogenesis is similar to mRNAs including transcriptional control and post-transcriptional processing. Diverse molecular and biological roles were assigned to lncRNAs although most of them probably did not acquire a detectable biological role. Since some lncRNAs serve as precursors for small noncoding regulatory RNAs in RNA silencing pathways, both types of noncoding RNA are reviewed together.

Long non-coding RNAs in osteosarcoma

Long non-coding RNAs (lncRNAs) with more than 200 nuleotides, have been explored to participate in various cancer types including osteosarcoma (OS), which is the most common kind of primary bone tumors with high morbidity in infants and adolescents. These oncogenic or tumor suppressive lncRNAs regulate OS pathogenesis, such as cell growth, proliferation, invasion, migration, metastasis and cell apoptosis, serve as independent prognostic biomarkers or play a significant role in multidrug resistance (MDR) in OS cells. In this review, we attempt to dissect the participation of lncRNAs in pathogenesis of OS and their potential clinical values, and also provide an outlook for viable biomarkers and therapeutic targets in OS.

LncRNA TRERNA1 Function as an Enhancer of SNAI1 Promotes Gastric Cancer Metastasis by Regulating Epithelial-Mesenchymal Transition

Long noncoding RNA (lncRNA) has been implicated in cancer, but little is known about the role of lncRNAs as regulators of tumor metastasis. In the present study, we demonstrate that lncRNA TRERNA1 acts like an enhancer of SNAI1 to promote cell invasion and migration and to contribute to metastasis of gastric cancer (GC). TRERNA1 is significantly unregulated in GCs and GC cell lines. Increased TRERNA1 is positively correlated with lymph node metastasis of GCs. RNA immunoprecipitation (RIP) and chromatin immunoprecipitation (ChIP) assays revealed that TRERNA1 functions as a scaffold to recruit EZH2 to epigenetically silence epithelial-mesenchymal transition marker CDH1 by H3K27me3 of its promoter region. TRERNA1 knockdown markedly reduced GC cell migration, invasion, tumorigenicity, and metastasis. Depletion of TRERNA1 reduced cell metastasis of GCs in vivo. Taken together, our findings indicated that TRERNA1 serves as a critical effector in GC progression by regulating CDH1 at the transcription level. It is implied that TRERNA1/CDH1 is a new potential target for GC therapy.

Systematic Identification and Molecular Characteristics of Long Noncoding RNAs in Pig Tissues

Long noncoding RNAs (lncRNAs) are non-protein-coding RNAs that are involved in a variety of biological processes. The pig is an important farm animal and an ideal biomedical model. In this study, we performed a genome-wide scan for lncRNAs in multiple tissue types from pigs. A total of 118 million paired-end 90 nt clean reads were obtained via strand-specific RNA sequencing, 80.4% of which were aligned to the pig reference genome. We developed a stringent bioinformatics pipeline to identify 2,139 high-quality multiexonic lncRNAs. The characteristic analysis revealed that the novel lncRNAs showed relatively shorter transcript length, fewer exons, and lower expression levels in comparison with protein-coding genes (PCGs). The guanine-cytosine (GC) content of the protein-coding exons and introns was significantly higher than that of the lncRNAs. Moreover, the single nucleotide polymorphism (SNP) density of lncRNAs was significantly higher than that of PCGs. Conservation analysis revealed that most lncRNAs were evolutionarily conserved among pigs, humans, and mice, such as CUFF.253988.1, which shares homology with human long noncoding RNA MALAT1. The findings of our study significantly increase the number of known lncRNAs in pigs.

Genome-wide computational analysis of potential long noncoding RNA mediated DNA:DNA:RNA triplexes in the human genome

BACKGROUND

Only a handful of long noncoding RNAs have been functionally characterized. They are known to modulate regulation through interacting with other biomolecules in the cell: DNA, RNA and protein. Though there have been detailed investigations on lncRNA-miRNA and lncRNA-protein interactions, the interaction of lncRNAs with DNA have not been studied extensively. In the present study, we explore whether lncRNAs could modulate genomic regulation by interacting with DNA through the formation of highly stable DNA:DNA:RNA triplexes.

METHODS

We computationally screened 23,898 lncRNA transcripts as annotated by GENCODE, across the human genome for potential triplex forming sequence stretches (PTS). The PTS frequencies were compared across 5'UTR, CDS, 3'UTR, introns, promoter and 1000 bases downstream of the transcription termination sites. These regions were annotated by mapping to experimental regulatory regions, classes of repeat regions and transcription factors. We validated few putative triplex mediated interactions where lncRNA-gene pair interaction is via pyrimidine triplex motif using biophysical methods.

RESULTS

We identified 20,04,034 PTS sites to be enriched in promoter and intronic regions across human genome. Additional analysis of the association of PTS with core promoter elements revealed a systematic paucity of PTS in all regulatory regions, except TF binding sites. A total of 25 transcription factors were found to be associated with PTS. Using an interaction network, we showed that a subset of the triplex forming lncRNAs, have a positive association with gene promoters. We also demonstrated an in vitro interaction of one lncRNA candidate with its predicted gene target promoter regions.

CONCLUSIONS

Our analysis shows that PTS are enriched in gene promoter and largely associated with simple repeats. The current study suggests a major role of a subset of lncRNAs in mediating chromatin organization modulation through CTCF and NSRF proteins.

Long intergenic non-coding RNA GALMD3 in chicken Marek's disease

Long intergenic non-coding RNAs (lincRNAs) are transcribed from non-coding DNA sequences. Studies have revealed that aberrant expressions of lincRNAs are associated with various types of cancers and neurological disorders. Marek's disease (MD) is a highly contagious T-cell lymphoid neoplasia of chicken induced by Marek's disease virus (MDV). In this study, we first identified and validated linc-GALMD3 highly expressed in MDV-infected CD4+ T cells by RNA-Seq and qRT-PCR. By RNA-Seq analysis in MDCC-MSB1 cells after loss of function of linc-GALMD3 by shRNA, we found that linc-GALMD3 could positively cis-regulate its downstream gga-miR-223 gene expression. In contrast, it could trans-regulate the 748 differentially expressed genes (FDR < 0.01) that were mainly enriched into mitochondrial structure and cell cycle processes using GO analysis. Of these, the most significantly expressed gene EPYC might cause iris lesion in MD. The other eight genes, NDUFA4, NDUFB6, NDUFV1, NDUFS8, SDHB, UQCRC1, UQCRC2, and COX7A2, actively participated in oxidative phosphorylation in mitochondrial dysfunction and cell death. Most importantly, we found that the MDV replication was repressed when linc-GALMD3 was knocked down in CEF cells. Our results suggested that linc-GALMD3 might be a critical regulator in chicken MD and could be used as a candidate-promising mark for MD prevention, diagnosis, and treatment.

The FANTOM5 collection, a data series underpinning mammalian transcriptome atlases in diverse cell types

The latest project from the FANTOM consortium, an international collaborative effort initiated by RIKEN, generated atlases of transcriptomes, in particular promoters, transcribed enhancers, and long-noncoding RNAs, across a diverse set of mammalian cell types. Here, we introduce the FANTOM5 collection, bringing together data descriptors, articles and analyses of FANTOM5 data published across the Nature Research journals. Associated data are openly available for reuse by all.

FANTOM5 CAGE profiles of human and mouse samples

In the FANTOM5 project, transcription initiation events across the human and mouse genomes were mapped at a single base-pair resolution and their frequencies were monitored by CAGE (Cap Analysis of Gene Expression) coupled with single-molecule sequencing. Approximately three thousands of samples, consisting of a variety of primary cells, tissues, cell lines, and time series samples during cell activation and development, were subjected to a uniform pipeline of CAGE data production. The analysis pipeline started by measuring RNA extracts to assess their quality, and continued to CAGE library production by using a robotic or a manual workflow, single molecule sequencing, and computational processing to generate frequencies of transcription initiation. Resulting data represents the consequence of transcriptional regulation in each analyzed state of mammalian cells. Non-overlapping peaks over the CAGE profiles, approximately 200,000 and 150,000 peaks for the human and mouse genomes, were identified and annotated to provide precise location of known promoters as well as novel ones, and to quantify their activities.

LncRNA-TP53TG1 Participated in the Stress Response Under Glucose Deprivation in Glioma

Gliomas are the most common brain tumors of the center nervous system. And long non-coding RNAs (lncRNAs) are non-protein coding transcripts, which have been considered as one type of gene expression regulator for cancer development. In this study, we investigated the role of lncRNA-TP53TG1 in response to glucose deprivation in human gliomas. The expression levels of TP53TG1 in glioma tissues and cells were analyzed by qRT-PCR. In addition, the influence of TP53TG1 on glucose metabolism related genes at the mRNA level during both high and low glucose treatment was detected by qRT-PCR. MTT, clonogenicity assays, and flow cytometry were performed to detect the cell proliferation and cell apoptosis. Furthermore, the migration of glioma cells was examined by Transwell assays. The expression of TP53TG1 was significantly higher in human glioma tissues or cell lines compared with normal brain tissue or NHA. Moreover, TP53TG1 and some tumor glucose metabolism related genes, such as GRP78, LDHA, and IDH1 were up-regulated significantly in U87 and LN18 cells under glucose deprivation. In addition, knockdown of TP53TG1 decreased cell proliferation and migration and down-regulated GRP78 and IDH1 expression levels and up-regulated PKM2 levels in U87 cells under glucose deprivation. However, over-expression of TP53TG1 showed the opposite tendency. Moreover, the effects of TP53TG1 were more remarkable in low glucose than that in high glucose. Our data showed that TP53TG1 under glucose deprivation may promote cell proliferation and migration by influencing the expression of glucose metabolism related genes in glioma. J. Cell. Biochem. 118: 4897-4904, 2017. © 2017 Wiley Periodicals, Inc.

Association of long non-coding RNA H19 and microRNA-21 expression with the biological features and prognosis of non-small cell lung cancer

In recent years, long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) have been shown to play important roles in tumor biological function. The aim of this study was to investigate the diagnostic and prognostic value of lncRNA H19 and miR-21 expression in non-small-cell lung cancer (NSCLC). H19 and miR-21 expression was measured in tumor tissues and corresponding non-tumor lung tissues from 200 patients by quantitative reverse transcription polymerase chain reaction. Moreover, the in vitro and in vivo effects of H19 knock out in A549 cells were investigated. Expression of both H19 and miR-21 was significantly higher in lung tissues from patients with NSCLC than in normal lung tissues. Increased expression of H19 and miR-21 was positively correlated with advanced tumor-node-metastasis stage and tumor size. miR-21 expression was highest in stage I and II NSCLC, whereas H19 expression was highest in stage III and IV NSCLC. Knockout of H19 significantly inhibited NSCLC cell proliferation both in vitro and in vivo. The results show that H19 may mainly contributes to the progression of NSCLC, and its expression levels can reflect the invasive and metastatic status to some extent. miR-21 expression more likely plays an important role in early stage NSCLC. Moreover, H19 and miR-21 interact in the regulation of NSCLC, and with greater expression of both H19 and miR-21, overall survival decreased. The combination of H19 and miR-21 may have diagnostic value in NSCLC and represent a target for new NSCLC treatments.

A resource for functional profiling of noncoding RNA in the yeast Saccharomyces cerevisiae

Eukaryotic genomes are extensively transcribed, generating many different RNAs with no known function. We have constructed 1502 molecular barcoded ncRNA gene deletion strains encompassing 443 ncRNAs in the yeast Saccharomyces cerevisiae as tools for ncRNA functional analysis. This resource includes deletions of small nuclear RNAs (snRNAs), transfer RNAs (tRNAs), small nucleolar RNAs (snoRNAs), and other annotated ncRNAs as well as the more recently identified stable unannotated transcripts (SUTs) and cryptic unstable transcripts (CUTs) whose functions are largely unknown. Specifically, deletions have been constructed for ncRNAs found in the intergenic regions, not overlapping genes or their promoters (i.e., at least 200 bp minimum distance from the closest gene start codon). The deletion strains carry molecular barcodes designed to be complementary with the protein gene deletion collection enabling parallel analysis experiments. These strains will be useful for the numerous genomic and molecular techniques that utilize deletion strains, including genome-wide phenotypic screens under different growth conditions, pooled chemogenomic screens with drugs or chemicals, synthetic genetic array analysis to uncover novel genetic interactions, and synthetic dosage lethality screens to analyze gene dosage. Overall, we created a valuable resource for the RNA community and for future ncRNA research.

Gene Evolution: Getting Something from Nothing

New genes arise from pre-existing genes, but some de novo origin from non-genic sequence also seems plausible. A new study has surprisingly concluded that 25% of random DNA sequences yield beneficial products when expressed in bacteria.

A Comprehensive Mouse Transcriptomic BodyMap across 17 Tissues by RNA-seq

The mouse has been widely used as a model organism for studying human diseases and for evaluating drug safety and efficacy. Many diseases and drug effects exhibit tissue specificity that may be reflected by tissue-specific gene-expression profiles. Here we construct a comprehensive mouse transcriptomic BodyMap across 17 tissues of six-weeks old C57BL/6JJcl mice using RNA-seq. We find different expression patterns between protein-coding and non-coding genes. Liver expressed the least complex transcriptomes, that is, the smallest number of genes detected in liver across all 17 tissues, whereas testis and ovary harbor more complex transcriptomes than other tissues. We report a comprehensive list of tissue-specific genes across 17 tissues, along with a list of 4,781 housekeeping genes in mouse. In addition, we propose a list of 27 consistently and highly expressed genes that can be used as reference controls in expression-profiling analysis. Our study provides a unique resource of mouse gene-expression profiles, which is helpful for further biomedical research.

HOTAIR: a key regulator in gynecologic cancers

Long non-coding RNAs (lncRNAs) play critical roles in the initiation and progression of human cancers. HOX transcript antisense RNA (HOTAIR) is an lncRNA localized to the mammalian HOXC gene cluster; it can interact with polycomb repressive complex 2 and the lysine-specific histone demethylase/CoREST/REST complex, and it manipulates the expression of various genes. HOTAIR promotes tumor invasion and metastasis by silencing tumor suppressors, and activating oncogenes and signaling pathways. HOTAIR is deregulated in many human cancers; despite its critical roles in health and disease, the underlying mechanisms governing HOTAIR function are unknown. In this review, we summarize the recent findings on the roles of HOTAIR in gynecologic cancers.

Relatively frequent switching of transcription start sites during cerebellar development

BACKGROUND

Alternative transcription start site (TSS) usage plays important roles in transcriptional control of mammalian gene expression. The growing interest in alternative TSSs and their role in genome diversification spawned many single-gene studies on differential usages of tissue-specific or temporal-specific alternative TSSs. However, exploration of the switching usage of alternative TSS usage on a genomic level, especially in the central nervous system, is largely lacking.

RESULTS

In this study, We have prepared a unique set of time-course data for the developing cerebellum, as part of the FANTOM5 consortium ( http://fantom.gsc.riken.jp/5/ ) that uses their innovative capturing of 5' ends of all transcripts followed by Helicos next generation sequencing. We analyzed the usage of all transcription start sites (TSSs) at each time point during cerebellar development that provided information on multiple RNA isoforms that emerged from the same gene. We developed a mathematical method that systematically compares the expression of different TSSs of a gene to identify temporal crossover and non-crossover switching events. We identified 48,489 novel TSS switching events in 5433 genes during cerebellar development. This includes 9767 crossover TSS switching events in 1511 genes, where the dominant TSS shifts over time.

CONCLUSIONS

We observed a relatively high prevalence of TSS switching in cerebellar development where the resulting temporally-specific gene transcripts and protein products can play important regulatory and functional roles.

Long Noncoding RNAs as Biomarkers in Cancer

Long noncoding RNAs (lncRNAs) are a relatively well-characterized class of noncoding RNA (ncRNA) molecules, involved in the regulation of various cell processes, including transcription, intracellular trafficking, and chromosome remodeling. Their deregulation has been associated with the development and progression of various cancer types, the fact which makes them suitable as biomarkers for cancer diagnosis and prognosis. In recent years, detection of cancer-associated lncRNAs in body fluids of cancer patients has proven itself as an especially valuable method to effectively diagnose cancer. Cancer diagnosis and prognosis employing circulating lncRNAs are preferential when compared to classical biopsies of tumor tissues, especially due to their noninvasiveness, and have great potential for routine usage in clinical practice. Thus, this review focuses on summarizing the perspectives of lncRNAs as biomarkers in cancer, based on evaluating their expression profiles determined in body fluids of cancer patients.

Identification and Characterization of a Class of MALAT1-like Genomic Loci

The MALAT1 (Metastasis-Associated Lung Adenocarcinoma Transcript 1) gene encodes a noncoding RNA that is processed into a long nuclear retained transcript (MALAT1) and a small cytoplasmic tRNA-like transcript (mascRNA). Using an RNA sequence- and structure-based covariance model, we identified more than 130 genomic loci in vertebrate genomes containing the MALAT1 3' end triple-helix structure and its immediate downstream tRNA-like structure, including 44 in the green lizard Anolis carolinensis. Structural and computational analyses revealed a co-occurrence of components of the 3' end module. MALAT1-like genes in Anolis carolinensis are highly expressed in adult testis, thus we named them testis-abundant long noncoding RNAs (tancRNAs). MALAT1-like loci also produce multiple small RNA species, including PIWI-interacting RNAs (piRNAs), from the antisense strand. The 3' ends of tancRNAs serve as potential targets for the PIWI-piRNA complex. Thus, we have identified an evolutionarily conserved class of long noncoding RNAs (lncRNAs) with similar structural constraints, post-transcriptional processing, and subcellular localization and a distinct function in spermatocytes.

The Dimensions, Dynamics, and Relevance of the Mammalian Noncoding Transcriptome

The combination of pervasive transcription and prolific alternative splicing produces a mammalian transcriptome of great breadth and diversity. The majority of transcribed genomic bases are intronic, antisense, or intergenic to protein-coding genes, yielding a plethora of short and long non-protein-coding regulatory RNAs. Long noncoding RNAs (lncRNAs) share most aspects of their biogenesis, processing, and regulation with mRNAs. However, lncRNAs are typically expressed in more restricted patterns, frequently from enhancers, and exhibit almost universal alternative splicing. These features are consistent with their role as modular epigenetic regulators. We describe here the key studies and technological advances that have shaped our understanding of the dimensions, dynamics, and biological relevance of the mammalian noncoding transcriptome.

Ectopically expressed Slc34a2a sense-antisense transcripts cause a cerebellar phenotype in zebrafish embryos depending on RNA complementarity and Dicer

Natural antisense transcripts (NATs) are complementary to protein coding genes and potentially regulate their expression. Despite widespread occurrence of NATs in the genomes of higher eukaryotes, their biological role and mechanism of action is poorly understood. Zebrafish embryos offer a unique model system to study sense-antisense transcript interplay at whole organism level. Here, we investigate putative antisense transcript-mediated mechanisms by ectopically co-expressing the complementary transcripts during early zebrafish development. In zebrafish the gene Slc34a2a (Na-phosphate transporter) is bi-directionally transcribed, the NAT predominantly during early development up to 48 hours after fertilization. Declining levels of the NAT, Slc34a2a(as), coincide with an increase of the sense transcript. At that time, sense and antisense transcripts co-localize in the endoderm at near equal amounts. Ectopic expression of the sense transcript during embryogenesis leads to specific failure to develop a cerebellum. The defect is RNA-mediated and dependent on sense-antisense complementarity. Overexpression of a Slc34a2a paralogue (Slc34a2b) or the NAT itself had no phenotypic consequences. Knockdown of Dicer rescued the brain defect suggesting that RNA interference is required to mediate the phenotype. Our results corroborate previous reports of Slc34a2a-related endo-siRNAs in two days old zebrafish embryos and emphasize the importance of coordinated expression of sense-antisense transcripts. Our findings suggest that RNAi is involved in gene regulation by certain natural antisense RNAs.

Cilia- and Flagella-Associated Protein 69 Regulates Olfactory Transduction Kinetics in Mice

Animals detect odorous chemicals through specialized olfactory sensory neurons (OSNs) that transduce odorants into neural electrical signals. We identified a novel and evolutionarily conserved protein, cilia- and flagella-associated protein 69 (CFAP69), in mice that regulates olfactory transduction kinetics. In the olfactory epithelium, CFAP69 is enriched in OSN cilia, where olfactory transduction occurs. Bioinformatic analysis suggests that a large portion of CFAP69 can form Armadillo-type α-helical repeats, which may mediate protein-protein interactions. OSNs lacking CFAP69, remarkably, displayed faster kinetics in both the on and off phases of electrophysiological responses at both the neuronal ensemble level as observed by electroolfactogram and the single-cell level as observed by single-cell suction pipette recordings. In single-cell analysis, OSNs lacking CFAP69 showed faster response integration and were able to fire APs more faithfully to repeated odor stimuli. Furthermore, both male and female mutant mice that specifically lack CFAP69 in OSNs exhibited attenuated performance in a buried food pellet test when a background of the same odor to the food pellet was present even though they should have better temporal resolution of coding olfactory stimulation at the peripheral. Therefore, the role of CFAP69 in the olfactory system seems to be to allow the olfactory transduction machinery to work at a precisely regulated range of response kinetics for robust olfactory behavior.SIGNIFICANCE STATEMENT Sensory receptor cells are generally thought to evolve to respond to sensory cues as fast as they can. This idea is consistent with mutational analyses in various sensory systems, where mutations of sensory receptor cells often resulted in reduced response size and slowed response kinetics. Contrary to this idea, we have found that there is a kinetic "damper" present in the olfactory transduction cascade of the mouse that slows down the response kinetics and, by doing so, it reduces the peripheral temporal resolution in coding odor stimuli and allows for robust olfactory behavior. This study should trigger a rethinking of the significance of the intrinsic speed of sensory transduction and the pattern of the peripheral coding of sensory stimuli.

The Novel Landscape of Long Non-Coding RNAs in Response to Human Foamy Virus Infection Characterized by RNA-Seq

Human foamy virus (HFV) is a complex and unique retrovirus with the longest genomes among retroviruses that are used as vectors for gene therapy. Long non-coding RNAs (lncRNAs) are regarded as key regulators that are involved in diverse biological processes during viral infection. However, the role of lncRNAs in HFV infection remains unknown. In this study, we utilized next-generation sequencing to first characterize lncRNAs in 293T cells after HFV infection, evaluating length distribution, exon number distribution, volcano picture, and lncRNA class distribution. We identified 11,336 lncRNAs (4,729 upregulated lncRNAs and 6,588 downregulated lncRNAs) and 61,367 mRNAs (30,133 upregulated mRNAs and 31,220 downregulated mRNAs), which were differentially expressed in the HFV-infected 293T cells. Subsequently, six differentially expressed lncRNAs characterized using RNA-seq were confirmed by quantitative real-time polymerase chain reaction assays. Interestingly, Gene Ontology (GO)/Gene Ontology Tree Machine (GOTM) and Kyoto Encyclopedia of Gene and Genomes (KEGG) pathway analyses indicated that positive regulation of interleukin 8 (IL8) production and cytokine-cytokine receptor interaction might be involved in the functional enrichment of lncRNAs. Moreover, cis-acting and trans-acting regulatory networks show that NR_028036 may target the fas gene in a cis-acting manner and that ENST00000354838 may target the IL18 gene in a trans-acting manner. Overall, these results not only provide novel insights into the relationship between HFV and lncRNAs in the host response to infection but also have implications for the future wider application of HFV as a vector.

LncRNA-Six1 Encodes a Micropeptide to Activate Six1 in Cis and Is Involved in Cell Proliferation and Muscle Growth

Long non-coding RNAs (lncRNAs) play important roles in epigenetic regulation of skeletal muscle development. In our previous RNA-seq study (accession number GSE58755), we found that lncRNA-Six1 is an lncRNA that is differentially expressed between White Recessive Rock (WRR) and Xinghua (XH) chicken. In this study, we have further demonstrated that lncRNA-Six1 is located 432 bp upstream of the gene encoding the protein Six homeobox 1 (Six1). A dual-luciferase reporter assay identified that lncRNA-Six1 overlaps the Six1 proximal promoter. In lncRNA-Six1, a micropeptide of about 7.26 kDa was found to play an important role in the lncRNA-Six1 in cis activity. Overexpression of lncRNA-Six1 promoted the mRNA and protein expression level of the Six1 gene, while knockdown of lncRNA-Six1 inhibited Six1 expression. Moreover, tissue expression profiles showed that both the lncRNA-Six1 and the Six1 mRNA were highly expressed in chicken breast tissue. LncRNA-Six1 overexpression promoted cell proliferation and induced cell division. Conversely, its loss of function inhibited cell proliferation and reduced cell viability. Similar effects were observed after overexpression or knockdown of the Six1 gene. In addition, overexpression or knockdown of Six1 promoted or inhibited, respectively, the expression levels of muscle-growth-related genes, such as MYOG, MYHC, MYOD, IGF1R, and INSR. Taken together, these data demonstrate that lncRNA-Six1 carries out cis-acting regulation of the protein-encoding Six1 gene, and encodes a micropeptide to activate Six1 gene, thus promoting cell proliferation and being involved in muscle growth.

Non-nuclear Pool of Splicing Factor SFPQ Regulates Axonal Transcripts Required for Normal Motor Development

Recent progress revealed the complexity of RNA processing and its association to human disorders. Here, we unveil a new facet of this complexity. Complete loss of function of the ubiquitous splicing factor SFPQ affects zebrafish motoneuron differentiation cell autonomously. In addition to its nuclear localization, the protein unexpectedly localizes to motor axons. The cytosolic version of SFPQ abolishes motor axonal defects, rescuing key transcripts, and restores motility in the paralyzed sfpq null mutants, indicating a non-nuclear processing role in motor axons. Novel variants affecting the conserved coiled-coil domain, so far exclusively found in fALS exomes, specifically affect the ability of SFPQ to localize in axons. They broadly rescue morphology and motility in the zebrafish mutant, but alter motor axon morphology, demonstrating functional requirement for axonal SFPQ. Altogether, we uncover the axonal function of the splicing factor SFPQ in motor development and highlight the importance of the coiled-coil domain in this process.

Video Abstract

•

SFPQ splicing factor is present in motor axons

•

Non-nuclear SFPQ is able to drive axon maturation and connectivity

•

Loss of axonal SFPQ affects axonal morphology

•

Coiled-coil domain of the protein is important for non-nuclear localization

Defining Trends in Global Gene Expression in Arabian Horses with Cerebellar Abiotrophy

Equine cerebellar abiotrophy (CA) is a hereditary neurodegenerative disease that affects the Purkinje neurons of the cerebellum and causes ataxia in Arabian foals. Signs of CA are typically first recognized either at birth to any time up to 6 months of age. CA is inherited as an autosomal recessive trait and is associated with a single nucleotide polymorphism (SNP) on equine chromosome 2 (13074277G>A), located in the fourth exon of TOE1 and in proximity to MUTYH on the antisense strand. We hypothesize that unraveling the functional consequences of the CA SNP using RNA-seq will elucidate the molecular pathways underlying the CA phenotype. RNA-seq (100 bp PE strand-specific) was performed in cerebellar tissue from four CA-affected and five age-matched unaffected horses. Three pipelines for differential gene expression (DE) analysis were used (Tophat2/Cuffdiff2, Kallisto/EdgeR, and Kallisto/Sleuth) with 151 significant DE genes identified by all three pipelines in CA-affected horses. TOE1 (Log2(foldchange) = 0.92, p = 0.66) and MUTYH (Log2(foldchange) = 1.13, p = 0.66) were not differentially expressed. Among the major pathways that were differentially expressed, genes associated with calcium homeostasis and specifically expressed in Purkinje neurons, CALB1 (Log2(foldchange) = -1.7, p < 0.01) and CA8 (Log2(foldchange) = -0.97, p < 0.01), were significantly down-regulated, confirming loss of Purkinje neurons. There was also a significant up-regulation of markers for microglial phagocytosis, TYROBP (Log2(foldchange) = 1.99, p < 0.01) and TREM2 (Log2(foldchange) = 2.02, p < 0.01). These findings reaffirm a loss of Purkinje neurons in CA-affected horses along with a potential secondary loss of granular neurons and activation of microglial cells.

Potentials of Long Noncoding RNAs (LncRNAs) in Sarcoma: From Biomarkers to Therapeutic Targets

Sarcoma includes some of the most heterogeneous tumors, which make the diagnosis, prognosis and treatment of these rare yet diverse neoplasms especially challenging. Long noncoding RNAs (lncRNAs) are important regulators of cancer initiation and progression, which implies their potential as neoteric prognostic and diagnostic markers in cancer, including sarcoma. A relationship between lncRNAs and sarcoma pathogenesis and progression is emerging. Recent studies demonstrate that lncRNAs influence sarcoma cell proliferation, metastasis, and drug resistance. Additionally, lncRNA expression profiles are predictive of sarcoma prognosis. In this review, we summarize contemporary advances in the research of lncRNA biogenesis and functions in sarcoma. We also highlight the potential for lncRNAs to become innovative diagnostic and prognostic biomarkers as well as therapeutic targets in sarcoma.

Genome-wide identification of conserved intronic non-coding sequences using a Bayesian segmentation approach

Background

Computational identification of non-coding RNAs (ncRNAs) is a challenging problem. We describe a genome-wide analysis using Bayesian segmentation to identify intronic elements highly conserved between three evolutionarily distant vertebrate species: human, mouse and zebrafish. We investigate the extent to which these elements include ncRNAs (or conserved domains of ncRNAs) and regulatory sequences.

Results

We identified 655 deeply conserved intronic sequences in a genome-wide analysis. We also performed a pathway-focussed analysis on genes involved in muscle development, detecting 27 intronic elements, of which 22 were not detected in the genome-wide analysis. At least 87% of the genome-wide and 70% of the pathway-focussed elements have existing annotations indicative of conserved RNA secondary structure. The expression of 26 of the pathway-focused elements was examined using RT-PCR, providing confirmation that they include expressed ncRNAs. Consistent with previous studies, these elements are significantly over-represented in the introns of transcription factors.

Conclusions

This study demonstrates a novel, highly effective, Bayesian approach to identifying conserved non-coding sequences. Our results complement previous findings that these sequences are enriched in transcription factors. However, in contrast to previous studies which suggest the majority of conserved sequences are regulatory factor binding sites, the majority of conserved sequences identified using our approach contain evidence of conserved RNA secondary structures, and our laboratory results suggest most are expressed.

Functional roles at DNA and RNA levels are not mutually exclusive, and many of our elements possess evidence of both. Moreover, ncRNAs play roles in transcriptional and post-transcriptional regulation, and this may contribute to the over-representation of these elements in introns of transcription factors. We attribute the higher sensitivity of the pathway-focussed analysis compared to the genome-wide analysis to improved alignment quality, suggesting that enhanced genomic alignments may reveal many more conserved intronic sequences.

Electronic supplementary material

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

Epigenetic Regulations in Diabetic Nephropathy

Diabetic nephropathy (DN) is a chronic complication of diabetes and the most common cause of end-stage kidney disease. It has been reported that multiple factors are involved in the pathogenesis of DN, while the molecular mechanisms that lead to DN are still not fully understood. Numerous risk factors for the development of diabetic nephropathy have been proposed, including ethnicity and inherited genetic differences. Recently, with the development of high-throughput technologies, there is emerging evidence that suggests the important role of epigenetic mechanisms in the pathogenesis of DN. Epigenetic regulations, including DNA methylation, noncoding RNAs, and histone modifications, play a pivotal role in DN pathogenesis by a second layer of gene regulation. All these findings can contribute to developing novel therapies for DN.

Human Pancreatic β Cell lncRNAs Control Cell-Specific Regulatory Networks

Recent studies have uncovered thousands of long non-coding RNAs (lncRNAs) in human pancreatic β cells. β cell lncRNAs are often cell type specific and exhibit dynamic regulation during differentiation or upon changing glucose concentrations. Although these features hint at a role of lncRNAs in β cell gene regulation and diabetes, the function of β cell lncRNAs remains largely unknown. In this study, we investigated the function of β cell-specific lncRNAs and transcription factors using transcript knockdowns and co-expression network analysis. This revealed lncRNAs that function in concert with transcription factors to regulate β cell-specific transcriptional networks. We further demonstrate that the lncRNA PLUTO affects local 3D chromatin structure and transcription of PDX1, encoding a key β cell transcription factor, and that both PLUTO and PDX1 are downregulated in islets from donors with type 2 diabetes or impaired glucose tolerance. These results implicate lncRNAs in the regulation of β cell-specific transcription factor networks.

Rpph1 Upregulates CDC42 Expression and Promotes Hippocampal Neuron Dendritic Spine Formation by Competing with miR-330-5p

Alzheimer's disease (AD) is a heterogeneous neurodegenerative disease. Recent studies employing microRNA-seq and genome-wide sequencing have identified some non-coding RNAs that are influentially involved in AD pathogenesis. Non-coding RNAs can compete with other endogenous RNAs by microRNA response elements (MREs) and manipulate biological processes, such as tumorigenesis. However, only a few non-coding RNAs have been reported in the pathogenesis of AD. In this study, we constructed the first competing endogenous RNA (ceRNA) network leveraging whole transcriptome sequencing and a previously studied microRNA-seq of APPswe/PS1ΔE9 transgenic mice. The underlying mechanisms for the involvement of ceRNA in AD were validated using the Dual Luciferase Reporter Assay, detection of transcription levels by quantitative RT-PCR and translation levels by Western blotting, and morphological examination in primary cultured neurons. In the ceRNA network, four lncRNAs (C030034L19Rik, Rpph1, A830012C17Rik, and Gm15477) and five miRNAs (miR-182-5p, miR-330-5p, miR-326-3p, miR-132-3p, and miR-484) are enriched in nine pathways and an AD-related gene pool. Among them, Ribonuclease P RNA component H1 (Rpph1) is upregulated in the cortex of APPswe/PS1ΔE9 mice compared to wild type controls. Rpph1 binds to miR326-3p/miR-330-5p and causes the release of their downstream target Cdc42, which leads to CDC42 upregulation. This effect was disrupted upon mutation of the MRE on Rpph1. Moreover, overexpression of Rpph1 increased dendritic spine density in primary cultured hippocampal pyramidal neurons, whereas knocking down of Rpph1 had the reverse effect. In conclusion, Rpph1 modulates CDC42 expression level in a ceRNA-dependent manner, which may represent a compensatory mechanism in the early stage of the AD pathogenesis.

Molecular Diversity of Midbrain Development in Mouse, Human, and Stem Cells

Understanding human embryonic ventral midbrain is of major interest for Parkinson's disease. However, the cell types, their gene expression dynamics, and their relationship to commonly used rodent models remain to be defined. We performed single-cell RNA sequencing to examine ventral midbrain development in human and mouse. We found 25 molecularly defined human cell types, including five subtypes of radial glia-like cells and four progenitors. In the mouse, two mature fetal dopaminergic neuron subtypes diversified into five adult classes during postnatal development. Cell types and gene expression were generally conserved across species, but with clear differences in cell proliferation, developmental timing, and dopaminergic neuron development. Additionally, we developed a method to quantitatively assess the fidelity of dopaminergic neurons derived from human pluripotent stem cells, at a single-cell level. Thus, our study provides insight into the molecular programs controlling human midbrain development and provides a foundation for the development of cell replacement therapies.

Tissue resolved, gene structure refined equine transcriptome

BACKGROUND

Transcriptome interpretation relies on a good-quality reference transcriptome for accurate quantification of gene expression as well as functional analysis of genetic variants. The current annotation of the horse genome lacks the specificity and sensitivity necessary to assess gene expression especially at the isoform level, and suffers from insufficient annotation of untranslated regions (UTR) usage. We built an annotation pipeline for horse and used it to integrate 1.9 billion reads from multiple RNA-seq data sets into a new refined transcriptome.

RESULTS

This equine transcriptome integrates eight different tissues from 59 individuals and improves gene structure and isoform resolution, while providing considerable tissue-specific information. We utilized four levels of transcript filtration in our pipeline, aimed at producing several transcriptome versions that are suitable for different downstream analyses. Our most refined transcriptome includes 36,876 genes and 76,125 isoforms, with 6474 candidate transcriptional loci novel to the equine transcriptome.

CONCLUSIONS

We have employed a variety of descriptive statistics and figures that demonstrate the quality and content of the transcriptome. The equine transcriptomes that are provided by this pipeline show the best tissue-specific resolution of any equine transcriptome to date and are flexible for several downstream analyses. We encourage the integration of further equine transcriptomes with our annotation pipeline to continue and improve the equine transcriptome.