Search

Gene selection is an indispensable step for analyzing noisy and high-dimensional single-cell RNA-seq (scRNA-seq) data. Compared with the commonly used variance-based methods, by mimicking the human maker selection in the 2D visualization of cells, a new feature selection method called HRG (Highly Regional Genes) is proposed to find the informative genes, which show regional expression patterns in the cell-cell similarity network. We mathematically find the optimal expression patterns that can maximize the proposed scoring function. In comparison with several unsupervised methods, HRG shows high accuracy and robustness, and can increase the performance of downstream cell clustering and gene correlation analysis. Also, it is applicable for selecting informative genes of sequencing-based spatial transcriptomic data.

Human experts are good at summarizing explicit strong patterns from small samples, while deep learning models can learn implicit weak patterns from big data. Biologists have traditionally described the sequence patterns of promoters via transcription factor binding sites (TFBSs), while the flanking sequences among TFBSs, which can also significantly influence promoter activity and function, remain largely uncharacterized. Thus, current synthetic promoters are mainly designed by the manipulation of TFBSs, while the flanking sequence is often chosen arbitrarily or by previous experience due to a lack of well-summarized optimization criteria. Here, we introduced an AI-aided promoter design framework, DeepSEED, that employs expert knowledge and deep learning methods to efficiently design synthetic promoters that have various desirable functions. DeepSEED incorporates the user-defined cis-regulatory sequences as ‘seeds’ and generates flanking sequences that match the ‘seeds’. We showed that DeepSEED can automatically capture k-mer frequencies and DNA shape features from active promoters in the training set and efficiently optimize the flanking sequences to better match desired properties in synthetic promoters. We validated the effectiveness of this framework for diverse synthetic promoter design tasks in both prokaryotic and eukaryotic cells. DeepSEED successfully designedE. coliconstitutive, isopropyl-beta-D-1-thiogalactopyranoside (IPTG)-inducible, and mammalian cell doxycycline (Dox)-inducible promoters with significant performance improvements, suggesting that DeepSEED has potential as an efficient AI-aided flanking sequence optimization approach for promoter design that may be of great benefit in synthetic biology applications.Availabilityhttps://github.com/WangLabTHU/deepseed

Discovering DNA regulatory sequence motifs and their relative positions is vital to understanding the mechanisms of gene expression regulation. Although deep convolutional neural networks (CNNs) have achieved great success in predicting cis-regulatory elements, the discovery of motifs and their combinatorial patterns from these CNN models has remained difficult. We show that the main difficulty is due to the problem of multifaceted neurons which respond to multiple types of sequence patterns. Since existing interpretation methods were mainly designed to visualize the class of sequences that can activate the neuron, the resulting visualization will correspond to a mixture of patterns. Such a mixture is usually difficult to interpret without resolving the mixed patterns. We propose the NeuronMotif algorithm to interpret such neurons. Given any convolutional neuron (CN) in the network, NeuronMotif first generates a large sample of sequences capable of activating the CN, which typically consists of a mixture of patterns. Then, the sequences are “demixed” in a layer-wise manner by backward clustering of the feature maps of the involved convolutional layers. NeuronMotif can output the sequence motifs, and the syntax rules governing their combinations are depicted by position weight matrices organized in tree structures. Compared to existing methods, the motifs found by NeuronMotif have more matches to known motifs in the JASPAR database. The higher-order patterns uncovered for deep CNs are supported by the literature and ATAC-seq footprinting. Overall, NeuronMotif enables the deciphering of cis-regulatory codes from deep CNs and enhances the utility of CNN in genome interpretation.

Heterochromatin remodeling is critical for various cell processes. In particular, the “loss of heterochromatin” phenotype in cellular senescence is associated with the process of aging and age-related disorders. Although biological processes of senescent cells, including senescence-associated heterochromatin foci (SAHF) formation, chromosome compaction, and redistribution of key proteins, have been closely associated with high-order chromatin structure, the relationship between the high-order chromatin reorganization and the loss of heterochromatin phenotype during senescence has not been fully understood. By using senescent and deep senescent fibroblasts induced by DNA damage harboring the “loss of heterochromatin” phenotype, we observed progressive 3D reorganization of heterochromatin during senescence. Facultative and constitutive heterochromatin marked by H3K27me3 and H3K9me3, respectively, show different alterations. Facultative heterochromatin tends to switch from the repressive B-compartment to the active A-compartment, whereas constitutive heterochromatin shows no significant changes at the compartment level but enhanced interactions between themselves. Both types of heterochromatin show increased chromatin accessibility and gene expression leakage during senescence. Furthermore, increased chromatin accessibility in potential CTCF binding sites accompanies the establishment of novel loops in constitutive heterochromatin. Finally, we also observed aberrant expression of repetitive elements, including LTR (long terminal repeat) and satellite classes. Overall, facultative and constitutive heterochromatin show both similar and distinct multiscale alterations in the 3D map, chromatin accessibility, and gene expression leakage. This study provides an epigenomic map of heterochromatin reorganization during senescence.

