Systematic mapping of organism-scale gene-regulatory networks in aging using population asynchrony.

In aging, physiologic networks decline in function at rates that differ between individuals, producing a wide distribution of lifespan. Though 70% of human lifespan variance remains unexplained by heritable factors, little is known about the intrinsic sources of physiologic heterogeneity in aging. To understand how complex physiologic networks generate lifespan variation, new methods are needed. Here, we present Asynch-seq, an approach that uses gene-expression heterogeneity within isogenic populations to study the processes generating lifespan variation. By collecting thousands of single-individual transcriptomes, we capture the Caenorhabditis elegans "pan-transcriptome"-a highly resolved atlas of non-genetic variation. We use our atlas to guide a large-scale perturbation screen that identifies the decoupling of total mRNA content between germline and soma as the largest source of physiologic heterogeneity in aging, driven by pleiotropic genes whose knockdown dramatically reduces lifespan variance. Our work demonstrates how systematic mapping of physiologic heterogeneity can be applied to reduce inter-individual disparities in aging.

An atlas of cells in the human tonsil.

Palatine tonsils are secondary lymphoid organs (SLOs) representing the first line of immunological defense against inhaled or ingested pathogens. We generated an atlas of the human tonsil composed of >556,000 cells profiled across five different data modalities, including single-cell transcriptome, epigenome, proteome, and immune repertoire sequencing, as well as spatial transcriptomics. This census identified 121 cell types and states, defined developmental trajectories, and enabled an understanding of the functional units of the tonsil. Exemplarily, we stratified myeloid slan-like subtypes, established a BCL6 enhancer as locally active in follicle-associated T and B cells, and identified SIX5 as putative transcriptional regulator of plasma cell maturation. Analyses of a validation cohort confirmed the presence, annotation, and markers of tonsillar cell types and provided evidence of age-related compositional shifts. We demonstrate the value of this resource by annotating cells from B cell-derived mantle cell lymphomas, linking transcriptional heterogeneity to normal B cell differentiation states of the human tonsil.

Identity Noise and Adipogenic Traits Characterize Dermal Fibroblast Aging.

During aging, stromal functions are thought to be impaired, but little is known whether this stems from changes of fibroblasts. Using population- and single-cell transcriptomics, as well as long-term lineage tracing, we studied whether murine dermal fibroblasts are altered during physiological aging under different dietary regimes that affect longevity. We show that the identity of old fibroblasts becomes undefined, with the fibroblast states present in young skin no longer clearly demarcated. In addition, old fibroblasts not only reduce the expression of genes involved in the formation of the extracellular matrix, but also gain adipogenic traits, paradoxically becoming more similar to neonatal pro-adipogenic fibroblasts. These alterations are sensitive to systemic metabolic changes: long-term caloric restriction reversibly prevents them, whereas a high-fat diet potentiates them. Our results therefore highlight loss of cell identity and the acquisition of adipogenic traits as a mechanism underlying cellular aging, which is influenced by systemic metabolism.

Machine Learning Identifies Stemness Features Associated with Oncogenic Dedifferentiation.

Cancer progression involves the gradual loss of a differentiated phenotype and acquisition of progenitor and stem-cell-like features. Here, we provide novel stemness indices for assessing the degree of oncogenic dedifferentiation. We used an innovative one-class logistic regression (OCLR) machine-learning algorithm to extract transcriptomic and epigenetic feature sets derived from non-transformed pluripotent stem cells and their differentiated progeny. Using OCLR, we were able to identify previously undiscovered biological mechanisms associated with the dedifferentiated oncogenic state. Analyses of the tumor microenvironment revealed unanticipated correlation of cancer stemness with immune checkpoint expression and infiltrating immune cells. We found that the dedifferentiated oncogenic phenotype was generally most prominent in metastatic tumors. Application of our stemness indices to single-cell data revealed patterns of intra-tumor molecular heterogeneity. Finally, the indices allowed for the identification of novel targets and possible targeted therapies aimed at tumor differentiation.

A Landscape of Pharmacogenomic Interactions in Cancer.

