Differentially Expressed mRNAs Research Articles

Whole transcriptome sequencing of testis and epididymis reveals genes associated with sperm development in roosters

BackgroundChickens play a crucial role as the primary global source of eggs and poultry, and the quality of rooster semen significantly impacts poultry reproductive efficiency. Therefore, it is imperative to comprehend the regulatory mechanisms underlying sperm development.ResultsIn this study, we established transcriptome profiles of lncRNAs, miRNAs, and mRNAs in 3 testis tissues and 3 epididymis tissues from “Jing Hong No.1” roosters at 24, 35, and 64 weeks of age. Using the data, we conducted whole transcriptome analysis and constructed a ceRNA network. We detected 10 differentially expressed mRNAs (DEmRNAs), 33 differentially expressed lncRNAs (DElncRNAs), and 10 differentially expressed miRNAs (DEmiRNAs) in the testis, as well as 149 DEmRNAs, 12 DElncRNAs, and 10 DEmiRNAs in the epididymis. These genes were found to be involved in cell differentiation and development, as well as various signaling pathways such as GnRH, MAPK, TGF-β, mTOR, VEGF, and calcium ion pathways. Subsequently, we constructed two competing endogenous RNA (ceRNA) networks comprising DEmRNAs, DElncRNAs, and DEmiRNAs. Furthermore, we identified four crucial lncRNA-mRNA-miRNA interactions that govern specific biological processes in the chicken reproductive system: MSTRG.2423.1-gga-miR-1563-PPP3CA and MSTRG.10064.2-gga-miR-32-5p-GPR12 regulating sperm motility in the testis; MSTRG.152556.1-gga-miR-9-3p-GREM1/THYN1 governing immunomodulation in the epididymis; and MSTRG.124708.1-gga-miR-375-NDUFB9/YBX1 controlling epididymal sperm maturation and motility.ConclusionsWhole transcriptome sequencing of chicken testis and epididymis screened several key genes and ceRNA regulatory networks, which may be involved in the regulation of epididymal immunity, spermatogenesis and sperm viability through the pathways of MAPK, TGF-β, mTOR, and calcium ion. These findings contribute to our comprehensive understanding of the intricate molecular processes underlying rooster spermatogenesis, maturation and motility.

KLF15 Transcription Activated Tachykinin Precursor 1 to Induce Perineural Invasion in Head and Neck Squamous Cell Carcinoma.

In this study, we aimed to identify novel biomarkers related to Peripheral Neural Invasion (PNI) in head and neck squamous cell carcinoma (HNSCC). The PNI-related differentially expressed mRNAs (DE-mRNAs) in HNSCC were identified to construct a PNI-related risk score model. The expression level and ROC curve for Tachykinin Precursor 1 (TAC1) were calculated. Additionally, two kinds of in vitro models of PNI were established for investigation, including the Matrigel-PNI model and the Transwell-PNI model. Furthermore, the transcription factor of the TAC1 was predicted and verified by qRTPCR. A total of 139 DE-mRNAs were identified in PNI positive and negative groups of HNSCC patients. The risk-score marker model incorporating 20 PNI-related DE-mRNAs was established. The TAC1 was identified as a potential highly expressed PNI marker, which exhibited good performance in predicting PNI events. Patients with higher TAC1 expressions demonstrated significantly shorter survival rates compared to those with lower TAC1 expressions in HNSCC. Besides, the knockdown of TAC1 significantly repressed neural invasion in HNSCC cells in vitro, according to the Matrigel-PNI model and Transwell-PNI model. Furthermore, KLF15 was predicted and verified as a transcription activator of TAC1 in HNSCC. This study highlights that the activation of KLF15 transcription of TAC1 promotes PNI in HNSCC cells, which provides guidance regarding the molecular diagnosis of PNI in HNSCC cells.

Bioinformatics analysis of immune infiltration in human diabetic retinopathy and identification of immune-related hub genes and their ceRNA networks

