Integrated analysis of microRNA and transcription factors in the bone marrow of patients with acute monocytic leukemia

A novel LGALS1-depended and immune-associated fatty acid metabolism risk model in acute myeloid leukemia stem cells

Leukemia stem cells (LSCs) are recognized as the root cause of leukemia initiation, relapse, and drug resistance. Lipid species are highly abundant and essential component of human cells, which often changed in tumor microenvironment. LSCs remodel lipid metabolism to sustain the stemness. However, there is no useful lipid related biomarker has been approved for clinical practice in AML prediction and treatment. Here, we constructed and verified fatty acid metabolism-related risk score (LFMRS) model based on TCGA database via a series of bioinformatics analysis, univariate COX regression analysis, and multivariate COX regression analysis, and found that the LFMRS model could be an independent risk factor and predict the survival time of AML patients combined with age. Moreover, we revealed that Galectin-1 (LGALS1, the key gene of LFMRS) was highly expressed in LSCs and associated with poor prognosis of AML patients, and LGALS1 repression inhibited AML cell and LSC proliferation, enhanced cell apoptosis, and decreased lipid accumulation in vitro. LGALS1 repression curbed AML progression, lipid accumulation, and CD8+ T and NK cell counts in vivo. Our study sheds light on the roles of LFMRS (especially LGALS1) model in AML, and provides information that may help clinicians improve patient prognosis and develop personalized treatment regimens for AML.

Data-driven modeling of core gene regulatory network underlying leukemogenesis in IDH mutant AML

Acute myeloid leukemia (AML) is characterized by uncontrolled proliferation of poorly differentiated myeloid cells, with a heterogenous mutational landscape. Mutations in IDH1 and IDH2 are found in 20% of the AML cases. Although much effort has been made to identify genes associated with leukemogenesis, the regulatory mechanism of AML state transition is still not fully understood. To alleviate this issue, here we develop a new computational approach that integrates genomic data from diverse sources, including gene expression and ATAC-seq datasets, curated gene regulatory interaction databases, and mathematical modeling to establish models of context-specific core gene regulatory networks (GRNs) for a mechanistic understanding of tumorigenesis of AML with IDH mutations. The approach adopts a new optimization procedure to identify the top network according to its accuracy in capturing gene expression states and its flexibility to allow sufficient control of state transitions. From GRN modeling, we identify key regulators associated with the function of IDH mutations, such as DNA methyltransferase DNMT1, and network destabilizers, such as E2F1. The constructed core regulatory network and outcomes of in-silico network perturbations are supported by survival data from AML patients. We expect that the combined bioinformatics and systems-biology modeling approach will be generally applicable to elucidate the gene regulation of disease progression.

The miR-146b-3p/TNFAIP2 axis regulates cell differentiation in acute myeloid leukaemia

Our purpose is to verify that miR-146b-3p targets the downstream transcript TNFAIP2 in order to reveal the machinery underlying the miR-146b-3p/TNFAIP2 axis regulating acute myeloid leukaemia (AML) differentiation. Bioinformatics analyses were performed using multiple databases and R packages. The CD11b+ and CD14+ cell frequencies were detected using flow cytometry and immunofluorescence staining. The TNFAIP2 protein expression was evaluated using western blotting, immunocytochemistry and immunofluorescence staining. The qRT-PCR was conducted to detect the expression of TNFAIP2 and miR-146b-3p. TNFAIP2 and its correlated genes were enriched in multiple cell differentiation pathways. TNFAIP2 was upregulated upon leukaemic cell differentiation. miR-146b-3p directly targeted TNFAIP2, resulting in a decrease in TNFAIP2 expression. Forced expression of TNFAIP2 or knockdown of miR-146b-3p significantly induced the differentiation of MOLM-13 cells. In this study, we demonstrated that TNFAIP2 is a critical driver in inducing differentiation and that the miR-146b-3p/TNFAIP2 axis involves in regulating cell differentiation in AML.

Molecular mechanisms in the miR-33a/LPPR4 pathway regulating unsaturated fatty acid synthesis in bovine mammary epithelial cells

