Search

We sought to better understand the immune response during the immediate post-diagnosis phase of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by identifying molecular associations with longitudinal disease outcomes. Multi-omic analyses identified differences in immune cell composition, cytokine levels, and cell subset-specific transcriptomic and epigenomic signatures between individuals on a more serious disease trajectory (Progressors) as compared to those on a milder course (Non-progressors). Higher levels of multiple cytokines were observed in Progressors, with IL-6 showing the largest difference. Blood monocyte cell subsets were also skewed, showing a comparative decrease in non-classical CD14 - CD16 + and intermediate CD14 + CD16 + monocytes. In lymphocytes, the CD8 + T effector memory cells displayed a gene expression signature consistent with stronger T cell activation in Progressors. These early stage observations could serve as the basis for the development of prognostic biomarkers of disease risk and interventional strategies to improve the management of severe COVID-19. BACKGROUND: Much of the literature on immune response post-SARS-CoV-2 infection has been in the acute and post-acute phases of infection. TRANSLATIONAL SIGNIFICANCE: We found differences at early time points of infection in approximately 160 participants. We compared multi-omic signatures in immune cells between individuals progressing to needing more significant medical intervention and non-progressors. We observed widespread evidence of a state of increased inflammation associated with progression, supported by a range of epigenomic, transcriptomic, and proteomic signatures. The signatures we identified support other findings at later time points and serve as the basis for prognostic biomarker development or to inform interventional strategies., (Copyright © 2023 Verily Life Sciences LLC. Published by Elsevier Inc. All rights reserved.)

Physiologic laminar shear stress (LSS) induces an endothelial gene expression profile that is vasculo-protective. In this report, we delineate how LSS mediates changes in the epigenetic landscape to promote this beneficial response. We show that under LSS, KLF4 interacts with the SWI/SNF nucleosome remodeling complex to increase accessibility at enhancer sites that promote the expression of homeostatic endothelial genes. By combining molecular and computational approaches we discover enhancers that loop to promoters of KLF4- and LSS-responsive genes that stabilize endothelial cells and suppress inflammation, such as BMPR2, SMAD5, and DUSP5. By linking enhancers to genes that they regulate under physiologic LSS, our work establishes a foundation for interpreting how non-coding DNA variants in these regions might disrupt protective gene expression to influence vascular disease., (© 2022. The Author(s).)

Macrophages use an array of innate immune sensors to detect intracellular pathogens and to tailor effective antimicrobial responses. In addition, extrinsic activation with the cytokine IFN-γ is often required as well to tip the scales of the host-pathogen balance toward pathogen restriction. However, little is known about how host-pathogen sensing impacts the antimicrobial IFN-γ-activated state. It was observed that in the absence of IRF3, a key downstream component of pathogen sensing pathways, IFN-γ-primed macrophages more efficiently restricted the intracellular bacterium Legionella pneumophila and the intracellular protozoan parasite Trypanosoma cruzi. This effect did not require IFNAR, the receptor for Type I IFNs known to be induced by IRF3, nor the sensing adaptors MyD88/TRIF, MAVS, or STING. This effect also did not involve differential activation of STAT1, the major signaling protein downstream of both Type 1 and Type 2 IFN receptors. IRF3-deficient macrophages displayed a significantly altered IFN-γ-induced gene expression program, with up-regulation of microbial restriction factors such as Nos2. Finally, we found that IFN-γ-primed but not unprimed macrophages largely excluded the activated form of IRF3 from the nucleus following bacterial infection. These data are consistent with a relationship of mutual inhibition between IRF3 and IFN-γ-activated programs, possibly as a component of a partially reversible mechanism for modulating the activity of potent innate immune effectors (such as Nos2) in the context of intracellular infection., (©2021 Society for Leukocyte Biology.)

Dysfunction in T cells limits the efficacy of cancer immunotherapy. We profiled the epigenome, transcriptome, and enhancer connectome of exhaustion-prone GD2-targeting HA-28z chimeric antigen receptor (CAR) T cells and control CD19-targeting CAR T cells, which present less exhaustion-inducing tonic signaling, at multiple points during their ex vivo expansion. We found widespread, dynamic changes in chromatin accessibility and three-dimensional (3D) chromosome conformation preceding changes in gene expression, notably at loci proximal to exhaustion-associated genes such as PDCD1 , CTLA4 , and HAVCR2 , and increased DNA motif access for AP-1 family transcription factors, which are known to promote exhaustion. Although T cell exhaustion has been studied in detail in mice, we find that the regulatory networks of T cell exhaustion differ between species and involve distinct loci of accessible chromatin and cis-regulated target genes in human CAR T cell exhaustion. Deletion of exhaustion-specific candidate enhancers of PDCD1 suppress the expression of PD-1 in an in vitro model of T cell dysfunction and in HA-28z CAR T cells, suggesting enhancer editing as a path forward in improving cancer immunotherapy., Competing Interests: Competing interest statement: H.Y.C. is a cofounder of Accent Therapeutics and Boundless Bio and is an advisor to 10× Genomics, Arsenal Bioscience, and Spring Discovery. C.L.M. is a cofounder of Lyell Immunopharma. R.C.L. is employed by and E.W.W. is a consultant for Lyell Immunopharma. A.T.S. is a cofounder of Immunai., (Copyright © 2021 the Author(s). Published by PNAS.)

