Systematic identification of long noncoding RNAs expressed during zebrafish embryogenesis

Biomolecules Interacting with Long Noncoding RNAs

This review explores the complex interactions between long noncoding RNAs (lncRNAs) and other biomolecules, highlighting their pivotal roles in gene regulation and cellular function. LncRNAs, defined as RNA transcripts exceeding 200 nucleotides without encoding proteins, are involved in diverse biological processes, from embryogenesis to pathogenesis. They interact with DNA through mechanisms like triplex structure formation, influencing chromatin organization and gene expression. LncRNAs also modulate RNA-mediated processes, including mRNA stability, translational control, and splicing regulation. Their versatility stems from their forming of complex structures that enable interactions with various biomolecules. This review synthesizes current knowledge on lncRNA functions, discusses emerging roles in development and disease, and evaluates potential applications in diagnostics and therapeutics. By examining lncRNA interactions, it provides insights into the intricate regulatory networks governing cellular processes, underscoring the importance of lncRNAs in molecular biology. Unlike the majority of previous reviews that primarily focused on individual aspects of lncRNA biology, this comprehensive review uniquely integrates structural, functional, and mechanistic perspectives on lncRNA interactions across diverse biomolecules. Additionally, this review critically evaluates cutting-edge methodologies for studying lncRNA interactions, bridges fundamental molecular mechanisms with potential clinical applications, and highlights their potential.

Syntenic lncRNA locus exhibits DNA regulatory functions with sequence evolution

Syntenic long non-coding RNAs (lncRNAs) often show limited sequence conservation across species, prompting concern in the field. This study delves into functional signatures of syntenic lncRNAs between humans and zebrafish. Syntenic lncRNAs are highly expressed in zebrafish, with ∼90 % located near protein-coding genes, either in sense or antisense orientation. During early zebrafish development and in human embryonic stem cells (H1-hESC), syntenic lncRNA loci are enriched with cis-regulatory repressor signatures, influencing the expression of development-associated genes. In later zebrafish developmental stages and specific human cell lines, these syntenic lncRNA loci function as enhancers or transcription start sites (TSS) for protein-coding genes. Analysis of transposable elements (TEs) in syntenic lncRNA sequences revealed intriguing patterns: human lncRNAs are enriched in simple repeat elements, while their zebrafish counterparts show enrichment in LTR elements. This sequence evolution likely arises from post-rearrangement mutations that enhance DNA elements or cis-regulatory functions. It may also contribute to vertebrate innovation by creating novel transcription factor binding sites within the locus. This study highlights the conserved functionality of syntenic lncRNA loci through DNA elements, emphasizing their conserved roles across species despite sequence divergence.

The miR3367-lncRNA67-GhCYP724B module regulates male sterility by modulating brassinosteroid biosynthesis and interacting with Aorf27 in Gossypium hirsutum

Cytoplasmic male sterile (CMS) lines play a crucial role in utilization of heterosis in crop plants. However, the mechanism underlying the manipulation of male sterility in cotton by long non-coding RNA (lncRNA) and brassinosteroids (BRs) remains elusive. Here, using an integrative approach combining lncRNA transcriptomic profiles with virus-induced gene silencing experiments, we identify a flower bud-specific lncRNA in the maintainer line 2074B, lncRNA67, negatively modulating with male sterility in upland cotton (Gossypium hirsutum). lncRNA67 positively regulates cytochrome P274B (GhCYP724B), which acted as an eTM (endogenous target mimic) for miR3367. The suppression of GhCYP724B induced symptoms of BR deficiency and male semi-sterility in upland cotton as well as in tobacco, which resulted from a reduction in the endogenous BR contents. GhCYP724B regulates BRs synthesis by interacting with GhDIM and GhCYP90B, two BRs biosynthesis proteins. Additionally, GhCYP724B suppressed a unique chimeric open reading frame (Aorf27) in 2074A mitochondrial genome. Ectopic expression of Aorf27 in yeast inhibited cellular growth, and over expression of Aorf27 in tobacco showed male sterility. Overall, the results proved that the miR3367-lncRNA67-GhCYP724B module positively regulates male sterility by modulating BRs biosynthesis. The findings uncovered the function of lncRNA67-GhCYP724B in male sterility, providing a new mechanism for understanding male sterility in upland cotton.

Localization is the key to action: regulatory peculiarities of lncRNAs

To understand the transcriptomic profile of an individual cell in a multicellular organism, we must comprehend its surrounding environment and the cellular space where distinct molecular stimuli responses are located. Contradicting the initial perception that RNAs were nonfunctional and that only a few could act in chromatin remodeling, over the last few decades, research has revealed that they are multifaceted, versatile regulators of most cellular processes. Among the various RNAs, long non-coding RNAs (LncRNAs) regulate multiple biological processes and can even impact cell fate. In this sense, the subcellular localization of lncRNAs is the primary determinant of their functions. It affects their behavior by limiting their potential molecular partner and which process it can affect. The fine-tuned activity of lncRNAs is also tissue-specific and modulated by their cis and trans regulation. Hence, the spatial context of lncRNAs is crucial for understanding the regulatory networks by which they influence and are influenced. Therefore, predicting a lncRNA's correct location is not just a technical challenge but a critical step in understanding the biological meaning of its activity. Hence, examining these peculiarities is crucial to researching and discussing lncRNAs. In this review, we debate the spatial regulation of lncRNAs and their tissue-specific roles and regulatory mechanisms. We also briefly highlight how bioinformatic tools can aid research in the area.

Long non-coding RNAs in cancer: multifaceted roles and potential targets for immunotherapy

Cancer remains a major global health concern with high mortality rates mainly due to late diagnosis and poor prognosis. Long non-coding RNAs (lncRNAs) are emerging as key regulators of gene expression in human cancer, functioning through various mechanisms including as competing endogenous RNAs (ceRNAs) and indirectly regulating miRNA expression. LncRNAs have been found to have both oncogenic and tumor-suppressive roles in cancer, with the former promoting cancer cell proliferation, migration, invasion, and poor prognosis. Recent research has shown that lncRNAs are expressed in various immune cells and are involved in cancer cell immune escape and the modulation of the tumor microenvironment, thus highlighting their potential as targets for cancer immunotherapy. Targeting lncRNAs in cancer or immune cells could enhance the anti-tumor immune response and improve cancer immunotherapy outcomes. However, further research is required to fully understand the functional roles of lncRNAs in cancer and the immune system and their potential as targets for cancer immunotherapy. This review offers a comprehensive examination of the multifaceted roles of lncRNAs in human cancers, with a focus on their potential as targets for cancer immunotherapy. By exploring the intricate mechanisms underlying lncRNA-mediated regulation of cancer cell proliferation, invasion, and immune evasion, we provide insights into the diverse therapeutic applications of these molecules.

Transcriptome Analysis Reveals the lncRNA-mRNA Co-expression Network Regulating the Aestivation of Sea Cucumber

LncRNAs are long non-coding RNAs that are widely recognized as crucial regulators of gene expression and metabolic control, involved in numerous dormancy-related processes. Aestivation is a common hypometabolism strategy of sea cucumber (Apostichopus japonicus) in response to high-temperature conditions and is typically characterized by the degradation of the intestine and respiratory tree. Although the aestivation process has been extensively studied in sea cucumbers, the role of lncRNAs in the context of aestivation states remains a conspicuous knowledge gap. Here, we identified and characterized 14,711 lncRNAs in A. japonicus and analyzed their differential expression patterns during the aestivation process in the intestine and respiratory tree. The results revealed the physiological differences, especially the metabolic processes, between the intestine and respiratory tree during the aestivation. The co-expression network of lncRNA-mRNA suggested the dominant role of lncRNA in regulating the differential response of the intestine and respiratory trees. Differentially co-expressed factors were significantly enriched in the deep-aestivation stage-specific modules. Conserved co-expressed factors included several transcription factors known to be involved in rhythm regulation, such as Klf2 and Egr1. Furthermore, a specific trans-acting lncRNA (lncrna.1393.1) was identified as a potential regulator of Klf2 and Egr1. Overall, the systematic identification, characterization, and expression analysis of lncRNAs in A. japonicus enhanced our knowledge of long non-coding regulation of aestivation in sea cucumber and provided new clues for understanding the common "toolkit" of dormancy regulatory mechanisms.

