miRNA-mRNA Relationships Research Articles

Identification of key genes and functions in lung metastasis of osteosarcoma based on bioinformatics

To screen the differentially expressed genes of lung metastasis of osteosarcoma by bioinformatics, and explore their functions and regulatory networks. The data set of GSE14359 was screened from GEO database(http://www.ncbi.nlm.nih.gov/gds) and the differentially expressed gene(DEG) was identified using GEO2R online tool. Download osteosarcoma disease related miRNAs from the online HMMD database(http://www.cuilab.cn/hmdd) and then FunRich software was used to predict the target gene, intersects with DEG to obtains the target gene. The miRNA-mRNA relationship pairs were formed according to the targeted joints, then the data was imported into Cytoscape for visualization, DAVID was used to performe GO and KEGG analysis on target genes, STRING was used to construct PPI network, Cytoscape visualization, CytoHubba plug-in screening central genes and online website for expression and survival analysis. Total 704 DEGs were identified, consisting of 477 up-regulated genes and 227 down regulated genes. FunRich predicted 7 888 mRNAs and 343 target genes were obtained through intersection of the two. KEGG analysis showed that it was mainly involved in focal adhesion, ECM receptor interaction, TNF signal pathway, PI3K-Akt signal pathway, IL-17 signal pathway and MAPK signal pathway. Ten central genes (CCNB1, CHEK1, AURKA, DTL, RRM2, MELK, CEP55, FEN1, KPNA2, TYMS) were identified as potential key genes. Among them, CCNB1, DTL, MELK were highly correlated with poor prognosis. The key genes and functional pathways identified in this study may be helpful to understand the molecular mechanism of the occurrence and progression of lung metastases from osteosarcoma, and provide potential therapeutic targets.

Exploring the molecular mechanisms by which per- and polyfluoroalkyl substances induce polycystic ovary syndrome through in silico toxicogenomic data mining

The pathogeny of polycystic ovary syndrome (PCOS) is intricate, with endocrine disruptors (EDCs) being acknowledged as significant environmental factors. Research has shown a link between exposure to per- and polyfluoroalkyl substances (PFAS) and the development and progression of PCOS, although the precise mechanism is not fully understood. This study utilized toxicogenomics and comparative toxicogenomics databases to analyze data and investigate how PFAS mixtures may contribute to the development of PCOS. The results indicated that 74 genes are associated with both PFAS exposure and PCOS progression. Enrichment analysis suggested that cell cycle regulation and steroid hormone synthesis may be crucial pathways through which PFAS mixtures participate in the development of PCOS, involving important genes such as CCNB1 and SRD5A1. Furthermore, the study identified transcription factors (TFs) and miRNAs that may be involved in the onset and progression of PCOS, constructing regulatory networks encompassing TFs-mRNA interactions and miRNA-mRNA relationships to elucidate their regulatory roles in gene expression. By utilizing data mining techniques based on toxicogenomic databases, this study provides relatively comprehensive insights into the association between exposure factors and diseases compared to traditional toxicology studies. These findings offer new perspectives for further in vivo or in vitro investigations and contribute to understanding the pathogenesis of PCOS, thereby providing valuable references for identifying clinical treatment targets.

Systematic identification and characterization of microRNAs with target genes involved in high night temperature stress at the filling stage of rice.

High night temperature stress is one of the main environmental factors affecting rice yield and quality. More and more evidence shows that microRNA (miRNA) plays an important role in various abiotic stresses. However, the molecular network of miRNA regulation on rice tolerance to high night temperatures remains unclear. Here, small RNA, transcriptome and degradome sequencing were integrated to identify differentially expressed miRNAs, genes, and key miRNA-target gene pairs in rice heat-sensitive and heat-tolerant lines at the filling stage suffering from high night temperature stress. It was discovered that there were notable differences in the relative expression of 102 miRNAs between the two rice lines under stress. Meanwhile, 5263 and 5405 mRNAs were differentially expressed in the heat-sensitive line and heat-tolerant line, and functional enrichment analysis revealed that these genes were involved in heat-related processes and pathways. The miRNAs-mRNAs target relationship was further verified by degradome sequencing. Eventually, 49 miRNAs-222 mRNAs target pairs with reverse expression patterns showed significant relative expression changes between the heat-tolerant and the heat-sensitive line, being suggested to be responsible for the heat tolerance difference of these two rice lines. Functional analysis of these 222 mRNA transcripts showed that high night temperature-responsive miRNAs targeted these mRNAs involved in many heat-related biological processes, such as transcription regulation, chloroplast regulation, mitochondrion regulation, protein folding, hormone regulation and redox process. This study identified possible miRNA-mRNA regulation relationships in response to high night temperature stress in rice and potentially contributed to heat resistance breeding of rice in the future.

Identification of key genes in sepsis-induced cardiomyopathy based on integrated bioinformatical analysis and experiments in vitro and in vivo.

Sepsis is a life-threatening disease that damages multiple organs and induced by the host's dysregulated response to infection with high morbidity and mortality. Heart remains one of the most vulnerable targets of sepsis-induced organ damage, and sepsis-induced cardiomyopathy (SIC) is an important factor that exacerbates the death of patients. However, the underlying genetic mechanism of SIC disease needs further research. The transcriptomic dataset, GSE171564, was downloaded from NCBI for further analysis. Gene expression matrices for the sample group were obtained by quartile standardization and log2 logarithm conversion prior to analysis. The time series, protein-protein interaction (PPI) network, and functional enrichment analysis via Gene Ontology and KEGG Pathway Databases were used to identify key gene clusters and their potential interactions. Predicted miRNA-mRNA relationships from multiple databases facilitated the construction of a TF-miRNA-mRNA regulatory network. In vivo experiments, along with qPCR and western blot assays, provided experimental validation. The transcriptome data analysis between SIC and healthy samples revealed 221 down-regulated, and 342 up-regulated expressed genes across two distinct clusters. Among these, Tpt1, Mmp9 and Fth1 were of particular significance. Functional analysis revealed their role in several biological processes and pathways, subsequently, in vivo experiments confirmed their overexpression in SIC samples. Notably, we found TPT1 play a pivotal role in the progression of SIC, and silencing TPT1 showed a protective effect against LPS-induced SIC. In our study, we demonstrated that Tpt1, Mmp9 and Fth1 have great potential to be biomarker of SIC. These findings will facilitated to understand the occurrence and development mechanism of SIC.

Expression Profiling of Adipogenic and Anti-Adipogenic MicroRNA Sequences following Methylmercury Exposure in Caenorhabditis elegans.

MicroRNA (miRNA) are important regulators of gene expression that respond not only to developmental and pathological cues, but also to environmental stimuli. Dyslipidemia is a hallmark of metabolic conditions and has been shown to significantly affect the expression of circulating miRNA sequences. Recently, our lab has shown that the environmental toxicant methylmercury (MeHg) causes dyslipidemia in the Caenorhabditis elegans model organism. While 10 and 20 μM MeHg increases the expression of adipogenic transcription factors and lipid-binding proteins in worms, there is limited information on how the toxicant affects the miRNA regulators of these genes. We hypothesized that MeHg would increase the expression of adipogenic miRNA sequences and/or decrease the expression of anti-adipogenic miRNA sequences. We further hypothesized that the target mRNA sequences for the miRNAs affected by MeHg would be consequently altered. We selected three potentially adipogenic (mir-34, mir-124, and mir-355) and three potentially anti-adipogenic (mir-240, mir-786, and let-7) miRNA sequences homologous to known human miRNA sequences altered in obesity, and quantified their levels 24 h and 48 h post MeHg treatment. At 24 h post exposure, MeHg significantly increased expression of both the adipogenic and anti-adipogenic miRNA sequences 1.5-3x above untreated control. By 48 h post exposure, only the adipogenic miRNA sequences were elevated, while the anti-adipogenic miRNA sequences were decreased by 50% compared to untreated control. These data suggest that there are developmental changes in miRNA expression over time following MeHg exposure. We next selected one target mRNA sequence for each miRNA sequence based on miRNA-mRNA relationships observed in humans. MeHg altered the gene expression of all the target genes assayed. Except for mir-34, all the tested miRNA-mRNA sequences showed a conserved relationship between nematode and humans. To determine whether the selected miRNA sequences were involved in lipid accumulation in response to MeHg, lipid storage was investigated in transgenic worm strains that lacked the specific miRNA strains. Of the six strains investigated, only the mir-124 and let-7 mutant worms had lipid storage levels that were statistically different from wild type, suggesting that these two sequences can be potential mediators of MeHg-induced lipid dysregulation.

Open a new epoch of arsenic trioxide investigation: ATOdb

Arsenic trioxide (ATO) is a great discovery in the treatment of acute promyelocytic leukemia (APL), which has been used in an increasing number of malignant diseases. Systematic integrative analysis will help to precisely understand the mechanism of ATO and find new combined drugs. Therefore, we developed a one-stop comprehensive database of ATO named ATOdb by manually compiling a wealth of experimentally supported ATO-related data from 3479 articles, and integrated analysis tools. The current version of ATOdb contains 8373 associations among 2300 ATO targets, 80 conditions and 262 combined drugs. Each entry in ATOdb contains detailed information on ATO targets, therapeutic/side effects, systems, cell names, cell types, regulations, detection methods, brief descriptions, references, etc. Furthermore, ATOdb also provides data visualization and analysis results such as the drug similarities, protein-protein interactions, and miRNA-mRNA relationships. An easy-to-use web interface was deployed in ATOdb for users to easily browse, search and download the data. In conclusion, ATOdb will serve as a valuable resource for in-depth study of the mechanism of ATO, discovery of new drug combination strategies, promotion of rational drug use and individualized treatments. ATOdb is freely available at http://bio-bigdata.hrbmu.edu.cn/ATOdb/index.jsp.

Construction of lncRNA-miRNA-mRNA regulatory network in severe asthmatic bronchial epithelial cells: A bioinformatics study.

Asthma is a chronic respiratory disease caused by environment-host interactions. Bronchial epithelial cells (BECs) are the first line of defense against environmental toxins. However, the mechanisms underlying the role of BECs in severe asthma (SA) are not yet fully understood. Long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) have been shown to play important roles in the regulation of gene expression in the pathogenesis of SA. In this study, bioinformatics was used for the first time to reveal the lncRNA-miRNA-mRNA regulatory network of BECs in SA. Five mRNA datasets of bronchial brushing samples from patients with SA and healthy controls (HC) were downloaded from the Gene Expression Omnibus (GEO) database. A combination of the Venn diagram and robust rank aggregation (RRA) method was used to identify core differentially expressed genes (DEGs). Protein-protein interaction (PPI) analysis of core DEGs was performed to screen hub genes. The miRDB, miRWalk, and ENCORI databases were used to predict the miRNA-mRNA relationships, and the ENCORI and starBase v2.0 databases were used to predict the upstream lncRNAs of the miRNA-mRNA relationships. Four core DEGs were identified: carcinoembryonic antigen-related cell adhesion molecule 5 (CEACAM5), interleukin-1 receptor type 2 (IL1R2), trefoil factor 3 (TFF3), and vascular endothelial growth factor A (VEGFA). These 4 core DEGs indicated that SA was not significantly associated with sex. Enrichment analysis showed that the MAPK, Rap1, Ras, PI3K-Akt and Calcium signaling pathways may serve as the principal pathways of BECs in SA. A lncRNA-miRNA-mRNA regulatory network of the severe asthmatic bronchial epithelium was constructed. The top 10 competing endogenous RNAs (ceRNAs) were FGD5 antisense RNA 1 (FGD5-AS1), metastasis associated lung adenocarcinoma transcript 1 (MALAT1), X inactive specific transcript (XIST), HLA complex group 18 (HCG18), small nucleolar RNA host gene 16 (SNHG16), has-miR-20b-5p, has-miR-106a-5p, hsa-miR-106b-5p, has-miR-519d-3p and Fms related receptor tyrosine kinase 1 (FLT1). Our study revealed a potential mechanism for the lncRNA-miRNA-mRNA regulatory network in BECs in SA.