SARS-CoV-2 is the cause of a pandemic with growing global mortality. Using comprehensive identification of RNA-binding proteins by mass spectrometry (ChIRP-MS), we identified 309 host proteins that bind the SARS-CoV-2 RNA during active infection. Integration of this data with ChIRP-MS data from three other RNA viruses defined viral specificity of RNA-host protein interactions. Targeted CRISPR screens revealed that the majority of functional RNA-binding proteins protect the host from virus-induced cell death, and comparative CRISPR screens across seven RNA viruses revealed shared and SARS-specific antiviral factors. Finally, by combining the RNA-centric approach and functional CRISPR screens, we demonstrated a physical and functional connection between SARS-CoV-2 and mitochondria, highlighting this organelle as a general platform for antiviral activity. Altogether, these data provide a comprehensive catalog of functional SARS-CoV-2 RNA-host protein interactions, which may inform studies to understand the host-virus interface and nominate host pathways that could be targeted for therapeutic benefit., Competing Interests: Declaration of interests A.T.S. is a scientific co-founder of Immunai and receives research funding from Arsenal Biosciences and 10x Genomics. K.R.P., H.Y.C., and A.T.S. are co-founders of Cartography Biosciences. H.Y.C. is a co-founder of Accent Therapeutics, Boundless Bio, and an advisor for 10x Genomics, Arsenal Biosciences, and Spring Discovery. Yale University (C.B.W.) has a patent pending related to this work entitled: “Compounds and Compositions for Treating, Ameliorating, and/or Preventing SARS-CoV-2 Infection and/or Complications Thereof.” Yale University has committed to rapidly executable non-exclusive royalty-free licenses to intellectual property rights for the purpose of making and distributing products to prevent, diagnose, and treat COVID-19 infection during the pandemic and for a short period thereafter., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

Genome-wide association studies of neurological diseases have identified thousands of variants associated with disease phenotypes. However, most of these variants do not alter coding sequences, making it difficult to assign their function. Here, we present a multi-omic epigenetic atlas of the adult human brain through profiling of single-cell chromatin accessibility landscapes and three-dimensional chromatin interactions of diverse adult brain regions across a cohort of cognitively healthy individuals. We developed a machine-learning classifier to integrate this multi-omic framework and predict dozens of functional SNPs for Alzheimer's and Parkinson's diseases, nominating target genes and cell types for previously orphaned loci from genome-wide association studies. Moreover, we dissected the complex inverted haplotype of the MAPT (encoding tau) Parkinson's disease risk locus, identifying putative ectopic regulatory interactions in neurons that may mediate this disease association. This work expands understanding of inherited variation and provides a roadmap for the epigenomic dissection of causal regulatory variation in disease.