SCAR-6 elncRNA locus epigenetically regulates PROZ and modulates coagulation and vascular function

In this study, we characterize a novel lncRNA-producing gene locus that we name Syntenic Cardiovascular Conserved Region-Associated lncRNA-6 (scar-6) and functionally validate its role in coagulation and cardiovascular function. A 12-bp deletion of the scar-6 locus in zebrafish (scar-6gib007Δ12/Δ12) results in cranial hemorrhage and vascular permeability. Overexpression, knockdown and rescue with the scar-6 lncRNA modulates hemostasis in zebrafish. Molecular investigation reveals that the scar-6 lncRNA acts as an enhancer lncRNA (elncRNA), and controls the expression of prozb, an inhibitor of factor Xa, through an enhancer element in the scar-6 locus. The scar-6 locus suppresses loop formation between prozb and scar-6 sequences, which might be facilitated by the methylation of CpG islands via the prdm14-PRC2 complex whose binding to the locus might be stabilized by the scar-6 elncRNA transcript. Binding of prdm14 to the scar-6 locus is impaired in scar-6gib007Δ12/Δ12 zebrafish. Finally, activation of the PAR2 receptor in scar-6gib007Δ12/Δ12 zebrafish triggers NF-κB-mediated endothelial cell activation, leading to vascular dysfunction and hemorrhage. We present evidence that the scar-6 locus plays a role in regulating the expression of the coagulation cascade gene prozb and maintains vascular homeostasis.

Evolutionary comparison of lncRNAs in four cotton species and functional identification of LncR4682-PAS2-KCS19 module in fiber elongation

Long non-coding RNAs (lncRNAs) play an important role in various biological processes in plants. However, there have been few reports on the evolutionary signatures of lncRNAs in closely related cotton species. The lncRNA transcription patterns in two tetraploid cotton species and their putative diploid ancestors were compared in this paper. By performing deep RNA sequencing, we identified 280 429 lncRNAs from 21 tissues in four cotton species. lncRNA transcription evolves more rapidly than mRNAs, and exhibits more severe turnover phenomenon in diploid species compared to that in tetraploid species. Evolutionarily conserved lncRNAs exhibit higher expression levels, and lower tissue specificity compared with species-specific lncRNAs. Remarkably, tissue expression of homologous lncRNAs in Gossypium hirsutum and G. barbadense exhibited similar patterns, suggesting that these lncRNAs may be functionally conserved and selectively maintained during domestication. An orthologous lncRNA, lncR4682, was identified and validated in fibers of G. hirsutum and G. barbadense with the highest conservatism and expression abundance. Through virus-induced gene silencing in upland cotton, we found that lncR4682 and its target genes GHPAS2 and GHKCS19 positively regulated fiber elongation. In summary, the present study provides a systematic analysis of lncRNAs in four closely related cotton species, extending the understanding of transcriptional conservation of lncRNAs across cotton species. In addition, LncR4682-PAS2-KCS19 contributes to cotton fiber elongation by participating in the biosynthesis of very long-chain fatty acids.

Systematic profiling and analysis of growth and development responsive DE-lncRNAs in cluster bean (Cyamopsis tetragonoloba)

Long non-coding RNAs (lncRNAs) play crucial role in regulating genes involved in various processes including growth & development, flowering, and stress response in plants. The study aims to identify and characterize tissue-specific, growth & development and floral responsive differentially expressed lncRNAs (DE-lncRNAs) in cluster bean from a high-throughput RNA sequencing data. We have identified 3309 DE-lncRNAs, with an average length of 818 bp. Merely, around 4 % of DE-lncRNAs across the tissues were found to be conserved as rate of evolution of lncRNAs is high. Among the identified DE-lncRNAs, 204 were common in leaf vs. shoot, leaf vs. flower and flower vs. shoot. A total of 60 DE-lncRNAs targeted 10 protein-coding genes involved in flower development and initiation processes. We investigated 179 tissue-specific DE-lncRNAs based on tissue specificity index. Three DE-lncRNAs: Cb_lnc_0820, Cb_lnc_0430, Cb_lnc_0260 and their target genes show their involvement in floral development and stress mechanisms, which were validated by Quantitative real-time PCR (qRT-PCR). The identified DE-lncRNAs were expressed higher in flower bud than in leaf and similar expression pattern was observed in both RNA-seq data and qRT-PCR analyses. Notably, 362 DE-lncRNAs were predicted as eTM-lncRNAs with the participation of 84 miRNAs. Whereas 46 DE-lncRNAs were predicted to possess the internal ribosomal entry sites (IRES) and can encode for small peptides. The regulatory networks established between DE-lncRNAs, mRNAs and miRNAs have provided an insight into their association with plant growth & development, flowering, and stress mechanisms. Comprehensively, the characterization of DE-lncRNAs in various tissues of cluster bean shed a light on interactions among lncRNAs, miRNAs and mRNAs and help understand their involvement in growth & development and floral initiation processes. The information retrieved from the analyses was shared in the public domain in the form of a database: Cb-DElncRNAdb, and made available at http://backlin.cabgrid.res.in/Cb-DElncRNA/index.php, which may be useful for the scientific community engaged cluster bean research.

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.

The effects of exercise on epigenetic modifications: focus on DNA methylation, histone modifications and non-coding RNAs

Physical activity on a regular basis has been shown to bolster the overall wellness of an individual; research is now revealing that these changes are accompanied by epigenetic modifications. Regular exercise has been proven to make intervention plans more successful and prolong adherence to them. When it comes to epigenetic changes, there are four primary components. This includes changes to the DNA, histones, expression of particular non-coding RNAs and DNA methylation. External triggers, such as physical activity, can lead to modifications in the epigenetic components, resulting in changes in the transcription process. This report pays attention to the current knowledge that pertains to the epigenetic alterations that occur after exercise, the genes affected and the resulting characteristics.

Evolution of diapause in the African turquoise killifish by remodeling the ancient gene regulatory landscape

Suspended animation states allow organisms to survive extreme environments. The African turquoise killifish has evolved diapause as a form of suspended development to survive a complete drought. However, the mechanisms underlying the evolution of extreme survival states are unknown. To understand diapause evolution, we performed integrative multi-omics (gene expression, chromatin accessibility, and lipidomics) in the embryos of multiple killifish species. We find that diapause evolved by a recent remodeling of regulatory elements at very ancient gene duplicates (paralogs) present in all vertebrates. CRISPR-Cas9-based perturbations identify the transcription factors REST/NRSF and FOXOs as critical for the diapause gene expression program, including genes involved in lipid metabolism. Indeed, diapause shows a distinct lipid profile, with an increase in triglycerides with very-long-chain fatty acids. Our work suggests a mechanism for the evolution of complex adaptations and offers strategies to promote long-term survival by activating suspended animation programs in other species.

A long noncoding RNA functions in pumpkin fruit development through S-adenosyl-L-methionine synthetase

Long noncoding RNAs (lncRNAs) play important roles in various biological processes. However, the regulatory roles of lncRNAs underlying fruit development have not been extensively studied. The pumpkin (Cucurbita spp.) is a preferred model for understanding the molecular mechanisms regulating fruit development because of its variable shape and size and large inferior ovary. Here, we performed strand-specific transcriptome sequencing on pumpkin (Cucurbita maxima "Rimu") fruits at 6 developmental stages and identified 5,425 reliably expressed lncRNAs. Among the 332 lncRNAs that were differentially expressed during fruit development, the lncRNA MSTRG.44863.1 was identified as a negative regulator of pumpkin fruit development. MSTRG.44863.1 showed a relatively high expression level and an obvious period-specific expression pattern. Transient overexpression and silencing of MSTRG.44863.1 significantly increased and decreased the content of 1-aminocyclopropane carboxylic acid (a precursor of ethylene) and ethylene production, respectively. RNA pull-down and microscale thermophoresis assays further revealed that MSTRG.44863.1 can interact with S-adenosyl-L-methionine synthetase (SAMS), an enzyme in the ethylene synthesis pathway. Considering that ethylene negatively regulates fruit development, these results indicate that MSTRG.44863.1 plays an important role in the regulation of pumpkin fruit development, possibly through interacting with SAMS and affecting ethylene synthesis. Overall, our findings provide a rich resource for further study of fruit-related lncRNAs while offering insights into the regulation of fruit development in plants.

