EryDB: A Transcriptomic Profile Database for Erythropoiesis and Erythroid-related Diseases.

Erythropoiesis is a finely regulated and complex process that involves multiple transformations from hematopoietic stem cells to mature red blood cells at hematopoietic sites from the embryonic to the adult stages. Investigations into its molecular mechanisms have generated a wealth of expression data, including bulk and single-cell RNA sequencing data. A comprehensively integrated and well-curated erythropoiesis-specific database will greatly facilitate the mining of gene expression data and enable large-scale research of erythropoiesis and erythroid-related diseases. Here, we present EryDB, an open-access and comprehensive database dedicated to the collection, integration, analysis, and visualization of transcriptomic data for erythropoiesis and erythroid-related diseases. Currently, the database includes expertly curated quality-assured data of 3803 samples and 1,187,119 single cells derived from 107 public studies of three species (Homo sapiens, Mus musculus, and Danio rerio), nine tissue types, and five diseases. EryDB provides users with the ability to not only browse the molecular features of erythropoiesis between tissues and species, but also perform computational analyses of single-cell and bulk RNA sequencing data, thus serving as a convenient platform for customized queries and analyses. EryDB v1.0 is freely accessible at https://ngdc.cncb.ac.cn/EryDB/home.

Transcriptomic and genomic characteristics of intrahepatic metastases of primary liver cancer.

BACKGROUND: Patients with primary multifocal hepatocellular carcinoma (HCC) have a poor prognosis and often experience a high rate of treatment failure. Multifocal HCC is mainly caused by intrahepatic metastasis (IM), and though portal vein tumor thrombosis (PVTT) is considered a hallmark of IM, the molecular mechanism by which primary HCC cells invade the portal veins remains unclear. Therefore, it is necessary to recognize the early signs of metastasis of HCC to arrange better treatment for patients. RESULTS: To determine the differential molecular features between primary HCC with and without phenotype of metastasis, we used the CIBERSORTx software to deconvolute cell types from bulk RNA-Seq based on a single-cell transcriptomic dataset. According to the relative abundance of tumorigenic and metastatic hepatoma cells, VEGFA+ macrophages, effector memory T cells, and natural killer cells, HCC samples were divided into five groups: Pro-T, Mix, Pro-Meta, NKC, and MemT, and the transcriptomic and genomic features of the first three groups were analyzed. We found that the Pro-T group appeared to retain native hepatic metabolic activity, whereas the Pro-Meta group underwent dedifferentiation. Genes highly expressed in the group Pro-Meta often signify a worse outcome. CONCLUSIONS: The HCC cohort can be well-typed and prognosis predicted according to tumor microenvironment components. Primary hepatocellular carcinoma may have obtained corresponding molecular features before metastasis occurred.

Single-cell RNA sequencing reveals the cellular and molecular changes that contribute to the progression of lung adenocarcinoma.

Pure ground glass nodules (GGNs) and solid nodules (SNs) represent early and relatively late stages of lung adenocarcinoma (LUAD) in radiology, respectively. The cellular and molecular characteristics of pure GGNs and SNs have not been comprehensively elucidated. Additionally, the mechanism driving the progression of lung adenocarcinoma from pure GGN to SN in radiology is also elusive. In this study, by analyzing the single-cell transcriptomic profiles of 76,762 cells from four pure GGNs, four SNs, and four normal tissues, we found that anti-tumor immunity mediated by NK and CD8+T cells gradually weakened with the progression of LUAD and humoral immunity mediated by plasma B cells was more active in SNs. Additionally, the proliferation ability of some special epithelial cell increased during the progression process from pure GGN to SN. Furthermore, stromal cells and M2 macrophages could assist the progression of LUAD. Through comprehensive analyses, we revealed dynamic changes in cellular components and intercellular interactions during the progression of LUAD. These findings could facilitate our understanding of LUAD and discovery of novel therapeutic targets.

Mapping the single-cell landscape of acral melanoma and analysis of the molecular regulatory network of the tumor microenvironments.

