Higher-Order Inter-chromosomal Hubs Shape 3D Genome Organization in the Nucleus.

Eukaryotic genomes are packaged into a 3-dimensional structure in the nucleus. Current methods for studying genome-wide structure are based on proximity ligation. However, this approach can fail to detect known structures, such as interactions with nuclear bodies, because these DNA regions can be too far apart to directly ligate. Accordingly, our overall understanding of genome organization remains incomplete. Here, we develop split-pool recognition of interactions by tag extension (SPRITE), a method that enables genome-wide detection of higher-order interactions within the nucleus. Using SPRITE, we recapitulate known structures identified by proximity ligation and identify additional interactions occurring across larger distances, including two hubs of inter-chromosomal interactions that are arranged around the nucleolus and nuclear speckles. We show that a substantial fraction of the genome exhibits preferential organization relative to these nuclear bodies. Our results generate a global model whereby nuclear bodies act as inter-chromosomal hubs that shape the overall packaging of DNA in the nucleus.

HIV-1-induced cytokines deplete homeostatic innate lymphoid cells and expand TCF7-dependent memory NK cells.

Human immunodeficiency virus 1 (HIV-1) infection is associated with heightened inflammation and excess risk of cardiovascular disease, cancer and other complications. These pathologies persist despite antiretroviral therapy. In two independent cohorts, we found that innate lymphoid cells (ILCs) were depleted in the blood and gut of people with HIV-1, even with effective antiretroviral therapy. ILC depletion was associated with neutrophil infiltration of the gut lamina propria, type 1 interferon activation, increased microbial translocation and natural killer (NK) cell skewing towards an inflammatory state, with chromatin structure and phenotype typical of WNT transcription factor TCF7-dependent memory T cells. Cytokines that are elevated during acute HIV-1 infection reproduced the ILC and NK cell abnormalities ex vivo. These results show that inflammatory cytokines associated with HIV-1 infection irreversibly disrupt ILCs. This results in loss of gut epithelial integrity, microbial translocation and memory NK cells with heightened inflammatory potential, and explains the chronic inflammation in people with HIV-1.

Systematic Functional Dissection of Common Genetic Variation Affecting Red Blood Cell Traits.

Genome-wide association studies (GWAS) have successfully identified thousands of associations between common genetic variants and human disease phenotypes, but the majority of these variants are non-coding, often requiring genetic fine-mapping, epigenomic profiling, and individual reporter assays to delineate potential causal variants. We employ a massively parallel reporter assay (MPRA) to simultaneously screen 2,756 variants in strong linkage disequilibrium with 75 sentinel variants associated with red blood cell traits. We show that this assay identifies elements with endogenous erythroid regulatory activity. Across 23 sentinel variants, we conservatively identified 32 MPRA functional variants (MFVs). We used targeted genome editing to demonstrate endogenous enhancer activity across 3 MFVs that predominantly affect the transcription of SMIM1, RBM38, and CD164. Functional follow-up of RBM38 delineates a key role for this gene in the alternative splicing program occurring during terminal erythropoiesis. Finally, we provide evidence for how common GWAS-nominated variants can disrupt cell-type-specific transcriptional regulatory pathways.

Ebola Virus Glycoprotein with Increased Infectivity Dominated the 2013-2016 Epidemic.

The magnitude of the 2013-2016 Ebola virus disease (EVD) epidemic enabled an unprecedented number of viral mutations to occur over successive human-to-human transmission events, increasing the probability that adaptation to the human host occurred during the outbreak. We investigated one nonsynonymous mutation, Ebola virus (EBOV) glycoprotein (GP) mutant A82V, for its effect on viral infectivity. This mutation, located at the NPC1-binding site on EBOV GP, occurred early in the 2013-2016 outbreak and rose to high frequency. We found that GP-A82V had heightened ability to infect primate cells, including human dendritic cells. The increased infectivity was restricted to cells that have primate-specific NPC1 sequences at the EBOV interface, suggesting that this mutation was indeed an adaptation to the human host. GP-A82V was associated with increased mortality, consistent with the hypothesis that the heightened intrinsic infectivity of GP-A82V contributed to disease severity during the EVD epidemic.

RNA-RNA interactions enable specific targeting of noncoding RNAs to nascent Pre-mRNAs and chromatin sites.