Cell-type-specific mRNA transcription and degradation kinetics in zebrafish embryogenesis from metabolically labeled single-cell RNA-seq

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the relative contributions of mRNA transcription and degradation to shaping those profiles remains challenging, especially within embryos with diverse cellular identities. Here, we combine single-cell RNA-Seq and metabolic labeling to capture temporal cellular transcriptomes of zebrafish embryos where newly-transcribed (zygotic) and pre-existing (maternal) mRNA can be distinguished. We introduce kinetic models to quantify mRNA transcription and degradation rates within individual cell types during their specification. These models reveal highly varied regulatory rates across thousands of genes, coordinated transcription and destruction rates for many transcripts, and link differences in degradation to specific sequence elements. They also identify cell-type-specific differences in degradation, namely selective retention of maternal transcripts within primordial germ cells and enveloping layer cells, two of the earliest specified cell types. Our study provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.

V. S, Kaushik K, Sabharwal A, Lalwani MK, K. S, Singh N, Scaria V, Sivasubbu S. Forward genetic screen using a gene-breaking trap approach identifies a novel role of grin2bb-associated RNA transcript (grin2bbART) in zebrafish heart function

LncRNA-based control affects cardiac pathophysiologies like myocardial infarction, coronary artery disease, hypertrophy, and myotonic muscular dystrophy. This study used a gene-break transposon (GBT) to screen zebrafish (Danio rerio) for insertional mutagenesis. We identified three insertional mutants where the GBT captured a cardiac gene. One of the adult viable GBT mutants had bradycardia (heart arrhythmia) and enlarged cardiac chambers or hypertrophy; we named it "bigheart." Bigheart mutant insertion maps to grin2bb or N-methyl D-aspartate receptor (NMDAR2B) gene intron 2 in reverse orientation. Rapid amplification of adjacent cDNA ends analysis suggested a new insertion site transcript in the intron 2 of grin2bb. Analysis of the RNA sequencing of wild-type zebrafish heart chambers revealed a possible new transcript at the insertion site. As this putative lncRNA transcript satisfies the canonical signatures, we called this transcript grin2bb associated RNA transcript (grin2bbART). Using in situ hybridization, we confirmed localized grin2bbART expression in the heart, central nervous system, and muscles in the developing embryos and wild-type adult zebrafish atrium and bulbus arteriosus. The bigheart mutant had reduced Grin2bbART expression. We showed that bigheart gene trap insertion excision reversed cardiac-specific arrhythmia and atrial hypertrophy and restored grin2bbART expression. Morpholino-mediated antisense downregulation of grin2bbART in wild-type zebrafish embryos mimicked bigheart mutants; this suggests grin2bbART is linked to bigheart. Cardiovascular tissues use Grin2bb as a calcium-permeable ion channel. Calcium imaging experiments performed on bigheart mutants indicated calcium mishandling in the heart. The bigheart cardiac transcriptome showed differential expression of calcium homeostasis, cardiac remodeling, and contraction genes. Western blot analysis highlighted Camk2d1 and Hdac1 overexpression. We propose that altered calcium activity due to disruption of grin2bbART, a putative lncRNA in bigheart, altered the Camk2d-Hdac pathway, causing heart arrhythmia and hypertrophy in zebrafish.

From Beats to Metabolism: the Heart at the Core of Interorgan Metabolic Cross Talk

The heart, once considered a mere blood pump, is now recognized as a multifunctional metabolic and endocrine organ. Its function is tightly regulated by various metabolic processes, at the same time it serves as an endocrine organ, secreting bioactive molecules that impact systemic metabolism. In recent years, research has shed light on the intricate interplay between the heart and other metabolic organs, such as adipose tissue, liver, and skeletal muscle. The metabolic flexibility of the heart and its ability to switch between different energy substrates play a crucial role in maintaining cardiac function and overall metabolic homeostasis. Gaining a comprehensive understanding of how metabolic disorders disrupt cardiac metabolism is crucial, as it plays a pivotal role in the development and progression of cardiac diseases. The emerging understanding of the heart as a metabolic and endocrine organ highlights its essential contribution to whole body metabolic regulation and offers new insights into the pathogenesis of metabolic diseases, such as obesity, diabetes, and cardiovascular disorders. In this review, we provide an in-depth exploration of the heart's metabolic and endocrine functions, emphasizing its role in systemic metabolism and the interplay between the heart and other metabolic organs. Furthermore, emerging evidence suggests a correlation between heart disease and other conditions such as aging and cancer, indicating that the metabolic dysfunction observed in these conditions may share common underlying mechanisms. By unraveling the complex mechanisms underlying cardiac metabolism, we aim to contribute to the development of novel therapeutic strategies for metabolic diseases and improve overall cardiovascular health.

eIF4E1b is a non-canonical eIF4E protecting maternal dormant mRNAs

Maternal mRNAs are essential for protein synthesis during oogenesis and early embryogenesis. To adapt translation to specific needs during development, maternal mRNAs are translationally repressed by shortening the polyA tails. While mRNA deadenylation is associated with decapping and degradation in somatic cells, maternal mRNAs with short polyA tails are stable. Here we report that the germline-specific eIF4E paralog, eIF4E1b, is essential for zebrafish oogenesis. eIF4E1b localizes to P-bodies in zebrafish embryos and binds to mRNAs with reported short or no polyA tails, including histone mRNAs. Loss of eIF4E1b results in reduced histone mRNA levels in early gonads, consistent with a role in mRNA storage. Using mouse and human eIF4E1Bs (in vitro) and zebrafish eIF4E1b (in vivo), we show that unlike canonical eIF4Es, eIF4E1b does not interact with eIF4G to initiate translation. Instead, eIF4E1b interacts with the translational repressor eIF4ENIF1, which is required for eIF4E1b localization to P-bodies. Our study is consistent with an important role of eIF4E1b in regulating mRNA dormancy and provides new insights into fundamental post-transcriptional regulatory principles governing early vertebrate development.

Identification of Dezhou donkey muscle development-related genes and long non-coding RNA based on differential expression analysis

Long non-coding RNAs (lncRNAs) play a critical role in the development of muscles. However, the role of lncRNAs in regulating skeletal muscle development has not been studied systematically in the donkey. In this study, we performed the RNA sequencing for different stages of muscles in donkeys, and investigate their expression profile, which showed that 3215 mRNAs (p-adjust <0.05) and 471 lncRNAs (p-value <0.05) were significantly differently expressed (DE) verified by RT-qPCR. GO and KEGG enrichment analysis indicated that DE genes and target genes of DE lncRNAs were associated with muscle development in the donkey. We also found these four target genes (DCN, ITM2A, MUSTN1, ARRDC2) involved in skeletal muscle growth and development. Combined with transcriptome data, network, and RT-qPCR results showed that four co-expression networks of DCN and lnc-008278, ITM2A and lnc_017247, MUSTN1 and lnc_030153, and ARRDC2 and lnc_033914, which may play an important role in the formation and development of muscle in the donkey.

The role of microRNAs in the pathophysiology of human central nervous system: A focus on neurodegenerative diseases