Diabetic retinopathy (DR) is the most common microvascular complication in diabetic patients, and recent studies have shown that immune regulatory mechanisms are closely associated with retinal damage in DR. Therefore, this study focused on exploring immune cells and immune-related genes (IRGs) in DR and gaining insight into the ceRNA mechanisms by which IRGs regulate DR progression. Four datasets from human DR model retinal tissues were obtained from the Gene Expression Omnibus (GEO) database. R software was first used to identify differentially expressed mRNAs (DE-mRNAs) in the dataset GSE160306-mRNAs, then the distribution of immune cells in the gene matrix was analyzed by xCell and ImmuCellAI, ImmPort and InnateDB database were used to obtain immune-related hub genes (IRHGs) in the DR, and finally the STRING online tool and Cytoscape to construct the immune-related ceRNA network. The datasets GSE102485, GSE160308 and GSE160306-lncRNAs were used to validate the results of the ceRNA network further. The results of immune cell infiltration analysis showed that macrophages are important immune cells in DR; immune-related gene screening showed that FCGR2B is an IRHG in DR, and 2 immune-related ceRNA networks of IRHG were obtained: DDN-AS1/miR-10a-5p/FCGR2B and LINC01515/miR-10a-5p/FCGR2B. Our study suggests that infiltration of immune cells, especially the immune role of macrophages, is an important component of DR progression; the immune-related hub gene FCGR2B and its ceRNA network may be a key regulatory network for DR progression. The discovery of key immune cells, IRHG and ceRNA networks in this study may provide new prospects for early intervention and targeted treatment of DR.

Deoxynivalenol-Induced Spleen Toxicity in Mice: Inflammation, Endoplasmic Reticulum Stress, Macrophage Polarization, and the Dysregulation of LncRNA Expression.

The spleen is a primary target of deoxynivalenol (DON) toxicity, but its underlying molecular mechanisms remain unclear. This study investigates the effects of DON on inflammation, splenic macrophage polarization, endoplasmic reticulum (ER) stress, and transcriptome changes (mRNA and lncRNAs) in mouse spleen. We found that DON exposure at doses of 2.5 or 5 mg/kg BW significantly induced inflammation and polarized splenic macrophages towards the M1 phenotype. Additionally, DON activated PERK-eIF2α-ATF4-mediated ER stress and upregulated apoptosis-related proteins (caspase-12, caspase-3). The ER stress inhibitor, 4-Phenylbutyric acid, significantly alleviated DON-induced ER stress, apoptosis, and the M1 polarization of splenic macrophages. Transcriptome analysis identified 1968 differentially expressed (DE) lncRNAs and 2664 DE mRNAs in mouse spleen following DON exposure. Functional enrichment analysis indicated that the upregulated genes were involved in pathways associated with immunity, including Th17 cell differentiation, TNF signaling, and IL-17 signaling, while downregulated mRNAs were linked to cell survival and growth pathways. Furthermore, 370 DE lncRNAs were predicted to target 255 DE target genes associated with immune processes, including the innate immune response, interferon-beta response, cytokine production regulation, leukocyte apoptosis, and NF-κB signaling genes. This study provides new insights into the mechanisms underlying DON toxicity and its effects on the immune system.

Whole-transcriptome analysis reveals the effect of retinoic acid on small intestinal mucosal injury in cage-stressed young laying ducks

Retinoic acid (RA) is an active derivative of vitamin A and is involved in a variety of physiological processes, including cell growth, antioxidant, and inflammation. However, the role of RA in intestinal oxidative stress injury in caged-stressed laying ducks is unknown. In this study, we analyzed the effect and underlying mechanism of RA supplementation on intestinal damage in cage-stressed young laying ducks. One hundred and sixty laying ducks were divided into 5 treatment groups, including a control group (CR) and 4 treatment groups exposed to different RA concentrations (2,500, 5,000, 7,500 and 10,000 IU/kg, TG1 to TG4). The experimental period comprised a 7-d prefeeding period and a 10-d experimental feeding period, for a total of 17 d. Phenotypic analysis revealed that compared with the control group, RA addition increased the intestinal villus height and the villus-to-crypt ratio; decreased the crypt depth (P < 0.01); decreased the serum diamine oxidase and D-lactate concentrations (P < 0.05); increased the serum antioxidant capacity and intestinal antioxidant gene expression levels (P < 0.05); and increased the expression levels of tight junction-related genes, with the greatest effect observed in TG2 group. Our further whole-transcriptome analysis of duodenum tissues from CR and TG2 ducks revealed 706 differentially expressed mRNAs (DEmRNAs), 357 differentially expressed lncRNAs (DElncRNAs), 14 differentially expressed circRNAs (DEcircRNAs), and 4 differentially expressed miRNAs (DEmiRNAs). These DEGs are involved in calcium signaling, NOD-like receptor signaling, pyruvate metabolism, Jak-STAT signaling, Wnt signaling, riboflavin metabolism, and the adherens junction and tight junction pathways. The results of omics and marker gene expression analysis suggested that RA treatment may play a role in endoplasmic reticulum stress (ERS) and apoptosis. In conclusion, the addition of RA to the diet improved intestinal injury by improving the redox homeostasis of intestinal cells associated with ERS, enhancing the intestinal tight junction structure and alleviating the apoptosis of intestinal epithelial cells; moreover, 5,000 IU/kg RA was determined to be the most appropriate concentration for supplementation.

