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

LncRNA H19 accelerates renal fibrosis by negatively regulating the let-7b-5p/TGF-βR1/COL1A1 axis

BACKGROUND

Transforming growth factor-beta1 (TGF-β1)-mediated renal fibrosis is a critical pathological process of chronic kidney disease worsening to end-stage renal disease. Recent studies have shown that long noncoding RNA H19 (lncRNA H19) is widely involved in the formation and progression of fibrosis in multiple organs. However, its molecular events in renal fibrosis remain to be elucidated.

METHODS

Rats were treated with adenine intragastrically and HK-2 cells were induced by TGF-β1 to construct renal fibrosis models in vivo and in vitro, respectively. Renal histopathological examination was performed using HE and Masson staining. Gene expression levels of interleukin-1beta (IL-1β), tumor necrosis factor-alpha (TNF-α), TGF-β1, fibronectin (Fn), alpha-smooth muscle actin (α-SMA), H19, let-7b-5p, TGF-β receptor 1 (TGF-βR1), and type I collagen (COL1A1) were detected by qRT-PCR. Immunohistochemistry, immunofluorescence, and western blot analysis were used to evaluate the expression of renal fibrosis biomarkers. Dual-luciferase reporter assay was used to verify the presence of binding sites between H19 and let-7b-5p, and between let-7b-5p and TGF-βR1 and COL1A1.

RESULTS

H19 was overexpressed in both in vivo and in vitro renal fibrosis models. H19 knockdown significantly reversed TGF-β1-induced upregulation of fibronectin, COL1A1, and α-SMA and downregulation of E-cadherin in HK-2 cells, accompanied by an increase in let-7b-5p. Let-7b-5p was bound to H19 in HK-2 cells, and its overexpression inhibited TGF-β1-induced HK-2 cell fibrosis. Further experiments determined that let-7b-5p directly targets TGF-βR1 and COL1A1 in HK-2 cells. In addition, inhibition of let-7b-5p reversed the reduction in HK-2 cell fibrosis induced by H19 knockdown. Finally, knockdown of H19 alleviated renal fibrosis in vivo and was associated with regulation of the let-7b-5p/TGF-βR1/COL1A1 axis.

CONCLUSION

Our results indicate that knockdown of H19 inhibits renal tubular epithelial fibrosis by negatively regulating the let-7b-5p/TGF-βR1/COL1A1 axis, which may provide new mechanistic insights into CRF progression.

Plant lncRNA-miRNA Interaction Prediction Based on Counterfactual Heterogeneous Graph Attention Network

Identifying interactions between long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) provides a new perspective for understanding regulatory relationships in plant life processes. Recently, computational methods based on graph neural networks (GNNs) have been widely employed to predict lncRNA-miRNA interactions (LMIs), which compensate for the inadequacy of biological experiments. However, the low-semantic and noise of graph limit the performance of existing GNN-based methods. In this paper, we develop a novel Counterfactual Heterogeneous Graph Attention Network (CFHAN) to improve the robustness to against the noise and the prediction of plant LMIs. Firstly, we construct a real-world based lncRNA-miRNA (L-M) heterogeneous network. Secondly, CFHAN utilizes the node-level attention, the semantic-level attention, and the counterfactual links to enhance the node embeddings learning. Finally, these embeddings are used as inputs for Multilayer Perceptron (MLP) to predict the interactions between lncRNAs and miRNAs. Evaluating our method on a benchmark dataset of plant LMIs, CFHAN outperforms five state-of-the-art methods, and achieves an average AUC and average ACC of 0.9953 and 0.9733, respectively. This demonstrates CFHAN's ability to predict plant LMIs and exhibits promising cross-species prediction ability, offering valuable insights for experimental LMI researches.

ETENLNC: An end to end lncRNA identification and analysis framework to facilitate construction of known and novel lncRNA regulatory networks

Long non-coding RNAs (lncRNAs) play crucial roles in the regulation of gene expression and maintenance of genomic integrity through various interactions with DNA, RNA, and proteins. The availability of large-scale sequence data from various high-throughput platforms has opened possibilities to identify, predict, and functionally annotate lncRNAs. As a result, there is a growing demand for an integrative computational framework capable of identifying known lncRNAs, predicting novel lncRNAs, and inferring the downstream regulatory interactions of lncRNAs at the genome-scale. We present ETENLNC (End-To-End-Novel-Long-NonCoding), a user-friendly, integrative, open-source, scalable, and modular computational framework for identifying and analyzing lncRNAs from raw RNA-Seq data. ETENLNC employs six stringent filtration steps to identify novel lncRNAs, performs differential expression analysis of mRNA and lncRNA transcripts, and predicts regulatory interactions between lncRNAs, mRNAs, miRNAs, and proteins. We benchmarked ETENLNC against six existing tools and optimized it for desktop workstations and high-performance computing environments using data from three different species. ETENLNC is freely available on GitHub: https://github.com/EvolOMICS-TU/ETENLNC.

Differential abundance of microRNAs in seminal plasma extracellular vesicles (EVs) in Sahiwal cattle bull related to male fertility

Sahiwal cattle, known for their high milk yield, are propagated through artificial insemination (AI) using male germplasm, largely contingent on semen quality. Spermatozoa, produced in the testes, carry genetic information and molecular signals essential for successful fertilization. Seminal plasma, in addition to sperm, contains nano-sized lipid-bound extracellular vesicles (SP-EVs) that carry key biomolecules, including fertility-related miRNAs, which are essential for bull fertility. The current study focused on miRNA profiling of SP-EVs from high-fertile (HF) and low-fertile (LF) Sahiwal bulls. SP-EVs were isolated using size exclusion chromatography (SEC) and characterized by dynamic light scattering (DLS) and nanoparticle tracking analysis (NTA). Western blotting detected the EV-specific protein markers TSG101 and CD63. The DLS analysis showed SP-EV sizes of 170-180 nm in HF and 130-140 nm in LF samples. The NTA revealed particle concentrations of 5.76 × 1010 to 5.86 × 1011 particles/mL in HF and 5.31 × 1010 to 2.70 × 1011 particles/mL in LF groups, with no significant differences in size and concentration between HF and LF. High-throughput miRNA sequencing identified 310 miRNAs in SP-EVs from both groups, with 61 upregulated and 119 downregulated in HF bull. Further analysis identified 41 miRNAs with significant fold changes and p-values, including bta-miR-1246, bta-miR-195, bta-miR-339b, and bta-miR-199b, which were analyzed for target gene prediction. Gene Ontology (GO) and KEGG pathway analyses indicated that these miRNAs target genes involved in transcription regulation, ubiquitin-dependent endoplasmic reticulum-associated degradation (ERAD) pathways, and signalling pathways. Functional exploration revealed that these genes play roles in spermatogenesis, motility, acrosome reactions, and inflammatory responses. qPCR analysis showed that bta-miR-195 had 80% higher expression in HF spermatozoa compared to LF, suggesting its association with fertility status (p < 0.05). In conclusion, this study elucidates the miRNA cargoes in SP-EVs as indicators of Sahiwal bull fertility, highlighting bta-miR-195 as a potential fertility factor among the various miRNAs identified.

Modelling cell type-specific lncRNA regulatory network in autism with Cycle

BACKGROUND

Autism spectrum disorder (ASD) is a class of complex neurodevelopment disorders with high genetic heterogeneity. Long non-coding RNAs (lncRNAs) are vital regulators that perform specific functions within diverse cell types and play pivotal roles in neurological diseases including ASD. Therefore, exploring lncRNA regulation would contribute to deciphering ASD molecular mechanisms. Existing computational methods utilize bulk transcriptomics data to identify lncRNA regulation in all of samples, which could reveal the commonalities of lncRNA regulation in ASD, but ignore the specificity of lncRNA regulation across various cell types.

RESULTS

Here, we present Cycle (Cell type-specific lncRNA regulatory network) to construct the landscape of cell type-specific lncRNA regulation in ASD. We have found that each ASD cell type is unique in lncRNA regulation, and more than one-third and all cell type-specific lncRNA regulatory networks are characterized as scale-free and small-world, respectively. Across 17 ASD cell types, we have discovered 19 rewired and 11 stable modules, along with eight rewired and three stable hubs within the constructed cell type-specific lncRNA regulatory networks. Enrichment analysis reveals that the discovered rewired and stable modules and hubs are closely related to ASD. Furthermore, more similar ASD cell types tend to be connected with higher strength in the constructed cell similarity network. Finally, the comparison results demonstrate that Cycle is a potential method for uncovering cell type-specific lncRNA regulation.

CONCLUSION

Overall, these results illustrate that Cycle is a promising method to model the landscape of cell type-specific lncRNA regulation, and provides insights into understanding the heterogeneity of lncRNA regulation between various ASD cell types.

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

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

Comparative Tissue Identification and Characterization of Long Non-Coding RNAs in the Globally Distributed Blue Shark Prionace glauca

Long non-coding RNAs (lncRNAs) are involved in numerous biological processes and serve crucial regulatory functions in both animals and plants. Nevertheless, there is limited understanding of lncRNAs and their patterns of expression and roles in sharks. In the current study, we systematically identified and characterized lncRNAs in the blue shark (Prionace glauca) from four tissues (liver, spleen, muscle, and kidney) using high-throughput sequencing and bioinformatics tools. A total of 21,932 high-confidence lncRNAs were identified, with 8984 and 3067 stably and tissue-specific expressed lncRNAs, respectively. In addition, a total of 45,007 differentially expressed (DE) lncRNAs were obtained among tissues, with kidney versus muscle having the largest numbers across tissues. DE lncRNAs trans target protein-coding genes were predicted, and functional gene ontology enrichment of these genes showed GO terms such as muscle system processes, cellular/metabolic processes, and stress and immune responses, all of which correspond with the specific biological functions of each tissue analyzed. These results advance our knowledge of lncRNAs in sharks and present novel data on tissue-specific lncRNAs, providing key information to support future functional shark investigations.

