Search

Eosinophilic esophagitis (EoE) is a rare atopic disorder associated with esophageal dysfunction, including difficulty swallowing, food impaction, and inflammation, that develops in a small subset of people with food allergies. Genome-wide association studies (GWASs) have identified 9 independent EoE risk loci reaching genome-wide significance (p < 5 × 10−8) and 27 additional loci of suggestive significance (5 × 10−8 < p < 1 × 10−5). In the current study, we perform linkage disequilibrium (LD) expansion of these loci to nominate a set of 531 variants that are potentially causal. To systematically interrogate the gene regulatory activity of these variants, we designed a massively parallel reporter assay (MPRA) containing the alleles of each variant within their genomic sequence context cloned into a GFP reporter library. Analysis of reporter gene expression in TE-7, HaCaT, and Jurkat cells revealed cell-type-specific gene regulation. We identify 32 allelic enhancer variants, representing 6 genome-wide significant EoE loci and 7 suggestive EoE loci, that regulate reporter gene expression in a genotype-dependent manner in at least one cellular context. By annotating these variants with expression quantitative trait loci (eQTL) and chromatin looping data in related tissues and cell types, we identify putative target genes affected by genetic variation in individuals with EoE. Transcription factor enrichment analyses reveal possible roles for cell-type-specific regulators, including GATA3. Our approach reduces the large set of EoE-associated variants to a set of 32 with allelic regulatory activity, providing functional insights into the effects of genetic variation in this disease. [Display omitted] Shook et al. developed and applied a massively parallel reporter assay for the allergic disease eosinophilic esophagitis (EoE). They identified common polymorphisms associated with EoE risk that also impacted transcriptional regulation. This study is an important step toward understanding how EoE genetic risk variants lead to genotype-dependent biology. [ABSTRACT FROM AUTHOR]

The most recent genome-wide association study (GWAS) of cutaneous melanoma identified 54 risk-associated loci, but functional variants and their target genes for most have not been established. Here, we performed massively parallel reporter assays (MPRAs) by using malignant melanoma and normal melanocyte cells and further integrated multi-layer annotation to systematically prioritize functional variants and susceptibility genes from these GWAS loci. Of 1,992 risk-associated variants tested in MPRAs, we identified 285 from 42 loci (78% of the known loci) displaying significant allelic transcriptional activities in either cell type (FDR < 1%). We further characterized MPRA-significant variants by motif prediction, epigenomic annotation, and statistical/functional fine-mapping to create integrative variant scores, which prioritized one to six plausible candidate variants per locus for the 42 loci and nominated a single variant for 43% of these loci. Overlaying the MPRA-significant variants with genome-wide significant expression or methylation quantitative trait loci (eQTLs or meQTLs, respectively) from melanocytes or melanomas identified candidate susceptibility genes for 60% of variants (172 of 285 variants). CRISPRi of top-scoring variants validated their cis -regulatory effect on the eQTL target genes, MAFF (22q13.1) and GPRC5A (12p13.1). Finally, we identified 36 melanoma-specific and 45 melanocyte-specific MPRA-significant variants, a subset of which are linked to cell-type-specific target genes. Analyses of transcription factor availability in MPRA datasets and variant-transcription-factor interaction in eQTL datasets highlighted the roles of transcription factors in cell-type-specific variant functionality. In conclusion, MPRAs along with variant scoring effectively prioritized plausible candidates for most melanoma GWAS loci and highlighted cellular contexts where the susceptibility variants are functional. Long et al. used massively parallel reporter assays combined with variant scoring to identify functional variants from 78% of known melanoma GWAS loci, including those specific to cell of origin versus cancer contexts. Linking prioritized functional variants to eQTLs identified target genes as validated by CRISPRi. [ABSTRACT FROM AUTHOR]

A meta-genome-wide association study across eight psychiatric disorders has highlighted the genetic architecture of pleiotropy in major psychiatric disorders. However, mechanisms underlying pleiotropic effects of the associated variants remain to be explored. We conducted a massively parallel reporter assay to decode the regulatory logic of variants with pleiotropic and disorder-specific effects. Pleiotropic variants differ from disorder-specific variants by exhibiting chromatin accessibility that extends across diverse cell types in the neuronal lineage and by altering motifs for transcription factors with higher connectivity in protein-protein interaction networks. We mapped pleiotropic and disorder-specific variants to putative target genes using functional genomics approaches and CRISPR perturbation. In vivo CRISPR perturbation of a pleiotropic and a disorder-specific gene suggests that pleiotropy may involve the regulation of genes expressed broadly across neuronal cell types and with higher network connectivity., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 Elsevier Inc. All rights reserved.)