Motivation Intermediately methylated regions occupy a significant fraction of the human genome and are closely associated with epigenetic regulations or cell-type deconvolution of bulk data. However, these regions show distinct methylation patterns, corresponding to different biological mechanisms. Although there have been some metrics developed for investigating these regions, the high noise sensitivity limits the utility for distinguishing distinct methylation patterns. Results We proposed a method named MeConcord to measure local methylation concordance across reads and CpG sites, respectively. MeConcord showed the most stable performance in distinguishing distinct methylation patterns (‘identical’, ‘uniform’ and ‘disordered’) compared with other metrics. Applying MeConcord to the whole genome data across 25 cell lines or primary cells or tissues, we found that distinct methylation patterns were associated with different genomic characteristics, such as CTCF binding or imprinted genes. Further, we showed the differences of CpG island hypermethylation patterns between senescence and tumorigenesis by using MeConcord. MeConcord is a powerful method to study local read-level methylation patterns for both the whole genome and specific regions of interest. Availability and implementation MeConcord is available at https://github.com/WangLabTHU/MeConcord. Supplementary information Supplementary data are available at Bioinformatics online.

Early detection is crucial to improve breast cancer (BC) patients’ outcomes and survival. Mammogram and ultrasound adopting the Breast Imaging Reporting and Data System (BI-RADS) categorization are widely used for BC early detection, while suffering high false-positive rate leading to unnecessary biopsy, especially in BI-RADS category-4 patients. Plasma cell-free DNA (cfDNA) carrying on DNA methylation information has emerged as a non-invasive approach for cancer detection. Here we present a prospective multi-center study with whole-genome bisulfite sequencing data to address the clinical utility of cfDNA methylation markers from 203 female patients with breast lesions suspected for malignancy. The cfDNA is enriched with hypo-methylated genomic regions. A practical computational framework was devised to excavate optimal cfDNA-rich DNA methylation markers, which significantly improved the early diagnosis of BI-RADS category-4 patients (AUC from 0.78–0.79 to 0.93–0.94). As a proof-of-concept study, we performed the first blood-based whole-genome DNA methylation study for detecting early-stage breast cancer from benign tumors at single-base resolution, which suggests that combining the liquid biopsy with the traditional diagnostic imaging can improve the current clinical practice, by reducing the false-positive rate and avoiding unnecessary harms. Supplementary Information The online version contains supplementary material available at 10.1186/s12943-021-01330-w.

Hematopoietic stem cell (HSC) aging correlates with an increasing risk of myeloproliferative disease and immunosenescence. In this study, we show that aging-related inflammation promotes HSC aging through tumor necrosis factor-α (TNF-α)→ERK→ETS1→interleukin27Ra (IL27Ra) pathway. TNF-α, a well-known biomarker of inflammation, increases during aging and induces the expression of IL27Ra on HSCs via ERK-ETS1 signaling. Deletion of IL27Ra rescues the functional decline and myeloid bias of HSCs and also reverses the inhibitory effect of TNF-α on HSCs. Aged IL27Ra−/− mice had a reduced proportion of myeloid-biased HSCs and did not display the biased myeloid differentiation that occurs in aged wild-type mice. IL27Ra+ HSCs exhibit impaired reconstitution capacity and myeloid-bias compared with IL27Ra− HSCs and serve as a myeloid-recovery pool upon inflammatory insult. Inflammation-related genes were enriched in IL27Ra+ HSCs and this enrichment increases with aging. Our study demonstrates that age-induced IL27Ra signaling impairs HSCs and raises the possibility that interfering with IL27Ra signaling can counter the physiologically deleterious effect of aging on hematopoietic capacity.