Genome-wide characterization of drought-responsive long non-coding RNAs in sorghum (Sorghum bicolor)

Drought stress strongly affects crop yield. Although knowledge of long non-coding RNAs (lncRNAs) has been updated continuously and rapidly, information about lncRNAs in drought resistance regulation is extremely limited in sorghum. Here, lncRNA-sequencing was performed with seedlings of a sorghum cultivar (Jinza29) under three water control treatments to investigate the mechanism of lncRNAs responsible for drought resistance in sorghum. A total of 377 differentially expressed lncRNAs (DElncRNAs) were identified. We also predicted 4322 and 2827 transcripts as potential cis-target and trans-target genes for drought-responsive lncRNAs, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that those target genes exhibited marked enrichment into "oxidoreductase activity", "signal transducer activity", "DNA repair", "photosynthesis", "glutathione metabolism", and "phenylpropanoid biosynthesis" and other terms associated with abiotic stress resistance. Moreover, several lncRNAs were estimated to modulate the expression of other genes related to stress response and photosynthetic carbon metabolism. Additionally, we found 107 DElncRNAs that might be candidate target mimics for 56 miRNAs. LncRNAs play important roles in drought adaptation of sorghum through interacting with protein-encoding genes. The obtained results provided novel insights into the biological characteristics of lncRNAs and offered potential regulatory factors for genetically enhancing drought resistance in sorghum.

Source leaves are regulated by sink strengths through non-coding RNAs and alternative polyadenylation in cucumber (Cucumis sativus L.)

BACKGROUND

The yield of major crops is generally limited by sink capacity and source strength. Cucumber is a typical raffinose family oligosaccharides (RFOs)-transporting crop. Non-coding RNAs and alternative polyadenylation (APA) play important roles in the regulation of growth process in plants. However, their roles on the sink‒source regulation have not been demonstrated in RFOs-translocating species.

RESULTS

Here, whole-transcriptome sequencing was applied to compare the leaves of cucumber under different sink strength, that is, no fruit-carrying leaves (NFNLs) and fruit-carrying leaves (FNLs) at 12th node from the bottom. The results show that 1101 differentially expressed (DE) mRNAs, 79 DE long non-coding RNAs (lncRNAs) and 23 DE miRNAs were identified, which were enriched in photosynthesis, energy production and conversion, plant hormone signal transduction, starch and carbohydrate metabolism and protein synthesis pathways. Potential co-expression networks like, DE lncRNAs-DE mRNAs/ DE miRNAs-DE mRNAs, and competing endogenous RNA (ceRNA) regulation models (DE lncRNAs-DE miRNAs-DE mRNAs) associated with sink‒source allocation, were constructed. Furthermore, 37 and 48 DE genes, which enriched in MAPK signaling and plant hormone signal transduction pathway, exist differentially APA, and SPS (CsaV3_2G033300), GBSS1 (CsaV3_5G001560), ERS1 (CsaV3_7G029600), PNO1 (CsaV3_3G003950) and Myb (CsaV3_3G022290) may be regulated by both ncRNAs and APA between FNLs and NFNLs, speculating that ncRNAs and APA are involved in the regulation of gene expression of cucumber sink‒source carbon partitioning.

CONCLUSIONS

These results reveal a comprehensive network among mRNAs, ncRNAs, and APA in cucumber sink-source relationships. Our findings also provide valuable information for further research on the molecular mechanism of ncRNA and APA to enhance cucumber yield.

Identification and Characterization of miRNAs and lncRNAs Associated with Salinity Stress in Rice Panicles

Salinity is a common abiotic stress that limits crop productivity. Although there is a wealth of evidence suggesting that miRNA and lncRNA play important roles in the response to salinity in rice seedlings and reproductive stages, the mechanism by which competing endogenous RNAs (ceRNAs) influence salt tolerance and yield in rice has been rarely reported. In this study, we conducted full whole-transcriptome sequencing of rice panicles during the reproductive period to clarify the role of ceRNAs in the salt stress response and yield. A total of 214 lncRNAs, 79 miRNAs, and 584 mRNAs were identified as differentially expressed RNAs under salt stress. Functional analysis indicates that they play important roles in GO terms such as response to stress, biosynthesis processes, abiotic stimuli, endogenous stimulus, and response to stimulus, as well as in KEGG pathways such as secondary metabolite biosynthesis, carotenoid biosynthesis, metabolic pathways, and phenylpropanoid biosynthesis. A ceRNA network comprising 95 lncRNA-miRNA-mRNA triplets was constructed. Two lncRNAs, MSTRG.51634.2 and MSTRG.48576.1, were predicted to bind to osa-miR172d-5p to regulate the expression of OsMYB2 and OsMADS63, which have been reported to affect salt tolerance and yield, respectively. Three lncRNAs, MSTRG.30876.1, MSTRG.44567.1, and MSTRG.49308.1, may bind to osa-miR5487 to further regulate the expression of a stress protein (LOC_Os07g48460) and an aquaporin protein (LOC_Os02g51110) to regulate the salt stress response. This study is helpful for understanding the underlying molecular mechanisms of ceRNA that drive the response of rice to salt stress and provide new genetic resources for salt-resistant rice breeding.

Transcriptomic analyses of differentially expressed human genes, micro RNAs and long-non-coding RNAs in severe, symptomatic and asymptomatic malaria infection

Malaria transmission and endemicity in Africa remains hugely disproportionate compared to the rest of the world. The complex life cycle of P. falciparum (Pf) between the vertebrate human host and the anopheline vector results in differential expression of genes within and between hosts. An in-depth understanding of Pf interaction with various human genes through regulatory elements will pave way for identification of newer tools in the arsenal for malaria control. Therefore, the regulatory elements (REs) involved in the over- or under-expression of various host immune genes hold the key to elucidating alternative control measures that can be applied for disease surveillance, prompt diagnosis and treatment. We carried out an RNAseq analysis to identify differentially expressed genes and network elucidation of non-coding RNAs and target genes associated with immune response in individuals with different clinical outcomes. Raw RNAseq datasets, retrieved for analyses include individuals with severe (Gambia-20), symptomatic (Burkina Faso-15), asymptomatic (Mali-16) malaria as well as uninfected controls (Tanzania-20; Mali-36). Of the total 107 datasets retrieved, we identified 5534 differentially expressed genes (DEGs) among disease and control groups. A peculiar pattern of DEGs was observed, with individuals presenting with severe/symptomatic malaria having the highest and most diverse upregulated genes, while a reverse phenomenon was recorded among asymptomatic and uninfected individuals. In addition, we identified 141 differentially expressed micro RNA (miRNA), of which 78 and 63 were upregulated and downregulated respectively. Interactome analysis revealed a moderate interaction between DEGs and miRNAs. Of all identified miRNA, five were unique (hsa-mir-32, hsa-mir-25, hsa-mir-221, hsa-mir-29 and hsa-mir-148) because of their connectivity to several genes, including hsa-mir-221 connected to 16 genes. Six-hundred and eight differentially expressed long non coding RNA (lncRNA) were also identified, including SLC7A11, LINC01524 among the upregulated ones. Our study provides important insight into host immune genes undergoing differential expression under different malaria conditions. It also identified unique miRNAs and lncRNAs that modify and/or regulate the expression of various immune genes. These regulatory elements we surmise, have the potential to serve a diagnostic purpose in discriminating between individuals with severe/symptomatic malaria and those with asymptomatic infection or uninfected, following further clinical validation from field isolates.

Leishmania infantum infection modulates messenger RNA, microRNA and long non-coding RNA expression in human neutrophils in vitro

In the Americas, L. infantum (syn. chagasi) is the main cause of human visceral leishmaniasis. The role of neutrophils as part of the innate response to Leishmania spp. infection is dubious and varies according to the species causing the infection. Global expression of coding RNAs, microRNAs and long non-coding RNAs changes as part of the immune response against pathogens. Changes in mRNA and non-coding RNA expression resulting from infection by Leishmania spp. are widely studied in macrophages, but scarce in neutrophils, the first cell to encounter the trypanosomatid, especially following infection by L. infantum. Herein, we aimed to understand the expression patterns of coding and non-coding transcripts during acute in vitro infection of human neutrophils by L. infantum. We isolated neutrophils from whole blood of healthy male donors (n = 5) and split into groups: 1) infected with L. infantum (MOI = 5:1), and 2) uninfected controls. After 3 hours of exposure of infected group to promastigotes of L. infantum, followed by 17 hours of incubation, total RNA was extracted and total RNA-Seq and miRNA microarray were performed. A total of 212 genes were differentially expressed in neutrophils following RNA-Seq analysis (log2(FC)±0.58, FDR≤0.05). In vitro infection with L. infantum upregulated the expression of 197 and reduced the expression of 92 miRNAs in human neutrophils (FC±2, FDR≤0.01). Lastly, 5 downregulated genes were classified as lncRNA, and of the 10 upregulated genes, there was only 1 lncRNA. Further bioinformatic analysis indicated that changes in the transcriptome and microtranscriptome of neutrophils, following in vitro infection with L. infantum, may impair phagocytosis, apoptosis and decrease nitric oxide production. Our work sheds light on several mechanisms used by L. infantum to control neutrophil-mediated immune response and identifies several targets for future functional studies, aiming at the development of preventive or curative treatments for this prevalent zoonosis.

