Allele-specific transcription factor binding across human brain regions offers mechanistic insight into eQTLs.

Transcription factors (TFs) regulate gene expression by facilitating or disrupting the formation of transcription initiation machinery at particular genomic loci. Because TF occupancy is driven in part by recognition of DNA sequence, genetic variation can influence TF-DNA associations and gene regulation. To identify variants that impact TF binding in human brain tissues, we assessed allele-specific binding (ASB) at heterozygous variants for 94 TFs in nine brain regions from two donors. Leveraging graph genomes constructed from phased genomic sequence data, we compared ChIP-seq signals between alleles at heterozygous variants within each brain region and identified thousands of variants exhibiting ASB for at least one TF. ASB reproducibility was measured by comparisons between independent experiments both within and between donors. We found that rare alleles in the general population more frequently led to reduced TF binding, whereas common alleles had an equal likelihood of increasing or decreasing binding. Further, for ASB variants in predicted binding motifs, the favored allele tended to be the one with the stronger expected motif match, but this concordance was not observed within highly occupied sites. We also found that neuron-specific cis-regulatory elements (cCREs), in contrast with oligodendrocyte-specific cCREs, showed depletion of ASB variants. We identified 2670 ASB variants associated with evidence for allele-specific gene expression in the brain from GTEx data and observed increasing eQTL effect direction concordance as ASB significance increases. These results provide a valuable and unique resource for mechanistic analysis of cis-regulatory variation in human brain tissue.

Multiomic profiling of transcription factor binding and function in human brain.

Transcription factors (TFs) orchestrate gene expression programs crucial for brain function, but we lack detailed information about TF binding in human brain tissue. We generated a multiomic resource (ChIP-seq, ATAC-seq, RNA-seq, DNA methylation) on bulk tissues and sorted nuclei from several postmortem brain regions, including binding maps for more than 100 TFs. We demonstrate improved measurements of TF activity, including motif recognition and gene expression modeling, upon identification and removal of high TF occupancy regions. Further, predictive TF binding models demonstrate a bias for these high-occupancy sites. Neuronal TFs SATB2 and TBR1 bind unique regions depleted for such sites and promote neuronal gene expression. Binding sites for TFs, including TBR1 and PKNOX1, are enriched for risk variants associated with neuropsychiatric disorders, predominantly in neurons. This work, titled BrainTF, is a powerful resource for future studies seeking to understand the roles of specific TFs in regulating gene expression in the human brain.

Probabilistic association of differentially expressed genes with cis-regulatory elements.

Differential gene expression in response to perturbations is mediated at least in part by changes in binding of transcription factors (TFs) and other proteins at specific genomic regions. Association of these cis-regulatory elements (CREs) with their target genes is a challenging task that is essential to address many biological and mechanistic questions. Many current approaches rely on chromatin conformation capture techniques or single-cell correlational methods to establish CRE-to-gene associations. These methods can be effective but have limitations, including resolution, gaps in detectable association distances, and cost. As an alternative, we have developed DegCre, a nonparametric method that evaluates correlations between measurements of perturbation-induced differential gene expression and differential regulatory signal at CREs to score possible CRE-to-gene associations. It has several unique features, including the ability to use any type of CRE activity measurement, yield probabilistic scores for CRE-to-gene pairs, and assess CRE-to-gene pairings across a wide range of sequence distances. We apply DegCre to six data sets, each using different perturbations and containing a variety of regulatory signal measurements, including chromatin openness, histone modifications, and TF occupancy. To test their efficacy, we compare DegCre associations to Hi-C loop calls and CRISPR-validated CRE-to-gene associations, establishing good performance by DegCre that is comparable or superior to competing methods. DegCre is a novel approach to the association of CREs to genes from a perturbation-differential perspective, with strengths that are complementary to existing approaches and allow for new insights into gene regulation.

Neuronal MAPT expression is mediated by long-range interactions with cis-regulatory elements.

