LncTar: a tool for predicting the RNA targets of long noncoding RNAs

Pan-cancer survey of lncRNA rewiring and functional alternation in tumor-infiltrating T cell by scLNC

Long non-coding RNAs (lncRNAs) have been reported to involve in diverse biological processes, including tumor immunity. Since lncRNAs are expressed with high cell-type specificity, investigation of lncRNAs at the single-cell level will unveil the cell-type-specific functions of lncRNAs. However, at the single-cell level, a systematic pan-cancer analysis of lncRNA functions in tumor immune microenvironments (TIMEs) remains lacking. Here, we performed pan-cancer single-cell profiling of lncRNA functions in TIMEs and developed a tool, scLNC, tailored for lncRNA functional characterization at the single-cell level. scLNC enabled the comparison of lncRNA function from the levels of lncRNA-mRNA pairs, lncRNA regulatory unit activity and unit function in a cell-type-specific manner. Applying scLNC, our analysis depicted the cross-tumor and tumor-specific lncRNA regulatory profiles in the T cell subtypes and revealed the new regulatory units that lncRNAs established in tumor-infiltrating T cells, particularly in the tumor-enriched T cells. We further characterized the activity and functional alternations of lncRNAs through their regulatory units. Overall, our findings suggested that lncRNAs played an important role in the regulation of cytokine production, cell activation and migration in tumor-enriched T cells and further in immunotherapy.

Co-expression analysis of lncRNA and mRNA suggests a role for ncRNA-mediated regulation of host-parasite interactions in primary skin lesions of patients with American tegumentary leishmaniasis

Leishmaniasis, caused by different Leishmania species, manifests as cutaneous or visceral forms. In the American continent, the cutaneous form is called American tegumentary leishmaniasis (ATL) and is primarily caused by Leishmania (Viannia) braziliensis. Mucosal leishmaniasis (ML), the most severe form of ATL, arises in approximately 20% of patients from a primary cutaneous lesion. Evidence indicates changes in overall expression patterns of mRNAs and lncRNAs of the host in response to Leishmania infection, with the parasite capable of modulating host immune response, which may contribute to disease progression. We evaluated whether the co-expression of lncRNAs and their putative target mRNAs in primary cutaneous lesions of patients with ATL could be associated with the development of ML. Previously available public RNA-Seq data from primary skin lesions of patients infected with L. braziliensis was employed. We identified 579 mRNAs and 46 lncRNAs differentially expressed in the primary lesion that subsequently progressed to mucosal disease. Co-expression analysis revealed 1,324 significantly correlated lncRNA-mRNA pairs. Among these, we highlight the positive correlation and trans-action between lncRNA SNHG29 and mRNA S100A8, both upregulated in the ML group. S100A8 and its heterodimeric partner S100A9 form a pro-inflammatory complex expressed by immune cells and seems to participate in host innate immune response processes of infection. These findings expand the knowledge of the Leishmania-host interaction and indicate that the expression of lncRNAs in the primary cutaneous lesion could regulate mRNAs and play roles in disease progression.

Interactome of Long Non-Coding RNAs: Transcriptomic Expression Patterns and Shaping Cancer Cell Phenotypes

RNA biology has gained extensive recognition in the last two decades due to the identification of novel transcriptomic elements and molecular functions. Cancer arises, in part, due to the accumulation of mutations that greatly contribute to genomic instability. However, the identification of differential gene expression patterns of wild-type loci has exceeded the boundaries of mutational study and has significantly contributed to the identification of molecular mechanisms that drive carcinogenic transformation. Non-coding RNA molecules have provided a novel avenue of exploration, providing additional routes for evaluating genomic and epigenomic regulation. Of particular focus, long non-coding RNA molecule expression has been demonstrated to govern and direct cellular activity, thus evidencing a correlation between aberrant long non-coding RNA expression and the pathological transformation of cells. lncRNA classification, structure, function, and therapeutic utilization have expanded cancer studies and molecular targeting, and understanding the lncRNA interactome aids in defining the unique transcriptomic signatures of cancer cell phenotypes.

GKLOMLI: a link prediction model for inferring miRNA-lncRNA interactions by using Gaussian kernel-based method on network profile and linear optimization algorithm

BACKGROUND

The limited knowledge of miRNA-lncRNA interactions is considered as an obstruction of revealing the regulatory mechanism. Accumulating evidence on Human diseases indicates that the modulation of gene expression has a great relationship with the interactions between miRNAs and lncRNAs. However, such interaction validation via crosslinking-immunoprecipitation and high-throughput sequencing (CLIP-seq) experiments that inevitably costs too much money and time but with unsatisfactory results. Therefore, more and more computational prediction tools have been developed to offer many reliable candidates for a better design of further bio-experiments.

METHODS

In this work, we proposed a novel link prediction model based on Gaussian kernel-based method and linear optimization algorithm for inferring miRNA-lncRNA interactions (GKLOMLI). Given an observed miRNA-lncRNA interaction network, the Gaussian kernel-based method was employed to output two similarity matrixes of miRNAs and lncRNAs. Based on the integrated matrix combined with similarity matrixes and the observed interaction network, a linear optimization-based link prediction model was trained for inferring miRNA-lncRNA interactions.

RESULTS

To evaluate the performance of our proposed method, k-fold cross-validation (CV) and leave-one-out CV were implemented, in which each CV experiment was carried out 100 times on a training set generated randomly. The high area under the curves (AUCs) at 0.8623 ± 0.0027 (2-fold CV), 0.9053 ± 0.0017 (5-fold CV), 0.9151 ± 0.0013 (10-fold CV), and 0.9236 (LOO-CV), illustrated the precision and reliability of our proposed method.

CONCLUSION

GKLOMLI with high performance is anticipated to be used to reveal underlying interactions between miRNA and their target lncRNAs, and deciphers the potential mechanisms of the complex diseases.

Long non-coding RNA HIF1A-As2 and MYC form a double-positive feedback loop to promote cell proliferation and metastasis in KRAS-driven non-small cell lung cancer

Lung cancer is the leading cause of cancer-related deaths worldwide. KRAS is the main oncogenic driver in lung cancer that can be activated by gene mutation or amplification, but whether long non-coding RNAs (lncRNAs) regulate its activation remains unknown. Through gain and loss of function approaches, we identified that lncRNA HIF1A-As2, a KRAS-induced lncRNA, is required for cell proliferation, epithelial-mesenchymal transition (EMT) and tumor propagation in non-small cell lung cancer (NSCLC) in vitro and in vivo. Integrative analysis of HIF1A-As2 transcriptomic profiling reveals that HIF1A-As2 modulates gene expression in trans, particularly regulating transcriptional factor genes including MYC. Mechanistically, HIF1A-As2 epigenetically activates MYC by recruiting DHX9 on MYC promoter, consequently stimulating the transcription of MYC and its target genes. In addition, KRAS promotes HIF1A-As2 expression via the induction of MYC, suggesting HIF1A-As2 and MYC form a double-regulatory loop to strengthen cell proliferation and tumor metastasis in lung cancer. Inhibition of HIF1A-As2 by LNA GapmeR antisense oligonucleotides (ASO) significantly improves sensitization to 10058-F4 (a MYC-specific inhibitor) and cisplatin treatment in PDX and KRAS^LSLG12D-driven lung tumors, respectively.

Cannabinoid enhancement of lncRNA MMP25-AS1/MMP25 interaction reduces neutrophil infiltration and intestinal epithelial injury in HIV/SIV infection

Intestinal epithelial barrier dysfunction, a hallmark of HIV/SIV infection, persists despite viral suppression by combination antiretroviral therapy (cART). Emerging evidence suggests a critical role for long noncoding RNAs (lncRNAs) in maintaining epithelial homeostasis. We simultaneously profiled lncRNA/mRNA expression exclusively in colonic epithelium (CE) of SIV-infected rhesus macaques (RMs) administered vehicle (VEH) or Δ-9-tetrahydrocannabinol (THC). Relative to controls, fewer lncRNAs were up- or downregulated in CE of THC/SIV compared with VEH/SIV RMs. Importantly, reciprocal expression of the natural antisense lncRNA MMP25-AS1 (up 2.3-fold) and its associated protein-coding gene MMP25 (attracts neutrophils by inactivating alpha-1 anti-trypsin/SERPINA1) (down 2.2-fold) was detected in CE of THC/SIV RMs. Computational analysis verified 2 perfectly matched complementary regions and an energetically stable (normalized binding free energy = -0.2626) MMP25-AS1/MMP25 duplex structure. MMP25-AS1 overexpression blocked IFN-γ-induced MMP25 mRNA and protein expression in vitro. Elevated MMP25 protein expression in CE of VEH/SIV but not THC/SIV RMs was associated with increased infiltration by myeloperoxidase/CD11b++ neutrophils (transendothelial migration) and epithelial CD47 (transepithelial migration) expression. Interestingly, THC administered in combination with cART increased MMP25-AS1 and reduced MMP25 mRNA/protein expression in jejunal epithelium of SIV-infected RMs. Our findings demonstrate that MMP25-AS1 is a potentially unique epigenetic regulator of MMP25 and that low-dose THC can reduce neutrophil infiltration and intestinal epithelial injury potentially by downregulating MMP25 expression through modulation of MMP25-AS1.

A hybrid approach to assess the structural impact of long noncoding RNA mutations uncovers key NEAT1 interactions in colorectal cancer