Extracellular vesicle-packaged lncRNA from cancer-associated fibroblasts promotes immune evasion by downregulating HLA-A in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is characterised by immune evasion that contribute to poor prognosis. Cancer-associated fibroblasts (CAFs) play a pivotal role in orchestrating the PDAC tumour microenvironment. We investigated the role of CAF-derived extracellular vesicle (EV)-packaged long non-coding RNAs (lncRNAs) in immune evasion and explored gene therapy using engineered EVs loading small interfering RNAs (siRNAs) as a potential therapeutic strategy. Our findings highlight the significance of EV-packaged lncRNA RP11-161H23.5 from CAF in promoting PDAC immune evasion by downregulating HLA-A expression, a key component of antigen presentation. Mechanistically, RP11-161H23.5 forms a complex with CNOT4, a subunit of the mRNA deadenylase CCR4-NOT complex, enhancing the degradation of HLA-A mRNA by shortening its poly(A) tail. This immune evasion mechanism compromises the anti-tumour immune response. To combat this, we propose an innovative approach utilising engineered EVs as natural and biocompatible nanocarriers for siRNA-based gene therapy and this strategy holds promise for enhancing the effectiveness of immunotherapy in PDAC. Overall, our study sheds light on the critical role of CAF-derived EV-packaged lncRNA RP11-161H23.5/CNOT4/HLA-A axis in PDAC immune evasion and presents a novel avenue for therapeutic intervention.

Full-length transcriptome sequencing of Arabidopsis plants provided new insights into the autophagic regulation of photosynthesis

Autophagy is a highly conserved eukaryotic pathway and plays a crucial role in cell survival under stress conditions. Here, we applied a full-length transcriptome approach to study an Arabidopsis autophagy mutant (atg5-1) subjected to nitrogen-starvation, using Oxford Nanopore Technologies. A total of 39,033 transcripts were identified, including 11,356 new transcripts. In addition, alternative splicing (AS) events and lncRNAs were also detected between Col-0 (WT) and atg5-1. Differentially expressed transcript enrichment showed that autophagy upregulates the expression of many stress-responsive genes and inhibits the transcription of photosynthesis-associated genes. The qRT-PCR results showed that the expression patterns of photosynthesis-related genes in the atg5-1 differed under the conditions of nitrogen starvation and carbon starvation. Under nitrogen starvation treatment, many genes related to photosynthesis also exhibited AS. Chlorophyll fluorescence images revealed that the Fv/Fm and ΦPSII of old atg5-1 leaves were significantly reduced after nitrogen starvation treatment, but the Y(NPQ) indices were significantly increased compared to those of the WT plants. The results of qRT-PCR suggest that autophagy appears to be involved in the degradation of genes related to photodamage repair in PSII. Taken together, the full-length transcriptiome sequencing provide new insights into how new transcripts, lncRNAs and alternative splicing (AS) are involved in plant autophagy through full-length transcriptome sequencing and suggest a new potential link between autophagy and photosynthesis.

Sex-biased gene expression precedes sexual dimorphism in the agonadal annelid Platynereis dumerilii

Gametogenesis is the process by which germ cells differentiate into mature sperm and oocytes, cells essential for sexual reproduction. The sex-specific molecular programs that drive spermatogenesis and oogenesis can also serve as sex identification markers. Platynereis dumerilii is a research organism that has been studied in many areas of developmental biology. However investigations often disregard sex, as P. dumerilii juveniles lack sexual dimorphism. The molecular mechanisms of gametogenesis in the segmented worm P. dumerilii are also largely unknown. In this study, we used RNA sequencing to investigate the transcriptomic profiles of gametogenesis in P. dumerilii juveniles. Our analysis revealed that sex-biased gene expression becomes increasingly pronounced during the advanced developmental stages, particularly during the meiotic phases of gametogenesis. We identified conserved genes associated with spermatogenesis, such as dmrt1, and a novel gene psmt, that is associated with oogenesis. Additionally, putative long non-coding RNAs were upregulated in both male and female gametogenic programs. This study provides a foundational resource for germ cell research in P. dumerilii, markers for sex identification, and offers comparative data to enhance our understanding of the evolution of gametogenesis mechanisms across species.

Characterisation of sRNAs enriched in outer membrane vesicles of pathogenic Flavobacterium psychrophilum causing Bacterial Cold Water Disease in rainbow trout

Flavobacterium psychrophilum (Fp) causes Bacterial Cold Water Disease in salmonids. During host-pathogen interactions, gram-negative bacteria, such as Fp, release external membrane vesicles (OMVs) harbouring cargos, such as DNA, RNA and virulence factors. This study aimed to characterise the potential role of the OMVs' small RNAs (sRNAs) in the Fp-rainbow trout host-pathogen interactions. sRNAs carried within OMVs were isolated from Fp. RNA-Seq datasets from whole-cell Fp and their isolated OMVs indicated substantial enrichment of specific sRNAs in the OMVs compared to the parent cell. Many of the OMV-packaged sRNAs were located in the pathogenicity islands of Fp. Conservation of sRNAs in 65 strains with variable degrees of virulence was reported. Dual RNA-Seq of host and pathogen transcriptomes on day 5 post-infection of Fp -resistant and -susceptible rainbow trout genetic lines revealed correlated expression of OMV-packaged sRNAs and their predicted host's immune gene targets. In vitro, treatment of the rainbow trout epithelial cell line RTgill-W1 with OMVs showed signs of cytotoxicity accompanied by dynamic changes in the expression of host genes when profiled 24 h following treatment. The OMV-treated cells, similar to the Fp -resistant fish, showed downregulated expression of the suppressor of cytokine signalling 1 (SOCS1) gene, suggesting induction of phagosomal maturation. Other signs of modulating the host gene expression following OMV-treatment include favouring elements from the phagocytic, endocytic and antigen presentation pathways in addition to HSP70, HSP90 and cochaperone proteins, which provide evidence for a potential role of OMVs in boosting the host immune response. In conclusion, the study identified novel microbial targets and inherent characteristics of OMVs that could open up new avenues of treatment and prevention of fish infections.

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

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

Update on functional analysis of long non-coding RNAs in common crops

With the rapid advances in next-generation sequencing technology, numerous non-protein-coding transcripts have been identified, including long noncoding RNAs (lncRNAs), which are functional RNAs comprising more than 200 nucleotides. Although lncRNA-mediated regulatory processes have been extensively investigated in animals, there has been considerably less research on plant lncRNAs. Nevertheless, multiple studies on major crops showed lncRNAs are involved in crucial processes, including growth and development, reproduction, and stress responses. This review summarizes the progress in the research on lncRNA roles in several major crops, presents key strategies for exploring lncRNAs in crops, and discusses current challenges and future prospects. The insights provided in this review will enhance our comprehension of lncRNA functions in crops, with potential implications for improving crop genetics and breeding.

Regulatory roles of long non-coding RNAs in short-term heat stress in adult worker bees

Long non-coding RNAs (lncRNAs) are crucial modulators of post-transcriptional gene expression regulation, cell fate determination, and disease development. However, lncRNA functions during short-term heat stress in adult worker bees are poorly understood. Here, we performed deep sequencing and bioinformatic analyses of honeybee lncRNAs. RNA interference was performed by using siRNA targeting the most highly expressed lncRNA. The silencing effect on lncRNA and the relative expression levels of seven heat shock protein (HSP) genes, were subsequently examined. Overall, 7,842 lncRNAs and 115 differentially expressed lncRNAs (DELs) were identified in adult worker bees following heat stress exposure. Structural analysis revealed that the overall expression abundance, length of transcripts, exon number, and open reading frames of lncRNAs were lower than those of mRNAs. GO analysis revealed that the target genes were mainly involved in "metabolism," "protein folding," "response to stress," and "signal transduction" pathways. KEGG analysis indicated that the "protein processing in endoplasmic reticulum" and "longevity regulating pathway-multiple species" pathways were most enriched. Quantitative real-time polymerase chain reaction (qRT-PCR) detection of the selected DELs confirmed the reliability of the sequencing data. Moreover, the siRNA experiment indicated that feeding siRNA yielded a silencing efficiency of 77.51% for lncRNA MSTRG.9645.5. Upon silencing this lncRNA, the expression levels of three HSP genes were significantly downregulated (p < 0.05), whereas those of three other HSP genes were significantly upregulated (p < 0.05). Our results provide a new perspective for understanding the regulatory mechanisms of lncRNAs in adult worker bees under short-term heat stress.

Transcriptome Analysis of Sesame (Sesamum indicum L.) Reveals the LncRNA and mRNA Regulatory Network Responding to Low Nitrogen Stress

