Identification of Rhythmically Expressed LncRNAs in the Zebrafish Pineal Gland and Testis

SUDAZFLNC - a curated and searchable online database for zebrafish lncRNAs, mRNAs, miRNAs, and circadian expression profiles

The zebrafish (Danio rerio) has emerged as a model organism for investigating lncRNAs-driven fundamental biological processes, such as circadian rhythms, physiology, metabolism, and various diseases. While state-of-the-art sequencing technologies have identified an increasing number of lncRNAs in zebrafish, their annotations are far from complete. In this study, we collect 28,925 lncRNAs from both the published studies and our own RNA-seq analyses and establish a novel webserver-based database called SUDAZFLNC (https://sudarna.website/). The database, containing 28,925 lncRNAs, 25,432 mRNAs, and 368 miRNAs, provides several crucial features and annotations for the zebrafish RNAs, such as sequence identifiers (IDs), sequence length, hexamer score, coding probabilities, GO and KEGG annotations, and micropeptides. SUDAZFLNC also includes time-course expression profiles of 3288 lncRNAs, 25,432 mRNAs, and 342 miRNAs generated from our RNA-seq experiments, and 149, 4407, and 43 rhythmically expressed lncRNAs, mRNAs, and miRNAs, respectively. Based on the peak expression patterns, we classified these RNAs into morning RNAs, evening RNAs, and night RNAs. Users of the database can access the RNA sequences and their expression profiles by searching the corresponding IDs from the Graphical User Interface (GUI) of the database. The database supports several features to investigate RNA sequences and expression profiles, including BLAST, search of sequence and data, ID conversion, and RNA-RNA interaction prediction. This is the largest curated database of zebrafish RNAs and their expression profiles to date.

Thousands of oscillating LncRNAs in the mouse testis

The long noncoding RNAs (lncRNAs) are involved in numerous fundamental biological processes, including circadian regulation. Although recent studies have revealed insights into the functions of lncRNAs, how the lncRNAs regulate circadian rhythms still requires a deeper investigation. In this study, we generate two datasets of RNA-seq profiles of the mouse (Mus musculus) testis under light-dark (LD) cycle. The first dataset included 18,613 unannotated transcripts measured at 12 time points, each with duplicate samples, under LD conditions; while the second dataset included 21,414 unannotated transcripts measured at six time points, each with three replicates, under desynchronized and control conditions. We identified 5964 testicular lncRNAs in each dataset by BLASTing these transcripts against the known mouse lncRNAs from the NONCODE database. MetaCycle analyses were performed to identify 519, 475, and 494 rhythmically expressed mouse testicular lncRNAs in the 12-time-point dataset, the six-time-point control dataset, and the six-time-point desynchronized dataset, respectively. A comparison of the expression profiles of the lncRNAs under desynchronized and control conditions revealed that 427 rhythmically expressed lncRNAs from the control condition became arrhythmic under the desynchronized condition, suggesting a possible loss of rhythmicity. In contrast, 446 arrhythmic lncRNAs from the control condition became rhythmic under the desynchronized condition, suggesting a possible gain of rhythmicity. Interestingly, 48 lncRNAs were rhythmically expressed under both desynchronized and control conditions. These oscillating lncRNAs were divided into morning lncRNAs, evening lncRNAs, and night lncRNAs based on their time-course expression patterns. We interrogated the promoter regions of these rhythmically expressed mouse testicular lncRNAs to predict their possible regulation by the E-box, D-box, or RORE promoter motifs. GO and KEGG analyses were performed to identify the possible biological functions of these rhythmically expressed mouse testicular lncRNAs. Further, we conducted conservation analyses of the rhythmically expressed mouse testicular lncRNAs with lncRNAs from humans, rats, and zebrafish, and uncovered three mouse testicular lncRNAs conserved across these four species. Finally, we computationally predicted the conserved lncRNA-encoded peptides and their 3D structures from each of the four species. Taken together, our study revealed thousands of rhythmically expressed lncRNAs in the mouse testis, setting the stage for further computational and experimental validations.

Novel Insights into the Circadian Rhythms Based on Long Noncoding and Circular RNA Profiling

Circadian rhythm disorders pose major risks to human health and animal production activity, and the hypothalamus is the center of circadian rhythm regulation. However, the epigenetic regulation of circadian rhythm based on farm animal models has been poorly investigated. We collected chicken hypothalamus samples at seven time points in one light/dark cycle and performed long noncoding RNA (lncRNA), circular RNA (circRNA), and mRNA sequencing to detect biomarkers associated with circadian rhythm. We enhanced the comprehensive expression profiling of ncRNAs and mRNAs in the hypothalamus and found two gene sets (circadian rhythm and retinal metabolism) associated with the light/dark cycle. Noncoding RNA networks with circadian expression patterns were identified by differential expression and circadian analysis was provided that included 38 lncRNAs, 15 circRNAs, and 200 candidate genes. Three lncRNAs (ENSGALT00000098661, ENSGALT00000100816, and MSTRG.16980.1) and one circRNA (novel_circ_010168) in the ncRNA-mRNA regulatory network were identified as key molecules influencing circadian rhythm by regulating AOX1 in retinal metabolism. These ncRNAs were predicted to be related to pernicious anemia, gonadal, eye disease and other disorders in humans. Together, the findings of this study provide insights into the epigenetic mechanisms of circadian rhythm and reveal AOX1 as a promising target of circadian rhythm regulation.

Day-night and seasonal variation of human gene expression across tissues

Circadian and circannual cycles trigger physiological changes whose reflection on human transcriptomes remains largely uncharted. We used the time and season of death of 932 individuals from GTEx to jointly investigate transcriptomic changes associated with those cycles across multiple tissues. Overall, most variation across tissues during day-night and among seasons was unique to each cycle. Although all tissues remodeled their transcriptomes, brain and gonadal tissues exhibited the highest seasonality, whereas those in the thoracic cavity showed stronger day-night regulation. Core clock genes displayed marked day-night differences across multiple tissues, which were largely conserved in baboon and mouse, but adapted to their nocturnal or diurnal habits. Seasonal variation of expression affected multiple pathways, and it was enriched among genes associated with the immune response, consistent with the seasonality of viral infections. Furthermore, they unveiled cytoarchitectural changes in brain regions. Altogether, our results provide the first combined atlas of how transcriptomes from human tissues adapt to major cycling environmental conditions. This atlas may have multiple applications; for example, drug targets with day-night or seasonal variation in gene expression may benefit from temporally adjusted doses.

Long-stranded non-coding RNAs temporal-specific expression profiles reveal longissimus dorsi muscle development and intramuscular fat deposition in Tianzhu white yak

The process of muscle development and intramuscular fat (IMF) deposition is quite complex and controlled by both mRNAs and ncRNAs. Long-stranded non-coding RNAs (LncRNAs) are involved in various biological processes in mammals while also playing a critical role in muscle development and fat deposition. In the present study, RNA-Seq was used to comprehensively study the expression of lncRNAs and mRNAs during muscle development and intramuscular fat deposition in postnatal Tianzhu white yaks at three stages, including 6 mo of age (calve, n = 6), 30 mo of age (young cattle, n = 6) and 54 mo of age (adult cattle, n = 6). The results indicated that a total of 2,101 lncRNAs and 20,855 mRNAs were screened across the three stages, of which the numbers of differential expression (DE) lncRNAs and DE mRNAs were 289 and 1,339, respectively, and DE lncRNAs were divided into eight different expression patterns based on expression trends. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that some DE mRNAs overlapped with target genes of lncRNAs, such as NEDD4L, SCN3B, AGT, HDAC4, DES, MYH14, KLF15 (muscle development), ACACB, PCK2, LIPE, PIK3R1, PNPLA2, and MGLL (intramuscular fat deposition). These DE mRNAs were significantly enriched in critical muscle development and IMF deposition-related pathways and GO terms, such as AMPK signaling pathway, PI3K-Akt signaling pathway, PPAR signaling pathway, etc. In addition, lncRNA-mRNA co-expression network analysis revealed that six lncRNAs (MSTRG.20152.2, MSTRG.20152.3, XR_001351700.1, MSTRG.8190.1, MSTRG.4827.1, and MSTRG.11486.1) may play a major role in Tianzhu white yak muscle development and lipidosis deposition. Therefore, this study enriches the database of yak lncRNAs and could help to further explore the functions and roles of lncRNAs in different stages of muscle development and intramuscular fat deposition in the Tianzhu white yak.

Day-night and seasonal variation of human gene expression across tissues

Circadian and circannual cycles trigger physiological changes whose reflection on human transcriptomes remains largely uncharted. We used the time and season of death of 932 individuals from GTEx to jointly investigate transcriptomic changes associated with those cycles across multiple tissues. Overall, most variation across tissues during day-night and among seasons was unique to each cycle. Although all tissues remodeled their transcriptomes, brain and gonadal tissues exhibited the highest seasonality, whereas those in the thoracic cavity showed stronger day-night regulation. Core clock genes displayed marked day-night differences across multiple tissues, which were largely conserved in baboon and mouse, but adapted to their nocturnal or diurnal habits. Seasonal variation of expression affected multiple pathways and it was enriched among genes associated with the immune response, consistent with the seasonality of viral infections. Furthermore, they unveiled cytoarchitectural changes in brain regions. Altogether, our results provide the first combined atlas of how transcriptomes from human tissues adapt to major cycling environmental conditions.

Hundreds of LncRNAs Display Circadian Rhythmicity in Zebrafish Larvae

Long noncoding RNAs (lncRNAs) have been shown to play crucial roles in various life processes, including circadian rhythms. Although next generation sequencing technologies have facilitated faster profiling of lncRNAs, the resulting datasets require sophisticated computational analyses. In particular, the regulatory roles of lncRNAs in circadian clocks are far from being completely understood. In this study, we conducted RNA-seq-based transcriptome analysis of zebrafish larvae under both constant darkness (DD) and constant light (LL) conditions in a circadian manner, employing state-of-the-art computational approaches to identify approximately 3220 lncRNAs from zebrafish larvae, and then uncovered 269 and 309 lncRNAs displaying circadian rhythmicity under DD and LL conditions, respectively, with 30 of them are coexpressed under both DD and LL conditions. Subsequently, GO, COG, and KEGG pathway enrichment analyses of all these circadianly expressed lncRNAs suggested their potential involvement in numerous biological processes. Comparison of these circadianly expressed zebrafish larval lncRNAs, with rhythmically expressed lncRNAs in the zebrafish pineal gland and zebrafish testis, revealed that nine (DD) and twelve (LL) larval lncRNAs are coexpressed in the zebrafish pineal gland and testis, respectively. Intriguingly, among peptides encoded by these coexpressing circadianly expressed lncRNAs, three peptides (DD) and one peptide (LL) were found to have the known domains from the Protein Data Bank. Further, the conservation analysis of these circadianly expressed zebrafish larval lncRNAs with human and mouse genomes uncovered one lncRNA and four lncRNAs shared by all three species under DD and LL conditions, respectively. We also investigated the conserved lncRNA-encoded peptides and found one peptide under DD condition conserved in these three species and computationally predicted its 3D structure and functions. Our study reveals that hundreds of lncRNAs from zebrafish larvae exhibit circadian rhythmicity and should help set the stage for their further functional studies.

LINC00662: A new oncogenic lncRNA with great potential

LINC00662 is located on chromosome 19q11 and is 2085 bp long. It is a long noncoding RNA (lncRNA) newly discovered. LINC00662 expression is upregulated in at least 14 tumors. In addition, the upregulation of LINC00662 expression is also closely related to the poor prognosis of cancer patients and resistance to radiotherapy and chemotherapy. LINC00662 can act as a ceRNA of at least 8 miRNAs. By regulating these miRNAs and their downstream genes, LINC00662 participates in the regulation of four signaling pathways, including the extracellular signal-regulated kinase (ERK) signaling pathway, the Wnt/β-catenin signaling pathway, the Hippo signaling pathway, and the SMD signaling pathway. In addition, the abnormal upregulation of LINC00662 can promote the stem-like features of lung cancer cells. LINC00662 can reduce the promoter methylation level of s-adenosylmethionine (SAM)-dependent hepatocellular carcinoma (HCC)-promoting genes by regulating the MAT1A/SAM and AHCY/SAH axes, thereby promoting the activation of oncogenes. This article summarizes the molecular regulation mechanism of LINC00662 in cancer and the diagnostic and prognostic value of LINC00662 in cancer.