Genome-wide identification and developmental expression profiling of long noncoding RNAs during Drosophila metamorphosis

Expression of De Novo Open Reading Frames in Natural Populations of Drosophila melanogaster

De novo genes, which originate from noncoding DNA, are known to have a high rate of turnover over short evolutionary timescales, such as within a species. Thus, their expression is often lineage- or genetic background-specific. However, little is known about their levels and breadth of expression as populations of a species diverge. In this study, we utilized publicly available RNA-seq data to examine the expression of newly evolved open reading frames (neORFs) in comparison to non- and protein-coding genes in Drosophila melanogaster populations from the derived species range in Europe and the ancestral range in sub-Saharan Africa. Our datasets included two adult tissue types as well as whole bodies at two temperatures for both sexes and three larval/prepupal developmental stages in a single tissue and sex, which allowed us to examine neORF expression and divergence across multiple sample types as well as sex and population. We detected a relatively large proportion (approximately 50%) of annotated neORFs as expressed in the population samples, with neORFs often showing greater expression divergence between populations than non- or protein-coding genes. However, differential expression of neORFs between populations tended to occur in a sample type-specific manner. On the other hand, neORFs displayed less sex-biased expression than the other two gene classes, with the majority of sex-biased neORFs detected in whole bodies, which may be attributable to the presence of the gonads. We also found that neORFs shared among multiple lines in the original set of inbred lines in which they were first detected were more likely to be both expressed and differentially expressed in the new population samples, suggesting that neORFs at a higher frequency (i.e. present in more individuals) within a species are more likely to be functional.

Transcriptome-wide analysis uncovers regulatory elements of the antennal transcriptome repertoire of bumblebee at different life stages

Bumblebees are crucial pollinators, providing essential ecosystem services and global food production. The success of pollination services relies on the interaction between sensory organs and the environment. The antenna functions as a versatile multi-sensory organ, pivotal in mediating chemosensory/olfactory information, and governs adaptive responses to environmental changes. Despite an increasing number of RNA-sequencing studies on insect antenna, comprehensive antennal transcriptome studies at the different life stages were not elucidated systematically. Here, we quantified the expression profile and dynamics of coding/microRNA genes of larval head and antennal tissues from early- and late-stage pupa to the adult of Bombus terrestris as suitable model organism among pollinators. We further performed Pearson correlation analyses on the gene expression profiles of the antennal transcriptome from larval head tissue to adult stages, exploring both positive and negative expression trends. The positively correlated coding genes were primarily enriched in sensory perception of chemical stimuli, ion transport, transmembrane transport processes and olfactory receptor activity. Negatively correlated genes were mainly enriched in organic substance biosynthesis and regulatory mechanisms underlying larval body patterning and the formation of juvenile antennal structures. As post-transcriptional regulators, miR-1000-5p, miR-13b-3p, miR-263-5p and miR-252-5p showed positive correlations, whereas miR-315-5p, miR-92b-3p, miR-137-3p, miR-11-3p and miR-10-3p exhibited negative correlations in antennal tissue. Notably, based on the inverse expression relationship, positively and negatively correlated microRNA (miRNA)-mRNA target pairs revealed that differentially expressed miRNAs predictively targeted genes involved in antennal development, shaping antennal structures and regulating antenna-specific functions. Our data serve as a foundation for understanding stage-specific antennal transcriptomes and large-scale comparative analysis of transcriptomes in different insects.

Long non-coding RNAs involved in Drosophila development and regeneration

The discovery of functional long non-coding RNAs (lncRNAs) changed their initial concept as transcriptional noise. LncRNAs have been identified as regulators of multiple biological processes, including chromatin structure, gene expression, splicing, mRNA degradation, and translation. However, functional studies of lncRNAs are hindered by the usual lack of phenotypes upon deletion or inhibition. Here, we used Drosophila imaginal discs as a model system to identify lncRNAs involved in development and regeneration. We examined a subset of lncRNAs expressed in the wing, leg, and eye disc development. Additionally, we analyzed transcriptomic data from regenerating wing discs to profile the expression pattern of lncRNAs during tissue repair. We focused on the lncRNA CR40469, which is upregulated during regeneration. We generated CR40469 mutant flies that developed normally but showed impaired wing regeneration upon cell death induction. The ability of these mutants to regenerate was restored by the ectopic expression of CR40469. Furthermore, we found that the lncRNA CR34335 has a high degree of sequence similarity with CR40469 and can partially compensate for its function during regeneration in the absence of CR40469. Our findings point to a potential role of the lncRNA CR40469 in trans during the response to damage in the wing imaginal disc.

A Combinatorial Single-Molecule Real-Time and Illumina Sequencing Analysis of Postembryonic Gene Expression in the Asian Citrus Psyllid Diaphorina citri

Huanglongbing (HLB) is a systemic plant disease caused by 'Candidatus Liberibacter asiaticus (CLas)' and transmitted by Diaphorina citri. D. citri acquires the CLas bacteria in the nymph stage and transmits it in the adult stage, indicating that molting from the nymph to adult stages is crucial for HLB transmission. However, the available D. citri reference genomes are incomplete, and gene function studies have been limited to date. In the current research, PacBio single-molecule real-time (SMRT) and Illumina sequencing were performed to investigate the transcriptome of D. citri nymphs and adults. In total, 10,641 full-length, non-redundant transcripts (FLNRTs), 594 alternative splicing (AS) events, 4522 simple sequence repeats (SSRs), 1086 long-coding RNAs (lncRNAs), 281 transcription factors (TFs), and 4459 APA sites were identified. Furthermore, 3746 differentially expressed genes (DEGs) between nymphs and adults were identified, among which 30 DEGs involved in the Hippo signaling pathway were found. Reverse transcription-quantitative PCR (RT-qPCR) further validated the expression levels of 12 DEGs and showed a positive correlation with transcriptome data. Finally, the spatiotemporal expression pattern of genes involved in the Hippo signaling pathway exhibited high expression in the D. citri testis, ovary, and egg. Silencing of the D. citri transcriptional co-activator (DcYki) gene significantly increased D. citri mortality and decreased the cumulative molting. Our results provide useful information and a reliable data resource for gene function research of D. citri.

Regulatory roles of long non-coding RNAs in short-term heat stress in adult worker bees

Long non-coding RNAs (lncRNAs) are crucial modulators of post-transcriptional gene expression regulation, cell fate determination, and disease development. However, lncRNA functions during short-term heat stress in adult worker bees are poorly understood. Here, we performed deep sequencing and bioinformatic analyses of honeybee lncRNAs. RNA interference was performed by using siRNA targeting the most highly expressed lncRNA. The silencing effect on lncRNA and the relative expression levels of seven heat shock protein (HSP) genes, were subsequently examined. Overall, 7,842 lncRNAs and 115 differentially expressed lncRNAs (DELs) were identified in adult worker bees following heat stress exposure. Structural analysis revealed that the overall expression abundance, length of transcripts, exon number, and open reading frames of lncRNAs were lower than those of mRNAs. GO analysis revealed that the target genes were mainly involved in "metabolism," "protein folding," "response to stress," and "signal transduction" pathways. KEGG analysis indicated that the "protein processing in endoplasmic reticulum" and "longevity regulating pathway-multiple species" pathways were most enriched. Quantitative real-time polymerase chain reaction (qRT-PCR) detection of the selected DELs confirmed the reliability of the sequencing data. Moreover, the siRNA experiment indicated that feeding siRNA yielded a silencing efficiency of 77.51% for lncRNA MSTRG.9645.5. Upon silencing this lncRNA, the expression levels of three HSP genes were significantly downregulated (p < 0.05), whereas those of three other HSP genes were significantly upregulated (p < 0.05). Our results provide a new perspective for understanding the regulatory mechanisms of lncRNAs in adult worker bees under short-term heat stress.

