Search

Long COVID (LC) occurs after at least 10% of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections, yet its etiology remains poorly understood. We used omic assays and serology to deeply characterize the global and SARS-CoV-2-specific immunity in the blood of individuals with clear LC and non-LC clinical trajectories, 8 months postinfection. We found that LC individuals exhibited systemic inflammation and immune dysregulation. This was evidenced by global differences in T cell subset distribution implying ongoing immune responses, as well as by sex-specific perturbations in cytolytic subsets. LC individuals displayed increased frequencies of CD4+ T cells poised to migrate to inflamed tissues and exhausted SARS-CoV-2-specific CD8+ T cells, higher levels of SARS-CoV-2 antibodies and a mis-coordination between their SARS-CoV-2-specific T and B cell responses. Our analysis suggested an improper crosstalk between the cellular and humoral adaptive immunity in LC, which can lead to immune dysregulation, inflammation and clinical symptoms associated with this debilitating condition.

Apolipoprotein E4 (APOE4) is the strongest genetic risk factor for late-onset Alzheimers disease (LOAD), leading to earlier age of clinical onset and exacerbating pathologies. There is a critical need to identify protective targets. Recently, a rare APOE variant, APOE3-R136S (Christchurch), was found to protect against early-onset AD in a PSEN1-E280A carrier. In this study, we sought to determine if the R136S mutation also protects against APOE4-driven effects in LOAD. We generated tauopathy mouse and human iPSC-derived neuron models carrying human APOE4 with the homozygous or heterozygous R136S mutation. We found that the homozygous R136S mutation rescued APOE4-driven Tau pathology, neurodegeneration and neuroinflammation. The heterozygous R136S mutation partially protected against APOE4-driven neurodegeneration and neuroinflammation but not Tau pathology. Single-nucleus RNA sequencing revealed that the APOE4-R136S mutation increased disease-protective and diminished disease-associated cell populations in a gene dose-dependent manner. Thus, the APOE-R136S mutation protects against APOE4-driven AD pathologies, providing a target for therapeutic development against AD.

Apolipoprotein E4 (APOE4) is the strongest known genetic risk factor for late-onset Alzheimer's disease (AD). Conditions of stress or injury induce APOE expression within neurons, but the role of neuronal APOE4 in AD pathogenesis is still unclear. Here we report the characterization of neuronal APOE4 effects on AD-related pathologies in an APOE4-expressing tauopathy mouse model. The selective genetic removal of APOE4 from neurons led to a significant reduction in tau pathology, gliosis, neurodegeneration, neuronal hyperexcitability and myelin deficits. Single-nucleus RNA-sequencing revealed that the removal of neuronal APOE4 greatly diminished neurodegenerative disease-associated subpopulations of neurons, oligodendrocytes, astrocytes and microglia whose accumulation correlated to the severity of tau pathology, neurodegeneration and myelin deficits. Thus, neuronal APOE4 plays a central role in promoting the development of major AD pathologies and its removal can mitigate the progressive cellular and tissue alterations occurring in this model of APOE4-driven tauopathy.

BackgroundEndometriosis is a chronic, estrogen-dependent disorder where inflammation contributes to disease-associated symptoms of pelvic pain and infertility. Immune dysfunction includes insufficient immune lesion clearance, a pro-inflammatory endometrial environment, and systemic inflammation. Comprehensive understanding of endometriosis immune pathophysiology in different hormonal milieu and disease severity has been hampered by limited direct characterization of immune populations in endometrium, blood, and lesions. Simultaneous deep phenotyping at single-cell resolution of complex tissues has transformed our understanding of the immune system and its role in many diseases. Herein, we report mass cytometry and high dimensional analyses to study immune cell phenotypes, abundance, activation states, and functions in endometrium and blood of women with and without endometriosis in different cycle phases and disease stages.MethodsA case-control study was designed. Endometrial biopsies and blood (n = 60 total) were obtained from women with (n = 20, n = 17, respectively) and without (n = 14, n = 9) endometriosis in the proliferative and secretory cycle phases of the menstrual cycle. Two mass cytometry panels were designed: one broad panel and one specific for mononuclear phagocytic cells (MPC), and all samples were multiplexed to characterize both endometrium and blood immune composition at unprecedented resolution. We combined supervised and unsupervised analyses to finely define the immune cell subsets with an emphasis on MPC. Then, association between cell types, protein expression, disease status, and cycle phase were performed.ResultsThe broad panel highlighted a significant modification of MPC in endometriosis; thus, they were studied in detail with an MPC-focused panel. Endometrial CD91+ macrophages overexpressed SIRPα (phagocytosis inhibitor) and CD64 (associated with inflammation) in endometriosis, and they were more abundant in mild versus severe disease. In blood, classical and intermediate monocytes were less abundant in endometriosis, whereas plasmacytoid dendritic cells and non-classical monocytes were more abundant. Non-classical monocytes were higher in severe versus mild disease.ConclusionsA greater inflammatory phenotype and decreased phagocytic capacity of endometrial macrophages in endometriosis are consistent with defective clearance of endometrial cells shed during menses and in tissue homeostasis, with implications in endometriosis pathogenesis and pathophysiology. Different proportions of monocytes and plasmacytoid dendritic cells in blood from endometriosis suggest systemically aberrant functionality of the myeloid system opening new venues for the study of biomarkers and therapies for endometriosis.

