Construction of a TF-miRNA-gene feed-forward loop network predicts biomarkers and potential drugs for myasthenia gravis

Inflammatory Gene Signature Identified by Machine Algorithms Reveals Novel Biomarkers of Coronary Artery Disease

Purpose

Inflammatory activation of immune cells plays a pivotal role in the development of coronary artery diseases (CAD). This study aims to investigate the immune responses of peripheral blood mononuclear cells (PBMCs) in CAD and to identify novel signature genes and biomarkers using machine learning algorithms.

Methods

The GSE113079 dataset was analyzed and differentially expressed genes (DEGs) were identified between CAD and normal samples. The intersection of DEGs with inflammation-related genes was used to identify the differentially expressed inflammation-related genes (DIRGs). Then, the receiver operating characteristic (ROC) curves were plotted for each DIRG, and those with an area under the curve (AUC) greater than 0.9 were selected for subsequent analysis. Furthermore, machine learning algorithms were employed to identify biomarkers. A nomogram was developed based on these biomarkers. The CIBERSORT algorithm and Wilcoxon test method were used to analyze the differences in immune cells between the CAD and normal samples. The identified biomarkers were validated in PBMCs from patients with CAD and in atherosclerotic aortas from ApoE-/- mice.

Results

A total of 574 DEGs were identified between CAD and normal samples. From this intersection, 29 DIRGs were identified, of which 14 DIRGs (PTGER1, IL17RC, KLKB1, GPR32, ADM, NUPR1, SCN9A, IL17B, CX3CL1, FFAR3, PYDC2, SYT11, RORA, and GPR31) exhibited high diagnostic efficacy. Four biomarkers (ADM, NUPR1, PTGER1, and PYDC2) were identified using Support Vector Machine (SVM). Ten types of immune cells, including CD8+ T cells, regulatory T cells (Tregs), and resting NK cells, showed significant differences between the CAD and normal groups. Furthermore, increased levels of ADM, NUPR1, PTGER1, and PYDC2 were validated in PBMCs isolated from CAD patients. In addition, ADM, NUPR1, and PTGER1 were upregulated in the mouse atherosclerotic aorta.

Conclusion

Our findings revealed novel inflammatory gene signatures of CAD that could be potential biomarkers for the early diagnosis of CAD.

Unveiling cell-type-specific microRNA networks through alternative polyadenylation in glioblastoma

BACKGROUND

Glioblastoma multiforme (GBM) is characterized by its cellular complexity, with a microenvironment consisting of diverse cell types, including oligodendrocyte precursor cells (OPCs) and neoplastic CD133 + radial glia-like cells. This study focuses on exploring the distinct cellular transitions in GBM, emphasizing the role of alternative polyadenylation (APA) in modulating microRNA-binding and post-transcriptional regulation.

RESULTS

Our research identified unique APA profiles that signify the transitional phases between neoplastic cells and OPCs, underscoring the importance of APA in cellular identity and transformation in GBM. A significant finding was the disconnection between differential APA events and gene expression alterations, indicating that APA operates as an independent regulatory mechanism. We also highlighted the specific genes in neoplastic cells and OPCs that lose microRNA-binding sites due to APA, which are crucial for maintaining stem cell characteristics and DNA repair, respectively. The constructed networks of microRNA-transcription factor-target genes provide insights into the cellular mechanisms influencing cancer cell survival and therapeutic resistance.

CONCLUSIONS

This study elucidates the APA-driven regulatory framework within GBM, spotlighting its influence on cell state transitions and microRNA network dynamics. Our comprehensive analysis using single-cell RNA sequencing data to investigate the microRNA-binding sites altered by APA profiles offers a robust foundation for future research, presenting a novel approach to understanding and potentially targeting the complex molecular interplay in GBM.

Advances in the genetics of myasthenia gravis: insights from cutting-edge neuroscience research