Nitrogen is one of the important factors restricting the development of sesame planting and industry in China. Cultivating sesame varieties tolerant to low nitrogen is an effective way to solve the problem of crop nitrogen deficiency. To date, the mechanism of low nitrogen tolerance in sesame has not been elucidated at the transcriptional level. In this study, two sesame varieties Zhengzhi HL05 (ZZ, nitrogen efficient) and Burmese prolific (MD, nitrogen inefficient) in low nitrogen were used for RNA-sequencing. A total of 3964 DEGs (differentially expressed genes) and 221 DELs (differentially expressed lncRNAs) were identified in two sesame varieties at 3d and 9d after low nitrogen stress. Among them, 1227 genes related to low nitrogen tolerance are mainly located in amino acid metabolism, starch and sucrose metabolism and secondary metabolism, and participate in the process of transporter activity and antioxidant activity. In addition, a total of 209 pairs of lncRNA-mRNA were detected, including 21 pairs of trans and 188 cis. WGCNA (weighted gene co-expression network analysis) analysis divided the obtained genes into 29 modules; phenotypic association analysis identified three low-nitrogen response modules; through lncRNA-mRNA co-expression network, a number of hub genes and cis/trans-regulatory factors were identified in response to low-nitrogen stress including GS1-2 (glutamine synthetase 1-2), PAL (phenylalanine ammonia-lyase), CHS (chalcone synthase, CHS), CAB21 (chlorophyll a-b binding protein 21) and transcription factors MYB54, MYB88 and NAC75 and so on. As a trans regulator, lncRNA MSTRG.13854.1 affects the expression of some genes related to low nitrogen response by regulating the expression of MYB54, thus responding to low nitrogen stress. Our research is the first to provide a more comprehensive understanding of DEGs involved in the low nitrogen stress of sesame at the transcriptome level. These results may reveal insights into the molecular mechanisms of low nitrogen tolerance in sesame and provide diverse genetic resources involved in low nitrogen tolerance research.

Whole transcriptome mapping reveals the lncRNA regulatory network of TFP5 treatment in diabetic nephropathy

BACKGROUND

TFP5 is a Cdk5 inhibitor peptide, which could restore insulin production. However, the role of TFP5 in diabetic nephropathy (DN) is still unclear.

OBJECTIVE

This study aims to characterize the transcriptome profiles of mRNA and lncRNA in TFP5-treated DN mice to mine key lncRNAs associated with TFP5 efficacy.

METHODS

We evaluated the role of TFP5 in DN pathology and performed RNA sequencing in C57BL/6J control mice, C57BL/6J db/db model mice, and TFP5 treatment C57BL/6J db/db model mice. The differentially expressed lncRNAs (DElncRNAs) and mRNAs (DEmRNAs) were analyzed. WGCNA was used to screen hub-gene of TFP5 in treatment of DN.

RESULTS

Our results showed that TFP5 therapy ameliorated renal tubular injury in DN mice. In addition, compared with the control group, the expression profile of lncRNAs in the model group was significantly disordered, while TFP5 alleviated the abnormal expression of lncRNAs. A total of 67 DElncRNAs shared among the three groups, 39 DElncRNAs showed a trend of increasing in the DN group and decreasing after TFP treatment, while the remaining 28 showed the opposite trend. DElncRNAs were enriched in glycosphingolipid biosynthesis signaling pathways, NF-κB signaling pathways, and complement activation signaling pathways. There were 1028 up-regulated and 1117 down-regulated DEmRNAs in the model group compared to control group, and 123 up-regulated and 153 down-regulated DEmRNAs in the TFP5 group compared to the model group. The DEmRNAs were involved in PPAR and MAPK signaling pathway. We confirmed that MSTRG.28304.1 is a key DElncRNA for TFP5 treatment of DN. TFP5 ameliorated DN maybe by inhibiting MSTRG.28304.1 through regulating the insulin resistance and PPAR signaling pathway. The qRT-PCR results confirmed the reliability of the sequencing data through verifying the expression of ENSMUST00000211209, MSTRG.31814.5, MSTRG.28304.1, and MSTRG.45642.14.

CONCLUSION

Overall, the present study provides novel insights into molecular mechanisms of TFP5 treatment in DN.

Exploring the regulatory role of non-coding RNAs in fiber development and direct regulation of GhKCR2 in the fatty acid metabolic pathway in upland cotton

Cotton fiber holds immense importance as the primary raw material for the textile industry. Consequently, comprehending the regulatory mechanisms governing fiber development is pivotal for enhancing fiber quality. Our study aimed to construct a regulatory network of competing endogenous RNAs (ceRNAs) and assess the impact of non-coding RNAs on gene expression throughout fiber development. Through whole transcriptome data analysis, we identified differentially expressed genes (DEGs) regulated by non-coding RNA (ncRNA) that were predominantly enriched in phenylpropanoid biosynthesis and the fatty acid elongation pathway. This analysis involved two contrasting phenotypic materials (J02-508 and ZRI015) at five stages of fiber development. Additionally, we conducted a detailed analysis of genes involved in fatty acid elongation, including KCS, KCR, HACD, ECR, and ACOT, to unveil the factors contributing to the variation in fatty acid elongation between J02-508 and ZRI015. Through the integration of histochemical GUS staining, dual luciferase assay experiments, and correlation analysis of expression levels during fiber development stages for lncRNA MSTRG.44818.23 (MST23) and GhKCR2, we elucidated that MST23 positively regulates GhKCR2 expression in the fatty acid elongation pathway. This identification provides valuable insights into the molecular mechanisms underlying fiber development, emphasizing the intricate interplay between non-coding RNAs and protein-coding genes.

Long Non-Coding RNAs as Determinants of Thyroid Cancer Phenotypes: Investigating Differential Gene Expression Patterns and Novel Biomarker Discovery

Thyroid Cancer (TC) is the most common endocrine malignancy, with increasing incidence globally. Papillary thyroid cancer (PTC), a differentiated form of TC, accounts for approximately 90% of TC and occurs predominantly in women of childbearing age. Although responsive to current treatments, recurrence of PTC by middle age is common and is much more refractive to treatment. Undifferentiated TC, particularly anaplastic thyroid cancer (ATC), is the most aggressive TC subtype, characterized by it being resistant and unresponsive to all therapeutic and surgical interventions. Further, ATC is one of the most aggressive and lethal malignancies across all cancer types. Despite the differences in therapeutic needs in differentiated vs. undifferentiated TC subtypes, there is a critical unmet need for the identification of molecular biomarkers that can aid in early diagnosis, prognosis, and actionable therapeutic targets for intervention. Advances in the field of cancer genomics have enabled for the elucidation of differential gene expression patterns between tumors and healthy tissue. A novel category of molecules, known as non-coding RNAs, can themselves be differentially expressed, and extensively contribute to the up- and downregulation of protein coding genes, serving as master orchestrators of regulated and dysregulated gene expression patterns. These non-coding RNAs have been identified for their roles in driving carcinogenic patterns at various stages of tumor development and have become attractive targets for study. The identification of specific genes that are differentially expressed can give insight into mechanisms that drive carcinogenic patterns, filling the gaps of deciphering molecular and cellular processes that modulate TC subtypes, outside of well-known driver mutations.

Genome-wide characterization of post-transcriptional processes related to wood formation in Dalbergia odorifera

BACKGROUND

Alternative polyadenylation (APA), alternative splicing (AS), and long non-coding RNAs (lncRNAs) play regulatory roles in post-transcriptional processes in plants. However, little is known about their involvement in xylem development in Dalbergia odorifera, a valuable rosewood species with medicinal and commercial significance. We addressed this by conducting Isoform Sequencing (Iso-Seq) using PacBio's SMRT technology and combined it with RNA-seq analysis (RNA sequencing on Illumina platform) after collecting xylem samples from the transition zone and the sapwood of D. odorifera.

RESULTS

We identified 14,938 full-length transcripts, including 9,830 novel isoforms, which has updated the D. odorifera genome annotation. Our analysis has revealed that 4,164 genes undergo APA, whereas 3,084 genes encounter AS. We have also annotated 118 lncRNAs. Furthermore, RNA-seq analysis identified 170 differential alternative splicing (DAS) events, 344 genes with differential APA site usage (DE-APA), and 6 differentially expressed lncRNAs in the transition zone when compared to the sapwood. AS, APA, and lncRNAs are differentially regulated during xylem development. Differentially expressed APA genes were enriched for terpenoid and flavonoid metabolism, indicating their role in the heartwood formation. Additionally, DE-APA genes were associated with cell wall biosynthesis and terpenoid metabolism, implying an APA's role in wood formation. A DAS gene (involved in chalcone accumulation) with a significantly greater inclusion of the last exon in the transition zone than in the sapwood was identified. We also found that differentially expressed lncRNAs targeted the genes related to terpene synthesis.

CONCLUSIONS

This study enhances our understanding of the molecular regulatory mechanisms underlying wood formation in D. odorifera, and provides valuable genetic resources and insights for its molecular-assisted breeding.

High-quality faba bean reference transcripts generated using PacBio and Illumina RNA-seq data

The genome of faba bean was first published in 2023. To promote future molecular breeding studies, we improved the quality of the faba genome based on high-density genetic maps and the Illumina and Pacbio RNA-seq datasets. Two high-density genetic maps were used to conduct the scaffold ordering and orientation of faba bean, culminating in an increased length (i.e., 14.28 Mbp) of chromosomes and a decrease in the number of scaffolds by 45. In gene model mining and optimisation, the PacBio and Illumina RNA-seq datasets from 37 samples allowed for the identification and correction 121,606 transcripts, and the data facilitated a prediction of 15,640 alternative splicing events, 2,148 lncRNAs, and 1,752 fusion transcripts, thus allowing for a clearer understanding of the gene structures underlying the faba genome. Moreover, a total of 38,850 new genes including 56,188 transcripts were identified compared with the reference genome. Finally, the genetic data of the reference genome was integrated and a comprehensive and complete faba bean transcriptome sequence of 103,267 transcripts derived from 54,753 uni-genes was formed.