Transposable elements (TEs) comprise nearly half of the human genome and are often transcribed or exhibit cis-regulatory properties with unknown function in specific processes such as heart development. In the case of endogenous retroviruses (ERVs), a TE subclass, experimental interrogation is constrained as many are primate-specific or human-specific. Here, we use primate pluripotent stem-cell-derived cardiomyocytes that mimic fetal cardiomyocytes in vitro to discover hundreds of ERV transcripts from the primate-specific MER41 family, some of which are regulated by the cardiogenic transcription factor TBX5. The most significant of these are located within BANCR, a long non-coding RNA (lncRNA) exclusively expressed in primate fetal cardiomyocytes. Functional studies reveal that BANCR promotes cardiomyocyte migration in vitro and ventricular enlargement in vivo. We conclude that recently evolved TE loci such as BANCR may represent potent de novo developmental regulatory elements that can be interrogated with species-matching pluripotent stem cell models., Competing Interests: Declaration of Interests H.Y.C. is a co-founder of Accent Therapeutics and Boundless Bio and an advisor to 10x Genomics, Arsenal Biosciences, and Spring Discovery. J.C.W. is a co-founder of Khloris Biosciences but has no competing interests, as the work presented here is completely independent. The other authors declare no competing interests., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Genetic variations underlying susceptibility to complex autoimmune and allergic diseases are concentrated within noncoding regulatory elements termed enhancers 1 . The functions of a large majority of disease-associated enhancers are unknown, in part owing to their distance from the genes they regulate, a lack of understanding of the cell types in which they operate, and our inability to recapitulate the biology of immune diseases in vitro. Here, using shared synteny to guide loss-of-function analysis of homologues of human enhancers in mice, we show that the prominent autoimmune and allergic disease risk locus at chromosome 11q13.5 2-7 contains a distal enhancer that is functional in CD4 + regulatory T (T reg ) cells and required for T reg -mediated suppression of colitis. The enhancer recruits the transcription factors STAT5 and NF-κB to mediate signal-driven expression of Lrrc32, which encodes the protein glycoprotein A repetitions predominant (GARP). Whereas disruption of the Lrrc32 gene results in early lethality, mice lacking the enhancer are viable but lack GARP expression in Foxp3 + T reg cells, which are unable to control colitis in a cell-transfer model of the disease. In human T reg cells, the enhancer forms conformational interactions with the promoter of LRRC32 and enhancer risk variants are associated with reduced histone acetylation and GARP expression. Finally, functional fine-mapping of 11q13.5 using CRISPR-activation (CRISPRa) identifies a CRISPRa-responsive element in the vicinity of risk variant rs11236797 capable of driving GARP expression. These findings provide a mechanistic basis for association of the 11q13.5 risk locus with immune-mediated diseases and identify GARP as a potential target in their therapy.

Systemic sclerosis (SSc) is an autoimmune disease that shows one of the highest mortality rates among rheumatic diseases. We perform a large genome-wide association study (GWAS), and meta-analysis with previous GWASs, in 26,679 individuals and identify 27 independent genome-wide associated signals, including 13 new risk loci. The novel associations nearly double the number of genome-wide hits reported for SSc thus far. We define 95% credible sets of less than 5 likely causal variants in 12 loci. Additionally, we identify specific SSc subtype-associated signals. Functional analysis of high-priority variants shows the potential function of SSc signals, with the identification of 43 robust target genes through HiChIP. Our results point towards molecular pathways potentially involved in vasculopathy and fibrosis, two main hallmarks in SSc, and highlight the spectrum of critical cell types for the disease. This work supports a better understanding of the genetic basis of SSc and provides directions for future functional experiments.

Understanding complex tissues requires single-cell deconstruction of gene regulation with precision and scale. Here, we assess the performance of a massively parallel droplet-based method for mapping transposase-accessible chromatin in single cells using sequencing (scATAC-seq). We apply scATAC-seq to obtain chromatin profiles of more than 200,000 single cells in human blood and basal cell carcinoma. In blood, application of scATAC-seq enables marker-free identification of cell type-specific cis- and trans-regulatory elements, mapping of disease-associated enhancer activity and reconstruction of trajectories of cellular differentiation. In basal cell carcinoma, application of scATAC-seq reveals regulatory networks in malignant, stromal and immune cells in the tumor microenvironment. Analysis of scATAC-seq profiles from serial tumor biopsies before and after programmed cell death protein 1 blockade identifies chromatin regulators of therapy-responsive T cell subsets and reveals a shared regulatory program that governs intratumoral CD8 + T cell exhaustion and CD4 + T follicular helper cell development. We anticipate that scATAC-seq will enable the unbiased discovery of gene regulatory factors across diverse biological systems.

Modular domains of long non-coding RNAs can serve as scaffolds to bring distant regions of the linear genome into spatial proximity. Here, we present HiChIRP, a method leveraging bio-orthogonal chemistry and optimized chromosome conformation capture conditions, which enables interrogation of chromatin architecture focused around a specific RNA of interest down to approximately ten copies per cell. HiChIRP of three nuclear RNAs reveals insights into promoter interactions (7SK), telomere biology (telomerase RNA component) and inflammatory gene regulation (lincRNA-EPS).

