Identification and properties of 1,119 candidate lincRNA loci in the Drosophila melanogaster genome

Untangling the gordian knot: The intertwining interactions between developmental hormone signaling and epigenetic mechanisms in insects

Hormones control developmental and physiological processes, often by regulating the expression of multiple genes simultaneously or sequentially. Crosstalk between hormones and epigenetics is pivotal to dynamically coordinate this process. Hormonal signals can guide the addition and removal of epigenetic marks, steering gene expression. Conversely, DNA methylation, histone modifications and non-coding RNAs can modulate regional chromatin structure and accessibility and regulate the expression of numerous (hormone-related) genes. Here, we provide a review of the interplay between the classical insect hormones, ecdysteroids and juvenile hormones, and epigenetics. We summarize the mode-of-action and roles of these hormones in post-embryonic development, and provide a general overview of epigenetic mechanisms. We then highlight recent advances on the interactions between these hormonal pathways and epigenetics, and their involvement in development. Furthermore, we give an overview of several 'omics techniques employed in the field. Finally, we discuss which questions remain unanswered and possible avenues for future research.

Antennapedia: The complexity of a master developmental transcription factor

Hox genes encode transcription factors that play an important role in establishing the basic body plan of animals. In Drosophila, Antennapedia is one of the five genes that make up the Antennapedia complex (ANT-C). Antennapedia determines the identity of the second thoracic segment, known as the mesothorax. Misexpression of Antennapedia at different developmental stages changes the identity of the mesothorax, including the muscles, nervous system, and cuticle. In Drosophila, Antennapedia has two distinct promoters highly regulated throughout development by several transcription factors. Antennapedia proteins are found with other transcription factors in different ANTENNAPEDIA transcriptional complexes to regulate multiple subsets of target genes. In this review, we describe the different mechanisms that regulate the expression and function of Antennapedia and the role of this Hox gene in the development of Drosophila.

Identification of long noncoding RNAs of silkworm at the early stage of Bombyx mori bidensovirus infection

Bombyx mori bidensovirus (BmBDV) is one of the most important pathogens of silkworm. It mainly infects midgut cells of silkworm and causes losses to the sericulture industry. Long noncoding RNAs (lncRNAs) have been reported to play an important role in the regulation of antiviral immune response in silkworm. To explore whether lncRNAs are involved in BmBDV infection and immune response of silkworm, we performed a comparative transcriptome analysis to identify the lncRNAs and mRNAs between the BmBDV infected and noninfected silkworm larvae at the early stage. A total of 16,069 genes and 974 candidate lncRNAs were identified, among which 142 messenger RNA (mRNAs) and four lncRNAs were differentially expressed (DE). Target gene prediction revealed that 142 DEmRNAs were coexpressed with four DElncRNAs, suggesting that the expression of mRNA is mainly affected through trans-regulation activities. A regulatory network of DElncRNAs and DEmRNAs was constructed, showing that many genes targeted by different DElncRNAs are involved in metabolism and immunity, which implies that these genes and lncRNAs play an important role in the replication of BmBDV. Our results will help us to improve our understanding of lncRNA-mediated regulatory roles in BmBDV infection, providing a new perspective for further exploring the interaction between host and BmBDV.

Comparative RNA Genomics

Over the last quarter of a century it has become clear that RNA is much more than just a boring intermediate in protein expression. Ancient RNAs still appear in the core information metabolism and comprise a surprisingly large component in bacterial gene regulation. A common theme with these types of mostly small RNAs is their reliance of conserved secondary structures. Large-scale sequencing projects, on the other hand, have profoundly changed our understanding of eukaryotic genomes. Pervasively transcribed, they give rise to a plethora of large and evolutionarily extremely flexible non-coding RNAs that exert a vastly diverse array of molecule functions. In this chapter we provide a-necessarily incomplete-overview of the current state of comparative analysis of non-coding RNAs, emphasizing computational approaches as a means to gain a global picture of the modern RNA world.

The pleiotropic functions of Pri smORF peptides synchronize leg development regulators

The last decade witnesses the emergence of the abundant family of smORF peptides, encoded by small ORF (<100 codons), whose biological functions remain largely unexplored. Bioinformatic analyses here identify hundreds of putative smORF peptides expressed in Drosophila imaginal leg discs. Thanks to a functional screen in leg, we found smORF peptides involved in morphogenesis, including the pioneer smORF peptides Pri. Since we identified its target Ubr3 in the epidermis and pri was known to control leg development through poorly understood mechanisms, we investigated the role of Ubr3 in mediating pri function in leg. We found that pri plays several roles during leg development both in patterning and in cell survival. During larval stage, pri activates independently of Ubr3 tarsal transcriptional programs and Notch and EGFR signaling pathways, whereas at larval pupal transition, Pri peptides cooperate with Ubr3 to insure cell survival and leg morphogenesis. Our results highlight Ubr3 dependent and independent functions of Pri peptides and their pleiotropy. Moreover, we reveal that the smORF peptide family is a reservoir of overlooked developmental regulators, displaying distinct molecular functions and orchestrating leg development.

LincRNA_XR209691.3 could promote Bombyx mori nucleopolyhedrovirus replication by interacting with BmHSP70

Long non-coding RNAs (lncRNAs), a class of transcripts exceeding 200 nucleotides and lacking protein coding potential, have been proven to play important roles in viral infection and host immunity. Bombyx mori nucleopolyhedrovirus (BmNPV) is an important pathogen, which causes the silkworm disease and leads to a huge challenge to the sericultural industry. At present, research on the roles of insect lncRNAs in host-virus interaction are relatively few. In this study, we explored the function of lincRNA_XR209691.3 that was significantly up-regulated in the silkworm fat body upon BmNPV infection. Firstly, the subcellular localization experiment confirmed that lincRNA_XR209691.3 was present in both the nucleus and cytoplasm. Enhancing the expression of lincRNA_XR209691.3 in BmN cells could promote the proliferation of BmNPV, while inhibition of lincRNA_XR209691.3 by RNA interference suppresses the proliferation of BmNPV. Combining RNA pull-down and mass spectrometry, we identified the host and BmNPV proteins that could interact with lincRNA_XR209691.3. Next, by using truncation experiment and RNA immunoprecipitation (RIP) assay, it was found that lincRNA_XR209691.3 could bind to the Actin domain of BmHSP70. Subsequently, overexpression of lncRNA_XR209691.3 in BmN cells promoted the expression of BmHSP70, while knockdown of BmHsp70 suppressed the replication of BmNPV. Based on the above results, it is speculated that lincRNA_XR209691.3 could promote the proliferation of BmNPV through interaction with BmHSP70, possibly by improving the stability of BmHSP70 and thereby enhancing the expression of BmHSP70. Our results shed light on the lncRNA function in insect-pathogen interactions and provide a new clue to elucidate the molecular mechanism of BmNPV infection.

Recent Advances and Future Potential of Long Non-Coding RNAs in Insects

Over the last decade, long non-coding RNAs (lncRNAs) have witnessed a steep rise in interest amongst the scientific community. Because of their functional significance in several biological processes, i.e., alternative splicing, epigenetics, cell cycle, dosage compensation, and gene expression regulation, lncRNAs have transformed our understanding of RNA's regulatory potential. However, most knowledge concerning lncRNAs comes from mammals, and our understanding of the potential role of lncRNAs amongst insects remains unclear. Technological advances such as RNA-seq have enabled entomologists to profile several hundred lncRNAs in insect species, although few are functionally studied. This article will review experimentally validated lncRNAs from different insects and the lncRNAs identified via bioinformatic tools. Lastly, we will discuss the existing research challenges and the future of lncRNAs in insects.

Analysis of Long Non-Coding RNA-Mediated Regulatory Networks of Plutella xylostella in Response to Metarhizium anisopliae Infection