Bioinformatics analysis of a disease-specific lncRNA-miRNA-mRNA regulatory network in recurrent spontaneous abortion (RSA)

BACKGROUND

This study investigated the molecular mechanisms of long non-coding RNAs (lncRNAs) in RSA using the lncRNA-miRNA-mRNA regulatory network.

METHODS

The present study obtained expression datasets of long non-coding RNAs (lncRNAs), messenger RNAs (mRNAs), and microRNAs (miRNAs) from blood samples of individuals with unexplained recurrent spontaneous abortion (RSA) and healthy controls. Differentially expressed lncRNAs (DELs), mRNAs (DEMs), and miRNAs (DEmiRs) were identified. A regulatory network comprising lncRNA, miRNA, and mRNA was constructed, and Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to analyze the biological functions of DEM. Also, a protein-protein interaction (PPI) network was made and key genes were identified.

RESULTS

A total of 57 DELs, 212 DEmiRs, and 301 DEMs regarding RSA were identified. Later analysis revealed a lncRNA-miRNA-mRNA network comprising nine lncRNAs, 14 miRNAs, and 65 mRNAs. Then, the ceRNA network genes were subjected to functional enrichment and pathway analysis, which showed their association with various processes, such as cortisol and thyroid hormone synthesis and secretion, human cytomegalovirus infection, and parathyroid hormone synthesis. In addition, ten hub genes (ITGB3, GNAI2, GNAS, SRC, PLEC, CDC42, RHOA, RAC1, CTNND1, and FN1) were identified based on the PPI network results.

CONCLUSION

In summary, the outcomes of our study provided some data regarding the alteration genes involved in RSA pathogenic mechanism via the lncRNA-miRNA-mRNA network and reveal the possibility of identifying new lncRNAs and miRNAs as promising molecular biomarkers.

Screening and identification of lncRNAs in preadipocyte differentiation in sheep

Studies of preadipocyte differentiation and fat deposition in sheep have mainly focused on functional genes, and with no emphasis placed on the role that long non-coding RNAs (lncRNAs) may have on the activity of those genes. Here, the expression profile of lncRNAs in ovine preadipocyte differentiation was investigated and the differentially expressed lncRNAs were screened on day 0 (D0), day 2(D2) and day 8(D8) of ovine preadipocyte differentiation, with their target genes being predicted. The competing endogenous RNA (ceRNA) regulatory network was constructed by GO and KEGG enrichment analysis for functional annotation, and some differentially expressed lncRNAs were randomly selected to verify the RNA-Seq results by RT-qPCR. In the study, a total of 2517 novel lncRNAs and 3943 known lncRNAs were identified from ovine preadipocytes at the three stages of differentiation, with the highest proportion being intergenic lncRNAs. A total of 3455 lncRNAs were expressed at all three stages of preadipocyte differentiation, while 214, 226 and 228 lncRNAs were uniquely expressed at day 0, day 2 and day 8, respectively. By comparing the expression of the lncRNAs between the three stages of differentiation stages, a total of 405, 272 and 359 differentially expressed lncRNAs were found in D0-vs-D2, D0-vs-D8, and D2-vs-D8, respectively. Functional analysis revealed that the differentially expressed lncRNAs were enriched in signaling pathways related to ovine preadipocyte differentiation, such as mitogen-activated protein kinase (MAPK) pathway, the phosphoinositide 3-kinase protein kinase B (PI3K-Akt) pathway, and the transforming growth factor beta (TGF-β) pathway. In summary, lncRNAs from preadipocytes at different stages of differentiation in sheep were identified and screened using RNA-Seq technology, and the regulatory mechanisms of lncRNAs in preadipocyte differentiation and lipid deposition were explored. This study provides a theoretical reference for revealing the roles of lncRNAs in ovine preadipocyte differentiation and also offers a theoretical basis for further understanding the regulatory mechanisms of ovine preadipocyte differentiation.

Long noncoding RNA transcriptome analysis reveals novel lncRNAs in Morus alba 'Yu-711' response to drought stress

Drought stress has been a key environmental factor affecting plant growth and development. The plant genome is capable of producing long noncoding RNAs (lncRNAs). To better understand white mulberry (Morus alba L.) drought response mechanism, we conducted a comparative transcriptome study comparing two treatments: drought-stressed (EG) and well-watered (CK) plants. A total of 674 differentially expressed lncRNAs (DElncRNAs) were identified. In addition, 782 differentially expressed messenger RNAs (DEmRNAs) were identified. We conducted Gene Ontology (GO) and KEGG enrichment analyses focusing on the differential lncRNAs cis-target genes. The target genes of the DElncRNAs were most significantly involved in the biosynthesis of secondary metabolites. Gene regulatory networks of the target genes involving DElncRNAs-mRNAs-DEmRNAs and DElncRNA-miRNA-DEmRNA were constructed. In the DElncRNAs-DEmRNAs network, 30 DEmRNAs involved in the biosynthesis of secondary metabolites are collocated with 46 DElncRNAs. The interaction between DElncRNAs and candidate genes was identified using LncTar. In summary, quantitative real-time polymerase chain reaction (qRT-PCR) validated nine candidate genes and seven target lncRNAs including those identified by LncTar. We predicted that the DElncRNAs-DEmRNAs might recruit microRNAs (miRNAs) to interact with gene regulatory networks under the drought stress response in mulberry. The findings will contribute to our understanding of the regulatory functions of lncRNAs under drought stress and will shed new light on the mulberry-drought stress interactions.

Multi-omics revealed rumen microbiota metabolism and host immune regulation in Tibetan sheep of different ages

The rumen microbiota and metabolites play an important role in energy metabolism and immune regulation of the host. However, the regulatory mechanism of rumen microbiota and metabolite interactions with host on Tibetan sheep's plateau adaptability is still unclear. We analyzed the ruminal microbiome and metabolome, host transcriptome and serum metabolome characteristics of Tibetan sheep at different ages. Biomarkers Butyrivibrio, Lachnospiraceae_XPB1014_group, Prevotella, and Rikenellaceae_RC9_gut_group were found in 4 months, 1.5 years, 3.5 years, and 6 years Tibetan sheep, respectively. The rumen microbial metabolites were mainly enriched in galactose metabolism, unsaturated fatty acid biosynthesis and fatty acid degradation pathways, and had significant correlation with microbiota. These metabolites further interact with mRNA, and are co-enriched in arginine and proline metabolism, metabolism of xenobiotics by cytochrome P450, propanoate metabolism, starch and sucrose metabolism, gap junction pathway. Meanwhile, serum metabolites also have a similar function, such as chemical carcinogenesis - reactive oxygen species, limonene and pinene degradation, and cutin, suberine and wax biosynthesis, thus participating in the regulation of the body's immune and energy-related metabolic processes. This study systematically revealed that rumen microbiota, metabolites, mRNA and serum metabolites of Tibetan sheep were involved in the regulation of fermentation metabolic function and immune level of Tibetan sheep at different ages, which provided a new perspective for plateau adaptability research of Tibetan sheep at different ages.

Identifying long non-coding RNAs involved in heat stress response during wheat pollen development

Introduction

Wheat is a staple food crop for over one-third of the global population. However, the stability of wheat productivity is threatened by heat waves associated with climate change. Heat stress at the reproductive stage can result in pollen sterility and failure of grain development.

Methods

This study used transcriptome data analysis to explore the specific expression of long non-coding RNAs (lncRNAs) in response to heat stress during pollen development in four wheat cultivars.

Results and discussion

We identified 11,054 lncRNA-producing loci, of which 5,482 lncRNAs showed differential expression in response to heat stress. Heat-responsive lncRNAs could target protein-coding genes in cis and trans and in lncRNA-miRNA-mRNA regulatory networks. Gene ontology analysis predicted that target protein-coding genes of lncRNAs regulate various biological processes such as hormonal responses, protein modification and folding, response to stress, and biosynthetic and metabolic processes. We also noted some paired lncRNA/protein-coding gene modules and some lncRNA-miRNA-mRNA regulatory modules shared in two or more wheat cultivars. These modules were related to regulating plant responses to heat stress, such as heat-shock proteins and transcription factors, and protein domains, such as MADS-box, Myc-type, and Alpha crystallin/Hsp20 domain.

Conclusion

Our results provide the basic knowledge and molecular resources for future functional studies investigating wheat reproductive development under heat stress.

Identification of lncRNA and weighted gene coexpression network analysis of germinating Rhizopus delemar causing mucormycosis

Rhizopus delemar, an opportunistic fungal pathogen, causes a highly fatal disease, mucormycosis. Spore germination is a crucial mechanism for disease pathogenesis. Thus, exploring the molecular mechanisms of fungal germination would underpin our knowledge of such transformation and, in turn, help control mucormycosis. To gain insight into the developmental process particularly associated with cell wall modification and synthesis, weighted gene co-expression network analysis (WGCNA) was performed including both coding and non-coding transcripts identified in the current study, to find out the module of interest in the germination stages. The module-trait relationship identified a particular module to have a high correlation only at the resting phase and further analysis revealed the module to be enriched for protein phosphorylation, carbohydrate metabolic process, and cellular response to stimulus. Moreover, co-expression network analysis of highly connected nodes revealed cell wall modifying enzymes, especially those involved in mannosylation, chitin-glucan crosslinking, and polygalacturonase activities co-expressing and interacting with the novel lncRNAs among which some of them predicted to be endogenous target mimic (eTM) lncRNAs. Hence, the present study provides an insight into the onset of spore germination and the information on the novel non-coding transcripts with key cell wall-related enzymes as potential targets against mucormycosis.

