Differentially Expressed lncRNAs Research Articles

Changes of mRNA, miRNA and lncRNA expression contributing to skeletal muscle differences between fetus and adult Mongolian horses

The growth and development of myofibers, as the fundamental units comprising muscle tissue, and their composition type are indeed among the most crucial factors influencing skeletal muscle types. Muscle fiber adaptation is closely associated with alterations in physiological conditions. Muscle fiber types undergo dynamic changes in fetus and adult horses. Our aim is to investigate the mechanisms influencing the differences in muscle fiber types between fetal and adult stages of Mongolian horses. The study investigated the distribution of muscle fiber types within longissimus dorsi muscle of fetus and adult Mongolian horses. A total of 652 differentially expressed genes (DEGs), 476 Differentially expressed lncRNAs (DELs), and 174 Differentially expressed miRNAs (DEMIRs) were identified using deep RNA-seq analysis. The results of functional analysis reveal the transformations in muscle fiber type from the fetal to adult stage in Mongolian horses. The up-regulated DEGs were implicated in the development and differentiation of muscle fibers, while down-regulated DEGs were associated with muscle fiber contraction, transformation, and metabolism. Additionally, connections between non-coding RNA and mRNA landscapes were identified based on their functional alterations, some non-coding RNA target genes may be associated with immunity. These data have broadened our understanding of the specific roles and interrelationships among regulatory molecules involved in Mongolian horse development, this provides new perspectives for selecting and breeding superior individuals and for disease prevention.

LncRNA cis- and trans-regulation provides new insight into drought stress responses in wild barley.

Drought is one of the most common abiotic stresses that affect barley productivity. Long noncoding RNA (lncRNA) has been reported to be widely involved in abiotic stress, however, its function in the drought stress response in wild barley remains unclear. In this study, RNA sequencing was performed to identify differentially expressed lncRNAs (DElncRNA) among two wild barley and two cultivated barley genotypes. Then, the cis-regulatory networks were according to the chromosome position and the expression level correction. The GO annotation indicates that these cis-target genes are mainly involved in "ion transport transporter activity" and "metal ion transport transporter activity". Through weighted gene co-expression network analysis (WGCNA), 10 drought-related modules were identified to contract trans-regulatory networks. The KEGG annotation demonstrated that these trans-target genes were enriched for photosynthetic physiology, brassinosteroid biosynthesis, and flavonoid metabolism. In addition, we constructed the lncRNA-mediated ceRNA regulatory network by predicting the microRNA response elements (MREs). Furthermore, the expressions of lncRNAs were verified by RT-qPCR. Functional verification of a candidate lncRNA, MSTRG.32128, demonstrated its positive role in drought response and root growth and development regulation. Hormone content analysis provided insights into the regulatory mechanisms of MSTRG.32128 in root development, revealing its involvement in auxin and ethylene signal transduction pathways. These findings advance our understanding of lncRNA-mediated regulatory mechanisms in barley under drought stress. Our results will provide new insights into the functions of lncRNAs in barley responding to drought stress.

Transcriptome analysis revealed the characteristics and functions of long non-coding RNAs in the hypothalamus during sexual maturation in goats

The hypothalamus is an essential neuroendocrine area in animals that regulates sexual development. Long non-coding RNAs (lncRNAs) are hypothesized to regulate physiological processes related to animal reproduction. However, the regulatory mechanism by which lncRNAs participate in sexual maturity in goats is poorly known, particularly from birth to sexual maturation. In this study, RNAseq analysis was conducted on the hypothalamus of four developmental stages (1day (D1, n = 5), 2 months (M2, n = 5), 4 months (M4, n = 5), and 6 months (M6, n = 5)) of Jining grey goats. The results showed that a total of 237 differentially expressed lncRNAs (DELs) were identified in the hypothalamus. Among these, 221 DELs exhibited cis-regulatory effects on 693 target genes, while 24 DELs demonstrated trans-regulatory effects on 63 target genes. The target genes of these DELs are mainly involved in biological processes related to energy metabolism, signal transduction and hormone secretion, such as sphingolipid signaling pathway, adipocytokine signaling pathway, neurotrophic signaling pathway, glutamatergic synapse, P53 signaling pathway and GnRH signaling pathway. In addition, XR_001918477.1, TCONS_00077463, XR_001918760.1, and TCONS_00029048 and their potential target genes may play a crucial role in the process of goat sexual maturation. This study advances our understanding of lncRNA in hypothalamic tissue during sexual maturation in goats and will give a theoretical foundation for improving goat reproductive features.

