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.

Immunobiology of Long Noncoding RNAs.

The discovery of long noncoding RNAs (lncRNA) has provided a new perspective on gene regulation in diverse biological contexts. lncRNAs are remarkably versatile molecules that interact with RNA, DNA, or proteins to promote or restrain the expression of protein-coding genes. Activation of immune cells is associated with dynamic changes in expression of genes, the products of which combat infectious microorganisms, initiate repair, and resolve inflammatory responses in cells and tissues. Recent evidence indicates that lncRNAs play important roles in directing the development of diverse immune cells and controlling the dynamic transcriptional programs that are a hallmark of immune cell activation. The importance of these molecules is underscored by their newly recognized roles in inflammatory diseases. In this review, we discuss the contribution of lncRNAs in the development and activation of immune cells and their roles in immune-related diseases. We also discuss challenges faced in identifying biological functions for this large and complex class of genes.

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.

Hormonal gatekeeping via the blood-brain barrier governs caste-specific behavior in ants.

Here, we reveal an unanticipated role of the blood-brain barrier (BBB) in regulating complex social behavior in ants. Using scRNA-seq, we find localization in the BBB of a key hormone-degrading enzyme called juvenile hormone esterase (Jhe), and we show that this localization governs the level of juvenile hormone (JH3) entering the brain. Manipulation of the Jhe level reprograms the brain transcriptome between ant castes. Although ant Jhe is retained and functions intracellularly within the BBB, we show that Drosophila Jhe is naturally extracellular. Heterologous expression of ant Jhe into the Drosophila BBB alters behavior in fly to mimic what is seen in ants. Most strikingly, manipulation of Jhe levels in ants reprograms complex behavior between worker castes. Our study thus uncovers a remarkable, potentially conserved role of the BBB serving as a molecular gatekeeper for a neurohormonal pathway that regulates social behavior.

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.

Chromosome-level organization of the regulatory genome in the Drosophila nervous system.

Previous studies have identified topologically associating domains (TADs) as basic units of genome organization. We present evidence of a previously unreported level of genome folding, where distant TAD pairs, megabases apart, interact to form meta-domains. Within meta-domains, gene promoters and structural intergenic elements present in distant TADs are specifically paired. The associated genes encode neuronal determinants, including those engaged in axonal guidance and adhesion. These long-range associations occur in a large fraction of neurons but support transcription in only a subset of neurons. Meta-domains are formed by diverse transcription factors that are able to pair over long and flexible distances. We present evidence that two such factors, GAF and CTCF, play direct roles in this process. The relative simplicity of higher-order meta-domain interactions in Drosophila, compared with those previously described in mammals, allowed the demonstration that genomes can fold into highly specialized cell-type-specific scaffolds that enable megabase-scale regulatory associations.

Linking chromatin acylation mark-defined proteome and genome in living cells.

A generalizable strategy with programmable site specificity for in situ profiling of histone modifications on unperturbed chromatin remains highly desirable but challenging. We herein developed a single-site-resolved multi-omics (SiTomics) strategy for systematic mapping of dynamic modifications and subsequent profiling of chromatinized proteome and genome defined by specific chromatin acylations in living cells. By leveraging the genetic code expansion strategy, our SiTomics toolkit revealed distinct crotonylation (e.g., H3K56cr) and ß-hydroxybutyrylation (e.g., H3K56bhb) upon short chain fatty acids stimulation and established linkages for chromatin acylation mark-defined proteome, genome, and functions. This led to the identification of GLYR1 as a distinct interacting protein in modulating H3K56cr's gene body localization as well as the discovery of an elevated super-enhancer repertoire underlying bhb-mediated chromatin modulations. SiTomics offers a platform technology for elucidating the "metabolites-modification-regulation" axis, which is widely applicable for multi-omics profiling and functional dissection of modifications beyond acylations and proteins beyond histones.

Whole-genome sequencing reveals a complex African population demographic history and signatures of local adaptation.

We conduct high coverage (>30×) whole-genome sequencing of 180 individuals from 12 indigenous African populations. We identify millions of unreported variants, many predicted to be functionally important. We observe that the ancestors of southern African San and central African rainforest hunter-gatherers (RHG) diverged from other populations >200 kya and maintained a large effective population size. We observe evidence for ancient population structure in Africa and for multiple introgression events from "ghost" populations with highly diverged genetic lineages. Although currently geographically isolated, we observe evidence for gene flow between eastern and southern Khoesan-speaking hunter-gatherer populations lasting until ∼12 kya. We identify signatures of local adaptation for traits related to skin color, immune response, height, and metabolic processes. We identify a positively selected variant in the lightly pigmented San that influences pigmentation in vitro by regulating the enhancer activity and gene expression of PDPK1.