LncRTPred: Predicting RNA-RNA mode of interaction mediated by lncRNA

Long non-coding RNAs (lncRNAs) play a significant role in various biological processes. Hence, it is utmost important to elucidate their functions in order to understand the molecular mechanism of a complex biological system. This versatile RNA molecule has diverse modes of interaction, one of which constitutes lncRNA-mRNA interaction. Hence, identifying its target mRNA is essential to understand the function of an lncRNA explicitly. Existing lncRNA target prediction tools mainly adopt thermodynamics approach. Large execution time and inability to perform real-time prediction limit their usage. Further, lack of negative training dataset has been a hindrance in the path of developing machine learning (ML) based lncRNA target prediction tools. In this work, we have developed a ML-based lncRNA-mRNA target prediction model- 'LncRTPred'. Here we have addressed the existing problems by generating reliable negative dataset and creating robust ML models. We have identified the non-interacting lncRNA and mRNAs from the unlabelled dataset using BLAT. It is further filtered to get a reliable set of outliers. LncRTPred provides a cumulative_model_score as the final output against each query. In terms of prediction accuracy, LncRTPred outperforms other popular target prediction protocols like LncTar. Further, we have tested its performance against experimentally validated disease-specific lncRNA-mRNA interactions. Overall, performance of LncRTPred is heavily dependent on the size of the training dataset, which is highly reflected by the difference in its performance for human and mouse species. Its performance for human species shows better as compared to that for mouse when applied on an unknown data due to smaller size of the training dataset in case of mouse compared to that of human. Availability of increased number of lncRNA-mRNA interaction data for mouse will improve the performance of LncRTPred in future. Both webserver and standalone versions of LncRTPred are available. Web server link: http://bicresources.jcbose.ac.in/zhumur/lncrtpred/index.html. Github Link: https://github.com/zglabDIB/LncRTPred.

LncHLEF promotes hepatic lipid synthesis through miR-2188-3p/GATA6 axis and encoding peptides and enhances intramuscular fat deposition via exosome

Long noncoding RNAs (lncRNAs) have emergingly been implicated in mammalian lipid metabolism. However, their biological functions and regulatory mechanisms underlying adipogenesis remain largely elusive in chicken. Here, we systematically characterized the genome-wide full-length lncRNAs in the livers of pre- and peak-laying hens, and identified a novel intergenic lncRNA, lncHLEF, an RNA macromolecule with a calculated molecular weight of 433 kDa. lncHLEF was primarily distributed in cytoplasm of chicken hepatocyte and significantly up-regulated in livers of peak-laying hens. Functionally, lncHLEF could promote hepatocyte lipid droplet formation, triglycerides and total cholesterol contents. Mechanistically, lncHLEF could not only serve as a competitive endogenous RNA to modulate miR-2188-3p/GATA6 axis, but also encode three small functional polypeptides that directly interact with ACLY protein to enable its stabilization. Importantly, adeno-associated virus-mediated liver-specific lncHLEF overexpression resulted in increased hepatic lipid synthesis and intramuscular fat (IMF) deposition, but did not alter abdominal fat (AbF) deposition. Furthermore, hepatocyte lncHLEF could be delivered into intramuscular and abdominal preadipocytes via hepatocyte-secreted exosome to enhance intramuscular preadipocytes differentiation without altering abdominal preadipocytes differentiation. In conclusion, this study revealed that the lncHLEF could promote hepatic lipid synthesis through two independent regulatory mechanisms, and could enhance IMF deposition via hepatocyte-adipocyte communications mediated by exosome.

The novel RNA-RNA activation of H19 on MyoD transcripts promoting myogenic differentiation of goat muscle satellite cells

The elaborate interplay of coding and noncoding factors governs muscle growth and development. Here, we reported a mutual activation between long noncoding RNA (lncRNA) H19 and MyoD (myogenic determination gene number 1) in the muscle process. We successfully cloned the two isoforms of goat H19, which were significantly enriched and positively correlated with MyoD transcripts in skeletal muscles or differentiating muscle satellite cells (MuSCs). To systematically screen genes altered by H19, we performed RNA-seq using cDNA libraries of differentiating H19-deficiency MuSCs and consequently anchored MyoD as the critical genes in mediating H19 function. Intriguingly, some transcripts of MyoD and H19 overlapped in the cytoplasm, which was dramatically damaged when the core complementary nucleotides were mutated. Meanwhile, MyoD RNA successfully pulled down H19 in MS2-RIP experiments. Furthermore, HuR could bind both H19 and MyoD transcripts, while H19 or its truncated mutants successfully stabilized MyoD mRNA, with or without HuR deficiency. In turn, novel functional MyoD protein-binding sites were identified in the promoter and exons of the H19 gene. Our results suggest that MyoD activates H19 transcriptionally, and RNA-RNA hybridization is critical for H19-promoted MyoD expression, which extends our knowledge of the hierarchy of regulatory networks in muscle growth.

LncRNA MIR200CHG inhibits EMT in gastric cancer by stabilizing miR-200c from target-directed miRNA degradation

Gastric cancer (GC) is a heterogeneous disease, threatening millions of lives worldwide, yet the functional roles of long non-coding RNAs (lncRNAs) in different GC subtypes remain poorly characterized. Microsatellite stable (MSS)/epithelial-mesenchymal transition (EMT) GC is the most aggressive subtype associated with a poor prognosis. Here, we apply integrated network analysis to uncover lncRNA heterogeneity between GC subtypes, and identify MIR200CHG as a master regulator mediating EMT specifically in MSS/EMT GC. The expression of MIR200CHG is silenced in MSS/EMT GC by promoter hypermethylation, associated with poor prognosis. MIR200CHG reverses the mesenchymal identity of GC cells in vitro and inhibits metastasis in vivo. Mechanistically, MIR200CHG not only facilitates the biogenesis of its intronic miRNAs miR-200c and miR-141, but also protects miR-200c from target-directed miRNA degradation (TDMD) through direct binding to miR-200c. Our studies reveal a landscape of a subtype-specific lncRNA regulatory network, providing clinically relevant biological insights towards MSS/EMT GC.

Co-expression analysis suggests lncRNA-mRNA interactions enhance antiviral immune response during acute Chikungunya fever in whole blood of pediatric patients

Chikungunya virus is an arbovirus that causes the neglected tropical disease chikungunya fever, common in tropical areas worldwide. There is evidence that arboviruses alter host transcriptome and modulate immune response; this modulation may involve transcriptional and post-transcriptional control mechanisms mediated by long non-coding RNAs (lncRNAs). Herein, we employed bioinformatic analysis to evaluate co-expression of lncRNAs and their putative target mRNAs in whole blood during natural Chikungunya infection in adolescent boys. Sequencing data from GSE99992 was uploaded to the Galaxy web server, where data was aligned with HISAT2, gene counts were estimated with HTSeq-count, and differential expression was run with DESeq2. After gene classification with Biomart, Pearson's correlation was applied to identify potential interactions between lncRNAs and mRNAs, which were later classified into cis and trans according to genomic location (FEELnc) and binding potential (LncTar), respectively. We identified 1,975 mRNAs and 793 lncRNAs that were differentially expressed between the acute and convalescent stages of infection in the blood. Of the co-expressed lncRNAs and mRNAs, 357 potentially interact in trans and 9 in cis; their target mRNAs enriched pathways related to immune response and viral infections. Out of 52 enriched KEGG pathways, the RIG-I like receptor signaling is enriched by the highest number of target mRNAs. This pathway starts with the recognition of viral pathogens, leading to innate immune response mediated by the production of IFN-I and inflammatory cytokines. Our findings indicate that alterations in lncRNA expression in adolescent boys, induced by acute Chikungunya infection, potentially modulate mRNAs that contribute to antiviral immune responses.

The expression landscape and pangenome of long non-coding RNA in the fungal wheat pathogen Zymoseptoria tritici

Long non-coding RNAs (lncRNAs) are regulatory molecules interacting in a wide array of biological processes. lncRNAs in fungal pathogens can be responsive to stress and play roles in regulating growth and nutrient acquisition. Recent evidence suggests that lncRNAs may also play roles in virulence, such as regulating pathogenicity-associated enzymes and on-host reproductive cycles. Despite the importance of lncRNAs, only a few model fungi have well-documented inventories of lncRNA. In this study, we apply a recent computational pipeline to predict high-confidence lncRNA candidates in Zymoseptoria tritici, an important global pathogen of wheat impacting global food production. We analyse genomic features of lncRNAs and the most likely associated processes through analyses of expression over a host infection cycle. We find that lncRNAs are frequently expressed during early infection, before the switch to necrotrophic growth. They are mostly located in facultative heterochromatic regions, which are known to contain many genes associated with pathogenicity. Furthermore, we find that lncRNAs are frequently co-expressed with genes that may be involved in responding to host defence signals, such as oxidative stress. Finally, we assess pangenome features of lncRNAs using four additional reference-quality genomes. We find evidence that the repertoire of expressed lncRNAs varies substantially between individuals, even though lncRNA loci tend to be shared at the genomic level. Overall, this study provides a repertoire and putative functions of lncRNAs in Z. tritici enabling future molecular genetics and functional analyses in an important pathogen.

Transcriptome analysis reveals a lncRNA-miRNA-mRNA regulatory network in OsRpp30-mediated disease resistance in rice

BACKGROUND