Intermolecular RNA-RNA interactions are used by many noncoding RNAs (ncRNAs) to achieve their diverse functions. To identify these contacts, we developed a method based on RNA antisense purification to systematically map RNA-RNA interactions (RAP-RNA) and applied it to investigate two ncRNAs implicated in RNA processing: U1 small nuclear RNA, a component of the spliceosome, and Malat1, a large ncRNA that localizes to nuclear speckles. U1 and Malat1 interact with nascent transcripts through distinct targeting mechanisms. Using differential crosslinking, we confirmed that U1 directly hybridizes to 5' splice sites and 5' splice site motifs throughout introns and found that Malat1 interacts with pre-mRNAs indirectly through protein intermediates. Interactions with nascent pre-mRNAs cause U1 and Malat1 to localize proximally to chromatin at active genes, demonstrating that ncRNAs can use RNA-RNA interactions to target specific pre-mRNAs and genomic sites. RAP-RNA is sensitive to lower abundance RNAs as well, making it generally applicable for investigating ncRNAs.

High-Resolution Mapping of Multiway Enhancer-Promoter Interactions Regulating Pathogen Detection.

Eukaryotic gene expression regulation involves thousands of distal regulatory elements. Understanding the quantitative contribution of individual enhancers to gene expression is critical for assessing the role of disease-associated genetic risk variants. Yet, we lack the ability to accurately link genes with their distal regulatory elements. To address this, we used 3D enhancer-promoter (E-P) associations identified using split-pool recognition of interactions by tag extension (SPRITE) to build a predictive model of gene expression. Our model dramatically outperforms models using genomic proximity and can be used to determine the quantitative impact of enhancer loss on gene expression in different genetic backgrounds. We show that genes that form stable E-P hubs have less cell-to-cell variability in gene expression. Finally, we identified transcription factors that regulate stimulation-dependent E-P interactions. Together, our results provide a framework for understanding quantitative contributions of E-P interactions and associated genetic variants to gene expression.

Local regulation of gene expression by lncRNA promoters, transcription and splicing.

Mammalian genomes are pervasively transcribed to produce thousands of long non-coding RNAs (lncRNAs). A few of these lncRNAs have been shown to recruit regulatory complexes through RNA-protein interactions to influence the expression of nearby genes, and it has been suggested that many other lncRNAs can also act as local regulators. Such local functions could explain the observation that lncRNA expression is often correlated with the expression of nearby genes. However, these correlations have been challenging to dissect and could alternatively result from processes that are not mediated by the lncRNA transcripts themselves. For example, some gene promoters have been proposed to have dual functions as enhancers, and the process of transcription itself may contribute to gene regulation by recruiting activating factors or remodelling nucleosomes. Here we use genetic manipulation in mouse cell lines to dissect 12 genomic loci that produce lncRNAs and find that 5 of these loci influence the expression of a neighbouring gene in cis. Notably, none of these effects requires the specific lncRNA transcripts themselves and instead involves general processes associated with their production, including enhancer-like activity of gene promoters, the process of transcription, and the splicing of the transcript. Furthermore, such effects are not limited to lncRNA loci: we find that four out of six protein-coding loci also influence the expression of a neighbour. These results demonstrate that cross-talk among neighbouring genes is a prevalent phenomenon that can involve multiple mechanisms and cis-regulatory signals, including a role for RNA splice sites. These mechanisms may explain the function and evolution of some genomic loci that produce lncRNAs and broadly contribute to the regulation of both coding and non-coding genes.

The Xist lncRNA interacts directly with SHARP to silence transcription through HDAC3.

Many long non-coding RNAs (lncRNAs) affect gene expression, but the mechanisms by which they act are still largely unknown. One of the best-studied lncRNAs is Xist, which is required for transcriptional silencing of one X chromosome during development in female mammals. Despite extensive efforts to define the mechanism of Xist-mediated transcriptional silencing, we still do not know any proteins required for this role. The main challenge is that there are currently no methods to comprehensively define the proteins that directly interact with a lncRNA in the cell. Here we develop a method to purify a lncRNA from cells and identify proteins interacting with it directly using quantitative mass spectrometry. We identify ten proteins that specifically associate with Xist, three of these proteins--SHARP, SAF-A and LBR--are required for Xist-mediated transcriptional silencing. We show that SHARP, which interacts with the SMRT co-repressor that activates HDAC3, is not only essential for silencing, but is also required for the exclusion of RNA polymerase II (Pol II) from the inactive X. Both SMRT and HDAC3 are also required for silencing and Pol II exclusion. In addition to silencing transcription, SHARP and HDAC3 are required for Xist-mediated recruitment of the polycomb repressive complex 2 (PRC2) across the X chromosome. Our results suggest that Xist silences transcription by directly interacting with SHARP, recruiting SMRT, activating HDAC3, and deacetylating histones to exclude Pol II across the X chromosome.