Myocardial RNA Sequencing Reveals New Potential Therapeutic Targets in Heart Failure with Preserved Ejection Fraction.

Heart failure with preserved ejection fraction (HFpEF) represents a global health challenge, with limited therapies proven to enhance patient outcomes. This makes the elucidation of disease mechanisms and the identification of novel potential therapeutic targets a priority. Here, we performed RNA sequencing on ventricular myocardial biopsies from patients with HFpEF, prospecting to discover distinctive transcriptomic signatures. A total of 306 differentially expressed mRNAs (DEG) and 152 differentially expressed microRNAs (DEM) were identified and enriched in several biological processes involved in HF. Moreover, by integrating mRNA and microRNA expression data, we identified five potentially novel miRNA-mRNA relationships in HFpEF: the upregulated hsa-miR-25-3p, hsa-miR-26a-5p, and has-miR4429, targeting HAPLN1; and NPPB mRNA, targeted by hsa-miR-26a-5p and miR-140-3p. Exploring the predicted miRNA-mRNA interactions experimentally, we demonstrated that overexpression of the distinct miRNAs leads to the downregulation of their target genes. Interestingly, we also observed that microRNA signatures display a higher discriminative power to distinguish HFpEF sub-groups over mRNA signatures. Our results offer new mechanistic clues, which can potentially translate into new HFpEF therapies.

