Integrated bioinformatic analysis of microarray data reveals shared gene signature between MDS and AML

Deciphering molecular heterogeneity in pediatric AML using a cancer vs. normal transcriptomic approach

BACKGROUND

Still 30-40% of pediatric acute myeloid leukemia (pedAML) patients relapse. Delineation of the transcriptomic profile of leukemic subpopulations could aid in a better understanding of molecular biology and provide novel biomarkers.

METHODS

Using microarray profiling and quantitative PCR validation, transcript expression was measured in leukemic stem cells (LSC, n = 24) and leukemic blasts (L-blast, n = 25) from pedAML patients in comparison to hematopoietic stem cells (HSCs, n = 19) and control myeloblasts (C-blast, n = 20) sorted from healthy subjects. Gene set enrichment analysis was performed to identify relevant gene set enrichment signatures, and functional protein associations were identified by STRING analysis.

RESULTS

Highly significantly overexpressed genes in LSC and L-blast were identified with a vast majority not studied in AML. CDKN1A, CFP, and CFD (LSC) and HOMER3, CTSA, and GADD45B (L-blast) represent potentially interesting biomarkers and therapeutic targets. Eleven LSC downregulated targets were identified that potentially qualify as tumor suppressor genes, with MYCT1, PBX1, and PTPRD of highest interest. Inflammatory and immune dysregulation appeared to be perturbed biological networks in LSC, whereas dysregulated metabolic profiles were observed in L-blast.

CONCLUSION

Our study illustrates the power of taking into account cell population heterogeneity and reveals novel targets eligible for functional evaluation and therapy in pedAML.

IMPACT

Novel transcriptional targets were discovered showing a significant differential expression in LSCs and blasts from pedAML patients compared to their normal counterparts from healthy controls. Deregulated pathways, including immune and metabolic dysregulation, were addressed for the first time in children, offering a deeper understanding of the molecular pathogenesis. These novel targets have the potential of acting as biomarkers for risk stratification, follow-up, and targeted therapy. Multiple LSC-downregulated targets endow tumor suppressor roles in other cancer entities, and further investigation whether hypomethylating therapy could result into LSC eradication in pedAML is warranted.

Cell lineage-specific methylome and genome alterations in gout

In this study, we examined data from 69 gout patients and 1,455 non-gout controls using a MethylationEPIC BeadChip assay and Illumina HiSeq platform to identify lineage-specific epigenetic alterations and associated genetic factors that contributed to gouty inflammation. Cell lineage-specific differentially methylated sites were identified using CellDMC after adjusting for sex, age, alcohol drinking, smoking status, and smoking history (total pack-years). Different cell lineages displayed distinct differential methylation. Ingenuity Pathway Analysis and NetworkAnalyst indicated that many differential methylated sites were associated with interleukin-1β expression in monocytes. On the UCSC Genome Browser and WashU Epigenome Browser, metabolic trait, cis-methylation quantitative trait loci, genetic, and functional annotation analyses identified nine methylation loci located in interleukin-1β-regulating genes (PRKCZ, CIDEC, VDAC1, CPT1A, BIRC2, BRCA1, STK11, and NLRP12) that were associated specifically with gouty inflammation. All nine sites mapped to active regulatory elements in monocytes. MoLoTool and ReMap analyses indicated that the nine methylation loci overlapped with binding sites of several transcription factors that regulated interleukin-1β production and gouty inflammation. Decreases in PRKCZ and STK11 methylation were also associated with higher numbers of first-degree relatives who also had gout. The gouty-inflammation specific methylome and genome alterations could potentially aid in the identification of novel therapeutic targets.

IGF‑IR promotes clonal cell proliferation in myelodysplastic syndromes via inhibition of the MAPK pathway

Type 1 insulin-like growth factor receptor (IGF-IR) signaling is considered to serve a key role in the development of cancer. However, the effects of IGF-IR on the malignant characteristics of myelodysplastic syndrome (MDS) clonal cells remains to be determined. In the present study it was demonstrated that knockdown of IGF-IR reduced the proliferation and increased the apoptosis of MDS/leukemia cells. Integrated analysis of gene expression profiles using bioinformatics identified the MAPK signaling pathway as a critical downstream factor of IGF-IR, and this was confirmed in vitro using western blotting which revealed that IGF-IR knockdown significantly increased the expression of activated MAPK. Furthermore, IGF-IR signaling was inhibited to investigate the potential of IGF-IR as a therapeutic target of MDS. The results revealed that the IGF-IR inhibitor picropodophyllin (PPP) inhibited cell proliferation, promoted cell apoptosis and arrested the cell cycle at the G2/M phase in MDS/leukemia cells. Similar to the effects of IGF-IR knockdown, PPP treatment also increased MAPK signaling in vitro. In conclusion, IGF-IR may serve as a potential therapeutic target of MDS.

Circulating Small Noncoding RNAs Have Specific Expression Patterns in Plasma and Extracellular Vesicles in Myelodysplastic Syndromes and Are Predictive of Patient Outcome

Myelodysplastic syndromes (MDS) are hematopoietic stem cell disorders with large heterogeneity at the clinical and molecular levels. As diagnostic procedures shift from bone marrow biopsies towards less invasive techniques, circulating small noncoding RNAs (sncRNAs) have become of particular interest as potential novel noninvasive biomarkers of the disease. We aimed to characterize the expression profiles of circulating sncRNAs of MDS patients and to search for specific RNAs applicable as potential biomarkers. We performed small RNA-seq in paired samples of total plasma and plasma-derived extracellular vesicles (EVs) obtained from 42 patients and 17 healthy controls and analyzed the data with respect to the stage of the disease, patient survival, response to azacitidine, mutational status, and RNA editing. Significantly higher amounts of RNA material and a striking imbalance in RNA content between plasma and EVs (more than 400 significantly deregulated sncRNAs) were found in MDS patients compared to healthy controls. Moreover, the RNA content of EV cargo was more homogeneous than that of total plasma, and different RNAs were deregulated in these two types of material. Differential expression analyses identified that many hematopoiesis-related miRNAs (e.g., miR-34a, miR-125a, and miR-150) were significantly increased in MDS and that miRNAs clustered on 14q32 were specifically increased in early MDS. Only low numbers of circulating sncRNAs were significantly associated with somatic mutations in the SF3B1 or DNMT3A genes. Survival analysis defined a signature of four sncRNAs (miR-1237-3p, U33, hsa_piR_019420, and miR-548av-5p measured in EVs) as the most significantly associated with overall survival (HR = 5.866, p < 0.001). In total plasma, we identified five circulating miRNAs (miR-423-5p, miR-126-3p, miR-151a-3p, miR-125a-5p, and miR-199a-3p) whose combined expression levels could predict the response to azacitidine treatment. In conclusion, our data demonstrate that circulating sncRNAs show specific patterns in MDS and that their expression changes during disease progression, providing a rationale for the potential clinical usefulness of circulating sncRNAs in MDS prognosis. However, monitoring sncRNA levels in total plasma or in the EV fraction does not reflect one another, instead, they seem to represent distinctive snapshots of the disease and the data should be interpreted circumspectly with respect to the type of material analyzed.