Chromatin loops enable transcription-factor-bound distal enhancers to interact with their target promoters to regulate transcriptional programs. Although developmental transcription factors such as active forms of Notch can directly stimulate transcription by activating enhancers, the effect of their oncogenic subversion on the 3D organization of cancer genomes is largely undetermined. By mapping chromatin looping genome-wide in Notch-dependent triple-negative breast cancer and B cell lymphoma, we show that beyond the well-characterized role of Notch as an activator of distal enhancers, Notch regulates its direct target genes by instructing enhancer repositioning. Moreover, a large fraction of Notch-instructed regulatory loops form highly interacting enhancer and promoter spatial clusters termed "3D cliques." Loss- and gain-of-function experiments show that Notch preferentially targets hyperconnected 3D cliques that regulate the expression of crucial proto-oncogenes. Our observations suggest that oncogenic hijacking of developmental transcription factors can dysregulate transcription through widespread effects on the spatial organization of cancer genomes., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The vast majority of polymorphisms for human dermatologic diseases fall in noncoding DNA regions, leading to difficulty interpreting their functional significance. Recent work using chromosome conformation capture technology in combination with chromatin immunoprecipitation (ChIP) has provided a systematic means of linking noncoding variants in active enhancer loci to putative gene targets. Here, we apply H3K27ac HiChIP high-resolution contact maps, generated from primary human T-cell subsets (CD4 + naïve, T helper type 17, and regulatory T cells), to 21 dermatologic conditions associated with single nucleotide polymorphisms from 106 genome-wide association studies. This "enhancer connectome" identified 1,492 HiChIP gene targets from 542 noncoding SNPs (P ≤ 5.0 × 10 -8 ). SNP-containing enhancers from inflammatory skin conditions were significantly enriched at the HLA locus and also targeted several key factors from the JAK-STAT signaling pathway, but nonimmune conditions were not. A focused profiling of systemic lupus erythematosus HiChIP genes identified enhancer interactions with factors important for effector CD4 + T-cell differentiation and function, including IRF8 and members of the Ikaros family of zinc-finger proteins. Our results show the ability of the enhancer connectome to nominate functionally relevant candidates from genome-wide association study-identified variants, representing a powerful tool to guide future studies into the genomic regulatory mechanisms underlying dermatologic diseases., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)

Regulatory T (Treg) cells maintain immune tolerance through the master transcription factor forkhead box P3 (FOXP3), which is crucial for Treg cell function and homeostasis. We identified an IPEX (immune dysregulation polyendocrinopathy enteropathy X-linked) syndrome patient with a FOXP3 mutation in the domain swap interface of the protein. Recapitulation of this Foxp3 variant in mice led to the development of an autoimmune syndrome consistent with an unrestrained T helper type 2 (Th2) immune response. Genomic analysis of Treg cells by RNA-sequencing, Foxp3 chromatin immunoprecipitation followed by high-throughput DNA sequencing (ChIP-sequencing), and H3K27ac-HiChIP revealed a specific de-repression of the Th2 transcriptional program leading to the generation of Th2-like Treg cells that were unable to suppress extrinsic Th2 cells. Th2-like Treg cells showed increased intra-chromosomal interactions in the Th2 locus, leading to type 2 cytokine production. These findings identify a direct role for Foxp3 in suppressing Th2-like Treg cells and implicate additional pathways that could be targeted to restrain Th2 trans-differentiated Treg cells., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Here, we present Perturb-ATAC, a method that combines multiplexed CRISPR interference or knockout with genome-wide chromatin accessibility profiling in single cells based on the simultaneous detection of CRISPR guide RNAs and open chromatin sites by assay of transposase-accessible chromatin with sequencing (ATAC-seq). We applied Perturb-ATAC to transcription factors (TFs), chromatin-modifying factors, and noncoding RNAs (ncRNAs) in ∼4,300 single cells, encompassing more than 63 genotype-phenotype relationships. Perturb-ATAC in human B lymphocytes uncovered regulators of chromatin accessibility, TF occupancy, and nucleosome positioning and identified a hierarchy of TFs that govern B cell state, variation, and disease-associated cis-regulatory elements. Perturb-ATAC in primary human epidermal cells revealed three sequential modules of cis-elements that specify keratinocyte fate. Combinatorial deletion of all pairs of these TFs uncovered their epistatic relationships and highlighted genomic co-localization as a basis for synergistic interactions. Thus, Perturb-ATAC is a powerful strategy to dissect gene regulatory networks in development and disease., (Copyright © 2018 Elsevier Inc. All rights reserved.)

