Gene Regulatory Circuits in Innate and Adaptive Immune Cells.

Cell identity and function largely rely on the programming of transcriptomes during development and differentiation. Signature gene expression programs are orchestrated by regulatory circuits consisting of cis-acting promoters and enhancers, which respond to a plethora of cues via the action of transcription factors. In turn, transcription factors direct epigenetic modifications to revise chromatin landscapes, and drive contacts between distal promoter-enhancer combinations. In immune cells, regulatory circuits for effector genes are especially complex and flexible, utilizing distinct sets of transcription factors and enhancers, depending on the cues each cell type receives during an infection, after sensing cellular damage, or upon encountering a tumor. Here, we review major players in the coordination of gene regulatory programs within innate and adaptive immune cells, as well as integrative omics approaches that can be leveraged to decipher their underlying circuitry.

A TCF4-dependent gene regulatory network confers resistance to immunotherapy in melanoma.

To better understand intrinsic resistance to immune checkpoint blockade (ICB), we established a comprehensive view of the cellular architecture of the treatment-naive melanoma ecosystem and studied its evolution under ICB. Using single-cell, spatial multi-omics, we showed that the tumor microenvironment promotes the emergence of a complex melanoma transcriptomic landscape. Melanoma cells harboring a mesenchymal-like (MES) state, a population known to confer resistance to targeted therapy, were significantly enriched in early on-treatment biopsies from non-responders to ICB. TCF4 serves as the hub of this landscape by being a master regulator of the MES signature and a suppressor of the melanocytic and antigen presentation transcriptional programs. Targeting TCF4 genetically or pharmacologically, using a bromodomain inhibitor, increased immunogenicity and sensitivity of MES cells to ICB and targeted therapy. We thereby uncovered a TCF4-dependent regulatory network that orchestrates multiple transcriptional programs and contributes to resistance to both targeted therapy and ICB in melanoma.

Massively parallel in vivo Perturb-seq reveals cell-type-specific transcriptional networks in cortical development.

Leveraging AAVs' versatile tropism and labeling capacity, we expanded the scale of in vivo CRISPR screening with single-cell transcriptomic phenotyping across embryonic to adult brains and peripheral nervous systems. Through extensive tests of 86 vectors across AAV serotypes combined with a transposon system, we substantially amplified labeling efficacy and accelerated in vivo gene delivery from weeks to days. Our proof-of-principle in utero screen identified the pleiotropic effects of Foxg1, highlighting its tight regulation of distinct networks essential for cell fate specification of Layer 6 corticothalamic neurons. Notably, our platform can label >6% of cerebral cells, surpassing the current state-of-the-art efficacy at <0.1% by lentivirus, to achieve analysis of over 30,000 cells in one experiment and enable massively parallel in vivo Perturb-seq. Compatible with various phenotypic measurements (single-cell or spatial multi-omics), it presents a flexible approach to interrogate gene function across cell types in vivo, translating gene variants to their causal function.

Transcriptional Control of Dendritic Cell Development.

The dendritic cells (DCs) of the immune system function in innate and adaptive responses by directing activity of various effector cells rather than serving as effectors themselves. DCs and closely related myeloid lineages share expression of many surface receptors, presenting a challenge in distinguishing their unique in vivo functions. Recent work has taken advantage of unique transcriptional programs to identify and manipulate murine DCs in vivo. This work has assigned several nonredundant in vivo functions to distinct DC lineages, consisting of plasmacytoid DCs and several subsets of classical DCs that promote different immune effector modules in response to pathogens. In parallel, a correspondence between human and murine DC subsets has emerged, underlying structural similarities for the DC lineages between these species. Recent work has begun to unravel the transcriptional circuitry that controls the development and diversification of DCs from common progenitors in the bone marrow.

The Innate Lymphoid Cell Precursor.

The discovery of tissue-resident innate lymphoid cell populations effecting different forms of type 1, 2, and 3 immunity; tissue repair; and immune regulation has transformed our understanding of mucosal immunity and allergy. The emerging complexity of these populations along with compounding issues of redundancy and plasticity raise intriguing questions about their precise lineage relationship. Here we review advances in mapping the emergence of these lineages from early lymphoid precursors. We discuss the identification of a common innate lymphoid cell precursor characterized by transient expression of the transcription factor PLZF, and the lineage relationships of innate lymphoid cells with conventional natural killer cells and lymphoid tissue inducer cells. We also review the rapidly growing understanding of the network of transcription factors that direct the development of these lineages.

Cooperative assembly confers regulatory specificity and long-term genetic circuit stability.