The possible molecular mechanism underlying the involvement of the variable shear factor QKI in the epithelial-mesenchymal transformation of oesophageal cancer.

Based on the GEO, TCGA and GTEx databases, we reveal the possible molecular mechanism of the variable shear factor QKI in epithelial mesenchymal transformation (EMT) of oesophageal cancer. Based on the TCGA and GTEx databases, the differential expression of the variable shear factor QKI in oesophageal cancer samples was analysed, and functional enrichment analysis of QKI was performed based on the TCGA-ESCA dataset. The percent-spliced in (PSI) data of oesophageal cancer samples were downloaded from the TCGASpliceSeq database, and the genes and variable splicing types that were significantly related to the expression of the variable splicing factor QKI were screened out. We further identified the significantly upregulated circRNAs and their corresponding coding genes in oesophageal cancer, screened the EMT-related genes that were significantly positively correlated with QKI expression, predicted the circRNA-miRNA binding relationship through the circBank database, predicted the miRNA-mRNA binding relationship through the TargetScan database, and finally obtained the circRNA-miRNA-mRNA network through which QKI promoted the EMT process. Compared with normal control tissue, QKI expression was significantly upregulated in tumour tissue samples of oesophageal cancer patients. High expression of QKI may promote the EMT process in oesophageal cancer. QKI promotes hsa_circ_0006646 and hsa_circ_0061395 generation by regulating the variable shear of BACH1 and PTK2. In oesophageal cancer, QKI may promote the production of the above two circRNAs by regulating variable splicing, and these circRNAs further competitively bind miRNAs to relieve the targeted inhibition of IL-11, MFAP2, MMP10, and MMP1 and finally promote the EMT process. Variable shear factor QKI promotes hsa_circ_0006646 and hsa_circ_0061395 generation, and downstream related miRNAs can relieve the targeted inhibition of EMT-related genes (IL11, MFAP2, MMP10, MMP1) and promote the occurrence and development of oesophageal cancer, providing a new theoretical basis for screening prognostic markers of oesophageal cancer patients.