The regulatory mechanisms governing metabolism of fatty acids in cow mammary gland are crucial for establishing relationships between milk quality and fatty acid content. Both, microRNAs (miRNAs) and protein-coding genes are important factors involved in the regulation of milk fat synthesis. In this study, high-throughput sequencing of miRNAs and mRNAs in bovine mammary gland tissue was performed during peak lactation (3 samples) and late lactation (3 samples) periods to characterize expression profiles. Differential expression (DE) analyses of miRNA and mRNA and miRNA-mRNA regulatory pathway screening were performed. Two-hundred eighty regulatory miRNA-mRNA pairs were identified, including the miR-33a-lipid phosphate phosphatase-related protein type 4 (LPPR4) pathway. Bioinformatics prediction, dual-luciferase reporter system detection, qRT-PCR, and Western blotting revealed that miR-33a can directly target LPPR4 and inhibit its expression. Experiments also revealed that miR-33a promotes the synthesis of triglycerides and increases the content of unsaturated fatty acids (UFAs) in bovine mammary epithelial cells (BMECs). These results indicate that miR-33a via LPPR4 plays an important role in the regulation of milk fat synthesis and UFA levels.

Data-driven modeling of core gene regulatory network underlying leukemogenesis in IDH mutant AML

Acute myeloid leukemia (AML) is characterized by uncontrolled proliferation of poorly differentiated myeloid cells, with a heterogenous mutational landscape. Mutations in IDH1 and IDH2 are found in 20% of the AML cases. Although much effort has been made to identify genes associated with leukemogenesis, the regulatory mechanism of AML state transition is still not fully understood. To alleviate this issue, here we develop a new computational approach that integrates genomic data from diverse sources, including gene expression and ATAC-seq datasets, curated gene regulatory interaction databases, and mathematical modeling to establish models of context-specific core gene regulatory networks (GRNs) for a mechanistic understanding of tumorigenesis of AML with IDH mutations. The approach adopts a novel optimization procedure to identify the optimal network according to its accuracy in capturing gene expression states and its flexibility to allow sufficient control of state transitions. From GRN modeling, we identify key regulators associated with the function of IDH mutations, such as DNA methyltransferase DNMT1, and network destabilizers, such as E2F1. The constructed core regulatory network and outcomes of in-silico network perturbations are supported by survival data from AML patients. We expect that the combined bioinformatics and systems-biology modeling approach will be generally applicable to elucidate the gene regulation of disease progression.

Cervical Cancer Stem-Like Cell Transcriptome Profiles Predict Response to Chemoradiotherapy

Cervical cancer (CC) represents a major global health issue, particularly impacting women from resource constrained regions worldwide. Treatment refractoriness to standard chemoradiotheraphy has identified cancer stem cells as critical coordinators behind the biological mechanisms of resistance, contributing to CC recurrence. In this work, we evaluated differential gene expression in cervical cancer stem-like cells (CCSC) as biomarkers related to intrinsic chemoradioresistance in CC. A total of 31 patients with locally advanced CC and referred to Mário Penna Institute (Belo Horizonte, Brazil) from August 2017 to May 2018 were recruited for the study. Fluorescence-activated cell sorting was used to enrich CD34+/CD45- CCSC from tumor biopsies. Transcriptome was performed using ultra-low input RNA sequencing and differentially expressed genes (DEGs) using Log2 fold differences and adjusted p-value < 0.05 were determined. The analysis returned 1050 DEGs when comparing the Non-Responder (NR) (n=10) and Responder (R) (n=21) groups to chemoradiotherapy. These included a wide-ranging pattern of underexpressed coding genes in the NR vs. R patients and a panel of lncRNAs and miRNAs with implications for CC tumorigenesis. A panel of biomarkers was selected using the rank-based AUC (Area Under the ROC Curve) and pAUC (partial AUC) measurements for diagnostic sensitivity and specificity. Genes overlapping between the 21 highest AUC and pAUC loci revealed seven genes with a strong capacity for identifying NR vs. R patients (ILF2, RBM22P2, ACO16722.1, AL360175.1 and AC092354.1), of which four also returned significant survival Hazard Ratios. This study identifies DEG signatures that provide potential biomarkers in CC prognosis and treatment outcome, as well as identifies potential alternative targets for cancer therapy.