Long noncoding RNAs (lncRNAs) are emerging players in cancer and they entail potential as prognostic biomarkers or therapeutic targets. Earlier studies have identified somatic mutations in lncRNAs that are associated with tumor relapse after therapy, but the underlying mechanisms behind these associations remain unknown. Given the relevance of secondary structure for the function of some lncRNAs, some of these mutations may have a functional impact through structural disturbance. Here, we examined the potential structural and functional impact of a novel A > G point mutation in NEAT1 that has been recurrently observed in tumors of colorectal cancer patients experiencing relapse after treatment. Here, we used the nextPARS structural probing approach to provide first empirical evidence that this mutation alters NEAT1 structure. We further evaluated the potential effects of this structural alteration using computational tools and found that this mutation likely alters the binding propensities of several NEAT1-interacting miRNAs. Differential expression analysis on these miRNA networks shows upregulation of Vimentin, consistent with previous findings. We propose a hybrid pipeline that can be used to explore the potential functional effects of lncRNA somatic mutations.

Adaptation of the Porcine Pituitary Transcriptome, Spliceosome and Editome during Early Pregnancy

The physiological mechanisms of the porcine reproduction are relatively well-known. However, transcriptomic changes and the mechanisms accompanying transcription and translation processes in various reproductive organs, as well as their dependence on hormonal status, are still poorly understood. The aim of this study was to gain a principal understanding of alterations within the transcriptome, spliceosome and editome occurring in the pituitary of the domestic pig (Sus scrofa domestica L.), which controls basic physiological processes in the reproductive system. In this investigation, we performed extensive analyses of data obtained by high-throughput sequencing of RNA from the gilts' pituitary anterior lobes during embryo implantation and the mid-luteal phase of the estrous cycle. During analyses, we obtained detailed information on expression changes of 147 genes and 43 long noncoding RNAs, observed 784 alternative splicing events and also found the occurrence of 8729 allele-specific expression sites and 122 RNA editing events. The expression profiles of the selected 16 phenomena were confirmed by PCR or qPCR techniques. As a final result of functional meta-analysis, we acquired knowledge regarding intracellular pathways that induce changes in the processes accompanying transcription and translation regulation, which may induce modifications in the secretory activity of the porcine adenohypophyseal cells.

Full-length transcriptome profiling for fruit development in Diospyros oleifera using nanopore sequencing

OBJECTIVES

Diospyros oleifera, one of the most economically important Diospyros species, is an ideal model for studying the fruit development of persimmon. While, the lack of whole-transcriptome has hindered the complex transcriptional regulation mechanisms of sugar and tannin during fruit development.

DATA DESCRIPTION

We applied Oxford Nanopore Technologies to six developmental stage of fruit from D. oleifera for use in transcriptome sequencing. As a result of full-length transcriptome sequencing, 55.87 Gb of clean data were generated. After mapping onto the reference genome of D. oleifera, 51,588 full-length collapsing transcripts, including 2,727 new gene loci and 43,223 transcripts, were obtained. Comprehensively annotated, 38,086 of new transcripts were functional annotation, and 972 lncRNAs, 7,159 AS events were predicted. Here, we released the transcriptome database of D. oleifera at different stage of fruit development,which will provide a fundamention of to investigatethe transcript structure, variants and evolution of persimmon.

The potential role of N6-methyladenosine modification of LncRNAs in contributing to the pathogenesis of chronic glomerulonephritis

BACKGROUND

Increasing evidence indicates that N6-methyladenosine (m6A) modification of mRNAs has been shown to play a critical role in the occurrence and development of many diseases, while little is known about m6A modification in long non-coding RNAs (LncRNAs). Our study aims to investigate the potential functions of LncRNA m6A modifications in lipopolysaccharide (LPS)-induced mouse mesangial cells (MMCs), providing us with a new perspective on the molecular mechanisms of chronic glomerulonephritis (CGN) pathogenesis.

METHODS

Differentially methylated LncRNAs were identified by Methylated RNA immunoprecipitation sequencing (MeRIP-seq). LncRNA-mRNA and LncRNA-associated LncRNA-miRNA-mRNA (CeRNA) networks were constructed by bioinformatics analysis. Furthermore, we utilized gene ontology (GO) and pathway enrichment analyses (KEGG) to explore target genes from co-expression networks. In addition, the total level of m6A RNA methylation and expression of methyltransferase and pro-inflammatory cytokines were detected by the colorimetric quantification method and western blot, respectively. Cell viability and cell cycle stage were detected by cell counting kit-8 (CCK-8) and flow cytometry.

RESULTS

In total, 1141 differentially m6A-methylated LncRNAs, including 529 hypermethylated LncRNAs and 612 hypomethylated LncRNAs, were determined by MeRIP-seq. The results of GO and KEGG analysis revealed that the target mRNAs were mainly enriched in signal pathways, such as the NF-kappa B signaling pathway, MAPK signaling pathway, Toll-like receptor signaling pathway, and apoptosis signaling pathway. In addition, higher METTL3 expression was found in CGN kidney tissues using the GEO database. METTL3 knockdown in MMC cells drastically reduced the levels of m6A RNA methylation, pro-inflammatory cytokines IL6 and TNF-α, and inhibited cell proliferation and cycle progression.

CONCLUSIONS

Our findings provide a basis and novel insight for further investigations of m6A modifications in LncRNAs for the pathogenesis of CGN.

Comparative transcriptome analysis of adult worker bees under short-term heat stress

High temperature affects behavior, physiology, survival, and the expression of related genes in adult honeybees. Apis mellifera is the common pollinator in greenhouse and is susceptible to high temperature stress. To further explore the molecular basis related to heat stress, we compared the transcriptome profiles of adult worker bees at 25 and 45°C, and detected the expression patterns of some differentially expressed genes (DEGs) in different tissues by q RT-PCR. Differential expression analysis showed that 277 DEGs were identified, including 167 genes upregulated and 110 genes downregulated after heat stress exposure in adult worker bees. In GO enrichment analysis, DEGs were mostly enriched for protein folding, unfold protein binding, and heme binding terms. Protein processing in endoplasmic reticulum and longevity regulating pathway-multiple species were significantly enriched in KEGG. The expression levels of 16 DEGs were consistent with the transcriptome results. The expression patterns of 9 DEGs in different tissues revealed high levels in the thorax, which was supposed that the thorax may be the most important part in the response to heat stress. This study provided valuable data for exploring the function of heat resistance-related genes.

Genome-wide analysis of long non-coding RNAs in sugar beet (Beta vulgaris L.) under drought stress

Drought stress is one of the most severe abiotic stresses that restrict global crop production. Long non-coding RNAs (lncRNAs) have been proved to play a key role in response to drought stress. However, genome-wide identification and characterization of drought-responsive lncRNAs in sugar beet is still lacking. Thus, the present study focused on analyzing lncRNAs in sugar beet under drought stress. We identified 32017 reliable lncRNAs in sugar beet by strand-specific high-throughput sequencing. A total of 386 differentially expressed lncRNAs (DElncRNAs) were found under drought stress. The most significantly upregulated and downregulated lncRNAs were TCONS_00055787 (upregulated by more than 6000 fold) and TCONS_00038334 (downregulated by more than 18000 fold), respectively. Quantitative real-time PCR results exhibited a high concordance with RNA sequencing data, which conformed that the expression patterns of lncRNAs based on RNA sequencing were highly reliable. In addition, we predicted 2353 and 9041 transcripts that were estimated to be the cis- and trans-target genes of the drought-responsive lncRNAs. As revealed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the target genes of DElncRNAs were significantly enriched in organelle subcompartment, thylakoid, endopeptidase activity, catalytic activity, developmental process, lipid metabolic process, RNA polymerase activity, transferase activity, flavonoid biosynthesis and several other terms associated with abiotic stress tolerance. Moreover, 42 DElncRNAs were predicted as potential miRNA target mimics. LncRNAs have important effects on plant adaptation to drought conditions through the interaction with protein-encoding genes. The present study leads to greater insights into lncRNA biology and offers candidate regulators for improving the drought tolerance of sugar beet cultivars at the genetic level.

Co-expression analysis of lncRNA and mRNA identifies potential adipogenesis regulatory non-coding RNAs involved in the transgenerational effects of tributyltin

The obesity epidemic is considered a global public health crisis, with an increase in caloric intake, sedentary lifestyles and/or genetic predispositions as contributing factors. Although the positive energy balance is one of the most significant causes of obesity, recent research has linked early exposure to Endocrine-Disrupting Chemicals (EDCs) such as the obesogen tributyltin (TBT) to the disease epidemic. In addition to their actions on the hormonal profile, EDCs can induce long-term changes in gene expression, possibly due to changes in epigenetic patterns. Long non-coding RNAs (lncRNAs) are epigenetic mediators that play important regulatory roles in several biological processes, through regulation of gene transcription and/or translation. In this study, we explored the differential expression of lncRNAs in gonadal white adipose tissue samples from adult male C57BL/6J F4 generation, female C57BL/6J offspring exposed (F0 generation) to 50 nM TBT or 0.1% DMSO (control of vehicle) via drinking water provided during pregnancy and lactation, analyzing RNA-seq data from a publicly available dataset (GSE105051). A total of 74 lncRNAs were differentially expressed (DE), 22 were up-regulated and 52 were down-regulated in the group whose F4 ancestor was exposed in utero to 50nM TBT when compared to those exposed to 0.1% DMSO (control). Regulation of DE lncRNAs and their potential partner genes in gonadal white adipose tissue of mice ancestrally exposed to EDC TBT may be related to the control of adipogenesis, as pathway enrichment analyses showed that these gene partners are mainly involved in the metabolism of lipids and glucose and in insulin-related pathways, which are essential for obesity onset and control.