Exploratory integrated analysis of circulating exosomal miRNA and tissue mRNA related to long-term physical activity for more than 25years: a bioinformatics study.

Physical activity exerts various positive effects on both physical and mental health. Although the comprehensive expression profiles of each microRNA (miRNA) or messenger RNA (mRNA) related to physical activity have already been reported, the association between miRNA and mRNA remains unclear. Here, the integrated study was conducted to comprehensively explore the potential miRNA-mRNA relationships related to long-term physical activity over 25years. Genome-wide public deposited mRNA expression data of adipose tissue (GSE20536) from six same-sex twin pairs (no information regarding gender) and of skeletal muscle tissue (GSE20319) from ten same-sex twin pairs (four female twin pairs) were used, and differentially expressed mRNAs (DEMs) related to discordant leisure-time physical activity for 30years were identified using GEO2R. Overlapped mRNAs between DEMs and predicted possible target mRNAs, based on a previous study and TargetScan tool, were then identified and used as long-term physical activity-related mRNAs targeted by miRNAs. In adipose tissue, 36 mRNAs and 42 mRNAs were identified as upregulated or downregulated DEMs, respectively. Based on the results of the overlapped analysis between DEMs and predicted possible target mRNAs targeted by miRNAs, 15 upregulated mRNAs, including NDRG4, FAM13A, ST3GAL6, and AFF1, and 10 downregulated mRNAs, including RPL14, LBP, and GLRX, were identified. In muscle tissue, three downregulated mRNAs overlapped with the predicted target mRNAs targeted by miRNAs. Fifteen upregulated mRNAs in adipose tissue showed a tendency to enrich in "Cardiovascular" in GAD_DISEASE_CLASS category. Potential miRNA-mRNA relationships related to long-term physical activity over 25years were identified through bioinformatics analysis.