Transcriptional changes during isoproterenol-induced cardiac fibrosis in mice

Introduction: β-adrenergic stimulation using β-agonists such as isoproterenol has been routinely used to induce cardiac fibrosis in experimental animal models. Although transcriptome changes in surgical models of cardiac fibrosis such as transverse aortic constriction (TAC) and coronary artery ligation (CAL) are well-studied, transcriptional changes during isoproterenol-induced cardiac fibrosis are not well-explored. Methods: Cardiac fibrosis was induced in male C57BL6 mice by administration of isoproterenol for 4, 8, or 11 days at 50 mg/kg/day dose. Temporal changes in gene expression were studied by RNA sequencing. Results and discussion: We observed a significant alteration in the transcriptome profile across the different experimental groups compared to the saline group. Isoproterenol treatment caused upregulation of genes associated with ECM organization, cell-cell contact, three-dimensional structure, and cell growth, while genes associated with fatty acid oxidation, sarcoplasmic reticulum calcium ion transport, and cardiac muscle contraction are downregulated. A number of known long non-coding RNAs (lncRNAs) and putative novel lncRNAs exhibited differential regulation. In conclusion, our study shows that isoproterenol administration leads to the dysregulation of genes relevant to ECM deposition and cardiac contraction, and serves as an excellent alternate model to the surgical models of heart failure.

Genome-Wide Identification and Involvement in Response to Biotic and Abiotic Stresses of lncRNAs in Turbot (Scophthalmus maximus)

Long non-coding RNAs (lncRNAs) play crucial roles in a variety of biological processes, including stress response. However, the number, characteristics and stress-related expression of lncRNAs in turbot are still largely unknown. In this study, a total of 12,999 lncRNAs were identified at the genome-wide level of turbot for the first time using 24 RNA-seq datasets. Sequence characteristic analyses of transcripts showed that lncRNA transcripts were shorter in average length, lower in average GC content and in average expression level as compared to the coding genes. Expression pattern analyses of lncRNAs in 12 distinct tissues showed that lncRNAs, especially lincRNA, exhibited stronger tissue-specific expression than coding genes. Moreover, 612, 1351, 1060, 875, 420 and 1689 differentially expressed (DE) lncRNAs under Vibrio anguillarum, Enteromyxum scophthalmi, and Megalocytivirus infection and heat, oxygen, and salinity stress conditions were identified, respectively. Among them, 151 and 62 lncRNAs showed differential expression under various abiotic and biotic stresses, respectively, and 11 lncRNAs differentially expressed under both abiotic and biotic stresses were selected as comprehensive stress-responsive lncRNA candidates. Furthermore, expression pattern analysis and qPCR validation both verified the comprehensive stress-responsive functions of these 11 lncRNAs. In addition, 497 significantly co-expressed target genes (correlation coefficient (R) > 0.7 and q-value < 0.05) for these 11 comprehensive stress-responsive lncRNA candidates were identified. Finally, GO and KEGG enrichment analyses indicated that these target genes were enriched mainly in molecular function, such as cytokine activity and active transmembrane transporter activity, in biological processes, such as response to stimulus and immune response, and in pathways, such as protein families: signaling and cellular processes, transporters and metabolism. These findings not only provide valuable reference resources for further research on the molecular basis and function of lncRNAs in turbot but also help to accelerate the progress of molecularly selective breeding of stress-resistant turbot strains or varieties.

Genome-wide analysis of long noncoding RNAs and their association in regulating the metamorphosis of the Sarcophaga peregrina (Diptera: Sarcophagidae)

BACKGROUND

The flesh fly, Sarcophaga peregrina (Diptera: Sarcophagidae), is an important hygiene pest, that causes myiasis in humans and other mammals, typically livestock, and as a vector for various parasitic agents, including bacteria, viruses, and parasites. The role of long non-coding RNAs (lncRNAs) in regulating gene expression during metamorphosis of the flesh fly has not been well established.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we performed genome-wide identification and characterization of lncRNAs from the early pupal stage (1-days pupae), mid-term pupal stage (5-days pupae), and late pupal stage (9-days pupae) of S. peregrina by RNA-seq, and a total of 6921 lncRNAs transcripts were identified. RT-qPCR and enrichment analyses revealed the differentially expressed lncRNAs (DE lncRNAs) that might be associated with insect metamorphosis development. Furthermore, functional analysis revealed that the DE lncRNA (SP_lnc5000) could potentially be involved in regulating the metamorphosis of S. peregrina. RNA interference of SP_lnc5000 caused reduced expression of metamorphosis-related genes in 20-hydroxyecdysone (20E) signaling (Br-c, Ftz-F1), cuticle tanning pathway (TH, DOPA), and chitin related pathway (Cht5). Injection of dsSP_lnc5000 in 3rd instar larvae of S. peregrina resulted in deformed pupae, stagnation of pupal-adult metamorphosis, and a decrease in development time of pupal, pupariation rates and eclosion rates. Hematoxylin-eosin staining (H&E), scanning electron microscope (SEM) observation and cuticle hydrocarbons (CHCs) analysis indicated that SP_lnc5000 had crucial roles in the metamorphosis developmental by modulating pupal cuticular development.

CONCLUSIONS/SIGNIFICANCE

We established that the lncRNA SP_lnc5000 potentially regulates the metamorphosis of S. peregrina by putatively affecting the structure and composition of the pupal cuticle. This study enhances our understanding of lncRNAs as regulators of metamorphosis in S. peregrina, and provide valuable insights into the identification of potential targets for vector control and the development of effective strategies for controlling the spread of myiasis and parasitic diseases.

Interaction Analysis of lncRNA and mRNA Based on the Full-Length Transcriptome of the Nymph-to-Adult Developmental Transition of Sogatella furcifera

Little is known on how long noncoding RNAs (lncRNAs) and mRNAs cooperatively participate in regulating the nymph-to-adult development transition of Sogatella furcifera. Herein, lncRNA and mRNA libraries were constructed in three different developmental stages of S. furcifera, namely, prior to (PE), during (DE), and after (AE) ecdysis. Overall, 4649 lncRNAs were identified and divided into intergenic (53.90%), intronic (1.33%), sense (8.99%), antisense (21.75%), and bidirectional (3.94%) lncRNAs. Moreover, 795 differentially expressed lncRNAs were identified. Specifically, upon comparing PE and DE, 2719 target mRNAs were predicted for 574 lncRNAs. Upon comparing PE and AE, 2816 target mRNAs were predicted for 627 lncRNAs. Finally, upon comparing DE and AE, 51 target mRNAs were predicted for 35 lncRNAs. Kyoto Encyclopedia of Genes and Genome functional enrichment analysis revealed that the target genes of 795 lncRNAs were enriched in metabolic pathways, amino sugar and nucleotide sugar metabolism, and fatty acid metabolism. Subsequently, interaction analysis indicated that MSTRG.16086.1, MSTRG.16087.1, and MSTRG.2447.1 were functionally associated with cuticle protein and chitin biosynthesis. Finally, 11 differentially expressed lncRNAs were significantly enriched in 3rd and 4th instar nymphs. Our findings suggest that lncRNAs play a critical regulatory role during the molting of S. furcifera.