BACKGROUND: GATA4 (GATA-binding protein 4), a zinc finger-containing, DNA-binding transcription factor, is essential for normal cardiac development and homeostasis in mice and humans, and mutations in this gene have been reported in human heart defects. Defects in alternative splicing are associated with many heart diseases, yet relatively little is known about how cell type- or cell state-specific alternative splicing is achieved in the heart. Here, we show that GATA4 regulates cell type-specific splicing through direct interaction with RNA and the spliceosome in human induced pluripotent stem cell-derived cardiac progenitors. METHODS: We leveraged a combination of unbiased approaches including affinity purification of GATA4 and mass spectrometry, enhanced cross-linking with immunoprecipitation, electrophoretic mobility shift assays, in vitro splicing assays, and unbiased transcriptomic analysis to uncover GATA4s novel function as a splicing regulator in human induced pluripotent stem cell-derived cardiac progenitors. RESULTS: We found that GATA4 interacts with many members of the spliceosome complex in human induced pluripotent stem cell-derived cardiac progenitors. Enhanced cross-linking with immunoprecipitation demonstrated that GATA4 also directly binds to a large number of mRNAs through defined RNA motifs in a sequence-specific manner. In vitro splicing assays indicated that GATA4 regulates alternative splicing through direct RNA binding, resulting in functionally distinct protein products. Correspondingly, knockdown of GATA4 in human induced pluripotent stem cell-derived cardiac progenitors resulted in differential alternative splicing of genes involved in cytoskeleton organization and calcium ion import, with functional consequences associated with the protein isoforms. CONCLUSIONS: This study shows that in addition to its well described transcriptional function, GATA4 interacts with members of the spliceosome complex and regulates cell type-specific alternative splicing via sequence-specific interactions with RNA. Several genes that have splicing regulated by GATA4 have functional consequences and many are associated with dilated cardiomyopathy, suggesting a novel role for GATA4 in achieving the necessary cardiac proteome in normal and stress-responsive conditions.

Endothelial and erythropoietic lineages arise from a common developmental progenitor. Etv2 is a master transcriptional regulator required for the development of both lineages. However, the mechanisms through which Etv2 initiates the gene-regulatory networks (GRNs) for endothelial and erythropoietic specification and how the two GRNs diverge downstream of Etv2 remain incompletely understood. Here, by analyzing a hypomorphic Etv2 mutant, we demonstrate different threshold requirements for initiation of the downstream GRNs for endothelial and erythropoietic development. We show that Etv2 functions directly in a coherent feedforward transcriptional network for vascular endothelial development, and a low level of Etv2 expression is sufficient to induce and sustain the endothelial GRN. In contrast, Etv2 induces the erythropoietic GRN indirectly via activation of Tal1, which requires a significantly higher threshold of Etv2 to initiate and sustain erythropoietic development. These results provide important mechanistic insight into the divergence of the endothelial and erythropoietic lineages.

BackgroundThe regenerative capacity of the heart after myocardial infarction is limited. Our previous study showed that ectopic introduction of 4 cell cycle factors (4F; CDK1 [cyclin-dependent kinase 1], CDK4 [cyclin-dependent kinase 4], CCNB [cyclin B1], and CCND [cyclin D1]) promotes cardiomyocyte proliferation in 15% to 20% of infected cardiomyocytes in vitro and in vivo and improves cardiac function after myocardial infarction in mice.MethodsUsing temporal single-cell RNA sequencing, we aimed to identify the necessary reprogramming stages during the forced cardiomyocyte proliferation with 4F on a single cell basis. Using rat and pig models of ischemic heart failure, we aimed to start the first preclinical testing to introduce 4F gene therapy as a candidate for the treatment of ischemia-induced heart failure.ResultsTemporal bulk and single-cell RNA sequencing and further biochemical validations of mature human induced pluripotent stem cell-derived cardiomyocytes treated with either LacZ or 4F adenoviruses revealed full cell cycle reprogramming in 15% of the cardiomyocyte population at 48 hours after infection with 4F, which was associated mainly with sarcomere disassembly and metabolic reprogramming (n=3/time point/group). Transient overexpression of 4F, specifically in cardiomyocytes, was achieved using a polycistronic nonintegrating lentivirus (NIL) encoding 4F; each is driven by a TNNT2 (cardiac troponin T isoform 2) promoter (TNNT2-4Fpolycistronic-NIL). TNNT2-4Fpolycistronic-NIL or control virus was injected intramyocardially 1 week after myocardial infarction in rats (n=10/group) or pigs (n=6-7/group). Four weeks after injection, TNNT2-4Fpolycistronic-NIL-treated animals showed significant improvement in left ventricular ejection fraction and scar size compared with the control virus-treated animals. At 4 months after treatment, rats that received TNNT2-4Fpolycistronic-NIL still showed a sustained improvement in cardiac function and no obvious development of cardiac arrhythmias or systemic tumorigenesis (n=10/group).ConclusionsThis study provides mechanistic insights into the process of forced cardiomyocyte proliferation and advances the clinical feasibility of this approach by minimizing the oncogenic potential of the cell cycle factors owing to the use of a novel transient and cardiomyocyte-specific viral construct.

