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.

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.

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.

Deleterious, protein-altering variants in the transcriptional coregulator ZMYM3 in 27 individuals with a neurodevelopmental delay phenotype.

Neurodevelopmental disorders (NDDs) result from highly penetrant variation in hundreds of different genes, some of which have not yet been identified. Using the MatchMaker Exchange, we assembled a cohort of 27 individuals with rare, protein-altering variation in the transcriptional coregulator ZMYM3, located on the X chromosome. Most (n = 24) individuals were males, 17 of which have a maternally inherited variant; six individuals (4 male, 2 female) harbor de novo variants. Overlapping features included developmental delay, intellectual disability, behavioral abnormalities, and a specific facial gestalt in a subset of males. Variants in almost all individuals (n = 26) are missense, including six that recurrently affect two residues. Four unrelated probands were identified with inherited variation affecting Arg441, a site at which variation has been previously seen in NDD-affected siblings, and two individuals have de novo variation resulting in p.Arg1294Cys (c.3880C>T). All variants affect evolutionarily conserved sites, and most are predicted to damage protein structure or function. ZMYM3 is relatively intolerant to variation in the general population, is widely expressed across human tissues, and encodes a component of the KDM1A-RCOR1 chromatin-modifying complex. ChIP-seq experiments on one variant, p.Arg1274Trp, indicate dramatically reduced genomic occupancy, supporting a hypomorphic effect. While we are unable to perform statistical evaluations to definitively support a causative role for variation in ZMYM3, the totality of the evidence, including 27 affected individuals, recurrent variation at two codons, overlapping phenotypic features, protein-modeling data, evolutionary constraint, and experimentally confirmed functional effects strongly support ZMYM3 as an NDD-associated gene.

Phosphoinositide species and filamentous actin formation mediate engulfment by senescent tumor cells.

Cancer cells survive chemotherapy and cause lethal relapse by entering a senescent state that facilitates expression of many phagocytosis/macrophage-related genes that engender a novel cannibalism phenotype. We used biosensors and live-cell imaging to reveal the basic steps and mechanisms of engulfment by senescent human and mouse tumor cells. We show filamentous actin in predator cells was localized to the prey cell throughout the process of engulfment. Biosensors to various phosphoinositide (PI) species revealed increased concentration and distinct localization of predator PI(4) P and PI(4,5)P2 at the prey cell during early stages of engulfment, followed by a transient burst of PI(3) P before and following internalization. PIK3C2B, the kinase responsible for generating PI(3)P, was required for complete engulfment. Inhibition or knockdown of Clathrin, known to associate with PIK3C2B and PI(4,5)P2, severely impaired engulfment. In sum, our data reveal the most fundamental cellular processes of senescent cell engulfment, including the precise localizations and dynamics of actin and PI species throughout the entire process.

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.

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.

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.

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.

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.

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.

Navigating Hidden Hunger: An Exploratory Analysis of the Lived Experience of Food Insecurity among College Students.

College students are a vulnerable population to food insecurity (FI), which has significant implications for academic and health outcomes. The aims of this study were to explore the meaning of FI and its impact on students' lived experiences and food decisions, facilitators and barriers to food access as a student, and students' proposed solutions to address FI. Semi-structured, qualitative interviews were conducted with thirty students from a large, public land grant university in the Southeast United States. Grounded theory methodology was utilized with a constant comparative coding strategy to guide thematic analysis. Nine main themes emerged. Themes included the perceived meaning of FI, students' lived experience with FI, and food related coping strategies and decisions. Facilitators to food access were found to be social-networks and on-campus resources, while barriers to food access included financial burden of higher education, and stigma and social comparison. Proposed solutions to FI aligned with two main themes: food access solutions and information access solutions. Both of these themes included multiple subthemes that provided specific suggestions to address food insecurity for students. The findings aid in understanding the complex lived experience of FI and can inform future efforts to center student experiences, perceptions, and feedback into institutional frameworks to best meet student needs.

An Adenovirus early region 4 deletion mutant induces G2/M arrest via ATM activation and reduces expression of the mitotic marker phosphorylated (ser10) histone 3.

Adenovirus (Ad) type 5 (Ad5) early region 4 (E4) proteins inhibit the DNA damage response (DDR) including activation of the DDR kinase ATM and its substrates, which can induce G2/M cell cycle arrest. Infection with Ad5 or the E4 deletion mutant H5dl1007 (1007) resulted in the accumulation of post G1 cells with > 2 N cellular DNA content. A greater fraction of cells with 4 N DNA content was observed in 1007 infections compared to Ad5; this population was dependent on activation of ATM. G2/M checkpoint kinases, phosphorylated Chk2 (pChk2), and phosphorylated Cdk1 (pCdk1) were upregulated in 1007 infections, and 1007 showed reduced levels of the mitosis marker phosphorylated (Ser10) histone 3 compared to Ad5. Our results show that E4 mutant activation of ATM induces G2/M arrest via activation of checkpoint kinases, thereby contributing to viral-mediated regulation of the cell cycle.

Breast cancer cells survive chemotherapy by activating targetable immune-modulatory programs characterized by PD-L1 or CD80.