Host miRNA and mRNA profiles during in DEF and duck after DHAV-1 infection

DHAV-1 is a highly infectious pathogen that can cause acute hepatitis in ducklings. MicroRNA (miRNA) plays an essential regulatory role in virus response. We characterized and compared miRNA and mRNA expression profiles in duck embryonic fibroblasts (DEF) and the liver of ducklings infected with DHAV-1. DHAV-1 infected DEF was divided into infection group (D group) and blank group (M group), and DHAV-1 infected duckling group was divided into infection group (H group) and blank group (N group). D vs. M have 130 differentially expressed (DE) miRNA (DEM) and 2204 differentially expressed (DE) mRNA (DEG), H vs. N have 72 DEM and 1976 DEG. By the intersection of D vs. M and H vs. N comparisons, 15 upregulated DEM, 5 downregulated DEM, 340 upregulated DEG and 50 downregulated DEG were found with both in vivo and in vitro DHAV-1 infection. In particular, we identified the same DE miRNA target genes and functional annotations of DE mRNA. We enriched with multiple gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, which may have important roles in viral virulence, host immunity, and metabolism. We selected miR-155, which is co-upregulated, and found that miR-155 targets SOCS1 to inhibit DHVA-1 replication.

UID-Dual Transcriptome Sequencing Analysis of the Molecular Interactions between Streptococcus agalactiae ATCC 27956 and Mammary Epithelial Cells.

Streptococcus agalactiae ATCC 27956 is a highly contagious Gram-positive bacterium that causes mastitis, has a high infectivity for mammary epithelial cells, and becomes challenging to treat. However, the molecular interactions between it and mammary epithelial cells remain poorly understood. This study analyzed differential gene expression in mammary epithelial cells with varying levels of S. agalactiae infection using UID-Dual transcriptome sequencing and bioinformatics tools. This study identified 211 differentially expressed mRNAs (DEmRNAs) and 452 differentially expressed lncRNAs (DElncRNAs) in host cells, primarily enriched in anti-inflammatory responses, immune responses, and cancer-related processes. Additionally, 854 pathogen differentially expressed mRNAs (pDEmRNAs) were identified, mainly enriched in protein metabolism, gene expression, and biosynthesis processes. Mammary epithelial cells activate pathways, such as the ERK1/2 pathway, to produce reactive oxygen species (ROS) to eliminate bacteria. The bacteria disrupt the host's innate immune mechanisms by interfering with the alternative splicing processes of mammary epithelial cells. Specifically, the bacterial genes of tsf, prfB, and infC can interfere with lncRNAs targeting RUNX1 and BCL2L11 in mammary epithelial cells, affecting the alternative splicing of target genes and altering normal molecular regulation.

Exploration and validation of ceRNA regulatory networks in colorectal cancer based on associations whole transcriptome sequencing

Colorectal cancer (CRC) is a wide-spread gastrointestinal cancer that is associated with augmented morbidity and mortality, and we do not yet have a deep understanding of its epidemiology and carcinogenicity. The transcriptome can reveal the complexity and heterogeneity of tumors and uncover new biomarkers or treatment options. In this study, we identified messenger RNAs (mRNAs), long non-coding RNAs (lncRNAs), round RNAs (circRNAs), and microRNAs (miRNAs) using whole-transcriptome sequencing and generated competing endogenous RNA (ceRNA) modulatory axes. We conducted whole transcriptome sequencing on 10 CRC and para-cancer (CRCP) samples and discovered 2465 differentially expressed (DE) mRNAs (DEmRNAs), 77 DE miRNAs (DEmiRNAs). 2852 DE lncRNAs (DElncRNAs) and 1477 DE circRNAs (DEcircRNAs). In addition, utilizing co-DE analysis, we generated the ceRNA axis. Subsequently, we employed the ceRNA axis to identify essential genes and corresponding associations with lncRNAs, circRNAs, and miRNAs in CRC. ceRNA regulatory network including mRNA-miRNA-lncRNA and mRNA-miRNA-circRNA. These modulatory axes potentially modulate the positive regulation of smooth muscle contraction, melanosome, plasma membrane, integral plasma membrane component and so on. Finally, the results of RNA sequencing (RNA-SEQ) were combined with the TCGA and GEO databases, and the DEGs strongly correlated with the TCGA-COAD overall survival (OS) as estimated by univariate cox and logarithmic rank analyses were cross-analyzed, and the co-upregulated DEGs were screened. Among the many DEs, KPNA2 was chosen for additional analysis. Using invitro experimentations, western blot, CCK8, EdU and other experiments were performed to verify the results. We found siRNA-based KPNA2 depletion reduces bladder cancer cells’ viability, migratory, and proliferative activities, which showed that the DEmRNA profiles were comparable to the sequencing information, confirming that the sequencing data were very reliable. These evidences highlight the ceRNA regulatory mechanisms in CRC and will aid future research into the molecular mechanisms behind colorectal cancer prevention and treatment.