We present the genome-wide chromatin accessibility profiles of 410 tumor samples spanning 23 cancer types from The Cancer Genome Atlas (TCGA). We identify 562,709 transposase-accessible DNA elements that substantially extend the compendium of known cis-regulatory elements. Integration of ATAC-seq (the assay for transposase-accessible chromatin using sequencing) with TCGA multi-omic data identifies a large number of putative distal enhancers that distinguish molecular subtypes of cancers, uncovers specific driving transcription factors via protein-DNA footprints, and nominates long-range gene-regulatory interactions in cancer. These data reveal genetic risk loci of cancer predisposition as active DNA regulatory elements in cancer, identify gene-regulatory interactions underlying cancer immune evasion, and pinpoint noncoding mutations that drive enhancer activation and may affect patient survival. These results suggest a systematic approach to understanding the noncoding genome in cancer to advance diagnosis and therapy., (Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Distinguishing classical dendritic cells from other myeloid cell types is complicated by the shared expression of cell surface markers. ZBTB46 is a zinc finger and BTB domain-containing transcription factor, which is expressed by dendritic cells and committed dendritic cell precursors, but not by plasmacytoid dendritic cells, monocytes, macrophages, or other immune cell populations. In this study, we demonstrate that expression of ZBTB46 identifies human dendritic cell neoplasms. We examined ZBTB46 expression in a range of benign and malignant histiocytic disorders and found that ZBTB46 is able to clearly define the dendritic cell identity of many previously unclassified histiocytic disease subtypes. In particular, all examined cases of Langerhans cell histiocytosis and histiocytic sarcoma expressed ZBTB46, while all cases of blastic plasmacytoid dendritic cell neoplasm, chronic myelomonocytic leukemia, juvenile xanthogranuloma, Rosai-Dorfman disease, and Erdheim-Chester disease failed to demonstrate expression of ZBTB46. Moreover, ZBTB46 expression clarified the identity of diagnostically challenging neoplasms, such as cases of indeterminate cell histiocytosis, classifying a fraction of these entities as dendritic cell malignancies. These findings clarify the lineage origins of human histiocytic disorders and distinguish dendritic cell disorders from all other myeloid neoplasms.

In this Letter, analysis of steady-state regulatory T (Treg) cell percentages from Il2ra enhancer deletion (EDEL) and wild-type (WT) mice revealed no differences between them (Extended Data Fig. 9d). This analysis included two mice whose genotypes were incorrectly assigned. Even after correction of the genotypes, no significant differences in Treg cell percentages were seen when data across experimental cohorts were averaged (as was done in Extended Data Fig. 9d). However, if we normalize the corrected data to account for variation among experimental cohorts, a subtle decrease in EDEL Treg cell percentages is revealed and, using the corrected and normalized data, we have redrawn Extended Data Fig. 9d in Supplementary Fig. 1. The Supplementary Information to this Amendment contains the corrected and reanalysed Extended Data Fig. 9d. The sentence "This enhancer deletion (EDEL) strain also had no obvious T cell phenotypes at steady state (Extended Data Fig. 9)." should read: "This enhancer deletion (EDEL) strain had a small decrease in the percentage of Treg cells (Extended Data Fig. 9).". This error does not affect any of the main figures in the Letter or the data from mice with the human autoimmune-associated single nucleotide polymorphism (SNP) knocked in or with a 12-base-pair deletion at the site (12DEL). In addition, we stated in the Methods that we observed consistent immunophenotypes of EDEL mice across three founders, but in fact, we observed consistent phenotypes in mice from two founders. This does not change any of our conclusions and the original Letter has not been corrected.

Noncoding mutations in cancer genomes are frequent but challenging to interpret. PVT1 encodes an oncogenic lncRNA, but recurrent translocations and deletions in human cancers suggest alternative mechanisms. Here, we show that the PVT1 promoter has a tumor-suppressor function that is independent of PVT1 lncRNA. CRISPR interference of PVT1 promoter enhances breast cancer cell competition and growth in vivo. The promoters of the PVT1 and the MYC oncogenes, located 55 kb apart on chromosome 8q24, compete for engagement with four intragenic enhancers in the PVT1 locus, thereby allowing the PVT1 promoter to regulate pause release of MYC transcription. PVT1 undergoes developmentally regulated monoallelic expression, and the PVT1 promoter inhibits MYC expression only from the same chromosome via promoter competition. Cancer genome sequencing identifies recurrent mutations encompassing the human PVT1 promoter, and genome editing verified that PVT1 promoter mutation promotes cancer cell growth. These results highlight regulatory sequences of lncRNA genes as potential disease-associated DNA elements., (Copyright © 2018 Elsevier Inc. All rights reserved.)

T cells create vast amounts of diversity in the genes that encode their T cell receptors (TCRs), which enables individual clones to recognize specific peptide-major histocompatibility complex (MHC) ligands. Here we combined sequencing of the TCR-encoding genes with assay for transposase-accessible chromatin with sequencing (ATAC-seq) analysis at the single-cell level to provide information on the TCR specificity and epigenomic state of individual T cells. By using this approach, termed transcript-indexed ATAC-seq (T-ATAC-seq), we identified epigenomic signatures in immortalized leukemic T cells, primary human T cells from healthy volunteers and primary leukemic T cells from patient samples. In peripheral blood CD4 + T cells from healthy individuals, we identified cis and trans regulators of naive and memory T cell states and found substantial heterogeneity in surface-marker-defined T cell populations. In patients with a leukemic form of cutaneous T cell lymphoma, T-ATAC-seq enabled identification of leukemic and nonleukemic regulatory pathways in T cells from the same individual by allowing separation of the signals that arose from the malignant clone from the background T cell noise. Thus, T-ATAC-seq is a new tool that enables analysis of epigenomic landscapes in clonal T cells and should be valuable for studies of T cell malignancy, immunity and immunotherapy.

The challenge of linking intergenic mutations to target genes has limited molecular understanding of human diseases. Here we show that H3K27ac HiChIP generates high-resolution contact maps of active enhancers and target genes in rare primary human T cell subtypes and coronary artery smooth muscle cells. Differentiation of naive T cells into T helper 17 cells or regulatory T cells creates subtype-specific enhancer-promoter interactions, specifically at regions of shared DNA accessibility. These data provide a principled means of assigning molecular functions to autoimmune and cardiovascular disease risk variants, linking hundreds of noncoding variants to putative gene targets. Target genes identified with HiChIP are further supported by CRISPR interference and activation at linked enhancers, by the presence of expression quantitative trait loci, and by allele-specific enhancer loops in patient-derived primary cells. The majority of disease-associated enhancers contact genes beyond the nearest gene in the linear genome, leading to a fourfold increase in the number of potential target genes for autoimmune and cardiovascular diseases.

Chromosome conformation is an important feature of metazoan gene regulation; however, enhancer-promoter contact remodeling during cellular differentiation remains poorly understood. To address this, genome-wide promoter capture Hi-C (CHi-C) was performed during epidermal differentiation. Two classes of enhancer-promoter contacts associated with differentiation-induced genes were identified. The first class ('gained') increased in contact strength during differentiation in concert with enhancer acquisition of the H3K27ac activation mark. The second class ('stable') were pre-established in undifferentiated cells, with enhancers constitutively marked by H3K27ac. The stable class was associated with the canonical conformation regulator cohesin, whereas the gained class was not, implying distinct mechanisms of contact formation and regulation. Analysis of stable enhancers identified a new, essential role for a constitutively expressed, lineage-restricted ETS-family transcription factor, EHF, in epidermal differentiation. Furthermore, neither class of contacts was observed in pluripotent cells, suggesting that lineage-specific chromatin structure is established in tissue progenitor cells and is further remodeled in terminal differentiation.

We present Omni-ATAC, an improved ATAC-seq protocol for chromatin accessibility profiling that works across multiple applications with substantial improvement of signal-to-background ratio and information content. The Omni-ATAC protocol generates chromatin accessibility profiles from archival frozen tissue samples and 50-μm sections, revealing the activities of disease-associated DNA elements in distinct human brain structures. The Omni-ATAC protocol enables the interrogation of personal regulomes in tissue context and translational studies.

The majority of genetic variants associated with common human diseases map to enhancers, non-coding elements that shape cell-type-specific transcriptional programs and responses to extracellular cues. Systematic mapping of functional enhancers and their biological contexts is required to understand the mechanisms by which variation in non-coding genetic sequences contributes to disease. Functional enhancers can be mapped by genomic sequence disruption, but this approach is limited to the subset of enhancers that are necessary in the particular cellular context being studied. We hypothesized that recruitment of a strong transcriptional activator to an enhancer would be sufficient to drive target gene expression, even if that enhancer was not currently active in the assayed cells. Here we describe a discovery platform that can identify stimulus-responsive enhancers for a target gene independent of stimulus exposure. We used tiled CRISPR activation (CRISPRa) to synthetically recruit a transcriptional activator to sites across large genomic regions (more than 100 kilobases) surrounding two key autoimmunity risk loci, CD69 and IL2RA. We identified several CRISPRa-responsive elements with chromatin features of stimulus-responsive enhancers, including an IL2RA enhancer that harbours an autoimmunity risk variant. Using engineered mouse models, we found that sequence perturbation of the disease-associated Il2ra enhancer did not entirely block Il2ra expression, but rather delayed the timing of gene activation in response to specific extracellular signals. Enhancer deletion skewed polarization of naive T cells towards a pro-inflammatory T helper (T H 17) cell state and away from a regulatory T cell state. This integrated approach identifies functional enhancers and reveals how non-coding variation associated with human immune dysfunction alters context-specific gene programs.

Genome conformation is central to gene control but challenging to interrogate. Here we present HiChIP, a protein-centric chromatin conformation method. HiChIP improves the yield of conformation-informative reads by over 10-fold and lowers the input requirement over 100-fold relative to that of ChIA-PET. HiChIP of cohesin reveals multiscale genome architecture with greater signal-to-background ratios than those of in situ Hi-C.

icSHAPE (in vivo click selective 2-hydroxyl acylation and profiling experiment) captures RNA secondary structure at a transcriptome-wide level by measuring nucleotide flexibility at base resolution. Living cells are treated with the icSHAPE chemical NAI-N3 followed by selective chemical enrichment of NAI-N3-modified RNA, which provides an improved signal-to-noise ratio compared with similar methods leveraging deep sequencing. Purified RNA is then reverse-transcribed to produce cDNA, with SHAPE-modified bases leading to truncated cDNA. After deep sequencing of cDNA, computational analysis yields flexibility scores for every base across the starting RNA population. The entire experimental procedure can be completed in ∼5 d, and the sequencing and bioinformatics data analysis take an additional 4-5 d with no extensive computational skills required. Comparing in vivo and in vitro icSHAPE measurements can reveal in vivo RNA-binding protein imprints or facilitate the dissection of RNA post-transcriptional modifications. icSHAPE reactivities can additionally be used to constrain and improve RNA secondary structure prediction models.

Shade from neighboring plants limits light for photosynthesis; as a consequence, plants have a variety of strategies to avoid canopy shade and compete with their neighbors for light. Collectively the response to foliar shade is called the shade avoidance syndrome (SAS). The SAS includes elongation of a variety of organs, acceleration of flowering time, and additional physiological responses, which are seen throughout the plant life cycle. However, current mechanistic knowledge is mainly limited to shade-induced elongation of seedlings. Here we use phenotypic profiling of seedling, leaf, and flowering time traits to untangle complex SAS networks. We used over-representation analysis (ORA) of shade-responsive genes, combined with previous annotation, to logically select 59 known and candidate novel mutants for phenotyping. Our analysis reveals shared and separate pathways for each shade avoidance response. In particular, auxin pathway components were required for shade avoidance responses in hypocotyl, petiole, and flowering time, whereas jasmonic acid pathway components were only required for petiole and flowering time responses. Our phenotypic profiling allowed discovery of seventeen novel shade avoidance mutants. Our results demonstrate that logical selection of mutants increased success of phenotypic profiling to dissect complex traits and discover novel components.

Protein expression is regulated by the production and degradation of messenger RNAs (mRNAs) and proteins, but their specific relationships remain unknown. We combine measurements of protein production and degradation and mRNA dynamics so as to build a quantitative genomic model of the differential regulation of gene expression in lipopolysaccharide-stimulated mouse dendritic cells. Changes in mRNA abundance play a dominant role in determining most dynamic fold changes in protein levels. Conversely, the preexisting proteome of proteins performing basic cellular functions is remodeled primarily through changes in protein production or degradation, accounting for more than half of the absolute change in protein molecules in the cell. Thus, the proteome is regulated by transcriptional induction for newly activated cellular functions and by protein life-cycle changes for remodeling of preexisting functions., (Copyright © 2015, American Association for the Advancement of Science.)

Pseudouridine is the most abundant RNA modification, yet except for a few well-studied cases, little is known about the modified positions and their function(s). Here, we develop Ψ-seq for transcriptome-wide quantitative mapping of pseudouridine. We validate Ψ-seq with spike-ins and de novo identification of previously reported positions and discover hundreds of unique sites in human and yeast mRNAs and snoRNAs. Perturbing pseudouridine synthases (PUS) uncovers which pseudouridine synthase modifies each site and their target sequence features. mRNA pseudouridinylation depends on both site-specific and snoRNA-guided pseudouridine synthases. Upon heat shock in yeast, Pus7p-mediated pseudouridylation is induced at >200 sites, and PUS7 deletion decreases the levels of otherwise pseudouridylated mRNA, suggesting a role in enhancing transcript stability. rRNA pseudouridine stoichiometries are conserved but reduced in cells from dyskeratosis congenita patients, where the PUS DKC1 is mutated. Our work identifies an enhanced, transcriptome-wide scope for pseudouridine and methods to dissect its underlying mechanisms and function., (Copyright © 2014 Elsevier Inc. All rights reserved.)

N6-methyladenosine (m6A) is a common modification of mRNA with potential roles in fine-tuning the RNA life cycle. Here, we identify a dense network of proteins interacting with METTL3, a component of the methyltransferase complex, and show that three of them (WTAP, METTL14, and KIAA1429) are required for methylation. Monitoring m6A levels upon WTAP depletion allowed the definition of accurate and near single-nucleotide resolution methylation maps and their classification into WTAP-dependent and -independent sites. WTAP-dependent sites are located at internal positions in transcripts, topologically static across a variety of systems we surveyed, and inversely correlated with mRNA stability, consistent with a role in establishing "basal" degradation rates. WTAP-independent sites form at the first transcribed base as part of the cap structure and are present at thousands of sites, forming a previously unappreciated layer of transcriptome complexity. Our data shed light on the proteomic and transcriptional underpinnings of this RNA modification., (Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.)

N(6)-methyladenosine (m(6)A) is the most ubiquitous mRNA base modification, but little is known about its precise location, temporal dynamics, and regulation. Here, we generated genomic maps of m(6)A sites in meiotic yeast transcripts at nearly single-nucleotide resolution, identifying 1,308 putatively methylated sites within 1,183 transcripts. We validated eight out of eight methylation sites in different genes with direct genetic analysis, demonstrated that methylated sites are significantly conserved in a related species, and built a model that predicts methylated sites directly from sequence. Sites vary in their methylation profiles along a dense meiotic time course and are regulated both locally, via predictable methylatability of each site, and globally, through the core meiotic circuitry. The methyltransferase complex components localize to the yeast nucleolus, and this localization is essential for mRNA methylation. Our data illuminate a conserved, dynamically regulated methylation program in yeast meiosis and provide an important resource for studying the function of this epitranscriptomic modification., (Copyright © 2013 Elsevier Inc. All rights reserved.)

The presence of 5-methylcytidine (m(5)C) in tRNA and rRNA molecules of a wide variety of organisms was first observed more than 40 years ago. However, detection of this modification was limited to specific, abundant, RNA species, due to the usage of low-throughput methods. To obtain a high resolution, systematic, and comprehensive transcriptome-wide overview of m(5)C across the three domains of life, we used bisulfite treatment on total RNA from both gram positive (B. subtilis) and gram negative (E. coli) bacteria, an archaeon (S. solfataricus) and a eukaryote (S. cerevisiae), followed by massively parallel sequencing. We were able to recover most previously documented m(5)C sites on rRNA in the four organisms, and identified several novel sites in yeast and archaeal rRNAs. Our analyses also allowed quantification of methylated m(5)C positions in 64 tRNAs in yeast and archaea, revealing stoichiometric differences between the methylation patterns of these organisms. Molecules of tRNAs in which m(5)C was absent were also discovered. Intriguingly, we detected m(5)C sites within archaeal mRNAs, and identified a consensus motif of AUCGANGU that directs methylation in S. solfataricus. Our results, which were validated using m(5)C-specific RNA immunoprecipitation, provide the first evidence for mRNA modifications in archaea, suggesting that this mode of post-transcriptional regulation extends beyond the eukaryotic domain., Competing Interests: The authors have declared that no competing interests exist.

High throughput phenotyping (phenomics) is a powerful tool for linking genes to their functions (see review and recent examples). Leaves are the primary photosynthetic organ, and their size and shape vary developmentally and environmentally within a plant. For these reasons studies on leaf morphology require measurement of multiple parameters from numerous leaves, which is best done by semi-automated phenomics tools. Canopy shade is an important environmental cue that affects plant architecture and life history; the suite of responses is collectively called the shade avoidance syndrome (SAS). Among SAS responses, shade induced leaf petiole elongation and changes in blade area are particularly useful as indices. To date, leaf shape programs (e.g. SHAPE, LAMINA, LeafAnalyzer, LEAFPROCESSOR) can measure leaf outlines and categorize leaf shapes, but can not output petiole length. Lack of large-scale measurement systems of leaf petioles has inhibited phenomics approaches to SAS research. In this paper, we describe a newly developed ImageJ plugin, called LeafJ, which can rapidly measure petiole length and leaf blade parameters of the model plant Arabidopsis thaliana. For the occasional leaf that required manual correction of the petiole/leaf blade boundary we used a touch-screen tablet. Further, leaf cell shape and leaf cell numbers are important determinants of leaf size. Separate from LeafJ we also present a protocol for using a touch-screen tablet for measuring cell shape, area, and size. Our leaf trait measurement system is not limited to shade-avoidance research and will accelerate leaf phenotyping of many mutants and screening plants by leaf phenotyping.