Recent advances in predicting lncRNA-disease associations based on computational methods

Mutations in and dysregulation of long non-coding RNAs (lncRNAs) are closely associated with the development of various human complex diseases, but only a few lncRNAs have been experimentally confirmed to be associated with human diseases. Predicting new potential lncRNA-disease associations (LDAs) will help us to understand the pathogenesis of human diseases and to detect disease markers, as well as in disease diagnosis, prevention and treatment. Computational methods can effectively narrow down the screening scope of biological experiments, thereby reducing the duration and cost of such experiments. In this review, we outline recent advances in computational methods for predicting LDAs, focusing on LDA databases, lncRNA/disease similarity calculations, and advanced computational models. In addition, we analyze the limitations of various computational models and discuss future challenges and directions for development.

Identification of birch lncRNAs and mRNAs responding to salt stress and characterization of functions of lncRNA

Long noncoding RNAs (lncRNAs) are important in abiotic stress tolerance. Here, we identified salt-responsive genes and lncRNAs in the roots and leaves of Betula platyphylla Suk. (birch), and characterized their lncRNAs functions. In total, 2660 mRNAs and 539 lncRNAs responding to salt treatment were identified using RNA-seq. The salt-responsive genes were substantially enriched in 'cell wall biogenesis' and 'wood development' in the roots and were enriched in 'photosynthesis' and 'response to stimulus' in the leaves. Meanwhile, the potential target genes of the salt-responsive lncRNAs in roots and leaves were both enriched in 'nitrogen compound metabolic process' and 'response to stimulus'. We further built a method for quickly identifying abiotic stress tolerance of lncRNAs, which employed transient transformation for overexpression and knock-down of the lncRNA, enabling gain- and loss-of-function analysis. Using this method, 11 randomly selected salt-responsive lncRNAs were characterized. Among them, six lncRNAs confer salt tolerance, two lncRNAs confer salt sensitivity, and the other three lncRNAs are not involved in salt tolerance. In addition, a lncRNA, LncY1, was further characterized, which improves salt tolerance by regulating two transcription factors, BpMYB96 and BpCDF3. Taken together, our results suggested that lncRNAs play important roles in the salt response of birch plants.

Novel lncRNA-prader willi/angelman region RNA, SNRPN neighbour (PWARSN) aggravates tubular epithelial cell pyroptosis by regulating TXNIP via dual way in diabetic kidney disease

OBJECTIVES

Elevated thioredoxin-interacting protein (TXNIP)-induced pyroptosis contributes to the pathology of diabetic kidney disease (DKD). However, the molecular mechanisms in dysregulated TXNIP in DKD remain largely unclear.

MATERIALS AND METHODS

Transcriptomic analysis identified a novel long noncoding RNA-Prader Willi/Angelman region RNA, SNRPN neighbour (PWARSN)-which was highly expressed in a proximal tubular epithelial cell (PTEC) under high glucose conditions. We focused on revealing the functions of PWARSN in regulating TXNIP-mediated pyroptosis in PTECs by targeting PWARSN expression via lentivirus-mediated overexpression and CRISPR-Cas9-based knockout in vitro and overexpressing PWARSN in the renal cortex by AAV-9 targeted injection in vivo. A number of molecular techniques disclosed the mechanisms of PWARSN in regulating TXNIP induced-pyroptosis in DKD.

RESULTS

TXNIP-NOD-like receptor thermal protein domain associated protein 3 (NLRP3) inflammasome and PTEC pyroptosis were activated in the renal tubules of patients with DKD and in diabetic mice. Then we explored that PWARSN enhanced TXNIP-driven PTECs pyroptosis in vitro and in vivo. Mechanistically, cytoplasmic PWARSN sponged miR-372-3p to promote TXNIP expression. Moreover, nuclear PWARSN interacted and facilitated RNA binding motif protein X-linked (RBMX) degradation through ubiquitination, resulting in the initiation of TXNIP transcription by reducing H3K9me3-enrichment at the TXNIP promoter. Further analysis indicated that PWARSN might be a potential biomarker for DKD.

CONCLUSIONS

These findings illustrate distinct dual molecular mechanisms for PWARSN-modulated TXNIP and PTECs pyroptosis in DKD, presenting PWARSN as a promising therapeutic target for DKD.

Comprehensive Analysis of Whole-Transcriptome Profiles in Response to Acute Hypersaline Challenge in Chinese Razor Clam Sinonovacula constricta

The Chinese razor clam (Sinonovacula constricta) is an important for Chinese aquaculture marine bivalve that naturally occurs across intertidal and estuarine areas subjected to significant changes in salinity level. However, the information on the molecular mechanisms related to high salinity stress in the species remain limited. In this study, nine gill samples of S. constricta treated with 20, 30, and 40 ppt salinity for 24 h were used for whole-transcriptome RNA sequencing, and a regulatory network of competing endogenous RNAs (ceRNAs) was constructed to better understand the mechanisms responsible for adaptation of the species to high salinity. A total of 83,262 lncRNAs, 52,422 mRNAs, 2890 circRNAs, and 498 miRNAs were identified, and 4175 of them displayed differential expression pattern among the three groups examined. The KEGG analyses of differentially expressed RNAs evidenced that amino acid synthesis and membrane transport were the dominant factors involved in the adaptation of the Chinese razor clam to acute salinity increase, while lipid metabolism and signaling played only a supporting role. In addition, lncRNA/circRNA-miRNA-mRNA regulatory networks (ceRNA network) showed clearly regulatory relationships among different RNAs. Moreover, the expression of four candidate genes, including tyrosine aminotransferase (TAT), hyaluronidase 4 (HYAL4), cysteine sulfinic acid decarboxylase (CSAD), and ∆¹-pyrroline-5-carboxylate synthase (P5CS) at different challenge time were detected by qRT-PCR. The expression trend of TAT and HYAL4 was consistent with that of the ceRNA network, supporting the reliability of established network. The expression of TAT, CSAD, and P5CS were upregulated in response to increased salinity. This might be associated with increased amino acid synthesis rate, which seems to play an essential role in adaptation of the species to high salinity stress. In contrast, the expression level of HYAL4 gene decreased in response to elevated salinity level, which is associated with reduction Hyaluronan hydrolysis to help maintain water in the cell. Our findings provide a very rich reference for understanding the important role of ncRNAs in the salinity adaptation of shellfish. Moreover, the acquired information may be useful for optimization of the artificial breeding of the Chinese razor clam under aquaculture conditions.

Full-length transcriptome sequencing reveals the molecular mechanism of monoterpene and sesquiterpene biosynthesis in Cinnamomum burmannii

Essential oil of Cinnamomum burmannii is rich in monoterpenes and sesquiterpenes and is widely used in cosmetics and medicines. Knowledge about the enzymes that catalyze the formation of monoterpenes and sesquiterpenes in C. burmannii is insufficient. Therefore, anatomy observation of C. burmannii at the four developmental stages (7 days, CBS1; 14 days, CBS2; 21 days, CBS3, and 28 days, CBS4) were conducted to elucidate the origins of essential oil production. Twelve full-length transcriptomes of C. burmannii leaves at the four stages were generated using Oxford Nanopore Technologies. GC-MS analysis revealed 15 monoterpene and sesquiterpenes dramatically increased from CBS1 to CBS4. A weighted correlation network analysis (WGCNA) in association and differentially expressed genes across four developmental stages were performed. A total of 44 differentially expressed genes (DEGs) were involved in terpenoid syntheses during leaf development. Among them, the DEGs of the mevalonate acid (MVA) pathway were predominantly expressed at CBS1, while those of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway showed increased expression from CBS2 to CBS4. Besides, fourteen genes were associated with monoterpene synthesis and nine with sesquiterpene synthesis. Functions of these DEGs were further predicted with regard to gene expression profile and phylogenetic relationship with those characterized in previous studies. In addition, 922 long noncoding RNAs (lncRNAs) were detected, of which twelve were predicted to regulate monoterpene and sesquiterpene biosynthesis. The present study provided new insights the molecular mechanisms of monoterpenoid and sesquiterpenoid syntheses of C. burmannii.

Genome-wide screening and characterization of long noncoding RNAs involved in flowering/bolting of Lactuca sativa

BACKGROUND

Lettuce (Lactuca sativa L.) is considered the most important vegetable in the leafy vegetable group. However, bolting affects quality, gives it a bitter taste, and as a result makes it inedible. Bolting is an event induced by the coordinated effects of various environmental factors and endogenous genetic components. Although bolting/flowering responsive genes have been identified in most sensitive and non-sensitive species, non-coding RNA molecules like long non-coding RNAs (lncRNAs) have not been investigated in lettuce. Hence, in this study, potential long non-coding RNAs that regulate flowering /bolting were investigated in two lettuce strains S24 (resistant strain) and S39 (susceptible strain) in different flowering times to better understand the regulation of lettuce bolting mechanism. For this purpose, we used two RNA-seq datasets to discover the lncRNA transcriptome profile during the transition from vegetative to reproductive phase.

RESULTS