Integrative transcriptomic profiling of ncRNAs and mRNAs in developing mouse lens.

In recent years, burgeoning research has underscored the pivotal role of non-coding RNA in orchestrating the growth, development, and pathogenesis of various diseases across organisms. However, despite these advances, our understanding of the specific contributions of long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) to lens development remains notably limited. Clarifying the intricate gene regulatory networks is imperative for unraveling the molecular underpinnings of lens-related disorders. In this study, we aimed to address this gap by conducting a comprehensive analysis of the expression profiles of messenger RNAs (mRNAs), lncRNAs, and circRNAs at critical developmental time points of the mouse lens, encompassing both embryonic (E10.5, E12.5, and E16.5) and postnatal stages (P0.5, P10.5, and P60). Leveraging RNA-sequencing technology, we identified key transcripts pivotal to lens development. Our analysis revealed differentially expressed (DE) mRNAs, lncRNAs, and circRNAs across various developmental stages. Particularly noteworthy, there were 1831 co-differentially expressed (CO-DE) mRNAs, 150 CO-DE lncRNAs, and 13 CO-DE circRNAs identified during embryonic stages. Gene Ontology (GO) enrichment analysis unveiled associations primarily related to lens development, DNA conformational changes, and angiogenesis among DE mRNAs and lncRNAs. Furthermore, employing protein-protein interaction networks, mRNA-lncRNA co-expression networks, and circRNA-microRNA-mRNA networks, we predicted candidate key molecules implicated in lens development. Our findings underscore the pivotal roles of lncRNAs and circRNAs in this process, offering fresh insights into the pathogenesis of lens-related disorders and paving the way for future exploration in this field.

Comprehensive analysis of the expression profiles of mRNA, lncRNA, circRNA, and miRNA in primary hair follicles of coarse sheep fetal skin

BackgroundThe Qinghai Tibetan sheep, a local breed renowned for its long hair, has experienced significant deterioration in wool characteristics due to the absence of systematic breeding practices. Therefore, it is imperative to investigate the molecular mechanisms underlying follicle development in order to genetically enhance wool-related traits and safeguard the sustainable utilization of valuable germplasm resources. However, our understanding of the regulatory roles played by coding and non-coding RNAs in hair follicle development remains largely elusive.ResultsA total of 20,874 mRNAs, 25,831 circRNAs, 4087 lncRNAs, and 794 miRNAs were annotated. Among them, we identified 58 DE lncRNAs, 325 DE circRNAs, 924 DE mRNAs, and 228 DE miRNAs during the development of medullary primary hair follicle development. GO and KEGG functional enrichment analyses revealed that the JAK-STAT, TGF-β, Hedgehog, PPAR, cGMP-PKG signaling pathway play crucial roles in regulating fibroblast and epithelial development during skin and hair follicle induction. Furthermore, the interactive network analysis additionally identified several crucial mRNA, circRNA, and lncRNA molecules associated with the process of primary hair follicle development. Ultimately, by investigating DEmir’s role in the ceRNA regulatory network mechanism, we identified 113 circRNA–miRNA pairs and 14 miRNA–mRNA pairs, including IGF2BP1-miR-23-x-novel-circ-01998-MSTRG.7111.3, DPT-miR-370-y-novel-circ-005802-MSTRG.14857.1 and TSPEAR-oar-miR-370-3p-novel-circ-005802- MSTRG.10527.1.ConclusionsOur study offers novel insights into the distinct expression patterns of various transcription types during hair follicle morphogenesis, establishing a solid foundation for unraveling the molecular mechanisms that drive hair development and providing a scientific basis for selectively breeding desirable wool-related traits in this specific breed.

Identifying autophagy-related mRNAs and potential ceRNA networks in meniscus degeneration based on RNA sequencing and experimental validation