Genome-wide identification of lncRNAs associated with viral infection in Spodoptera frugiperda

The role of lncRNAs in immune defence has been demonstrated in many multicellular and unicellular organisms. However, investigation of the identification and characterization of long non-coding RNAs (lncRNAs) involved in the insect immune response is still limited. In this study, we used RNA sequencing (RNA-seq) to investigate the expression profiles of lncRNAs and mRNAs in the fall armyworm Spodoptera frugiperda in response to virus infection. To assess the tissue- and virus-specificity of lncRNAs, we analysed and compared their expression profiles in haemocytes and fat body of larvae infected with two entomopathogenic viruses with different lifestyles, i.e. the polydnavirus HdIV (Hyposoter didymator IchnoVirus) and the densovirus JcDV (Junonia coenia densovirus). We identified 1883 candidate lncRNAs, of which 529 showed differential expression following viral infection. Expression profiles differed considerably between samples, indicating that many differentially expressed (DE) lncRNAs showed virus- and tissue-specific expression patterns. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and target prediction analyses indicated that DE-LncRNAs were mainly enriched in metabolic process, DNA replication and repair, immune response, metabolism of insect hormone and cell adhesion. In addition, we identified three DE-lncRNAs potentially acting as microRNA host genes, suggesting that they participate in gene regulation by producing miRNAs in response to virus infection. This study provides a catalogue of lncRNAs expressed in two important immune tissues and potential insight into their roles in the antiviral defence in S. frugiperda. The results may help future in-depth functional studies to better understand the biological function of lncRNAs in interaction between viruses and the fall armyworm.

Noncoding RNA Regulation of Hormonal and Metabolic Systems in the Fruit Fly Drosophila

The importance of RNAs is commonly recognised thanks to protein-coding RNAs, whereas non-coding RNAs (ncRNAs) were conventionally regarded as 'junk'. In the last decade, ncRNAs' significance and roles are becoming noticeable in various biological activities, including those in hormonal and metabolic regulation. Among the ncRNAs: microRNA (miRNA) is a small RNA transcript with ~20 nucleotides in length; long non-coding RNA (lncRNA) is an RNA transcript with >200 nucleotides; and circular RNA (circRNA) is derived from back-splicing of pre-mRNA. These ncRNAs can regulate gene expression levels at epigenetic, transcriptional, and post-transcriptional levels through various mechanisms in insects. A better understanding of these crucial regulators is essential to both basic and applied entomology. In this review, we intend to summarise and discuss the current understanding and knowledge of miRNA, lncRNA, and circRNA in the best-studied insect model, the fruit fly Drosophila.

Pleiotropic fitness effects of the lncRNA Uhg4 in Drosophila melanogaster

BACKGROUND

Long noncoding RNAs (lncRNAs) are a diverse class of RNAs that are critical for gene regulation, DNA repair, and splicing, and have been implicated in development, stress response, and cancer. However, the functions of many lncRNAs remain unknown. In Drosophila melanogaster, U snoRNA host gene 4 (Uhg4) encodes an antisense long noncoding RNA that is host to seven small nucleolar RNAs (snoRNAs). Uhg4 is expressed ubiquitously during development and in all adult tissues, with maximal expression in ovaries; however, it has no annotated function(s).

RESULTS

We used CRISPR-Cas9 germline gene editing to generate multiple deletions spanning the promoter region and first exon of Uhg4. Females showed arrested egg development and both males and females were sterile. In addition, Uhg4 deletion mutants showed delayed development and decreased viability, and changes in sleep and responses to stress. Whole-genome RNA sequencing of Uhg4 deletion flies and their controls identified co-regulated genes and genetic interaction networks associated with Uhg4. Gene ontology analyses highlighted a broad spectrum of biological processes, including regulation of transcription and translation, morphogenesis, and stress response.

CONCLUSION

Uhg4 is a lncRNA essential for reproduction with pleiotropic effects on multiple fitness traits.

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.

Systematic identification and characterization of long noncoding RNAs (lncRNAs) during Aedes albopictus development

Background

Aedes albopictus originated in the tropical forests of Southeast Asia and can currently be found on all continents. As one of the main arboviral vectors, the control of Ae. albopictus requires novel strategies, informed by a deep knowledge of its biology. Little is known regarding mosquito long noncoding RNAs (lncRNAs), which are transcripts longer than 200 nucleotides that lack protein-coding potential and have roles in developmental regulation.

Results

Based on RNA-seq data from five developmental time points, eggs, early larvae, late larvae, pupae, and adults (female and male) of Ae. albopictus, 21,414 lncRNAs were characterized in this study. Differential expression analysis revealed that lncRNAs exhibited developmental stage specificity. The expression of most lncRNAs was upregulated at the onset of metamorphosis developmental stages. More differentially expressed lncRNAs were observed between eggs and early larvae. Weighted gene co-expression network analysis (WGCNA) further confirmed that the expression patterns of lncRNAs were obviously correlated with specific developmental time points. Functional annotation using co-expression analysis revealed that lncRNAs may be involved in the regulation of metamorphic developmental transitions of Ae. albopictus. The hub lncRNAs and hub gene clusters were identified for each module that were highly associated with specific developmental time points.

Conclusions

The results of this study will facilitate future researches to elucidate the regulatory mechanisms of lncRNAs in the development of Ae. albopictus and utilize lncRNAs to assist with mosquito control.

The long non-coding RNA landscape of Candida yeast pathogens

Long non-coding RNAs (lncRNAs) constitute a poorly studied class of transcripts with emerging roles in key cellular processes. Despite efforts to characterize lncRNAs across a wide range of species, these molecules remain largely unexplored in most eukaryotic microbes, including yeast pathogens of the Candida clade. Here, we analyze thousands of publicly available sequencing datasets to infer and characterize the lncRNA repertoires of five major Candida pathogens: Candida albicans, Candida tropicalis, Candida parapsilosis, Candida auris and Candida glabrata. Our results indicate that genomes of these species encode hundreds of lncRNAs that show levels of evolutionary constraint intermediate between those of intergenic genomic regions and protein-coding genes. Despite their low sequence conservation across the studied species, some lncRNAs are syntenic and are enriched in shared sequence motifs. We find co-expression of lncRNAs with certain protein-coding transcripts, hinting at potential functional associations. Finally, we identify lncRNAs that are differentially expressed during infection of human epithelial cells for four of the studied species. Our comprehensive bioinformatic analyses of Candida lncRNAs pave the way for future functional characterization of these transcripts.

Long non-coding RNAs (lncRNAs) play roles in key cellular processes, but remain largely unexplored in fungal pathogens such as Candida. Here, Hovhannisyan and Gabaldón analyze thousands of sequencing datasets to infer and characterize the lncRNA repertoires of five Candida species, paving the way for their future functional characterization.

Genome-wide analysis of long non-coding RNAs and their association with wing development in Aphis citricidus (Hemiptera: Aphididae)

