Identification and evolutionary analysis of long non-coding RNAs in zebra finch

Genome-wide transcriptome analysis reveals the diversity and function of long non-coding RNAs in dinoflagellates

Dinoflagellates are a diverse group of phytoplankton, ranging from harmful bloom-forming microalgae to photosymbionts of coral reefs. Genome-scale data from dinoflagellates reveal atypical genomic features, extensive genomic divergence, and lineage-specific innovation of gene functions. Long non-coding RNAs (lncRNAs), known to regulate gene expression in eukaryotes, are largely unexplored in dinoflagellates. Here, using high-quality genome and transcriptome data, we identified 48039 polyadenylated lncRNAs in three dinoflagellate species: the coral symbionts Cladocopium proliferum and Durusdinium trenchii, and the bloom-forming species, Prorocentrum cordatum. These lncRNAs have fewer introns and lower G+C content than protein-coding sequences; 37 768 (78.6%) are unique with respect to sequence similarity. We classified all lncRNAs based on conserved motifs (k-mers) into distinct clusters, following properties of protein-binding and/or subcellular localisation. Interestingly, 3708 (7.7%) lncRNAs are differentially expressed under heat stress, algal lifestyle, and/or growth phase, and share co-expression patterns with protein-coding genes. Based on inferred triplex interactions between lncRNA and putative promoter regions, we identified 19 460 putative gene targets for 3721 lncRNAs; 907 genes exhibit differential expression under heat stress. These results reveal, for the first time, the diversity of lncRNAs in dinoflagellates and how lncRNAs may regulate gene expression as a heat-stress response in these ecologically important microbes.

Long non-coding RNA profile in banana shrimp, Fenneropenaeus merguiensis and the potential role of lncPV13 in vitellogenesis

The long non-coding RNAs (lncRNAs) have been known to play important roles in several biological processes as well as in reproduction. This study aimed to identify lncRNA in ovary female banana shrimp, Fenneropenaeus merguiensis, and investigate the potential role of lncPV13 in the vitellogenesis. After the in silico identification of the ovarian transcriptome, a total of 24,733 putative lncRNAs were obtained, and only 147 putative lncRNAs were significantly differentially expressed among the ovarian development stages. To validate the in silico identification of lncRNAs, the 16 lncRNAs with the highest differential expression in the transcriptome analysis were evaluated by RT-qPCR. The 6 lncRNAs showed higher expression levels in the mature stage than in the previtellogenic stage and were found in several tissues such as in eyestalks, brains, thoracic ganglia, gills, and muscle. Furthermore, most candidate lncRNAs were amplifiable in Litopenaeus vannamei's and Penaeus monodon's DNA but not in Macrobrachium rosenbergii's DNA, suggesting some lncRNAs are expressed in a species-specific manner among penaeid shrimp. In this study, the lncPV13 was investigated for its vitellogenin regulating function by RNA interference. The result indicates that the lncPV13 expression was suppressed in the ovary on day 7 after the injection of double-stranded RNA specific to lncPV13 (dslncPV13), while vitellogenin (Vg) expression was significantly decreased. In contrast, the gonad inhibiting hormone (GIH) expression was significantly increased in the lncPV13 knockdown shrimp. However, the oocyte proliferation was not significantly different between control and lncPV13 knockdown shrimp. This suggests that lncPV13 regulate Vg synthesis through GIH inhibition. Finally, our findings provide lncRNA information and potential lncRNAs involved in the vitellogenesis of female banana shrimp.

In Silico Analysis of Common Long Noncoding RNAs in Schistosoma mansoni and Schistosoma haematobium

Background

Schistosomiasis caused by Schistosoma parasites is one of the most common parasitic infections worldwide. Genetic regulation of the genus Schistosoma, which has different developmental stages throughout its life, is quite complex. In these parasites, thousands of long noncoding RNAs (lncRNAs) estimated to be functional were identified. Identifying the transcripts expressed in common and detecting their functions for better understanding of the role of these lncRNAs require a comparative study.

Methods

Assembled RNA-seq datasets belonging to S. mansoni and S. haematobium were obtained from the National Center for Biotechnology. A basic local alignment search tool (BLASTN) analysis was conducted against previously constructed lncRNA library to identify the common lncRNAs between two species. LncRNAs target genes and their gene ontology annotation was performed.

Results