Myasthenia gravis (MG) is an autoimmune disorder involving complex interactions between genetic and environmental factors. Genome-wide association studies (GWAS), transcriptome-wide association studies (TWAS), and other methods have identified multiple novel susceptibility loci and genes, providing crucial insights into the genetic etiology of MG. Moreover, the pivotal roles of epigenetic mechanisms, such as DNA methylation, histone modifications, and non-coding RNAs, in the pathogenesis of MG are gradually being unveiled. This review comprehensively summarizes the latest advances in MG genetic research, focusing on the discovery and validation of susceptibility genes, genetic heterogeneity and subtype-specific genetic factors, gene-environment interactions, epigenetic mechanisms, and progress in genetics-based diagnostic and prognostic biomarkers.

The intricate dance of non-coding RNAs in myasthenia gravis pathogenesis and treatment

Myasthenia gravis (MG) stands as a perplexing autoimmune disorder affecting the neuromuscular junction, driven by a multitude of antibodies targeting postsynaptic elements. However, the mystery of MG pathogenesis has yet to be completely uncovered, and its heterogeneity also challenges diagnosis and treatment. Growing evidence shows the differential expression of non-coding RNAs (ncRNAs) in MG has played an essential role in the development of MG in recent years. Remarkably, these aberrantly expressed ncRNAs exhibit distinct profiles within diverse clinical subgroups and among patients harboring various antibody types. Furthermore, they have been implicated in orchestrating the production of inflammatory cytokines, perturbing the equilibrium of T helper 1 cells (Th1), T helper 17 cells (Th17), and regulatory T cells (Tregs), and inciting B cells to generate antibodies. Studies have elucidated that certain ncRNAs mirror the clinical severity of MG, while others may hold therapeutic significance, showcasing a propensity to return to normal levels following appropriate treatments or potentially foretelling the responsiveness to immunosuppressive therapies. Notably, the intricate interplay among these ncRNAs does not follow a linear trajectory but rather assembles into a complex network, with competing endogenous RNA (ceRNA) emerging as a prominent hub in some cases. This comprehensive review consolidates the landscape of dysregulated ncRNAs in MG, briefly delineating their pivotal role in MG pathogenesis. Furthermore, it explores their promise as prospective biomarkers, aiding in the elucidation of disease subtypes, assessment of disease severity, monitoring therapeutic responses, and as novel therapeutic targets.

Investigating the link between miR-34a-5p and TLR6 signaling in sepsis-induced ARDS

Acute lung injury (ALI)/acute respiratory distress syndrome (ARDS) are lung complications diagnosed by impaired gaseous exchanges leading to mortality. From the diverse etiologies, sepsis is a prominent contributor to ALI/ARDS. In the present study, we retrieved sepsis-induced ARDS mRNA expression profile and identified 883 differentially expressed genes (DEGs). Next, we established an ARDS-specific weighted gene co-expression network (WGCN) and picked the blue module as our hub module based on highly correlated network properties. Later we subjected all hub module DEGs to form an ARDS-specific 3-node feed-forward loop (FFL) whose highest-order subnetwork motif revealed one TF (STAT6), one miRNA (miR-34a-5p), and one mRNA (TLR6). Thereafter, we screened a natural product library and identified three lead molecules that showed promising binding affinity against TLR6. We then performed molecular dynamics simulations to evaluate the stability and binding free energy of the TLR6-lead molecule complexes. Our results suggest these lead molecules may be potential therapeutic candidates for treating sepsis-induced ALI/ARDS. In-silico studies on clinical datasets for sepsis-induced ARDS indicate a possible positive interaction between miR-34a and TLR6 and an antagonizing effect on STAT6 to promote inflammation. Also, the translational study on septic mice lungs by IHC staining reveals a hike in the expression of TLR6. We report here that miR-34a actively augments the effect of sepsis on lung epithelial cell apoptosis. This study suggests that miR-34a promotes TLR6 to heighten inflammation in sepsis-induced ALI/ARDS.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03700-1.

Major regulators of the multi-step metastatic process are potential therapeutic targets for breast cancer management