For identifying unannotated transcripts in these datasets, a 7-step pipeline was employed to filter out these transcripts and terminate with 293 novel lncRNAs predicted by PLncPRO and CREMA. These transcripts were then utilized to predict cis and trans flowering-associated targets and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Computational predictions of target gene function showed the involvement of putative flowering-related genes and enrichment of the floral regulators FLC, CO, FT, and SOC1 in both datasets. Finally, 17 and 18 lncRNAs were proposed as competing endogenous target mimics (eTMs) for novel and known lncRNA miRNAs, respectively.

CONCLUSION

Overall, this study provides new insights into lncRNAs that control the flowering time of plants known for bolting, such as lettuce, and opens new windows for further study.

A long non-coding RNA PelncRNA1 is involved in Phyllostachys edulis response to UV-B stress

Phyllostachys edulis (moso bamboo) is China's most widespread bamboo species, with significant economic and ecological values. Long non-coding RNA (lncRNA) is a type of regulatory RNA that is longer than 200 nucleotides and incapable of encoding proteins, and is frequently involved in regulating biotic and abiotic stress and plant development. However, the biological functions of lncRNA in moso bamboo are unknown. In this study, a lncRNA (named PelncRNA1) differentially expressed following UV-B treatment was discovered in the whole transcriptome sequencing database of moso bamboo. The target genes were filtered and defined by correlation analysis of PelncRNA1 and gene expression pattern. The expression levels of PelncRNA1 and its target genes were verified using qRT-PCR. The results demonstrated that the expression levels of PelncRNA1 and its target genes increased during UV-B treatment. In Arabidopsis transgenic seedlings and moso bamboo protoplasts, PelncRNA1 was discovered to influence the expression of its target genes when overexpressed. In addition, transgenic Arabidopsis showed higher tolerance to UV-B stress. These results suggest that PelncRNA1 and its target genes are involved in the response of moso bamboo to UV-B stress. The novel findings would contribute to our understanding of how lncRNAs regulate the response to abiotic stresses in moso bamboo.

Rumen Epithelial Development- and Metabolism-Related Genes Regulate Their Micromorphology and VFAs Mediating Plateau Adaptability at Different Ages in Tibetan Sheep

The rumen is an important hallmark organ of ruminants and plays an important role in the metabolism and immune barrier of Tibetan sheep on the Plateau. However, there are few studies on rumen development and metabolism regulation in Tibetan sheep at different ages. Here, we comprehensively analyzed the immune function, fermentation function, rumen epithelial micromorphology and transcriptome profile of Tibetan sheep at different ages. The results showed that the concentration of IgG decreased and the concentration of IgM increased with age (p < 0.05), and the highest concentration of IgA was observed at 1.5 and 3.5 years of age. In terms of rumen fermentation characteristics, VFAs of 4-month-old lambs were the highest, followed by VFAs and NH3-N of Tibetan sheep at 3.5 years of age. Hematoxylin-eosin staining and transmission electron microscopy section examination of rumen epithelial tissue showed that the rumen papilla width increased with age (p < 0.001), the thickness of the stratum corneum decreased, the cells in the stratum corneum showed accelerated migration and the thickness of the rumen muscle layer increased (p < 0.001). Desmosomal junctions between the layers of rumen epithelium increased at 1.5 and 3.5 years old, forming a compact barrier structure, and the basal layer had more mitochondria involved in the regulation of energy metabolism. RNA-seq analysis revealed that a total of 1006 differentially expressed genes (DEGs) were identified at four ages. The DEGs of Tibetan sheep aged 4 months and 6 years were mainly enriched in the oxidation−reduction process and ISG15-protein conjugation pathway. The 1.5 and 3.5-year-olds were mainly enriched in skeletal muscle thin filament assembly, mesenchyme migration and the tight junction pathway. WGCNA showed that DEGs related to rumen microbiota metabolite VFAs and epithelial morphology were enriched in “Metabolism of xenobiotics by cytochrome P450, PPAR signaling pathway, Butanoate metabolism pathways” and participated in the regulation of rumen epithelial immune and fermentation metabolism functions of Tibetan sheep at different ages. This study systematically revealed the regulatory mechanism of rumen epithelial development and metabolism in the plateau adaptation of Tibetan sheep, providing a new approach for the study of plateau adaptation.

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.

Identification and differential analysis of noncoding RNAs in response to drought in Phyllostachys aureosulcata f. spectabilis

The role of noncoding RNAs (ncRNAs) in plant resistance to abiotic stresses is increasingly being discovered. Drought stress is one of the most common stresses that affecting plant growth, and high intensity drought has a significant impact on the normal growth of plants. In this study, a high-throughput sequencing was performed on plant tissue samples of Phyllostachys aureosulcata f. spectabilis C. D. Chu et C. S. Chao by drought treatment for 0, 2, 4 and 6 days. The sequencing results were analysed bioinformatically. We detected 336,946 RNAs among all 12 samples, including 192,098 message RNAs (mRNAs), 142,761 long noncoding RNAs (lncRNAs), 1,670 circular RNAs (circRNAs), and 417 microRNAs (miRNAs). We detected 2,419 differentially expressed (DE) ncRNAs, including 213 DE circRNAs, 2,088 DE lncRNAs and 118 DE miRNAs. Then, we used Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) to functionally predict DE ncRNAs. The results showed that most DE ncRNAs are involved in the response to drought stress, mainly in biochemical reactions involved in some metabolites, as well as in organelle activities. In addition, we validated two random circRNAs and demonstrated their circularity. We also found a stable internal reference gene available for Phyllostachys aureosulcata f. spectabilis and validated the accuracy of this experiment by quantitative real-time polymerase chain reaction (qRT-PCR).

Response of Ruminal Microbiota-Host Gene Interaction to High-Altitude Environments in Tibetan Sheep

Altitude is the main external environmental pressure affecting the production performance of Tibetan sheep, and the adaptive evolution of many years has formed a certain response mechanism. However, there are few reports on the response of ruminal microbiota and host genomes of Tibetan sheep to high-altitude environments. Here, we conducted an integrated analysis of volatile fatty acids (VFAs), microbial diversity (16S rRNA), epithelial morphology, and epithelial transcriptome in the rumen of Tibetan sheep at different altitudes to understand the changes in ruminal microbiota−host interaction in response to high altitude. The differences in the nutritional quality of forage at different altitudes, especially the differences in fiber content (ADF/NDF), led to changes in rumen VFAs of Tibetan sheep, in which the A/P value (acetic acid/propionic acid) was significantly decreased (p < 0.05). In addition, the concentrations of IgA and IgG in Middle-altitude (MA) and High-altitude Tibetan sheep (HA) were significantly increased (p < 0.05), while the concentrations of IgM were significantly increased in MA (p < 0.05). Morphological results showed that the width of the rumen papilla and the thickness of the basal layer increased significantly in HA Tibetan sheep (p < 0.05). The 16S rRNA analysis found that the rumen microbial diversity of Tibetan sheep gradually decreased with increasing altitude, and there were some differences in phylum- and genus-level microbes at the three altitudes. RDA analysis found that the abundance of the Rikenellaceae RC9 gut group and the Ruminococcaceae NK4A214 group increased with altitudes. Furthermore, a functional analysis of the KEGG microbial database found the “lipid metabolism” function of HA Tibetan sheep to be significantly enriched. WGCNA revealed that five gene modules were enriched in “energy production and conversion”, “lipid transport and metabolism”, and “defense mechanisms”, and cooperated with microbiota to regulate rumen fermentation and epithelial immune barrier function, so as to improve the metabolism and immune level of Tibetan sheep at high altitude.

Noncoding RNAs responsive to nitric oxide and their protein-coding gene targets shed light on root hair formation in Arabidopsis thaliana

An overview of the total Arabidopsis thaliana transcriptome, described previously by our research group, pointed some noncoding RNA (ncRNA) as participants in the restoration of hair-root phenotype in A. thaliana rhd6 mutants, leading us to a deeper investigation. A transcriptional gene expression profiling of seedling roots was performed aiming to identify ncRNA responsive to nitric oxide (GSNO) and auxin (IAA), and their involvement in root hair formation in the rhd6 null mutant. We identified 3,631 ncRNAs, including new ones, in A. thaliana and differential expression (DE) analysis between the following: 1) GSNO-treated rhd6 vs. untreated rhd6, 2) IAA-treated rhd6 vs. untreated rhd6, 3) GSNO-treated rhd6 vs. IAA-treated rhd6, and 4) WS-2 vs. untreated rhd6 detected the greatest number of DE genes in GSNO-treated rhd6. We detected hundreds of in silico interactions among ncRNA and protein-coding genes (PCGs), highlighting MIR5658 and MIR171 precursors highly upregulated in GSNO-treated rhd6 and wild type, respectively. Those ncRNA interact with many DE PCGs involved in hormone signaling, cell wall development, transcription factors, and root hair formation, becoming candidate genes in cell wall modulation and restoration of root hair phenotype by GSNO treatment. Our data shed light on how GSNO modulates ncRNA and their PCG targets in A. thaliana root hair formation.

Population transcriptomic analysis identifies the comprehensive lncRNAs landscape of spike in wheat (Triticum aestivum L.)

BACKGROUND

Long noncoding RNAs (lncRNAs) are emerging as the important regulators involving in growth and development as well as stress response in plants. However, current lncRNA studies were mainly performed at the individual level and the significance of it is not well understood in wheat.

RESULTS