Congenital heart disease (CHD) is present in 1% of live births, yet identification of causal mutations remains challenging. We hypothesized that genetic determinants for CHDs may lie in the protein interactomes of transcription factors whose mutations cause CHDs. Defining the interactomes of two transcription factors haplo-insufficient in CHD, GATA4 and TBX5, within human cardiac progenitors, and integrating the results with nearly 9,000 exomes from proband-parent trios revealed an enrichment of de novo missense variants associated with CHD within the interactomes. Scoring variants of interactome members based on residue, gene, and proband features identified likely CHD-causing genes, including the epigenetic reader GLYR1. GLYR1 and GATA4 widely co-occupied and co-activated cardiac developmental genes, and the identified GLYR1 missense variant disrupted interaction with GATA4, impairing in vitro and in vivo function in mice. This integrative proteomic and genetic approach provides a framework for prioritizing and interrogating genetic variants in heart disease.

Differentiation proceeds along a continuum of increasingly fate-restricted intermediates, referred to as canalization1,2. Canalization is essential for stabilizing cell fate, but the mechanisms that underlie robust canalization are unclear. Here we show that the BRG1/BRM-associated factor (BAF) chromatin-remodelling complex ATPase gene Brm safeguards cell identity during directed cardiogenesis of mouse embryonic stem cells. Despite the establishment of a well-differentiated precardiac mesoderm, Brm-/- cells predominantly became neural precursors, violating germ layer assignment. Trajectory inference showed a sudden acquisition of a non-mesodermal identity in Brm-/- cells. Mechanistically, the loss of Brm prevented de novo accessibility of primed cardiac enhancers while increasing the expression of neurogenic factor POU3F1, preventing the binding of the neural suppressor REST and shifting the composition of BRG1 complexes. The identity switch caused by the Brm mutation was overcome by increasing BMP4 levels during mesoderm induction. Mathematical modelling supports these observations and demonstrates that Brm deletion affects cell fate trajectory by modifying saddle-node bifurcations2. In the mouse embryo, Brm deletion exacerbated mesoderm-deleted Brg1-mutant phenotypes, severely compromising cardiogenesis, and reveals an in vivo role for Brm. Our results show that Brm is a compensable safeguard of the fidelity of mesoderm chromatin states, and support a model in which developmental canalization is not a rigid irreversible path, but a highly plastic trajectory.

Answer ALS is a biological and clinical resource of patient-derived, induced pluripotent stem (iPS) cell lines, multi-omic data derived from iPS neurons and longitudinal clinical and smartphone data from over 1,000 patients with ALS. This resource provides population-level biological and clinical data that may be employed to identify clinical-molecular-biochemical subtypes of amyotrophic lateral sclerosis (ALS). A unique smartphone-based system was employed to collect deep clinical data, including fine motor activity, speech, breathing and linguistics/cognition. The iPS spinal neurons were blood derived from each patient and these cells underwent multi-omic analytics including whole-genome sequencing, RNA transcriptomics, ATAC-sequencing and proteomics. The intent of these data is for the generation of integrated clinical and biological signatures using bioinformatics, statistics and computational biology to establish patterns that may lead to a better understanding of the underlying mechanisms of disease, including subgroup identification. A web portal for open-source sharing of all data was developed for widespread community-based data analytics.

T and natural killer (NK) cells are effector cells with key roles in anti-HIV immunity, including in lymphoid tissues, the major site of HIV persistence. However, little is known about the features of these effector cells from people living with HIV (PLWH), particularly from those who initiated antiretroviral therapy (ART) during acute infection. Our study design was to use 42-parameter CyTOF to conduct deep phenotyping of paired blood- and lymph node (LN)-derived T and NK cells from three groups of HIV+ aviremic individuals: elite controllers (N = 5), and ART-suppressed individuals who had started therapy during chronic (N = 6) vs. acute infection (N = 8), the latter of which is associated with better outcomes. We found that acute-treated individuals are enriched for specific subsets of T and NK cells, including blood-derived CD56-CD16+ NK cells previously associated with HIV control, and LN-derived CD4+ T follicular helper cells with heightened expansion potential. An in-depth comparison of the features of the cells from blood vs. LNs of individuals from our cohort revealed that T cells from blood were more activated than those from LNs. By contrast, LNs were enriched for follicle-homing CXCR5+ CD8+ T cells, which expressed increased levels of inhibitory receptors and markers of survival and proliferation as compared to their CXCR5- counterparts. In addition, a subset of memory-like CD56brightTCF1+ NK cells was enriched in LNs relative to blood. These results together suggest unique T and NK cell features in acute-treated individuals, and highlight the importance of examining effector cells not only in blood but also the lymphoid tissue compartment, where the reservoir mostly persists, and where these cells take on distinct phenotypic features.