Long non-coding RNAs (lncRNAs) represent a diverse class of RNAs that are structurally similar to messenger RNAs (mRNAs) but do not encode proteins. Growing evidence suggests that in response to biotic and abiotic stresses, the lncRNAs play crucial regulatory roles in plants and animals. However, the potential role of lncRNAs during fungal infection has yet to be characterized in Plutella xylostella, a devastating pest of cruciferous crops. In the current study, we performed a strand-specific RNA sequencing of Metarhizium anisopliae-infected (Px36hT, Px72hT) and uninfected (Px36hCK, Px72hCK) P. xylostella fat body tissues. Comprehensive bioinformatic analysis revealed a total of 5665 and 4941 lncRNAs at 36 and 72-h post-infection (hpi), including 563 (Px36hT), 532 (Px72hT) known and 5102 (Px36hT), 4409 (Px72hT) novel lncRNA transcripts. These lncRNAs shared structural similarities with their counterparts in other species, including shorter exon and intron length, fewer exon numbers, and a lower expression profile than mRNAs. LncRNAs regulate the expression of neighboring protein-coding genes by acting in a cis and trans manner. Functional annotation and pathway analysis of cis-acting lncRNAs revealed their role in several immune-related genes, including Toll, serpin, transferrin, βGRP etc. Furthermore, we identified multiple lncRNAs acting as microRNA (miRNA) precursors. These miRNAs can potentially regulate the expression of mRNAs involved in immunity and development, suggesting a crucial lncRNA-miRNA-mRNA complex. Our findings will provide a genetic resource for future functional studies of lncRNAs involved in P. xylostella immune responses to M. anisopliae infection and shed light on understanding insect host-pathogen interactions.

Identification of the stress granule transcriptome via RNA-editing in single cells and in vivo

Stress granules are phase-separated assemblies formed around RNAs. So far, the techniques available to identify these RNAs are not suitable for single cells and small tissues displaying cell heterogeneity. Here, we used TRIBE (target of RNA-binding proteins identified by editing) to profile stress granule RNAs. We used an RNA-binding protein (FMR1) fused to the catalytic domain of an RNA-editing enzyme (ADAR), which coalesces into stress granules upon oxidative stress. RNAs colocalized with this fusion are edited, producing mutations that are detectable by VASA sequencing. Using single-molecule FISH, we validated that this purification-free method can reliably identify stress granule RNAs in bulk and single S2 cells and in Drosophila neurons. Similar to mammalian cells, we find that stress granule mRNAs encode ATP binding, cell cycle, and transcription factors. This method opens the possibility to identify stress granule RNAs and other RNA-based assemblies in other single cells and tissues.

Genomic and functional conservation of lncRNAs: lessons from flies

Over the last decade, the increasing interest in long non-coding RNAs (lncRNAs) has led to the discovery of these transcripts in multiple organisms. LncRNAs tend to be specifically, and often lowly, expressed in certain tissues, cell types and biological contexts. Although lncRNAs participate in the regulation of a wide variety of biological processes, including development and disease, most of their functions and mechanisms of action remain unknown. Poor conservation of the DNA sequences encoding for these transcripts makes the identification of lncRNAs orthologues among different species very challenging, especially between evolutionarily distant species such as flies and humans or mice. However, the functions of lncRNAs are unexpectedly preserved among different species supporting the idea that conservation occurs beyond DNA sequences and reinforcing the potential of characterising lncRNAs in animal models. In this review, we describe the features and roles of lncRNAs in the fruit fly Drosophila melanogaster, focusing on genomic and functional comparisons with human and mouse lncRNAs. We also discuss the current state of advances and limitations in the study of lncRNA conservation and future perspectives.

Identification and the potential roles of long non-coding RNAs in regulating acetyl-CoA carboxylase ACC transcription in spirotetramat-resistant Aphis gossypii

Long non-coding RNAs (lncRNAs) represent the largest class of non-coding transcripts. They act a pivotal part in various insect developmental processes and stress responses. However, the investigation of lncRNA functions in insecticide resistant remains at an early phase. Herein, we conducted whole-transcriptome RNA sequencing for two cotton aphid (Aphis gossypii Glover) strains, i.e., insecticide-susceptible (SS) and spirotetramat-resistant (SR). We discovered 6059 lncRNAs in the RNA-Seq data, and 874 lncRNAs showed differential expression. In addition, 5 lncRNAs among 874 lncRNAs were predicted as targets of acetyl-CoA carboxylase (ACC). Reverse transcription real-time quantitative PCR (RT-qPCR) combined with RNA interference (RNAi) confirmed that selected ACC lncRNA was related to the expression of ACC. Moreover, we also identified two transcription factors, i.e., C/EBP and C/EBPzeta, that regulate the transcription level of ACC lncRNA. These results provide a good basis for the study of cotton aphid lncRNA functions in insecticide resistance development.

Emerging role of non-coding RNA in health and disease

Human diseases have always been a significant turf of concern since the origin of mankind. It is cardinal to know the cause, treatment, and cure for every disease condition. With the advent and advancement in technology, the molecular arena at the microscopic level to study the mechanism, progression, and therapy is more rational and authentic pave than a macroscopic approach. Non-coding RNAs (ncRNAs) have now emerged as indispensable players in the diagnosis, development, and therapeutics of every abnormality concerning physiology, pathology, genetics, epigenetics, oncology, and developmental diseases. This is a comprehensive attempt to collate all the existing and proven strategies, techniques, mechanisms of genetic disorders including Silver Russell Syndrome, Fascio- scapula humeral muscular dystrophy, cardiovascular diseases (atherosclerosis, cardiac fibrosis, hypertension, etc.), neurodegenerative diseases (Spino-cerebral ataxia type 7, Spino-cerebral ataxia type 8, Spinal muscular atrophy, Opitz-Kaveggia syndrome, etc.) cancers (cervix, breast, lung cancer, etc.), and infectious diseases (viral) studied so far. This article encompasses discovery, biogenesis, classification, and evolutionary prospects of the existence of this junk RNA along with the integrated networks involving chromatin remodelling, dosage compensation, genome imprinting, splicing regulation, post-translational regulation and proteomics. In conclusion, all the major human diseases are discussed with a facilitated technology transfer, advancements, loopholes, and tentative future research prospects have also been proposed.

Identification and characterization of long noncoding RNAs and their association with acquisition of blood meal in Culex quinquefasciatus

The Southern house mosquito, Culex quinquefasciatus (Cx. quinquefasciatus) is an important vector that transmit multiple diseases including West Nile encephalitis, Japanese encephalitis, St. Louis encephalitis and lymphatic filariasis. Long noncoding RNAs (lncRNAs) involve in many biological processes such as development, infection, and virus-host interaction. However, there is no systematic identification and characterization of lncRNAs in Cx. quinquefasciatus. Here, we report the first lncRNA identification in Cx. quinquefasciatus. By using 31 public RNA-seq datasets, a total of 4763 novel lncRNA transcripts were identified, of which 3591, 569, and 603 were intergenic, intronic, and antisense respectively. Examination of genomic features revealed that Cx. quinquefasciatus shared similar characteristics with other species such as short in length, low GC content, low sequence conservation, and low coding potential. Furthermore, compared to protein-coding genes, Cx. quinquefasciatus lncRNAs had lower expression values, and tended to be expressed in temporally specific fashion. In addition, weighted correlation network and functional annotation analyses showed that lncRNAs may have roles in blood meal acquisition of adult female Cx. quinquefasciatus mosquitoes. This study presents the first systematic identification and analysis of Cx. quinquefasciatus lncRNAs and their association with blood feeding. Results generated from this study will facilitate future investigation on the function of Cx. quinquefasciatus lncRNAs.

Expression profile of long non-coding RNAs in porcine lymphnode response to porcine circovirus type 2 infection