Systematic studies of cancer genomes have provided unprecedented insights into the molecular nature of cancer. Using this information to guide the development and application of therapies in the clinic is challenging. Here, we report how cancer-driven alterations identified in 11,289 tumors from 29 tissues (integrating somatic mutations, copy number alterations, DNA methylation, and gene expression) can be mapped onto 1,001 molecularly annotated human cancer cell lines and correlated with sensitivity to 265 drugs. We find that cell lines faithfully recapitulate oncogenic alterations identified in tumors, find that many of these associate with drug sensitivity/resistance, and highlight the importance of tissue lineage in mediating drug response. Logic-based modeling uncovers combinations of alterations that sensitize to drugs, while machine learning demonstrates the relative importance of different data types in predicting drug response. Our analysis and datasets are rich resources to link genotypes with cellular phenotypes and to identify therapeutic options for selected cancer sub-populations.

An Adenine Code for DNA: A Second Life for N6-Methyladenine.

DNA N6-methyladenine (6mA) protects against restriction enzymes in bacteria. However, isolated reports have suggested additional activities and its presence in other organisms, such as unicellular eukaryotes. New data now find that 6mA may have a gene regulatory function in green alga, worm, and fly, suggesting m6A as a potential "epigenetic" mark.

The emerging landscape of spatial profiling technologies.

Improved scale, multiplexing and resolution are establishing spatial nucleic acid and protein profiling methods as a major pillar for cellular atlas building of complex samples, from tissues to full organisms. Emerging methods yield omics measurements at resolutions covering the nano- to microscale, enabling the charting of cellular heterogeneity, complex tissue architectures and dynamic changes during development and disease. We present an overview of the developing landscape of in situ spatial genome, transcriptome and proteome technologies, exemplify their impact on cell biology and translational research, and discuss current challenges for their community-wide adoption. Among many transformative applications, we envision that spatial methods will map entire organs and enable next-generation pathology.

Metastatic recurrence in colorectal cancer arises from residual EMP1+ cells.

Around 30-40% of patients with colorectal cancer (CRC) undergoing curative resection of the primary tumour will develop metastases in the subsequent years1. Therapies to prevent disease relapse remain an unmet medical need. Here we uncover the identity and features of the residual tumour cells responsible for CRC relapse. An analysis of single-cell transcriptomes of samples from patients with CRC revealed that the majority of genes associated with a poor prognosis are expressed by a unique tumour cell population that we named high-relapse cells (HRCs). We established a human-like mouse model of microsatellite-stable CRC that undergoes metastatic relapse after surgical resection of the primary tumour. Residual HRCs occult in mouse livers after primary CRC surgery gave rise to multiple cell types over time, including LGR5+ stem-like tumour cells2-4, and caused overt metastatic disease. Using Emp1 (encoding epithelial membrane protein 1) as a marker gene for HRCs, we tracked and selectively eliminated this cell population. Genetic ablation of EMP1high cells prevented metastatic recurrence and mice remained disease-free after surgery. We also found that HRC-rich micrometastases were infiltrated with T cells, yet became progressively immune-excluded during outgrowth. Treatment with neoadjuvant immunotherapy eliminated residual metastatic cells and prevented mice from relapsing after surgery. Together, our findings reveal the cell-state dynamics of residual disease in CRC and anticipate that therapies targeting HRCs may help to avoid metastatic relapse.

Dietary palmitic acid promotes a prometastatic memory via Schwann cells.

Fatty acid uptake and altered metabolism constitute hallmarks of metastasis1,2, yet evidence of the underlying biology, as well as whether all dietary fatty acids are prometastatic, is lacking. Here we show that dietary palmitic acid (PA), but not oleic acid or linoleic acid, promotes metastasis in oral carcinomas and melanoma in mice. Tumours from mice that were fed a short-term palm-oil-rich diet (PA), or tumour cells that were briefly exposed to PA in vitro, remained highly metastatic even after being serially transplanted (without further exposure to high levels of PA). This PA-induced prometastatic memory requires the fatty acid transporter CD36 and is associated with the stable deposition of histone H3 lysine 4 trimethylation by the methyltransferase Set1A (as part of the COMPASS complex (Set1A/COMPASS)). Bulk, single-cell and positional RNA-sequencing analyses indicate that genes with this prometastatic memory predominantly relate to a neural signature that stimulates intratumoural Schwann cells and innervation, two parameters that are strongly correlated with metastasis but are aetiologically poorly understood3,4. Mechanistically, tumour-associated Schwann cells secrete a specialized proregenerative extracellular matrix, the ablation of which inhibits metastasis initiation. Both the PA-induced memory of this proneural signature and its long-term boost in metastasis require the transcription factor EGR2 and the glial-cell-stimulating peptide galanin. In summary, we provide evidence that a dietary metabolite induces stable transcriptional and chromatin changes that lead to a long-term stimulation of metastasis, and that this is related to a proregenerative state of tumour-activated Schwann cells.