SUMMARYSingle-cell omics data can characterize multifaceted features of massive cells and bring significant insights to biomedical researches. The accumulation of single-cell data provides growing resources for constructing atlases for all cells of a human organ or the whole body. The true assembly of a cell atlas should be cell-centric rather than file-centric. We proposed a unified information framework enabling seamless cell-centric data assembly and developed a human Ensemble Cell Atlas (hECA) as an instance. hECA version 1.0 assembled scRNA-seq data across multiple studies into one orchestrated data repository. It contains 1,093,299 labeled cells and metadata from 116 published human single-cell studies, covering 38 human organs and 11 systems. We invented three methods of applications based on the cell-centric assembly: “In data” cell sorting enables targeted data retrieval in the full atlas with customizable logic expressions; The “quantitative portraiture” system provides a multi-view presentation of biological entities (organs, cell types, and genes) of multiple granularities; The customizable reference creation allows users to use the cell-centric assembly to generate references for their own cell type annotations. Case studies on agile construction of user-defined sub-atlases and “in data” investigation of CAR-T off-targets in multiple organs showed the great potential of cell-centric atlas assembly.

Functional genetic elements are one of the most essential units for synthetic biology. However, both knowledge-driven and data-driven methodology can hardly accomplish the complicated task of genetic elements design efficiently due to the lack of explicit regulatory logics and training samples. Here, we proposed a knowledge-constraint deep learning model named PccGEO to automatically design functional genetic elements with high success rate and efficiency. PccGEO utilized a novel “fill-in-the-flank” strategy with a conditional generative adversarial network structure to optimize the flanking regions of known functional sequences derived from the biological prior knowledge, which can efficiently capture the implicit patterns with a reduced searching space. We applied PccGEO in the design of Escherichia coli promoters, and found that the implicit patterns in flanking regions matter to the properties of promoters such as the expression level. The PccGEO-designed constitutive and inducible promoters showed more than 91.6% chance of success by in vivo validation. We further utilized PccGEO by setting a limited frequency of nucleotide modifications and surprisingly found that the expression level of E. coli sigma 70 promoters could show up to a 159.3-fold increase with only 10-bp nucleotide modifications. The results supported that the implicit patterns are important in the design of functional gene elements and validated the strong capacity of our method in the efficient design of functional genetic elements.Availabilityhttps://github.com/WangLabTHU/PccGEO

Hematopoietic stem cells (HSC) give rise to the cells of the blood system over the whole lifespan. N6-methyladenosine (m6A), the most prevalent RNA modification, modulates gene expression via the processes of "writing" and "reading". Recent studies showed that m6A "writer" genes (Mettl3 and Mettl14) play an essential role in HSC. However, which reader deciphers the m6A modification to modulate HSC remains unknown. In this study, we observed that dysfunction of Ythdf3 and Ccnd1 severely impaired the reconstitution capacity of HSC, which phenocopies Mettl3-deficient HSC. Dysfunction of Ythdf3 and Mettl3 results in a translational defect of Ccnd1. Ythdf3 and Mettl3 regulate HSC by transmitting m6A RNA methylation on the 5' untranslated region of Ccnd1. Enforced Ccnd1 expression completely rescued the defect of Ythdf3-/- HSC and partially rescued Mettl3-compromised HSC. Taken together, this study identified, for the first time, that Ccnd1 is the target of METTL3 and YTHDF3 to transmit the m6A RNA methylation signal and thereby regulate the reconstitution capacity of HSC.

Detecting cancer signals in cell-free DNA (cfDNA) high-throughput sequencing data is emerging as a novel noninvasive cancer detection method. Due to the high cost of sequencing, it is crucial to make robust and precise predictions with low-depth cfDNA sequencing data. Here we propose a novel approach named DISMIR, which can provide ultrasensitive and robust cancer detection by integrating DNA sequence and methylation information in plasma cfDNA whole-genome bisulfite sequencing (WGBS) data. DISMIR introduces a new feature termed as ‘switching region’ to define cancer-specific differentially methylated regions, which can enrich the cancer-related signal at read-resolution. DISMIR applies a deep learning model to predict the source of every single read based on its DNA sequence and methylation state and then predicts the risk that the plasma donor is suffering from cancer. DISMIR exhibited high accuracy and robustness on hepatocellular carcinoma detection by plasma cfDNA WGBS data even at ultralow sequencing depths. Further analysis showed that DISMIR tends to be insensitive to alterations of single CpG sites’ methylation states, which suggests DISMIR could resist to technical noise of WGBS. All these results showed DISMIR with the potential to be a precise and robust method for low-cost early cancer detection.