Congenital heart disease (CHD) is present in 1% of live births, yet despite large-scale genomic sequencing efforts, identification of causal mutations remains a challenge. We hypothesized that genetic determinants for CHDs may lie in the protein interactomes of GATA4 and TBX5, two transcription factors whose mutation cause CHDs. Defining the GATA4 or TBX5 interactomes in human cardiac progenitors via affinity purification-mass spectrometry and integrating the results with genetic data from the Pediatric Cardiac Genomic Consortium revealed an enrichment of de novo variants associated with CHD. A consolidative score that prioritized interactome members based on variant, gene, and proband features identified likely CHD-causing genes, including the epigenetic reader GLYR1. GLYR1 and GATA4 widely co-occupied and co-activated cardiac developmental genes, and the GLYR1 missense variant identified disrupted interaction with GATA4 and impaired transcriptional co-regulation in cardiomyocyte differentiation in vitro and cardiogenesis in vivo. This integrative proteomic and genetic approach provides a framework for prioritizing and interrogating the contribution of genetic variants in disease.

Extrinsic inhibitors at sites of blood-brain barrier disruption and neurovascular damage contribute to remyelination failure in neurological diseases. However, therapies to overcome the extrinsic inhibition of remyelination are not widely available and the dynamics of glial progenitor niche remodelling at sites of neurovascular dysfunction are largely unknown. By integrating in vivo two-photon imaging co-registered with electron microscopy and transcriptomics in chronic neuroinflammatory lesions, we found that oligodendrocyte precursor cells clustered perivascularly at sites of limited remyelination with deposition of fibrinogen, a blood coagulation factor abundantly deposited in multiple sclerosis lesions. By developing a screen (OPC-X-screen) to identify compounds that promote remyelination in the presence of extrinsic inhibitors, we showed that known promyelinating drugs did not rescue the extrinsic inhibition of remyelination by fibrinogen. In contrast, bone morphogenetic protein type I receptor blockade rescued the inhibitory fibrinogen effects and restored a promyelinating progenitor niche by promoting myelinating oligodendrocytes, while suppressing astrocyte cell fate, with potent therapeutic effects in chronic models of multiple sclerosis. Thus, abortive oligodendrocyte precursor cell differentiation by fibrinogen is refractory to known promyelinating compounds, suggesting that blockade of the bone morphogenetic protein signalling pathway may enhance remyelinating efficacy by overcoming extrinsic inhibition in neuroinflammatory lesions with vascular damage.

Diverse gene products contribute to the pathogenesis of Alzheimer's disease (AD). Experimental models have helped elucidate their mechanisms and impact on brain functions. Human amyloid precursor protein (hAPP) transgenic mice from line J20 (hAPP-J20 mice) are widely used to simulate key aspects of AD. However, they also carry an insertional mutation in noncoding sequence of one Zbtb20 allele, a gene involved in neural development. We demonstrate that heterozygous hAPP-J20 mice have reduced Zbtb20 expression in some AD-relevant brain regions, but not others, and that Zbtb20 levels are higher in hAPP-J20 mice than heterozygous Zbtb20 knock-out (Zbtb20 +/-) mice. Whereas hAPP-J20 mice have premature mortality, severe deficits in learning and memory, other behavioral alterations, and prominent nonconvulsive epileptiform activity, Zbtb20 +/- mice do not. Thus, the insertional mutation in hAPP-J20 mice does not ablate the affected Zbtb20 allele and is unlikely to account for the AD-like phenotype of this model.

Diverse gene products contribute to the pathogenesis of Alzheimer's disease (AD). Experimental models have helped elucidate their mechanisms and impact on brain functions. Human amyloid precursor protein (hAPP) transgenic mice from line J20 (hAPP-J20 mice) are widely used to simulate key aspects of AD. However, they also carry an insertional mutation in noncoding sequence of one Zbtb20 allele, a gene involved in neural development. We demonstrate that heterozygous hAPP-J20 mice have reduced Zbtb20 expression in some AD-relevant brain regions, but not others, and that Zbtb20 levels are higher in hAPP-J20 mice than heterozygous Zbtb20 knock-out (Zbtb20+/-) mice. Whereas hAPP-J20 mice have premature mortality, severe deficits in learning and memory, other behavioral alterations, and prominent nonconvulsive epileptiform activity, Zbtb20+/- mice do not. Thus, the insertional mutation in hAPP-J20 mice does not ablate the affected Zbtb20 allele and is unlikely to account for the AD-like phenotype of this model.