Metastasis is a multi-step process that leads to the dissemination of tumor cells to new sites and, consequently, to multi-organ neoplasia. Although most lethal breast cancer cases are related to metastasis occurrence, little is known about the dysregulation of each step, and clinicians still lack reliable therapeutic targets for metastasis impairment. To fill these gaps, we constructed and analyzed gene regulatory networks for each metastasis step (cell adhesion loss, epithelial-to-mesenchymal transition, and angiogenesis). Through topological analysis, we identified E2F1, EGR1, EZH2, JUN, TP63, and miR-200c-3p as general hub-regulators, FLI1 for cell-adhesion loss specifically, and TRIM28, TCF3, and miR-429 for angiogenesis. Applying the FANMOD algorithm, we identified 60 coherent feed-forward loops regulating metastasis-related genes associated with distant metastasis-free survival prediction. miR-139-5p, miR-200c-3p, miR-454-3p, and miR-1301-3p, among others, were the FFL's mediators. The expression of the regulators and mediators was observed to impact overall survival and to go along with metastasis occurrence. Lastly, we selected 12 key regulators and observed that they are potential therapeutic targets for canonical and candidate antineoplastics and immunomodulatory drugs, like trastuzumab, goserelin, and calcitriol. Our results highlight the relevance of miRNAs in mediating feed-forward loops and regulating the expression of metastasis-related genes. Altogether, our results contribute to understanding the multi-step metastasis complexity and identifying novel therapeutic targets and drugs for breast cancer management.

Identification of functional TF-miRNA-hub gene regulatory network associated with ovarian endometriosis

Endometriosis (EMs), one of the most common gynecological diseases, seriously affects the health and wellness of women; however, the underlying pathogenesis remains unclear. This study focused on dysregulated genes and their predicted transcription factors (TFs) and miRNAs, which may provide ideas for further mechanistic research. The microarray expression dataset GSE58178, which included six ovarian endometriosis (OE) samples and six control samples, was downloaded from the Gene Expression Omnibus (GEO) to identify differentially expressed genes (DEGs). Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to study the cellular and organism-level functions of DEGs. The protein-protein interaction (PPI) network was built and visualized using Cytoscape, and modules and hub genes were explored using various algorithms. Furthermore, we predicted miRNAs and TFs of hub genes using online databases, and constructed the TF-miRNA-hub gene network. There were 124 upregulated genes and 66 downregulated genes in EMs tissues. GO enrichment analysis showed that DEGs were concentrated in reproductive structure development and collagen-containing extracellular matrix, while KEGG pathway analysis showed that glycolysis/gluconeogenesis and central carbon metabolism in cancer require further exploration. Subsequently, HIF1A, LDHA, PGK1, TFRC, and CD9 were identified as hub genes, 22 miRNAs and 34 TFs were predicted to be upstream regulators of hub genes, and these molecules were pooled together. In addition, we found three key feedback loops in the network, MYC-miR-34a-5p-LDHA, YY1-miR-155-5p-HIF1A, and RELA-miR-93-5p-HIF1A, which may be closely related to OE development. Taken together, our study structured a TF-miRNA-hub gene network to decipher the molecular mechanism of OE, which may provide novel insights for clinical diagnosis and treatment.

Evolution of miRNA-Binding Sites and Regulatory Networks in Cichlids

The divergence of regulatory regions and gene regulatory network (GRN) rewiring is a key driver of cichlid phenotypic diversity. However, the contribution of miRNA-binding site turnover has yet to be linked to GRN evolution across cichlids. Here, we extend our previous studies by analyzing the selective constraints driving evolution of miRNA and transcription factor (TF)-binding sites of target genes, to infer instances of cichlid GRN rewiring associated with regulatory binding site turnover. Comparative analyses identified increased species-specific networks that are functionally associated to traits of cichlid phenotypic diversity. The evolutionary rewiring is associated with differential models of miRNA- and TF-binding site turnover, driven by a high proportion of fast-evolving polymorphic sites in adaptive trait genes compared with subsets of random genes. Positive selection acting upon discrete mutations in these regulatory regions is likely to be an important mechanism in rewiring GRNs in rapidly radiating cichlids. Regulatory variants of functionally associated miRNA- and TF-binding sites of visual opsin genes differentially segregate according to phylogeny and ecology of Lake Malawi species, identifying both rewired, for example, clade-specific and conserved network motifs of adaptive trait associated GRNs. Our approach revealed several novel candidate regulators, regulatory regions, and three-node motifs across cichlid genomes with previously reported associations to known adaptive evolutionary traits.