Porcine circovirus type 2 (PCV2) can cause various clinical diseases in pigs, resulting in huge losses for the pig farms all over the world. In order to develop a new strategy to control PCV2, it is essential to understand its mechanisms firstly, especially PCV2 interferes with the host's innate immunity. In the present study, lncRNA and mRNA expression profiles in porcine lymphnode response to PCV2 infection were deeply sequenced and analyzed. 3271 novel lncRNAs were identified in all. 1898 mRNAs and 282 lncRNAs showed differential expression between control and PCV2-infected groups. The bioinformatics analysis including lncRNA-mRNA co-expression network construction, as well as GO and KEGG pathway analysis focused on the DEGs was carried out. The results indicated that lncRNAs might participate in PCV2 infection-induced the pathogenesis of immunosuppression through regulating the host's immune responses, biological regulation, response to stimulus, cellular component organization or biogenesis and metabolism. And these differentially expressed lncRNAs might play important roles in response to PCV2 infection in the host's innate immune system. These findings provided a large-scale survey of dysregulated lncRNAs after PCV2 infection, especially the lncRNAs responded to host's innate immune within the lymphnode. This study will provide a novel insight into the lncRNAs' functions and the possible immunosuppressive mechanism induced by PCV2 infection. However, further research will be required to verify the characteristic function of the dysregulated lncRNAs.

Transcriptome analysis reveals potential function of long non-coding RNAs in 20-hydroxyecdysone regulated autophagy in Bombyx mori

Background

20-hydroxyecdysone (20E) plays important roles in insect molting and metamorphosis. 20E-induced autophagy has been detected during the larval–pupal transition in different insects. In Bombyx mori, autophagy is induced by 20E in the larval fat body. Long non-coding RNAs (lncRNAs) function in various biological processes in many organisms, including insects. Many lncRNAs have been reported to be potential for autophagy occurrence in mammals, but it has not been investigated in insects.

Results

RNA libraries from the fat body of B. mori dissected at 2 and 6 h post-injection with 20E were constructed and sequenced, and comprehensive analysis of lncRNAs and mRNAs was performed. A total of 1035 lncRNAs were identified, including 905 lincRNAs and 130 antisense lncRNAs. Compared with mRNAs, lncRNAs had longer transcript length and fewer exons. 132 lncRNAs were found differentially expressed at 2 h post injection, compared with 64 lncRNAs at 6 h post injection. Thirty differentially expressed lncRNAs were common at 2 and 6 h post-injection, and were hypothesized to be associated with the 20E response. Target gene analysis predicted 6493 lncRNA-mRNA cis pairs and 42,797 lncRNA-mRNA trans pairs. The expression profiles of LNC_000560 were highly consistent with its potential target genes, Atg4B, and RNAi of LNC_000560 significantly decreased the expression of LNC_000560 and Atg4B. These results indicated that LNC_000560 was potentially involved in the 20E-induced autophagy of the fat body by regulating Atg4B.

Conclusions

This study provides the genome-wide identification and functional characterization of lncRNAs associated with 20E-induced autophagy in the fat body of B. mori. LNC_000560 and its potential target gene were identified to be related to 20-regulated autophagy in B. mori. These results will be helpful for further studying the regulatory mechanisms of lncRNAs in autophagy and other biological processes in this insect model.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07692-1.

Genome-wide identification of long non-coding RNAs (lncRNAs) associated with malathion resistance in Bactrocera dorsalis

BACKGROUND

Long non-coding RNAs (lncRNAs) play important roles in the regulation of biological processes and have been identified in many species including insects. However, the association between lncRNAs and pesticide resistance in insect species such as Bactrocera dorsalis is unknown.

RESULTS

RNA-seq was performed on malathion resistant (MR1) and susceptible (MS) strains of B. dorsalis and a total of 6171 lncRNAs transcripts were identified. These included 3728 lincRNAs, 653 antisense lncRNAs, 1402 intronic lncRNAs, and 388 sense lncRNAs. A total of 40 and 52 upregulated lncRNAs were found in females and males of the MR1 strain compared to 54 and 49 in the same sexes of the MS strain, respectively. Twenty-seven of these lncRNAs showed the same trend of expression in both females and males in the MR1 strain, in which 15 lncRNAs were upregulated and 12 were downregulated. RT-qPCR results indicated that the differentially expressed lncRNAs were associated with malathion resistance. The lnc15010.10 and lnc3774.2 were highly expressed in the cuticle of the MR1 strain, indicating that these two lncRNAs may be related to malathion resistance. RNAi of lnc3774.2 and a bioassay showed that malathion resistance was possibly influenced by changes in the B. dorsalis cuticle.

CONCLUSION

LncRNAs of B. dorsalis potentially related to the malathion resistance were identified. Two lncRNAs appear to influence malathion resistance via modulating the structure, or components, of the cuticle. © 2021 Society of Chemical Industry.

Systematic and computational identification of Androctonus crassicauda long non-coding RNAs

The potential function of long non-coding RNAs in regulating neighbor protein-coding genes has attracted scientists' attention. Despite the important role of lncRNAs in biological processes, a limited number of studies focus on non-model animal lncRNAs. In this study, we used a stringent step-by-step filtering pipeline and machine learning-based tools to identify the specific Androctonus crassicauda lncRNAs and analyze the features of predicted scorpion lncRNAs. 13,401 lncRNAs were detected using pipeline in A. crassicauda transcriptome. The blast results indicated that the majority of these lncRNAs sequences (12,642) have no identifiable orthologs even in closely related species and those considered as novel lncRNAs. Compared to lncRNA prediction tools indicated that our pipeline is a helpful approach to distinguish protein-coding and non-coding transcripts from RNA sequencing data of species without reference genomes. Moreover, analyzing lncRNA characteristics in A. crassicauda uncovered that lower protein-coding potential, lower GC content, shorter transcript length, and less number of isoform per gene are outstanding features of A. crassicauda lncRNAs transcripts.

Filtering the Junk: Assigning Function to the Mosquito Non-Coding Genome

The portion of the mosquito genome that does not code for proteins contains regulatory elements that likely underlie variation for important phenotypes including resistance and susceptibility to infection with arboviruses and Apicomplexan parasites. Filtering the non-coding genome to uncover these functional elements is an expanding area of research, though identification of non-coding regulatory elements is challenging due to the lack of an amino acid-like code for the non-coding genome and a lack of sequence conservation across species. This review focuses on three types of non-coding regulatory elements: (1) microRNAs (miRNAs), (2) long non-coding RNAs (lncRNAs), and (3) enhancers, and summarizes current advances in technical and analytical approaches for measurement of each of these elements on a genome-wide scale. The review also summarizes and highlights novel findings following application of these techniques in mosquito-borne disease research. Looking beyond the protein-coding genome is essential for understanding the complexities that underlie differential gene expression in response to arboviral or parasite infection in mosquito disease vectors. A comprehensive understanding of the regulation of gene and protein expression will inform transgenic and other vector control methods rooted in naturally segregating genetic variation.

The landscape of lncRNAs in Cydia pomonella provides insights into their signatures and potential roles in transcriptional regulation

Background

Long noncoding RNAs (lncRNAs) have emerged as an important class of transcriptional regulators in cellular processes. The past decades have witnessed great progress in lncRNA studies in a variety of organisms. The codling moth (Cydia pomonella L.) is an important invasive insect in China. However, the functional impact of lncRNAs in this insect remains unclear. In this study, an atlas of codling moth lncRNAs was constructed based on publicly available RNA-seq datasets.

Results

In total, 9875 lncRNA transcripts encoded by 9161 loci were identified in the codling moth. As expected, the lncRNAs exhibited shorter transcript lengths, lower GC contents, and lower expression levels than protein-coding genes (PCGs). Additionally, the lncRNAs were more likely to show tissue-specific expression patterns than PCGs. Interestingly, a substantial fraction of the lncRNAs showed a testis-biased expression pattern. Additionally, conservation analysis indicated that lncRNA sequences were weakly conserved across insect species, though additional lncRNAs with homologous relationships could be identified based on synteny, suggesting that synteny could be a more reliable approach for the cross-species comparison of lncRNAs. Furthermore, the correlation analysis of lncRNAs with neighbouring PCGs indicated a stronger correlation between them, suggesting potential cis-acting roles of these lncRNAs in the regulation of gene expression.

