Shortcut barcoding and early pooling for scalable multiplex single-cell reduced-representation CpG methylation sequencing at single nucleotide resolution.

DNA methylation is essential for a wide variety of biological processes, yet the development of a highly efficient and robust technology remains a challenge for routine single-cell analysis. We developed a multiplex scalable single-cell reduced representation bisulfite sequencing (msRRBS) technology. It allows cell-specific barcoded DNA fragments of individual cells to be pooled before bisulfite conversion, free of enzymatic modification or physical capture of the DNA ends, and achieves read mapping rates of 62.5 ± 3.9%, covering 60.0 ± 1.4% of CpG islands and 71.6 ± 1.6% of promoters in K562 cells. Its reproducibility is shown in duplicates of bulk cells with close to perfect correlation (R = 0.97-0.99). At a low 1 Mb of clean reads, msRRBS provides highly consistent coverage of CpG islands and promoters, outperforming the conventional methods with orders of magnitude reduction in cost. Here, we use this method to characterize the distinct methylation patterns and cellular heterogeneity of six cell lines, plus leukemia and hepatocellular carcinoma models. Taking 4 h of hands-on time, msRRBS offers a unique, highly efficient approach for dissecting methylation heterogeneity in a variety of multicellular systems.

The Genetics Analysis of Molecular Pathogenesis for Alzheimer's Disease.

INTRODUCTION: The molecular pathogenesis of Alzheimer's disease (AD) is still not clear, and the relationship between gene expression profile for different brain regions has not been studied. OBJECTIVE: Bioinformatic analysis at the genetic level has become the best way for the pathogenesis research of AD, which can analyze the abovementioned relationship. METHODS: In this study, the datasets of AD were obtained from the Gene Expression Omnibus (GEO), and Qlucore Omics Explorer (QOE) software was used to screen differentially expressed genes of GSE36980 and GSE9770 and verify gene expression of GSE63060. The Gene Ontology (GO) function enrichment analysis of these selected genes was conducted by Database for Annotation, Visualization, and Integrated Discovery (DAVID), and then the gene/protein interaction network was established by STRING to find the related proteins. R language was used for drafting maps and plots. RESULTS: There were 20 differentially expressed genes related to AD selected from GSE36980 (p = 6.2e-6, q = 2.9422e-4) and GSE9770 (p = 3.3e-4, q = 0.016606). Their expression levels of the AD group were lower than those in the control group and varied among different brain regions. Cellular morphogenesis and establishment or maintenance of cell polarity were enriched, and LRRTM1 and RASAL1 were identified by the integration network. Moreover, the analysis of GSE63060 verified the expression level of LRRTM1 and RASAL1 in Alzheimer's patients, which was much lower than that in normal people aged >65 years. CONCLUSIONS: The pathogenesis of AD at molecular levels may link to cell membrane structures and signal transduction; hence, a list of 20 genes, including LRRTM1 and RASAL1,potentially are important for the discovery of treatment target or molecular marker of AD.

Identification of cystic fibrosis transmembrane conductance regulator as a prognostic marker for juvenile myelomonocytic leukemia via the whole-genome bisulfite sequencing of monozygotic twins and data mining.

Background: Linked deoxyribonucleic acid (DNA) hypermethylation investigations of promoter methylation levels of candidate genes may help to increase the progressiveness and mortality rates of juvenile myelomonocytic leukemia (JMML), which is a unique myelodysplastic/myeloproliferative neoplasm caused by excessive monocyte and granulocyte proliferation in infancy/early childhood. However, the roles of hypermethylation in this malignant disease are uncertain. Methods: Bone marrow samples from a JMML patient and peripheral blood samples from a healthy monozygotic twin and an unrelated healthy donor were collected with the informed consent of the participant's parents. Whole-genome bisulfite sequencing (WGBS) was then performed. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to analyze specific differentially methylated region (DMG) related genes. The target genes were screened with Cytoscape and Search Tool for the Retrieval of Interacting Genes/Proteins (STRING), which are gene/protein interaction databases. A data mining platform was used to examine the expression level data of the healthy control and JMML patient tissues in Gene Expression Omnibus data sets, and a survival analysis was performed for all the JMML patients. Results: The STRING analysis revealed that the red node [i.e., the cystic fibrosis transmembrane conductance regulator (CFTR)] was the gene of interest. The gene-expression microarray data set analysis suggested that the CFTR expression levels did not differ significantly between the JMML patients and healthy controls (P=0.81). A statistically significant difference was observed in the CFTR promoter methylation level but not in the CFTR gene body methylation level. The overall survival analysis demonstrated that a high level of CFTR expression was associated with a worse survival rate in patients with JMML (P=0.039). Conclusions: CFTR promoter hypermethylation may be a novel biomarker for the diagnosis, monitoring of disease progression, and prognosis of JMML.