microRNAs (miRNAs) are suggested to play substantial roles in regulating the development and various physiologic functions of the central nervous system (CNS). These include neurogenesis, cell fate and differentiation, morphogenesis, formation of dendrites, and targeting non-neural mRNAs. Notably, deregulation of an increasing number of miRNAs is associated with several neurodegenerative diseases including Alzheimer's disease, Parkinson's disease, multiple sclerosis, amyotrophic lateral sclerosis and CNS tumors. They are particularly known to affect the amyloid β (Aβ) cleavage and accumulation, tau protein homeostasis, and expression of alpha-synuclein (α-syn), Parkin, PINK1, and brain-derived neurotrophic factor (BDNF) that play pivotal roles in the pathogenesis of neurodegenerative diseases. These include miR-16, miR-17-5p, miR-20a, miR-106a, miR-106b, miR-15a, miR-15b, miR-103, miR-107, miR-298, miR-328, miR-195, miR-485, and miR-29. In CNS tumors, several miRNAs, including miR-31, miR-16, and miR-21 have been identified to modulate tumorigenesis through impacting tumor invasion and apoptosis. In this review article, we have a look at the recent advances on our knowledge about the role of miRNAs in human brain development and functions, neurodegenerative diseases, and their clinical potentials.

LncRNA MAFG-AS1 is involved in human cancer progression

Long noncoding RNAs (lncRNAs) refer to a type of non-protein-coding transcript of more than 200 nucleotides. LncRNAs play fundamental roles in disease development and progression, and lncRNAs are dysregulated in many pathophysiological processes. Thus, lncRNAs may have potential value in clinical applications. The lncRNA, MAF BZIP Transcription Factor G (MAFG)-AS1, is dysregulated in several cancer, including breast cancer, lung cancer, liver cancer, bladder cancer, colorectal cancer, gastric cancer, esophagus cancer, prostate cancer, pancreatic cancer, ovarian cancer, and glioma. Altered MAFG-AS1 levels are also associated with diverse clinical characteristics and patient outcomes. Mechanistically, MAFG-AS1 mediates a variety of cellular processes via the regulation of target gene expression. Therefore, the diagnostic, prognostic, and therapeutic aspects of MAFG-AS1 have been widely explored. In this review, we discuss the expression, major roles, and molecular mechanisms of MAFG-AS1, the relationship between MAFG-AS1 and clinical features of diseases, and the clinical applications of MAFG-AS1.

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.

Deciphering shared attributes of plant long non-coding RNAs through a comparative computational approach

Over the past decade, long non-coding RNA (lncRNA), which lacks protein-coding potential, has emerged as an essential regulator of the genome. The present study examined 13,599 lncRNAs in Arabidopsis thaliana, 11,565 in Oryza sativa, and 32,397 in Zea mays for their characteristic features and explored the associated genomic and epigenomic features. We found lncRNAs were distributed throughout the chromosomes and the Helitron family of transposable elements (TEs) enriched, while the terminal inverted repeat depleted in lncRNA transcribing regions. Our analyses determined that lncRNA transcribing regions show rare or weak signals for most epigenetic marks except for H3K9me2 and cytosine methylation in all three plant species. LncRNAs showed preferential localization in the nucleus and cytoplasm; however, the distribution ratio in the cytoplasm and nucleus varies among the studied plant species. We identified several conserved endogenous target mimic sites in the lncRNAs among the studied plants. We found 233, 301, and 273 unique miRNAs, potentially targeting the lncRNAs of A. thaliana, O. sativa, and Z. mays, respectively. Our study has revealed that miRNAs, which interact with lncRNAs, target genes that are involved in a diverse array of biological and molecular processes. The miRNA-targeted lncRNAs displayed a strong affinity for several transcription factors, including ERF and BBR-BPC, mutually present in all three plants, advocating their conserved functions. Overall, the present study showed that plant lncRNAs exhibit conserved genomic and epigenomic characteristics and potentially govern the growth and development of plants.

Genome-wide analysis of long noncoding RNA profiles in pseudorabies-virus-infected PK15 cells

Recently, an increasing number of studies have shown that long noncoding RNAs (lncRNAs) are involved in host metabolism after infection with pseudorabies virus (PRV). In our study, via RNA sequencing analysis, a total of 418 mRNAs, 137 annotated lncRNAs, and 312 new lncRNAs were found to be differentially expressed. These lncRNAs were closely associated with metabolic regulation and immunity-related signalling pathways, including the T-cell receptor signalling pathway, chemokine signalling pathway, mitogen-activated protein kinase (MAPK) signalling pathway, TNF signalling pathway, Ras signalling pathway, calcium signalling pathway, and phosphatidylinositol signalling system. Real-time PCR indicated that several mRNAs and lncRNAs involved in the regulation of the immune effector process, T-cell receptor signalling pathway, TNF signalling pathway, MAPK signalling pathway, and chemokine signalling pathways were significantly expressed. These mRNAs and lncRNAs might play a role in PRV infection.

Epigenetic mechanisms of lncRNA in response to thermal stress during embryogenesis of allotetraploid Cyprinus carpio

Embryogenesis and epigenetic mechanisms of lncRNA may play an important role in the formation of temperature tolerance in allotetraploid Cyprinus carpio. To investigate the response of lncRNA to thermal stress during embryogenesis of C. carpio, transcriptome sequencing was performed on 81 embryo or larva samples from different early development stages and temperatures. We identified 45,097 lncRNAs and analyzed transcriptome variation during embryogenesis. Stage-specific and temperature-specific DE lncRNAs and DEGs were screened. GO and KEGG analysis identified numerous pathways involved in thermal stress. Temperature-specific regulation of cis-/trans-/antisense lncRNAs was analyzed. Interaction network analysis identified 6 hub lncRNAs and many hub genes, such as cdk1 and hsf1. Decreased expression of many essential genes regulated by lncRNAs may lead to the death of embryos at 33 °C. Our findings provide new insights into the regulation of lncRNA in thermal stress response during embryogenesis and contribute to the understanding of environmental adaptation of allotetraploid species.

Discovery of putative long non-coding RNAs expressed in the eyes of Astyanax mexicanus (Actinopterygii: Characidae)

Astyanax mexicanus is a well-known model species, that has two morphotypes, cavefish, from subterranean rivers and surface fish, from surface rivers. They are morphologically distinct due to many troglomorphic traits in the cavefish, such as the absence of eyes. Most studies on A. mexicanus are focused on eye development and protein-coding genes involved in the process. However, lncRNAs did not get the same attention and very little is known about them. This study aimed to fill this knowledge gap, identifying, describing, classifying, and annotating lncRNAs expressed in the embryo's eye tissue of cavefish and surface fish. To do so, we constructed a concise workflow to assemble and evaluate transcriptomes, annotate protein-coding genes, ncRNAs families, predict the coding potential, identify putative lncRNAs, map them and predict interactions. This approach resulted in the identification of 33,069 and 19,493 putative lncRNAs respectively mapped in cavefish and surface fish. Thousands of these lncRNAs were annotated and identified as conserved in human and several species of fish. Hundreds of them were validated in silico, through ESTs. We identified lncRNAs associated with genes related to eye development. This is the case of a few lncRNAs associated with sox2, which we suggest being isomorphs of the SOX2-OT, a lncRNA that can regulate the expression of sox2. This work is one of the first studies to focus on the description of lncRNAs in A. mexicanus, highlighting several lncRNA targets and opening an important precedent for future studies focusing on lncRNAs expressed in A. mexicanus.

Natural antisense transcription of presenilin in sea urchin reveals a possible role for natural antisense transcription in the general control of gene expression during development

One presenilin gene (PSEN) is expressed in the sea urchin embryo, in the vegetal pole of the gastrula and then mainly in cilia cells located around the digestive system of the pluteus, as we recently have reported. PSEN expression must be accurately regulated for correct execution of these two steps of development. While investigating PSEN expression changes in embryos after expansion of endoderm with LiCl or of ectoderm with Zn2+ by whole-mount in situ hybridization (WISH) and quantitative PCR (qPCR), we detected natural antisense transcription of PSEN. We then found that Endo16 and Wnt5, markers of endo-mesoderm, and of Hnf6 and Gsc, markers of ectoderm, are also sense and antisense transcribed. We discuss that general gene expression could depend on both sense and antisense transcription. This mechanism, together with the PSEN gene, should be included in gene regulatory networks (GRNs) that theorize diverse processes in this species. We suggest that it would also be relevant to investigate natural antisense transcription of PSEN in the field of Alzheimer's disease (AD) where the role of human PSEN1 and PSEN2 is well known.

