The Tetraodon nigroviridis reference transcriptome: developmental transition, length retention and microsynteny of long non-coding RNAs in a compact vertebrate genome

High Diversity of Long Terminal Repeat Retrotransposons in Compact Vertebrate Genomes: Insights from Genomes of Tetraodontiformes

This study aimed to investigate the evolutionary profile (including diversity, activity, and abundance) of retrotransposons (RTNs) with long terminal repeats (LTRs) in ten species of Tetraodontiformes. These species, Arothron firmamentum, Lagocephalus sceleratus, Pao palembangensis, Takifugu bimaculatus, Takifugu flavidus, Takifugu ocellatus, Takifugu rubripes, Tetraodon nigroviridis, Mola mola, and Thamnaconus septentrionalis, are known for having the smallest genomes among vertebrates. Data mining revealed a high diversity and wide distribution of LTR retrotransposons (LTR-RTNs) in these compact vertebrate genomes, with varying abundances among species. A total of 819 full-length LTR-RTN sequences were identified across these genomes, categorized into nine families belonging to four different superfamilies: ERV (Orthoretrovirinae and Epsilon retrovirus), Copia, BEL-PAO, and Gypsy (Gmr, Mag, V-clade, CsRN1, and Barthez). The Gypsy superfamily exhibited the highest diversity. LTR family distribution varied among species, with Takifugu bimaculatus, Takifugu flavidus, Takifugu ocellatus, and Takifugu rubripes having the highest richness of LTR families and sequences. Additionally, evidence of recent invasions was observed in specific tetraodontiform genomes, suggesting potential transposition activity. This study provides insights into the evolution of LTR retrotransposons in Tetraodontiformes, enhancing our understanding of their impact on the structure and evolution of host genomes.

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.

Identification and analysis of long non-coding RNAs that are involved in response to GCRV infection in grass carp (Ctenopharyngodon idella)

Long noncoding RNAs (lncRNAs) play important roles in many biological processes including the immune response against virus infection. However, their roles in grass carp reovirus (GCRV) pathogenicity are largely unknown. In this study, the next-generation sequencing (NGS) technology was used to analyze the profiles of lncRNAs in GCRV-infected and mock-infected grass carp kidney (CIK) cells. Our results showed that 37 lncRNAs and 1039 mRNA transcripts exhibited differential expression in CIK cells after GCRV infection compared with the mock infection. Functional analysis through the gene ontology and Kyoto Encyclopedia of Genes and Genomes databases (KEGG) indicated that target genes of the differentially expressed lncRNAs were mainly enriched in the biological processes - biological regulation, cellular process, metabolic process and regulation of the biological process, such as MAPK signaling pathway and Notch signaling. Furthermore, we observed that the lncRNA3076 (ON693852) was markedly upregulated after the GCRV infection. In addition, silencing lncRNA3076 decreased the GCRV replication, which indicates that it might play an important role in the replication of GCRV.

The potential regulatory role of the lncRNA-miRNA-mRNA axis in teleost fish

Research over the past two decades has confirmed that noncoding RNAs (ncRNAs), which are abundant in cells from yeast to vertebrates, are no longer "junk" transcripts but functional regulators that can mediate various cellular and physiological processes. The dysregulation of ncRNAs is closely related to the imbalance of cellular homeostasis and the occurrence and development of various diseases. In mammals, ncRNAs, such as long noncoding RNAs (lncRNAs) and microRNAs (miRNAs), have been shown to serve as biomarkers and intervention targets in growth, development, immunity, and disease progression. The regulatory functions of lncRNAs on gene expression are usually mediated by crosstalk with miRNAs. The most predominant mode of lncRNA-miRNA crosstalk is the lncRNA-miRNA-mRNA axis, in which lncRNAs act as competing endogenous RNAs (ceRNAs). Compared to mammals, little attention has been given to the role and mechanism of the lncRNA-miRNA-mRNA axis in teleost species. In this review, we provide current knowledge about the teleost lncRNA-miRNA-mRNA axis, focusing on its physiological and pathological regulation in growth and development, reproduction, skeletal muscle, immunity to bacterial and viral infections, and other stress-related immune responses. Herein, we also explored the potential application of the lncRNA-miRNA-mRNA axis in the aquaculture industry. These findings contribute to an enhanced understanding of ncRNA and ncRNA-ncRNA crosstalk in fish biology to improve aquaculture productivity, fish health and quality.

Ancestrally Duplicated Conserved Noncoding Element Suggests Dual Regulatory Roles of HOTAIR in cis and trans