Long non-coding RNAs (lncRNAs) play critical roles in various biological processes in plants. Extensive studies utilizing high-throughput RNA sequencing have revealed that many lncRNAs are involved in plant disease resistance. Oryza sativa RNase P protein 30 (OsRpp30) has been identified as a positive regulator of rice immunity against fungal and bacterial pathogens. Nevertheless, the specific functions of lncRNAs in relation to OsRpp30-mediated disease resistance in rice remain elusive.

RESULTS

We conducted a comprehensive analysis of lncRNAs, miRNAs, and mRNAs expression patterns in wild type (WT), OsRpp30 overexpression (OsRpp30-OE), and OsRpp30 knockout (OsRpp30-KO) rice plants. In total, we identified 91 differentially expressed lncRNAs (DElncRNAs), 1671 differentially expressed mRNAs (DEmRNAs), and 41 differentially expressed miRNAs (DEmiRNAs) across the different rice lines. To gain further insights, we investigated the interaction between DElncRNAs and DEmRNAs, leading to the discovery of 10 trans- and 27 cis-targeting pairs specific to the OsRpp30-OE and OsRpp30-KO samples. In addition, we constructed a competing endogenous RNA (ceRNA) network comprising differentially expressed lncRNAs, miRNAs, and mRNAs to elucidate their intricate interplay in rice disease resistance. The ceRNA network analysis uncovered a set of gene targets regulated by lncRNAs and miRNAs, which were found to be involved in pathogen recognition, hormone pathways, transcription factor activation, and other biological processes related to plant immunity.

CONCLUSIONS

Our study provides a comprehensive expression profiling of lncRNAs, miRNAs, and mRNAs in a collection of defense mutants in rice. To decipher the putative functional significance of lncRNAs, we constructed trans- and cis-targeting networks involving differentially expressed lncRNAs and mRNAs, as well as a ceRNA network incorporating differentially expressed lncRNAs, miRNAs, and mRNAs. Together, the findings from this study provide compelling evidence supporting the pivotal roles of lncRNAs in OsRpp30-mediated disease resistance in rice.

Regulation of cashmere fineness traits by noncoding RNA in Jiangnan cashmere goats

BACKGROUND

Cashmere has long been used as the raw material for wool textiles. The diameter of the cashmere fibre determines its quality and economic value. However, the regulatory role of noncoding RNAs (ncRNAs) in cashmere fineness remains unclear, especially regarding the interaction between ncRNAs and coding RNAs.

RESULTS

Transcriptome sequencing was used to identify the expression profiles of long noncoding RNAs (lncRNAs), circular RNAs (circRNAs) and microRNAs (miRNAs) in the skin tissues of Jiangnan cashmere goats with different cashmere fineness levels. Integration analysis of ncRNA and coding RNA was performed in combination with previous research results. The results showed that 16,437 lncRNAs, 2234 circRNAs, and 1322 miRNAs were identified in 8 skin samples of cashmere goats. A total of 403 differentially expressed (DE) lncRNAs, 62 DE circRNAs and 30 DE miRNAs were identified in the skin tissues of the fine groups (Fe) and coarse groups (Ce). We predicted the target gene of DE lncRNA, the target gene of DE miRNA and the host gene of DE circRNA. Based on functional annotation and enrichment analysis of target genes, we found that DE lncRNAs could be involved in regulating the fineness traits of cashmere. The most potential lncRNAs were MSTRG.42054.1, MSTRG.18602.3, and MSTRG.2199.13.

CONCLUSIONS

The data from this study enriched the cashmere goat noncoding RNA database and helped to supplement the annotation of the goat genome. The results provided a new direction for the breeding of cashmere characters.

Synergistic Responses of Tibetan Sheep Rumen Microbiota, Metabolites, and the Host to the Plateau Environment

Plateau adaptation in animals involves genetic mechanisms as well as coevolutionary mechanisms of the microbiota and metabolome of the animal. Therefore, the characteristics of the rumen microbiome and metabolome, transcriptome, and serum metabolome of Tibetan sheep at different altitudes (4500 m, 3500 m, and 2500 m) were analyzed. The results showed that the rumen differential metabolites at 3500 m and 4500 m were mainly enriched in amino acid metabolism, lipid metabolism, and carbohydrate metabolism, and there was a significant correlation with microbiota. The differentially expressed genes and metabolites at middle and high altitudes were coenriched in asthma, arachidonic acid metabolism, and butanoate and propanoate metabolism. In addition, the serum differential metabolites at 3500 m and 4500 m were mainly enriched in amino acid metabolism, lipid metabolism, and metabolism of xenobiotics by cytochrome P450, and they were also related to microbiota. Further analysis revealed that rumen metabolites accounted for 7.65% of serum metabolites. These common metabolites were mainly enriched in metabolic pathways and were significantly correlated with host genes (p < 0.05). This study found that microbiota, metabolites, and epithelial genes were coenriched in pathways related to lipid metabolism, energy metabolism, and immune metabolism, which may be involved in the regulation of Tibetan sheep adaptation to plateau environmental changes.

Orchestration of a blood-derived and ADARB1-centred network in Alzheimer's and Parkinson's disease

The peripheral immune system is thought to influence the pathogenesis of the central nervous system in Alzheimer's disease (AD) and Parkinson's disease (PD). This study aimed to investigate the characteristics of peripheral leukocytes in AD and PD by comprehensively analyzing the transcriptomic and metabolic features in the blood (NCONTROL = 15; NAD = 11; NPD = 13). The study found an ADARB1-centered module that was associated with diagnosis, phenethylamine, and glutamate. The module consisted of ADARB1, a vital RNA-editing enzyme, LINC02960, and 109 miRNAs. The study also predicted that the ADARB1 and involved regulators were targeted by miRNAs in the ADARB1 module. The integrated analysis of transcriptome and metabolic panel revealed a central role of ADARB1, miR-199b-5p, miR-26a, miR-450b-5p, miR-34c-5p, glutamate and phenethylamine in the regulatory relationships. The study highlights a set of synergetic non-coding RNA related to ADARB1 in the blood ecosystem of AD and PD with dynamic glutamate and phenethylamine, providing new insights into the pathogenesis of these diseases.

Integrative time-serial networks for genome-wide lncRNA-mRNA interactions reveal interferon-inducible antiviral and T-cell receptor regulations against PRRSV infection

Porcine reproductive and respiratory syndrome virus (PRRSV) infection severely affects the swine industry each year. Although the host mechanisms against PRRSV infection have been identified in key target tissues through whole transcriptome sequencing, specific molecular regulators have not been elucidated. Long non-coding RNA (lncRNA) expression is highly specific and could thus be used to effectively identify PRRSV-specific candidates. Here, we identified novel lncRNAs in lungs, bronchial lymph nodes, and tonsils after PRRSV infection and constructed phenotype-based integrative co-expression networks using time-series differentially expressed (DE) lncRNAs and mRNAs. After the analyses, a total of 309 lncRNA-mRNA interactions were identified. During early host innate signalling, interferon-inducible and interferon genes were positively regulated by specific lncRNA. Moreover, T-cell receptor genes in lung adaptive immune signalling were negatively regulated by specific lncRNA. Collectively, our findings provide insights into the genome-wide lncRNA-mRNA interactions and dynamic regulation of lncRNA-mediated mechanisms against PRRSV infection.

A New Understanding of Long Non-Coding RNA in Hepatocellular Carcinoma-From m6A Modification to Blood Biomarkers

With recent advancements in biological research, long non-coding RNAs (lncRNAs) with lengths exceeding 200 nucleotides have emerged as pivotal regulators of gene expression and cellular phenotypic modulation. Despite initial skepticism due to their low sequence conservation and expression levels, their significance in various biological processes has become increasingly apparent. We provided an overview of lncRNAs and discussed their defining features and modes of operation. We then explored their crucial function in the hepatocarcinogenesis process, elucidating their complex involvement in hepatocellular carcinoma (HCC). The influential role of lncRNAs within the HCC tumor microenvironment is emphasized, illustrating their potential as key modulators of disease dynamics. We also investigated the significant influence of N6-methyladenosine (m6A) modification on lncRNA function in HCC, enhancing our understanding of both their roles and their upstream regulators. Additionally, the potential of lncRNAs as promising biomarkers was discussed in liver cancer diagnosis, suggesting a novel avenue for future research and clinical application. Finally, our work underscored the dual potential of lncRNAs as both contributors to HCC pathogenesis and innovative tools for its diagnosis. Existing challenges and prospective trajectories in lncRNA research are also discussed, emphasizing their potential in advancing liver cancer research.

Integrated co-expression analysis of regulatory elements (miRNA, lncRNA, and TFs) in bovine monocytes induced by Str. uberis

Non-coding RNAs, including long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), together with transcription factors, are critical pre-, co-, and post-transcriptional regulators. In addition to their criteria as ideal biomarkers, they have great potential in disease prognosis, diagnosis, and treatment of complex diseases. Investigation of regulatory mechanisms in the context of bovine mastitis, as most common and economic disease in the dairy industry, to identify elements influencing the expression of candidate genes as key regulators of the mammary immune response is not yet fully understood. Transcriptome profiles (50 RNA-Seq and 50 miRNA-Seq samples) of bovine monocytes induced by Str. uberis were used for co-expression module detection and preservation analysis using the weighted gene co-expression network analysis (WGCNA) approach. Assigned mi-, lnc-, and m-modules used to construct the integrated regulatory networks and miRNA-lncRNA-mRNA regulatory sub-networks. Remarkably, we have identified 18 miRNAs, five lncRNAs, and seven TFs as key regulators of str. uberis-induced mastitis. Most of the genes introduced here, mainly involved in immune response, inflammation, and apoptosis, were new to mastitis. These findings may help to further elucidate the underlying mechanisms of bovine mastitis, and the discovered genes may serve as signatures for early diagnosis and treatment of the disease.