Esrrb Regulates Specific Feed-Forward Loops to Transit From Pluripotency Into Early Stages of Differentiation

Characterization of pluripotent states, in which cells can both self-renew or differentiate, with the irreversible loss of pluripotency, are important research areas in developmental biology. Although microRNAs (miRNAs) have been shown to play a relevant role in cellular differentiation, the role of miRNAs integrated into gene regulatory networks and its dynamic changes during these early stages of embryonic stem cell (ESC) differentiation remain elusive. Here we describe the dynamic transcriptional regulatory circuitry of stem cells that incorporate protein-coding and miRNA genes based on miRNA array expression and quantitative sequencing of short transcripts upon the downregulation of the Estrogen Related Receptor Beta (Esrrb). The data reveals how Esrrb, a key stem cell transcription factor, regulates a specific stem cell miRNA expression program and integrates dynamic changes of feed-forward loops contributing to the early stages of cell differentiation upon its downregulation. Together these findings provide new insights on the architecture of the combined transcriptional post-transcriptional regulatory network in embryonic stem cells.

Potential Small Molecules for Therapy of Lupus Nephritis Based on Genetic Effect and Immune Infiltration

Lupus nephritis (LN) is the most common and significant complication of systemic lupus erythematosus (SLE) due to its poor prognosis and mortality rates in SLE patients. There is a critical need for new drugs as the pathogenesis of LN remains to be elucidated and immunosuppressive therapy comes with many deficiencies. In this study, 23 hub genes (IFI6, PLSCR1, XAF1, IFI16, IFI44, MX1, IFI44L, IFIT3, IFIT2, IFI27, DDX58, EIF2AK2, IFITM1, RTP4, IFITM3, TRIM22, PARP12, IFIH1, OAS1, HERC6, RSAD2, DDX60, and MX2) were identified through bioinformatics and network analysis and are closely related to interferon production and function. Interestingly, immune cell infiltration analysis and correlation analysis demonstrate a positive correlation between the expression of 23 hub genes and monocyte infiltration in glomeruli and M2 macrophage infiltration in the tubulointerstitium of LN patients. Additionally, the CTD database, DsigDB database, and DREIMT database were used to explore the bridging role of genes in chemicals and LN as well as the potential influence of these chemicals on immune cells. After comparison and discussion, six small molecules (Acetohexamide, Suloctidil, Terfenadine, Prochlorperazine, Mefloquine, and Triprolidine) were selected for their potential ability in treating lupus nephritis.

Construction Immune Related Feed-Forward Loop Network Reveals Angiotensin II Receptor Blocker as Potential Neuroprotective Drug for Ischemic Stroke