HOTAIR was proposed to regulate either HoxD cluster genes in trans or HoxC cluster genes in cis, a mechanism that remains unclear. We have identified a 32-nucleotide conserved noncoding element (CNE) as HOTAIR ancient sequence that likely originated at the root of vertebrate. The second round of whole-genome duplication resulted in one copy of the CNE within HOTAIR and another copy embedded in noncoding transcript of HOXD11. Paralogous CNEs underwent compensatory mutations, exhibit sequence complementarity with respect to transcripts directionality, and have high affinity in vitro. The HOTAIR CNE resembled a poised enhancer in stem cells and an active enhancer in HOTAIR-expressing cells. HOTAIR expression is positively correlated with HOXC11 in cis and negatively correlated with HOXD11 in trans. We propose a dual modality of HOTAIR regulation where transcription of HOTAIR and its embedded enhancer regulates HOXC11 in cis and sequence complementarity between paralogous CNEs suggests HOXD11 regulation in trans.

Identification and Evaluation of Long Noncoding RNAs in Response to Handling Stress in Red Cusk-Eel (Genypterus chilensis) via RNA-seq

The red cusk-eel (Genypterus chilensis) is a native species with strong potential to support Chilean aquaculture diversification. Under commercial conditions, fish are exposed to several stressors. To date, little is known about the mechanism involved in the stress response of red cusk-eel, and there is no information related to the regulation mediated by long noncoding RNAs (lncRNAs). The objective of this work was to identify for the first time the lncRNAs in the transcriptome of G. chilensis and to evaluate the differential expression levels of lncRNAs in the liver, head kidney, and skeletal muscle in response to handling stress. We used previously published transcriptome data to identify the lncRNAs by applying a series of filters based on annotation information in several databases to discard coding sequences. We identified a total of 14,614 putative lncRNAs in the transcriptome of red cusk-eel, providing a useful lncRNA reference resource to be used in future studies. We evaluated their differential expression in response to handling stress in the liver, head kidney, and skeletal muscle, identifying 112, 323, and 108 differentially expressed lncRNAs, respectively. The results suggest that handling stress in red cusk-eel generate an altered metabolic status in liver, altered immune response in head kidney, and skeletal muscle atrophy through an important coding and noncoding gene network. This is the first study that identifies lncRNAs in Genypterus genus and that evaluates the relation between handling stress and lncRNAs in teleost fish, thereby providing valuable information regarding noncoding responses to stress in Genypterus species.

Immunity-associated long non-coding RNA and expression in response to bacterial infection in large yellow croaker (Larimichthys crocea)

Long non-coding RNA refers to an RNA transcript of a non-coding protein with a sequence length greater than 200 bp. More and more reports indicated that lncRNA was involved in the regulation of gene expression as a signalling molecule, an inducing molecule, a leader molecule and a scaffold molecule. Previous studies have sequenced the draft genome and several transcriptome data sets for protein-coding genes of the large yellow croaker (Larimichthys crocea), but little is known about the expression and function of lncRNAs in this species. In order to obtain a catalogue of lncRNAs for this croaker, Vibrio parahaemolyticus infection challenge experiment was conducted and long non-coding RNA sequences were obtained. Using high-throughput sequencing of lncRNA, a total of 73,233 high-confidence transcripts were reconstructed in 32,726 loci, recovering most of the expressed reference transcripts, and 6473 novel expressed loci were identified. The tissue expression profile revealed that most lacunas were specifically enriched in distinct tissues. A set of 163 lncRNAs were identified as being specifically expressed in the spleen and may be involved in the immune response. It is the first time to identify specific lncRNAs in the L. crocea systematically in this croaker, aiming to benefit the future genomic study of this species.

Embryonic development of the bullseye puffer Sphoeroides annulatus (Tetraodontidae): A morphofunctional approach to ontogenetic steps

The embryonic development of the bullseye puffer, Sphoeroides annulatus, was characterized on the basis of the theory of saltatory ontogeny. This theory predicts a correlative relationship between the ontogeny-type in an altricial-precocial spectrum and the habitat that a species occupies within an unstable-stable environmental spectrum. Because S. annulatus inhabits a variety of unstable environments along a wide latitudinal range, the hypothesis that this species presents one of the most altricial embryonic developments among tetraodontids was tested. Based on major developmental events that marked the ontogenetic thresholds nine embryonic steps were identified. Developmental features such as small adhesives eggs, lack of vitelline circulation, small free embryos swimming up at hatching guided by positive phototaxis, and small first-feeding larvae actively swam in the water column, suggest that S. annulatus belongs to the reproductive guild of the nonguarders-lithopelagophils. Moreover, a comparative analysis of the developmental sequences, egg size, and first-feeding larvae size between tetraodontids confirms the hypothesis of this study and supports the evolutionary principle of the altricial-precocial spectrum postulated in the theory of saltatory ontogeny.