LncRNA Functional Screening in Organismal Development

Controversy continues over the functional prevalence of long non-coding RNAs (lncRNAs) despite their being widely investigated in all kinds of cells and organisms. In animals, lncRNAs have aroused general interest from exponentially increasing transcriptomic repertoires reporting their highly tissue-specific and developmentally dynamic expression, and more importantly, from growing experimental evidence supporting their functionality in facilitating organogenesis and individual fitness. In mammalian testes, while a great multitude of lncRNA species are identified, only a minority of them have been shown to be useful, and even fewer have been demonstrated as true requirements for male fertility using knockout models to date. This noticeable gap is attributed to the virtual existence of a large number of junk lncRNAs, the lack of an ideal germline culture system, difficulty in loss-of-function interrogation, and limited screening strategies. Facing these challenges, in this review, we discuss lncRNA functionality in organismal development and especially in mouse testis, with a focus on lncRNAs with functional screening.

LncRNA109897-JrCCR4-JrTLP1b forms a positive feedback loop to regulate walnut resistance against anthracnose caused by Colletotrichum gloeosporioides

Walnut anthracnose induced by Colletotrichum gloeosporioides is a disastrous disease that severely restricts the development of the walnut industry in China. Long non-coding RNAs (lncRNAs) are involved in adaptive responses to disease, but their roles in the regulation of walnut anthracnose resistance response are not well defined. In this study, transcriptome analysis demonstrated that a C. gloeosporioides-induced lncRNA, lncRNA109897, located upstream from the target gene JrCCR4, upregulated the expression of JrCCR4. JrCCR4 interacted with JrTLP1b and promoted its transcriptional activity. In turn, JrTLP1b induced the transcription of lncRNA109897 to promote its expression. Meanwhile, transient expression in walnut leaves and stable transformation of Arabidopsis thaliana further proved that lncRNA, JrCCR4, and JrTLP1b improve the resistance of C. gloeosporioides. Collectively, these findings provide insights into the mechanism by which the lncRNA109897-JrCCR4-JrTLP1b transcriptional cascade regulates the resistance of walnut to anthracnose.

Transcriptomics analysis reveals key lncRNAs and genes related to the infection of feline kidney cell line by panleukopenia virus

Feline panleukopenia virus (FPV) can cause a viral disease and is responsible for severe leukopenia, gastroenteritis, and nervous signs with significant economic losses. Biochemically long non-coding RNAs (lncRNAs) can regulate the expression of mRNA in different ways, thereby causing the functional changes in host cells in response to viral infection. However, no attention has been paid until now to investigate the link between FPV pathogenesis and lncRNA. Here, through RNA sequencing, we performed a comprehensive analysis of lncRNA and mRNA in F81 cells after FPV-BJ04 strain infection. Consistent with previous studies, our data showed that lncRNAs have distinct features from mRNA. A total of 291 lncRNAs and 873 mRNAs were differentially expressed in F81 cells after FPV-BJ04 infection. GO and KEGG enrichment analysis showed that the differentially upregulated lncRNAs target genes were mainly involved in the positive regulation of transcription by RNA polymerase II and MAPK signaling pathway. The differentially downregulated lncRNAs target genes were mainly involved in the mRNA splicing and endocytosis. In addition, the differentially expressed immune pathway related genes that are targeted by lncRNA were also screened out to construct a lncRNA-miRNA-mRNA axes as a potential novel biomarkers in regulating the immune response of feline against FPV infection. Our results contribute to understand the basic role of lncRNA in F81 cells during FPV infection and lay the foundation for following research.

Single-cell temporal dynamics reveals the relative contributions of transcription and degradation to cell-type specific gene expression in zebrafish embryos

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the underlying regulation of mRNA transcription and degradation remains a challenge, especially within whole embryos with diverse cellular identities. Here, we collect temporal cellular transcriptomes of zebrafish embryos, and decompose them into their newly-transcribed (zygotic) and pre-existing (maternal) mRNA components by combining single-cell RNA-Seq and metabolic labeling. We introduce kinetic models capable of quantifying regulatory rates of mRNA transcription and degradation within individual cell types during their specification. These reveal different regulatory rates between thousands of genes, and sometimes between cell types, that shape spatio-temporal expression patterns. Transcription drives most cell-type restricted gene expression. However, selective retention of maternal transcripts helps to define the gene expression profiles of germ cells and enveloping layer cells, two of the earliest specified cell-types. Coordination between transcription and degradation restricts expression of maternal-zygotic genes to specific cell types or times, and allows the emergence of spatio-temporal patterns when overall mRNA levels are held relatively constant. Sequence-based analysis links differences in degradation to specific sequence motifs. Our study reveals mRNA transcription and degradation events that control embryonic gene expression, and provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.

Diagnostic and therapeutic potential of LINC01929 as an oncogenic LncRNA in human cancers

Recent findings have shown the significant role of long non-coding RNAs in the pathogenesis of various cancers. In this regard, the variation in the expression of LINC01929 was explored in various cancers to explore its impact on the development of diverse malignancies and cancers. The data of the cancer genome atlas (TCGA) were utilized to evaluate the changes in the expression of LINC01929 in various cancers, as well as its relationship with the patients' survival rate. The co-expression of the genes and data merging of TCGA were utilized to identify the LINC01929-associated pathways. The samples of colorectal, gastric, and breast cancers were also examined by the RT-qPCR to confirm the results and evaluate the expression of LINC01929 in the mentioned cancers. In silico investigations indicated a remarkable enhancement in the expression of LINC01929 within the tumor tissues compared to normal samples in 10 types of cancer. Based on the survival results, the increase in the LINC01929 expression is linked to poor prognosis of bladder, breast, colorectal, kidney, and liver cancers. The gene co-expression network showed the strong co-expression of LINC01929 with genes involved in the metastatic pathways including COL5A1. RT-qPCR findings showed a remarkable increment in the expression level of LINC01929 in the colorectal, gastric, and breast tumor tissues versus the adjacent normal tissues. A significant and strong relationship was also found between the expression of LINC01929 and COL5A1. This study indicated a significant enhancement in the expression level of LINC01929 in various cancer types, accompanied by the mortality rate. Moreover, LINC01929 exhibited a strong co-expression with the metastatic genes such as COL5A1. As an oncogene and regulator of the metastatic pathways, LINC01929 can be a proper candidate for diagnostic and therapeutic purposes.

SLAMseq resolves the kinetics of maternal and zygotic gene expression during early zebrafish embryogenesis

The maternal-to-zygotic transition (MZT) is a key developmental process in metazoan embryos that involves the activation of zygotic transcription (ZGA) and degradation of maternal transcripts. We employed metabolic mRNA sequencing (SLAMseq) to deconvolute the compound embryonic transcriptome in zebrafish. While mitochondrial zygotic transcripts prevail prior to MZT, we uncover the spurious transcription of hundreds of short and intron-poor genes as early as the 2-cell stage. Upon ZGA, most zygotic transcripts originate from thousands of maternal-zygotic (MZ) genes that are transcribed at rates comparable to those of hundreds of purely zygotic genes and replenish maternal mRNAs at distinct timescales. Rapid replacement of MZ transcripts involves transcript decay features unrelated to major maternal degradation pathways and promotes de novo synthesis of the core gene expression machinery by increasing poly(A)-tail length and translation efficiency. SLAMseq hence provides insights into the timescales, molecular features, and regulation of MZT during zebrafish embryogenesis.

Two Novel lncRNAs Regulate Primordial Germ Cell Development in Zebrafish