In this study, the lncRNA landscape of wheat spike was characterized through analysing a total of 186 spike RNA-seq datasets from 93 wheat genotypes. A total of 35,913 lncRNAs as well as 1,619 lncRNA-mRNA pairs comprised of 443 lncRNAs and 464 mRNAs were obtained. Compared to coding genes, these lncRNAs displayed rather low conservation among wheat and other gramineous species. Based on re-sequencing data, the genetic variations of these lncRNA were investigated and obvious genetic bottleneck were found on them during wheat domestication process. Furthermore, 122 lncRNAs were found to act as ceRNA to regulate endogenous competition. Finally, association and co-localization analysis of the candidate lncRNA-mRNA pairs identified 170 lncRNAs and 167 target mRNAs significantly associated with spike-related traits, including lncRNA.127690.1/TraesCS2A02G518500.1 (PMEI) and lncRNA.104854.1/TraesCS6A02G050300.1 (ATG5) associated with heading date and spike length, respectively.

CONCLUSIONS

This study reported the lncRNA landscape of wheat spike through the population transcriptome analysis, which not only contribute to better understand the wheat evolution from the perspective of lncRNA, but also lay the foundation for revealing roles of lncRNA playing in spike development.

Full-length transcriptome revealed the accumulation of polyunsaturated fatty acids in developing seeds of Plukenetia volubilis

Background

Plukenetia volubilis is cultivated as a valuable oilseed crop, and its mature seeds are rich in polyunsaturated fatty acids (FAs), which are widely used in food and pharmaceutical industries. Recently, next-generation sequencing (NGS) transcriptome studies in P. volubilis indicated that some candidate genes were involved in oil biosynthesis. The NGS were inaccuracies in assembly of some candidate genes, leading to unknown errors in date analyses. However, single molecular real-time (SMRT) sequencing can overcome these assembled errors. Unfortunately, this technique has not been reported in P. volubilis.

Methods

The total oil content of P. volubilis seed (PVS) was determined using Soxhlet extraction system. The FA composition were analyzed by gas chromatography. Combining PacBio SMRT and Illumina technologies, the transcriptome analysis of developing PVS was performed. Functional annotation and differential expression were performed by BLAST software (version 2.2.26) and RSEM software (version 1.2.31), respectively. The lncRNA-targeted transcripts were predicted in developing PVS using LncTar tool.

Results

By Soxhlet extraction system, the oil content of superior plant-type (SPT) was 13.47% higher than that of inferior plant-type (IPT) at mature PVS. The most abundant FAs were C18:2 and C18:3, among which C18:3 content of SPT was 1.11-fold higher than that of IPT. Combined with PacBio and Illumina platform, 68,971 non-redundant genes were obtained, among which 7,823 long non-coding RNAs (lncRNAs) and 7,798 lncRNA-targeted genes were predicted. In developing seed, the expressions of 57 TFs showed a significantly positive correlation with oil contents, including WRI1-like1, LEC1-like1, and MYB44-like. Comparative analysis of expression profiles between SPT and IPT implied that orthologs of FAD3, PDCT, PDAT, and DAGT2 were possibly important for the accumulation of polyunsaturated FAs. Together, these results provide a reference for oil biosynthesis of P. volubilis and genetic improvement of oil plants.

Integrative Investigation of Root-Related mRNAs, lncRNAs and circRNAs of “Muscat Hamburg” (Vitis vinifera L.) Grapevine in Response to Root Restriction through Transcriptomic Analyses

Root restriction is a physical and ecological cultivation mode which restricts plant roots into a limited container to regulate vegetative and reproduction growth by reshaping root architecture. However, little is known about related molecular mechanisms. To uncover the root-related regulatory network of endogenous RNAs under root restriction cultivation (referred to RR), transcriptome-wide analyses of mRNAs, long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) involved in root development were performed. During root development, RR treatment had a positive effect on root weight, typically, young roots were significantly higher than conventional cultivation (referred to NR) treatment, suggesting that root architecture reconstruction under RR was attributed to the vigorous induction into lateral roots. Furthermore, a total of 26,588 mRNAs, 1971 lncRNAs, and 2615 circRNAs were identified in root of annual “Muscat Hamburg” grapevine by the transcriptomic analyses. The expression profile of mRNAs, lncRNAs and circRNA were further confirmed by the quantitative real-time PCR (RT-qPCR). Gene ontology enrichment analysis showed that a majority of the differentially expressed mRNAs, lncRNAs and circRNAs were enriched into the categories of cellular process, metabolic process, cell part, binding, and catalytic activity. In addition, the regulatory network of endogenous RNAs was then constructed by the prediction of lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA network, implying that these RNAs play significant regulatory roles for root architecture shaping in response to root restriction. Our results, for the first time, the regulatory network of competitive endogenous RNAs (ceRNAs) functions of lncRNA and circRNA was integrated, and a basis for studying the potential functions of non-coding RNAs (ncRNAs) during root development of grapevine was provided.

LncRNA-mRNA modules involved in goat rumen development: Insights from genome-wide transcriptome profiling

The rumen is an essential digestive and absorption organ of ruminants. During fetal life, lactation, and post-weaning period, goat rumen undergoes drastic morphological and metabolic-functional changes triggered by potential regulated genes and non-coding RNA molecules. As the essential regulatory factors, long non-coding RNAs (lncRNAs) have vital functions in various biological activities. However, their roles during rumen development are still poorly explored in ruminants. To explore the genome-wide expression profiles of lncRNAs and mRNAs in the goat rumens, we generated 5,007 lncRNAs and 19,738 mRNAs identified during the fetal and prepubertal stages by the high-throughput RNA sequencing. Notably, 365 lncRNAs and 2,877 mRNAs were considered to be differentially expressed. The weighted gene co-expression network analysis and functional analysis were performed to explore the regulatory roles of those differentially expressed molecules. The cis-and trans-target genes of differently expressed lncRNAs were enriched for pathways related to focal adhesion, cGMP-PKG signaling pathway, alpha-linolenic acid metabolism, arachidonic acid metabolism, and fat digestion and absorption. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes analyses showed that the differently expressed genes mainly participated in mitotic cytokinesis, desmosome, fatty acid degradation, cell adhesion molecules, and fatty acid metabolism. The prediction of lncRNA-mRNA interaction networks further revealed transcripts potentially involved in rumen development. The present study profiles a global overview of lncRNAs and mRNAs during rumen development. Our findings provide valuable resources for genetic regulation and molecular mechanisms of rumen development in ruminants.

Transcriptome sequencing reveals differences between leydig cells and sertoli cells of yak

In this study, we detected the expression of mRNAs, lncRNAs, and miRNAs in primary cultured leydig cells (LCs) and sertoli cells (SCs) of yak by RNA sequencing technology. A total of 84 differently expression mRNAs (DEmRNAs) (LCs vs. SCs: 15 up and 69 down), 172 differently expression lncRNAs (DElncRNAs) (LCs vs. SCs: 36 up and 136 down), and 90 differently expression miRNAs (DEmiRNAs) (LCs vs. SCs: 72 up and 18 down) were obtained between the two types of cells. GO enrichment and KEGG analysis indicated that the differential expression genes (DEGs) were more enriched in the regulation of actin cytoskeleton, Rap1/MAPK signaling pathway, steroid biosynthesis, focal adhesion, and pathways associated with metabolism. Targeted regulation relationship pairs of 3β-HSD and MSTRG.54630.1, CNTLN and MSTRG.19058.1, BRCA2 and MSTRG.28299.4, CA2 and novel-miR-148, and ceRNA network of LAMC3-MSTRG.68870.1- bta-miR-7862/novel-miR-151/novel-miR-148 were constructed by Cytoscape software. In conclusion, the differences between LCs and SCs were mainly reflected in steroid hormone synthesis, cell proliferation and metabolism, and blood-testicular barrier (BTB) dynamic regulation, and 3β-HSD, CNTLN, BRCA2, CA2, and LAMC3 may be the key factors causing these differences, which may be regulated by ncRNAs. This study provides a basic direction for exploring the differential regulation of LCs and SCs by ncRNAs.

B, Iquebal MA, Bhattacharya RC, Rai A, Kumar D. Genome wide identification of lncRNAs and circRNAs having regulatory role in fruit shelf life in health crop cucumber (Cucumis sativus L.)

Cucumber is an extremely perishable vegetable; however, under room conditions, the fruits become unfit for consumption 2-3 days after harvesting. One natural variant, DC-48 with an extended shelf-life was identified, fruits of which can be stored up to 10-15 days under room temperature. The genes involved in this economically important trait are regulated by non-coding RNAs. The study aims to identify the long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) by taking two contrasting genotypes, DC-48 and DC-83, at two different fruit developmental stages. The upper epidermis of the fruits was collected at 5 days and 10 days after pollination (DAP) for high throughput RNA sequencing. The differential expression analysis was performed to identify differentially expressed (DE) lncRNAs and circRNAs along with the network analysis of lncRNA, miRNA, circRNA, and mRNA interactions. A total of 97 DElncRNAs were identified where 18 were common under both the developmental stages (8 down regulated and 10 upregulated). Based on the back-spliced reads, 238 circRNAs were found to be distributed uniformly throughout the cucumber genomes with the highest numbers (71) in chromosome 4. The majority of the circRNAs (49%) were exonic in origin followed by inter-genic (47%) and intronic (4%) origin. The genes related to fruit firmness, namely, polygalacturonase, expansin, pectate lyase, and xyloglucan glycosyltransferase were present in the target sites and co-localized networks indicating the role of the lncRNA and circRNAs in their regulation. Genes related to fruit ripening, namely, trehalose-6-phosphate synthase, squamosa promoter binding protein, WRKY domain transcription factors, MADS box proteins, abscisic stress ripening inhibitors, and different classes of heat shock proteins (HSPs) were also found to be regulated by the identified lncRNA and circRNAs. Besides, ethylene biosynthesis and chlorophyll metabolisms were also found to be regulated by DElncRNAs and circRNAs. A total of 17 transcripts were also successfully validated through RT PCR data. These results would help the breeders to identify the complex molecular network and regulatory role of the lncRNAs and circRNAs in determining the shelf-life of cucumbers.