Inference of Developmental Gene Regulatory Networks Beyond Classical Model Systems: New Approaches in the Post-genomic Era

The advent of high-throughput sequencing (HTS) technologies has revolutionized the way we understand the transformation of genetic information into morphological traits. Elucidating the network of interactions between genes that govern cell differentiation through development is one of the core challenges in genome research. These networks are known as developmental gene regulatory networks (dGRNs) and consist largely of the functional linkage between developmental control genes, cis-regulatory modules, and differentiation genes, which generate spatially and temporally refined patterns of gene expression. Over the last 20 years, great advances have been made in determining these gene interactions mainly in classical model systems, including human, mouse, sea urchin, fruit fly, and worm. This has brought about a radical transformation in the fields of developmental biology and evolutionary biology, allowing the generation of high-resolution gene regulatory maps to analyze cell differentiation during animal development. Such maps have enabled the identification of gene regulatory circuits and have led to the development of network inference methods that can recapitulate the differentiation of specific cell-types or developmental stages. In contrast, dGRN research in non-classical model systems has been limited to the identification of developmental control genes via the candidate gene approach and the characterization of their spatiotemporal expression patterns, as well as to the discovery of cis-regulatory modules via patterns of sequence conservation and/or predicted transcription-factor binding sites. However, thanks to the continuous advances in HTS technologies, this scenario is rapidly changing. Here, we give a historical overview on the architecture and elucidation of the dGRNs. Subsequently, we summarize the approaches available to unravel these regulatory networks, highlighting the vast range of possibilities of integrating multiple technical advances and theoretical approaches to expand our understanding on the global gene regulation during animal development in non-classical model systems. Such new knowledge will not only lead to greater insights into the evolution of molecular mechanisms underlying cell identity and animal body plans, but also into the evolution of morphological key innovations in animals.

Non-coding RNAs Function as Immune Regulators in Teleost Fish

Non-coding RNAs (ncRNAs) are functional RNA molecules that are transcribed from DNA but not translated into proteins. ncRNAs function as key regulators of gene expression and chromatin modification. Recently, the functional role of ncRNAs in teleost fish has been extensively studied. Teleost fish are a highly diverse group among the vertebrate lineage. Fish are also important in terms of aquatic ecosystem, food production and human life, being the source of animal proteins worldwide and models of biomedical research. However, teleost fish always suffer from the invasion of infectious pathogens including viruses and bacteria, which has resulted in a tremendous economic loss to the fishing industry worldwide. Emerging evidence suggests that ncRNAs, especially miRNAs and lncRNAs, may serve as important regulators in cytokine and chemokine signaling, antigen presentation, complement and coagulation cascades, and T cell response in teleost fish. In this review, we summarize current knowledge and understanding of the roles of both miRNAs and lncRNAs in immune regulation in teleost fish. Molecular mechanism insights into the function of ncRNAs in fish immune response may contribute to the development of potential biomarkers and therapeutic targets for the prevention and treatment of fish diseases.

Tissue- and Cell Type-Specific Expression of the Long Noncoding RNA Klhl14-AS in Mouse

lncRNAs are acquiring increasing relevance as regulators in a wide spectrum of biological processes. The extreme heterogeneity in the mechanisms of action of these molecules, however, makes them very difficult to study, especially regarding their molecular function. A novel lncRNA has been recently identified as the most enriched transcript in mouse developing thyroid. Due to its genomic localization antisense to the protein-encoding Klhl14 gene, we named it Klhl14-AS. In this paper, we highlight that mouse Klhl14-AS produces at least five splicing variants, some of which have not been previously described. Klhl14-AS is expressed with a peculiar pattern, characterized by diverse relative abundance of its isoforms in different mouse tissues. We examine the whole expression level of Klhl14-AS in a panel of adult mouse tissues, showing that it is expressed in the thyroid, lung, kidney, testis, ovary, brain, and spleen, although at different levels. In situ hybridization analysis reveals that, in the context of each organ, Klhl14-AS shows a cell type-specific expression. Interestingly, databases report a similar expression profile for human Klhl14-AS. Our observations suggest that this lncRNA could play cell type-specific roles in several organs and pave the way for functional characterization of this gene in appropriate biological contexts.