Haploinsufficiency of transcriptional regulators causes human congenital heart disease (CHD); however, the underlying CHD gene regulatory network (GRN) imbalances are unknown. Here, we define transcriptional consequences of reduced dosage of the CHD transcription factor, TBX5, in individual cells during cardiomyocyte differentiation from human induced pluripotent stem cells (iPSCs). We discovered highly sensitive dysregulation of TBX5-dependent pathways-including lineage decisions and genes associated with heart development, cardiomyocyte function, and CHD genetics-in discrete subpopulations of cardiomyocytes. Spatial transcriptomic mapping revealed chamber-restricted expression for many TBX5-sensitive transcripts. GRN analysis indicated that cardiac network stability, including vulnerable CHD-linked nodes, is sensitive to TBX5 dosage. A GRN-predicted genetic interaction between Tbx5 and Mef2c, manifesting as ventricular septation defects, was validated in mice. These results demonstrate exquisite and diverse sensitivity to TBX5 dosage in heterogeneous subsets of iPSC-derived cardiomyocytes and predicts candidate GRNs for human CHDs, with implications for quantitative transcriptional regulation in disease.

Microglial surveillance is a key feature of brain physiology and disease. Here, we found that Gi-dependent microglial dynamics prevent neuronal network hyperexcitability. By generating MgPTX mice to genetically inhibit Gi in microglia, we show that sustained reduction of microglia brain surveillance and directed process motility induced spontaneous seizures and increased hypersynchrony after physiologically evoked neuronal activity in awake adult mice. Thus, Gi-dependent microglia dynamics may prevent hyperexcitability in neurological diseases.

The latent reservoir is a major barrier to HIV cure. As latently infected cells cannot be phenotyped directly, the features of the in vivo reservoir have remained elusive. Here, we describe a method that leverages high-dimensional phenotyping using CyTOF to trace latently infected cells reactivated ex vivo to their original pre-activation states. Our results suggest that, contrary to common assumptions, the reservoir is not randomly distributed among cell subsets, and is remarkably conserved between individuals. However, reservoir composition differs between tissues and blood, as do cells successfully reactivated by different latency reversing agents. By selecting 8-10 of our 39 original CyTOF markers, we were able to isolate highly purified populations of unstimulated in vivo latent cells. These purified populations were highly enriched for replication-competent and intact provirus, transcribed HIV, and displayed clonal expansion. The ability to isolate unstimulated latent cells from infected individuals enables previously impossible studies on HIV persistence.

To discover regulatory elements driving the specificity of gene expression in different cell types and regions of the developing human brain, we generated an atlas of open chromatin from nine dissected regions of the mid-gestation human telencephalon, as well as microdissected upper and deep layers of the prefrontal cortex. We identified a subset of open chromatin regions (OCRs), termed predicted regulatory elements (pREs), that are likely to function as developmental brain enhancers. pREs showed temporal, regional, and laminar differences in chromatin accessibility and were correlated with gene expression differences across regions and gestational ages. We identified two functional de novo variants in a pRE for autism risk gene SLC6A1, and using CRISPRa, demonstrated that this pRE regulates SCL6A1. Additionally, mouse transgenic experiments validated enhancer activity for pREs proximal to FEZF2 and BCL11A. Thus, this atlas serves as a resource for decoding neurodevelopmental gene regulation in health and disease.

Oxidative stress is a central part of innate immune-induced neurodegeneration. However, the transcriptomic landscape of central nervous system (CNS) innate immune cells contributing to oxidative stress is unknown, and therapies to target their neurotoxic functions are not widely available. Here, we provide the oxidative stress innate immune cell atlas in neuroinflammatory disease and report the discovery of new druggable pathways. Transcriptional profiling of oxidative stress-producing CNS innate immune cells identified a core oxidative stress gene signature coupled to coagulation and glutathione-pathway genes shared between a microglia cluster and infiltrating macrophages. Tox-seq followed by a microglia high-throughput screen and oxidative stress gene network analysis identified the glutathione-regulating compound acivicin, with potent therapeutic effects that decrease oxidative stress and axonal damage in chronic and relapsing multiple sclerosis models. Thus, oxidative stress transcriptomics identified neurotoxic CNS innate immune populations and may enable discovery of selective neuroprotective strategies.