Phloridzin Reveals New Treatment Strategies for Liver Fibrosis

Liver fibrosis is an urgent public health problem which is difficult to resolve. However, various drugs for the treatment of liver fibrosis in clinical practice have their own problems during use. In this study, we used phloridzin to treat hepatic fibrosis in the CCl4-induced C57/BL6N mouse model, which was extracted from lychee core, a traditional Chinese medicine. The therapeutic effect was evaluated by biochemical index detections and ultrasound detection. Furthermore, in order to determine the mechanism of phloridzin in the treatment of liver fibrosis, we performed high-throughput sequencing of mRNA and lncRNA in different groups of liver tissues. The results showed that compared with the model group, the phloridzin-treated groups revealed a significant decrease in collagen deposition and decreased levels of serum alanine aminotransferase, aspartate aminotransferase, laminin, and hyaluronic acid. GO and KEGG pathway enrichment analysis of the differential mRNAs was performed and revealed that phloridzin mainly affects cell ferroptosis. Gene co-expression analysis showed that the target genes of lncRNA were obvious in cell components such as focal adhesions, intercellular adhesion, and cell–substrate junctions and in metabolic pathways such as carbon metabolism. These results showed that phloridizin can effectively treat liver fibrosis, and the mechanism may involve ferroptosis, carbon metabolism, and related changes in biomechanics.

In Silico Identification of lncRNAs Regulating Sperm Motility in the Turkey (Meleagris gallopavo L.)

Long non-coding RNAs (lncRNAs) are transcripts not translated into proteins with a length of more than 200 bp. LncRNAs are considered an important factor in the regulation of countless biological processes, mainly through the regulation of gene expression and interactions with proteins. However, the detailed mechanism of interaction as well as functions of lncRNAs are still unclear and therefore constitute a serious research challenge. In this study, for the first time, potential mechanisms of lncRNA regulation of processes related to sperm motility in turkey were investigated and described. Customized bioinformatics analysis was used to detect and identify lncRNAs, and their correlations with differentially expressed genes and proteins were also investigated. Results revealed the expression of 863 new/unknown lncRNAs in ductus deferens, testes and epididymis of turkeys. Moreover, potential relationships of the lncRNAs with the coding mRNAs and their products were identified in turkey reproductive tissues. The results obtained from the OMICS study may be useful in describing and characterizing the way that lncRNAs regulate genes and proteins as well as signaling pathways related to sperm motility.

Whole Transcriptome Sequencing Reveals Drought Resistance-Related Genes in Upland Cotton

China, particularly the cotton-growing province of Xinjiang, is experiencing acute agricultural water shortages, stifling the expansion of the cotton sector. Discovering drought resistance genes in cotton and generating high-quality, drought-resistant cotton varieties through molecular breeding procedures are therefore critical to the cotton industry’s success. The drought-resistant cotton variety Xinluzhong No. 82 and the drought-sensitive cotton variety Kexin No. 1 were utilised in this study to uncover a batch of drought-resistant candidate genes using whole transcriptome sequencing. The following are the key research findings: A competing endogenous RNA network (ceRNA) was built using complete transcriptional sequencing to screen the core genes in the core pathway, and two drought-related candidate genes were discovered. It was found that γ-aminobutyric acid aminotransferase (GhGABA-T, Gohir.A11G156000) was upregulated at 0 h vs. 12 h and downregulated at 12 h vs. 24 h. L-Aspartate oxidase (GhAO, Gohir.A07G220600) was downregulated at 0 h vs. 12 h and upregulated at 12 h vs. 24 h. GABA-T is analogous to a pyridoxal phosphate-dependent transferase superfamily protein (POP2) in Arabidopsis thaliana and influences plant drought resistance by controlling γ-aminobutyric acid (GABA) concentration. The analogue of GhAO in A. thaliana is involved in the early steps of nicotinamide adenine dinucleotide (NAD) production as well as in plant antioxidant responses. This study revealed that gene expression regulatory networks can be used for rapid screening of reliable drought resistance genes and then utilised to validate gene function.

Full-Length Transcriptional Analysis of the Same Soybean Genotype With Compatible and Incompatible Reactions to Heterodera glycines Reveals Nematode Infection Activating Plant Defense Response

Full-length transcriptome sequencing with long reads is a powerful tool to analyze transcriptional and post-transcriptional events; however, it has not been applied on soybean (Glycine max). Here, a comparative full-length transcriptome analysis was performed on soybean genotype 09-138 infected with soybean cyst nematode (SCN, Heterodera glycines) race 4 (SCN4, incompatible reaction) and race 5 (SCN5, compatible reaction) using Oxford Nanopore Technology. Each of 9 full-length samples collected 8 days post inoculation with/without nematodes generated an average of 6.1 GB of clean data and a total of 65,038 transcript sequences. After redundant transcripts were removed, 1,117 novel genes and 41,096 novel transcripts were identified. By analyzing the sequence structure of the novel transcripts, a total of 28,759 complete open reading frame (ORF) sequences, 5,337 transcription factors, 288 long non-coding RNAs, and 40,090 novel transcripts with function annotation were predicted. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of differentially expressed genes (DEGs) revealed that growth hormone, auxin-activated signaling pathway and multidimensional cell growth, and phenylpropanoid biosynthesis pathway were enriched by infection with both nematode races. More DEGs associated with stress response elements, plant-hormone signaling transduction pathway, and plant-pathogen interaction pathway with more upregulation were found in the incompatible reaction with SCN4 infection, and more DEGs with more upregulation involved in cell wall modification and carbohydrate bioprocess were detected in the compatible reaction with SCN5 infection when compared with each other. Among them, overlapping DEGs with a quantitative difference was triggered. The combination of protein-protein interaction with DEGs for the first time indicated that nematode infection activated the interactions between transcription factor WRKY and VQ (valine-glutamine motif) to contribute to soybean defense. The knowledge of the SCN-soybean interaction mechanism as a model will present more understanding of other plant-nematode interactions.

Identification of Long Noncoding RNAs Involved in Eyelid Pigmentation of Hereford Cattle

Several ocular pathologies in cattle, such as ocular squamous cell carcinoma and infectious keratoconjunctivitis, have been associated with low pigmentation of the eyelids. The main objective of this study was to analyze the transcriptome of eyelid skin in Hereford cattle using strand-specific RNA sequencing technology to characterize and identify long noncoding RNAs (lncRNAs). We compared the expression of lncRNAs between pigmented and unpigmented eyelids and analyzed the interaction of lncRNAs and putative target genes to reveal the genetic basis underlying eyelid pigmentation in cattle. We predicted 4,937 putative lncRNAs mapped to the bovine reference genome, enriching the catalog of lncRNAs in Bos taurus. We found 27 differentially expressed lncRNAs between pigmented and unpigmented eyelids, suggesting their involvement in eyelid pigmentation. In addition, we revealed potential links between some significant differentially expressed lncRNAs and target mRNAs involved in the immune response and pigmentation. Overall, this study expands the catalog of lncRNAs in cattle and contributes to a better understanding of the biology of eyelid pigmentation.

Multi-Omics Reveals That the Rumen Transcriptome, Microbiome, and Its Metabolome Co-regulate Cold Season Adaptability of Tibetan Sheep

Tibetan sheep can maintain a normal life and reproduce in harsh environments under extreme cold and lack of nutrition. However, the molecular and metabolic mechanisms underlying the adaptability of Tibetan sheep during the cold season are still unclear. Hence, we conducted a comprehensive analysis of rumen epithelial morphology, epithelial transcriptomics, microbiology and metabolomics in a Tibetan sheep model. The results showed that morphological structure of rumen epithelium of Tibetan sheep in cold season had adaptive changes. Transcriptomics analysis showed that the differential genes were primarily enriched in the PPAR signaling pathway (ko03320), legionellosis (ko05134), phagosome (ko04145), arginine and proline metabolism (ko00330), and metabolism of xenobiotics by cytochrome P450 (ko00980). Unique differential metabolites were identified in cold season, such as cynaroside A, sanguisorbin B and tryptophyl-valine, which were mainly enriched in arachidonic acid metabolism, arachidonic acid metabolism and linolenic acid metabolism pathways, and had certain correlation with microorganisms. Integrated transcriptome-metabolome-microbiome analysis showed that epithelial gene-GSTM3 expression was upregulated in the metabolism of xenobiotics by the cytochrome P450 pathway during the cold season, leading to the downregulation of some harmful metabolites; TLR5 gene expression was upregulated and CD14 gene expression was downregulated in the legionellosis pathway during the cold season. This study comprehensively described the interaction mechanism between the rumen host and microbes and their metabolites in grazing Tibetan sheep during the cold season. Rumen epithelial genes, microbiota and metabolites act together in some key pathways related to cold season adaptation.