Whole-transcriptome sequencing analysis to identify key circRNAs, miRNAs, and mRNAs in the development of yak testes

BackgroundThe Testis is an important reproductive organ in male mammals and the site for spermatogenesis, androgen synthesis, and secretion. Non-coding RNAs (ncRNAs) play an important regulatory role in various biological processes. However, the regulatory role of ncRNAs in the development of yak testes and spermatogenesis remains largely unclear.ResultIn this study, we compared the expression profiles of circular RNAs (circRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs) in yak testicular tissue samples collected at 6 months (Y6M), 18 months (Y18M), and 4 years (Y4Y). Using RNA sequencing (RNA-Seq), we observed a significant difference in the expression patterns of ncRNAs in the samples collected at different testicular development stages. Twenty-two differentially expressed (DE) circRNAs, 69 DE miRNAs, and 64 DE mRNAs were detected in Y6M, Y18M, and Y4Y testicular samples, respectively. The results of gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that the source genes of DE circRNAs, predicted target genes of DE miRNAs, and DE mRNAs were specifically associated with signaling pathways and GO terms that were related to sperm synthesis, sperm vitality, and testicular development, such as cell cycle, Wnt signaling pathway, MAPK signaling pathway, GnRH signaling pathway, and spermatogenesis. The analysis of the circRNA-miRNA-mRNA network revealed that some DE ncRNAs, including miR-574, miR-449a, CDC42, and CYP11A1, among others, may be involved in testicular spermatogenesis. Concurrently, various circRNA-miRNA interaction pairs were observed.ConclusionOur findings provide a database of circRNAs, miRNAs, and mRNAs expression profiles in testicular tissue of yaks at different developmental stages and a detailed understanding of the regulatory network of ncRNAs in yak testicular development and provide data that can help elucidate the molecular mechanisms underlying yak testicular development.

Analysis of the Long Non-Coding and Messenger RNA Expression Profiles in the Skin Tissue of Super Merino and Small-Tailed Han Sheep.

Wool quality and yield are two important economic livestock traits. However, there are relatively few molecular studies on lncRNA for improving sheep wool, so these require further exploration. In this study, we examined skin tissue from the upper scapula of Super Merino (SM) and Small-Tailed Han (STH) sheep during the growing period. The apparent difference was verified via histological examination. High-throughput RNA sequencing identified differentially expressed (DE) long non-coding (lncRNAs) and messenger RNAs (mRNAs). The target gene of DE lncRNA and DE genes were enrichment analyzed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). A Reverse Transcription quantitative Polymerase Chain Reaction (RT-qPCR) was used to verify randomly selected DE lncRNAs and mRNAs. Finally, the DE, RAC2, WNT11, and FZD2 genes, which were enriched in the Wnt signaling pathway, were detected via immunohistochemistry. The results showed that a total of 20,888 lncRNAs and 31,579 mRNAs were identified in the skin tissues of the two sheep species. Among these, 56 lncRNAs and 616 mRNAs were differentially expressed. Through qRT-PCR, the trends in the randomly selected DE genes' expression were confirmed to be aligned with the RNA-seq results. GO and KEGG enrichment analysis showed that DE lncRNA target genes were enriched in GO terms as represented by epidermal and skin development and keratin filature and in KEGG terms as represented by PI3K-Akt, Ras, MAPK, and Wnt signaling pathways, which were related to hair follicle growth and development. Finally, immunohistochemistry staining results indicated that RAC2, WNT11, and FZD2 were expressed in dermal papilla (DP). The lncRNAs MSTRG.9225.1 and MSTRG.98769.1 may indirectly participate in the regulation of hair follicle growth, development, and fiber traits by regulating their respective target genes, LOC114113396(KRTAP15-1), FGF1, and IGF1. In addition, MSTRG.84658.1 may regulate the Wnt signaling pathway involved in the development of sheep hair follicles by targeting RAC2. This study provides a theoretical reference for improving sheep breeding in the future and lays a foundation for further research on the effects of MSTRG.84658.1 and the target gene RAC2 on dermal papilla cells (DPC).