Long non-coding RNAs (lncRNAs) play critical roles in the various physiological processes of insects. The wing is a successful adaptation allowing insects to escape from unfavorable environments, while information on lncRNAs related to wing development is limited. In this study, we constructed 12 libraries from two RNA-seq comparisons: 4th instar winged nymphs versus winged adults and 4th instar wingless nymphs versus wingless adults in the brown citrus aphid Aphis citricidus, to identify the wing development-associated lncRNAs. A total of 2914 lncRNAs were identified and 50 lncRNAs were differentially expressed during the 4th instar winged nymphs to winged adults transition, and 28 lncRNAs changed during the 4th instar wingless nymphs to wingless adults transition. The differentially expressed lncRNAs were grouped into six clusters according to the expression patterns in the combined two-winged morphs. lncRNA Ac_lnc54106.1 was up-regulated during 4th instar winged nymphs to winged adults transition, but a lack of change during the 4th instar wingless nymphs to wingless adults transition implied a critical role in the specific regulation of wing development. RNA interference of Ac_lnc54106.1 resulted in malformed wings. Targets prediction, expression patterns, and RNAi assay results showed that Ac_lnc54106.1 may target the PiggyBac transposable element-derived protein 4 (PGBD4) gene, decrease expression of the canonical wing development-related genes, and finally regulate wing development. The systematic identification of lncRNAs in an aphid increases our understanding of how non-coding RNA mediates the wing plasticity of insects.

lnc94638 is a testis-specific long non-coding RNA involved in spermatozoa formation in Zeugodacus cucurbitae (Coquillett)

Long non-coding RNAs (lncRNAs) generally display tissue-specific distributions, and testis-specific lncRNAs form the highest proportion of lncRNAs in many species. Here, we presented a detailed analysis of testis-specific lncRNAs in the melon fly, Zeugodacus cucurbitae, a highly destructive insect pest of cucurbitaceous and other related crops. Most testis-specific lncRNAs were found to be long intergenic non-coding RNAs (lincRNA). The size distribution of these lncRNAs ranged between 600 and 1000 nucleotides. Testis-specific lncRNAs that harboured one isoform number and two exons were the most abundant. Compared to other male tissues, the testis had more highly expressed lncRNAs. The quantitative real-time polymerase chain reaction results of 10 randomly selected testis-specific lncRNAs showed expression patterns consistent with RNA-seq data. Further analysis of the most highly expressed testis-specific lncRNA, lnc94638, was undertaken. Fluorescent in situ hybridization assays localized lnc94638 to the apical region of the testis that contains mature spermatozoa. RNA interference-mediated knockdown of lnc94638 expression reduced spermatozoa numbers and impaired the fertility of Z. cucurbitae male. This study provides a catalogue of testis-specific lncRNAs, shows that the testis-specific lnc94638 is involved in spermatogenesis and has the potential to be used for treating male sterility.

Genome-wide identification and characterization of long non-coding RNAs in Tribolium castaneum

Long non-coding RNAs (lncRNAs) are poorly understood in insects. In this study, we performed genome-wide analysis of lncRNAs in Tribolium castaneum by RNA-seq. In total, 4516 lncRNA transcripts corresponding to 3917 genes were identified from late embryos, early larvae, late larvae, early pupae, late pupae and early adults of T. castaneum, including 3152 novel lncRNAs and 1364 known lncRNAs. These lncRNAs have few exons and transcripts, and are short in length. During development, they exhibited nine different expression patterns. Functionally, they can act either by targeting messenger RNAs (1813 lncRNAs) and lncRNAs (45 lncRNAs) or as micro RNA (miRNA) precursors (46 lncRNAs). LncRNAs were observed to target the metabolic enzymes of glycolysis, TCA cycle and amino acids, demonstrating that lncRNAs control metabolism by regulating metabolic enzymes. Moreover, lncRNAs were shown to participate in cell differentiation and development via their targets. As miRNA precursors, lncRNAs could participate in the ecdysone signaling pathway. This study provides comprehensive information for lncRNAs of T. castaneum, and will promote functional analysis and target identification of lncRNAs in the insect.

Genome-Wide Characterization and Identification of Long Non-Coding RNAs during the Molting Process of a Spider Mite, Panonychus citri

Molting is essential for arthropods to grow. As one of the important arthropod pests in agriculture, key spider mite species (Tetranychus and Panonychus) can normally molt three times from the larva to adult stage within a week. This physiological strategy results in the short lifecycle of spider mites and difficulties in their control in the field. Long non-coding RNAs (lncRNAs) regulate transcriptional editing, cellular function, and biological processes. Thus, analysis of the lncRNAs in the spider mite molting process may provide new insights into their roles in the molting mechanism. For this purpose, we used high-throughput RNA-seq to examine the expression dynamics of lncRNAs and mRNAs in the molting process of different development stages in Panonychus citri. We identified 9199 lncRNAs from 18 transcriptomes. Analysis of the lncRNAs suggested that they were shorter and had fewer exons and transcripts than mRNAs. Among these, 356 lncRNAs were differentially expressed during three molting processes: late larva to early protonymph, late protonymph to early deutonymph, and late deutonymph to early adult. A time series profile analysis of differentially expressed lncRNAs showed that 77 lncRNAs were clustered into two dynamic expression profiles (Pattern a and Pattern c), implying that lncRNAs were involved in the molting process of spider mites. Furthermore, the lncRNA-mRNA co-expression networks showed that several differentially expressed hub lncRNAs were predicted to be functionally associated with typical molting-related proteins, such as cuticle protein and chitin biosynthesis. These data reveal the potential regulatory function of lncRNAs in the molting process and provide datasets for further analysis of lncRNAs and mRNAs in spider mites.

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.

Genome-wide identification of Aedes albopictus long noncoding RNAs and their association with dengue and Zika virus infection

The Asian tiger mosquito, Aedes albopictus (Ae. albopictus), is an important vector that transmits arboviruses such as dengue (DENV), Zika (ZIKV) and Chikungunya virus (CHIKV). Long noncoding RNAs (lncRNAs) are known to regulate various biological processes. Knowledge on Ae. albopictus lncRNAs and their functional role in virus-host interactions are still limited. Here, we identified and characterized the lncRNAs in the genome of an arbovirus vector, Ae. albopictus, and evaluated their potential involvement in DENV and ZIKV infection. We used 148 public datasets, and identified a total of 10, 867 novel lncRNA transcripts, of which 5,809, 4,139, and 919 were intergenic, intronic and antisense respectively. The Ae. albopictus lncRNAs shared many characteristics with other species such as short length, low GC content, and low sequence conservation. RNA-sequencing of Ae. albopictus cells infected with DENV and ZIKV showed that the expression of lncRNAs was altered upon virus infection. Target prediction analysis revealed that Ae. albopictus lncRNAs may regulate the expression of genes involved in immunity and other metabolic and cellular processes. To verify the role of lncRNAs in virus infection, we generated mutations in lncRNA loci using CRISPR-Cas9, and discovered that two lncRNA loci mutations, namely XLOC_029733 (novel lncRNA transcript id: lncRNA_27639.2) and LOC115270134 (known lncRNA transcript id: XR_003899061.1) resulted in enhancement of DENV and ZIKV replication. The results presented here provide an important foundation for future studies of lncRNAs and their relationship with virus infection in Ae. albopictus.

Ae. albopictus is an important vector of arboviruses such as dengue and Zika viruses. Studies on virus-host interaction at gene expression and molecular level are crucial especially in devising methods to inhibit virus replication in Aedes mosquitoes. Previous reports have shown that, besides protein-coding genes, noncoding RNAs such as lncRNAs are also involved in virus-host interaction. In this study, we report a comprehensive catalog of novel lncRNA transcripts in the genome of Ae. albopictus. We also show that the expression of lncRNAs was altered upon infection with dengue and Zika. Additionally, depletion of certain lncRNAs resulted in increased replication of dengue and Zika; hence, suggesting potential association of lncRNAs in virus infection. Results of this study provide a new avenue to the investigation of mosquito-virus interactions, especially in the aspect of noncoding genes.

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.