A ubiquitous feature of eukaryotic transcriptional regulation is cooperative self-assembly between transcription factors (TFs) and DNA cis-regulatory motifs. It is thought that this strategy enables specific regulatory connections to be formed in gene networks between otherwise weakly interacting, low-specificity molecular components. Here, using synthetic gene circuits constructed in yeast, we find that high regulatory specificity can emerge from cooperative, multivalent interactions among artificial zinc-finger-based TFs. We show that circuits "wired" using the strategy of cooperative TF assembly are effectively insulated from aberrant misregulation of the host cell genome. As we demonstrate in experiments and mathematical models, this mechanism is sufficient to rescue circuit-driven fitness defects, resulting in genetic and functional stability of circuits in long-term continuous culture. Our naturally inspired approach offers a simple, generalizable means for building high-fidelity, evolutionarily robust gene circuits that can be scaled to a wide range of host organisms and applications.

Human prefrontal cortex gene regulatory dynamics from gestation to adulthood at single-cell resolution.

Human brain development is underpinned by cellular and molecular reconfigurations continuing into the third decade of life. To reveal cell dynamics orchestrating neural maturation, we profiled human prefrontal cortex gene expression and chromatin accessibility at single-cell resolution from gestation to adulthood. Integrative analyses define the dynamic trajectories of each cell type, revealing major gene expression reconfiguration at the prenatal-to-postnatal transition in all cell types followed by continuous reconfiguration into adulthood and identifying regulatory networks guiding cellular developmental programs, states, and functions. We uncover links between expression dynamics and developmental milestones, characterize the diverse timing of when cells acquire adult-like states, and identify molecular convergence from distinct developmental origins. We further reveal cellular dynamics and their regulators implicated in neurological disorders. Finally, using this reference, we benchmark cell identities and maturation states in organoid models. Together, this captures the dynamic regulatory landscape of human cortical development.

The Chinese pine genome and methylome unveil key features of conifer evolution.

Conifers dominate the world's forest ecosystems and are the most widely planted tree species. Their giant and complex genomes present great challenges for assembling a complete reference genome for evolutionary and genomic studies. We present a 25.4-Gb chromosome-level assembly of Chinese pine (Pinus tabuliformis) and revealed that its genome size is mostly attributable to huge intergenic regions and long introns with high transposable element (TE) content. Large genes with long introns exhibited higher expressions levels. Despite a lack of recent whole-genome duplication, 91.2% of genes were duplicated through dispersed duplication, and expanded gene families are mainly related to stress responses, which may underpin conifers' adaptation, particularly in cold and/or arid conditions. The reproductive regulation network is distinct compared with angiosperms. Slow removal of TEs with high-level methylation may have contributed to genomic expansion. This study provides insights into conifer evolution and resources for advancing research on conifer adaptation and development.

Disrupting autorepression circuitry generates "open-loop lethality" to yield escape-resistant antiviral agents.

Across biological scales, gene-regulatory networks employ autorepression (negative feedback) to maintain homeostasis and minimize failure from aberrant expression. Here, we present a proof of concept that disrupting transcriptional negative feedback dysregulates viral gene expression to therapeutically inhibit replication and confers a high evolutionary barrier to resistance. We find that nucleic-acid decoys mimicking cis-regulatory sites act as "feedback disruptors," break homeostasis, and increase viral transcription factors to cytotoxic levels (termed "open-loop lethality"). Feedback disruptors against herpesviruses reduced viral replication >2-logs without activating innate immunity, showed sub-nM IC50, synergized with standard-of-care antivirals, and inhibited virus replication in mice. In contrast to approved antivirals where resistance rapidly emerged, no feedback-disruptor escape mutants evolved in long-term cultures. For SARS-CoV-2, disruption of a putative feedback circuit also generated open-loop lethality, reducing viral titers by >1-log. These results demonstrate that generating open-loop lethality, via negative-feedback disruption, may yield a class of antimicrobials with a high genetic barrier to resistance.

A modular master regulator landscape controls cancer transcriptional identity.

Despite considerable efforts, the mechanisms linking genomic alterations to the transcriptional identity of cancer cells remain elusive. Integrative genomic analysis, using a network-based approach, identified 407 master regulator (MR) proteins responsible for canalizing the genetics of individual samples from 20 cohorts in The Cancer Genome Atlas (TCGA) into 112 transcriptionally distinct tumor subtypes. MR proteins could be further organized into 24 pan-cancer, master regulator block modules (MRBs), each regulating key cancer hallmarks and predictive of patient outcome in multiple cohorts. Of all somatic alterations detected in each individual sample, >50% were predicted to induce aberrant MR activity, yielding insight into mechanisms linking tumor genetics and transcriptional identity and establishing non-oncogene dependencies. Genetic and pharmacological validation assays confirmed the predicted effect of upstream mutations and MR activity on downstream cellular identity and phenotype. Thus, co-analysis of mutational and gene expression profiles identified elusive subtypes and provided testable hypothesis for mechanisms mediating the effect of genetic alterations.