Comprehensive analysis of lncRNA-mediated ceRNA network in renal cell carcinoma based on GEO database.

Renal cell carcinoma (RCC) ranks among the leading causes of cancer-related mortality. Despite extensive research, the precise etiology and progression of RCC remain incompletely elucidated. Long noncoding RNA (lncRNA) has been identified as competitive endogenous RNA (ceRNA) capable of binding to microRNA (miRNA) sites, thereby modulating the expression of messenger RNAs (mRNA) and target genes. This regulatory network is known to exert a pivotal influence on cancer initiation and progression. However, the specific role and functional significance of the lncRNA-miRNA-mRNA ceRNA network in RCC remain poorly understood. The RCC transcriptome data was obtained from the gene expression omnibus database. The identification of differentially expressed long noncoding RNAs (DElncRNAs), differentially expressed miRNAs, and differentially expressed mRNAs (DEmRNAs) between RCC and corresponding paracancer tissues was performed using the "Limma" package in R 4.3.1 software. We employed a weighted gene co-expression network analysis to identify the key DElncRNAs that are most relevant to RCC. Subsequently, we utilized the encyclopedia of RNA interactomes database to predict the interactions between these DElncRNAs and miRNAs, and the miRDB database to predict the interactions between miRNAs and mRNAs. Therefore, key DElncRNAs were obtained to verify the expression of their related genes in the The Cancer Genome Atlas database and to analyze the prognosis. The construction of RCC-specific lncRNA-miRNA-mRNA ceRNA network was carried out using Cytoscape 3.7.0. A total of 286 DElncRNAs, 56 differentially expressed miRNAs, and 2065 DEmRNAs were identified in RCC. Seven key DElncRNAs (GAS6 antisense RNA 1, myocardial infarction associated transcript, long intergenic nonprotein coding RNA 921, MMP25 antisense RNA 1, Chromosome 22 Open Reading Frame 34, MIR34A host gene, MIR4435-2 host gene) were identified using weighted gene co-expression network analysis and encyclopedia of RNA interactomes databases. Subsequently, a network diagram comprising 217 nodes and 463 edges was constructed based on these key DElncRNAs. The functional analysis of DEmRNAs in the ceRNA network was conducted using Kyoto Encyclopedia of Genes and Genomes and gene ontology. We constructed RCC-specific ceRNA networks and identified the crucial lncRNAs associated with RCC using bioinformatics analysis, which will help us further understand the pathogenesis of this disease.

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

BackgroundThe 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.ResultsHere, 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.ConclusionsThese 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 functional characterization of differentially expressed circRNAs in high glucose treated endothelial cells: Construction of circRNA-miRNA-mRNA network

BackgroundEndothelial dysfunction is a complication of diabetes mellitus (DM), characterized by impaired endothelial function in both microvessels and macrovessels, closely linked to atherosclerosis (AS). Endothelial dysfunction, characterized by impaired endothelial cell (EC) function, is a pivotal factor in AS and DM. Circular RNAs (circRNAs) are endogenous non-coding RNAs that can act as competing endogenous RNAs (ceRNAs) and regulate gene expression. However, the role of circRNAs in ECs dysfunction and AS under high glucose (HG) condition remains elusive. MethodsWe performed high-throughput sequencing to identify differentially expressed (DE) circRNAs in human umbilical vein endothelial cells (HUVEC) exposed to HG, one risk factors of endothelial dysfunction and AS. We then validated eight candidate circRNAs by qRT-PCR and functional analysis, directing our attention to hsa_circ_0122319. Moreover, microarray analysis identified the differential expression profiles of miRNAs and mRNAs regulated by hsa_circ_0122319. Subsequently, the construction of the ceRNAs network employed bioinformatic analysis and Cytoscape software. Furthermore, the role of the PI3K-Akt signaling pathway in regulating ceRNAs was evaluated. ResultsWe detected 917 DE circRNAs in HG treated HUVEC. The parental genes of these circRNAs were enriched in cell cycle, cellular senescence and endocytosis related pathways. The differential expression of hsa_circ_0122319 was confirmed to be most obvious at the cellular level and in clinical samples by qPCR experiments. After overexpression of hsa_circ_0122319, 49 DE miRNAs and 459 DE mRNAs were identified using microarray analysis. Subsequently, a ceRNAs network was constructed, comprising hsa_circ_0122319, 8 miRNAs, and 41 mRNAs. ConclusionIn summary, our study delves into the role of circRNAs in endothelial dysfunction associated with DM and AS. Through high-throughput sequencing and validation, we identified hsa_circ_0122319 as a pivotal regulator of ECs function under HG conditions. It also showed that hsa_circ_0123319 has the potential to serve as a biomarker for DM and its vascular complications, and provides new evidence for future exploration of the intricate molecular mechanisms of endothelial dysfunction in the progression of DM and AS.