Genome-wide analysis of long non-coding RNAs in adult tissues of the melon fly, Zeugodacus cucurbitae (Coquillett)

Background

Long non-coding RNAs (lncRNAs) are involved in many fundamental biological processes, such as transcription regulation, protein degradation, and cell differentiation. Information on lncRNA in the melon fly, Zeugodacus cucurbitae (Coquillett) is currently limited.

Results

We constructed 24 RNA-seq libraries from eight tissues (midgut, Malpighian tubules, fat body, ovary, and testis) of Z. cucurbitae adults. A total of 3124 lncRNA transcripts were identified. Among those, 1464 were lincRNAs, 1037 were intronic lncRNAs, 301 were anti-sense lncRNAs, and 322 were sense lncRNAs. The majority of lncRNAs contained two exons and one isoform. Differentially expressed lncRNAs were analyzed between tissues, and Malpighian tubules versus testis had the largest number. Some lncRNAs exhibited strong tissue specificity. Specifically expressed lncRNAs were identified and filtered in tissues of female and male Z. cucurbitae based on their expression levels. Four midgut-specific lncRNAs were validated by quantitative real-time polymerase chain reaction (RT-qPCR), and the data were consistent with RNA-seq data. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of targets of midgut-specific lncRNAs indicated an enrichment of the metabolic process.

Conclusions

This was the first systematic identification of lncRNA in the melon fly. Expressions of lncRNAs in multiple adult tissues were evaluated by quantitative transcriptomic analysis. These qualitative and quantitative analyses of lncRNAs, especially the tissue-specific lncRNAs in Z. cucurbitae, provide useful data for further functional studies.

Coding and Non-coding RNAs: Molecular Basis of Forest-Insect Outbreaks

Insect population dynamics are closely related to 'human' ecological and economic environments, and a central focus of research is outbreaks. However, the lack of molecular-based investigations restricts our understanding of the intrinsic mechanisms responsible for insect outbreaks. In this context, the moth Dendrolimus punctatus Walker can serve as an ideal model species for insect population dynamics research because it undergoes periodic outbreaks. Here, high-throughput whole-transcriptome sequencing was performed using D. punctatus, sampled during latent and outbreak periods, to systemically explore the molecular basis of insect outbreaks and to identify the involved non-coding RNA (ncRNA) regulators, namely microRNAs, long non-coding RNAs, and circular RNAs. Differentially expressed mRNAs of D. punctatus from different outbreak periods were involved in developmental, reproductive, immune, and chemosensory processes; results that were consistent with the physiological differences in D. punctatus during differing outbreak periods. Targets analysis of the non-coding RNAs indicated that long non-coding RNAs could be the primary ncRNA regulators of D. punctatus outbreaks, while circular RNAs mainly regulated synapses and cell junctions. The target genes of differentially expressed microRNAs mainly regulated the metabolic and reproductive pathways during the D. punctatus outbreaks. Developmental, multi-organismal, and reproductive processes, as well as biological adhesion, characterized the competing endogenous RNA network. Chemosensory and immune genes closely related to the outbreak of D. punctatus were further analyzed in detail: from their ncRNA regulators' analysis, we deduce that both lncRNA and miRNA may play significant roles. This is the first report to examine the molecular basis of coding and non-coding RNAs' roles in insect outbreaks. The results provide potential biomarkers for control targets in forest insect management, as well as fresh insights into underlying outbreak-related mechanisms, which could be used for improving insect control strategies in the future.

Long noncoding RNA PAHAL modulates locust behavioural plasticity through the feedback regulation of dopamine biosynthesis

Some long noncoding RNAs (lncRNAs) are specifically expressed in brain cells, implying their neural and behavioural functions. However, how lncRNAs contribute to neural regulatory networks governing the precise behaviour of animals is less explored. Here, we report the regulatory mechanism of the nuclear-enriched lncRNA PAHAL for dopamine biosynthesis and behavioural adjustment in migratory locusts (Locusta migratoria), a species with extreme behavioral plasticity. PAHAL is transcribed from the sense (coding) strand of the gene encoding phenylalanine hydroxylase (PAH), which is responsible for the synthesis of dopamine from phenylalanine. PAHAL positively regulates PAH expression resulting in dopamine production in the brain. In addition, PAHAL modulates locust behavioral aggregation in a population density-dependent manner. Mechanistically, PAHAL mediates PAH transcriptional activation by recruiting serine/arginine-rich splicing factor 2 (SRSF2), a transcription/splicing factor, to the PAH proximal promoter. The co-activation effect of PAHAL requires the interaction of the PAHAL/SRSF2 complex with the promoter-associated nascent RNA of PAH. Thus, the data support a model of feedback modulation of animal behavioural plasticity by an lncRNA. In this model, the lncRNA mediates neurotransmitter metabolism through orchestrating a local transcriptional loop.

The neurotransmitter dopamine is crucial for the neuronal and behavioral response in animals. Phenylalanine hydroxylase (PAH) is involved in dopamine biosynthesis and behavioral regulation in the migratory locust. However, the molecular mechanism for the fine tuning of PAH expression in behavioral response remains ambiguous. Here we discovered a nuclear-enriched lncRNA PAHAL that is transcribed from the coding strand of the PAH gene in the locust (i.e., sense lncRNA). PAHAL positively regulated PAH expression and dopamine production in the brain. In addition, PAHAL modulated behavioral aggregation of the locust. Mechanistically, PAHAL mediated the transcriptional activation of PAH by recruiting SRSF2, a transcription/splicing factor, to the promoter-associated nascent RNA of PAH. These data support a model of feedback modulation of dopamine biosynthesis and behavioral plasticity via a sense lncRNA in the catecholamine metabolic pathway.

Insights into the Functions of LncRNAs in Drosophila

Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs longer than 200 nucleotides (nt). LncRNAs have high spatiotemporal specificity, and secondary structures have been preserved throughout evolution. They have been implicated in a range of biological processes and diseases and are emerging as key regulators of gene expression at the epigenetic, transcriptional, and post-transcriptional levels. Comparative analyses of lncRNA functions among multiple organisms have suggested that some of their mechanisms seem to be conserved. Transcriptome studies have found that some Drosophila lncRNAs have highly specific expression patterns in embryos, nerves, and gonads. In vivo studies of lncRNAs have revealed that dysregulated expression of lncRNAs in Drosophila may result in impaired embryo development, impaired neurological and gonadal functions, and poor stress resistance. In this review, we summarize the epigenetic, transcriptional, and post-transcriptional mechanisms of lncRNAs and mainly focus on recent insights into the transcriptome studies and biological functions of lncRNAs in Drosophila.

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.

Re-thinking adaptive immunity in the beetles: Evolutionary and functional trajectories of lncRNAs

Unlike vertebrate animals, invertebrates lack lymphocytes and therefore have historically been believed not to develop immune memory. A few studies have reported evidence of immune priming in insects; however, these studies lack the molecular mechanism and proposed it might be different among taxa. Since lncRNAs are known to regulate the immune response, we identified 10,120 lncRNAs in Tribolium castaneum genome-wide followed by transcriptome analysis of primed and unprimed larvae of different infectious status. A shift in lncRNA expression between Btt primed larvae and other treatment groups provides evidence of immune memory response. A few "priming" lncRNAs (n = 9) were uniquely regulated in Btt primed larvae. Evidence suggests these lncRNAs are likely controlling immune priming in Tribolium by regulating expression of genes involved in proteasomal machinery, Notch system, zinc metabolism, and methyltransferase activity, which are necessary to modulate phagocytosis. Our results support a conserved immune priming mechanism in a macrophage-dependent manner.