The contraceptive effectiveness of intrauterine devices (IUDs) has been attributed in part to a foreign body reaction in the endometrium. We performed this study to better understand mechanisms of action of contraceptives of by studying their effects on endometrial and cervical transcriptomes. We collected endometrial and cervical biopsies from women using the levonorgestrel-releasing intrauterine system (LNG-IUS, n = 11), copper intrauterine device (cu-IUD, n = 13) or levonorgestrel-containing combined oral contraceptives (COC, n = 12), and from women not using contraceptives (control group, n = 11). Transcriptional profiling was performed with Affymetrix arrays, Principal Component Analysis and the bioconductor package limma. In endometrial samples from cu-IUD users, there were no genes with statistically significant differential expression compared to controls. In LNG-IUS users, 2509 genes were differentially expressed and mapped predominantly onto immune and inflammatory pathways. The cervical samples showed no statistically significant differential gene expression compared to controls. Hormonal and copper IUDs have significantly different effects on the endometrial transcriptome, with the LNG-IUS transcriptome showing pronounced inflammation and immune activation compared to controls whereas the cu-IUD transcriptome was indistinguishable from luteal phase endometrium. These findings argue against a foreign body reaction as a common mechanism of action of IUDs.

Abstract Gene expression differs between cell types and regions within complex tissues such as the developing brain. To discover regulatory elements underlying this specificity, we generated genome-wide maps of chromatin accessibility in eleven anatomically-defined regions of the developing human telencephalon, including upper and deep layers of the prefrontal cortex. We predicted a subset of open chromatin regions (18%) that are most likely to be active enhancers, many of which are dynamic with 26% differing between early and late mid-gestation and 28% present in only one brain region. These region-specific predicted regulatory elements (pREs) are enriched proximal to genes with expression differences across regions and developmental stages and harbor distinct sequence motifs that suggest potential upstream regulators of regional and temporal transcription. We leverage this atlas to identify regulators of genes associated with autism spectrum disorder (ASD) including an enhancer of BCL11A , validated in mouse, and two functional de novo mutations in individuals with ASD in an enhancer of SLC6A1 , validated in neuroblastoma cells. These applications demonstrate the utility of this atlas for decoding neurodevelopmental gene regulation in health and disease. Summary To discover regulatory elements driving the specificity of gene expression in different cell types and regions of the developing human brain, we generated an atlas of open chromatin from eleven dissected regions of the mid-gestation human telencephalon, including upper and deep layers of the prefrontal cortex. We identified a subset of open chromatin regions (OCRs), termed predicted regulatory elements (pREs), that are likely to function as developmental brain enhancers. pREs showed regional differences in chromatin accessibility, including many specific to one brain region, and were correlated with gene expression differences across the same regions and gestational ages. pREs allowed us to map neurodevelopmental disorder risk genes to developing telencephalic regions, and we identified three functional de novo noncoding variants in pREs that alter enhancer function. In addition, transgenic experiments in mouse validated enhancer activity for a pRE proximal to BCL11A , showing how this atlas serves as a resource for decoding neurodevelopmental gene regulation in health and disease.

Chromatin remodeling complexes instruct cellular differentiation and lineage specific transcription. The BRG1/BRM-associated factor (BAF) complexes are important for several aspects of differentiation. We show that the catalytic subunit gene Brg1 has a specific role in cardiac precursors (CPs) to initiate cardiac gene expression programs and repress non-cardiac expression. Using immunopurification with mass spectrometry, we have determined the dynamic composition of BAF complexes during mammalian cardiac differentiation, identifying several cell-type specific subunits. We focused on the CP- and cardiomyocyte (CM)-enriched subunits BAF60c (SMARCD3) and BAF170 (SMARCC2). Baf60c and Baf170 co-regulate gene expression with Brg1 in CPs, and in CMs their loss results in broadly deregulated cardiac gene expression. BRG1, BAF60c and BAF170 modulate chromatin accessibility, to promote accessibility at activated genes while closing chromatin at repressed genes. BAF60c and BAF170 are required for proper BAF complex composition, and BAF170 loss leads to retention of BRG1 at CP-specific sites. Thus, dynamic interdependent BAF complex subunit assembly modulates chromatin states and thereby participates in directing temporal gene expression programs in cardiogenesis.

Ectopic expression of combinations of transcription factors (TFs) can drive direct lineage conversion, thereby reprogramming a somatic cell's identity. To determine the molecular mechanisms by which Gata4, Mef2c, and Tbx5 (GMT) induce conversion from a cardiac fibroblast toward an induced cardiomyocyte, we performed comprehensive transcriptomic, DNA-occupancy, and epigenomic interrogation throughout the reprogramming process. Integration of these datasets identified new TFs involved in cardiac reprogramming and revealed context-specific roles for GMT, including the ability of Mef2c and Tbx5 to independently promote chromatin remodeling at previously inaccessible sites. We also find evidence for cooperative facilitation and refinement of each TF's binding profile in a combinatorial setting. A reporter assay employing newly defined regulatory elements confirmed that binding of a single TF can be sufficient for gene activation, suggesting that co-binding events do not necessarily reflect synergy. These results shed light on fundamental mechanisms by which combinations of TFs direct lineage conversion.