Conclusions

Taken together, our work provides a valuable resource for the comparative and functional study of lncRNAs, which will facilitate the understanding of their mechanistic roles in transcriptional regulation.

Genome-wide characterization and developmental expression profiling of long non-coding RNAs in Sogatella furcifera

The genome-wide characterization of long non-coding RNA (lncRNA) in insects demonstrates their importance in fundamental biological processes. Essentially, an in-depth understanding of the functional repertoire of lncRNA in insects is pivotal to insect resources utilization and sustainable pest control. Using a custom bioinformatics pipeline, we identified 1861 lncRNAs encoded by 1852 loci in the Sogatella furcifera genome. We profiled lncRNA expression in different developmental stages and observed that the expression of lncRNAs is more highly temporally restricted compared to protein-coding genes. More up-regulated Sogatella furcifera lncRNA expressed in the embryo, 4th and 5th instars, suggesting that increased lncRNA levels may play a role in these developmental stages. We compared the relationship between the expression of Sogatella furcifera lncRNA and its nearest protein gene and found that lncRNAs were more correlated to their downstream coding neighbors on the opposite strand. Our genome-wide profiling of lncRNAs in Sogatella furcifera identifies exciting candidates for characterization of lncRNAs, and also provides information on lncRNA regulation during insect development.

An integrative proteogenomics approach reveals peptides encoded by annotated lincRNA in the mouse kidney inner medulla

Long noncoding RNAs (lncRNAs) are intracellular transcripts longer than 200 nucleotides and lack protein-coding information. A subclass of lncRNA known as long intergenic noncoding RNAs (lincRNAs) are transcribed from genomic regions that share no overlap with annotated protein-coding genes. Increasing evidence has shown that some annotated lincRNA transcripts do in fact contain open reading frames (ORFs) encoding functional short peptides in the cell. Few robust methods for lincRNA-encoded peptide identification have been reported, and the tissue-specific expression of these peptides has been largely unexplored. Here we propose an integrative workflow for lincRNA-encoded peptide discovery and test it on the mouse kidney inner medulla (IM). In brief, low molecular weight protein fractions were enriched from homogenate of IMs and trypsinized into shorter peptides, which were sequenced by high resolution liquid chromatography-tandem mass spectrometry (LC-MS/MS). To curate a hypothetical lincRNA-encoded peptide database for peptide-spectrum matching following LC-MS/MS, we performed RNA-Seq on IMs, computationally removed reads overlapping with annotated protein-coding genes, and remapped the remaining reads to a database of mouse noncoding transcripts to infer lincRNA expression. Expressed lincRNAs were searched for ORFs by an existing rule-based algorithm, and translated ORFs were used for peptide-spectrum matching. Peptides identified by LC-MS/MS were further evaluated by using several quality control criteria and bioinformatics methods. We discovered three novel lincRNA-encoded peptides, which are conserved in mouse, rat, and human. The workflow can be adapted for discovery of small protein-coding genes in any species or tissue where noncoding transcriptome information is available.

Characterization and analysis of full-length transcriptomes from two grasshoppers, Gomphocerus licenti and Mongolotettix japonicus

Acrididae are diverse in size, body shape, behavior, ecology and life history; widely distributed; easy to collect; and important to agriculture. They represent promising model candidates for functional genomics, but their extremely large genomes have hindered this research; establishing a reference transcriptome for a species is the primary means of obtaining genetic information. Here, two Acrididae species, Gomphocerus licenti and Mongolotettix japonicus, were selected for full-length (FL) PacBio transcriptome sequencing. For G. licenti and M. japonicus, respectively, 590,112 and 566,165 circular consensus sequences (CCS) were generated, which identified 458,131 and 428,979 full-length nonchimeric (FLNC) reads. After isoform-level clustering, next-generation sequencing (NGS) short sequences were used for error correction, and remove redundant sequences with CD-HIT, 17,970 and 16,766 unigenes were generated for G. licenti and M. japonicus. In addition, we obtained 17,495 and 16,373 coding sequences, 1,082 and 813 transcription factors, 11,840 and 10,814 simple sequence repeats, and 905 and 706 long noncoding RNAs by analyzing the transcriptomes of G. licenti and M. japonicus, respectively, and 15,803 and 14,846 unigenes were annotated in eight functional databases. This is the first study to sequence FL transcriptomes of G. licenti and M. japonicus, providing valuable genetic resources for further functional genomics research.

Comprehensive Identification and Characterization of Long Non-coding RNAs Associated With Rice Black-Streaked Dwarf Virus Infection in Laodelphax striatellus (Fallén) Midgut

Long non-coding RNAs (lncRNAs) are involved in a variety of biological functions through transcriptional and post-transcriptional regulation. However, little is known about their functions in the process of insect mediated virus transmission. In the present study, we analyzed using RNA-Seq, the lncRNAs that were differentially expressed in response to Rice black-streaked dwarf virus (RBSDV) infection in Laodelphax striatellus (Fallén) midgut. A total of 13,927 lncRNAs were identified and over 69% were assigned to intergenic regions. Among them, 176 lncRNAs were differentially expressed and predicted to target 168 trans-regulatory genes. Ten differentially expressed lncRNAs were selected and their expression changes were validated by RT-qPCR. KEGG analysis showed that these target genes were enriched in the essential biological process, such as Purine metabolism, Valine, leucine and isoleucine degradation, and Fatty acid elongation. The expression levels of the differentially expressed lncRNAs and the predicted target genes that were significantly enriched in the Human papillomavirus infection pathway were analyzed by RT-qPCR. The results showed that several lncRNAs were co-expressed with their target genes. One of the lncRNAs called MSTRG15394 and its target gene, encoding a secreted protease inhibitor (PI), were up-regulated at the transcriptional level after RBSDV infection. Knockdown of MSTRG15394 could down-regulate the PI expression at mRNA level. Inhibition of either MSTRG15394 or PI expression by RNA interference promoted RBSDV accumulation in L. striatellus midgut. Our finding provides new insights into the function of lncRNAs in regulating virus infection in an important insect vector.

Tunable self-cleaving ribozymes for modulating gene expression in eukaryotic systems

Advancements in the field of synthetic biology have been possible due to the development of genetic tools that are able to regulate gene expression. However, the current toolbox of gene regulatory tools for eukaryotic systems have been outpaced by those developed for simple, single-celled systems. Here, we engineered a set of gene regulatory tools by combining self-cleaving ribozymes with various upstream competing sequences that were designed to disrupt ribozyme self-cleavage. As a proof-of-concept, we were able to modulate GFP expression in mammalian cells, and then showed the feasibility of these tools in Drosophila embryos. For each system, the fold-reduction of gene expression was influenced by the location of the self-cleaving ribozyme/upstream competing sequence (i.e. 5' vs. 3' untranslated region) and the competing sequence used. Together, this work provides a set of genetic tools that can be used to tune gene expression across various eukaryotic systems.

Long non-coding RNAs in development and disease: conservation to mechanisms

Our genomes contain the blueprint of what makes us human and many indications as to why we develop disease. Until the last 10 years, most studies had focussed on protein-coding genes, more specifically DNA sequences coding for proteins. However, this represents less than 5% of our genomes. The other 95% is referred to as the 'dark matter' of our genomes, our understanding of which is extremely limited. Part of this 'dark matter' includes regions that give rise to RNAs that do not code for proteins. A subset of these non-coding RNAs are long non-coding RNAs (lncRNAs), which in particular are beginning to be dissected and their importance to human health revealed. To improve our understanding and treatment of disease it is vital that we understand the molecular and cellular function of lncRNAs, and how their misregulation can contribute to disease. It is not yet clear what proportion of lncRNAs is actually functional; conservation during evolution is being used to understand the biological importance of lncRNA. Here, we present key themes within the field of lncRNAs, emphasising the importance of their roles in both the nucleus and the cytoplasm of cells, as well as patterns in their modes of action. We discuss their potential functions in development and disease using examples where we have the greatest understanding. Finally, we emphasise why lncRNAs can serve as biomarkers and discuss their emerging potential for therapy. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Identification of LincRNA from Dermatophagoides farinae (Acari: Pyroglyphidae) for Potential Allergen-Related Targets