Acral melanoma (AM) exhibits a high incidence in Asian patients with melanoma, and it is not well treated with immunotherapy. However, little attention has been paid to the characteristics of the immune microenvironment in AM. Therefore, in this study, we collected clinical samples from Chinese patients with AM and conducted single-cell RNA sequencing to analyze the heterogeneity of its tumor microenvironments (TMEs) and the molecular regulatory network. Our analysis revealed that genes, such as TWIST1, EREG, TNFRSF9, and CTGF could drive the deregulation of various TME components. The molecular interaction relationships between TME cells, such as MIF-CD44 and TNFSF9-TNFRSF9, might be an attractive target for developing novel immunotherapeutic agents.

Acral melanoma is a type of cancer that affects the hands and feet. It tends to form on the palms, soles, and under the nails. It is rare in people of European descent, but in Asian populations it makes up more than half of all melanoma cases. Unlike other types of skin cancer, it does not respond well to immunotherapy, but scientists did not understand why. Historically, cancer research has focused on the genetics of whole tumors. But cancer is complicated. Malignant cells recruit other cells to help them survive and grow, and to protect them from attacks by the immune system. Together, they create their own ecosystem, called the tumor microenvironment. The exact makeup of the tumor microenvironment differs depending on the type of cancer and on the genetics of the individual. Investigating the cells that 'support' the tumor could help to explain how acral melanoma develops and why it does not respond to treatment. To address these questions, He et al. collected samples from six patients with acral melanoma and examined the genes used by more than 60,000 individual cells. This revealed nine different types of cells in the tumor microenvironment. Most were cancer cells, but there were also immune cells, blood vessel cells, skin cells, and a type of cell that makes connective tissue. He et al. also identified four genes that most likely shape the tumor microenvironment, and two gene pairs that may control some of the interactions between the cells. Investigating these early findings in more detail could open new treatment avenues for acral melanoma. The number of samples in this study was small, but it provides a starting point for future investigation. With more data, researchers could start to develop treatments that target the unique tumor microenvironment of this type of cancer.

A deep learning model for early risk prediction of heart failure with preserved ejection fraction by DNA methylation profiles combined with clinical features.

BACKGROUND: Heart failure with preserved ejection fraction (HFpEF), affected collectively by genetic and environmental factors, is the common subtype of chronic heart failure. Although the available risk assessment methods for HFpEF have achieved some progress, they were based on clinical or genetic features alone. Here, we have developed a deep learning framework, HFmeRisk, using both 5 clinical features and 25 DNA methylation loci to predict the early risk of HFpEF in the Framingham Heart Study Cohort. RESULTS: The framework incorporates Least Absolute Shrinkage and Selection Operator and Extreme Gradient Boosting-based feature selection, as well as a Factorization-Machine based neural network-based recommender system. Model discrimination and calibration were assessed using the AUC and Hosmer-Lemeshow test. HFmeRisk, including 25 CpGs and 5 clinical features, have achieved the AUC of 0.90 (95% confidence interval 0.88-0.92) and Hosmer-Lemeshow statistic was 6.17 (P = 0.632), which outperformed models with clinical characteristics or DNA methylation levels alone, published chronic heart failure risk prediction models and other benchmark machine learning models. Out of them, the DNA methylation levels of two CpGs were significantly correlated with the paired transcriptome levels (R < -0.3, P < 0.05). Besides, DNA methylation locus in HFmeRisk were associated with intercellular signaling and interaction, amino acid metabolism, transport and activation and the clinical variables were all related with the mechanism of occurrence of HFpEF. Together, these findings give new evidence into the HFmeRisk model. CONCLUSION: Our study proposes an early risk assessment framework for HFpEF integrating both clinical and epigenetic features, providing a promising path for clinical decision making.

Long noncoding RNA PCED1B-AS1 promotes erythroid differentiation coordinating with GATA1 and chromatin remodeling.

Erythropoiesis is a complex and sophisticated multi-stage process regulated by a variety of factors, including the transcription factor GATA1 and non-coding RNA. GATA1 is regarded as an essential transcriptional regulator promoting transcription of erythroid-specific genes-such as long non-coding RNAs (lncRNA). Here, we comprehensively screened lncRNAs that were potentially regulated by GATA1 in erythroid cells. We identified a novel lncRNA-PCED1B-AS1-and verified its role in promoting erythroid differentiation of K562 erythroid cells. We also predicted a model in which PCED1B-AS1 participates in erythroid differentiation via dynamic chromatin remodeling involving GATA1. The relationship between lncRNA and chromatin in the process of erythroid differentiation remains to be revealed, and in our study we have carried out preliminary explorations.