Insight into GATA1 transcriptional activity through interrogation of cis elements disrupted in human erythroid disorders.

Whole-exome sequencing has been incredibly successful in identifying causal genetic variants and has revealed a number of novel genes associated with blood and other diseases. One limitation of this approach is that it overlooks mutations in noncoding regulatory elements. Furthermore, the mechanisms by which mutations in transcriptionalcis-regulatory elements result in disease remain poorly understood. Here we used CRISPR/Cas9 genome editing to interrogate three such elements harboring mutations in human erythroid disorders, which in all cases are predicted to disrupt a canonical binding motif for the hematopoietic transcription factor GATA1. Deletions of as few as two to four nucleotides resulted in a substantial decrease (>80%) in target gene expression. Isolated deletions of the canonical GATA1 binding motif completely abrogated binding of the cofactor TAL1, which binds to a separate motif. Having verified the functionality of these three GATA1 motifs, we demonstrate strong evolutionary conservation of GATA1 motifs in regulatory elements proximal to other genes implicated in erythroid disorders, and show that targeted disruption of such elements results in altered gene expression. By modeling transcription factor binding patterns, we show that multiple transcription factors are associated with erythroid gene expression, and have created predictive maps modeling putative disruptions of their binding sites at key regulatory elements. Our study provides insight into GATA1 transcriptional activity and may prove a useful resource for investigating the pathogenicity of noncoding variants in human erythroid disorders.

Multiple cell and population-level interactions with mouse embryonic stem cell heterogeneity.

Much of development and disease concerns the generation of gene expression differences between related cells sharing similar niches. However, most analyses of gene expression only assess population and time-averaged levels of steady-state transcription. The mechanisms driving differentiation are buried within snapshots of the average cell, lacking dynamic information and the diverse regulatory history experienced by individual cells. Here, we use a quantitative imaging platform with large time series data sets to determine the regulation of developmental gene expression by cell cycle, lineage, motility and environment. We apply this technology to the regulation of the pluripotency gene Nanog in mouse embryonic stem cells. Our data reveal the diversity of cell and population-level interactions with Nanog dynamics and heterogeneity, and how this regulation responds to triggers of pluripotency. Cell cycles are highly heterogeneous and cycle time increases with Nanog reporter expression, with longer, more variable cycle times as cells approach ground-state pluripotency. Nanog reporter expression is highly stable over multiple cell generations, with fluctuations within cycles confined by an attractor state. Modelling reveals an environmental component to expression stability, in addition to any cell-autonomous behaviour, and we identify interactions of cell density with both cycle behaviour and Nanog. Rex1 expression dynamics showed shared and distinct regulatory effects. Overall, our observations of multiple partially overlapping dynamic heterogeneities imply complex cell and environmental regulation of pluripotent cell behaviour, and suggest simple deterministic views of stem cell states are inappropriate.

Sin3a is essential for the genome integrity and viability of pluripotent cells.

The Sin3a/HDAC co-repressor complex is a critical regulator of transcription networks that govern cell cycle control and apoptosis throughout development. Previous studies have identified Sin3a as essential for embryonic development around the time of implantation, during which the epiblast cell cycle is uniquely structured to achieve very rapid divisions with little tolerance of DNA damage. This study investigates the specific requirement for Sin3a in the early mouse embryo and shows that embryos lacking Sin3a suffer unresolved DNA damage and acute p53-independent apoptosis specifically in the E3.5-4.5 epiblast. Surprisingly, Myc and E2F targets in Sin3a-null ICMs are downregulated, suggesting a central but non-canonical role for Sin3a in regulating the pluripotent embryonic cell cycle. ES cells deleted for Sin3a mount a DNA damage response indicative of unresolved double-strand breaks, profoundly arrest at G2, and undergo apoptosis. These results indicate that Sin3a protects the genomic integrity of pluripotent embryonic cells and governs their unusual cell cycle.

Clustered DNA motifs mark X chromosomes for repression by a dosage compensation complex.