Long noncoding RNAs (lncRNAs), especially their important subclass of long intergenic noncoding RNAs (lincRNAs), have been identified in some insects. They play important roles in the regulation of biological processes, such as immune response or cell differentiation and as possible evolutionary precursors for protein coding genes. House dust mites (HDMs) are recognized as allergenic mites because allergens are found in their feces and bodies. Dermatophagoides farinae is one of the most important pyroglyphid mites because of its abundance in the household. To determine if lincRNAs can regulate allergen presentation in HDMs, we analyzed RNA-seq data for HDMs. We identified 11 lincRNAs that are related to mRNAs coding for allergens in HDMs. Using qRT-PCR, we amplified 10 lincRNAs and their putative target allergen-encoding mRNAs, confirming expression of these lincRNAs and allergen genes. The results suggest that lincRNAs might be involved in the regulation of allergen production in HDMs and might represent potential acaricidal candidates to inhibit mite allergen production.

Comparative analysis of dsRNA-induced lncRNAs in three kinds of insect species

Long noncoding RNAs (lncRNAs) that have immune responses to various stimuli have been identified in some insects. One type of pathogen-associated molecular pattern, double-stranded RNA (dsRNA), can trigger the RNA interference (RNAi) pathway and immune response. Interestingly, there has been no research into characterizing the relationship between lncRNA and dsRNA-induced RNAi pathways. In this study, dsRNA-induced lncRNAs were investigated in two species of lepidopteran insects, Helicoverpa armigera and Plutella xylostella, and one species of coleopteran insects, Tribolium castaneum. Between untreated group and dsRNA-induced group; 3,463 H. armigera, 6,245 P. xylostella, and 3,067 T. castaneum differentially expressed lncRNAs were identified while 156 H. armigera, 247 P. xylostella, 415 T. castaneum lncRNAs and their putative target genes showed consistent changes in gene expression. In T. castaneum, most target genes of the differentially expressed lncRNAs are enriched in the cyclic adenosine monophosphate signaling pathway, ABC transporters, and Janus kinase-signal transducers and activators of the transcription signaling pathway. Conversely, in H. armigera and P. xylostella, the differentially expressed lncRNAs were mainly enriched in the metabolic, digestive, and synthetic signaling pathways. This result indicates that dsRNA-induced lncRNA is species-dependent. We also found that both Dicer-2 and the lncRNA that targets Dicer-2 were significantly upregulated after dsRNA treatment in P. xylostella, indicating that some lncRNAs may be involved in the regulation of the core RNAi pathway in insects. Our results are the first to identify a relationship between lncRNAs and dsRNA in various insect species with different RNAi efficiencies. These results provide a reference for future study of the dsRNA-induced RNAi pathway and different RNAi efficiencies among insect species.

High-Quality Genome Assemblies Reveal Long Non-coding RNAs Expressed in Ant Brains

Ants are an emerging model system for neuroepigenetics, as embryos with virtually identical genomes develop into different adult castes that display diverse physiology, morphology, and behavior. Although a number of ant genomes have been sequenced to date, their draft quality is an obstacle to sophisticated analyses of epigenetic gene regulation. We reassembled de novo high-quality genomes for two ant species, Camponotus floridanus and Harpegnathos saltator. Using long reads enabled us to span large repetitive regions and improve genome contiguity, leading to comprehensive and accurate protein-coding annotations that facilitated the identification of a Gp-9-like gene as differentially expressed in Harpegnathos castes. The new assemblies also enabled us to annotate long non-coding RNAs in ants, revealing caste-, brain-, and developmental-stage-specific long non-coding RNAs (lncRNAs) in Harpegnathos. These upgraded genomes, along with the new gene annotations, will aid future efforts to identify epigenetic mechanisms of phenotypic and behavioral plasticity in ants.

Combining next-generation sequencing and single-molecule sequencing to explore brown plant hopper responses to contrasting genotypes of japonica rice

BACKGROUND

The brown plant hopper (BPH), Nilaparvata lugens, is one of the major pest of rice (Oryza sativa). Plant defenses against insect herbivores have been extensively studied, but our understanding of insect responses to host plants' resistance mechanisms is still limited. The purpose of this study is to characterize transcripts of BPH and reveal the responses of BPH insects to resistant rice at transcription level by using the advanced molecular techniques, the next-generation sequencing (NGS) and the single-molecule, real-time (SMRT) sequencing.

RESULTS

The current study obtained 24,891 collapsed isoforms of full-length transcripts, and 20,662 were mapped to known annotated genes, including 17,175 novel transcripts. The current study also identified 915 fusion genes, 1794 novel genes, 2435 long non-coding RNAs (lncRNAs), and 20,356 alternative splicing events. Moreover, analysis of differentially expressed genes (DEGs) revealed that genes involved in metabolic and cell proliferation processes were significantly enriched in up-regulated and down-regulated sets, respectively, in BPH fed on resistant rice relative to BPH fed on susceptible wild type rice. Furthermore, the FoxO signaling pathway was involved and genes related to BPH starvation response (Nlbmm), apoptosis and autophagy (caspase 8, ATG13, BNIP3 and IAP), active oxygen elimination (catalase, MSR, ferritin) and detoxification (GST, CarE) were up-regulated in BPH responses to resistant rice.

CONCLUSIONS

The current study provides the first demonstrations of the full diversity and complexity of the BPH transcriptome, and indicates that BPH responses to rice resistance, might be related to starvation stress responses, nutrient transformation, oxidative decomposition, and detoxification. The current result findings will facilitate further exploration of molecular mechanisms of interaction between BPH insects and host rice.

Systematic identification and characterization of Aedes aegypti long noncoding RNAs (lncRNAs)

Long noncoding RNAs (lncRNAs) play diverse roles in biological processes. Aedes aegypti (Ae. aegypti), a blood-sucking mosquito, is the principal vector responsible for replication and transmission of arboviruses including dengue, Zika, and Chikungunya virus. Systematic identification and developmental characterisation of Ae. aegypti lncRNAs are still limited. We performed genome-wide identification of lncRNAs, followed by developmental profiling of lncRNA in Ae. aegypti. We identified a total of 4,689 novel lncRNA transcripts, of which 2,064, 2,076, and 549 were intergenic, intronic, and antisense respectively. Ae. aegypti lncRNAs share many characteristics with other species including low expression, low GC content, short in length, and low conservation. Besides, the expression of Ae. aegypti lncRNAs tend to be correlated with neighbouring and antisense protein-coding genes. A subset of lncRNAs shows evidence of maternal inheritance; hence, suggesting potential role of lncRNAs in early-stage embryos. Additionally, lncRNAs show higher tendency to be expressed in developmental and temporal specific manner. The results from this study provide foundation for future investigation on the function of Ae. aegypti lncRNAs.

The role of long noncoding RNA in major human disease

BACKGROUND

Long noncoding RNAs (lncRNAs) are RNAs whose transcripts are longer than 200nt in length and lack the ability to encode proteins due to lack of specific open reading frames. lncRNAs were once thought to represent transcriptome noise or garbage sequences and a byproduct of RNA polymerase II (Pol II), and thereby ignored by researchers. In fact, lncRNA was involved in a wide variety of physiological and pathological processes in organisms. Comprehensive study of lncRNA does not only provide explanations to the physiological and pathological processes of living organisms, but also gives us new perspectives to the diagnosis, prevention and treatment of some clinical diseases. Therefore, the study of lncRNA is a very broad field of great research value and significance.

RESULTS

This article reviews the function of lncRNAs and their role in major human diseases.

CONCLUSIONS

Numerous studies show that lncRNA might serve as a biomarker for diagnosis and prognosis of various diseases. Compared to conventional biomarkers, lncRNA seems to have a higher diagnostic and prognostic values, not only because of their tissue and disease specific expression patterns, but also due to their highly stable physical and chemical properties.