In any given cell type, dozens of transcription factors (TFs) act in concert to control the activity of the genome by binding to specific DNA sequences in regulatory elements. Despite their considerable importance, we currently lack simple tools to directly measure the activity of many TFs in parallel. Massively parallel reporter assays (MPRAs) allow the detection of TF activities in a multiplexed fashion; however, we lack basic understanding to rationally design sensitive reporters for many TFs. Here, we use an MPRA to systematically optimize transcriptional reporters for 86 TFs and evaluate the specificity of all reporters across a wide array of TF perturbation conditions. We thus identified critical TF reporter design features and obtained highly sensitive and specific reporters for 62 TFs, many of which outperform available reporters. The resulting collection of "prime" TF reporters can be used to uncover TF regulatory networks and to illuminate signaling pathways. A record of this paper's transparent peer review process is included in the supplemental information., Competing Interests: Declaration of interests The authors declare that they have filed a patent application to secure intellectual property rights for the designed TF reporters. C.R. is a co-founder and shareholder of Metric Biotechnologies, Inc., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Our study investigated the underlying mechanism for the 14q24 renal cell carcinoma (RCC) susceptibility risk locus identified by a genome-wide association study (GWAS). The sentinel single-nucleotide polymorphism (SNP), rs4903064, at 14q24 confers an allele-specific effect on expression of the double PHD fingers 3 (DPF3) of the BAF SWI/SNF complex as assessed by massively parallel reporter assay, confirmatory luciferase assays, and eQTL analyses. Overexpression of DPF3 in renal cell lines increases growth rates and alters chromatin accessibility and gene expression, leading to inhibition of apoptosis and activation of oncogenic pathways. siRNA interference of multiple DPF3-deregulated genes reduces growth. Our results indicate that germline variation in DPF3 , a component of the BAF complex, part of the SWI/SNF complexes, can lead to reduced apoptosis and activation of the STAT3 pathway, both critical in RCC carcinogenesis. In addition, we show that altered DPF3 expression in the 14q24 RCC locus could influence the effectiveness of immunotherapy treatment for RCC by regulating tumor cytokine secretion and immune cell activation. [ABSTRACT FROM AUTHOR]

Using machine learning (ML), we interrogated the function of all human-chimpanzee variants in 2,645 human accelerated regions (HARs), finding 43% of HARs have variants with large opposing effects on chromatin state and 14% on neurodevelopmental enhancer activity. This pattern, consistent with compensatory evolution, was confirmed using massively parallel reporter assays in chimpanzee and human neural progenitor cells. The species-specific enhancer activity of HARs was accurately predicted from the presence and absence of transcription factor footprints in each species. Despite these striking cis effects, activity of a given HAR sequence was nearly identical in human and chimpanzee cells. This suggests that HARs did not evolve to compensate for changes in the trans environment but instead altered their ability to bind factors present in both species. Thus, ML prioritized variants with functional effects on human neurodevelopment and revealed an unexpected reason why HARs may have evolved so rapidly. • 31% of HARs act as enhancers in chimpanzee and human neural progenitor cells • 43% of HARs have human variants with large opposing effects on chromatin state • Transcription factor footprints predict human:chimpanzee enhancer activity • Rapid evolution of HAR sequences reflects compensatory tuning of enhancer activity Whalen et al. couple deep learning with functional assays in chimpanzee and human cells to interrogate the neurodevelopmental enhancer potential of 2,645 human accelerated regions (HARs). Activity is dominated by cis rather than trans effects, and compensatory changes are identified as a driver of rapid HAR evolution. [ABSTRACT FROM AUTHOR]

The notochord is a defining feature of all chordates. The transcription factors Zic and ETS regulate enhancer activity within the notochord. We conduct high-throughput screens of genomic elements within developing Ciona embryos to understand how Zic and ETS sites encode notochord activity. Our screen discovers an enhancer located near Lama, a gene critical for notochord development. Reversing the orientation of an ETS site within this enhancer abolishes expression, indicating that enhancer grammar is critical for notochord activity. Similarly organized clusters of Zic and ETS sites occur within mouse and human Lama1 introns. Within a Brachyury (Bra) enhancer, FoxA and Bra, in combination with Zic and ETS binding sites, are necessary and sufficient for notochord expression. This binding site logic also occurs within other Ciona and vertebrate Bra enhancers. Collectively, this study uncovers the importance of grammar within notochord enhancers and discovers signatures of enhancer logic and grammar conserved across chordates., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Gene expression is in part controlled by cis -regulatory elements (CREs) such as enhancers and repressive elements. Anecdotal evidence has indicated that a CRE and a promoter need to be biochemically compatible for promoter regulation to occur, but this compatibility has remained poorly characterized in mammalian cells. We used high-throughput combinatorial reporter assays to test thousands of CRE-promoter pairs from three Mb-sized genomic regions in mouse cells. This revealed that CREs vary substantially in their promoter compatibility, ranging from striking specificity to broad promiscuity. More than half of the tested CREs exhibit significant promoter selectivity. Housekeeping promoters tend to have similar CRE preferences, but other promoters exhibit a wide diversity of compatibilities. Higher-order transcription factors (TF) motif combinations may account for compatibility. CRE-promoter selectivity does not correlate with looping interactions in the native genomic context, suggesting that chromatin folding and compatibility are two orthogonal mechanisms that confer specificity to gene regulation. [Display omitted] • Intrinsic compatibility of thousands of enhancer-promoter combinations was tested • Compatibilities exhibit a broad spectrum, from promiscuous to highly specific • Enhancer-promoter compatibility and chromatin looping appear to be independent A major question in genome biology is how enhancers "choose" their target promoter. Using high-throughput reporter assays, Martinez-Ara and colleagues surveyed the intrinsic compatibilities of thousands of enhancer-promoter pairs in mouse cells. They found a wide diversity of specificities, mostly likely driven by a complex grammar of sequence motifs. [ABSTRACT FROM AUTHOR]