Tauopathies are a group of neurodegenerative diseases defined by abnormal aggregates of tau, a microtubule-associated protein encoded by MAPT. MAPT expression is near absent in neural progenitor cells (NPCs) and increases during differentiation. This temporally dynamic expression pattern suggests that MAPT expression could be controlled by transcription factors and cis-regulatory elements specific to differentiated cell types. Given the relevance of MAPT expression to neurodegeneration pathogenesis, identification of such elements is relevant to understanding disease risk and pathogenesis. Here, we performed chromatin conformation assays (HiC & Capture-C), single-nucleus multiomics (RNA-seq+ATAC-seq), bulk ATAC-seq, and ChIP-seq for H3K27ac and CTCF in NPCs and differentiated neurons to nominate candidate cis-regulatory elements (cCREs). We assayed these cCREs using luciferase assays and CRISPR interference (CRISPRi) experiments to measure their effects on MAPT expression. Finally, we integrated cCRE annotations into an analysis of genetic variation in neurodegeneration-affected individuals and control subjects. We identified both proximal and distal regulatory elements for MAPT and confirmed the regulatory function for several regions, including three regions centromeric to MAPT beyond the H1/H2 haplotype inversion breakpoint. We also found that rare and predicted damaging genetic variation in nominated CREs was nominally depleted in dementia-affected individuals relative to control subjects, consistent with the hypothesis that variants that disrupt MAPT enhancer activity, and thereby reduced MAPT expression, may be protective against neurodegenerative disease. Overall, this study provides compelling evidence for pursuing detailed knowledge of CREs for genes of interest to permit better understanding of disease risk.

Allele biased transcription factor binding across human brain regions gives mechanistic insight into eQTLs.

Transcription Factors (TFs) influence gene expression by facilitating or disrupting the formation of transcription initiation machinery at particular genomic loci. Because genomic localization of TFs is in part driven by TF recognition of DNA sequence, variation in TF binding sites can disrupt TF-DNA associations and affect gene regulation. To identify variants that impact TF binding in human brain tissues, we quantified allele bias for 93 TFs analyzed with ChIP-seq experiments of multiple structural brain regions from two donors. Using graph genomes constructed from phased genomic sequence data, we compared ChIP-seq signal between alleles at heterozygous variants within each tissue sample from each donor. Comparison of results from different brain regions within donors and the same regions between donors provided measures of allele bias reproducibility. We identified thousands of DNA variants that show reproducible bias in ChIP-seq for at least one TF. We found that alleles that are rarer in the general population were more likely than common alleles to exhibit large biases, and more frequently led to reduced TF binding. Combining ChIP-seq with RNA-seq, we identified TF-allele interaction biases with RNA bias in a phased allele linked to 6,709 eQTL variants identified in GTEx data, 3,309 of which were found in neural contexts. Our results provide insights into the effects of both common and rare variation on gene regulation in the brain. These findings can facilitate mechanistic understanding of cis-regulatory variation associated with biological traits, including disease.

Early-onset Alzheimer's disease explained by polygenic risk of late-onset disease?

Early-onset Alzheimer's disease (AD) is highly heritable, yet only 10% of cases are associated with known pathogenic mutations. For early-onset AD patients without an identified autosomal dominant cause, we hypothesized that their early-onset disease reflects further enrichment of the common risk-conferring single nucleotide polymorphisms associated with late-onset AD. We applied a previously validated polygenic hazard score for late-onset AD to 193 consecutive patients diagnosed at our tertiary dementia referral center with symptomatic early-onset AD. For comparison, we included 179 participants with late-onset AD and 70 healthy controls. Polygenic hazard scores were similar in early- versus late-onset AD. The polygenic hazard score was not associated with age-of-onset or disease biomarkers within early-onset AD. Early-onset AD does not represent an extreme enrichment of the common single nucleotide polymorphisms associated with late-onset AD. Further exploration of novel genetic risk factors of this highly heritable disease is warranted.Highlights: There is a unique genetic architecture of early- versus late-onset Alzheimer's disease (AD).Late-onset AD polygenic risk is not an explanation for early-onset AD.Polygenic risk of late-onset AD does not predict early-onset AD biology.Unique genetic architecture of early- versus late-onset AD parallels AD heterogeneity.

MAPT expression is mediated by long-range interactions with cis-regulatory elements.