Functional investigation of lncRNAs and target cytochrome P450 genes related to spirotetramat resistance in Aphis gossypii Glover

BACKGROUND

Spirotetramat is a tetramic acid derivative insecticide with novel modes of action for controlling Aphis gossypii Glover in the field. Previous studies have shown that long noncoding RNAs (lncRNAs) and cytochrome P450 monooxygenases (P450s) are involved in the detoxification process. However, the functions of lncRNAs in regulating P450 gene expression in spirotetramat resistance in A. gossypii are unknown.

RESULTS

In this study, we found CYP4CJ1, CYP6CY7 and CYP6CY21 expression levels to be significantly upregulated in a spirotetramat-resistant (SR) strain compared with a susceptible (SS) strain. Furthermore, knockdown of CYP4CJ1, CYP6CY7 and CYP6CY21 increased nymph and adult mortality in the SR strain following exposure to spirotetramat. Drosophila ectopically expressing CYP380C6, CYP4CJ1, CYP6DA2, CYP6CY7 and CYP6CY21 showed significantly decreased mortality after spirotetramat exposure, and CYP380C6, CYP4CJ1 and CYP6CY21 are putative targets of six lncRNAs. Silencing of lncRNAs MSTRG.36649.2/5 and MSTRG.71880.1 changed CYP6CY21 and CYP380C6 expression, altering the sensitivity of the SR strain to spirotetramat. Moreover, MSTRG.36649.2/5 did not compete for microRNA (miRNA) binding to regulate CYP6CY21 expression.

CONCLUSION

Our results confirm that CYP380C6, CYP4CJ1, CYP6DA2, CYP6CY7 and CYP6CY21 are potentially involved in the development of spirotetramat resistance in A. gossypii, and MSTRG.36649.2/5 and MSTRG.71880.1 probably regulate CYP6CY21 and CYP380C6 expression other than through the "sponge effect" of competing for miRNA binding. Our results provide a favorable molecular basis for studying cotton aphid P450 genes and lncRNA functions in spirotetramat resistance development.

LINC01431 Promotes Histone H4R3 Methylation to Impede HBV Covalently Closed Circular DNA Transcription by Stabilizing PRMT1

Covalently closed circular DNA (cccDNA) is the transcriptional template of hepatitis B virus (HBV), which interacts with both host and viral proteins to form minichromosome in the nucleus and is resistant to antiviral agents. Identification of host factors involved in cccDNA transcriptional regulation is expected to prove a new venue for HBV therapy. Recent evidence suggests the involvement of long noncoding RNAs (lncRNAs) in mediating the interaction of host factors with various viruses, however, lncRNAs that HBV targets and represses cccDNA transcription have not been fully elucidated. Here, the authors identified LINC01431 as a novel host restriction factor for HBV transcription. Mechanically, LINC01431 competitively bound with type I protein arginine methyltransferase (PRMT1) to block the HBx-mediated PRMT1 ubiquitination and degradation. Consequently, LINC01431 increased the occupancy of PRMT1 on cccDNA, leading to enhanced H4R3me2a modification and reduced acetylation of cccDNA-bound histones, thereby repressing cccDNA transcription. In turn, to facilitate viral replication, HBV transcriptionally repressed LINC01431 expression by HBx-mediated repression of transcription factor Zinc fingers and homeoboxes 2 (ZHX2). Collectively, the study demonstrates LINC01431 as a novel epigenetic regulator of cccDNA minichromosome and highlights a feedback loop of HBx-LINC01431-PRMT1 in HBV replication, which provides potential therapeutic targets for HBV treatment.

Research on lncRNA related to drought resistance of Shanlan upland rice

Background

Drought has become the major abiotic stress that causes losses in rice yields and consequently is one of the main environmental factors threatening food security. Long non-coding RNA (lncRNA) is known to play an important role in plant response to drought stress, while the mechanisms of competing endogenous RNA (ceRNA) in drought resistance in upland rice have been rarely reported.

Results

In our study, a total of 191 lncRNAs, 2115 mRNAs and 32 miRNAs (microRNAs) were found by strand-specific sequencing and small RNA sequencing to be differentially expressed in drought-stressed rice. Functional analysis of results indicate that they play important roles in hormone signal transduction, chlorophyll synthesis, protein synthesis and other pathways. Construction of a ceRNA network revealed that MSTRG.28732.3 may interact with miR171 in the chlorophyll biosynthesis pathway and affect the ability of plants to withstand drought stress by regulating Os02g0662700, Os02g0663100 and Os06g0105350. The accuracy of the regulatory network was verified by qRT-PCR.

Conclusion

Our results provide a theoretical basis for future studies on the potential function of lncRNA in plant drought resistance, and they provide new genetic resources for drought-resistant rice breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08546-0.

HIF1A promotes miR-210/miR-424 transcription to modulate the angiogenesis in HUVECs and HDMECs via sFLT1 under hypoxic stress

Angiogenesis is a critical process during human skin wound healing. However, hypoxia might lead to the dysfunction of the cellular interplay of endothelial cells and subcutaneous fibroblasts, resulting in the deregulation of angiogenesis. HIF1A is a key regulatory of the recovery of intracellular homeostasis under hypoxia. In the present study, the detailed role and mechanism of HIF1A in the angiogenesis under hypoxia were investigated. Via bioinformatic analyses on microarray profiles (GSE1041 and GSE17944), solube fms-related tyrosine kinase 1 (sFLT1, also known as sVEGFR1) and miR-210/miR-424 might be involved in HIF1A function on the angiogenesis under hypoxia in human umbilical vascular endothelium cells (HUVECs) and human dermal microvascular endothelial cells (HDMECs). In the present study, we identified sFLT1 as a downregulated gene in response to hypoxia and HIF1A overexpression in HUVECs and HDMECs. sFLT1 overexpression inhibited the capacity of migration and angiogenesis and significantly reversed the inducible effects of HIF1A on the migration and angiogenesis in both cell lines. miR-210 and miR-424 were upregulated by hypoxia and targeted sFLT1 3'-UTR to negatively modulate its expression. HIF1A modulated sFLT1 expression, VEGF signaling, and the migration and angiogenesis in HUVECs and HDMECs via miR-210/miR-424. Regarding the molecular mechanism, HIF1A bound the promoter region of miR-210 and miR-424 to activate their transcription, while miR-210/miR-424 bound sFLT1 3'-UTR to suppress its expression. In summary, HIF1A/miR-210/miR-424/sFLT1 axis modulates the angiogenesis in HUVECs and HDMECs upon hypoxic condition via VEGF signaling.

Integrative Analysis of Long- and Short-Read Transcriptomes Identify the Regulation of Terpenoids Biosynthesis Under Shading Cultivation in Oenanthe javanica

Water dropwort (Oenanthe javanica) is a popular vegetable with high nutritional value and distinctive flavor. The flavor is mainly correlate with the biosynthesis of terpenoids. Shading cultivation was used to improve the flavor in the production of water dropwort. However, the changes of terpenoids and the genes involved in terpenoids biosynthesis under shading treatment remains unclear. In this study, the long- and short-reads transcriptomes of water dropwort were constructed. In total, 57,743 non-redundant high-quality transcripts were obtained from the transcriptome. 28,514 SSRs were identified from non-redundant transcripts and the mono-nucleotide repeats were the most abundant SSRs. The lncRNAs of water dropwort were recognized and their target genes were predicted. The volatile compound contents in petioles and leaf blades of water dropwort were decreased after the shading treatment. The DEGs analysis was performed to identify the terpenoids biosynthesis genes. The results indicated that 5,288 DEGs were differentially expressed in petiole, of which 22 DEGs were enriched in the terpenoids backbone biosynthesis pathway. A total of 12 DEGs in terpenoids biosynthesis pathway were selected and further verified by qRT-PCR assay, demonstrating that the terpenoids biosynthesis genes were down-regulated under shading treatment. Here, the full-length transcriptome was constructed and the regulatory genes related to terpenoids biosynthesis in water dropwort were also investigated. These results will provide useful information for future researches on functional genomics and terpenoids biosynthesis mechanism in water dropwort.

Comparative analysis of long noncoding RNA and mRNA expression provides insights into adaptation to hypoxia in Tibetan sheep