[Genetic variability and phylogenetic analysis of 39 short tandem repeat loci in Beijing Han population].

In this study, we studied the genetic polymorphisms of short tandem repeat (STR) loci from 13 CODIS and 26 non-CODIS system in Beijing Han population for the first time, and established a database of 39 STR loci whose forensic parameters were further evaluated. Our results demonstrated no significant deviation from the Hardy-Weinberg equilibrium of 39 STR loci and no pairwise linkage disequilibrium between them. The power of discriminations, expected heterozygosity, polymorphic information content, and power of exclusion of 39 STR loci ranged from 0.7740-0.9818, 0.6000-0.9350, 0.5317-0.9047 and 0.2909-0.8673. The cumulated discrimination power and cumulative probability of exclusion were 0.999999999999999999999999999999999999999964971 and 0.999999999973878, respectively. Moreover, the genetic distance was calculated based on allele frequency and phylogenetic tree was built using STR loci data from Beijing Han and other 11 Chinese ethnic groups.This study provides important basic data for Chinese forensic DNA database and population genetics database, and has important significance in carrying out forensic individual identification, paternity testing, and population genetic study.

Molecular biomarkers screened by next-generation RNA sequencing for non-sentinel lymph node status prediction in breast cancer patients with metastatic sentinel lymph nodes.

BACKGROUND: Non-sentinel lymph node (NSLN) status prediction with molecular biomarkers may make some sentinel lymph node (SLN) positive breast cancer patients avoid the axillary lymph node dissection, but the available markers remain limited. METHODS: SLN positive patients with and without NSLN invasion were selected, and genes differentially expressed or fused in SLN metastasis were screened by next-generation RNA sequencing. RESULTS: Six candidates (all ER/PR+, HER2-, Ki-67 <20%) with metastatic SLNs selected from 305 patients were equally categorized as NSLN negative and positive. We identified 103 specifically expressed genes in the NSLN negative group and 47 in the NSLN positive group. Among them, FABP1 (negative group) and CYP2A13 (positive group) were the only 2 protein-encoding genes with expression levels in the 8th to 10th deciles. Using a false discovery rate threshold of <0.05, 62 up-regulated genes and 98 down-regulated genes were discovered in the NSLN positive group. Furthermore, 10 gene fusions were identified in this group with the most frequently fused gene being IGLL5. CONCLUSIONS: The biomarkers screened in present study may broaden our understanding of the mechanisms of breast cancer metastasis to the lymph nodes and contribute to the axillary surgery selection for SLN positive patients.

Insulin-like growth factor binding protein 5 (IGFBP5) functions as a tumor suppressor in human melanoma cells.

The insulin-like growth factor binding protein 5 (IGFBP5), which is often dysregulated in human cancers, plays a crucial role in carcinogenesis and cancer development. However, the function and underlying mechanism of IGFBP5 in tumor growth and metastasis has been elusive, particularly in malignant human melanoma. Here, we reported that IGFBP5 acts as an important tumor suppressor in melanoma tumorigenicity and metastasis by a series of experiments including transwell assay, xenograft model, in vivo tumor metastasis experiment, and RNA-Seq. Overexpression of IGFBP5 in A375, a typical human melanoma cell line, inhibited cell malignant behaviors significantly, including in vitro proliferation, anchorage-independent growth, migration and invasion, as well as in vivo tumor growth and pulmonary metastasis. In addition, overexpression of IGFBP5 suppressed epithelial-mesenchymal transition (EMT), and decreased the expression of E-cadherin and the key stem cell markers NANOG, SOX2, OCT4, KLF4, and CD133. Furthermore, IGFBP5 exerts its inhibitory activities by reducing the phosphorylation of IGF1R, ERK1/2, and p38-MAPK kinases and abating the expression of HIF1α and its target genes, VEGF and MMP9. All these findings were confirmed by IGFBP5 knockdown in human melanoma cell line A2058. Taken together, these results shed light on the mechanism of IGFBP5 as a potential tumor-suppressor in melanoma progression, indicating that IGFBP5 might be a novel therapeutic target for human melanoma.

Whole transcriptome RNA-seq analysis: tumorigenesis and metastasis of melanoma.