In S. mansoni and S. haematobium, 5132 and 3589 lncRNA transcripts were detected, respectively. These two species had 694 lncRNAs in common. A significant number of lncRNAs was determined to be transcribed from sex chromosomes. The frequently expressed lncRNAs appear to be involved in metabolic and biological regulation processes.

Conclusions

These two species share similar lncRNAs; thus, this finding is a clue that they might have similar functions. In sexual development, they especially might play important roles. Our results will provide important clues to further studies about interactions between human hosts and parasites and the infection mechanisms of Schistosoma parasites.

A systematic evaluation of bioinformatics tools for identification of long noncoding RNAs

High-throughput RNA sequencing unveiled the complexity of transcriptome and significantly increased the records of long noncoding RNAs (lncRNAs), which were reported to participate in a variety of biological processes. Identification of lncRNAs is a key step in lncRNA analysis, and a bunch of bioinformatics tools have been developed for this purpose in recent years. While these tools allow us to identify lncRNA more efficiently and accurately, they may produce inconsistent results, making selection a confusing issue. We compared the performance of 41 analysis models based on 14 software packages and different data sets, including high-quality data and low-quality data from 33 species. In addition, computational efficiency, robustness, and joint prediction of the models were explored. As a practical guidance, key points for lncRNA identification under different situations were summarized. In this investigation, no one of these models could be superior to others under all test conditions. The performance of a model relied to a great extent on the source of transcripts and the quality of assemblies. As general references, FEELnc_all_cl, CPC, and CPAT_mouse work well in most species while COME, CNCI, and lncScore are good choices for model organisms. Since these tools are sensitive to different factors such as the species involved and the quality of assembly, researchers must carefully select the appropriate tool based on the actual data. Alternatively, our test suggests that joint prediction could behave better than any single model if proper models were chosen. All scripts/data used in this research can be accessed at http://bioinfo.ihb.ac.cn/elit.

Feather Evolution from Precocial to Altricial Birds

Birds are the most abundant terrestrial vertebrates and their diversity is greatly shaped by the feathers. How avian evolution is linked to feather evolution has long been a fascinating question. Numerous excellent studies have shed light on this complex relationship by investigating feather diversity and its underlying molecular mechanisms. However, most have focused on adult domestic birds, and the contribution of feather diversity to environmental adaptation has not been well-studied. In this review, we described bird diversity using the traditional concept of the altricial-precocial spectrum in bird hatchlings. We combined the spectrum with a recently published avian phylogeny to profile the spectrum evolution. We then focused on the discrete diagnostic character of the spectrum, the natal down, and propose a hypothesis for the precocial-to-altricial evolution. For the underlying molecular mechanisms in feather diversity and bird evolution, we reviewed the literature and constructed the known mechanisms for feather tract definition and natal down development. Finally, we suggested some future directions for research on altricial-precocial divergence, which may expand our understanding of the relationship between natal down diversity and bird evolution.

Re-analysis of the coral Acropora digitifera transcriptome reveals a complex lncRNAs-mRNAs interaction network implicated in Symbiodinium infection

BACKGROUND

Being critically important to the ecosystem, the stability of coral reefs is directly related to the marine and surrounding terrestrial environments. However, coral reefs are now undergoing massive and accelerating devastation due to bleaching. The fact that the breakdown of symbiosis between coral and the dinoflagellate, zooxanthellae, has been well elucidated as the main cause of bleaching, implying the establishment of symbiosis with zooxanthellae plays a crucial role in maintaining coral survival. However, the relevant molecular and cellular mechanisms have not been well studied yet. In this study, based on the deep RNA-sequencing data derived from Mohamed, A. R. et al., an integrated transcriptome analysis was performed to deeply investigate global transcriptome changes of the coral Acropora digitifera in response to infection by dinoflagellate of the genus Symbiodinium.

RESULTS

The results revealed that compared to RefTranscriptome_v1.0 (A. digitifera transcriptome assembly v1.0), numerous novel transcripts and isoforms were identified, the Symbiodinium-infected coral larvae at 4 h generated the highest proportion of specific isoforms. Alternative splicing analysis showed that intron retention predominated in all alternative transcripts among six statuses. Additionally, 8117 lncRNAs were predicted via a stringent stepwise filtering pipeline. A complex lncRNAs-mRNAs network including 815 lncRNAs and 6395 mRNAs were established, in which 21 lncRNAs were differentially expressed at 4 h post infection. Functional clustering analysis for those differentially lncRNAs-coexpressed mRNAs demonstrated that several biological processes and pathways related to protein kinase activity, metabolic pathways, mitochondrion, ribosome, etc. were enriched.