Application of lncRNA-miRNA-mRNA ceRNA network analysis in the treatment of androgenic alopecia.

Long noncoding RNAs (lncRNAs) can be used as competitive endogenous RNAs (ceRNAs) to bind to microRNAs (miRNAs) to regulate gene expression. Previous studies have demonstrated that ceRNAs play an important role in the development of tumors. However, it is not clear whether the lncRNA-miRNA-mRNA ceRNA network plays a role in androgenic alopecia (AGA). The hair follicles of three AGA patients and three healthy individuals were collected for high-throughput whole transcriptome sequencing to screen for differentially expressed lncRNAs. Differentially expressed lncRNA target genes were subjected to databases to predict miRNA-mRNA and lncRNA-miRNA relationship pairs, and a ceRNA network was constructed using Cytoscape software. Relative expression was verified by real-time quantitative reverse transcription-polymerase chain reaction (qRT-PCR). 84 lncRNAs were significantly differentially expressed between the hair follicles of AGA patients and those of healthy individuals; 30 were upregulated, and 54 were downregulated. The top 10 upregulated lncRNAs were ENST00000501520, ENST00000448179, ENST00000318291, ENST00000568280, ENST00000561121, ENST00000376609, ENST00000602414, ENST00000573866, ENST00000513358, and ENST00000564194. The top 10 downregulated lncRNAs were ENST00000566804, ENST00000561973, ENST00000587680, ENST00000569927, ENST00000340444, ENST00000424345, ENST00000589787, NR_024344, NR_073026, and NR_110001. The qRT-PCR validation results and receiver-operating characteristic curve analysis indicated that one upregulated lncRNA, LOXL1-AS1 (ENST00000564194), had the most significant clinical diagnostic potential. After further analysis, it was concluded that LOXL1-AS1 could be used as a sponge to target hsa-miR-5193, thereby regulating TP53 expression. The ceRNA network-regulating AGA was constructed through high-throughput sequencing. Our study also identified a key lncRNA that is possibly related to the AGA pathological process.

The Effects of Transcranial Focused Ultrasound Stimulation of Nucleus Accumbens on Neuronal Gene Expression and Brain Tissue in High Alcohol-Preferring Rats.

We investigated the effect of low-intensity focused ultrasound (LIFU) on gene expression related to alcohol dependence and histological effects on brain tissue. We also aimed at determining the miRNA-mRNA relationship and their pathways in alcohol dependence-induced expression changes after focused ultrasound therapy. We designed a case-control study for 100days of observation to investigate differences in gene expression in the short-term stimulation group (STS) and long-term stimulation group (LTS) compared with the control sham group (SG). The study was performed in our Experimental Research Laboratory. 24 male high alcohol-preferring rats 63 to 79days old, weighing 270 to 300g, were included in the experiment. LTS received 50-day LIFU and STS received 10-day LIFU and 40-day sham stimulation, while the SG received 50-day sham stimulation. In miRNA expression analysis, it was found that LIFU caused gene expression differences in NAc. Significant differences were found between the groups for gene expression. Compared to the SG, the expression of 454 genes in the NAc region was changed in the STS while the expression of 382 genes was changed in the LTS. In the LTS, the expression of 32 genes was changed in total compared to STS. Our data suggest that LIFU targeted on NAc may assist in the treatment of alcohol dependence, especially in the long term possibly through altering gene expression. Our immunohistochemical studies verified that LIFU does not cause any tissue damage. These findings may lead to new studies in investigating the efficacy of LIFU for the treatment of alcohol dependence and also for other psychiatric disorders.

Construction of a transcription factor-miRNA-mRNA interactive network elucidates underlying pathogenesis for osteosarcoma and validation by qRT-PCR.