Genome-wide identification and characterization of long intergenic noncoding RNAs in the regenerative flatworm Macrostomum lignano

The free-living flatworm Macrostoma lignano (M. lignano) is an emerging model organism for aging and regeneration research. Long intergenic non-coding RNAs (lincRNAs) have important roles in many biological processes such as aging, stem cell maintenance and differentiation. However, to date, there is no systematic identification of lincRNAs in M. lignano. By using public RNA-seq data, we identified a total of 2547 lincRNA transcripts in M. lignano genome. We discovered that M. lignano lincRNAs shared many characteristics with other species such as shorter in length, lower GC content, and lower in expression compared to protein-coding genes. Unlike protein-coding genes, M. lignano lincRNAs showed higher tendency to be expressed in temporal and region-specific fashion. Additionally, co-expression network analysis and functional enrichment suggest that M. lignano lincRNAs have potential roles in regeneration. This study will provide important resources and pave the way for investigations on non-coding genes involved in aging and regeneration.

lncRNA profile of Apis mellifera and its possible role in behavioural transition from nurses to foragers

BACKGROUND

The behavioural transition from nurses to foragers in honey bees is known to be affected by intrinsic and extrinsic factors, including colony demography, hormone levels, brain chemistry and structure, and gene expression in the brain. However, the molecular mechanism underlying this behavioural transition of honey bees is still obscure.

RESULTS

Through RNA sequencing, we performed a comprehensive analysis of lncRNAs and mRNAs in honey bee nurses and foragers. Nurses and foragers from both typical colonies and single-cohort colonies were used to prepare six libraries to generate 49 to 100 million clear reads per sample. We obtained 6863 novel lncRNAs, 1480 differentially expressed lncRNAs between nurses and foragers, and 9308 mRNAs. Consistent with previous studies, lncRNAs showed features distinct from mRNAs, such as shorter lengths, lower exon numbers, and lower expression levels compared to mRNAs. Bioinformatic analysis showed that differentially expressed genes were mostly involved in the regulation of sensory-related events, such as olfactory receptor activity and odorant binding, and enriched Wnt and FoxO signaling pathways. Moreover, we found that lncRNAs TCONS_00356023, TCONS_00357367, TCONS_00159909 and mRNAs dop1, Kr-h1 and HR38 may play important roles in behavioural transition in honey bees.

CONCLUSION

This study characterized the expression profile of lncRNAs in nurses and foragers and provided a framework for further study of the role of lncRNAs in honey bee behavioural transition.

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.

Point centromere activity requires an optimal level of centromeric noncoding RNA

Budding yeast harbors a simple point centromere, which is originally believed to be sequence dependent without much epigenetic regulation and is transcription incompatible, as inserting a strong promoter upstream inactivates the centromere completely. Here, we demonstrate that an optimal level centromeric noncoding RNA is required for budding yeast centromere activity. Centromeric transcription is induced in S phase, coinciding with the assembly of new centromeric proteins. Too much or too little centromeric noncoding RNA leads to centromere malfunction. Overexpression of centromeric noncoding RNA reduces the protein levels and chromatin localization of inner centromere and kinetochore proteins, such as CENP-A, CENP-C, and the chromosome passenger complex. This work shows that point centromere is epigenetically regulated by noncoding RNA.

In budding yeast, which possesses simple point centromeres, we discovered that all of its centromeres express long noncoding RNAs (cenRNAs), especially in S phase. Induction of cenRNAs coincides with CENP-A^Cse4 loading time and is dependent on DNA replication. Centromeric transcription is repressed by centromere-binding factor Cbf1 and histone H2A variant H2A.Z^Htz1. Deletion of CBF1 and H2A.Z^HTZ1 results in an up-regulation of cenRNAs; an increased loss of a minichromosome; elevated aneuploidy; a down-regulation of the protein levels of centromeric proteins CENP-A^Cse4, CENP-A chaperone HJURP^Scm3, CENP-C^Mif2, Survivin^Bir1, and INCENP^Sli15; and a reduced chromatin localization of CENP-A^Cse4, CENP-C^Mif2, and Aurora B^Ipl1. When the RNA interference system was introduced to knock down all cenRNAs from the endogenous chromosomes, but not the cenRNA from the circular minichromosome, an increase in minichromosome loss was still observed, suggesting that cenRNA functions in trans to regulate centromere activity. CenRNA knockdown partially alleviates minichromosome loss in cbf1Δ, htz1Δ, and cbf1Δ htz1Δ in a dose-dependent manner, demonstrating that cenRNA level is tightly regulated to epigenetically control point centromere function.

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.

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.

The lncRNA male-specific abdominal plays a critical role in Drosophila accessory gland development and male fertility

Although thousands of long non-coding RNAs (lncRNA) have been identified in the genomes of higher eukaryotes, the precise function of most of them is still unclear. Here, we show that a >65 kb, male-specific, lncRNA, called male-specific abdominal (msa) is required for the development of the secondary cells of the Drosophila male accessory gland (AG). msa is transcribed from within the Drosophila bithorax complex and shares much of its sequence with another lncRNA, the iab-8 lncRNA, which is involved in the development of the central nervous system (CNS). Both lncRNAs perform much of their functions via a shared miRNA embedded within their sequences. Loss of msa, or of the miRNA it contains, causes defects in secondary cell morphology and reduces male fertility. Although both lncRNAs express the same miRNA, the phenotype in the secondary cells and the CNS seem to reflect misregulation of different targets in the two tissues.

In many animals, the male seminal fluid induces physiology changes in the mated female that increase a male’s reproductive success. These changes are often referred to as the post-mating response (PMR). In Drosophila, the seminal fluid proteins responsible for generating the PMR are made in a specialized gland, analogous to the mammalian seminal vesicle and prostate, called the accessory gland (AG). In this work, we show that a male-specific, long, non-coding RNA (lncRNA), called msa, plays a critical role in the development and function of this gland, primarily through a microRNA (miRNA) encoded within its sequence. This same miRNA had previously been shown to be expressed in the central nervous system (CNS) via an alternative promoter, where its ability to repress homeotic genes is required for both male and female fertility. Here, we present evidence that the targets of this miRNA in the AG are likely different from those found in the CNS. Thus, the same miRNA seems to have been selected to affect Drosophila fertility through two different mechanisms. Although many non-coding RNAs have now been identified, very few can be shown to have function. Our work highlights a lncRNA that has multiple biological functions, affecting cellular morphology and fertility.

Long non-coding regulatory RNAs in sponges and insights into the origin of animal multicellularity

How animals evolved from a single-celled ancestor over 700 million years ago is poorly understood. Recent transcriptomic and chromatin analyses in the sponge Amphimedon queenslandica, a morphologically-simple representative of one of the oldest animal phyletic lineages, have shed light on what innovations in the genome and its regulation underlie the emergence of animal multicellularity. Comparisons of the regulatory genome of this sponge with those of more complex bilaterian model species and even simpler unicellular relatives have revealed that fundamental changes in genome regulatory complexity accompanied the evolution of animal multicellularity. Here, we review and discuss the results of these recent investigations by specifically focusing on the contribution of long non-coding RNAs to the evolution of the animal regulatory genome.