Identification of a novel immune-related long noncoding RNA in carp primary macrophages associated with bisphenol A' s immunoregulatory effects

Increasing evidence suggests that long non-coding RNAs (lncRNAs) play pivotal roles in various biological processes. However, current studies on lncRNAs mostly focus on mammalian species, with little research on the functional roles of lncRNAs in teleost fish. Here, we identified a novel intergenic lncRNA (linc-93.2) in the head kidney primary macrophages of common carp (Cyprinus carpio) after exposure to a typical environmental endocrine disrupting chemical, bisphenol A (BPA). As a result, linc-93.2 was more than 3,619 bp in length and predominantly localized to the nucleus of primary macrophages other than cytoplasm, with the highest expression level in spleen followed by head kidney among different organs. Bioinformatic analysis predicted a cis-target gene, dennd1b, and 20 trans-target genes including hsp70, gna13 and rasgap, were potentially regulated by linc-93.2; NFκB and estrogen receptor (ERα) binding sites were located in the promoter region upstream of its transcription start site, which together suggested the involvement of linc-93.2 in immune and neurological functions in fish. Based on that, the expression level of linc-93.2 was determined in macrophages following acute lipopolysaccharide (LPS) and BPA treatments, both of which significantly induced linc-93.2 and IL-1β expression in cells. Moreover, a NF-κB inhibitor PDTC significantly reduced linc-93.2 expression in macrophages, but co-exposure of macrophages to PDTC with BPA or LPS could significantly rescue linc-93.2 expression, consistent with the observation on that LPS or BPA alone significantly induced both linc-93.2 and its target gene expression. Interestingly, linc-93.2 and its target gene expression was significantly suppressed by an ER antagonist ICI 182,780, however, the co-exposure of macrophages to ICI 182,780 with BPA failed to attenuate their declined expression. Overall, the current study demonstrated that linc-93.2, a novel immune-related lncRNA, may participate in the immune processes of common carp macrophages via the NF-κB and ER pathway. The results presented in this study enhance our understanding of the immunotoxin mechanisms of BPA in teleost fish.

VirhostlncR: A comprehensive database to explore lncRNAs and their targets in viral infections

Long non-coding-RNAs (lncRNAs) are an expanding set of cis-/trans-regulatory RNA genes that outnumber the protein-coding genes. Although being increasingly discovered, the functional role of the majority of lncRNAs in diverse biological conditions is undefined. Increasing evidence supports the critical role of lncRNAs in the emergence, regulation, and progression of various viral infections including influenza, hepatitis, coronavirus, and human immunodeficiency virus. Hence, the identification of signature lncRNAs would facilitate focused analysis of their functional roles accounting for their targets and regulatory mechanisms associated with infections. Towards this, we compiled 2803 lncRNAs identified to be modulated by 33 viral strains in various mammalian cell types and are provided through the resource named VirhostlncR (http://ciods.in/VirhostlncR/). The information on each of the viral strains, their multiplicity of infection, duration of infection, host cell name and cell types, fold change of lncRNA expression, and their specific identification methods are integrated into VirhostlncR. Based on the current datasets, we report 150 lncRNAs including differentiation antagonizing non-protein coding RNA (DANCR), metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), maternally expressed gene 3 (MEG3), nuclear paraspeckle assembly transcript 1 (NEAT1), and plasmacytoma variant translocation 1 (PVT1) to be perturbed by two or more viruses. Analysis of viral protein interactions with human transcription factors (TFs) or TF-containing protein complexes identified that distinct viruses can transcriptionally regulate many of these lncRNAs through multiple protein complexes. Together, we believe that the current dataset will enable priority selection of lncRNAs for identification of their targets and serve as an effective platform for the analysis of noncoding RNA-mediated regulations in viral infections.

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

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

Molecular regulatory mechanism of key LncRNAs in subclinical mastitic cows with folic acid supplementation

BACKGROUND

Folic acid is a water-soluble B vitamin (B9), which is closely related to the body's immune and other metabolic pathways. The folic acid synthesized by rumen microbes has been unable to meet the needs of high-yielding dairy cows. The incidence rate of subclinical mastitis in dairy herds worldwide ranged between 25%~65% with no obvious symptoms, but it significantly causes a decrease in lactation and milk quality. Therefore, this study aims at exploring the effects of folic acid supplementation on the expression profile of lncRNAs, exploring the molecular mechanism by which lncRNAs regulate immunity in subclinical mastitic dairy cows.

RESULTS

The analysis identified a total of 4384 lncRNA transcripts. Subsequently, differentially expressed lncRNAs in the comparison of two groups (SF vs. SC, HF vs. HC) were identified to be 84 and 55 respectively. Furthermore, the weighted gene co-expression network analysis (WGCNA) and the KEGG enrichment analysis result showed that folic acid supplementation affects inflammation and immune response-related pathways. The two groups have few pathways in common. One important lncRNA MSTRG.11108.1 and its target genes (ICAM1, CCL3, CCL4, etc.) were involved in immune-related pathways. Finally, through integrated analysis of lncRNAs with GWAS data and animal QTL database, we found that differential lncRNA and its target genes could be significantly enriched in SNPs and QTLs related to somatic cell count (SCC) and mastitis, such as MSTRG.11108.1 and its target gene ICAM1, CXCL3, GRO1.

CONCLUSIONS

For subclinical mastitic cows, folic acid supplementation can significantly affect the expression of immune-related pathway genes such as ICAM1 by regulating lncRNAs MSTRG.11108.1, thereby affecting related immune phenotypes. Our findings laid a ground foundation for theoretical and practical application for feeding folic acid supplementation in subclinical mastitic cows.

MALAT1 in Liquid Biopsy: The Diagnostic and Prognostic Promise for Colorectal Cancer and Adenomas?

Introduction

The development of colorectal cancer (CRC) is a multistep process accompanied by the accumulation of mutations that start from specific precancerous lesion - colorectal adenomas (CA). CRC incidence and mortality can be reduced by the early identification of these neoplasm. Colonoscopy is the most widely used screening method for CRC identification. Nowadays, clinical research interest is shifting to the use of liquid biopsy that may help with the early diagnosis of CA and CRC. In our previous study, we identified long non-coding RNA MALAT1 gene amplification associated with the development of CA.

Methods

This study aimed to describe the potential of MALAT1 expression levels in the adenoma tissue of patients used in the previous study by real-time qPCR. Furthermore, we analysed the plasma samples of an independent group of patients with CA (n=97), CRC (n=101), and cancer-free individuals (CFI, n=48).

Results

There was no difference in the MALAT1 expression level between CA patients with or without MALAT1 amplification. However, the plasma MALAT1 expression levels were significantly upregulated in patients with CRC and CA compared to CFI (for both p<0.001). Moreover, a correlation between MALAT1 expression and histological types of adenomas was identified- high-CRC-risk adenomas also displayed the highest MALAT1 expression levels. Furthermore, in CRC patients, MALAT1 levels were associated with a response to therapy.

Conclusion

MALAT1 expression levels could serve as a promising circulating biomarker for early CA and CRC diagnosis, and even as a predictor of therapy response in CRC patients.

Construction of chronic glomerulonephritis‑related lncRNA‑mRNA regulatory network and lncRNA‑-miRNA‑mRNA ceRNA network by bioinformatics analysis

Long non-coding RNAs (lncRNAs) are ncRNA transcripts >200 nucleotides that are important genetic regulators. LncRNAs can directly regulate mRNA through a lncRNA-mRNA regulatory mode and can also regulate mRNA through competitive binding to micro (mi)RNA, which is generally known as the competitive endogenous RNA (ceRNA) network. The present study evaluated the functional roles and regulatory networks of lncRNAs in chronic glomerulonephritis (CGN). The proliferative ability of mouse glomerular mesangial cells (GMCs) induced by different concentrations of lipopolysaccharide (LPS) was assessed using the Cell Counting Kit-8 assay, and RNA sequencing (RNA-seq) was performed to identify differentially expressed lncRNAs in LPS-induced GMCs. Based on the sequencing results, six lncRNAs were selected for validation using reverse transcription-quantitative PCR (RT-qPCR). Furthermore, the lncRNA-mRNA regulatory network and the lncRNA-miRNA-mRNA ceRNA network were constructed to assess the role and mechanism of CGN-related lncRNAs. To elucidate the biological functions of lncRNAs, Gene Ontology (GO) biological process term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on all mRNAs involved in the lncRNA-mRNA regulatory network and in the ceRNA network. A total of 1,532 differentially expressed lncRNAs, including 594 upregulated lncRNAs and 938 downregulated lncRNAs, were identified using RNA-seq. The results of RT-qPCR validation were consistent with RNA-seq results. An lncRNA-mRNA regulatory network, including 236 lncRNAs and 556 mRNAs, and a ceRNA network, including 6 lncRNAs, 18 miRNAs and 419 mRNAs, were successfully constructed. The GO biological process term enrichment and KEGG pathway enrichment analyses demonstrated that those lncRNAs were often related to inflammatory response and substance metabolism. The present study identified key CGN-related lncRNAs in LPS-induced GMCs, and further demonstrated a global view of the lncRNA-mRNA regulatory network and ceRNA network involved in CGN. These results offered novel insights into the roles of lncRNAs in the pathogenesis of CGN and identified potential diagnostic biomarkers.