Purpose:Osteosarcoma is characterized by features of rapid growth and early metastasis with a poor prognosis. The aim of our research is to investigate the potential transcription factor (TF)-miRNA-mRNA regulatory mechanism in osteosarcoma utilizing bioinformatics methods and validate by qRT-PCR.Methods:The microRNA (miRNA) expression profiling datasets (GSE28423 and GSE65071) and mRNA expression profiling dataset GSE33382 were collected from the Gene Expression Omnibus (GEO) database. Differentially expressed miRNAs (DEMs) and differentially expressed genes (DEGs) were screened using the limma package. Then, the TransmiR v2.0, miRDB, and Targetscan 7.2 database were applied for the acquisition of TF-miRNA and miRNA-mRNA interaction relationships, respectively. Finally, we built a TF-miRNA-mRNA interactive network. Furthermore, survival analysis was performed to identify sub-network with prognostic value and validate through qRT-PCR.Results:Eight overlapping DEMs and 682 DEGs were identified. Based on bioinformatics methods, 30 TF-miRNA interaction pairs and 25 miRNA-mRNA interaction pairs were screened. Finally, we constructed a TF-miRNA-mRNA regulatory network. Furthermore, laminin subunit gamma 1 (LAMC1) and thrombospondin-1 (THBS1), which involved in the network, were determined to have prognostic value and the corresponding subnetwork was identified. qRT-PCR results showed that LAMC1 mRNA expression was higher in osteosarcoma cells.Conclusion:Based on the survival analysis, a TF-miRNA–mRNA sub-network, that is TFs (SPI1, HEY1, and CEBPB)-hsa-miR-338-3p-target genes (LAMC1 and THBS1) was established. In conclusion, the construction of a potential TF-related regulatory network will help elucidate the underlying pathological mechanisms of osteosarcoma, and may provide novel insights for the diagnosis and treatment of osteosarcoma.

Transcriptional analysis of microRNAs related to unsaturated fatty acid synthesis by interfering bovine adipocyte ACSL1 gene.

Long-chain fatty acyl-CoA synthase 1 (ACSL1) plays a vital role in the synthesis and metabolism of fatty acids. The proportion of highly unsaturated fatty acids in beef not only affects the flavor and improves the meat’s nutritional value. In this study, si-ACSL1 and NC-ACSL1 were transfected in bovine preadipocytes, respectively, collected cells were isolated on the fourth day of induction, and then RNA-Seq technology was used to screen miRNAs related to unsaturated fatty acid synthesis. A total of 1,075 miRNAs were characterized as differentially expressed miRNAs (DE-miRNAs), of which the expressions of 16 miRNAs were upregulated, and that of 12 were downregulated. Gene ontology analysis indicated that the target genes of DE-miRNAs were mainly involved in biological regulation and metabolic processes. Additionally, KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis identified that the target genes of DE-miRNAs were mainly enriched in metabolic pathways, fatty acid metabolism, PI3K-Akt signaling pathway, glycerophospholipid metabolism, fatty acid elongation, and glucagon signaling pathway. Combined with the previous mRNA sequencing results, several key miRNA-mRNA targeting relationship pairs, i.e., novel-m0035-5p—ACSL1, novel-m0035-5p—ELOVL4, miR-9-X—ACSL1, bta-miR-677—ACSL1, miR-129-X—ELOVL4, and bta-miR-485—FADS2 were screened via the miRNA-mRNA interaction network. Thus, the results of this study provide a theoretical basis for further research on miRNA regulation of unsaturated fatty acid synthesis in bovine adipocytes.

Comprehensive analysis of fifteen hub genes to identify a promising diagnostic model, regulated networks, and immune cell infiltration in acute kidney injury.

Acute kidney injury is a common clinical problem with no sensitive and specific diagnostic biomarkers and definitive treatments. The underlying molecular mechanisms of acute kidney injury are unclear. Therefore, it is pivotal to explore the underlying mechanisms and screen for novel diagnostic biomarkers, and therapeutic targets. The present study identified 15 hub genes by WGCNA analysis. LASSO-based logistic regression analysis was used to select key features and construct a diagnostic model of AKI. In addition, GO and KEGG analyses were performed and TF-mRNA and miRNA-mRNA network analysis and immune infiltration analysis of hub genes were performed to reveal the underlying mechanisms of AKI. A diagnostic model was constructed by LASSO-based logistic regression analysis and was validated by RT-qPCR based on 15 hub genes. GO and KEGG analyses revealed DEGs were enriched in oxidation-reduction process, cell adhesion, proliferation, migration, and metabolic process. The enriched TFs were BRD2, EP300, ETS1, MYC, SPI1, and ZNF263. The enriched miRNAs were miR-181c-5p, miR-218-5p, miR-485-5p, miR-532-5p and miR-6884-5p. The immune infiltration analysis showed that Macrophages M2 was decreasing significantly revealing a protective factor for further AKI treatment. The present study identified 15 hub genes based on WGCNA. Development and validation of a potentially diagnostic model based on 15 hub genes. In addition, exploring the interaction between transcriptional factors and 15 hub genes, and miRNA-mRNA relationship pairs. Furthermore, immune infiltration analysis was performed by analyzing gene expression profiles of AKI. Our study provides some basis for further experimental studies.