Sponge Long Non-Coding RNAs Are Expressed in Specific Cell Types and Conserved Networks

Although developmental regulation by long non-coding RNAs (lncRNAs) appears to be a widespread feature amongst animals, the origin and level of evolutionary conservation of this mode of regulation remain unclear. We have previously demonstrated that the sponge Amphimedon queenslandica—a morphologically-simple animal—developmentally expresses an array of lncRNAs in manner akin to more complex bilaterians (insects + vertebrates). Here, we first show that Amphimedon lncRNAs are expressed in specific cell types in larvae, juveniles and adults. Thus, as in bilaterians, sponge developmental regulation involves the dynamic, cell type- and context-specific regulation of specific lncRNAs. Second, by comparing gene co-expression networks between Amphimedon queenslandica and Sycon ciliatum—a distantly-related calcisponge—we identify several putative co-expression modules that appear to be shared in sponges; these network-embedded sponge lncRNAs have no discernable sequence similarity. Together, these results suggest sponge lncRNAs are developmentally regulated and operate in conserved gene regulatory networks, as appears to be the case in more complex bilaterians.

Integration of lncRNA and mRNA Transcriptome Analyses Reveals Genes and Pathways Potentially Involved in Calf Intestinal Growth and Development during the Early Weeks of Life

A better understanding of the factors that regulate growth and immune response of the gastrointestinal tract (GIT) of calves will promote informed management practices in calf rearing. This study aimed to explore genomics (messenger RNA (mRNA)) and epigenomics (long non-coding RNA (lncRNA)) mechanisms regulating the development of the rumen and ileum in calves. Thirty-two calves (≈5-days-old) were reared for 96 days following standard procedures. Sixteen calves were humanely euthanized on experiment day 33 (D33) (pre-weaning) and another 16 on D96 (post-weaning) for collection of ileum and rumen tissues. RNA from tissues was subjected to next generation sequencing and 3310 and 4217 mRNAs were differentially expressed (DE) between D33 and D96 in ileum and rumen tissues, respectively. Gene ontology and pathways enrichment of DE genes confirmed their roles in developmental processes, immunity and lipid metabolism. A total of 1568 (63 known and 1505 novel) and 4243 (88 known and 4155 novel) lncRNAs were detected in ileum and rumen tissues, respectively. Cis target gene analysis identified BMPR1A, an important gene for a GIT disease (juvenile polyposis syndrome) in humans, as a candidate cis target gene for lncRNAs in both tissues. LncRNA cis target gene enrichment suggested that lncRNAs might regulate growth and development in both tissues as well as posttranscriptional gene silencing by RNA or microRNA processing in rumen, or disease resistance mechanisms in ileum. This study provides a catalog of bovine lncRNAs and set a baseline for exploring their functions in calf GIT development.

Investigating the genetic and epigenetic basis of big biological questions with the parthenogenetic marbled crayfish: A review and perspectives

In the last 15 years, considerable attempts have been undertaken to develop the obligately parthenogenetic marbled crayfish Procambarus virginalis as a new model in biology. Its main advantage is the production of large numbers of offspring that are genetically identical to the mother, making this crustacean particularly suitable for research in epigenetics. Now, a draft genome, transcriptome and genome-wide methylome are available opening new windows for research. In this article, I summarize the biological advantages and genomic and epigenetic features of marbled crayfish and, based on first promising data, discuss what this new model could contribute to answering of ''big'' biological questions. Genome mining is expected to reveal new insights into the genetic specificities of decapod crustaceans, the genetic basis of arthropod reproduction, moulting and immunity, and more general topics such as the genetic underpinning of adaptation to fresh water, omnivory, biomineralization, sexual system change, behavioural variation, clonal genome evolution, and resistance to cancer. Epigenetic investigations with the marbled crayfish can help clarifying the role of epigenetic mechanisms in gene regulation, tissue specification, adult stem cell regulation, cell ageing, organ regeneration and disease susceptibility. Marbled crayfish is further suitable to elucidate the relationship between genetic and epigenetic variation, the transgenerational inheritance of epigenetic signatures and the contribution of epigenetic phenotype variation to the establishment of social hierarchies, environmental adaptation and speciation. These issues can be tackled by experiments with highly standardized laboratory lineages, comparison of differently adapted wild populations and the generation of genetically and epigenetically edited strains.

From Heterochromatin to Long Noncoding RNAs in Drosophila: Expanding the Arena of Gene Function and Regulation

Recent years have witnessed a remarkable interest in exploring the significance of pervasive noncoding transcripts in diverse eukaryotes. Classical cytogenetic studies using the Drosophila model system unraveled the perplexing attributes and "functions" of the "gene"-poor heterochromatin. Recent molecular studies in the fly model are likewise revealing the very diverse and significant roles played by long noncoding RNAs (lncRNAs) in development, gene regulation, chromatin organization, cell and nuclear architecture, etc. There has been a rapid increase in the number of identified lncRNAs, although a much larger number still remains unknown. The diversity of modes of actions and functions of the limited number of Drosophila lncRNAs, which have been examined, already reflects the profound roles of such RNAs in generating and sustaining the biological complexities of eukaryotes. Several of the known Drosophila lncRNAs originate as independent sense or antisense transcripts from promoter or intergenic, intronic, or 5'/3'-UTR regions, while many of them are independent genes that produce only lncRNAs or coding as well as noncoding RNAs. The different lncRNAs affect chromatin organization (local or large-scale pan-chromosomal), transcription, RNA processing/stability, or translation either directly through interaction with their target DNA sequences or indirectly by acting as intermediary molecules for specific regulatory proteins or may act as decoys/sinks, or storage sites for specific proteins or groups of proteins, or may provide a structural framework for the assembly of substructures in nucleus/cytoplasm. It is interesting that many of the "functions" alluded to heterochromatin in earlier cytogenetic studies appear to find correlates with the known subtle as well as far-reaching actions of the different small and long noncoding RNAs. Further studies exploiting the very rich and powerful genetic and molecular resources available for the Drosophila model are expected to unravel the mystery underlying the long reach of ncRNAs.

Genome-wide identification of long non-coding RNA genes and their association with insecticide resistance and metamorphosis in diamondback moth, Plutella xylostella

Long non-coding RNA (lncRNA) is a class of noncoding RNA >200 bp in length that has essential roles in regulating a variety of biological processes. Here, we constructed a computational pipeline to identify lncRNA genes in the diamondback moth (Plutella xylostella), a major insect pest of cruciferous vegetables. In total, 3,324 lncRNAs corresponding to 2,475 loci were identified from 13 RNA-Seq datasets, including samples from parasitized, insecticide-resistant strains and different developmental stages. The identified P. xylostella lncRNAs had shorter transcripts and fewer exons than protein-coding genes. Seven out of nine randomly selected lncRNAs were validated by strand-specific RT-PCR. In total, 54-172 lncRNAs were specifically expressed in the insecticide resistant strains, among which one lncRNA was located adjacent to the sodium channel gene. In addition, 63-135 lncRNAs were specifically expressed in different developmental stages, among which three lncRNAs overlapped or were located adjacent to the metamorphosis-associated genes. These lncRNAs were either strongly or weakly co-expressed with their overlapping or neighboring mRNA genes. In summary, we identified thousands of lncRNAs and presented evidence that lncRNAs might have key roles in conferring insecticide resistance and regulating the metamorphosis development in P. xylostella.