Background: Tauopathies are a group of neurodegenerative diseases driven by abnormal aggregates of tau, a microtubule associated protein encoded by the MAPT gene. MAPT expression is absent in neural progenitor cells (NPCs) and increases during differentiation. This temporally dynamic expression pattern suggests that MAPT expression is controlled by transcription factors and cis-regulatory elements specific to differentiated cell types. Given the relevance of MAPT expression to neurodegeneration pathogenesis, identification of such elements is relevant to understanding genetic risk factors. Methods: We performed HiC, chromatin conformation capture (Capture-C), single-nucleus multiomics (RNA-seq+ATAC-seq), bulk ATAC-seq, and ChIP-seq for H3K27Ac and CTCF in NPCs and neurons differentiated from human iPSC cultures. We nominated candidate cis-regulatory elements (cCREs) for MAPT in human NPCs, differentiated neurons, and pure cultures of inhibitory and excitatory neurons. We then assayed these cCREs using luciferase assays and CRISPR interference (CRISPRi) experiments to measure their effects on MAPT expression. Finally, we integrated cCRE annotations into an analysis of genetic variation in AD cases and controls. Results: Using orthogonal genomics approaches, we nominated 94 cCREs for MAPT, including the identification of cCREs specifically active in differentiated neurons. Eleven regions enhanced reporter gene transcription in luciferase assays. Using CRISPRi, 5 of the 94 regions tested were identified as necessary for MAPT expression as measured by RT-qPCR and RNA-seq. Rare and predicted damaging genetic variation in both nominated and confirmed CREs was depleted in AD cases relative to controls (OR = 0.40, p = 0.004), consistent with the hypothesis that variants that disrupt MAPT enhancer activity, and thereby reduce MAPT expression, may be protective against neurodegenerative disease. Conclusions: We identified both proximal and distal regulatory elements for MAPT and confirmed the regulatory function for several regions, including three regions centromeric to MAPT beyond the well-described H1/H2 haplotype inversion breakpoint. This study provides compelling evidence for pursuing detailed knowledge of CREs for genes of interest to permit better understanding of disease risk.

Single nucleus multiomics identifies ZEB1 and MAFB as candidate regulators of Alzheimer's disease-specific cis-regulatory elements.

Cell type-specific transcriptional differences between brain tissues from donors with Alzheimer's disease (AD) and unaffected controls have been well documented, but few studies have rigorously interrogated the regulatory mechanisms responsible for these alterations. We performed single nucleus multiomics (snRNA-seq plus snATAC-seq) on 105,332 nuclei isolated from cortical tissues from 7 AD and 8 unaffected donors to identify candidate cis-regulatory elements (CREs) involved in AD-associated transcriptional changes. We detected 319,861 significant correlations, or links, between gene expression and cell type-specific transposase accessible regions enriched for active CREs. Among these, 40,831 were unique to AD tissues. Validation experiments confirmed the activity of many regions, including several candidate regulators of APP expression. We identified ZEB1 and MAFB as candidate transcription factors playing important roles in AD-specific gene regulation in neurons and microglia, respectively. Microglia links were globally enriched for heritability of AD risk and previously identified active regulatory regions.

Genetic associations with age at dementia onset in the PSEN1 E280A Colombian kindred.

INTRODUCTION: Genetic associations with Alzheimer's disease (AD) age at onset (AAO) could reveal genetic variants with therapeutic applications. We present a large Colombian kindred with autosomal dominant AD (ADAD) as a unique opportunity to discover AAO genetic associations. METHODS: A genetic association study was conducted to examine ADAD AAO in 340 individuals with the PSEN1 E280A mutation via TOPMed array imputation. Replication was assessed in two ADAD cohorts, one sporadic early-onset AD study and four late-onset AD studies. RESULTS: 13 variants had p<1×10-7 or p<1×10-5 with replication including three independent loci with candidate associations with clusterin including near CLU. Other suggestive associations were identified in or near HS3ST1, HSPG2, ACE, LRP1B, TSPAN10, and TSPAN14. DISCUSSION: Variants with suggestive associations with AAO were associated with biological processes including clusterin, heparin sulfate, and amyloid processing. The detection of these effects in the presence of a strong mutation for ADAD reinforces their potentially impactful role.

Environmental stress maintains trioecy in nematode worms.

Sex is determined by chromosomes in mammals but it can be influenced by the environment in many worms, crustaceans, and vertebrates. Despite this, there is little understanding of the relationship between ecology and the evolution of sexual systems. The nematode Auanema freiburgensis has a unique sex determination system in which individuals carrying one X chromosome develop into males while XX individuals develop into females in stress-free environments and self-fertile hermaphrodites in stressful environments. Theory predicts that trioecious populations with coexisting males, females, and hermaphrodites should be unstable intermediates in evolutionary transitions between mating systems. In this article, we study a mathematical model of reproductive evolution based on the unique life history and sex determination of A. freiburgensis. We develop the model in two scenarios, one where the relative production of hermaphrodites and females is entirely dependent on the environment and one based on empirical measurements of a population that displays incomplete, "leaky" environmental dependence. In the first scenario environmental conditions can push the population along an evolutionary continuum and result in the stable maintenance of multiple reproductive systems. The second "leaky" scenario results in the maintenance of three sexes for all environmental conditions. Theoretical investigations of reproductive system transitions have focused on the evolutionary costs and benefits of sex. Here, we show that the flexible sex determination system of A. freiburgensis may contribute to population-level resilience in the microscopic nematode's patchy, ephemeral natural habitat. Our results demonstrate that life history, ecology, and environment may play defining roles in the evolution of sexual systems.