A transcription factor atlas of directed differentiation.

Transcription factors (TFs) regulate gene programs, thereby controlling diverse cellular processes and cell states. To comprehensively understand TFs and the programs they control, we created a barcoded library of all annotated human TF splice isoforms (>3,500) and applied it to build a TF Atlas charting expression profiles of human embryonic stem cells (hESCs) overexpressing each TF at single-cell resolution. We mapped TF-induced expression profiles to reference cell types and validated candidate TFs for generation of diverse cell types, spanning all three germ layers and trophoblasts. Targeted screens with subsets of the library allowed us to create a tailored cellular disease model and integrate mRNA expression and chromatin accessibility data to identify downstream regulators. Finally, we characterized the effects of combinatorial TF overexpression by developing and validating a strategy for predicting combinations of TFs that produce target expression profiles matching reference cell types to accelerate cellular engineering efforts.

Repression and 3D-restructuring resolves regulatory conflicts in evolutionarily rearranged genomes.

Regulatory landscapes drive complex developmental gene expression, but it remains unclear how their integrity is maintained when incorporating novel genes and functions during evolution. Here, we investigated how a placental mammal-specific gene, Zfp42, emerged in an ancient vertebrate topologically associated domain (TAD) without adopting or disrupting the conserved expression of its gene, Fat1. In ESCs, physical TAD partitioning separates Zfp42 and Fat1 with distinct local enhancers that drive their independent expression. This separation is driven by chromatin activity and not CTCF/cohesin. In contrast, in embryonic limbs, inactive Zfp42 shares Fat1's intact TAD without responding to active Fat1 enhancers. However, neither Fat1 enhancer-incompatibility nor nuclear envelope-attachment account for Zfp42's unresponsiveness. Rather, Zfp42's promoter is rendered inert to enhancers by context-dependent DNA methylation. Thus, diverse mechanisms enabled the integration of independent Zfp42 regulation in the Fat1 locus. Critically, such regulatory complexity appears common in evolution as, genome wide, most TADs contain multiple independently expressed genes.

Transcription Factor Dynamics: One Molecule at a Time.

Cells must tightly regulate their gene expression programs and yet rapidly respond to acute biochemical and biophysical cues within their environment. This information is transmitted to the nucleus through various signaling cascades, culminating in the activation or repression of target genes. Transcription factors (TFs) are key mediators of these signals, binding to specific regulatory elements within chromatin. While live-cell imaging has conclusively proven that TF-chromatin interactions are highly dynamic, how such transient interactions can have long-term impacts on developmental trajectories and disease progression is still largely unclear. In this review, we summarize our current understanding of the dynamic nature of TF functions, starting with a historical overview of early live-cell experiments. We highlight key factors that govern TF dynamics and how TF dynamics, in turn, affect downstream transcriptional bursting. Finally, we conclude with open challenges and emerging technologies that will further our understanding of transcriptional regulation.

Kr-h1 maintains distinct caste-specific neurotranscriptomes in response to socially regulated hormones.

Behavioral plasticity is key to animal survival. Harpegnathos saltator ants can switch between worker and queen-like status (gamergate) depending on the outcome of social conflicts, providing an opportunity to study how distinct behavioral states are achieved in adult brains. Using social and molecular manipulations in live ants and ant neuronal cultures, we show that ecdysone and juvenile hormone drive molecular and functional differences in the brains of workers and gamergates and direct the transcriptional repressor Kr-h1 to different target genes. Depletion of Kr-h1 in the brain caused de-repression of "socially inappropriate" genes: gamergate genes were upregulated in workers, whereas worker genes were upregulated in gamergates. At the phenotypic level, loss of Kr-h1 resulted in the emergence of worker-specific behaviors in gamergates and gamergate-specific traits in workers. We conclude that Kr-h1 is a transcription factor that maintains distinct brain states established in response to socially regulated hormones.

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.

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.

Dynamic landscape of immune cell-specific gene regulation in immune-mediated diseases.

Genetic studies have revealed many variant loci that are associated with immune-mediated diseases. To elucidate the disease pathogenesis, it is essential to understand the function of these variants, especially under disease-associated conditions. Here, we performed a large-scale immune cell gene-expression analysis, together with whole-genome sequence analysis. Our dataset consists of 28 distinct immune cell subsets from 337 patients diagnosed with 10 categories of immune-mediated diseases and 79 healthy volunteers. Our dataset captured distinctive gene-expression profiles across immune cell types and diseases. Expression quantitative trait loci (eQTL) analysis revealed dynamic variations of eQTL effects in the context of immunological conditions, as well as cell types. These cell-type-specific and context-dependent eQTLs showed significant enrichment in immune disease-associated genetic variants, and they implicated the disease-relevant cell types, genes, and environment. This atlas deepens our understanding of the immunogenetic functions of disease-associated variants under in vivo disease conditions.