Controlled X-chromosome dynamics defines meiotic potential of female mouse in vitro germ cells.

The mammalian germline is characterized by extensive epigenetic reprogramming during its development into functional eggs and sperm. Specifically, the epigenome requires resetting before parental marks can be established and transmitted to the next generation. In the female germline, X-chromosome inactivation and reactivation are among the most prominent epigenetic reprogramming events, yet very little is known about their kinetics and biological function. Here, we investigate X-inactivation and reactivation dynamics using a tailor-made in vitro system of primordial germ cell-like cell (PGCLC) differentiation from mouse embryonic stem cells. We find that X-inactivation in PGCLCs in vitro and in germ cell-competent epiblast cells in vivo is moderate compared to somatic cells, and frequently characterized by escaping genes. X-inactivation is followed by step-wise X-reactivation, which is mostly completed during meiotic prophase I. Furthermore, we find that PGCLCs which fail to undergo X-inactivation or reactivate too rapidly display impaired meiotic potential. Thus, our data reveal fine-tuned X-chromosome remodelling as a critical feature of female germ cell development towards meiosis and oogenesis.

Single-cell multi-omics analysis of COVID-19 patients with pre-existing autoimmune diseases shows aberrant immune responses to infection.

In COVID-19, hyperinflammatory and dysregulated immune responses contribute to severity. Patients with pre-existing autoimmune conditions can therefore be at increased risk of severe COVID-19 and/or associated sequelae, yet SARS-CoV-2 infection in this group has been little studied. Here, we performed single-cell analysis of peripheral blood mononuclear cells from patients with three major autoimmune diseases (rheumatoid arthritis, psoriasis, or multiple sclerosis) during SARS-CoV-2 infection. We observed compositional differences between the autoimmune disease groups coupled with altered patterns of gene expression, transcription factor activity, and cell-cell communication that substantially shape the immune response under SARS-CoV-2 infection. While enrichment of HLA-DRlow CD14+ monocytes was observed in all three autoimmune disease groups, type-I interferon signaling as well as inflammatory T cell and monocyte responses varied widely between the three groups of patients. Our results reveal disturbed immune responses to SARS-CoV-2 in patients with pre-existing autoimmunity, highlighting important considerations for disease treatment and follow-up.

Comparative Analysis of Single-Cell RNA Sequencing Methods.

Single-cell RNA sequencing (scRNA-seq) offers new possibilities to address biological and medical questions. However, systematic comparisons of the performance of diverse scRNA-seq protocols are lacking. We generated data from 583 mouse embryonic stem cells to evaluate six prominent scRNA-seq methods: CEL-seq2, Drop-seq, MARS-seq, SCRB-seq, Smart-seq, and Smart-seq2. While Smart-seq2 detected the most genes per cell and across cells, CEL-seq2, Drop-seq, MARS-seq, and SCRB-seq quantified mRNA levels with less amplification noise due to the use of unique molecular identifiers (UMIs). Power simulations at different sequencing depths showed that Drop-seq is more cost-efficient for transcriptome quantification of large numbers of cells, while MARS-seq, SCRB-seq, and Smart-seq2 are more efficient when analyzing fewer cells. Our quantitative comparison offers the basis for an informed choice among six prominent scRNA-seq methods, and it provides a framework for benchmarking further improvements of scRNA-seq protocols.

A single-cell tumor immune atlas for precision oncology.

The tumor immune microenvironment is a main contributor to cancer progression and a promising therapeutic target for oncology. However, immune microenvironments vary profoundly between patients, and biomarkers for prognosis and treatment response lack precision. A comprehensive compendium of tumor immune cells is required to pinpoint predictive cellular states and their spatial localization. We generated a single-cell tumor immune atlas, jointly analyzing published data sets of >500,000 cells from 217 patients and 13 cancer types, providing the basis for a patient stratification based on immune cell compositions. Projecting immune cells from external tumors onto the atlas facilitated an automated cell annotation system. To enable in situ mapping of immune populations for digital pathology, we applied SPOTlight, combining single-cell and spatial transcriptomics data and identifying colocalization patterns of immune, stromal, and cancer cells in tumor sections. We expect the tumor immune cell atlas, together with our versatile toolbox for precision oncology, to advance currently applied stratification approaches for prognosis and immunotherapy.