Genome-Wide Identification of Long Non-Coding RNAs and Their Regulatory Networks Involved in Apis mellifera ligustica Response to Nosema ceranae Infection

Long non-coding RNAs (lncRNAs) are a diverse class of transcripts that structurally resemble mRNAs but do not encode proteins, and lncRNAs have been proven to play pivotal roles in a wide range of biological processes in animals and plants. However, knowledge of expression patterns and potential roles of honeybee lncRNA response to Nosema ceranae infection is completely unknown. Here, we performed whole transcriptome strand-specific RNA sequencing of normal midguts of Apis mellifera ligustica workers (Am7CK, Am10CK) and N. ceranae-inoculated midguts (Am7T, Am10T), followed by comprehensive analyses using bioinformatic and molecular approaches. A total of 6353 A. m. ligustica lncRNAs were identified, including 4749 conserved lncRNAs and 1604 novel lncRNAs. These lncRNAs had minimal sequence similarities with other known lncRNAs in other species; however, their structural features were similar to counterparts in mammals and plants, including shorter exon and intron length, lower exon number, and lower expression level, compared with protein-coding transcripts. Further, 111 and 146 N. ceranae-responsive lncRNAs were identified from midguts at 7-days post-inoculation (dpi) and 10 dpi compared with control midguts. Twelve differentially expressed lncRNAs (DElncRNAs) were shared by Am7CK vs. Am7T and Am10CK vs. Am10T comparison groups, while the numbers of unique DElncRNAs were 99 and 134, respectively. Functional annotation and pathway analysis showed that the DElncRNAs may regulate the expression of neighboring genes by acting in cis and trans fashion. Moreover, we discovered 27 lncRNAs harboring eight known miRNA precursors and 513 lncRNAs harboring 2257 novel miRNA precursors. Additionally, hundreds of DElncRNAs and their target miRNAs were found to form complex competitive endogenous RNA (ceRNA) networks, suggesting that these DElncRNAs may act as miRNA sponges. Furthermore, DElncRNA-miRNA-mRNA networks were constructed and investigated, the results demonstrated that a portion of the DElncRNAs were likely to participate in regulating the host material and energy metabolism as well as cellular and humoral immune host responses to N. ceranae invasion. Our findings revealed here offer not only a rich genetic resource for further investigation of the functional roles of lncRNAs involved in the A. m. ligustica response to N. ceranae infection, but also a novel insight into understanding the host-pathogen interaction during honeybee microsporidiosis.

lncRedibly versatile: biochemical and biological functions of long noncoding RNAs

Long noncoding RNAs (lncRNAs) are transcripts that do not code for proteins, but nevertheless exert regulatory effects on various biochemical pathways, in part via interactions with proteins, DNA, and other RNAs. LncRNAs are thought to regulate transcription and other biological processes by acting, for example, as guides that target proteins to chromatin, scaffolds that facilitate protein-protein interactions and complex formation, and orchestrators of phase-separated compartments. The study of lncRNAs has reached an exciting time, as recent advances in experimental and computational methods allow for genome-wide interrogation of biochemical and biological mechanisms of these enigmatic transcripts. A better appreciation for the biochemical versatility of lncRNAs has allowed us to begin closing gaps in our knowledge of how they act in diverse cellular and organismal contexts, including development and disease.

A gene expression atlas of embryonic neurogenesis in Drosophila reveals complex spatiotemporal regulation of lncRNAs

Cell type specification during early nervous system development in Drosophila melanogaster requires precise regulation of gene expression in time and space. Resolving the programs driving neurogenesis has been a major challenge owing to the complexity and rapidity with which distinct cell populations arise. To resolve the cell type-specific gene expression dynamics in early nervous system development, we have sequenced the transcriptomes of purified neurogenic cell types across consecutive time points covering crucial events in neurogenesis. The resulting gene expression atlas comprises a detailed resource of global transcriptome dynamics that permits systematic analysis of how cells in the nervous system acquire distinct fates. We resolve known gene expression dynamics and uncover novel expression signatures for hundreds of genes among diverse neurogenic cell types, most of which remain unstudied. We also identified a set of conserved long noncoding RNAs (lncRNAs) that are regulated in a tissue-specific manner and exhibit spatiotemporal expression during neurogenesis with exquisite specificity. lncRNA expression is highly dynamic and demarcates specific subpopulations within neurogenic cell types. Our spatiotemporal transcriptome atlas provides a comprehensive resource for investigating the function of coding genes and noncoding RNAs during crucial stages of early neurogenesis.

Summary: DIV-MARIS, an adapted technique for examining stage- and cell type-specific gene expression, reveals a complex network of mRNAs and lncRNAs expressed in specific cell types during early Drosophila embryonic nervous system development.

SSS-test: a novel test for detecting positive selection on RNA secondary structure

BACKGROUND

Long non-coding RNAs (lncRNAs) play an important role in regulating gene expression and are thus important for determining phenotypes. Most attempts to measure selection in lncRNAs have focused on the primary sequence. The majority of small RNAs and at least some parts of lncRNAs must fold into specific structures to perform their biological function. Comprehensive assessments of selection acting on RNAs therefore must also encompass structure. Selection pressures acting on the structure of non-coding genes can be detected within multiple sequence alignments. Approaches of this type, however, have so far focused on negative selection. Thus, a computational method for identifying ncRNAs under positive selection is needed.

RESULTS

We introduce the SSS-test (test for Selection on Secondary Structure) to identify positive selection and thus adaptive evolution. Benchmarks with biological as well as synthetic controls yield coherent signals for both negative and positive selection, demonstrating the functionality of the test. A survey of a lncRNA collection comprising 15,443 families resulted in 110 candidates that appear to be under positive selection in human. In 26 lncRNAs that have been associated with psychiatric disorders we identified local structures that have signs of positive selection in the human lineage.

CONCLUSIONS

It is feasible to assay positive selection acting on RNA secondary structures on a genome-wide scale. The detection of human-specific positive selection in lncRNAs associated with cognitive disorder provides a set of candidate genes for further experimental testing and may provide insights into the evolution of cognitive abilities in humans.

AVAILABILITY

The SSS-test and related software is available at: https://github.com/waltercostamb/SSS-test . The databases used in this work are available at: http://www.bioinf.uni-leipzig.de/Software/SSS-test/ .

Transcriptomic analyses reveal groups of co-expressed, syntenic lncRNAs in four species of the genus Caenorhabditis

Long non-coding RNAs (lncRNAs) are a heterogeneous class of genes that do not code for proteins. Since lncRNAs (or a fraction thereof) are expected to be functional, many efforts have been dedicated to catalog lncRNAs in numerous organisms, but our knowledge of lncRNAs in non vertebrate species remains very limited. Here, we annotated lncRNAs using transcriptomic data from the same larval stage of four Caenorhabditis species. The number of annotated lncRNAs in self-fertile nematodes was lower than in out-crossing species. We used a combination of approaches to identify putatively homologous lncRNAs: synteny, sequence conservation, and structural conservation. We classified a total of 1,532 out of 7,635 genes from the four species into families of lncRNAs with conserved synteny and expression at the larval stage, suggesting that a large fraction of the predicted lncRNAs may be species specific. Despite both sequence and local secondary structure seem to be poorly conserved, sequences within families frequently shared BLASTn hits and short sequence motifs, which were more likely to be unpaired in the predicted structures. We provide the first multi-species catalog of lncRNAs in nematodes and identify groups of lncRNAs with conserved synteny and expression, that share exposed motifs.

Identification of functional long non-coding RNAs in C. elegans

BACKGROUND

Functional characterisation of the compact genome of the model organism Caenorhabditis elegans remains incomplete despite its sequencing 20 years ago. The last decade of research has seen a tremendous increase in the number of non-coding RNAs identified in various organisms. While we have mechanistic understandings of small non-coding RNA pathways, long non-coding RNAs represent a diverse class of active transcripts whose function remains less well characterised.