A cis-regulatory atlas in maize at single-cell resolution.

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

Changes in genome architecture and transcriptional dynamics progress independently of sensory experience during post-natal brain development.

Both transcription and three-dimensional (3D) architecture of the mammalian genome play critical roles in neurodevelopment and its disorders. However, 3D genome structures of single brain cells have not been solved; little is known about the dynamics of single-cell transcriptome and 3D genome after birth. Here, we generated a transcriptome (3,517 cells) and 3D genome (3,646 cells) atlas of the developing mouse cortex and hippocampus by using our high-resolution multiple annealing and looping-based amplification cycles for digital transcriptomics (MALBAC-DT) and diploid chromatin conformation capture (Dip-C) methods and developing multi-omic analysis pipelines. In adults, 3D genome "structure types" delineate all major cell types, with high correlation between chromatin A/B compartments and gene expression. During development, both transcriptome and 3D genome are extensively transformed in the first post-natal month. In neurons, 3D genome is rewired across scales, correlated with gene expression modules, and independent of sensory experience. Finally, we examine allele-specific structure of imprinted genes, revealing local and chromosome (chr)-wide differences. These findings uncover an unknown dimension of neurodevelopment.

Innovation, conservation, and repurposing of gene function in root cell type development.

Plant species have evolved myriads of solutions, including complex cell type development and regulation, to adapt to dynamic environments. To understand this cellular diversity, we profiled tomato root cell type translatomes. Using xylem differentiation in tomato, examples of functional innovation, repurposing, and conservation of transcription factors are described, relative to the model plant Arabidopsis. Repurposing and innovation of genes are further observed within an exodermis regulatory network and illustrate its function. Comparative translatome analyses of rice, tomato, and Arabidopsis cell populations suggest increased expression conservation of root meristems compared with other homologous populations. In addition, the functions of constitutively expressed genes are more conserved than those of cell type/tissue-enriched genes. These observations suggest that higher order properties of cell type and pan-cell type regulation are evolutionarily conserved between plants and animals.

In vivo CD8+ T cell CRISPR screening reveals control by Fli1 in infection and cancer.

Improving effector activity of antigen-specific T cells is a major goal in cancer immunotherapy. Despite the identification of several effector T cell (TEFF)-driving transcription factors (TFs), the transcriptional coordination of TEFF biology remains poorly understood. We developed an in vivo T cell CRISPR screening platform and identified a key mechanism restraining TEFF biology through the ETS family TF, Fli1. Genetic deletion of Fli1 enhanced TEFF responses without compromising memory or exhaustion precursors. Fli1 restrained TEFF lineage differentiation by binding to cis-regulatory elements of effector-associated genes. Loss of Fli1 increased chromatin accessibility at ETS:RUNX motifs, allowing more efficient Runx3-driven TEFF biology. CD8+ T cells lacking Fli1 provided substantially better protection against multiple infections and tumors. These data indicate that Fli1 safeguards the developing CD8+ T cell transcriptional landscape from excessive ETS:RUNX-driven TEFF cell differentiation. Moreover, genetic deletion of Fli1 improves TEFF differentiation and protective immunity in infections and cancer.

A phosphate starvation response-centered network regulates mycorrhizal symbiosis.

To secure phosphorus (P) from soil, most land plants use a direct phosphate uptake pathway via root hairs and epidermis and an indirect phosphate uptake pathway via mycorrhizal symbiosis. The interaction between these two pathways is unclear. Here, we mapped a network between transcription factors and mycorrhizal symbiosis-related genes using Y1H. Intriguingly, this gene regulatory network is governed by the conserved P-sensing pathway, centered on phosphate starvation response (PHR) transcription factors. PHRs are required for mycorrhizal symbiosis and regulate symbiosis-related genes via the P1BS motif. SPX-domain proteins suppress OsPHR2-mediated induction of symbiosis-related genes and inhibit mycorrhizal infection. In contrast, plants overexpressing OsPHR2 show improved mycorrhizal infection and are partially resistant to P-mediated inhibition of symbiosis. Functional analyses of network nodes revealed co-regulation of hormonal signaling and mycorrhizal symbiosis. This network deciphers extensive regulation of mycorrhizal symbiosis by endogenous and exogenous signals and highlights co-option of the P-sensing pathway for mycorrhizal symbiosis.

Differential pre-malignant programs and microenvironment chart distinct paths to malignancy in human colorectal polyps.