The cellular landscape of the endochondral bone during the transition to extrauterine life.

The cellular complexity of the endochondral bone underlies its essential and pleiotropic roles during organismal life. While the adult bone has received significant attention, we still lack a deep understanding of the perinatal bone cellulome. Here, we have profiled the full composition of the murine endochondral bone at the single-cell level during the transition from fetal to newborn life and in comparison with the adult tissue, with particular emphasis on the mesenchymal compartment. The perinatal bone contains different fibroblastic clusters with blastema-like characteristics in organizing and supporting skeletogenesis, angiogenesis and hematopoiesis. Our data also suggest dynamic inter- and intra-compartment interactions, as well as a bone marrow milieu that seems prone to anti-inflammation, which we hypothesize is necessary to ensure the proper program of lymphopoiesis and the establishment of central and peripheral tolerance in early life. Our study provides an integrative roadmap for the future design of genetic and cellular functional assays to validate cellular interactions and lineage relationships within the perinatal bone.

SPOTlight: seeded NMF regression to deconvolute spatial transcriptomics spots with single-cell transcriptomes.

Spatially resolved gene expression profiles are key to understand tissue organization and function. However, spatial transcriptomics (ST) profiling techniques lack single-cell resolution and require a combination with single-cell RNA sequencing (scRNA-seq) information to deconvolute the spatially indexed datasets. Leveraging the strengths of both data types, we developed SPOTlight, a computational tool that enables the integration of ST with scRNA-seq data to infer the location of cell types and states within a complex tissue. SPOTlight is centered around a seeded non-negative matrix factorization (NMF) regression, initialized using cell-type marker genes and non-negative least squares (NNLS) to subsequently deconvolute ST capture locations (spots). Simulating varying reference quantities and qualities, we confirmed high prediction accuracy also with shallowly sequenced or small-sized scRNA-seq reference datasets. SPOTlight deconvolution of the mouse brain correctly mapped subtle neuronal cell states of the cortical layers and the defined architecture of the hippocampus. In human pancreatic cancer, we successfully segmented patient sections and further fine-mapped normal and neoplastic cell states. Trained on an external single-cell pancreatic tumor references, we further charted the localization of clinical-relevant and tumor-specific immune cell states, an illustrative example of its flexible application spectrum and future potential in digital pathology.

bigSCale: an analytical framework for big-scale single-cell data.

Single-cell RNA sequencing (scRNA-seq) has significantly deepened our insights into complex tissues, with the latest techniques capable of processing tens of thousands of cells simultaneously. Analyzing increasing numbers of cells, however, generates extremely large data sets, extending processing time and challenging computing resources. Current scRNA-seq analysis tools are not designed to interrogate large data sets and often lack sensitivity to identify marker genes. With bigSCale, we provide a scalable analytical framework to analyze millions of cells, which addresses the challenges associated with large data sets. To handle the noise and sparsity of scRNA-seq data, bigSCale uses large sample sizes to estimate an accurate numerical model of noise. The framework further includes modules for differential expression analysis, cell clustering, and marker identification. A directed convolution strategy allows processing of extremely large data sets, while preserving transcript information from individual cells. We evaluated the performance of bigSCale using both a biological model of aberrant gene expression in patient-derived neuronal progenitor cells and simulated data sets, which underlines the speed and accuracy in differential expression analysis. To test its applicability for large data sets, we applied bigSCale to assess 1.3 million cells from the mouse developing forebrain. Its directed down-sampling strategy accumulates information from single cells into index cell transcriptomes, thereby defining cellular clusters with improved resolution. Accordingly, index cell clusters identified rare populations, such as reelin (Reln)-positive Cajal-Retzius neurons, for which we report previously unrecognized heterogeneity associated with distinct differentiation stages, spatial organization, and cellular function. Together, bigSCale presents a solution to address future challenges of large single-cell data sets.