RESULTS

By analysing hundreds of published transcriptome datasets, we annotated 3392 potential lncRNAs including 143 multi-exonic loci that showed increased nucleotide conservation and GC content relative to other non-coding regions. Using CRISPR/Cas9 genome editing, we generated deletion mutants for ten long non-coding RNA loci. Using automated microscopy for in-depth phenotyping, we show that six of the long non-coding RNA loci are required for normal development and fertility. Using RNA interference-mediated gene knock-down, we provide evidence that for two of the long non-coding RNA loci, the observed phenotypes are dependent on the corresponding RNA transcripts.

CONCLUSIONS

Our results highlight that a large section of the non-coding regions of the C. elegans genome remains unexplored. Based on our in vivo analysis of a selection of high-confidence lncRNA loci, we expect that a significant proportion of these high-confidence regions is likely to have a biological function at either the genomic or the transcript level.

Systematic evaluation of C. elegans lincRNAs with CRISPR knockout mutants

Background

Long intergenic RNAs (lincRNAs) play critical roles in eukaryotic cells, but systematic analyses of the lincRNAs of an animal for phenotypes are lacking. We generate CRISPR knockout strains for Caenorhabditis elegans lincRNAs and evaluate their phenotypes.

Results

C. elegans lincRNAs demonstrate global features such as shorter length and fewer exons than mRNAs. For the systematic evaluation of C. elegans lincRNAs, we produce CRISPR knockout strains for 155 of the total 170 C. elegans lincRNAs. Mutants of 23 lincRNAs show phenotypes in 6 analyzed traits. We investigate these lincRNAs by phenotype for their gene expression patterns and potential functional mechanisms. Some C. elegans lincRNAs play cis roles to modulate the expression of their neighboring genes, and several lincRNAs play trans roles as ceRNAs against microRNAs. We also examine the regulation of lincRNA expression by transcription factors, and we dissect the pathway by which two transcription factors, UNC-30 and UNC-55, together control the expression of linc-73. Furthermore, linc-73 possesses a cis function to modulate the expression of its neighboring kinesin gene unc-104 and thus plays roles in C. elegans locomotion.

Conclusions

By using CRISPR/cas9 technology, we generate knockout strains of 155 C. elegans lincRNAs as valuable resources for studies in noncoding RNAs, and we provide biological insights for 23 lincRNAs with the phenotypes identified in this study.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1619-6) contains supplementary material, which is available to authorized users.

An exact test for comparing a fixed quantitative property between gene sets

Motivation

A significant difference in the distribution of a feature between two gene sets can provide insight into function or regulation. This statistical setting differs from much of hypothesis testing theory because the genome is often considered to be effectively fixed, finite and entirely known in commonly studied organisms, such as human. The Mann-Whitney U test is commonly employed in this scenario despite the assumptions of the test not being met, leading to unreliable and generally underpowered results. Permutation tests are also commonly employed for this purpose, but are computationally burdensome and are not tractable for obtaining small P values or for multiple comparisons.

Results

We present an exact test for the null hypothesis that gene set membership is independent of the quantitative gene feature of interest. We derive an analytic expression for the randomization distribution of the median of the quantitative feature under the null hypothesis. Efficient implementation permits calculation of precise P values of arbitrary magnitude and makes thousands of simultaneous tests of transcriptome-sized gene sets computationally tractable. The flexibility of the hypothesis testing framework presented permits extension to a variety of related tests commonly found in genomics. The exact test is used to identify signatures of translation control and protein function in the human genome.

Availability and implementation

The exact test presented here is implemented in R in the package kpmt available on CRAN.

Contact

map2085@med.cornell.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Identification and characterization of novel conserved RNA structures in Drosophila

BACKGROUND

Comparative genomics approaches have facilitated the discovery of many novel non-coding and structured RNAs (ncRNAs). The increasing availability of related genomes now makes it possible to systematically search for compensatory base changes - and thus for conserved secondary structures - even in genomic regions that are poorly alignable in the primary sequence. The wealth of available transcriptome data can add valuable insight into expression and possible function for new ncRNA candidates. Earlier work identifying ncRNAs in Drosophila melanogaster made use of sequence-based alignments and employed a sliding window approach, inevitably biasing identification toward RNAs encoded in the more conserved parts of the genome.

RESULTS

To search for conserved RNA structures (CRSs) that may not be highly conserved in sequence and to assess the expression of CRSs, we conducted a genome-wide structural alignment screen of 27 insect genomes including D. melanogaster and integrated this with an extensive set of tiling array data. The structural alignment screen revealed ∼30,000 novel candidate CRSs at an estimated false discovery rate of less than 10%. With more than one quarter of all individual CRS motifs showing sequence identities below 60%, the predicted CRSs largely complement the findings of sliding window approaches applied previously. While a sixth of the CRSs were ubiquitously expressed, we found that most were expressed in specific developmental stages or cell lines. Notably, most statistically significant enrichment of CRSs were observed in pupae, mainly in exons of untranslated regions, promotors, enhancers, and long ncRNAs. Interestingly, cell lines were found to express a different set of CRSs than were found in vivo. Only a small fraction of intergenic CRSs were co-expressed with the adjacent protein coding genes, which suggests that most intergenic CRSs are independent genetic units.

CONCLUSIONS

This study provides a more comprehensive view of the ncRNA transcriptome in fly as well as evidence for differential expression of CRSs during development and in cell lines.

Porcine endemic diarrhea virus infection regulates long noncoding RNA expression

Long noncoding RNAs (lncRNAs) have been implicated in various life processes. However, the lncRNA expression and potential functions in porcine endemic diarrhea virus (PEDV) infection and host defense are still poorly understood. In this study, we investigated the lncRNA expression profiles during PEDV infection in intestinal porcine epithelial cell-jejunum 2 (IPEC-J2) cell lines by next-generation sequencing and identified 6188 novel lncRNAs. The functional annotation analysis revealed that these lncRNAs might be associated with many immunity-related genes. We next selected candidate lncRNAs related to immune response pathways and further identified their differential expression in PEDV-infected IPEC-J2 cells and newborn piglets. Our results demonstrated that PEDV infection regulated lncRNA expression patterns in both the IPEC-J2 cell line and piglet ileum. These findings provide the first large-scale survey of lncRNAs associated with PEDV infection, specifically the lncRNAs responsible for the activation of the immune system within the ileum.

Non-coding RNA Expression, Function, and Variation during Drosophila Embryogenesis

Long non-coding RNAs (lncRNAs) can often function in the regulation of gene expression during development; however, their generality as essential regulators in developmental processes and organismal phenotypes remains unclear. Here, we performed a tailored investigation of lncRNA expression and function during Drosophila embryogenesis, interrogating multiple stages, tissue specificity, nuclear localization, and genetic backgrounds. Our results almost double the number of annotated lncRNAs expressed at these embryonic stages. lncRNA levels are generally positively correlated with those of their neighboring genes, with little evidence of transcriptional interference. Using fluorescent in situ hybridization, we report the spatiotemporal expression of 15 new lncRNAs, revealing very dynamic tissue-specific patterns. Despite this, deletion of selected lncRNA genes had no obvious developmental defects or effects on viability under standard and stressed conditions. However, two lncRNA deletions resulted in modest expression changes of a small number of genes, suggesting that they fine-tune expression of non-essential genes. Several lncRNAs have strain-specific expression, indicating that they are not fixed within the population. This intra-species variation across genetic backgrounds may thereby be a useful tool to distinguish rapidly evolving lncRNAs with as yet non-essential roles.

Genome-wide identification and differentially expression analysis of lncRNAs in tilapia

Background

Long noncoding RNAs (LncRNAs) play important roles in fundamental biological processes. However, knowledge about the genome-wide distribution and stress-related expression of lncRNAs in tilapia is still limited.

Results