Expression profile of long non-coding RNA in inner Mongolian cashmere goat with putative roles in hair follicles development.

The hair follicle is a complex skin accessory organ, which determines hair growth. Long non-coding RNAs (lncRNAs) have been proven to play an important role in hair follicle development, but their specific mechanism is still unclear. In this study, high-throughput sequencing was used to obtain the expression profiles of lncRNA in the hair follicles of Inner Mongolian cashmere goats at different embryonic stages (45, 55, 65, and 75 days), and a total of 6,630 lncRNA were identified. According to the rules of hair follicle development, we combined miRNA and mRNA databases (published) and predicted lncRNA-miRNA, miRNA-mRNA, and lncRNA-mRNA interaction pairs in the 45 vs. 75 comparison group. We obtained 516 lncRNA-mRNA, 1,011 lncRNA-miRNA, and 7,411 miRNA-mRNA relationship pairs. Finally, target genes were analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), and it was found that they were mainly enriched in the Wnt signaling pathway and PI3K-Akt signaling pathway related to hair follicle development, indicating that lncRNA may interact with miRNA/mRNA to directly or indirectly regulate the expression of genes related to hair follicle development. Dual-luciferase reporter gene analysis showed that lncRNA MSTRG.1705.1 could bind to Chi-miR-1, while lncRNA MSTRG.11809.1 had no binding site for Chi-miR-433. In conclusion, this study aims to further analyze the molecular regulation mechanism of hair follicle development and to lay a theoretical foundation for revealing the regulation mechanism of cashmere hair follicle growth.

Multiomics analysis of the impact of polychlorinated biphenyls on environmental liver disease in a mouse model

Exposure to high fat diet (HFD) and persistent organic pollutants including polychlorinated biphenyls (PCBs) is associated with liver injury in human populations and non-alcoholic fatty liver disease (NAFLD) and steatohepatitis (NASH) in animal models. Previously, exposure of HFD-fed male mice to the non-dioxin-like (NDL) PCB mixture Aroclor1260, dioxin-like (DL) PCB126, or Aroclor1260 + PCB126 co-exposure caused toxicant-associated steatohepatitis (TASH) and differentially altered the liver proteome. Here unbiased mRNA and miRNA sequencing (mRNA- and miRNA- seq) was used to identify biological pathways altered in these liver samples. Fewer transcripts and miRs were up- or down- regulated by PCB126 or Aroclor1260 compared to the combination, suggesting that crosstalk between the receptors activated by these PCBs amplifies changes in the transcriptome. Pathway enrichment analysis identified “positive regulation of Wnt/β-catenin signaling” and “role of miRNAs in cell migration, survival, and angiogenesis” for differentially expressed mRNAs and miRNAs, respectively. We evaluated the five miRNAs increased in human plasma with PCB exposure and suspected TASH and found that miR-192–5p was increased with PCB exposure in mouse liver. Although we observed little overlap between differentially expressed mRNA transcripts and proteins, biological pathway-relevant PCB-induced miRNA-mRNA and miRNA-protein inverse relationships were identified that may explain protein changes. These results provide novel insights into miRNA and mRNA transcriptome changes playing direct and indirect roles in the functional protein pathways in PCB-related hepatic lipid accumulation, inflammation, and fibrosis in a mouse model of TASH and its relevance to human liver disease in exposed populations.

MiR-27 and miR-124 target AR coregulators in prostate cancer: Bioinformatics and in vitro analysis.