Long noncoding RNAs (lncRNAs) are regulatory transcripts in various biological processes. However, the role of lncRNAs in germline development remains poorly understood, especially for fish primordial germ cell (PGC) development. In this study, the lncRNA profile of zebrafish PGC was revealed by single cell RNA-sequencing and bioinformatic prediction. We established the regulation network of lncRNA-mRNA associated with PGC development, from which we identified three novel lncRNAs-lnc172, lnc196, and lnc304-highly expressing in PGCs and gonads. Fluorescent in situ hybridization indicated germline-specific localization of lnc196 and lnc304 in the cytoplasm and nucleus of spermatogonia, spermatocyte, and occyte, and they were co-localized with vasa in the cytoplasm of the spermatogonia. By contrast, lnc172 was localized in the cytoplasm of male germline, myoid cells and ovarian somatic cells. Loss- and gain-of-function experiments demonstrated that knockdown and PGC-specific overexpression of lnc304 as well as universal overexpression of lnc172 significantly disrupted PGC development. In summary, the present study revealed the lncRNA profile of zebrafish PGC and identified two novel lncRNAs associated with PGC development, providing new insights for understanding the regulatory mechanism of PGC development.

Research progress of lncRNA and miRNA in hepatic ischemia-reperfusion injury

BACKGROUND

Hepatic ischemia-reperfusion injury (HIRI) is a common complication of liver surgeries, such as hepatectomy and liver transplantation. In recent years, several non-coding RNAs (ncRNAs) including long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) have been identified as factors involved in the pathological progression of HIRI. In this review, we summarized the latest research on lncRNAs, miRNAs and the lncRNA-miRNA regulatory networks in HIRI.

DATA SOURCES

The PubMed and Web of Science databases were searched for articles published up to December 2021 using the following keywords: "hepatic ischemia-reperfusion injury", "lncRNA", "long non-coding RNA", "miRNA" and "microRNA". The bibliography of the selected articles was manually screened to identify additional studies.

RESULTS

The mechanism of HIRI is complex, and involves multiple lncRNAs and miRNAs. The roles of lncRNAs such as AK139328, CCAT1, MALAT1, TUG1 and NEAT1 have been established in HIRI. In addition, numerous miRNAs are associated with apoptosis, autophagy, oxidative stress and cellular inflammation that accompany HIRI pathogenesis. Based on the literature, we conclude that four lncRNA-miRNA regulatory networks mediate the pathological progression of HIRI. Furthermore, the expression levels of some lncRNAs and miRNAs undergo significant changes during the progression of HIRI, and thus are potential prognostic markers and therapeutic targets.

CONCLUSIONS

Complex lncRNA-miRNA-mRNA networks regulate HIRI progression through mutual activation and antagonism. It is necessary to screen for more HIRI-associated lncRNAs and miRNAs in order to identify novel therapeutic targets.

Transcriptomic Analysis of Long Non-Coding RNA during Candida albicans Infection

Candida albicans is one of the most commonly found species in fungal infections. Due to its clinical importance, molecular aspects of the host immune defense against the fungus are of interest to biomedical sciences. Long non-coding RNAs (lncRNAs) have been investigated in different pathologies and gained widespread attention regarding their role as gene regulators. However, the biological processes in which most lncRNAs perform their function are still unclear. This study investigates the association between lncRNAs with host response to C. albicans using a public RNA-Seq dataset from lung samples of female C57BL/6J wild-type Mus musculus with induced C. albicans infection. The animals were exposed to the fungus for 24 h before sample collection. We selected lncRNAs and protein-coding genes related to the host immune response by combining the results from different computational approaches used for gene selection: differential expression gene analysis, co-expression genes network analysis, and machine learning-based gene selection. Using a guilt by association strategy, we inferred connections between 41 lncRNAs and 25 biological processes. Our results indicated that nine up-regulated lncRNAs were associated with biological processes derived from the response to wounding: 1200007C13Rik, 4833418N02Rik, Gm12840, Gm15832, Gm20186, Gm38037, Gm45774, Gm4610, Mir22hg, and Mirt1. Additionally, 29 lncRNAs were related to genes involved in immune response, while 22 lncRNAs were associated with processes related to reactive species production. These results support the participation of lncRNAs during C. albicans infection, and may contribute to new studies investigating lncRNA functions in the immune response.

Nanog organizes transcription bodies

The localization of transcriptional activity in specialized transcription bodies is a hallmark of gene expression in eukaryotic cells.^1-3 How proteins of the transcriptional machinery come together to form such bodies, however, is unclear. Here, we take advantage of two large, isolated, and long-lived transcription bodies that reproducibly form during early zebrafish embryogenesis to characterize the dynamics of transcription body formation. Once formed, these transcription bodies are enriched for initiating and elongating RNA polymerase II, as well as the transcription factors Nanog and Sox19b. Analyzing the events leading up to transcription, we find that Nanog and Sox19b cluster prior to transcription. The clustering of transcription factors is sequential; Nanog clusters first, and this is required for the clustering of Sox19b and the initiation of transcription. Mutant analysis revealed that both the DNA-binding domain as well as one of the two intrinsically disordered regions of Nanog are required to organize the two bodies of transcriptional activity. Taken together, our data suggest that the clustering of transcription factors dictates the formation of transcription bodies.

Long noncoding RNA lincsc5d regulates hepatic cholesterol synthesis by modulating sterol C5 desaturase in large yellow croaker

Although long noncoding RNA (lncRNA) plays a vital role in cholesterol metabolism, very little information is available in fish. Thus, a 10-week feeding experiment was performed to estimate the effects of lncRNA on cholesterol metabolism in large yellow croaker fed with fish oil (FO), soybean oil (SO), olive oil (OO), and palm oil (PO) diets. Results showed that fish fed with OO and PO diets had higher liver total cholesterol (TC) and cholesterol ester (CE) contents compared with fish fed with FO diets. Analysis of the KEGG pathway showed that the steroid biosynthesis pathway was enriched in comparisons FO vs SO, FO vs OO, and FO vs PO. Meanwhile, sterol C5 desaturase (SC5D), a cholesterol synthase, was up-regulated in the steroid biosynthesis pathway. SC5D was widely expressed in all tissues examined, and the highest expression of SC5D was detected in brain. More importantly, a novel lncRNA associated with sc5d gene was identified by RNA sequencing and named as lincsc5d. The tissue distribution of lincsc5d was similar to that of sc5d. A nuclear/cytoplasmic RNA separation assay showed that lincsc5d was a nucleus-enriched lncRNA. qRT-PCR results demonstrated that lincsc5d was markedly up-regulated in the SO, OO, and PO groups. Furthermore, the results of TC content and the lincsc5d and sc5d expression in hepatocytes agreed with in vivo results. In conclusion, this study indicated that vegetable oils, especially OO and PO, increased hepatic cholesterol levels by promoting cholesterol synthesis, and lncRNA lincsc5d and sc5d might be involved in cholesterol synthesis.

Throughout in vitro first spermatogenic wave: Next-generation sequencing gene expression patterns of fresh and cryopreserved prepubertal mice testicular tissue explants

INTRODUCTION

Suitable cryopreservation procedures of pre-pubertal testicular tissue associated with efficient culture conditions are crucial in the fields of fertility preservation and restoration. In vitro spermatogenesis remains a challenging technical procedure to undergo a complete spermatogenesis.The number of haploid cells and more specifically the spermatic yield produced in vitro in mice is still extremely low compared to age-matched in vivo controls and this procedure has never yet been successfully transferred to humans.

METHODS

To evaluate the impact of in vitro culture and freezing procedure, pre-pubertal testicular mice testes were directly cultured until day 4 (D4), D16 and D30 or cryopreserved by controlled slow freezing then cultured until D30. Testes composed of a panel of 6.5 dpp (days postpartum), 10.5 dpp, 22.5 dpp, and 36.5 dpp mice were used as in vivo controls. Testicular tissues were assessed by histological (HES) and immunofluorescence (stimulated by retinoic acid gene 8, STRA8) analyses. Moreover, a detailed transcriptome evaluation study has been carried out to study the gene expression patterns throughout the first in vitro spermatogenic wave.

RESULTS

Transcriptomic analyses reveal that cultured tissues expression profiles are almost comparable between D16 and D30; highlighting an abnormal kinetic throughout the second half of the first spermatogenesis during in vitro cultures. In addition, testicular explants have shown dysregulation of their transcriptomic profile compared to controls with genes related to inflammation response, insulin-like growth factor and genes involved in steroidogenesis.