Whole transcriptome sequencing revealed the gene regulatory network of hypoxic response in yak Sertoli cells

Yaks live in the Qinghai-Tibet Plateau for a long time where oxygen is scarce, but can ensure the smooth development of testis and spermatogenesis. The key lies in the functional regulation of the Sertoli cells under hypoxia. In this study, we sequenced yak Sertoli cells cultured in normal oxygen concentration (Normoxia) and treated with low oxygen concentration (Hypoxia) by whole transcriptomics, and screened out 194 differentially expressed mRNAs (DEmRNAs), 934 differentially expressed LncRNAs (DELncRNAs) and 129 differentially expressed miRNAs (DEmiRNAs). GO and KEGG analysis showed that these differential genes were mainly concentrated in PI3K-AKT, MAPK, RAS, and other signaling pathways, and were associated with glucose metabolism, tight junction, steroid hormone synthesis, cell fusion, and immunity of yak Sertoli cells. We constructed the gene interaction network of yak Sertoli cells in hypoxia and screened out the relationship pairs related to glucose metabolism and tight junction. The results suggested that the changes in energy metabolism, tight junction, and immune regulation of yak Sertoli cells under hypoxia might provide favorable conditions for spermatogenesis. This study provides data for further study on the role of non-coding RNA in testis development and spermatogenesis of yak.

Core competing endogenous RNA network based on mRNA and non-coding RNA expression profiles in chicken fatty liver.

Fatty liver disease is a common metabolic disease in chickens. This disease can lead to a decrease in egg production and increase the risk of death in chickens. Long non-coding RNAs (lncRNAs) are involved in fatty liver formation by directly targeting genes or regulating gene expression by competitively binding microRNAs. However, a large proportion of competing endogenous RNA (ceRNA) networks in fatty liver diseases are still unclear. The total of 300 Jingxing-Huang chickens were used for fatty liver model construction. Then, differentially expressed (DE) genes (DEGs) identified through whole-transcriptome sequencing from four chickens with fatty liver and four chickens without fatty liver were chosen from the F1 generation. A total of 953 DEGs were identified between the fatty liver group and the control group, including 26 DE micro (mi)RNAs and 56 DE lncRNAs. Differential expression heatmaps and volcano plots were obtained after clustering expression analysis. Gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed that these DEGs were involved in many biological processes and signaling pathways related to fatty acid metabolism and lipid synthesis. Furthermore, cytoscape was used to construct a ceRNA network of the DE miRNAs, DE mRNAs, and DE lncRNAs. Eleven DE lncRNAs, seven DE miRNAs, and 13 DE mRNAs were found to be associated with the pathogenesis of fatty liver disease. An lncRNA-miRNA-mRNA ceRNA network was constructed to elucidate the mechanisms of fatty liver diseases, and the ENSGALT00000079786-miR-140/miR-143/miR-1a/miR-22/miR-375 network was identified. These results provide a valuable resource for further elucidating the posttranscriptional regulatory mechanisms of chicken liver and adipose fat development or deposition.

Multi-omics integrative analysis provided new insights into alkaline stress in alfalfa

Saline-alkali stress is one of the main abiotic stresses that limits plant growth. Salt stress has been widely studied, but alkaline salt degradation caused by NaHCO3 has rarely been investigated. In the present study, the alfalfa cultivar ‘Zhongmu No. 1′ was treated with 50 mM NaHCO3 (0, 4, 8, 12 and 24 h) to study the resulting enzyme activity and changes in mRNA, miRNA and metabolites in the roots. The results showed that the enzyme activity changed significantly after alkali stress treatment. The genomic analysis revealed 14,970 differentially expressed mRNAs (DEMs), 53 differentially expressed miRNAs (DEMis), and 463 differentially accumulated metabolites (DAMs). Combined analysis of DEMs and DEMis revealed that 21 DEMis negatively regulated 42 DEMs. In addition, when combined with Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of DEMs and DAMs, we found that phenylpropanoid biosynthesis, flavonoid biosynthesis, starch and sucrose metabolism and plant hormone signal transduction played important roles in the alkali stress response. The results of this study further elucidated the regulatory mechanism underlying the plant response to alkali stress and provided valuable information for the breeding of new saline-alkaline tolerance plant varieties.