Colorectal cancers (CRCs) arise from precursor polyps whose cellular origins, molecular heterogeneity, and immunogenic potential may reveal diagnostic and therapeutic insights when analyzed at high resolution. We present a single-cell transcriptomic and imaging atlas of the two most common human colorectal polyps, conventional adenomas and serrated polyps, and their resulting CRC counterparts. Integrative analysis of 128 datasets from 62 participants reveals adenomas arise from WNT-driven expansion of stem cells, while serrated polyps derive from differentiated cells through gastric metaplasia. Metaplasia-associated damage is coupled to a cytotoxic immune microenvironment preceding hypermutation, driven partly by antigen-presentation differences associated with tumor cell-differentiation status. Microsatellite unstable CRCs contain distinct non-metaplastic regions where tumor cells acquire stem cell properties and cytotoxic immune cells are depleted. Our multi-omic atlas provides insights into malignant progression of colorectal polyps and their microenvironment, serving as a framework for precision surveillance and prevention of CRC.

Genome-wide CRISPR Screens Reveal Host Factors Critical for SARS-CoV-2 Infection.

Identification of host genes essential for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection may reveal novel therapeutic targets and inform our understanding of coronavirus disease 2019 (COVID-19) pathogenesis. Here we performed genome-wide CRISPR screens in Vero-E6 cells with SARS-CoV-2, Middle East respiratory syndrome CoV (MERS-CoV), bat CoV HKU5 expressing the SARS-CoV-1 spike, and vesicular stomatitis virus (VSV) expressing the SARS-CoV-2 spike. We identified known SARS-CoV-2 host factors, including the receptor ACE2 and protease Cathepsin L. We additionally discovered pro-viral genes and pathways, including HMGB1 and the SWI/SNF chromatin remodeling complex, that are SARS lineage and pan-coronavirus specific, respectively. We show that HMGB1 regulates ACE2 expression and is critical for entry of SARS-CoV-2, SARS-CoV-1, and NL63. We also show that small-molecule antagonists of identified gene products inhibited SARS-CoV-2 infection in monkey and human cells, demonstrating the conserved role of these genetic hits across species. This identifies potential therapeutic targets for SARS-CoV-2 and reveals SARS lineage-specific and pan-CoV host factors that regulate susceptibility to highly pathogenic CoVs.

Runx factors launch T cell and innate lymphoid programs via direct and gene network-based mechanisms.

Runx factors are essential for lineage specification of various hematopoietic cells, including T lymphocytes. However, they regulate context-specific genes and occupy distinct genomic regions in different cell types. Here, we show that dynamic Runx binding shifts in mouse early T cell development are mostly not restricted by local chromatin state but regulated by Runx dosage and functional partners. Runx cofactors compete to recruit a limited pool of Runx factors in early T progenitor cells, and a modest increase in Runx protein availability at pre-commitment stages causes premature Runx occupancy at post-commitment binding sites. This increased Runx factor availability results in striking T cell lineage developmental acceleration by selectively activating T cell-identity and innate lymphoid cell programs. These programs are collectively regulated by Runx together with other, Runx-induced transcription factors that co-occupy Runx-target genes and propagate gene network changes.

Metabolic regulation of T lymphocytes.

T cell activation leads to dramatic shifts in cell metabolism to protect against pathogens and to orchestrate the action of other immune cells. Quiescent T cells require predominantly ATP-generating processes, whereas proliferating effector T cells require high metabolic flux through growth-promoting pathways. Further, functionally distinct T cell subsets require distinct energetic and biosynthetic pathways to support their specific functional needs. Pathways that control immune cell function and metabolism are intimately linked, and changes in cell metabolism at both the cell and system levels have been shown to enhance or suppress specific T cell functions. As a result of these findings, cell metabolism is now appreciated as a key regulator of T cell function specification and fate. This review discusses the role of cellular metabolism in T cell development, activation, differentiation, and function to highlight the clinical relevance and opportunities for therapeutic interventions that may be used to disrupt immune pathogenesis.

Pathogenic Mechanisms of Somatic Mutation and Genome Mosaicism in Aging.

Age-related accumulation of postzygotic DNA mutations results in tissue genetic heterogeneity known as somatic mosaicism. Although implicated in aging as early as the 1950s, somatic mutations in normal tissue have been difficult to study because of their low allele fractions. With the recent emergence of cost-effective high-throughput sequencing down to the single-cell level, enormous progress has been made in our capability to quantitatively analyze somatic mutations in human tissue in relation to aging and disease. Here we first review how recent technological progress has opened up this field, providing the first broad sets of quantitative information on somatic mutations in vivo necessary to gain insight into their possible causal role in human aging and disease. We then propose three major mechanisms that can lead from accumulated de novo mutations across tissues to cell functional loss and human disease.