DISCUSSION

The present work first shows that cryopreservation had very little impact on gene expression in testicular tissue, either directly after thawing or after 30 days in culture. Transcriptomic analysis of testis tissue samples is highly informative due to the large number of expressed genes and identified isoforms. This study provides a very valuable basis for future studies concerning in vitro spermatogenesis in mice.

Genome-wide analysis of long non-coding RNAs under diel light exhibits role in floral development and the circadian clock in Arabidopsis thaliana

The circadian clock is regulated by signaling networks that enhance a plant's ability to coordinate internal events with the external environment. In this study, we examine the rhythmic expression of long non-coding RNAs (lncRNAs) using multiple transcriptomes of Arabidopsis thaliana in the diel light cycle and integrated this information to have a better understanding of the functions of lncRNAs in regulating the circadian clock. We identified 968, 1050, and 998 lncRNAs at 8 h light, 16 h light and 8 h dark conditions, respectively. Among these, 423, 486, and 417 lncRNAs were uniquely present at 8 h light, 16 h light, and 8 h dark, respectively, whereas 334 lncRNAs were common under the three conditions. The specificity of identified lncRNAs under different light conditions was verified using qRT-PCR. The identified lncRNAs were less GC-rich and expressed at a significantly lower level than the mRNAs of protein-coding genes. In addition, we identified enriched motifs in lncRNA transcribing regions that were associated with light-responsive genes (SORLREP and SORLIP), flower development (AGAMOUS), and circadian clock (CCA1) under all three light conditions. We identified 10 and 12 different lncRNAs targeting different miRNAs with perfect and interrupted complementarity (endogenous target mimic). These predicted lncRNA-interacting miRNAs govern the function of a set of genes involved in the developmental process, reproductive structure development, gene silencing and transcription regulation. We demonstrated that the lncRNA transcribing regions were enriched for epigenetic marks such as H3.3, H3K4me2, H3K4me3, H4K16ac, H3K36ac, H3K56ac and depleted for heterochromatic (H3K9me2 and H3K27me1) and repressive (H3K27me3) histone modifications. Further, we found that hypermethylated genomic regions negatively correlated with lncRNA transcribing regions. Overall, our study showed that lncRNAs expressed corresponding to the diel light cycle are implicated in regulating the circadian rhythm and governing the developmental stage-specific growth.

Casting CRISPR-Cas13d to fish for microprotein functions in animal development

Protein coding genes were originally identified with sequence-based definitions that included a 100-codon cutoff to avoid annotating irrelevant open reading frames. However, many active proteins contain less than 100 amino acids. Indeed, functional genetics, ribosome profiling, and proteomic profiling have identified many short, translated open reading frames, including those with biologically active peptide products (microproteins). Yet, functions for most of these peptide products remain unknown. Because microproteins often act as key signals or fine-tune processes, animal development has already revealed functions for a handful of microproteins and provides an ideal context to uncover additional microprotein functions. However, many mRNAs during early development are maternally provided and hinder targeted mutagenesis approaches to characterize developmental microprotein functions. The recently established, RNA-targeting CRISPR-Cas13d system in zebrafish overcomes this barrier and produces potent knockdown of targeted mRNA, including maternally provided mRNA, and enables flexible, efficient interrogation of microprotein functions in animal development.

Genome-wide identification of antisense lncRNAs and their association with susceptibility to Flavobacterium psychrophilum in rainbow trout

Eukaryotic genomes encode long noncoding natural antisense transcripts (lncNATs) that have been increasingly recognized as regulatory members of gene expression. Recently, we identified a few antisense transcripts correlating in expression with immune-related genes. However, a systematic genome-wide analysis of lncNATs in rainbow trout is lacking. This study used 134 RNA-Seq datasets from five different projects to identify antisense transcripts. A total of 13,503 lncNATs were identified genome-wide. About 75% of lncNATs showed multiple exons compared to 36.5% of the intergenic lncRNAs. RNA-Seq datasets from resistant, control, and susceptible rainbow trout genetic lines with significant differences in survival rate following Flavobacterium psychrophilum (Fp) infection were analyzed to investigate the potential role of the lncNATs during infection. Twenty-four pairwise comparisons between the different genetic lines, infectious status, and time points revealed 581 differentially expressed (DE) lncNATs and 179 differentially used exons (DUEs). Most of the DE lncNATs strongly and positively correlated in expression with their corresponding sense transcripts across 24 RNA-Seq datasets. LncNATs complementary to genes related to immunity, muscle contraction, proteolysis, and iron/heme metabolism were DE following infection. LncNATs complementary to hemolysis-related genes were DE in the resistant fish compared to susceptible fish on day 5 post-infection, suggesting enhanced clearance of free hemoglobin (Hb) and heme and increased erythropoiesis. LncNATs complementary to hepcidin, a master negative regulator of the plasma iron concentration, were the most downregulated lncNATs on day 5 of bacterial infection in the resistant fish. Ninety-four DE lncNAT, including five complementary to hepcidin, are located within 26 QTL regions previously identified in association with bacterial cold water disease (BCWD) in rainbow trout. Collectively, lncNATs are involved in the molecular architecture of fish immunity and should be further investigated for potential applications in genomic selection and genetic manipulation in aquaculture.

Single-molecule Real-time (SMRT) Sequencing Facilitates Transcriptome Research and Genome Annotation of the Fish Sillago sinica

As a newly described Sillaginidae species, Chinese sillago (Sillago sinica) needs a better understanding of gene annotation information. In this study, we reported the first full-length transcriptome data of S. sinica using the PacBio isoform sequencing Iso-seq and a description of transcriptome structure analysis. A total of 454,979 high-quality full-length transcripts were obtained by single-molecule real-time (SMRT) sequencing, which was corrected by Illumina sequencing data. After that, 66,948 non-redundant full-length transcripts were generated after mapping to the reference genome of S. sinica, including 49 fusion isoforms and 9,250 novel isoforms. 63,459 isoforms were successfully annotated by one of the Nr, Nt, SwissProt, Pfam, KOG, GO, and KEGG databases. Additionally, 30,987 alternative polyadenylation (APA) sites, 451,867 alternative splicing (AS) events, 21,928 long non-coding RNAs (lncRNAs) and 12,911 transcription factors (TFs) were identified. The full-length transcripts of S. sinica would provide a precious resource for characterizing the transcriptome of S. sinica and for the further study of gene function and regulatory mechanism of this species.

Functional analysis of long non-coding RNAs involved in alkaline stress responses in wheat

Saline-alkali soil is a major environmental problem affecting crop productivity. One of the most effective approaches to combat it is to breed stress-tolerant plants through genetic engineering. Shanrong No. 4 (SR4) is an alkaline-tolerant cultivar of bread wheat (Triticum aestivum) derived from asymmetric somatic hybridization between the common wheat cultivar Jinan 177 (JN177) and tall wheatgrass. In this study, we aimed to explore the structure and function of alkalinity stress-responsive long non-coding RNAs (lncRNAs) in wheat. Sequencing was employed to identify the lncRNAs associated with stress tolerance and their corresponding targets. Approximately 19 000 novel lncRNA sequences were detected in SR4 and JN177. Upon exposure to alkaline stress, SR4 differentially expressed 5691 lncRNAs, whilst JN177 differentially expressed 5932. We selected five of them (L0760, L6247, L0208, L2098, and L3065) and generated seedlings of transiently knocked down strains using the virus-induced gene-silencing method. Knockdown of L0760 and L2098 caused the plants to exhibit sensitivity to alkaline stress, whereas knockdown of L6247, L0208, and L3065 increased the ability of plants to tolerate alkaline stress. We constructed lncRNA-miRNA-target-mRNA networks and alkali-response-related lncRNA-target-mRNA association networks to analyse the functions of lncRNAs. Collectively, our results demonstrate that lncRNAs may perform different roles under alkaline stress conditions.