CONCLUSIONS

Our study not only refined A. digitifera transcriptome via identification of novel transcripts and isoforms, but also predicted a high-confidence dataset of lncRNAs. Functional study based on the construction of lncRNAs-mRNAs co-expression network has disclosed a complex lncRNA-mediated regulation in response to Symbiodinium infection exhibited in A. digitifera. Once validated, these lncRNAs could be good potential targets for treatment and prevention of bleaching in coral.

Identification and analysis of long non-coding RNAs in response to H5N1 influenza viruses in duck (Anas platyrhynchos)

BACKGROUND

Long non-coding RNAs (lncRNAs) are important component of mammalian genomes, where their numbers are even larger than that of protein-coding genes. For example, human (Homo sapiens) (96,308 vs. 20,376) and mouse (Mus musculus) (87,774 vs. 22,630) have more lncRNA genes than protein-coding genes in the NONCODEv5 database. Recently, mammalian lncRNAs were reported to play critical roles in immune response to influenza A virus infections. Such observation inspired us to identify lncRNAs related to immune response to influenza A virus in duck, which is the most important natural host of influenza A viruses.

RESULTS

We explored features of 62,447 lncRNAs from human, mouse, chicken, zebrafish and elegans, and developed a pipeline to identify lncRNAs using the identified features with transcriptomic data. We then collected 151,970 assembled transcripts from RNA-Seq data of 21 individuals from three tissues and annotated 4094 duck lncRNAs. Comparing to duck protein-coding transcripts, we found that 4094 lncRNAs had smaller number of exons (2.4 vs. 10.2) and longer length of transcripts (1903.0 bp vs. 1686.9 bp) on average. Among them, 3586 (87.6%) lncRNAs located in intergenic regions and 619 lncRNAs showed differential expression in ducks infected by H5N1 virus when compared to control individuals. 58 lncRNAs were involved into two co-expressional modules related to anti-influenza A virus immune response. Moreover, we confirmed that eight lncRNAs showed remarkably differential expression both in vivo (duck individuals) and in vitro (duck embryo fibroblast cells, DEF cells) after infected with H5N1 viruses, implying they might play important roles in response to influenza A virus infection.

CONCLUSIONS

This study presented an example to annotate lncRNA in new species based on model species using transcriptome data. These data and analysis provide information for duck lncRNAs' function in immune response to influenza A virus.

Genetic and Molecular Basis of Feather Diversity in Birds

Feather diversity is striking in many aspects. Although the development of feather has been studied for decades, genetic and genomic studies of feather diversity have begun only recently. Many questions remain to be answered by multidisciplinary approaches. In this review, we discuss three levels of feather diversity: Feather morphotypes, intraspecific variations, and interspecific variations. We summarize recent studies of feather evolution in terms of genetics, genomics, and developmental biology and provide perspectives for future research. Specifically, this review includes the following topics: 1) Diversity of feather morphotype; 2) feather diversity among different breeds of domesticated birds, including variations in pigmentation pattern, in feather length or regional identity, in feather orientation, in feather distribution, and in feather structure; and 3) diversity of feathers among avian species, including plumage color and morph differences between species and the regulatory differences in downy feather development between altricial and precocial birds. Finally, we discussed future research directions.

Genome-wide differential expression of long noncoding RNAs and mRNAs in ovarian follicles of two different chicken breeds

Bashang long-tail chickens are an indigenous breed with dual purpose in China (meat and eggs) but have low egg laying performance. To improve the low egg laying performance, a genome-wide analysis of mRNAs and long noncoding RNAs (lncRNAs) from Bashang long-tail chickens and Hy-Line brown layers was performed. A total of 16,354 mRNAs and 8691 lncRNAs were obtained from ovarian follicles. Between the breeds, 160 mRNAs and 550 lncRNAs were found to be significantly differentially expressed. Integrated network analysis suggested some differentially expressed genes were involved in ovarian follicular development through oocyte meiosis, progesterone-mediated oocyte maturation, and cell cycle. The impact of lncRNAs on cis and trans target genes, indicating some lncRNAs may play important roles in ovarian follicular development. The current results provided a catalog of chicken ovarian follicular lncRNAs and genes for further study to understand their roles in regulation of egg laying performance.