Start codon variant in LAG3 is associated with decreased LAG-3 expression and increased risk of autoimmune thyroid disease.

Autoimmune thyroid disease (AITD) is a common autoimmune disease. In a GWAS meta-analysis of 110,945 cases and 1,084,290 controls, 290 sequence variants at 225 loci are associated with AITD. Of these variants, 115 are previously unreported. Multiomics analysis yields 235 candidate genes outside the MHC-region and the findings highlight the importance of genes involved in T-cell regulation. A rare 5'-UTR variant (rs781745126-T, MAF = 0.13% in Iceland) in LAG3 has the largest effect (OR = 3.42, P = 2.2 × 10-16) and generates a novel start codon for an open reading frame upstream of the canonical protein translation initiation site. rs781745126-T reduces mRNA and surface expression of the inhibitory immune checkpoint LAG-3 co-receptor on activated lymphocyte subsets and halves LAG-3 levels in plasma among heterozygotes. All three homozygous carriers of rs781745126-T have AITD, of whom one also has two other T-cell mediated diseases, that is vitiligo and type 1 diabetes. rs781745126-T associates nominally with vitiligo (OR = 5.1, P = 6.5 × 10-3) but not with type 1 diabetes. Thus, the effect of rs781745126-T is akin to drugs that inhibit LAG-3, which unleash immune responses and can have thyroid dysfunction and vitiligo as adverse events. This illustrates how a multiomics approach can reveal potential drug targets and safety concerns.

eRNA co-expression network uncovers TF dependency and convergent cooperativity.

Enhancer RNAs (eRNAs) are non-coding RNAs produced by transcriptional enhancers that are highly correlated with their activity. Using a capped nascent RNA sequencing (PRO-cap) dataset in human lymphoblastoid cell lines across 67 individuals, we identified inter-individual variation in the expression of over 80 thousand transcribed transcriptional regulatory elements (tTREs), in both enhancers and promoters. Co-expression analysis of eRNAs from tTREs across individuals revealed how enhancers are associated with each other and with promoters. Mid- to long-range co-expression showed a distance-dependent decay that was modified by TF occupancy. In particular, we found a class of "bivalent" TFs, including Cohesin, that both facilitate and isolate the interaction between enhancers and/or promoters, depending on their topology. At short distances, we observed strand-specific correlations between nearby eRNAs in both convergent and divergent orientations. Our results support a cooperative model of convergent eRNAs, consistent with eRNAs facilitating adjacent enhancers rather than interfering with each other. Therefore, our approach to infer functional interactions from co-expression analyses provided novel insights into the principles of enhancer interactions as a function of distance, orientation, and binding landscapes of TFs.

Revealing eRNA interactions: TF dependency and convergent cooperativity.

Enhancer RNAs (eRNAs) are non-coding RNAs produced from transcriptional enhancers that are highly correlated with their activities. Using capped nascent RNA sequencing (PRO-cap) dataset in human lymphoblastoid cell lines across individuals, we identified inter-individual variation of expression in over 80 thousand transcribed transcriptional regulatory elements (tTREs), in both enhancers and promoters. Co-expression analysis of eRNAs from tTREs across individuals revealed how enhancers interact with each other and with promoters. Mid-to-long range interactions showed distance-dependent decay, which was modified by TF occupancy. In particular, we found a class of 'bivalent' TFs, including Cohesin, which both facilitates and insulates the interaction between enhancers and/or promoters depending on the topology. In short ranges, we observed strand specific interactions between nearby eRNAs in both convergent or divergent orientations. Our finding supports a cooperative convergent eRNA model, which is compatible with eRNA remodeling neighboring enhancers rather than interfering with each other. Therefore, our approach to infer functional interactions from co-expression analyses provided novel insights into the principles of enhancer interactions depending on the distance, orientation, and the binding landscapes of TFs.

Population-scale study of eRNA transcription reveals bipartite functional enhancer architecture.

Enhancer RNAs (eRNA) are unstable non-coding RNAs, transcribed bidirectionally from active regulatory sequences, whose expression levels correlate with enhancer activity. We use capped-nascent-RNA sequencing to efficiently capture bidirectional transcription initiation across several human lymphoblastoid cell lines (Yoruba population) and detect ~75,000 eRNA transcription sites with high sensitivity and specificity. The use of nascent-RNA sequencing sidesteps the confounding effect of eRNA instability. We identify quantitative trait loci (QTLs) associated with the level and directionality of eRNA expression. High-resolution analyses of these two types of QTLs reveal distinct positions of enrichment at the central transcription factor (TF) binding regions and at the flanking eRNA initiation regions, both of which are associated with mRNA expression QTLs. These two regions-the central TF-binding footprint and the eRNA initiation cores-define a bipartite architecture of enhancers, inform enhancer function, and can be used as an indicator of the significance of non-coding regulatory variants.