Melanoma is the most malignant cutaneous cancer and causes over 9000 deaths annually. Because fatality rates from malignant melanoma (MM) increase dramatically upon metastasis, we investigated tumorigenesis and metastasis of MM in transcriptome analyses of three distinct cell lines that correspond with the stages of MM pathogenesis: the normal stage (HEMn-LP), the onset of MM (A375), and the metastasis stage (A2058). Using next-generation sequencing (NGS) technology, we detected asymmetrical expression of genes among the three cell lines, notably on chromosomes 9, 11, 12, and 14, suggesting their involvement in tumorigenesis and metastasis of MM. These genes were clustered into 41 categories based on their expression patterns, and their biological functions were analyzed using Ingenuity Pathway Analysis. In the top cancer-associated category, HIF1A, IL8, TERT, ONECUT1, and FOXA1 directly interacted with either transcription factors or cytokines that are known to be involved in the tumorigenesis or metastasis of other malignant tumors. The present data suggest that cytokine regulatory pathways in macrophages predominate over other pathways during the pathogenesis of MM. This study provides new targets for the downstream mechanistic studies of the tumorigenesis and metastasis of MM and demonstrates a new strategy for studies of the progression of other malignant cancers.

Characterization of miRNomes in acute and chronic myeloid leukemia cell lines.

Myeloid leukemias are highly diverse diseases and have been shown to be associated with microRNA (miRNA) expression aberrations. The present study involved an in-depth miRNome analysis of two human acute myeloid leukemia (AML) cell lines, HL-60 and THP-1, and one human chronic myeloid leukemia (CML) cell line, K562, via massively parallel signature sequencing. mRNA expression profiles of these cell lines that were established previously in our lab facilitated an integrative analysis of miRNA and mRNA expression patterns. miRNA expression profiling followed by differential expression analysis and target prediction suggested numerous miRNA signatures in AML and CML cell lines. Some miRNAs may act as either tumor suppressors or oncomiRs in AML and CML by targeting key genes in AML and CML pathways. Expression patterns of cell type-specific miRNAs could partially reflect the characteristics of K562, HL-60 and THP-1 cell lines, such as actin filament-based processes, responsiveness to stimulus and phagocytic activity. miRNAs may also regulate myeloid differentiation, since they usually suppress differentiation regulators. Our study provides a resource to further investigate the employment of miRNAs in human leukemia subtyping, leukemogenesis and myeloid development. In addition, the distinctive miRNA signatures may be potential candidates for the clinical diagnosis, prognosis and treatment of myeloid leukemias.

Transcriptome dynamics during human erythroid differentiation and development.

To explore the mechanisms controlling erythroid differentiation and development, we analyzed the genome-wide transcription dynamics occurring during the differentiation of human embryonic stem cells (HESCs) into the erythroid lineage and development of embryonic to adult erythropoiesis using high throughput sequencing technology. HESCs and erythroid cells at three developmental stages: ESER (embryonic), FLER (fetal), and PBER (adult) were analyzed. Our findings revealed that the number of expressed genes decreased during differentiation, whereas the total expression intensity increased. At each of the three transitions (HESCs-ESERs, ESERs-FLERs, and FLERs-PBERs), many differentially expressed genes were observed, which were involved in maintaining pluripotency, early erythroid specification, rapid cell growth, and cell-cell adhesion and interaction. We also discovered dynamic networks and their central nodes in each transition. Our study provides a fundamental basis for further investigation of erythroid differentiation and development, and has implications in using ESERs for transfusion product in clinical settings.

Comprehensive characterization of erythroid-specific enhancers in the genomic regions of human Krüppel-like factors.

BACKGROUND: Mapping of DNase I hypersensitive sites (DHSs) is a powerful tool to experimentally identify cis-regulatory elements (CREs). Among CREs, enhancers are abundant and predominantly act in driving cell-specific gene expression. Krüppel-like factors (KLFs) are a family of eukaryotic transcription factors. Several KLFs have been demonstrated to play important roles in hematopoiesis. However, transcriptional regulation of KLFs via CREs, particularly enhancers, in erythroid cells has been poorly understood. RESULTS: In this study, 23 erythroid-specific or putative erythroid-specific DHSs were identified by DNase-seq in the genomic regions of 17 human KLFs, and their enhancer activities were evaluated using dual-luciferase reporter (DLR) assay. Of the 23 erythroid-specific DHSs, the enhancer activities of 15 DHSs were comparable to that of the classical enhancer HS2 in driving minimal promoter (minP). Fifteen DHSs, some overlapping those that increased minP activities, acted as enhancers when driving the corresponding KLF promoters (KLF-Ps) in erythroid cells; of these, 10 DHSs were finally characterized as erythroid-specific KLF enhancers. These 10 erythroid-specific KLF enhancers were further confirmed using chromatin immunoprecipitation coupled to sequencing (ChIP-seq) data-based bioinformatic and biochemical analyses. CONCLUSION: Our present findings provide a feasible strategy to extensively identify gene- and cell-specific enhancers from DHSs obtained by high-throughput sequencing, which will help reveal the transcriptional regulation and biological functions of genes in some specific cells.