PurposeThe intimate connection between long noncoding RNA (lncRNA) and autophagy has been established in cartilage degeneration. However, their roles in meniscal degeneration remain ambiguous. This study aimed to identify the key autophagy-related lncRNA and its associated regulatory network in meniscal degeneration in the context of osteoarthritis (OA). MethodsRNA sequencing was performed to identify differentially expressed lncRNAs (DELs) and mRNAs (DEMs), which were then conducted to enrichment analyses using the DAVID database and Metascape. Autophagy-related DEMs were identified by combining DEMs with data from the Human Autophagy Database. Three databases were used to predict miRNA, and the DIANA LncBase Predicted database was utilized to predict miRNA-lncRNA interactions. Based on these predictions, comprehensive competitive endogenous RNA (ceRNA) network were constructed. The expression levels of the classical autophagy markers and autophagy-related ceRNA network were validated. Additionally, Gene Set Enrichment Analysis (GSEA) was performed using autophagy-related DEMs. Results310 DELs and 320 DEMs were identified, with five upregulated and one downregulated autophagy-related DEMs. Through reverse prediction of miRNA, paired miRNA-lncRNA interactions, and verification using RT-qPCR, two lncRNAs (PCAT19, CLIP1-AS1), two miRNA (has-miR-3680-3p and has-miR-4795-3p) and two mRNAs (BAG3 and HSP90AB1) were included in the constructed ceRNA regulatory networks. GSEA indicated that the increased expression of autophagy-related mRNAs inhibited glycosaminoglycan biosynthesis in the degenerative meniscus. ConclusionThis study presented the first construction of regulatory ceRNA network involving autophagy-related lncRNA-miRNA-mRNA interactions in OA meniscus. These findings offered valuable insights into the mechanisms underlying meniscal degeneration and provided potential targets for therapeutic intervention.

Dynamic whole-transcriptome landscape of acute bilirubin encephalopathy in newborns

Hyperbilirubinemia in newborns may progress to acute bilirubin encephalopathy (ABE), posing short- and long-term health risks. Despite extensive research identifying numerous mRNAs, lncRNAs, circRNAs, and miRNAs associated with brain injury, their roles in neonatal bilirubin-induced brain injury remain elusive. This study employed whole-transcriptome sequencing to ascertain the differentially expressed (DE) RNA profiles in a newborn ABE rat model, followed by bioinformatic analysis. A time-series competing endogenous RNA (ceRNA) regulatory network was established, and the expression trends of 9 arbitrarily chosen RNAs were verified through quantitative real-time polymerase chain reaction(qRT-PCR). In comparison with the control group, we identified 595, 888, and 1448 DE mRNAs; 22, 37, and 37 DE miRNAs; 1945, 1869, and 1997 DE lncRNAs; and 31, 28, and 36 DE circRNAs at 6 h, 12 h, and 24 h, respectively. Predominantly, these DERNAs contribute to biological functions and pathways associated with inflammation, immunity, metabolism, cell death, and neurodevelopmental regulation. Moreover, we constructed ceRNA networks of DE lncRNA/circRNA-DE miRNA-DE mRNA based on time series. The qRT-PCR expression trends for the selected 9 RNAs were generally similar to the RNA-seq outcomes. This investigation uniquely delineated the temporal expression patterns of mRNA and non-coding RNA in ABE, establishing ceRNA networks and identifying potential molecular mechanisms underlying bilirubin-induced hippocampal damage. Nonetheless, further studies are warranted to corroborate these findings in humans.

Genome-wide analysis of long noncoding RNAs as cis-acting regulators of transcription factor-encoding genes in IgA nephropathy.