Ischemic stroke (IS) accounts for the leading cause of disability and mortality in China. Increasing researchers are studying the effects of neuroprotective agents on IS. However, the molecular mechanisms of feed-forward loops (FFLs) associated with neuroprotection in the pathogenesis of IS need to be further studied. A protein-protein interaction (PPI) network of IS immune genes was constructed to decipher the characters and excavate 3 hub genes (PI3K, IL6, and TNF) of immunity. Then, we identified two hub clusters of IS immune genes, and the cytokine-cytokine receptor interaction pathway was discovered on the pathway enrichment results of both clusters. Combined with GO enrichment analysis, the cytokines participate in the inflammatory response in the extracellular space of IS patients. Next, a transcription factor (TF)-miRNA-immune gene network (TMIGN) was established by extracting four regulatory pairs (TF-miRNA, TF-gene, miRNA-gene, and miRNA-TF). Then, we detected 3-node regulatory motif types in the TMIGN network. According to the criteria we set for defining 3-node motifs, the motif with the highest Z-score (3-node composite FFL) was picked as the statistically evident motif, which was merged to construct an immune-associated composite FFL motif-specific sub-network (IA-CFMSN), which contained 21 3-node FFLs composed of 13 miRNAs, 4 TFs, 9 immune genes, and 1 TF& immune gene, among which TP53 and VEGFA were prominent TF and immune gene, respectively. In addition, the immune genes in IA-CFMSN were used for identifying associated pathways and drugs to further clarify the immune regulation mechanism and neuroprotection after IS. As a result, 5 immune genes targeted by 20 drugs were identified and the Angiotensin II Receptor Blockers (ARBs) target AGTR1 was found to be a neuroprotective drug for IS. In the present study, the construction of IA-CFMSN provides IS immune-associated FFLs for further experimental studies, providing new prospects for the discovery of new biomarkers and potential drugs for IS.

Identification of the hub genes and transcription factor-miRNA axes involved in Helicobacter pylori-associated gastric cancer

It has been previously reported that transcription factor-microRNA (TF-miRNA) axes play a significant role in the carcinogenesis of several types of malignant tumor. However, there is a lack of research into the differences in the mechanism of Helicobacter pylori (HP)-positive [HP(+)] and HP-negative [HP(−)] gastric cancer. The aim of the present study was to identify the hub genes and TF-miRNA axes, and to determine the potential mechanisms involved in HP-associated gastric cancer. HP-associated mRNA and miRNA data, as well as the corresponding clinical information, was downloaded from The Cancer Genome Atlas database. Differentially expressed genes (DEGs) and DE miRNAs (DEMs) were then identified from the HP(+) and HP(−) cancer mRNA and miRNA datasets, respectively. Subsequently, gene set enrichment analysis and the protein-protein interaction (PPI) networks were investigated using the ClusterProfiler packages. Lastly, TF-miRNA-DEG networks were constructed using the miRWalk online tool. A total of 1,050 DEGs and 13 DEMs were identified from the normalized mRNA and miRNA expression datasets, respectively. In addition, 180 Gene Ontology terms and 30 Kyoto Encyclopedia of Genes and Genomes pathways were found to be enriched, while 6 hub genes were identified from the PPI analysis. Furthermore, 7 TF-miRNA interactions and 181 TF-miRNA-DEG axes were constructed using an integrated bioinformatics approach, while 2 TF-miRNA interactions (ZEB1-miRNA-144-3p and PAX2-miRNA-592) were confirmed using reverse transcription-quantitative PCR in samples from enrolled patients. Moreover, the ZEB1-miRNA-144-3p axis was further validated based on dual luciferase reporter assay results. In summary, an integrated bioinformatics approach was used to screen the significant molecular and regulatory axes, which may provide a novel direction to investigate the pathogenesis of gastric cancer associated with HP.

A review on the role of epidermal growth factor signaling in the development, progression and treatment of cervical cancer