When does hepatitis B virus meet long-stranded noncoding RNAs?

Hepatitis B virus (HBV) infection in humans and its associated diseases are long-standing problems. HBV can produce a large number of non-self-molecules during its life cycle, which acts as targets for innate immune recognition and initiation. Among these, interferon and its large number of downstream interferon-stimulated gene molecules are important early antiviral factors. However, the development of an effective antiviral immune response is not simple and depends not only on the delicate regulation of the immune response but also on the various mechanisms of virus-related immune escape and immune tolerance. Therefore, despite there being a relatively well-established consensus on the major pathways of the antiviral response and their component molecules, the complete clearance of HBV remains a challenge in both basic and clinical research. Long-noncoding RNAs (lncRNAs) are generally >200 bp in length and perform different functions in the RNA strand encoding the protein. As an important part of the IFN-inducible genes, interferon-stimulated lncRNAs are involved in the regulation of several HBV infection-related pathways. This review traces the basic elements of such pathways and characterizes the various recent targets of lncRNAs, which not only complement the regulatory mechanisms of pathways related to chronic HBV infection, fibrosis, and cancer promotion but also present with new potential therapeutic targets for controlling HBV infection and the malignant transformation of hepatocytes.

Noncoding RNAs in liver physiology and metabolic diseases

The liver holds central roles in detoxification, energy metabolism and whole-body homeostasis but can develop malignant phenotypes when being chronically overwhelmed with fatty acids and glucose. The global rise of metabolic-associated fatty liver disease (MAFLD) is already affecting a quarter of the global population. Pharmaceutical treatment options against different stages of MAFLD do not yet exist and several clinical trials against hepatic transcription factors and other proteins have failed. However, emerging roles of noncoding RNAs, including long (lncRNA) and short noncoding RNAs (sRNA), in various cellular processes pose exciting new avenues for treatment interventions. Actions of noncoding RNAs mostly rely on interactions with proteins, whereby the noncoding RNA fine-tunes protein function in a process termed riboregulation. The developmental stage-, disease stage- and cell type-specific nature of noncoding RNAs harbors enormous potential to precisely target certain cellular pathways in a spatio-temporally defined manner. Proteins interacting with RNAs can be categorized into canonical or non-canonical RNA binding proteins (RBPs) depending on the existence of classical RNA binding domains. Both, RNA- and RBP-centric methods have generated new knowledge of the RNA-RBP interface and added an additional regulatory layer. In this review, we summarize recent advances of how of RBP-lncRNA interactions and various sRNAs shape cellular physiology and the development of liver diseases such as MAFLD and hepatocellular carcinoma.

Genome-Wide Identification of Long Non-coding RNAs in Crassostrea gigas and Their Association with Heat Stress

Oysters face a complex and changeable environment in the intertidal zone. Heat stress is the main cause of their mass summer deaths. Several important genes are identified to be associated with heat response for oysters. However, regulation of these heat response genes in oysters remains largely unknown. In this study, 27 RNA-seq datasets are used to give a relatively comprehensive of long noncoding RNA (lncRNA) sets for C. gigas. Then, the differential expressed genes and lncRNAs are identified under heat stress. Among all the heat shock proteins (HSPs) and inhibitor of apoptosis proteins (IAPs) in the C. gigas genome, 25 heat shock proteins and 14 IAPs are differential expressed. The Gene Ontology analysis reveals that differential expressed genes (DEGs) are enriched in 6, 7, and 7 GO terms in cellular components, molecular function, and biological process, respectively. Within these terms, cellular response to stimulus is the most abundant term. Furthermore, the potential cis target of the differential expressed lncRNAs (DELs) are predicted to investigate their functions. Of the 394 DELs, there are 80 DELs being found to be corresponded to 113 protein-coding genes. Of them, eight HSPs are found to be regulated by their lncRNA regulators under heat stress. This work provides valuable resource of lncRNA and their regulatory roles under heat stress in C. gigas and gives new insights into adaptive evolution in marine mollusks.

Temporal Dynamic Analysis of Alternative Splicing During Embryonic Development in Zebrafish

Alternative splicing is pervasive in mammalian genomes and involved in embryo development, whereas research on crosstalk of alternative splicing and embryo development was largely restricted to mouse and human and the alternative splicing regulation during embryogenesis in zebrafish remained unclear. We constructed the alternative splicing atlas at 18 time-course stages covering maternal-to-zygotic transition, gastrulation, somitogenesis, pharyngula stages, and post-fertilization in zebrafish. The differential alternative splicing events between different developmental stages were detected. The results indicated that abundance alternative splicing and differential alternative splicing events are dynamically changed and remarkably abundant during the maternal-to-zygotic transition process. Based on gene expression profiles, we found splicing factors are expressed with specificity of developmental stage and largely expressed during the maternal-to-zygotic transition process. The better performance of cluster analysis was achieved based on the inclusion level of alternative splicing. The biological function analysis uncovered the important roles of alternative splicing during embryogenesis. The identification of isoform switches of alternative splicing provided a new insight into mining the regulated mechanism of transcript isoforms, which always is hidden by gene expression. In conclusion, we inferred that alternative splicing activation is synchronized with zygotic genome activation and discovered that alternative splicing is coupled with transcription during embryo development in zebrafish. We also unveiled that the temporal expression dynamics of splicing factors during embryo development, especially co-orthologous splicing factors. Furthermore, we proposed that the inclusion level of alternative splicing events can be employed for cluster analysis as a novel parameter. This work will provide a deeper insight into the regulation of alternative splicing during embryogenesis in zebrafish.

Comprehensive transcriptomic analysis of two RIL parents with contrasting salt responsiveness identifies polyadenylated and non-polyadenylated flower lncRNAs in chickpea

Salinity severely affects the yield of chickpea. Understanding the role of lncRNAs can shed light on chickpea salt tolerance mechanisms. However, because lncRNAs are encoded by multiple sites within the genome, their classification to reveal functional versatility at the transcriptional and the post-transcriptional levels is challenging. To address this, we deep sequenced 24 salt-challenged flower transcriptomes from two parental genotypes of a RIL population that significantly differ in salt tolerance ability. The transcriptomes for the first time included 12 polyadenylated and 12 non-polyadenylated RNA libraries to a sequencing depth of ~50 million reads. The ab initio transcriptome assembly comprised ~34 082 transcripts from three biological replicates of salt-tolerant (JG11) and salt-sensitive (ICCV2) flowers. A total of 9419 lncRNAs responding to salt stress were identified, 2345 of which were novel lncRNAs specific to chickpea. The expression of poly(A+) lncRNAs and naturally antisense transcribed RNAs suggest their role in post-transcriptional modification and gene silencing. Notably, 178 differentially expressed lncRNAs were induced in the tolerant genotype but repressed in the sensitive genotype. Co-expression network analysis revealed that the induced lncRNAs interacted with the FLOWERING LOCUS (FLC), chromatin remodelling and DNA methylation genes, thus inducing flowering during salt stress. Furthermore, 26 lncRNAs showed homology with reported lncRNAs such as COOLAIR, IPS1 and AT4, thus confirming the role of chickpea lncRNAs in controlling flowering time as a crucial salt tolerance mechanism in tolerant chickpea genotype. These robust set of differentially expressed lncRNAs provide a deeper insight into the regulatory mechanisms controlled by lncRNAs under salt stress.

Multiomic atlas with functional stratification and developmental dynamics of zebrafish cis-regulatory elements

Zebrafish, a popular organism for studying embryonic development and for modeling human diseases, has so far lacked a systematic functional annotation program akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created a central repository to store and process zebrafish developmental functional genomic data. Our data coordination center ( https://danio-code.zfin.org ) combines a total of 1,802 sets of unpublished and re-analyzed published genomic data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements throughout development, including classes with distinct features dependent on their activity in time and space. We delineated the distinct distance topology and chromatin features between regulatory elements active during zygotic genome activation and those active during organogenesis. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predicted functional relationships between them beyond sequence similarity, thus extending the utility of zebrafish developmental genomics to mammals.