IgA nephropathy (IgAN) is the most common form of primary glomerulonephritis in the world, but the disease pathogenesis noncoding is yet to be elucidated. Previous studies have revealed regulatory functions for long noncoding RNA (lncRNA) in various diseases; however, the roles of lncRNA in IgAN and regulation of transcription factors (TFs) have been scarcely investigated. Renal tissue samples (n = 5) from patients with IgAN and control samples (n = 4) were collected and RNA sequencing (RNA-seq) was performed. Four software programs were employed for lncRNA prediction. GO (Gene Ontology)/KEGG (Kyoto Encyclopedia of Genes and Genomes) were employed for analysis of the identified differentially expressed genes (DEGs). A regulatory network model of DE lncRNA-TF-DEG was developed, and the levels of expression of key lncRNAs, TFs, and corresponding target genes were assessed using qRT-PCR and immunofluorescence. The current study identified 674 upregulated and 1,011 downregulated DE mRNAs and 260 upregulated and 232 downregulated DE lncRNAs in IgAN samples compared with control samples. The upregulated DE mRNAs showed enrichment in cell adhesion and collagen glial fiber organization pathways. The DE lncRNAs-DE mRNAs showing co-expression are associated with transmembrane transport. A novel regulatory network model of lncRNA-TF-DEG has been developed. This study identified seven TFs that are cis-regulated by 6 DE lncRNAs, and show co-expression with 132 DEGs (correlation coefficient ≥ 0.8, P ≤ 0.01), generating 158 pairs that showed co-expression. The lncRNAs NQO1-DT and RP5-1057120.6 were found to be highly expressed in IgAN samples. The TFs vitamin D Receptor (VDR) and NFAT5, along with their target genes were also aberrantly expressed. Key lncRNAs and TFs centrally associated with IgAN have been identified in this study. A regulatory network model of lncRNA-TF-mRNA was constructed. Further studies on the genes identified herewith could provide insight into the pathogenesis of IgAN.

Systematic analysis of co-expression network and gene expression profile of differentially expressed lncRNA-mRNA one week after spinal cord injury

Spinal cord injury (SCI) is a severe event that occurs in the central nervous system (CNS). Long noncoding RNAs (lncRNAs) have already been proven to play key roles in different diseases of the CNS. Nevertheless, the roles of messenger RNAs (mRNAs) and lncRNAs one week after SCI remain to be elucidated. We established a rat model of SCI and examined the differential levels of mRNAs and lncRNAs, and then, the differentially expressed mRNAs (DE mRNAs) and lncRNAs (DE lncRNAs) were investigated. Gene Ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DE lncRNAs and mRNAs were subsequently carried out to explore crucial biological processes and signaling pathways. Moreover, we constructed protein–protein interaction (PPI) networks, the STRING database, and identified the 10 high degree core nodes. We also established a co-expression network by calculating the correlation between lncRNA expression profile and mRNA expression profile. The identified DE lncRNAs and DE mRNAs may potentially be critical targets for the diagnosis and treatment of SCI in the future. Abbreviations: SCI: spinal cord injury; CNS: central nervous system; YLDs: years of life lived with disability; lncRNAs: long noncoding RNAs; mRNAs: messenger RNAs; DE: differentially expressed; GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; GSEA: gene set enrichment analysis; RIN: Renewable Identification Number; PCA: principal component analysis; PPI: protein–protein interaction; MF, molecular function; BP: biological process; CC: cellular component; MCODE: molecular complex detection

Transcriptomic profiling reveals salt-responsive long non-coding RNAs and their putative target genes for improving salt tolerance in upland cotton (Gossypium hirsutum)

Salt stress continues to exert notable impacts on the growth and development of cotton plants. In various biological processes of plants and in response to abiotic stress, long noncoding RNAs (lncRNAs) play a crucial regulatory role. In this study, we conducted transcriptome sequencing of cotton plants grown under salt stress and normal water conditions. A total of 519 lncRNAs and 75,376 messenger RNAs (mRNAs) were identified. Compared with mRNAs, lncRNAs have shorter transcript lengths and open reading frames, fewer exons, and lower expression levels. By comparing and analyzing the differentially expressed lncRNAs (DELs) and differentially expressed coding genes (DEGs) after salt treatment, we identified a total of 60 DELs and 217 DEGs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that these target genes are enriched primarily in intracellular organelles and participate in metabolic regulation pathways. Ultimately, lncRNA XLOC_043671, along with its target gene Gohir.D02G125200, and the lncRNA XLOC_034831, along with its target gene Gohir.A12G114800 was identified as a participant in salt tolerance regulation in cotton. RNA-seq and RTqPCR revealed that both these lncRNAs and their target genes exhibited significantly enhanced expression under salt stress. Subsequently, utilizing virus-induced gene silencing (VIGS) technology, we successfully silenced the two lncRNAs in cotton plants, which resulted in reduced expression levels of the corresponding target genes, rendering the plants less salt tolerant than did the control plants. Similarly, we successfully silenced these two target genes in cotton seedlings, and the silenced plants demonstrated a pronounced salt-sensitive phenotype. In summary, our transcriptome analysis revealed a significant number of differentially expressed lncRNAs and their corresponding coding genes. These findings lay the foundation for further investigating the mechanisms of lncRNAs and their target genes in the salt stress response of cotton.

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.