The sub-committee constituted by the Indian Council of Medical Research (ICMR) for the management of cervical cancer (CC) detailed in the consensus document (2016) reported CC as a significant cause of morbidity and mortality in women. The incidence of an increase in CC and associated mortality in women is a major cause of cancer. To date, human papilloma viral (HPV) infection accounts for more than 99% of CC. However, there are individuals infected with HPV do not develop CC. There is a greater correlation between HPV infection and upregulation of the epidermal growth factor receptor (EGFR) signaling cascade during the initiation, sustenance, and progression of CC. Therefore, EGFR is often targeted to treat CC using tyrosine kinase inhibitors (TKIs) and monoclonal antibodies (mAB). The current review analyzed the existing clinical/pre-clinical studies and the significance of EGFR abundance using the Kaplan-Meier (KM) survival plot analysis for disease-free survival (DFS) and overall survival (OS). We performed a series of bioinformatics analyses to screen the crucial role of the EGFR gene in CC. Further, different transcription factors that are dysregulated due to EGFR abundance and their relevance were determined using computational tools in this review. Endogenous microRNAs (miRNA) that undergo changes due to alterations in EGFR during CC were identified using computational database and consolidated the information obtained with the published in the area of miRNA and EGFR with special reference to the initiation, sustenance and progression of CC. The current review aims to consolidate contemporary approaches for targeting CC using EGFR and highlight the current role of miRNA and genes that are differently regulated during CC involving EGFR mutations. Potential resistance to the available EGFR therapies such as TKIs and mABs and the need for better therapies are also extensively reviewed for the development of newer therapeutic molecules with better efficacy.

A Review on the Role of Non-Coding RNAs in the Pathogenesis of Myasthenia Gravis

Myasthenia gravis (MG) is an autoimmune condition related to autoantibodies against certain proteins in the postsynaptic membranes in the neuromuscular junction. This disorder has a multifactorial inheritance. The connection between environmental and genetic factors can be established by epigenetic factors, such as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs). XLOC_003810, SNHG16, IFNG-AS1, and MALAT-1 are among the lncRNAs with a possible role in the pathoetiology of MG. Moreover, miR-150-5p, miR-155, miR-146a-5p, miR-20b, miR-21-5p, miR-126, let-7a-5p, and let-7f-5p are among miRNAs whose roles in the pathogenesis of MG has been assessed. In the current review, we summarize the impact of miRNAs and lncRNAs in the development or progression of MG.

Identification of the regulatory role of lncRNA HCG18 in myasthenia gravis by integrated bioinformatics and experimental analyses

BACKGROUND

Long non-coding RNAs (lncRNAs), functioning as competing endogenous RNAs (ceRNAs), have been reported to play important roles in the pathogenesis of autoimmune diseases. However, little is known about the regulatory roles of lncRNAs underlying the mechanism of myasthenia gravis (MG). The aim of the present study was to explore the roles of lncRNAs as ceRNAs associated with the progression of MG.

METHODS

MG risk genes and miRNAs were obtained from public databases. Protein-protein interaction (PPI) network analysis and module analysis were performed. A lncRNA-mediated module-associated ceRNA (LMMAC) network, which integrated risk genes in modules, risk miRNAs and predicted lncRNAs, was constructed to systematically explore the regulatory roles of lncRNAs in MG. Through performing random walk with restart on the network, HCG18/miR-145-5p/CD28 ceRNA axis was found to play important roles in MG, potentially. The expression of HCG18 in MG patients was detected using RT-PCR. The effects of HCG18 knockdown on cell proliferation and apoptosis were determined by CCK-8 assay and flow cytometry. The interactions among HCG18, miR-145-5p and CD28 were explored by luciferase assay, RT-PCR and western blot assay.

RESULTS

Based on PPI network, we identified 9 modules. Functional enrichment analyses revealed these modules were enriched in immune-related signaling pathways. We then constructed LMMAC network, containing 25 genes, 50 miRNAs, and 64 lncRNAs. Through bioinformatics algorithm, we found lncRNA HCG18 as a ceRNA, might play important roles in MG. Further experiments indicated that HCG18 was overexpressed in MG patients and was a target of miR-145-5p. Functional assays illustrated that HCG18 suppressed Jurkat cell apoptosis and promoted cell proliferation. Mechanistically, knockdown of HCG18 inhibited the CD28 mRNA and protein expression levels in Jurkat cells, while miR-145-5p inhibitor blocked the reduction of CD28 expression induced by HCG18 suppression.

CONCLUSION

We have reported a novel HCG18/miR-145-5p/CD28 ceRNA axis in MG. Our findings will contribute to a deeper understanding of the molecular mechanism of and provide a novel potential therapeutic target for MG.