Spatiotemporal expression profiling of long intervening noncoding RNAs in Caenorhabditis elegans

To better understand the biological function of long noncoding RNAs, it is critical to determine their spatiotemporal expression patterns. We generated transgenic reporter strains for 149 out of the 170 annotated C. elegans long intervening noncoding RNAs (lincRNAs) and profiled their temporal activity. For the 68 lincRNAs with integrated reporter lines, we profiled their expression at the resolution of single cells in L1 larvae, and revealed that the expression of lincRNAs is more specific, heterogeneous and at lower level than transcription factors (TFs). These expression patterns can be largely attributed to transcriptional regulation because they were observed in assays using reporters of promoter activity. The spatial expression patterns of the 68 lincRNAs were further examined in 18 tissue categories throughout eight developmental stages. We compared the expression dynamics of lincRNAs, miRNAs and TFs during development. lincRNA and miRNA promoters are less active at embryo stage than those of TFs, but become comparable to TFs after embryogenesis. Finally, the lincRNA gene set shows a similar tissue distribution to that of miRNAs and TFs. We also generated a database, CELE, for the storage and retrieval of lincRNA reporter expression patterns and other relevant information. The data and strains described here will provide a valuable guide and resource for future functional exploration of C. elegans lincRNAs.

Genome-wide identification of lncRNAs associated with chlorantraniliprole resistance in diamondback moth Plutella xylostella (L.)

BACKGROUND

Long noncoding RNAs (lncRNAs) are now considered important regulatory factors, with a variety of biological functions in many species including insects. Some lncRNAs have the ability to show rapid responses to diverse stimuli or stress factors and are involved in responses to insecticide. However, there are no reports to date on the characterization of lncRNAs associated with chlorantraniliprole resistance in Plutella xylostella.

RESULTS

Nine RNA libraries constructed from one susceptible (CHS) and two chlorantraniliprole-resistant P. xylostella strains (CHR, ZZ) were sequenced, and 1309 lncRNAs were identified, including 877 intergenic lncRNAs, 190 intronic lncRNAs, 76 anti-sense lncRNAs and 166 sense-overlapping lncRNAs. Of the identified lncRNAs, 1059 were novel. Furthermore, we found that 64 lncRNAs were differentially expressed between CHR and CHS and 83 were differentially expressed between ZZ and CHS, of which 22 were differentially expressed in both CHR and ZZ. Most of the differentially expressed lncRNAs were hypothesized to be associated with chlorantraniliprole resistance in P. xylostella. The targets of lncRNAs via cis- (<10 kb upstream and downstream) or trans- (Pearson's correlation, r > 0.9 or < -0.9, P < 0.05) regulatory effects were also identified; many of the differently expressed lncRNAs were correlated with various important protein-coding genes involved in insecticide resistance, such as the ryanodine receptor, uridine diphosphate glucuronosyltransferase (UGTs), cytochrome P450, esterase and the ATP-binding cassette transporter.

CONCLUSIONS

Our results represent the first global identification of lncRNAs associated with chlorantraniliprole resistance in P. xylostella. These results will facilitate future studies of the regulatory mechanisms of lncRNAs in chlorantraniliprole and other insecticide resistance and in other biological processes in P. xylostella.

Origin and evolution of the metazoan non-coding regulatory genome

Animals rely on genomic regulatory systems to direct the dynamic spatiotemporal and cell-type specific gene expression that is essential for the development and maintenance of a multicellular lifestyle. Although it is widely appreciated that these systems ultimately evolved from genomic regulatory mechanisms present in single-celled stem metazoans, it remains unclear how this occurred. Here, we focus on the contribution of the non-coding portion of the genome to the evolution of animal gene regulation, specifically on recent insights from non-bilaterian metazoan lineages, and unicellular and colonial holozoan sister taxa. High-throughput next-generation sequencing, largely in bilaterian model species, has led to the discovery of tens of thousands of non-coding RNA genes (ncRNAs), including short, long and circular forms, and uncovered the central roles they play in development. Based on the analysis of non-bilaterian metazoan, unicellular holozoan and fungal genomes, the evolution of some ncRNAs, such as Piwi-interacting RNAs, correlates with the emergence of metazoan multicellularity, while others, including microRNAs, long non-coding RNAs and circular RNAs, appear to be more ancient. Analysis of non-coding regulatory DNA and histone post-translational modifications have revealed that some cis-regulatory mechanisms, such as those associated with proximal promoters, are present in non-animal holozoans, while others appear to be metazoan innovations, most notably distal enhancers. In contrast, the cohesin-CTCF system for regulating higher-order chromatin structure and enhancer-promoter long-range interactions appears to be restricted to bilaterians. Taken together, most bilaterian non-coding regulatory mechanisms appear to have originated before the divergence of crown metazoans. However, differential expansion of non-coding RNA and cis-regulatory DNA repertoires in bilaterians may account for their increased regulatory and morphological complexity relative to non-bilaterians.

BmncRNAdb: a comprehensive database of non-coding RNAs in the silkworm, Bombyx mori

BACKGROUND

Long non-coding RNAs (lncRNAs) may play critical roles in a wide range of developmental processes of higher organisms. Recently, lncRNAs have been widely identified across eukaryotes and many databases of lncRNAs have been developed for human, mouse, fruit fly, etc. However, there is rare information about them in the only completely domesticated insect, silkworm (Bombyx mori).

DESCRIPTION

In this study, we systematically scanned lncRNAs using the available silkworm RNA-seq data and public unigenes. Finally, we identified and collected 6281 lncRNAs in the silkworm. Besides, we also collected 1986 microRNAs (miRNAs) from previous studies. Then, we organized them into a comprehensive and web-based database, BmncRNAdb. This database offers a user-friendly interface for data browse and online analysis as well as the three online tools for users to predict the target genes of lncRNA or miRNA.

CONCLUSIONS

We have systematically identified and collected the silkworm lncRNAs and constructed a comprehensive database of the silkworm lncRNAs and miRNAs. This work gives a glimpse into lncRNAs of the silkworm and lays foundations for the ncRNAs study of the silkworm and other insects in the future. The BmncRNAdb is freely available at http://gene.cqu.edu.cn/BmncRNAdb/index.php .

Developmentally regulated long non-coding RNAs in Xenopus tropicalis

Advances in RNA sequencing technologies have led to the surprising discovery that a vast number of transcripts emanate from regions of the genome that are not part of coding genes. Although some of the smaller ncRNAs such as microRNAs have well-characterized functions, the majority of long ncRNA (lncRNA) functions remain poorly understood. Understanding the significance of lncRNAs is an important challenge facing biology today. A powerful approach to uncovering the function of lncRNAs is to explore temporal and spatial expression profiling. This may be particularly useful for classes of lncRNAs that have developmentally important roles as the expression of such lncRNAs will be expected to be both spatially and temporally regulated during development. Here, we take advantage of our ultra-high frequency (temporal) sampling of Xenopus embryos to analyze gene expression trajectories of lncRNA transcripts over the first 3 days of development. We computationally identify 5689 potential single- and multi-exon lncRNAs. These lncRNAs demonstrate clear dynamic expression patterns. A subset of them displays highly correlative temporal expression profiles with respect to those of the neighboring genes. We also identified spatially localized lncRNAs in the gastrula stage embryo. These results suggest that lncRNAs have regulatory roles during early embryonic development.