Senescence: An Identity Crisis Originating from Deep Within the Nucleus.

Cellular senescence is implicated in a wide range of physiological and pathological conditions throughout an organism's entire lifetime. In particular, it has become evident that senescence plays a causative role in aging and age-associated disorders. This is not due simply to the loss of function of senescent cells. Instead, the substantial alterations of the cellular activities of senescent cells, especially the array of secretory factors, impact the surrounding tissues or even entire organisms. Such non-cell-autonomous functionality is largely coordinated by tissue-specific genes, constituting a cell fate-determining state. Senescence can be viewed as a gain-of-function phenotype or a process of cell identity shift. Cellular functionality or lineage-specific gene expression is tightly linked to the cell type-specific epigenetic landscape, reinforcing the heterogeneity of senescence across cell types. Here, we aim to define the senescence cellular functionality and epigenetic features that may contribute to the gain-of-function phenotype.

Evaluating Enhancer Function and Transcription.

Cell-type- and condition-specific profiles of gene expression require coordination between protein-coding gene promoters and cis-regulatory sequences called enhancers. Enhancers can stimulate gene activity at great genomic distances from their targets, raising questions about how enhancers communicate with specific gene promoters and what molecular mechanisms underlie enhancer function. Characterization of enhancer loci has identified the molecular features of active enhancers that accompany the binding of transcription factors and local opening of chromatin. These characteristics include coactivator recruitment, histone modifications, and noncoding RNA transcription. However, it remains unclear which of these features functionally contribute to enhancer activity. Here, we discuss what is known about how enhancers regulate their target genes and how enhancers and promoters communicate. Further, we describe recent data demonstrating many similarities between enhancers and the gene promoters they control, and we highlight unanswered questions in the field, such as the potential roles of transcription at enhancers.

Transcription in Living Cells: Molecular Mechanisms of Bursting.

Transcription in several organisms from certain bacteria to humans has been observed to be stochastic in nature: toggling between active and inactive states. Periods of active nascent RNA synthesis known as bursts represent individual gene activation events in which multiple polymerases are initiated. Therefore, bursting is the single locus illustration of both gene activation and repression. Although transcriptional bursting was originally observed decades ago, only recently have technological advances enabled the field to begin elucidating gene regulation at the single-locus level. In this review, we focus on how biochemical, genomic, and single-cell data describe the regulatory steps of transcriptional bursts.

Dynamic Competition of Polycomb and Trithorax in Transcriptional Programming.

Predicting regulatory potential from primary DNA sequences or transcription factor binding patterns is not possible. However, the annotation of the genome by chromatin proteins, histone modifications, and differential compaction is largely sufficient to reveal the locations of genes and their differential activity states. The Polycomb Group (PcG) and Trithorax Group (TrxG) proteins are the central players in this cell type-specific chromatin organization. PcG function was originally viewed as being solely repressive and irreversible, as observed at the homeotic loci in flies and mammals. However, it is now clear that modular and reversible PcG function is essential at most developmental genes. Focusing mainly on recent advances, we review evidence for how PcG and TrxG patterns change dynamically during cell type transitions. The ability to implement cell type-specific transcriptional programming with exquisite fidelity is essential for normal development.

Chromatin Potential Identified by Shared Single-Cell Profiling of RNA and Chromatin.

Cell differentiation and function are regulated across multiple layers of gene regulation, including modulation of gene expression by changes in chromatin accessibility. However, differentiation is an asynchronous process precluding a temporal understanding of regulatory events leading to cell fate commitment. Here we developed simultaneous high-throughput ATAC and RNA expression with sequencing (SHARE-seq), a highly scalable approach for measurement of chromatin accessibility and gene expression in the same single cell, applicable to different tissues. Using 34,774 joint profiles from mouse skin, we develop a computational strategy to identify cis-regulatory interactions and define domains of regulatory chromatin (DORCs) that significantly overlap with super-enhancers. During lineage commitment, chromatin accessibility at DORCs precedes gene expression, suggesting that changes in chromatin accessibility may prime cells for lineage commitment. We computationally infer chromatin potential as a quantitative measure of chromatin lineage-priming and use it to predict cell fate outcomes. SHARE-seq is an extensible platform to study regulatory circuitry across diverse cells in tissues.