The inadequate efficacy of the current treatments for metastatic prostate cancer has directed efforts to the discovery of novel therapies. MicroRNAs (miRNAs) have been considered potential therapeutic agents due to their ability to control gene expression and cellular pathways. The accurate identification of genes and pathways which are targeted by a miRNA is the first step in the therapeutic use of these molecules. In this regard, there are multiple experimental and computational methods to predict and confirm the miRNA-mRNA relationships. The targeting the androgen receptor (AR) indirectly as the most important mediator of prostate cancer has been posited to both control the disease and prevent resistance to treatment. This study aimed to identify miRNAs targeting AR coregulators. For this purpose, we examined target genes by combining miRNA-mRNA computational and experimental data from various databases. miR-27a-3p and miR-124 displayed the highest scores and were selected as miRNAs with the potential to target candidate genes. Next, three cell lines of prostate cancer including PC3, LNCAP, and DU145 were transfected with plasmids which were expressed these selected miRNAs. Then, the gene expression and cell cycle analysis were performed. A decrease was observed in cell viability in all three cell lines than the cells transfected with backbone plasmid. Furthermore, the findings indicated that miR-27a-3p and miR-124 led to a significant decrease in the expression of all genes that were studied in PC3 cell line. In addition, miR-124 caused significant the cellular arrest in the G0/G1 stage, while for miR-27a-3p, this arrest occurred was in the G2/M stage. Our results indicated that the function of a unique miRNA could be different in different cell lines with particular cancer phenotype based on the cell line stage. These findings offer the possibility of employing the miR-124 and miR-27a-3p as therapeutic agents for prostate cancer treatment.

MirTarRnaSeq: An R/Bioconductor Statistical Package for miRNA-mRNA Target Identification and Interaction Analysis

We introduce mirTarRnaSeq, an R/Bioconductor package for quantitative assessment of miRNA-mRNA relationships within sample cohorts. mirTarRnaSeq is a statistical package to explore predicted or pre-hypothesized miRNA-mRNA relationships following target prediction.We present two use cases applying mirTarRnaSeq. First, to identify miRNA targets, we examined EBV miRNAs for interaction with human and virus transcriptomes of stomach adenocarcinoma. This revealed enrichment of mRNA targets highly expressed in CD105+ endothelial cells, monocytes, CD4+ T cells, NK cells, CD19+ B cells, and CD34 cells. Next, to investigate miRNA-mRNA relationships in SARS-CoV-2 (COVID-19) infection across time, we used paired miRNA and RNA sequenced datasets of SARS-CoV-2 infected lung epithelial cells across three time points (4, 12, and 24 hours post-infection). mirTarRnaSeq identified evidence for human miRNAs targeting cytokine signaling and neutrophil regulation immune pathways from 4 to 24 hours after SARS-CoV-2 infection. Confirming the clinical relevance of these predictions, three of the immune specific mRNA-miRNA relationships identified in human lung epithelial cells after SARS-CoV-2 infection were also observed to be differentially expressed in blood from patients with COVID-19. Overall, mirTarRnaSeq is a robust tool that can address a wide-range of biological questions providing improved prediction of miRNA-mRNA interactions.

TLR-9 agonist CpG activates macrophages and modulates expression of microRNAs involved in cancer drug resistance and various signaling pathways

Abstract TLR-9 agonist, CpG induces pro-inflammatory immune response. We studied the immune response elicited in C57Bl/6 female mice by repeated stimulation of TLR-9 through multiple doses of CpG administration. Complete blood count analysis after CpG treatment revealed peripheral pancytopenia and the spleen showed splenomegaly. Flow cytometry analyses of spleen and liver mononuclear cells disclosed macrophage activation, evidenced by substantial increase in F4/80+ cells in the treatment group compared to control. Analysis of serum by ELISA demonstrated enhanced levels of pro-inflammatory cytokines, IL12/23, IL6 and TNFα in the CpG-treated group whereas no significant change was observed in the amount of anti-inflammatory cytokine, IL10. The spleen cells were further subjected to microarray analysis in order to examine the modulation of microRNA expression pattern following CpG treatment. Investigation of miRNA-mRNA relationships using the ingenuity pathway analysis (IPA) revealed differential expression of miRNAs associated with cancer drug resistance, HOTAIR and BAG signaling pathways along with microRNAs involved in immunomodulation. The major microRNAs altered with CpG treatment included mir27a, mir181, mir130a, mir487a, mir29b1, mir130a, mir3473b, mir21a, mir142 and mir144 targeting genes involved in drug resistance, apoptosis, cell cycle, angiogenesis and major histocompatibility complex-related proteins. Together, our studies demonstrate that chronic macrophage activation may lead to uncontrolled activation and proliferation of immune cells. (Supported by NIH grants P01AT003961, P20GM103641, R01AI129788, R01ES030144, R01AI160896 and R01AI123947). Supported by NIH grants P01AT003961, P20GM103641, R01AI129788, R01ES030144, R01AI160896 and R01AI123947