Quantitative Trait Loci Identify Functional Noncoding Variation in Cancer.

The interpretation of noncoding alterations in cancer genomes presents an unresolved problem in cancer studies. While the impact of somatic variations in protein-coding regions is widely accepted, noncoding aberrations are mostly considered as passenger events. However, with the advance of genome-wide profiling strategies, alterations outside the coding context entered the focus, and multiple examples highlight the role of gene deregulation as cancer-driving events. This review describes the implication of noncoding alterations in oncogenesis and provides a theoretical framework for the identification of causal somatic variants using quantitative trait loci (QTL) analysis. Assuming that functional noncoding alterations affect quantifiable regulatory processes, somatic QTL studies constitute a valuable strategy to pinpoint cancer gene deregulation. Eventually, the comprehensive identification and interpretation of coding and noncoding alterations will guide our future understanding of cancer biology.

Epigenetic inactivation of the p53-induced long noncoding RNA TP53 target 1 in human cancer.

Long noncoding RNAs (lncRNAs) are important regulators of cellular homeostasis. However, their contribution to the cancer phenotype still needs to be established. Herein, we have identified a p53-induced lncRNA, TP53TG1, that undergoes cancer-specific promoter hypermethylation-associated silencing. In vitro and in vivo assays identify a tumor-suppressor activity for TP53TG1 and a role in the p53 response to DNA damage. Importantly, we show that TP53TG1 binds to the multifaceted DNA/RNA binding protein YBX1 to prevent its nuclear localization and thus the YBX1-mediated activation of oncogenes. TP53TG1 epigenetic inactivation in cancer cells releases the transcriptional repression of YBX1-targeted growth-promoting genes and creates a chemoresistant tumor. TP53TG1 hypermethylation in primary tumors is shown to be associated with poor outcome. The epigenetic loss of TP53TG1 therefore represents an altered event in an lncRNA that is linked to classical tumoral pathways, such as p53 signaling, but is also connected to regulatory networks of the cancer cell.

Metformin pharmacogenomics: a genome-wide association study to identify genetic and epigenetic biomarkers involved in metformin anticancer response using human lymphoblastoid cell lines.

Metformin is currently considered as a promising anticancer agent in addition to its anti-diabetic effect. To better individualize metformin therapy and explore novel molecular mechanisms in cancer treatment, we conducted a pharmacogenomic study using 266 lymphoblastoid cell lines (LCLs). Metformin cytotoxicity assay was performed using the MTS assay. Genome-wide association (GWA) analyses were performed in LCLs using 1.3 million SNPs, 485k DNA methylation probes, 54k mRNA expression probe sets, and metformin cytotoxicity (IC50s). Top candidate genes were functionally validated using siRNA screening, followed by MTS assay in breast cancer cell lines. Further study of one top candidate, STUB1, was performed to elucidate the mechanisms by which STUB1 might contribute to metformin action. GWA analyses in LCLs identified 198 mRNA expression probe sets, 12 SNP loci, and 5 DNA methylation loci associated with metformin IC50 with P-values <10−4 or <10−5. Integrated SNP/methylation loci-expression-IC50 analyses found 3 SNP loci or 5 DNA methylation loci associated with metformin IC50 through trans-regulation of expression of 11 or 26 genes with P-value <10−4. Functional validation of top 61 candidate genes in 4 IPA networks indicated down regulation of 14 genes significantly altered metformin sensitivity in two breast cancer cell lines. Mechanistic studies revealed that the E3 ubiquitin ligase, STUB1, could influence metformin response by facilitating proteasome-mediated degradation of cyclin A. GWAS using a genomic data-enriched LCL model system, together with functional and mechanistic studies using cancer cell lines, help us to identify novel genetic and epigenetic biomarkers involved in metformin anticancer response.

DNA methylation profiling in the clinic: applications and challenges.

Knowledge of epigenetic alterations in disease is rapidly increasing owing to the development of genome-wide techniques for their identification. The ever-growing number of genes that show epigenetic alterations in disease emphasizes the crucial role of these epigenetic alterations - particularly DNA methylation - for future diagnosis, prognosis and prediction of response to therapies. This Review focuses on epigenetic profiling, which has started to be of clinical value in cancer and may in the future be extended to other diseases, such as neurological and autoimmune disorders.