Hydroxyapatite nanoparticle improves ovine oocyte developmental capacity by alleviating oxidative stress in response to vitrification stimuli

The wide application of ovine oocyte vitrification is limited by its relatively low efficiency. Nanoparticle is potentially to be used in cryopreservation technology for its unique characteristics with high biocompatibility, potent antioxidant property as well as superiority in membrane permeation and heat transduction. However, the effect of nanoparticle on ovine oocyte cryopreservation as well as the underlying mechanism has not been systematically evaluated. The objective of this study was to investigate the impact of nanoparticles on ovine oocytes cryopreservation and further identify the underlying mechanism. Firstly, the effects of Hydroxyapatite (HA) and Fe3O4 nanoparticles on the developmental potential of vitrified ovine oocytes were determined, and the results showed that neither HA (VC = 85.95 ± 6.23 % vs. VH = 92.47 ± 8.11 %, P > 0.05) nor Fe3O4 (VC = 85.95 ± 6.23 % vs. VF = 89.39 ± 6.32 %, P > 0.05) had adverse effect on the survival rate of vitrified-thawed oocytes. Notably, both HA (VC = 77.78 ± 0.09 % vs. VH = 44.00 ± 0.09 %, P＜0.01) and Fe3O4 (VC = 77.78 ± 0.09 % vs. VF = 51.67 ± 0.15 %, P＜0.01) nanoparticles effectively reduced the level of oocyte apoptosis after freezing and thawing. What's more, HA could significantly improve the cleavage rate of frozen oocytes (VC = 33.79 ± 2.83 % vs. VH = 59.54 ± 4.13 %, P＜0.05). Moreover, reduced reactive oxygen species (ROS) level (VC = 13.66 ± 0.47 vs. VH = 12.61 ± 0.53, P < 0.05), increased glutathione (GSH) content (VC = 60.69 ± 7.89 vs. VH = 87.92 ± 1.05, P < 0.05) and elevated mitochondrial membrane potential (MMP) level (VC = 1.43 ± 0.04 vs. VH = 1.63 ± 0.01，P＜0.01) were observed in oocytes treated with HA nanoparticles when compared with that of the control group. Furthermore, Smart-RNA sequence technology was utilized to identify differentially expressed mRNAs (DEMs) induced by nanoparticles during cryopreservation. When compared with the control counterparts, a total of 721 DEMs (309 up-regulated and 412 down-regulated mRNAs) were identified in oocytes treated with HA, while 702 DEMs (480 up-regulated and 222 down-regulated mRNAs) were identified in oocytes treated with Fe3O4. A comparison of DEMs showed that total 692 mRNAs were expressed in oocytes treated with HA and Fe3O4. Notably, we discovered that 15 mRNAs were specially highly expressed in oocytes treated with HA, and Focal adhesion signaling pathway mainly contributed to the improved ovine oocyte quality after vitrification by alleviating oxidative stress.

Regulatory Roles of Long Non-Coding RNAs in Arterial Stiffness and Hypertension: Insights from Two African American Studies.

Arterial stiffness, commonly assessed via pulse wave velocity (PWV), is marked by reduced arterial elasticity and serves as a significant risk factor for cardiovascular disease and an early indicator of hypertension. This study investigated the regulatory roles of long non-coding RNAs (lncRNAs) in modulating mRNAs associated with arterial stiffness and hypertension, with a particular focus on African Americans, a population disproportionately impacted by hypertension. We utilized whole-blood transcriptome sequencing data from two African American (AA) cohorts with high hypertension prevalence: the GENE-FORECAST study (436 subjects) and the MH-GRID study (179 subjects). Our objectives were to: (1) identify lncRNAs and mRNAs differentially expressed (DE) between the upper and lower tertiles of PWV, (2) determine DE lncRNAs associated with the expression levels of each DE mRNA, and (3) link the lncRNA-modulated mRNAs to hypertension across both datasets. Differential expression analysis revealed 1,035 DE mRNAs and 31 DE lncRNAs between upper and lower PWV groups. Then lncRNA-mRNA pairs significantly associated were identified, involving 31 unique lncRNAs and 1,034 unique mRNAs. Finally, 22 of the lncRNA-modulated mRNAs initially linked to PWV were found associated with hypertension, in both datasets. Interestingly, 30 lncRNAs were linked to the expression of UCP2 (Uncoupling Protein 2), a gene implicated in oxidative stress and endothelial function. Our findings underscore the significant roles of lncRNAs in regulating gene expression associated with arterial stiffness and hypertension. The differential expression of UCP2 in relation to PWV and hypertension, along with its potential regulation by lncRNAs, offers valuable insights into the molecular mechanisms underlying arterial stiffness and its connection with hypertension.

