Scalable and unbiased sequence-informed embedding of single-cell ATAC-seq data with CellSpace.

Standard scATAC sequencing (scATAC-seq) analysis pipelines represent cells as sparse numeric vectors relative to an atlas of peaks or genomic tiles and consequently ignore genomic sequence information at accessible loci. Here we present CellSpace, an efficient and scalable sequence-informed embedding algorithm for scATAC-seq that learns a mapping of DNA k-mers and cells to the same space, to address this limitation. We show that CellSpace captures meaningful latent structure in scATAC-seq datasets, including cell subpopulations and developmental hierarchies, and can score transcription factor activities in single cells based on proximity to binding motifs embedded in the same space. Importantly, CellSpace implicitly mitigates batch effects arising from multiple samples, donors or assays, even when individual datasets are processed relative to different peak atlases. Thus, CellSpace provides a powerful tool for integrating and interpreting large-scale scATAC-seq compendia.

Discovery of Competent Chromatin Regions in Human Embryonic Stem Cells.

The mechanisms underlying the ability of embryonic stem cells (ESCs) to rapidly activate lineage-specific genes during differentiation remain largely unknown. Through multiple CRISPR-activation screens, we discovered human ESCs have pre-established transcriptionally competent chromatin regions (CCRs) that support lineage-specific gene expression at levels comparable to differentiated cells. CCRs reside in the same topological domains as their target genes. They lack typical enhancer-associated histone modifications but show enriched occupancy of pluripotent transcription factors, DNA demethylation factors, and histone deacetylases. TET1 and QSER1 protect CCRs from excessive DNA methylation, while HDAC1 family members prevent premature activation. This "push and pull" feature resembles bivalent domains at developmental gene promoters but involves distinct molecular mechanisms. Our study provides new insights into pluripotency regulation and cellular plasticity in development and disease. One sentence summary: We report a class of distal regulatory regions distinct from enhancers that confer human embryonic stem cells with the competence to rapidly activate the expression of lineage-specific genes.

Novel antigen-presenting cell imparts Treg-dependent tolerance to gut microbiota.

Establishing and maintaining tolerance to self-antigens or innocuous foreign antigens is vital for the preservation of organismal health. Within the thymus, medullary thymic epithelial cells (mTECs) expressing autoimmune regulator (AIRE) have a critical role in self-tolerance through deletion of autoreactive T cells and promotion of thymic regulatory T (Treg) cell development1-4. Within weeks of birth, a separate wave of Treg cell differentiation occurs in the periphery upon exposure to antigens derived from the diet and commensal microbiota5-8, yet the cell types responsible for the generation of peripheral Treg (pTreg) cells have not been identified. Here we describe the identification of a class of RORγt+ antigen-presenting cells called Thetis cells, with transcriptional features of both mTECs and dendritic cells, comprising four major sub-groups (TC I-TC IV). We uncover a developmental wave of Thetis cells within intestinal lymph nodes during a critical window in early life, coinciding with the wave of pTreg cell differentiation. Whereas TC I and TC III expressed the signature mTEC nuclear factor AIRE, TC IV lacked AIRE expression and was enriched for molecules required for pTreg generation, including the TGF-ß-activating integrin αvß8. Loss of either major histocompatibility complex class II (MHCII) or ITGB8 by Thetis cells led to a profound impairment in intestinal pTreg differentiation, with ensuing colitis. By contrast, MHCII expression by RORγt+ group 3 innate lymphoid cells (ILC3) and classical dendritic cells was neither sufficient nor required for pTreg generation, further implicating TC IV as the tolerogenic RORγt+ antigen-presenting cell with an essential function in early life. Our studies reveal parallel pathways for the establishment of tolerance to self and foreign antigens in the thymus and periphery, respectively, marked by the involvement of shared cellular and transcriptional programmes.

CRISPR screening uncovers a central requirement for HHEX in pancreatic lineage commitment and plasticity restriction.

The pancreas and liver arise from a common pool of progenitors. However, the underlying mechanisms that drive their lineage diversification from the foregut endoderm are not fully understood. To tackle this question, we undertook a multifactorial approach that integrated human pluripotent-stem-cell-guided differentiation, genome-scale CRISPR-Cas9 screening, single-cell analysis, genomics and proteomics. We discovered that HHEX, a transcription factor (TF) widely recognized as a key regulator of liver development, acts as a gatekeeper of pancreatic lineage specification. HHEX deletion impaired pancreatic commitment and unleashed an unexpected degree of cellular plasticity towards the liver and duodenum fates. Mechanistically, HHEX cooperates with the pioneer TFs FOXA1, FOXA2 and GATA4, shared by both pancreas and liver differentiation programmes, to promote pancreas commitment, and this cooperation restrains the shared TFs from activating alternative lineages. These findings provide a generalizable model for how gatekeeper TFs like HHEX orchestrate lineage commitment and plasticity restriction in broad developmental contexts.

Demographic and genetic factors influence the abundance of infiltrating immune cells in human tissues.

Despite infiltrating immune cells having an essential function in human disease and patients' responses to treatments, mechanisms influencing variability in infiltration patterns remain unclear. Here, using bulk RNA-seq data from 46 tissues in the Genotype-Tissue Expression project, we apply cell-type deconvolution algorithms to evaluate the immune landscape across the healthy human body. We discover that 49 of 189 infiltration-related phenotypes are associated with either age or sex (FDR < 0.1). Genetic analyses further show that 31 infiltration-related phenotypes have genome-wide significant associations (iQTLs) (P < 5.0 × 10-8), with a significant enrichment of same-tissue expression quantitative trait loci in suggested iQTLs (P < 10-5). Furthermore, we find an association between helper T cell content in thyroid tissue and a COMMD3/DNAJC1 regulatory variant (P = 7.5 × 10-10), which is associated with thyroiditis in other cohorts. Together, our results identify key factors influencing inter-individual variability of immune infiltration, to provide insights on potential therapeutic targets.