Gene expression in metazoans is regulated not only at the level of individual genes but also in a coordinated manner across large chromosomal domains (for example centromeres, telomeres and imprinted gene clusters) and along entire chromosomes (for example X-chromosome dosage compensation). The primary DNA sequence usually specifies the regulation of individual genes, but the nature of cis-acting information that controls genes over large regions has been elusive: higher-order DNA structure, specific histone modifications, subnuclear compartmentalization and primary DNA sequence are possibilities. One paradigm of chromosome-wide gene regulation is Caenorhabditis elegans dosage compensation in which a large dosage compensation complex (DCC) is targeted to both X chromosomes of hermaphrodites to repress transcript levels by half. This essential process equalizes X-linked gene expression between the sexes (XO males and XX hermaphrodites). Here we report the discovery and dissection of cis-acting sites that mark nematode X chromosomes as targets for gene repression by the DCC. These rex (recruitment element on X) sites are widely dispersed along X and reside in promoters, exons and intergenic regions. rex sites share at least two distinct motifs that act in combination to recruit the DCC. Mutating these motifs severely reduces or abolishes DCC binding in vivo, demonstrating the importance of primary DNA sequence in chromosome-wide regulation. Unexpectedly, the motifs are not enriched on X, but altering motif numbers within rex sites demonstrates that motif co-occurrence in unusually high densities is essential for optimal DCC recruitment. Thus, X-specific repression is established through sequences not specific to X. The distribution of common motifs provides the foundation for repression along an entire chromosome.

scRNA-seq of human vitiligo reveals complex networks of subclinical immune activation and a role for CCR5 in Treg function.

Vitiligo is an autoimmune skin disease characterized by the targeted destruction of melanocytes by T cells. Cytokine signaling between keratinocytes and T cells results in CD8+ T cell infiltration of vitiligo lesions, but the full scope of signals required to coordinate autoimmune responses is not completely understood. We performed single-cell RNA sequencing on affected and unaffected skin from patients with vitiligo, as well as healthy controls, to define the role of each cell type in coordinating autoimmunity during disease progression. We confirmed that type 1 cytokine signaling occupied a central role in disease, but we also found that this pathway was used by regulatory T cells (Tregs) to restrain disease progression in nonlesional skin. We determined that CCL5-CCR5 signaling served as a chemokine circuit between effector CD8+ T cells and Tregs, and mechanistic studies in a mouse model of vitiligo revealed that CCR5 expression on Tregs was required to suppress disease in vivo but not in vitro. CCR5 was not required for Treg recruitment to skin but appeared to facilitate Treg function by properly positioning these cells within the skin. Our data provide critical insights into the pathogenesis of vitiligo and uncover potential opportunities for therapeutic interventions.

Keeping things quiet: roles of NuRD and Sin3 co-repressor complexes during mammalian development.

Gene inactivation studies of mammalian histone and DNA-modifying proteins have demonstrated a role for many such proteins in embryonic development. Post-implantation embryonic lethality implies a role for epigenetic factors in differentiation and in development of specific lineages or tissues. However a handful of chromatin-modifying enzymes have been found to be required in pre- or peri-implantation embryos. This is significant as implantation is the time when inner cell mass cells of the blastocyst exit pluripotency and begin to commit to form the various lineages that will eventually form the adult animal. These observations indicate a critical role for chromatin-modifying proteins in the earliest lineage decisions of mammalian development, and/or in the formation of the first embryonic cell types. Recent work has shown that the two major class I histone deacetylase-containing co-repressor complexes, the NuRD and Sin3 complexes, are both required at peri-implantation stages of mouse development, demonstrating the importance of histone deacetylation in cell fate decisions. Over the past 10 years both genetic and biochemical studies have revealed surprisingly divergent roles for these two co-repressors in mammalian cells. In this review we will summarise the evidence that the two major class I histone deacetylase complexes in mammalian cells, the NuRD and Sin3 complexes, play important roles in distinct aspects of embryonic development.

Approaches for Understanding the Mechanisms of Long Noncoding RNA Regulation of Gene Expression.

Mammalian genomes encode tens of thousands of long noncoding RNAs (lncRNAs) that have been implicated in a diverse array of biological processes and human diseases. In recent years, the development of new tools for studying lncRNAs has enabled important progress in defining the mechanisms by which Xist and other lncRNAs function. This collective work provides a framework for how to define the mechanisms by which lncRNAs act. This includes defining lncRNA function, identifying and characterizing lncRNA-protein interactions, and lncRNA localization in the cell. In this review, we discuss various experimental approaches for deciphering lncRNA mechanisms and discuss issues and limitations in interpreting these results. We explore what these data can reveal about lncRNA function and mechanism as well as emerging insights into lncRNA biology that have been derived from these studies.