Tibetan sheep have lived on the Qinghai-Tibetan Plateau for thousands of years and have good adaptability to the hypoxic environment and strong disease resistance. However, the molecular mechanism by which Tibetan sheep adapt to this extreme environment, especially the role of genetic regulation, is still unknown. Emerging evidence suggests that long noncoding RNAs (lncRNAs) participate in the regulation of a diverse range of biological processes. To explore the potential lncRNAs involved in the adaptation to high-altitude hypoxia of Tibetan sheep, we analysed the expression profile of lncRNAs and mRNAs in the liver and lung tissues of sheep using comparative transcriptome analysis between four Tibetan sheep populations (high altitude) and one Hu sheep population (low altitude). The results showed a total of 7848 differentially expressed (DE) lncRNA transcripts, and 22,971 DE mRNA transcripts were detected by pairwise comparison. The expression patterns of selected mRNAs and lncRNAs were validated by qRT-PCR, and the results correlated well with the transcriptome data. Moreover, the functional annotation analysis based on the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases showed that DE mRNAs and the target genes of the lncRNAs were significantly enriched in organ morphogenesis, response to stimulus, haem binding, the immune system, arginine and proline metabolism, and fatty acid biosynthesis. The prediction of mRNA–mRNA and lncRNA–mRNA interaction networks further revealed transcripts potentially involved in adaptation to high-altitude hypoxia, and the hub genes DDX24, PDCD11, EIF4A3, NDUFA11, SART1, PRPF8 and TCONS_00306477, TCONS_00306029, TCONS_00139593, TCONS_00293272, and TCONS_00313398 were selected. Additionally, a set of target genes, PIK3R1, IGF1R, FZD6, IFNB2, ATF3, MB, CYP2B4, PSMD13, and TGFB1, were also identified as candidate genes associated with high-altitude hypoxia adaptation. In conclusion, a collection of novel expressed lncRNAs, a set of target genes and biological pathways known to be relevant for altitude adaptation were identified by comparative transcriptome analysis between Tibetan sheep and Hu sheep. Our results are the first to identify the characterization and expression profile of lncRNAs between Tibetan sheep and Hu sheep and provide insights into the genetic regulation mechanisms by which Tibetan sheep adapt to high-altitude hypoxic environments.

Temporal Dynamic Transcriptome Landscape Reveals Regulatory Network During the Early Differentiation of Female Strobilus Buds in Ginkgo biloba

Reproductive bud differentiation is one of the most critical events for the reproductive success of seed plants. Yet, our understanding of genetic basis remains limited for the development of the reproductive organ of gymnosperms, namely, unisexual strobilus or cone, leaving its regulatory network largely unknown for strobilus bud differentiation. In this study, we analyzed the temporal dynamic landscapes of genes, long non-coding RNAs (lncRNAs), and microRNAs (miRNAs) during the early differentiation of female strobilus buds in Ginkgo biloba based on the whole transcriptome sequencing. Results suggested that the functions of three genes, i.e., Gb_19790 (GbFT), Gb_13989 (GinNdly), and Gb_16301 (AG), were conserved in both angiosperms and gymnosperms at the initial differentiation stage. The expression of genes, lncRNAs, and miRNAs underwent substantial changes from the initial differentiation to the enlargement of ovule stalk primordia. Besides protein-coding genes, 364 lncRNAs and 15 miRNAs were determined to be functional. Moreover, a competing endogenous RNA (ceRNA) network comprising 10,248 lncRNA-miRNA-mRNA pairs was identified, which was highly correlated with the development of ovulate stalk primordia. Using the living fossil ginkgo as the study system, this study not only reveals the expression patterns of genes related to flowering but also provides novel insights into the regulatory networks of lncRNAs and miRNAs, especially the ceRNA network, paving the way for future studies concerning the underlying regulation mechanisms of strobilus bud differentiation.

Genome-Wide DNA Methylation and Its Effect on Gene Expression During Subclinical Mastitis in Water Buffalo

Subclinical mastitis (SCM) in buffalo is one of the most challenging paradoxes for the dairy sector with very significant milk production losses and poses an imminent danger to milch animal’s milk-producing ability. We present here the genome-wide methylation specific to SCM in water buffalo and its consequential effect on the gene expression landscape for the first time. Whole-genome DNA methylation profiles from peripheral blood lymphocytes and gene expression profiles from milk somatic cells of healthy and SCM cases were catalogued from the MeDIP-Seq and RNA-Seq data. The average methylation in healthy buffaloes was found to be higher than that in the SCM-infected buffaloes. DNA methylation was abundant in the intergenic region followed by the intronic region in both healthy control and SCM groups. A total of 3,950 differentially methylated regions (DMRs) were identified and annotated to 370 differentially methylated genes (DMGs), most of which were enriched in the promoter region. Several important pathways were activated due to hypomethylation and belonged to the Staphylococcus aureus infection, Th17 cell differentiation, and antigen processing and presentation pathways along with others of defense responses. DNA methylome was compared with transcriptome to understand the regulatory role of DNA methylation on gene expression specific to SCM in buffaloes. A total of 4,778 significant differentially expressed genes (DEGs) were extracted in response to SCM, out of which 67 DMGs were also found to be differentially expressed, suggesting that during SCM, DNA methylation could be one of the epigenetic regulatory mechanisms of gene expression. Genes like CSF2RB, LOC102408349, C3 and PZP like, and CPAMD8 were found to be downregulated in our study, which are known to be involved in the immune response to SCM. Association of DNA methylation with transposable elements, miRNAs, and lncRNAs was also studied. The present study reports a buffalo SCM web resource (BSCM2TDb) available at http://webtom.cabgrid.res.in/BSCM2TDb that catalogues all the mastitis-related information of the analyses results of this study in a single place. This will be of immense use to buffalo researchers to understand the host–pathogen interaction involving SCM, which is required in endeavors of mastitis control and management.

Long Non-Coding RNAs Might Regulate Phenotypic Switch of Vascular Smooth Muscle Cells Acting as ceRNA: Implications for In-Stent Restenosis

Coronary in-stent restenosis is a late complication of angioplasty. It is a multifactorial process that involves vascular smooth muscle cells (VSMCs), endothelial cells, and inflammatory and genetic factors. In this study, the transcriptomic landscape of VSMCs’ phenotypic switch process was assessed under stimuli resembling stent injury. Co-cultured contractile VSMCs and endothelial cells were exposed to a bare metal stent and platelet-derived growth factor (PDGF-BB) 20 ng/mL. Migratory capacity (wound healing assay), proliferative capacity, and cell cycle analysis of the VSMCs were performed. RNAseq analysis of contractile vs. proliferative VSMCs was performed. Gene differential expression (DE), identification of new long non-coding RNA candidates (lncRNAs), gene ontology (GO), and pathway enrichment (KEGG) were analyzed. A competing endogenous RNA network was constructed, and significant lncRNA–miRNA–mRNA axes were selected. VSMCs exposed to “stent injury” conditions showed morphologic changes, with proliferative and migratory capacities progressing from G0-G1 cell cycle phase to S and G2-M. RNAseq analysis showed DE of 1099, 509 and 64 differentially expressed mRNAs, lncRNAs, and miRNAs, respectively. GO analysis of DE genes showed significant enrichment in collagen and extracellular matrix organization, regulation of smooth muscle cell proliferation, and collagen biosynthetic process. The main upregulated nodes in the lncRNA-mediated ceRNA network were PVT1 and HIF1-AS2, with downregulation of ACTA2-AS1 and MIR663AHG. The PVT1 ceRNA axis appears to be an attractive target for in-stent restenosis diagnosis and treatment.

Seed Germination Mechanism of Carex rigescens Under Variable Temperature Determinded Using Integrated Single-Molecule Long-Read and Illumina Sequence Analysis

The turfgrass species Carex rigescens has broad development and utilization prospects in landscaping construction. However, seed dormancy and a low germination rate have inhibited its application. Furthermore, the molecular mechanisms of seed germination in C. rigescens have not been thoroughly studied. Therefore, in the present study, PacBio full-length transcriptome sequencing combined with Illumina sequencing was employed to elucidate the germination mechanism of C. rigescens seeds under variable temperatures. In general, 156,750 full-length non-chimeric sequences, including those for 62,086 high-quality transcripts, were obtained using single-molecule long read sequencing. In total, 40,810 high-quality non-redundant, 1,675 alternative splicing, 28,393 putative coding sequences, and 1,052 long non-coding RNAs were generated. Based on the newly constructed full-length reference transcriptome, 23,147 differentially expressed genes were identified. We screened four hub genes participating in seed germination using weighted gene co-expression network analysis. Combining these results with the physiological observations, the important roles of sucrose and starch metabolic pathways in germination are further discussed. In conclusion, we report the first full-length transcriptome of C. rigescens, and investigated the physiological and transcriptional mechanisms of seed germination under variable temperatures. Our results provide valuable information for future transcriptional analyses and gene function studies of C. rigescens.

Global identification and integrated analysis of heat-responsive long non-coding RNAs in contrasting rice cultivars

Identified 2743 rice lncRNAs LncRNAs in response to heat stress Function prediction of HRLs Network among HRLs, genes and miRNAs co-localization of HRLs with QTLs Significant motifs in HRL sequences Long non-coding RNAs (lncRNAs) play vital roles in plant responses to environmental challenges. A better understanding of the gene regulation mediated by lncRNAs and their systematic identification would provide great benefits for modern agriculture. In this study, we performed strand-specific RNA sequencing for two rice varieties, heat-tolerant ZS97B and heat-susceptible SYD2 under heat stress. In total, 2743 putative lncRNAs were identified, and their expression profiles in response to heat treatments were established. We identified 231 differentially expressed lncRNAs (DELs) under heat stress, including 31 DELs common to both varieties and 103 and 97 specific to ZS97B and SYD2, respectively, all defined as heat-responsive lncRNAs (HRLs). The target-coding genes of HRLs were predicted, and GO and KEGG annotations of HRL targets revealed functions in which HRLs might be involved. The interaction network between HRLs, target genes and relevant miRNAs was constructed. The HRLs and their targets were compared with publicly available QTLs for rice seedling growth under heat stimulus. Ten HRLs and twelve target genes were linked with five heat stress-relevant QTLs. Sequence analysis revealed several motifs significantly enriched within the 231 HRL sequences. Our findings provide a valuable resource for further characterization of lncRNAs in terms of heat response and plant heat tolerance improvement.