TED-Seq Identifies the Dynamics of Poly(A) Length during ER Stress.

Post-transcriptional RNA processing is a core mechanism of gene expression control in cell stress response. The poly(A) tail influences mRNA translation and stability, but it is unclear whether there are global roles of poly(A)-tail lengths in cell stress. To address this, we developed tail-end displacement sequencing (TED-seq) for an efficient transcriptome-wide profiling of poly(A) lengths and applied it to endoplasmic reticulum (ER) stress in human cells. ER stress induced increases in the poly(A) lengths of certain mRNAs, including known ER stress regulators, XBP1, DDIT3, and HSPA5. Importantly, the mRNAs with increased poly(A) lengths are both translationally de-repressed and stabilized. Furthermore, mRNAs in stress-induced RNA granules have shorter poly(A) tails than in the cytoplasm, supporting the view that RNA processing is compartmentalized. In conclusion, TED-seq reveals that poly(A) length is dynamically regulated upon ER stress, with potential consequences for both translation and mRNA turnover.

Systematic analysis of the Hmga2 3' UTR identifies many independent regulatory sequences and a novel interaction between distal sites.

The 3' untranslated regions (3' UTRs) of mRNAs regulate transcripts by serving as binding sites for regulatory factors, including microRNAs and RNA binding proteins. Binding of such trans-acting factors can control the rates of mRNA translation, decay, and other aspects of mRNA biology. To better understand the role of 3' UTRs in gene regulation, we performed a detailed analysis of a model mammalian 3' UTR, that of Hmga2, with the principal goals of identifying the complete set of regulatory elements within a single 3' UTR, and determining the extent to which elements interact with and affect one another. Hmga2 is an oncogene whose overexpression in cancers often stems from mutations that remove 3'-UTR regulatory sequences. We used reporter assays in cultured cells to generate maps of cis-regulatory information across the Hmga2 3' UTR at different resolutions, ranging from 50 to 400 nt. We found many previously unidentified regulatory sites, a large number of which were up-regulating. Importantly, the overall location and impact of regulatory sites was conserved between different species (mouse, human, and chicken). By systematically comparing the regulatory impact of 3'-UTR segments of different sizes we were able to determine that the majority of regulatory sequences function independently; only a very small number of segments showed evidence of any interactions. However, we discovered a novel interaction whereby terminal 3'-UTR sequences induced internal up-regulating elements to convert to repressive elements. By fully characterizing one 3' UTR, we hope to better understand the principles of 3'-UTR-mediated gene regulation.

Coordinate regulation of DNA methylation and H3K27me3 in mouse embryonic stem cells.

Chromatin is separated into functional domains distinguished by combinatorial patterns of post-translational histone modifications and DNA methylation. Recent studies examining multiple histone modifications have found numerous chromatin states with distinct profiles of chromatin marks and functional enrichments. There are data showing coordinate regulation between DNAme and H3K27me3, which are both involved in the establishment and maintenance of epigenetic gene silencing, but the data are conflicting. Multiple studies have presented evidence to support the theory that PRC2 and DNAme cooperate to achieve silencing, or alternatively that H3K27me3 and DNAme act antagonistically. Here we examine the effect loss of either PRC2 or DNA methyltransferase activity has on the placement of the reciprocal mark in mouse ES cells. We find that DNAme is acting globally to antagonize the placement of H3K27me3, in accordance with recently published results. At least 471,011 domains in the mouse genome acquire H3K27me3 when DNAme is diminished. Of these 466,563 have been shown to be fully methylated in wildtype ES cells, indicating the effects of DNAme on H3K27me3 are direct. In a reciprocal experiment, we examine the effect loss of PRC2 has on the placement of DNAme. In contrast to the global antagonism DNAme has on the placement of H3K27me3, loss of H3K27me3 has a modest effect on DNAme, with only 4% of genes undergoing changes in DNAme, including 861 showing increases and 552 showing losses of overall DNAme. We anticipate that integrating genomic datasets where the effect of loss of a particular epigenetic mark has on the placement of other marks will help elucidate the rules governing epigenetic regulation and what role coordinate regulation of epigenetic marks plays in development and disease.