Comparative Analysis of Immune Cells Reveals a Conserved Regulatory Lexicon.

Most well-characterized enhancers are deeply conserved. In contrast, genome-wide comparative studies of steady-state systems showed that only a small fraction of active enhancers are conserved. To better understand conservation of enhancer activity, we used a comparative genomics approach that integrates temporal expression and epigenetic profiles in an innate immune system. We found that gene expression programs diverge among mildly induced genes, while being highly conserved for strongly induced genes. The fraction of conserved enhancers varies greatly across gene expression programs, with induced genes and early-response genes, in particular, being regulated by a higher fraction of conserved enhancers. Clustering of conserved accessible DNA sequences within enhancers resulted in over 60 sequence motifs including motifs for known factors, as well as many with unknown function. We further show that the number of instances of these motifs is a strong predictor of the responsiveness of a gene to pathogen detection.

Xist recruits the X chromosome to the nuclear lamina to enable chromosome-wide silencing.

The Xist long noncoding RNA orchestrates X chromosome inactivation, a process that entails chromosome-wide silencing and remodeling of the three-dimensional (3D) structure of the X chromosome. Yet, it remains unclear whether these changes in nuclear structure are mediated by Xist and whether they are required for silencing. Here, we show that Xist directly interacts with the Lamin B receptor, an integral component of the nuclear lamina, and that this interaction is required for Xist-mediated silencing by recruiting the inactive X to the nuclear lamina and by doing so enables Xist to spread to actively transcribed genes across the X. Our results demonstrate that lamina recruitment changes the 3D structure of DNA, enabling Xist and its silencing proteins to spread across the X to silence transcription.

The Xist lncRNA exploits three-dimensional genome architecture to spread across the X chromosome.

Many large noncoding RNAs (lncRNAs) regulate chromatin, but the mechanisms by which they localize to genomic targets remain unexplored. We investigated the localization mechanisms of the Xist lncRNA during X-chromosome inactivation (XCI), a paradigm of lncRNA-mediated chromatin regulation. During the maintenance of XCI, Xist binds broadly across the X chromosome. During initiation of XCI, Xist initially transfers to distal regions across the X chromosome that are not defined by specific sequences. Instead, Xist identifies these regions by exploiting the three-dimensional conformation of the X chromosome. Xist requires its silencing domain to spread across actively transcribed regions and thereby access the entire chromosome. These findings suggest a model in which Xist coats the X chromosome by searching in three dimensions, modifying chromosome structure, and spreading to newly accessible locations.

Phylogenetic relationships in Gleditsia (Leguminosae) based on ITS sequences.

We used nucleotide sequences from the internal transcribed spacers and 5.8S gene of nuclear ribosomal DNA to test competing phylogenetic and biogeographic hypotheses in Gleditsia. Eleven of 13 Gleditsia species were sampled, along with two species of its sister genus, Gymnocladus. Analyses of ITS data and of a combined data set that included sequences of ITS and two chloroplast genes supported several conclusions that were interpreted in light of fossil data and current legume phylogeny. Gleditsia and Gymnocladus appear to have originated in eastern Asia during the Eocene. Eastern North American species of both genera most likely evolved from ancestors that migrated across the Bering land bridge, but the eastern Asian/eastern North American disjunction appears to be much older in Gymnocladus than in Gleditsia. Gleditsia amorphoides, from temperate South America, is sister to the rest of the genus, suggesting early long-distance dispersal from Asia. The remainder of Gleditsia is divided into three unresolved clades, possibly indicating a split early in the evolution of the genus. Two of those clades contain only Asian species, and one contains Asian and North American species. The North American species, Gleditsia triacanthos and Gleditsia aquatica, are polymorphic and paraphyletic with respect to their ITS and cpDNA sequences, which suggests recent diversification.

Recruitment and spreading of the C. elegans dosage compensation complex along X chromosomes.

To achieve X-chromosome dosage compensation, organisms must distinguish X chromosomes from autosomes. We identified multiple, cis-acting regions that recruit the Caenorhabditis elegans dosage compensation complex (DCC) through a search for regions of X that bind the complex when detached from X. The DCC normally assembles along the entire X chromosome, but not all detached regions recruit the complex, despite having genes known to be dosage compensated on the native X. Thus, the DCC binds first to recruitment sites, then spreads to neighboring X regions to accomplish chromosome-wide gene repression. From a large chromosomal domain, we defined a 793-base pair fragment that functions in vivo as an X-recognition element to recruit the DCC.