Whole-Transcriptome Analysis Sheds Light on the Biological Contexts of Intramuscular Fat Deposition in Ningxiang Pigs.

The quality of pork is significantly impacted by intramuscular fat (IMF). However, the regulatory mechanism of IMF depositions remains unclear. We performed whole-transcriptome sequencing of the longissimus dorsi muscle (IMF) from the high (5.1 ± 0.08) and low (2.9 ± 0.51) IMF groups (%) to elucidate potential mechanisms. In summary, 285 differentially expressed genes (DEGs), 14 differentially expressed miRNAs (DEMIs), 83 differentially expressed lncRNAs (DELs), and 79 differentially expressed circRNAs (DECs) were identified. DEGs were widely associated with IMF deposition and liposome differentiation. Furthermore, competing endogenous RNA (ceRNA) regulatory networks were constructed through co-differential expression analyses, which included circRNA-miRNA-mRNA (containing 6 DEMIs, 6 DEGs, 47 DECs) and lncRNA-miRNA-mRNA (containing 6 DEMIs, 6 DEGs, 36 DELs) regulatory networks. The circRNAs sus-TRPM7_0005, sus-MTUS1_0004, the lncRNAs SMSTRG.4269.1, and MSTRG.7983.2 regulate the expression of six lipid metabolism-related target genes, including PLCB1, BAD, and GADD45G, through the binding sites of 2-4068, miR-7134-3p, and miR-190a. For instance, MSTRG.4269.1 regulates its targets PLCB1 and BAD via miRNA 2_4068. Meanwhile, sus-TRPM7_0005 controls its target LRP5 through ssc-miR-7134-3P. These findings indicate molecular regulatory networks that could potentially be applied for the marker-assisted selection of IMF to enhance pork quality.

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.

Integrative analysis of plasma and substantia nigra in Parkinson's disease: unraveling biomarkers and insights from the lncRNA-miRNA-mRNA ceRNA network.

Parkinson's disease (PD) is a rapidly growing neurological disorder characterized by diverse movement symptoms. However, the underlying causes have not been clearly identified, and accurate diagnosis is challenging. This study aimed to identify potential biomarkers suitable for PD diagnosis and present an integrative perspective on the disease. We screened the GSE7621, GSE8397-GPL96, GSE8397-GPL97, GSE20163, and GSE20164 datasets in the NCBI GEO database to identify differentially expressed (DE) mRNAs in the substantia nigra (SN). We also screened the GSE160299 dataset from the NCBI GEO database to identify DE lncRNAs and miRNAs in plasma. We then constructed 2 lncRNA-miRNA-mRNA competing endogenous RNA (ceRNA) regulatory networks based on the ceRNA hypothesis. To understand the biological function, we performed Kyoto Encyclopedia of Genes and Genomes pathway and Gene Ontology analyses for each ceRNA network. The receiver operating characteristic analyses (ROC) was used to assess ceRNA results. We identified 7 upregulated and 29 downregulated mRNAs as common DE mRNAs in the 5 SN datasets. In the blood dataset, we identified 31 DE miRNAs (9 upregulated and 22 downregulated) and 332 DE lncRNAs (69 upregulated and 263 downregulated). Based on the determined interactions, 5 genes (P2RX7, HSPA1, SLCO4A1, RAD52, and SIRT4) appeared to be upregulated as a result of 10 lncRNAs sponging 4 miRNAs (miR-411, miR-1193, miR-301b, and miR-514a-2/3). Competing with 9 genes (ANK1, CBLN1, RGS4, SLC6A3, SYNGR3, VSNL1, DDC, KCNJ6, and SV2C) for miR-671, a total of 26 lncRNAs seemed to function as ceRNAs, influencing genes to be downregulated. In this study, we successfully constructed 2 novel ceRNA regulatory networks in patients with PD, including 36 lncRNAs, 5 miRNAs, and 14 mRNAs. Our results suggest that these plasma lncRNAs are involved in the pathogenesis of PD by sponging miRNAs and regulating gene expression in the SN of the brain. We propose that the upregulated and downregulated lncRNA-mediated ceRNA networks represent mechanisms of neuroinflammation and dopamine neurotransmission, respectively. Our ceRNA network, which was associated with PD, suggests the potential use of DE miRNAs and lncRNAs as body fluid diagnostic biomarkers. These findings provide an integrated view of the mechanisms underlying gene regulation and interactions in PD.