New potential diagnostic markers for verrucous hyperplasia and verrucous carcinoma based on RNA-sequencing data

Verrucous carcinoma (VC) is a rare subtype of squamous cell carcinoma (SCC) characterized by its histological presentation as a low-grade tumor with no potential for metastasis, setting it apart from invasive SCC. However, distinguishing VC from its benign counterpart, verrucous hyperplasia (VH), is challenging due to their clinical and morphological similarities. Despite the importance of accurate diagnosis for determining treatment strategies, diagnosis of VH and VC relied only on lesion recurrence after resection. To address this challenge, we generated RNA profiling data from tissue samples of VH and VC patients to identify novel diagnostic markers. We analyzed differentially expressed (DE) mRNA and long non-coding RNA (lncRNA) in tissue samples from VH and VC patients. Additionally, ChIP-X Enrichment Analysis 3 (ChEA3) was conducted to identify the top five transcription factors potentially regulating the expression of DE mRNAs in VH and VC. Our analysis of mRNA and lncRNA expression profiles in VH and VC provides insights into the underlying molecular characteristics of these diseases and offers potential new diagnostic markers. The identification of specific DE genes and lncRNAs may enable clinicians to more accurately differentiate between VH and VC, leading to better treatment choices.

Identification of key genes and long non‑coding RNA expression profiles in osteoporosis with rheumatoid arthritis based on bioinformatics analysis

BackgroundAlthough rheumatoid arthritis (RA) is a chronic systemic tissue disease often accompanied by osteoporosis (OP), the molecular mechanisms underlying this association remain unclear. This study aimed to elucidate the pathogenesis of RA and OP by identifying differentially expressed mRNAs (DEmRNAs) and long non-coding RNAs (lncRNAs) using a bioinformatics approach.MethodsExpression profiles of individuals diagnosed with OP and RA were retrieved from the Gene Expression Omnibus database. Differential expression analysis was conducted. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) pathway enrichment analyses were performed to gain insights into the functional categories and molecular/biochemical pathways associated with DEmRNAs. We identified the intersection of common DEmRNAs and lncRNAs and constructed a protein-protein interaction (PPI) network. Correlation analysis between the common DEmRNAs and lncRNAs facilitated the construction of a coding-non-coding network. Lastly, serum peripheral blood mononuclear cells (PBMCs) from patients with RA and OP, as well as healthy controls, were obtained for TRAP staining and qRT-PCR to validate the findings obtained from the online dataset assessments.ResultsA total of 28 DEmRNAs and 2 DElncRNAs were identified in individuals with both RA and OP. Chromosomal distribution analysis of the consensus DEmRNAs revealed that chromosome 1 had the highest number of differential expression genes. GO and KEGG analyses indicated that these DEmRNAs were primarily associated with " platelets (PLTs) degranulation”, “platelet alpha granules”, “platelet activation”, “tight junctions” and “leukocyte transendothelial migration”, with many genes functionally related to PLTs. In the PPI network, MT-ATP6 and PTGS1 emerged as potential hub genes, with MT-ATP6 originating from mitochondrial DNA. Co-expression analysis identified two key lncRNA-mRNA pairs: RP11 − 815J21.2 with MT − ATP6 and RP11 − 815J21.2 with PTGS1. Experimental validation confirmed significant differential expression of RP11-815J21.2, MT-ATP6 and PTGS1 between the healthy controls and the RA + OP groups. Notably, knockdown of RP11-815J21.2 attenuated TNF + IL-6-induced osteoclastogenesis.ConclusionsThis study successfully identified shared dysregulated genes and potential therapeutic targets in individuals with RA and OP, highlighting their molecular similarities. These findings provide new insights into the pathogenesis of RA and OP and suggest potential avenues for further research and targeted therapies.