Genome-wide identification of lncRNAs in the tilapia genome was carried out in this study using bioinformatics tools. 103 RNAseq datasets that generated in our laboratory or collected from NCBI database were analyzed. In total, 72,276 high-confidence lncRNAs were identified. The averaged positive correlation coefficient (r_mean = 0.286) between overlapped lncRNA and mRNA pairs showed significant differences with the values for all lncRNA-mRNA pairs (r_mean = 0.176, z statistics = − 2.45, p value = 0.00071) and mRNA-mRNA pairs (r_mean = 0.186, z statistics = − 2.23, p value = 0.0129). Weighted correlation network analysis of the lncRNA and mRNA datasets from 12 tissues identified 21 modules and many interesting mRNA genes that clustered with lncRNAs. Overrepresentation test indicated that these mRNAs enriched in many biological processes, such as meiosis (p = 0.00164), DNA replication (p = 0.00246), metabolic process (p = 0.000838) and in molecular function, e.g., helicase activity (p = 0.000102) and catalytic activity (p = 0.0000612). Differential expression (DE) analysis identified 99 stress-related lncRNA genes and 1955 tissue-specific DE lncRNA genes. MiRNA-lncRNA interaction analysis detected 72,267 lncRNAs containing motifs with sequence complementary to 458 miRNAs.

Conclusions

This study provides an invaluable resource for further studies on molecular bases of lncRNAs in tilapia genomes. Further function analysis of the lncRNAs will help to elucidate their roles in regulating stress-related adaptation in tilapia.

Electronic supplementary material

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

LncRNA SRA deregulation contributes to the development of atherosclerosis by causing dysfunction of endothelial cells through repressing the expression of adipose triglyceride lipase

It has been well established that long non‑coding RNAs (lncRNAs) are crucial mediators in a diverse range of diseases, including atherosclerosis. The present study aimed to examine the molecular mechanisms underlying the association between steroid receptor RNA activator (SRA) and atherosclerosis. Reverse transcription‑quantitative polymerase chain reaction analysis, western blot analysis and luciferase assays were performed to examine interactions among SRA, adipose triglyceride lipase (ATGL) and peroxisome proliferator‑activated receptor (PPARγ), and the effect of resveratrol (RSV) on the levels of SRA, ATGL and PPARγ. ELISA was performed to determine the effects of SRA and RSV on the production of inflammatory‑associated cytokines. The results showed that knockdown of the expression of SRA by transfecting HUVECs with short hairpin RNA‑SRA inhibited the production of ATGL and PPARγ. A plasmid coding SRA RNA, but not the SRAP protein, attenuated the luciferase activity of the ATGL promoter. PPARγ had no effect on the luciferase activity driven by the ATGL promoter in the absence of rosiglitazone, whereas the luciferase activity of the ATGL promoter was elevated in the presence of rosiglitazone. This effect was eliminated by SRA. SRA enhanced the production of inflammatory‑associated cytokines, including tumor necrosis factor‑α, interleukin‑6, monocyte chemotactic protein‑1 and intercellular adhesion molecule‑1; however, the promoting effect of SRA was eliminated by RSV. RSV increased the expression of ATGL and PPARγ, but not that of SRA. RSV distinctly and concentration‑dependently upregulated the luciferase activity of ATGL, compared with that in the cells without RSV treatment, whereas treating with rosiglitazone inhibited the effect of RSV on the luciferase activity of ATGL. The present study examined the roles of SRA in atherosclerosis, and the effects of changes in SRA and ATGL on inflammatory cytokines and HUVEC dysfunction.

Genome-wide identification of tissue-specific long non-coding RNA in three farm animal species

Background

Numerous long non-coding RNAs (lncRNAs) have been identified and their roles in gene regulation in humans, mice, and other model organisms studied; however, far less research has been focused on lncRNAs in farm animal species. While previous studies in chickens, cattle, and pigs identified lncRNAs in specific developmental stages or differentially expressed under specific conditions in a limited number of tissues, more comprehensive identification of lncRNAs in these species is needed. The goal of the FAANG Consortium (Functional Annotation of Animal Genomes) is to functionally annotate animal genomes, including the annotation of lncRNAs. As one of the FAANG pilot projects, lncRNAs were identified across eight tissues in two adult male biological replicates from chickens, cattle, and pigs.

Results

Comprehensive lncRNA annotations for the chicken, cattle, and pig genomes were generated by utilizing RNA-seq from eight tissue types from two biological replicates per species at the adult developmental stage. A total of 9393 lncRNAs in chickens, 7235 lncRNAs in cattle, and 14,429 lncRNAs in pigs were identified. Including novel isoforms and lncRNAs from novel loci, 5288 novel lncRNAs were identified in chickens, 3732 in cattle, and 4870 in pigs. These transcripts match previously known patterns of lncRNAs, such as generally lower expression levels than mRNAs and higher tissue specificity. An analysis of lncRNA conservation across species identified a set of conserved lncRNAs with potential functions associated with chromatin structure and gene regulation. Tissue-specific lncRNAs were identified. Genes proximal to tissue-specific lncRNAs were enriched for GO terms associated with the tissue of origin, such as leukocyte activation in spleen.

Conclusions

LncRNAs were identified in three important farm animal species using eight tissues from adult individuals. About half of the identified lncRNAs were not previously reported in the NCBI annotations for these species. While lncRNAs are less conserved than protein-coding genes, a set of positionally conserved lncRNAs were identified among chickens, cattle, and pigs with potential functions related to chromatin structure and gene regulation. Tissue-specific lncRNAs have potential regulatory functions on genes enriched for tissue-specific GO terms. Future work will include epigenetic data from ChIP-seq experiments to further refine these annotations.

Electronic supplementary material

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

Identification and comparison of long non-coding RNAs in the silk gland between domestic and wild silkworms

Under long-term artificial selection, the domestic silkworm (Bombyx mori) has increased its silk yield tremendously in comparison with its wild progenitor, Bombyx mandarina. However, the molecular mechanism of silk yield increase is still unknown. Comparative analysis of long non-coding RNAs (lncRNAs) may provide some insights into understanding this phenotypic variation. In this study, using RNA sequencing technology data of silk gland in domestic and wild silkworms, we identified 599 lncRNAs in the silk gland of the silkworm. Compared with protein-coding genes, the silk gland lncRNA genes tend to have fewer exon numbers, shorter transcript length and lower GC-content. Moreover, we found that three lncRNA genes are significantly and differentially expressed between domestic and wild silkworms. The potential targets of two differentially expressed lncRNAs (DELs) (dw4sg_0040 and dw4sg_0483) and the expression-correlated genes with the two DELs are mainly enriched in the related processes of silk protein translation. This implies that these DELs may affect the phenotypic variation in silk yield between the domestic and wild silkworms through the post-transcriptional regulation of silk protein.

Genome Wide Identification of Novel Long Non-coding RNAs and Their Potential Associations With Milk Proteins in Chinese Holstein Cows

Long non-coding RNAs (lncRNAs) have emerged as a novel class of regulatory molecules involved in various biological processes. However, their role in milk performance is unknown. Here, whole transcriptome RNA sequencing was used to generate the lncRNA transcriptome profiles in mammary tissue samples from 6 Chinese Holstein cows with 3 extremely high and 3 low milk protein percentage phenotypes. In this study, 6,450 lncRNA transcripts were identified through 5 stringent steps and filtration by coding potential. In total, 31 lncRNAs and 18 novel genes were identified to be differentially expressed in high milk protein samples (HP) relative to low milk protein samples (LP), respectively. Differentially expressed lncRNAs were selected to predict target genes through bioinformatics analysis, followed by the integration of differentially expressed mRNA data, gene function, gene ontology (GO) and pathway, genome wide association study (GWAS) and quantitative trait locus (QTL) information, as well as network analysis to further characterize potential interactions. Several lncRNAs were found (such as XLOC_059976) that could be used as candidate markers for milk protein content prediction. This is the first study to perform global expression profiling of lncRNAs and mRNAs related to milk protein traits in dairy cows. These results provide important information and insights into the synthesis of milk proteins, and potential targets for the future improvement of milk quality.