Dynamic transcriptomes of human myeloid leukemia cells.

To identify the mechanisms controlling chronic myeloid leukemia (CML) and acute myeloid leukemia (AML) in humans, we analyzed genome-wide transcription dynamics in three myeloid leukemia cell lines (K562, HL-60, and THP1) using high-throughput sequencing technology. Using KEGG analysis, we found that the ERK/MAPK, JAK-STAT and ErbB pathways promoted proliferation and metabolism in CML. However, in AML, differentiation and apoptosis blocking resulted in the accumulation of blast cells in marrow. In addition, each cell type had unique characteristics. K562 cells are an ideal model for studying erythroid differentiation and globin gene expression. The chemokine signaling pathway and Fc gamma R-mediated phagocytosis were markedly upregulated in HL-60 cells. In THP1 cells, highly expressed genes ensured strong phagocytosis by monocytes. Further, we provide a new insight into myeloid development. The abundant data sets and well-defined analysis methods will provide a resource and strategy for further investigation of myeloid leukemia.

Influence of carbon monoxide on growth and apoptosis of human umbilical artery smooth muscle cells and vein endothelial cells.

Carbon monoxide (CO) is a vasoactive molecule that is generated by vascular cells as a byproduct of heme catabolism and it plays an important physiological role in circulation system. In order to investigate whether exogenous CO can mediate the growth and proliferation of vascular cells, in this study, we used 250 parts per million (ppm) of CO to treat human umbilical artery smooth muscle cell (hUASMC) and human umbilical vein endothelial cell (HuVEC) and further evaluated the growth and apoptosis status of SMC and HuVEC. After SMC and HuVEC were exposed to CO for 7-day, the growth of SMC and HuVEC was significantly inhibited by CO in vitro on day 5 of CO exposure. And CO blocked cell cycle progress of SMC and HuVEC, more SMC and HuVEC stagnated at G0/G1 phase by flow cytometric analysis. Moreover, CO treatment inhibited SMC and HuVEC apoptosis caused by hydrogen peroxide through decreasing caspase 3 and 9 activities. To confirm the molecular mechanism of CO effect on SMC and HuVEC growth, we compared the gene expression profile in SMC and CO-treated SMC, HuVEC and CO-treated HuVEC. By microarray analysis, we found the expression level of some genes which are related to cell cycle regulation, cell growth and proliferation, and apoptosis were changed during CO exposure. We further identified that the down-regulated CDK2 contributed to arresting cell growth and the down-regulated Caspase 3 (CASP3) and Caspase 9 (CASP9) were associated with the inhibition of cell apoptosis. Therefore, CO exerts a certain growth arrest on SMC and HuVEC by inhibiting cell cycle transition from G0/G1 phase to S phase and has regulatory effect on cell apoptosis by regulating the expression of apoptosis-associated genes.

[Progress on Sp1-like and Krüppel-like factors].

Sp1-like and Krüppel-like factors (Sp1/KLFs) are a family of zinc-finger proteins that are important components of transcriptional machinery in eukaryotic cells. Sp1/KLFs have highly conserved DNA binding domains at the carboxyl terminus that have three tandem zinc-finger motifs. Their amino terminal regions vary greatly and contain transcription regulatory domains that interact with co-regulators. By regulating the expression of a large number of genes containing GC-rich or CACCC promoters in a tissue-, cell- and developmental stage-specific manner, Sp1/KLFs may be involved in many biological processes, including cell proliferation, differentiation, apoptosis and neoplastic transformation. In this article, we reviewed the structure, molecular mechanisms and biological functions of Sp1/KLFs.