Identification of different myofiber types in pigs muscles and construction of regulatory networks

BackgroundSkeletal muscle is composed of muscle fibers with different physiological characteristics, which plays an important role in regulating skeletal muscle metabolism, movement and body homeostasis. The type of skeletal muscle fiber directly affects meat quality. However, the transcriptome and gene interactions between different types of muscle fibers are not well understood.ResultsIn this paper, we selected 180-days-old Large White pigs and found that longissimus dorsi (LD) muscle was dominated by fast-fermenting myofibrils and soleus (SOL) muscle was dominated by slow-oxidizing myofibrils by frozen sections and related mRNA and protein assays. Here, we selected LD muscle and SOL muscle for transcriptomic sequencing, and identified 312 differentially expressed mRNA (DEmRs), 30 differentially expressed miRNA (DEmiRs), 183 differentially expressed lncRNA (DElRs), and 3417 differentially expressed circRNA (DEcRs). The ceRNA network included ssc-miR-378, ssc-miR-378b-3p, ssc-miR-24-3p, XR_308817, XR_308823, SMIM8, MAVS and FOS as multiple core nodes that play important roles in muscle development. Moreover, we found that different members of the miR-10 family expressed differently in oxidized and glycolytic muscle fibers, among which miR-10a-5p was highly expressed in glycolytic muscle fibers (LD) and could target MYBPH gene mRNA. Therefore, we speculate that miR-10a-5p may be involved in the transformation of muscle fiber types by targeting the MYHBP gene. In addition, PPI analysis of differentially expressed mRNA genes showed that ACTC1, ACTG2 and ACTN2 gene had the highest node degree, suggesting that this gene may play a key role in the regulatory network of muscle fiber type determination.ConclusionsWe can conclude that these genes play a key role in regulating muscle fiber type transformation. Our study provides transcriptomic profiles and ceRNA interaction networks for different muscle fiber types in pigs, providing reference for the transformation of pig muscle fiber types and the improvement of meat quality.

Key mRNAs and lncRNAs of pituitary that affect the reproduction of FecB + + small tail han sheep

BackgroundThe pituitary directly regulates the reproductive process through follicle-stimulating hormone (FSH) and luteinizing hormone (LH). Transcriptomic research on the pituitaries of ewes with different FecB (fecundity Booroola) genotypes has shown that some key genes and lncRNAs play an important role in pituitary function and sheep fecundity. Our previous study found that ewes with FecB + + genotypes (without FecB mutation) still had individuals with more than one offspring per birth. It is hoped to analyze this phenomenon from the perspective of the pituitary transcriptome.ResultsThe 12 Small Tail Han Sheep were equally divided into polytocous sheep in the follicular phase (PF), polytocous sheep in the luteal phase (PL), monotocous sheep in the follicular phase (MF), and monotocous sheep in the luteal phase (ML). Pituitary tissues were collected after estrus synchronous treatment for transcriptomic analysis. A total of 384 differentially expressed genes (DEGs) (182 in PF vs. MF and 202 in PL vs. ML) and 844 differentially expressed lncRNAs (DELs) (427 in PF vs. MF and 417 in PL vs. ML) were obtained from the polytocous-monotocous comparison groups in the two phases. Functional enrichment analysis showed that the DEGs in the two phases were enriched in signaling pathways known to play an important role in sheep fecundity, such as calcium ion binding and cAMP signaling pathways. A total of 1322 target relationship pairs (551 pairs in PF vs. MF and 771 pairs in PL vs. ML) were obtained for the target genes prediction of DELs, of which 29 DEL-DEG target relationship pairs (nine pairs in PF vs. MF and twenty pairs in PL vs. ML). In addition, the competing endogenous RNA (ceRNA) networks were constructed to explore the regulatory relationships of DEGs, and some important regulatory relationship pairs were obtained.ConclusionAccording to the analysis results, we hypothesized that the pituitary first receives steroid hormone signals from the ovary and uterus and that VAV3 (Vav Guanine Nucleotide Exchange Factor 3), GABRG1 (Gamma-Aminobutyric Acid A Receptor, Gamma 1), and FNDC1 (Fibronectin Type III Domain Containing 1) played an important role in this process. Subsequently, the reproductive process was regulated by gonadotropins, and IGFBP1 (Insulin-like Growth Factor Binding Protein 1) was directly involved in this process, ultimately affecting litter size. In addition, TGIF1 (Transforming Growth Factor-Beta-Induced Factor 1) and TMEFF2 (Transmembrane Protein With EGF Like And Two Follistatin Like Domains 2) compensated for the effect of the FecB mutation and function by acting on TGF-β/SMAD signaling pathway, an important pathway for sheep reproduction. These results provided a reference for understanding the mechanism of multiple births in Small Tail Han Sheep without FecB mutation.