Breast cancer cells must avoid intrinsic and extrinsic cell death to relapse following chemotherapy. Entering senescence enables survival from mitotic catastrophe, apoptosis and nutrient deprivation, but mechanisms of immune evasion are poorly understood. Here we show that breast tumors surviving chemotherapy activate complex programs of immune modulation. Characterization of residual disease revealed distinct tumor cell populations. The first population was characterized by interferon response genes, typified by Cd274, whose expression required chemotherapy to enhance chromatin accessibility, enabling recruitment of IRF1 transcription factor. A second population was characterized by p53 signaling, typified by CD80 expression. Treating mammary tumors with chemotherapy followed by targeting the PD-L1 and/or CD80 axes resulted in marked accumulation of T cells and improved response; however, even combination strategies failed to fully eradicate tumors in the majority of cases. Our findings reveal the challenge of eliminating residual disease populated by senescent cells expressing redundant immune inhibitory pathways and highlight the need for rational immune targeting strategies.

An Analysis of Factors That Influence Students to Pursue Immunology.

One considers many factors before choosing a career path, such as interest, accessibility of resources, academic ability, and social network support. As employment around the world in science, technology, engineering, and math (STEM) disciplines continues to increase, there is a need to understand why students select specific majors in an effort to increase overall enrollment and retention of STEM majors. The purpose of this study was to elucidate how undergraduate and graduate students were introduced to immunology, a STEM discipline, and how these experiences influenced their desire to pursue immunology as a major. The findings from this study show that a majority of both immunology and nonimmunology majors were initially exposed to immunology through an educational experience compared with a personal experience. Our data also indicate that the timing of the experience is critical, such that an educational experience at an advanced academic level, for example, in college, or a personal experience as a teen or young adult correlated with the decision to pursue an immunology degree. Moreover, graduate students studying immunology report that having research experiences and/or an experience with a mentor positively influenced their decision to pursue immunology. Overall, the findings from this research highlight the type and timing of exposures that influence individuals to major in the field of immunology, and these data can be used in the future to increase the number of immunology graduates.

Using real-world examples of the COVID-19 pandemic to increase student confidence in their scientific literacy skills.

Over the last few decades, there has been a shift in the classroom from lecture-based to active learning settings with the argument that students retain more information when they are involved in the learning process. This correlation is even stronger when the active learning setting incorporates a real-world or personal connection. Using active learning activities that develop students' ability to comprehend primary scientific literature is particularly important in the field of immunology, due to the rapid expansion of information in the field, which has been further accelerated due to the COVID-19 pandemic. By nature, immunology is interdisciplinary, requiring an integrated knowledge of concepts from several scientific disciplines to understand complex immune processes. Engaging undergraduate students through the use of primary literature can improve scientific literacy, develop critical thinking, and enhance understanding of complex topics. To explore this, we utilized a group learning activity in an introductory immunology course that incorporated both a coronavirus-related review and COVID-19 clinical research article. We found that this learning activity significantly enhanced student confidence in key scientific literacy skills: reading scientific literature, clearly explaining relevant points, and describing conclusions generated from the data. Moreover, all students reported that they enjoyed the activity and that it helped them understand more about the current COVID-19 pandemic in the context of the immune response.

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.

BH3 mimetics selectively eliminate chemotherapy-induced senescent cells and improve response in TP53 wild-type breast cancer.

TP53 wild-type breast tumors rarely undergo a complete pathological response after chemotherapy treatment. These patients have an extremely poor survival rate and studies show these tumors preferentially undergo senescence instead of apoptosis. These senescent cells persist after chemotherapy and secrete cytokines and chemokines comprising the senescence associated secretory phenotype, which promotes survival, proliferation, and metastasis. We hypothesized that eliminating senescent tumor cells would improve chemotherapy response and extend survival. Previous studies have shown "senolytic" agents selectively kill senescent normal cells, but their efficacy in killing chemotherapy-induced senescent cancer cells is unknown. We show that ABT-263, a BH3 mimetic that targets antiapoptotic proteins BCL2/BCL-XL/BCL-W, had no effect on proliferating cells, but rapidly and selectively induced apoptosis in a subset of chemotherapy-treated cancer cells, though sensitivity required days to develop. Low NOXA expression conferred resistance to ABT-263 in some cells, necessitating additional MCL1 inhibition. Gene editing confirmed breast cancer cells relied on BCL-XL or BCL-XL/MCL1 for survival in senescence. In a mouse model of breast cancer, ABT-263 treatment following chemotherapy led to apoptosis, greater tumor regression, and longer survival. Our results reveal cancer cells that have survived chemotherapy by entering senescence can be eliminated using BH3 mimetic drugs that target BCL-XL or BCL-XL/MCL1. These drugs could help minimize residual disease and extend survival in breast cancer patients that otherwise have a poor prognosis and are most in need of improved therapies.

IL-7 signalling represses Bcl-6 and the TFH gene program.

The transcriptional repressor Bcl-6 is linked to the development of both CD4(+) T follicular helper (TFH) and central memory T (TCM) cells. Here, we demonstrate that in response to decreased IL-2 signalling, T helper 1 (TH1) cells upregulate Bcl-6 and co-initiate TFH- and TCM-like gene programs, including expression of the cytokine receptors IL-6Rα and IL-7R. Exposure of this potentially bi-potent cell population to IL-6 favours the TFH gene program, whereas IL-7 signalling represses TFH-associated genes including Bcl6 and Cxcr5, but not the TCM-related genes Klf2 and Sell. Mechanistically, IL-7-dependent activation of STAT5 contributes to Bcl-6 repression. Importantly, antigen-specific IL-6Rα(+)IL-7R(+) CD4(+) T cells emerge from the effector population at late time points post influenza infection. These data support a novel role for IL-7 in the repression of the TFH gene program and evoke a divergent regulatory mechanism by which post-effector TH1 cells may contribute to long-term cell-mediated and humoral immunity.