3′ untranslated region (3′UTR) variants are strongly associated with human traits and diseases, yet few have been causally identified. We developed the massively parallel reporter assay for 3′UTRs (MPRAu) to sensitively assay 12,173 3′UTR variants. We applied MPRAu to six human cell lines, focusing on genetic variants associated with genome-wide association studies (GWAS) and human evolutionary adaptation. MPRAu expands our understanding of 3′UTR function, suggesting that simple sequences predominately explain 3′UTR regulatory activity. We adapt MPRAu to uncover diverse molecular mechanisms at base pair resolution, including an adenylate-uridylate (AU)-rich element of LEPR linked to potential metabolic evolutionary adaptations in East Asians. We nominate hundreds of 3′UTR causal variants with genetically fine-mapped phenotype associations. Using endogenous allelic replacements, we characterize one variant that disrupts a miRNA site regulating the viral defense gene TRIM14 and one that alters PILRB abundance, nominating a causal variant underlying transcriptional changes in age-related macular degeneration. [Display omitted] • Assayed thousands of GWAS and adaptation associated 3′UTR variants in 6 cell lines • Nominated hundreds of causal GWAS variants with functional evidence of activity • Characterized mechanistic regulatory motifs at base pair resolution • Used allelic replacement on causal variants for macular degeneration and viral defense Massively parallel reporter assay for 3′UTRs measures individual regulatory effects of over 12,000 3′UTR variants associated with human disease and evolutionary selection in many cell types, nominating functional genetic variation. [ABSTRACT FROM AUTHOR]

Hematopoietic stem cell gene therapy is a promising approach for treating disorders of the hematopoietic system. Identifying combinations of cis -regulatory elements that do not impede packaging or transduction efficiency when included in lentiviral vectors has proven challenging. In this study, we deploy LV-MPRA (lentiviral vector-based, massively parallel reporter assay), an approach that simultaneously analyzes thousands of synthetic DNA fragments in parallel to identify sequence-intrinsic and lineage-specific enhancer function at near-base-pair resolution. We demonstrate the power of LV-MPRA in elucidating the boundaries of previously unknown intrinsic enhancer sequences of the human β-globin locus control region. Our approach facilitated the rapid assembly of novel therapeutic β AS3 -globin lentiviral vectors harboring strong lineage-specific recombinant control elements capable of correcting a mouse model of sickle cell disease. LV-MPRA can be used to map any genomic locus for enhancer activity and facilitates the rapid development of therapeutic vectors for treating disorders of the hematopoietic system or other specific tissues and cell types., (© 2020 The Author(s).)

Alternative polyadenylation (APA) is a major driver of transcriptome diversity in human cells. Here, we use deep learning to predict APA from DNA sequence alone. We trained our model (APARENT, APA REgression NeT) on isoform expression data from over 3 million APA reporters. APARENT's predictions are highly accurate when tasked with inferring APA in synthetic and human 3′UTRs. Visualizing features learned across all network layers reveals that APARENT recognizes sequence motifs known to recruit APA regulators, discovers previously unknown sequence determinants of 3′ end processing, and integrates these features into a comprehensive, interpretable, cis -regulatory code. We apply APARENT to forward engineer functional polyadenylation signals with precisely defined cleavage position and isoform usage and validate predictions experimentally. Finally, we use APARENT to quantify the impact of genetic variants on APA. Our approach detects pathogenic variants in a wide range of disease contexts, expanding our understanding of the genetic origins of disease. • Trained a neural network to predict APA using data from over 3 million reporters • Visualized learned features to reveal a rich cis -regulatory code for APA • Developed and tested an algorithm to accurately engineer polyadenylation signals • Predicted and experimentally characterized over 12,000 human APA variants A deep neural network enables precise engineering of polyadenylation signals, identifies human genetic variants that act through mis-regulating APA, and learns a comprehensive model of the cis -regulatory APA code. [ABSTRACT FROM AUTHOR]

G protein-coupled receptors (GPCRs) are central to how mammalian cells sense and respond to chemicals. Mammalian olfactory receptors (ORs), the largest family of GPCRs, mediate the sense of smell through activation by small molecules, though for most bonafide ligands, they have not been identified. Here, we introduce a platform to screen large chemical panels against multiplexed GPCR libraries using next-generation sequencing of barcoded genetic reporters in stably engineered human cell lines. We mapped 39 mammalian ORs against 181 odorants and identified 79 interactions that have not been reported to our knowledge, including ligands for 15 previously orphaned receptors. This multiplexed receptor assay allows the cost-effective mapping of large chemical libraries to receptor repertoires at scale., (Copyright © 2019. Published by Elsevier Inc.)