LncRNA-mRNA Co-Expression and Regulation Analysis in Lung Fibroblasts from Idiopathic Pulmonary Fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a progressive lung disease marked by abnormal accumulation of extracellular matrix (ECM) due to dysregulated expression of various RNAs in pulmonary fibroblasts. This study utilized RNA-seq data meta-analysis to explore the regulatory network of hub long non-coding RNAs (lncRNAs) and messenger RNAs (mRNAs) in IPF fibroblasts. The meta-analysis unveiled 584 differentially expressed mRNAs (DEmRNA) and 75 differentially expressed lncRNAs (DElncRNA) in lung fibroblasts from IPF. Among these, BCL6, EFNB1, EPHB2, FOXO1, FOXO3, GNAI1, IRF4, PIK3R1, and RXRA were identified as hub mRNAs, while AC008708.1, AC091806.1, AL442071.1, FAM111A-DT, and LINC01989 were designated as hub lncRNAs. Functional characterization revealed involvement in TGF-β, PI3K, FOXO, and MAPK signaling pathways. Additionally, this study identified regulatory interactions between sequences of hub mRNAs and lncRNAs. In summary, the findings suggest that AC008708.1, AC091806.1, FAM111A-DT, LINC01989, and AL442071.1 lncRNAs can regulate BCL6, EFNB1, EPHB2, FOXO1, FOXO3, GNAI1, IRF4, PIK3R1, and RXRA mRNAs in fibroblasts bearing IPF and contribute to fibrosis by modulating crucial signaling pathways such as FoxO signaling, chemical carcinogenesis, longevity regulatory pathways, non-small cell lung cancer, and AMPK signaling pathways.

Whole-transcriptome profiling reveals potential biomarkers for the reversal of thymic epithelial cell senescence by umbilical cord mesenchymal stem cells.

Reduced numbers and dysfunction of thymic epithelial cells (TECs) are important factors of thymic degeneration. Previous studies have found that umbilical cord mesenchymal stem cells (UCMSCs) reverse the structure and function of the senescent thymus in vivo. However, the transcriptomic regulation mechanism is unclear. TECs were cultured with H2O2 for 72 hours to induce senescence. UCMSCs were cocultured with senescent TECs for 48 hours to detect SA-β-gal, P16 and Ki67. The cocultured TECs were collected for lncRNA, mRNA and miRNA sequencing to establish a competitive endogenous regulatory network (ceRNA). And RT-qPCR, immunofluorescence staining, and western blot were used to identified key genes. Our results showed that H2O2 induced TEC aging and that UCMSCs reversed these changes. Compared with those in aged TECs, 2260 DE mRNAs, 1033 DE lncRNAs and 67 DE miRNAs were differentially expressed, and these changes were reversed by coculturing the cells with UCMSCs. Differential mRNA enrichment analysis of ceRNA regulation revealed that the PI3K-AKT pathway was a significant signaling pathway. UCMSC coculture upregulated VEGFA, which is the upstream factor of the PI3K-AKT signaling pathway, and the expression of the key proteins PI3K and AKT. Thus, the expression of the cell cycle suppressor P27, which is downstream of the PI3K-AKT signaling pathway, was downregulated, while the expression of the cell cycle regulators CDK2 and CCNE was upregulated. UCMSC coculture upregulated the expression of VEGFA, activated the PI3K-AKT signaling pathway, increased the expression of CDK2 and CCNE, decreased the expression of P27, and promoted the proliferation of TECs.