Background: Small nuclear ribonucleoprotein polypeptide A (SNRPA1) is a splicing factor (SF) responsible for the processing of pre−mRNA into mRNA. The purpose of this study was to explore the clinicopathological characteristics and prognostic significance of SNRPA1 mRNA expression in liver cancer and its potential mechanism as an SF.Methods: A total of 418 RNA−Seq and clinical data were downloaded from The Cancer Genome Atlas (TCGA) database. Alternative splicing (AS) profiles were downloaded from TCGA SpliceSeq. Wilcoxon rank-sum tests were used to compare normal tissues with tumor tissues and analyze subgroups. We used the Kaplan−Meier analysis method to draw the survival curves. Univariate and multivariate Cox analyses were employed to estimate the prognostic value of SNRPA1. Gene set enrichment analysis (GSEA) was performed to identify the signaling pathways. Then we used univariate and Pearson's correlation tests to analyze the correlation between SFs and exon skip (ES) events. Wilcoxon rank-sum tests were applied to analyze the relationships between different spliceosome and cancers. We also analyzed the relationship between the expression of SCP2 (related to SNRPA1) and clinical symptoms. We then evaluated the expression levels of these genes with clinical samples and The Clinical Proteomic Tumor Analysis Consortium (CPTAC). Results: The level of SNRPA1 mRNA expression in liver cancer was significantly up−regulated in tumor tissues compared with normal tissues (p = 1.411e−27) in liver cancer and was positively correlated with survival status (p = 0.035). In addition, we found that SNRPA1 mRNA levels can reflect the prognosis of liver cancer (hazard ratio [HR] = 1.08, 95% confidence interval [CI]: 1.02–1.14, p = 0.005). The enriched KEGG pathway by GESA revealed a splicesome as the main pathway of SNRPA1. Alternative splicing events of SNRPA1 was related to the expression of SCP2. SNRPA1 exon 6 skip and SCP2 exon 12 skip correlated with many cancer types. Finally, SNRPA1 and SCP2 percent−spliced−in (PSI) values positively correlated with survival status (p = 3.022e−04 and p = 2.932e−03). Conclusions: Splicing Factor SNRPA1 is superior to its mRNA expression in predicting the prognosis of liver cancer and correlates with the SCP2 spliceosome, which is consistent with that of SNRPA1 in clinical samples.

Quantifying cell proportions, especially for rare cell types in some scenarios, is of great value in tracking signals associated with certain phenotypes or diseases. Although some methods have been proposed to infer cell proportions from multicomponent bulk data, they are substantially less effective for estimating the proportions of rare cell types which are highly sensitive to feature outliers and collinearity. Here we proposed a new deconvolution algorithm named ARIC to estimate cell type proportions from gene expression or DNA methylation data. ARIC employs a novel two-step marker selection strategy, including collinear feature elimination based on the component-wise condition number and adaptive removal of outlier markers. This strategy can systematically obtain effective markers for weighted $\upsilon$-support vector regression to ensure a robust and precise rare proportion prediction. We showed that ARIC can accurately estimate fractions in both DNA methylation and gene expression data from different experiments. We further applied ARIC to the survival prediction of ovarian cancer and the condition monitoring of chronic kidney disease, and the results demonstrate the high accuracy and robustness as well as clinical potentials of ARIC. Taken together, ARIC is a promising tool to solve the deconvolution problem of bulk data where rare components are of vital importance.

Cellular senescence is an important mechanism of autonomous tumor suppression, while its consequence such as the senescence-associated secretory phenotype (SASP) may drive tumorigenesis and age-related diseases. Therefore, controlling the cell fate optimally when encountering senescence stress is helpful for anti-cancer or anti-aging treatments. To identify genes essential for senescence establishment or maintenance, we carried out a CRISPR-based screen with a deliberately designed single-guide RNA (sgRNA) library. The library comprised of about 12,000 kinds of sgRNAs targeting 1378 senescence-associated genes selected by integrating the information of literature mining, protein-protein interaction network, and differential gene expression. We successfully detected a dozen gene deficiencies potentially causing senescence bypass, and their phenotypes were further validated with a high true positive rate. RNA-seq analysis showed distinct transcriptome patterns of these bypass cells. Interestingly, in the bypass cells, the expression of SASP genes was maintained or elevated with CHEK2, HAS1, or MDK deficiency; but neutralized with MTOR, CRISPLD2, or MORF4L1 deficiency. Pathways of some age-related neurodegenerative disorders were also downregulated with MTOR, CRISPLD2, or MORF4L1 deficiency. The results demonstrated that disturbing these genes could lead to distinct cell fates as a consequence of senescence bypass, suggesting that they may play essential roles in cellular senescence.