BackgroundPreviously, using microarrays and mRNA-Sequencing (mRNA-Seq) we found that occupational exposure to a range of benzene levels perturbed gene expression in peripheral blood mononuclear cells.ObjectivesIn the current study, we sought to identify gene expression biomarkers predictive of benzene exposure below 1 part per million (ppm), the occupational standard in the U.S.MethodsFirst, we used the nCounter platform to validate altered expression of 30 genes in 33 unexposed controls and 57 subjects exposed to benzene (0.7, p0.9 (p

BackgroundHyperexcitability of neuronal networks can lead to excessive release of the excitatory neurotransmitter glutamate, which in turn can cause neuronal damage by overactivating NMDA-type glutamate receptors and related signaling pathways. This process (excitotoxicity) has been implicated in the pathogenesis of many neurological conditions, ranging from childhood epilepsies to stroke and neurodegenerative disorders such as Alzheimer's disease (AD). Reducing neuronal levels of the microtubule-associated protein tau counteracts network hyperexcitability of diverse causes, but whether this strategy can also diminish downstream excitotoxicity is less clear.MethodsWe established a cell-based assay to quantify excitotoxicity in primary cultures of mouse hippocampal neurons and investigated the role of tau in exicitotoxicity by modulating neuronal tau expression through genetic ablation or transduction with lentiviral vectors expressing anti-tau shRNA or constructs encoding wildtype versus mutant mouse tau.ResultsWe demonstrate that shRNA-mediated knockdown of tau reduces glutamate-induced, NMDA receptor-dependent Ca2+ influx and neurotoxicity in neurons from wildtype mice. Conversely, expression of wildtype mouse tau enhances Ca2+ influx and excitotoxicity in tau-deficient (Mapt -/-) neurons. Reconstituting tau expression in Mapt -/- neurons with mutant forms of tau reveals that the tau-related enhancement of Ca2+ influx and excitotoxicity depend on the phosphorylation of tau at tyrosine 18 (pY18), which is mediated by the tyrosine kinase Fyn. These effects are most evident at pathologically elevated concentrations of glutamate, do not involve GluN2B-containing NMDA receptors, and do not require binding of Fyn to tau's major interacting PxxP motif or of tau to microtubules.ConclusionsAlthough tau has been implicated in diverse neurological diseases, its most pathogenic forms remain to be defined. Our study suggests that reducing the formation or level of pY18-tau can counteract excitotoxicity by diminishing NMDA receptor-dependent Ca2+ influx.

AbstractChromatin remodeling complexes instruct cellular differentiation and lineage specific transcription. The BRG1/BRM associated factor (BAF) complexes are important for several aspects of differentiation. We show that the catalytic subunit Brg1 has a specific role in cardiac precursors (CPs) to initiate cardiac gene expression programs and repress non-cardiac expression. Using immunoprecipitation with mass spectrometry (IP-MS), we determined the dynamic composition of BAF complexes during mammalian cardiac differentiation, and identified BAF60c (SMARCD3) and BAF170 (SMARCC2) as subunits enriched in CPs and cardiomyocytes (CM). Baf60c and Baf170 co-regulate gene expression with Brg1 in CPs, but in CMs control different gene expression programs, although still promoting a cardiac-specific gene set. BRG1, BAF60, and BAF170 all modulate chromatin accessibility, to either promote accessibility at activated genes, while closing up chromatin at repressed genes. BAF60c and BAF170 are required for proper BAF complex composition and stoichiometry, and promote BRG1 occupancy in CM. Additionally, BAF170 facilitates expulsion of BRG1-containing complexes in the transition from CP to CM. Thus, dynamic interdependent BAF complex subunit assembly modulates chromatin states and thereby directs temporal gene expression programs in cardiogenesis.Significance statementBRG1/BRM associated factors (BAF) form multi-subunit protein complexes that reorganize chromatin and regulate transcription. Specific BAF complex subunits have important roles during cell differentiation and development. We systematically identify BAF subunit composition and find temporal enrichment of subunits during cardiomyocyte differentiation. We find the catalytic subunit BRG1 has important contributions in initiating gene expression programs in cardiac progenitors along with cardiac-enriched subunits BAF60c and BAF170. Both these proteins regulated BAF subunit composition and chromatin accessibility and prevent expression of non-cardiac developmental genes during precursor to cardiomyocyte differentiation. Mechanistically, we find BAF170 destabilizes the BRG1 complex and expels BRG1 from cardiomyocyte-specific genes. Thus, our data shows synergies between diverse BAF subunits in facilitating temporal gene expression programs during cardiogenesis.

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is an important tool for studying gene regulatory proteins, such as transcription factors and histones. Peak calling is one of the first steps in the analysis of these data. Peak calling consists of two sub-problems: identifying candidate peaks and testing candidate peaks for statistical significance. We surveyed 30 methods and identified 12 features of the two sub-problems that distinguish methods from each other. We picked six methods GEM, MACS2, MUSIC, BCP, Threshold-based method (TM) and ZINBA] that span this feature space and used a combination of 300 simulated ChIP-seq data sets, 3 real data sets and mathematical analyses to identify features of methods that allow some to perform better than the others. We prove that methods that explicitly combine the signals from ChIP and input samples are less powerful than methods that do not. Methods that use windows of different sizes are more powerful than the ones that do not. For statistical testing of candidate peaks, methods that use a Poisson test to rank their candidate peaks are more powerful than those that use a Binomial test. BCP and MACS2 have the best operating characteristics on simulated transcription factor binding data. GEM has the highest fraction of the top 500 peaks containing the binding motif of the immunoprecipitated factor, with 50% of its peaks within 10 base pairs of a motif. BCP and MUSIC perform best on histone data. These findings provide guidance and rationale for selecting the best peak caller for a given application.

Enhancers frequently contain multiple binding sites for the same transcription factor. These homotypic binding sites often exhibit synergy, whereby the transcriptional output from two or more binding sites is greater than the sum of the contributions of the individual binding sites alone. Although this phenomenon is frequently observed, the mechanistic basis for homotypic binding site synergy is poorly understood. Here, we identify a bona fide cardiac-specific Prkaa2 enhancer that is synergistically activated by homotypic MEF2 binding sites. We show that two MEF2 sites in the enhancer function cooperatively due to bridging of the MEF2C-bound sites by the SAP domain-containing co-activator protein myocardin, and we show that paired sites buffer the enhancer from integration site-dependent effects on transcription in vivo Paired MEF2 sites are prevalent in cardiac enhancers, suggesting that this might be a common mechanism underlying synergy in the control of cardiac gene expression in vivo.

BackgroundBlood miRNAs are a new promising area of disease research, but variability in miRNA measurements may limit detection of true-positive findings. Here, we measured sources of miRNA variability and determine whether repeated measures can improve power to detect fold-change differences between comparison groups.MethodsBlood from healthy volunteers (N = 12) was collected at three time points. The miRNAs were extracted by a method predetermined to give the highest miRNA yield. Nine different miRNAs were quantified using different qPCR assays and analyzed using mixed models to identify sources of variability. A larger number of miRNAs from a publicly available blood miRNA microarray dataset with repeated measures were used for a bootstrapping procedure to investigate effects of repeated measures on power to detect fold changes in miRNA expression for a theoretical case-control study.ResultsTechnical variability in qPCR replicates was identified as a significant source of variability (P < 0.05) for all nine miRNAs tested. Variability was larger in the TaqMan qPCR assays (SD = 0.15-0.61) versus the qScript qPCR assays (SD = 0.08-0.14). Inter- and intraindividual and extraction variability also contributed significantly for two miRNAs. The bootstrapping procedure demonstrated that repeated measures (20%-50% of N) increased detection of a 2-fold change for approximately 10% to 45% more miRNAs.ConclusionStatistical power to detect small fold changes in blood miRNAs can be improved by accounting for sources of variability using repeated measures and choosing appropriate methods to minimize variability in miRNA quantification.ImpactThis study demonstrates the importance of including repeated measures in experimental designs for blood miRNA research. See all the articles in this CEBP Focus section, "Biomarkers, Biospecimens, and New Technologies in Molecular Epidemiology."

Benzene, a ubiquitous environmental pollutant, causes acute myeloid leukemia (AML). Recently, through transcriptome profiling of peripheral blood mononuclear cells (PBMC), we reported dose-dependent effects of benzene exposure on gene expression and biochemical pathways in 83 workers exposed across four airborne concentration ranges (from 10 ppm) compared with 42 subjects with non-workplace ambient exposure levels. Here, we further characterize these dose-dependent effects with continuous benzene exposure in all 125 study subjects. We estimated air benzene exposure levels in the 42 environmentally-exposed subjects from their unmetabolized urinary benzene levels. We used a novel non-parametric, data-adaptive model selection method to estimate the change with dose in the expression of each gene. We describe non-parametric approaches to model pathway responses and used these to estimate the dose responses of the AML pathway and 4 other pathways of interest. The response patterns of majority of genes as captured by mean estimates of the first and second principal components of the dose-response for the five pathways and the profiles of 6 AML pathway response-representative genes (identified by clustering) exhibited similar apparent supra-linear responses. Responses at or below 0.1 ppm benzene were observed for altered expression of AML pathway genes and CYP2E1. Together, these data show that benzene alters disease-relevant pathways and genes in a dose-dependent manner, with effects apparent at doses as low as 100 ppb in air. Studies with extensive exposure assessment of subjects exposed in the low-dose range between 10 ppb and 1 ppm are needed to confirm these findings.