Whole-Transcriptome Sequencing of Knee Joint Cartilage from Kashin-Beck Disease and Osteoarthritis Patients.

The aim of this study was to provide a comprehensive understanding of similarities and differences in mRNAs, lncRNAs, and circRNAs within cartilage for Kashin-Beck disease (KBD) compared to osteoarthritis (OA). We conducted a comparison of the expression profiles of mRNAs, lncRNAs, and circRNAs via whole-transcriptome sequencing in eight KBD and ten OA individuals. To facilitate functional annotation-enriched analysis for differentially expressed (DE) genes, DE lncRNAs, and DE circRNAs, we employed bioinformatic analysis utilizing Gene Ontology (GO) and KEGG. Additionally, using quantitative reverse transcriptase polymerase chain reaction (qRT-PCR), we validated the expression levels of four cartilage-related genes in chondrocytes. We identified a total of 43 DE mRNAs, 1451 DE lncRNAs, and 305 DE circRNAs in KBD cartilage tissue compared to OA (q value < 0.05; |log2FC| > 1). We also performed competing endogenous RNA network analysis, which identified a total of 65 lncRNA-mRNA interactions and 4714 miRNA-circRNA interactions. In particular, we observed that circRNA12218 had binding sites for three miRNAs targeting ACAN, while circRNA12487 had binding sites for seven miRNAs targeting COL2A1. Our results add a novel set of genes and non-coding RNAs that could potentially serve as candidate diagnostic biomarkers or therapeutic targets for KBD patients.

Discovery of long non-coding RNAs in naïve CD4+ T cells in response to initiating antiretroviral therapy at acute or chronic phase of HIV-1 infection

Long non-coding RNAs (lncRNAs) have gained prominence due to their involvement in various cellular processes, but their specific roles remain elusive. Dysregulation of lncRNAs has been implicated in the pathogenesis of several diseases. In this study, we aimed to shed light on the role of lncRNAs in individuals infected with human immunodeficiency virus type 1 (HIV-1) by examining their changes in the expression patterns related to the initiation of antiretroviral therapy (ART) during acute or chronic phases of infection, compared to healthy controls. We found 316 differentially expressed (DE) lncRNAs in patients receiving long-term ART, shedding light on their potential roles. We also observed interactions between these DE lncRNAs and specific microRNAs (miRNAs). Some of these miRNAs, such as hsa-miR-574-5p, hsa-miR-765, hsa-miR-6165, hsa-miR-1207-5p, and hsa-miR-378i, are associated with cancer progression or suppression, while others, including hsa-miR-328-5p, hsa-miR-4753-3p, and MiR-664, play roles in immune system regulation. Furthermore, our study revealed substantial enrichment in distinct Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, such as axon guidance, PI3K-Akt signaling, and MAPK signaling pathways. Although our results indicate possible molecular processes impacted by the discovered lncRNAs, we cannot explicitly establish causality or specific connections between lncRNAs and genes in these pathways, fostering more specific studies. Furthermore, Gene Ontology (GO) analysis highlighted terms such as cytoskeletal protein binding, ion channel function, synaptic processes, neuron projection, and the somatodendritic compartment, underscoring the relevance of lncRNAs in these cellular components within the context of HIV-1 infection and ART treatment. In conclusion, our study emphasizes the need for further exploration of lncRNAs as potential biomarkers and therapeutic targets in HIV-1-infected patients, with a particular focus on CD4+ T cells. Understanding the functions of lncRNAs in these contexts may pave the way for novel treatment strategies and improved patient outcomes, aligning with the broader goals of our research.