A combined analysis of bulk and single-cell sequencing data reveals that depleted extracellular matrix and enhanced immune processes co-contribute to fluorouracil beneficial responses in gastric cancer.

Fluorouracil, also known as 5-FU, is one of the most commonly used chemotherapy drugs in the treatment of advanced gastric cancer (GC). Whereas, the presence of innate or acquired resistance largely limits its survival benefit in GC patients. Although accumulated studies have demonstrated the involvement of tumor microenvironments (TMEs) in chemo-resistance induction, so far little is known about the relevance of GC TMEs in 5-FU resistance. To this end, in this study, we investigated the relationship between TME features and 5-FU responses in GC patients using a combined analysis involving both bulk sequencing data from the TCGA database and single-cell RNA sequencing data from the GEO database. We found that depleted extracellular matrix (ECM) components such as capillary/stroma cells and enhanced immune processes such as increased number of M1 polarized macrophages/Memory T cells/Natural Killer T cells/B cells and decreased number of regulatory T cells are two important features relating to 5-FU beneficial responses in GC patients, especially in diffuse-type patients. We further validated these two features in the tumor tissues of 5-FU-benefit GC patients using immunofluorescence staining experiments. Based on this finding, we also established a Pro (63 genes) and Con (199 genes) gene cohort that could predict 5-FU responses in GC with an AUC (area under curve) score of 0.90 in diffuse-type GC patients, and further proved the partial applicability of this gene panel pan-cancer-wide. Moreover, we identified possible communications mediated by heparanase and galectin-1 which could regulate ECM remodeling and tumor immune microenvironment (TIME) reshaping. Altogether, these findings deciphered the relationship between GC TMEs and 5-FU resistance for the first time, as well as provided potential therapeutic targets and predicting rationale to overcome this chemo-resistance, which could shed some light on developing novel precision treatment strategies in clinical practice.

[Prognostic significance and risk factors of minimal residual disease ≥1% on 19th day of induction chemotherapy in children with acute lymphoblastic leukemia].

OBJECTIVE: To assess the prognostic value of minimal residual disease on 19th day of induction chemotherapy (D19 MRD) and the risk factors of D19 MRD ≥ 1% in children with acute lymphoblastic leukemia (ALL) treated following the Chinese Children's Cancer Group ALL protocol. METHODS: We retrospectively analyzed the data of 243 children with ALL diagnosed between January 1, 2015 and December 31, 2018 in the Department of Pediatrics of Nanfang Hospital (Guangzhou China). Kaplan Meier-survival analysis was performed to compare the survival time between the patients with D19 MRD < 1% and those with D19 MRD ≥ 1%; logistic regression analyisis and Chi-square test were used to identify the risk factors of D19 MRD ≥ 1%. RESULTS: Compared with those with D19 MRD ≥ 1%, the children with D19 MRD < 1% had significantly better 3-year overall survival (100% vs 90.2%, P=0.004) and event-free survival (97.6% vs 71.6%, P < 0.001). Univariate analysis showed that the odds ratio (OR) for mediastinal invasion, T-cell immunophenotype, TEL/AML1 fusion gene and the presence of blasts in peripheral blood on the 5th day were 4.47 (95%CI: 0.275-72.968, P=0.034), 5.250 (95%CI: 1.950-14.133, P=0.02), 0.330 (95%CI: 0.112-0.970, P=0.036) and 4.407 (95%CI: 1.782-10.895, P=0.01), respectively. The initial risk stratification (P < 0.001), white blood cell grades (P=0.018) and its counts (P=0.027), and the number of blasts on the 5th day (P < 0.001) were significantly different between the two groups. Multivariate analysis showed that initial risk stratification as intermediate and high risks (OR=2.889, 95% CI: 1.193-6.996) and the presence of blasts in peripheral blood on the 5th day (OR=4.477, 95% CI: 1.692-11.843) were independent risk factors for poor early treatment response. CONCLUSIONS: D19 MRD ≥ 1% is a predictor of poor prognosis in children with ALL. Mediastinal invasion, T-cell immunophenotype and the presence of blasts in peripheral blood on the 5th day are all risk factors for poor early treatment response, while TEL/AML1 fusion gene is a protective factor; the initial risk stratification as intermediate to high risk and the presence of blasts in peripheral blood on the 5th day are independent risk factors for poor early treatment response of the patients.