3D promoter architecture re-organization during iPSC-derived neuronal cell differentiation implicates target genes for neurodevelopmental disorders

Dynamic chromatin architecture identifies new autoimmune-associated enhancers for IL2 and novel genes regulating CD4+ T cell activation

Genome-wide association studies (GWAS) have identified hundreds of genetic signals associated with autoimmune disease. The majority of these signals are located in non-coding regions and likely impact cis-regulatory elements (cRE). Because cRE function is dynamic across cell types and states, profiling the epigenetic status of cRE across physiological processes is necessary to characterize the molecular mechanisms by which autoimmune variants contribute to disease risk. We localized risk variants from 15 autoimmune GWAS to cRE active during TCR-CD28 co-stimulation of naïve human CD4+ T cells. To characterize how dynamic changes in gene expression correlate with cRE activity, we measured transcript levels, chromatin accessibility, and promoter-cRE contacts across three phases of naive CD4+ T cell activation using RNA-seq, ATAC-seq, and HiC. We identified ~1200 protein-coding genes physically connected to accessible disease-associated variants at 423 GWAS signals, at least one-third of which are dynamically regulated by activation. From these maps, we functionally validated a novel stretch of evolutionarily conserved intergenic enhancers whose activity is required for activation-induced IL2 gene expression in human and mouse, and is influenced by autoimmune-associated genetic variation. The set of genes implicated by this approach are enriched for genes controlling CD4+ T cell function and genes involved in human inborn errors of immunity, and we pharmacologically validated eight implicated genes as novel regulators of T cell activation. These studies directly show how autoimmune variants and the genes they regulate influence processes involved in CD4+ T cell proliferation and activation.

Implicating type 2 diabetes effector genes in relevant metabolic cellular models using promoter-focused Capture-C

AIMS/HYPOTHESIS

Genome-wide association studies (GWAS) have identified hundreds of type 2 diabetes loci, with the vast majority of signals located in non-coding regions; as a consequence, it remains largely unclear which 'effector' genes these variants influence. Determining these effector genes has been hampered by the relatively challenging cellular settings in which they are hypothesised to confer their effects.

METHODS

To implicate such effector genes, we elected to generate and integrate high-resolution promoter-focused Capture-C, assay for transposase-accessible chromatin with sequencing (ATAC-seq) and RNA-seq datasets to characterise chromatin and expression profiles in multiple cell lines relevant to type 2 diabetes for subsequent functional follow-up analyses: EndoC-BH1 (pancreatic beta cell), HepG2 (hepatocyte) and Simpson-Golabi-Behmel syndrome (SGBS; adipocyte).

RESULTS

The subsequent variant-to-gene analysis implicated 810 candidate effector genes at 370 type 2 diabetes risk loci. Using partitioned linkage disequilibrium score regression, we observed enrichment for type 2 diabetes and fasting glucose GWAS loci in promoter-connected putative cis-regulatory elements in EndoC-BH1 cells as well as fasting insulin GWAS loci in SGBS cells. Moreover, as a proof of principle, when we knocked down expression of the SMCO4 gene in EndoC-BH1 cells, we observed a statistically significant increase in insulin secretion.

CONCLUSIONS/INTERPRETATION

These results provide a resource for comparing tissue-specific data in tractable cellular models as opposed to relatively challenging primary cell settings.

DATA AVAILABILITY

Raw and processed next-generation sequencing data for EndoC-BH1, HepG2, SGBS_undiff and SGBS_diff cells are deposited in GEO under the Superseries accession GSE262484. Promoter-focused Capture-C data are deposited under accession GSE262496. Hi-C data are deposited under accession GSE262481. Bulk ATAC-seq data are deposited under accession GSE262479. Bulk RNA-seq data are deposited under accession GSE262480.

Variant-to-function mapping of late-onset Alzheimer's disease GWAS signals in human microglial cell models implicates RTFDC1 at the CASS4 locus

Late-onset Alzheimer's disease (LOAD) research has principally focused on neurons over the years due to their known role in the production of amyloid beta plaques and neurofibrillary tangles. In contrast, recent genomic studies of LOAD have implicated microglia as culprits of the prolonged inflammation exacerbating the neurodegeneration observed in patient brains. Indeed, recent LOAD genome-wide association studies (GWAS) have reported multiple loci near genes related to microglial function, including TREM2, ABI3, and CR1. However, GWAS alone cannot pinpoint underlying causal variants or effector genes at such loci, as most signals reside in non-coding regions of the genome and could presumably confer their influence frequently via long-range regulatory interactions. We elected to carry out a combination of ATAC-seq and high-resolution promoter-focused Capture-C in two human microglial cell models (iPSC-derived microglia and HMC3) in order to physically map interactions between LOAD GWAS-implicated candidate causal variants and their corresponding putative effector genes. Notably, we observed consistent evidence that rs6024870 at the GWAS CASS4 locus contacted the promoter of nearby gene, RTFDC1. We subsequently observed a directionallly consistent decrease in RTFDC1 expression with the the protective minor A allele of rs6024870 via both luciferase assays in HMC3 cells and expression studies in primary human microglia. Through CRISPR-Cas9-mediated deletion of the putative regulatory region harboring rs6024870 in HMC3 cells, we observed increased pro-inflammatory cytokine secretion and decreased DNA double strand break repair related, at least in part, to RTFDC1 expression levels. Our variant-to-function approach therefore reveals that the rs6024870-harboring regulatory element at the LOAD 'CASS4' GWAS locus influences both microglial inflammatory capacity and DNA damage resolution, along with cumulative evidence implicating RTFDC1 as a novel candidate effector gene.

Shared and unique 3D genomic features of substance use disorders across multiple cell types

Recent genome-wide association studies (GWAS) have revealed shared genetic components among alcohol, opioid, tobacco and cannabis use disorders. However, the extent of the underlying shared causal variants and effector genes, along with their cellular context, remain unclear. We leveraged our existing 3D genomic datasets comprising high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq and RNA-seq across >50 diverse human cell types to focus on genomic regions that coincide with GWAS loci. Using stratified LD regression, we determined the proportion of genomewide SNP heritability attributable to the features assayed across our cell types by integrating recent GWAS summary statistics for the relevant traits: alcohol use disorder (AUD), tobacco use disorder (TUD), opioid use disorder (OUD) and cannabis use disorder (CanUD). Statistically significant enrichments (P<0.05) were observed in 14 specific cell types, with heritability reaching 9.2-fold for iPSC-derived cortical neurons and neural progenitors, confirming that they are crucial cell types for further functional exploration. Additionally, several pancreatic cell types, notably pancreatic beta cells, showed enrichment for TUD, with heritability enrichments up to 4.8-fold, suggesting genomic overlap with metabolic processes. Further investigation revealed significant positive genetic correlations between T2D with both TUD and CanUD (FDR<0.05) and a significant negative genetic correlation with AUD. Interestingly, after partitioning the heritability for each cell type's cis-regulatory elements, the correlation between T2D and TUD for pancreatic beta cells was greater (r=0.2) than the global genetic correlation value. Our study provides new genomic insights into substance use disorders and implicates cell types where functional follow-up studies could reveal causal variant-gene mechanisms underpinning these disorders.

MotifbreakR v2: extended capability and database integration

MotifbreakR is a software tool that scans genetic variants against position weight matrices of transcription factors (TF) to determine the potential for the disruption of TF binding at the site of the variant. It leverages the Bioconductor suite of software packages and annotations to operate across a diverse array of genomes and motif databases. Initially developed to interrogate the effect of single nucleotide variants (common and rare SNVs) on potential TF binding sites, in motifbreakR v2, we have updated the functionality. New features include the ability to query other types of more complex genetic variants, such as short insertions and deletions (indels). This function allows modeling a more extensive array of variants that may have more significant effects on TF binding. Additionally, while TF binding is based partly on sequence preference, predictions of TF binding based on sequence preference alone can indicate many more potential binding events than observed. Adding information from DNA-binding sequencing datasets lends confidence to motif disruption prediction by demonstrating TF binding in cell lines and tissue types. Therefore, motifbreakR implements querying the ReMap2022 database for evidence that a TF matching the disrupted motif binds over the disrupting variant. Finally, in motifbreakR, in addition to the existing interface, we have implemented an R/Shiny graphical user interface to simplify and enhance access to researchers with different skill sets.

Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2

The ch12q13 locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via cis-regulation. We implicated rs7132908 as a putative causal variant by leveraging our in-house 3D genomic data and public domain datasets. Using a luciferase reporter assay, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. We generated isogenic human embryonic stem cell lines homozygous for either rs7132908 allele to assess changes in gene expression and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. The rs7132908 obesity risk allele influenced expression of FAIM2 and other genes and decreased the proportion of neurons produced by differentiation. We have functionally validated rs7132908 as a causal obesity variant that temporally regulates nearby effector genes and influences neurodevelopment and survival.

GWAS-informed data integration and non-coding CRISPRi screen illuminate genetic etiology of bone mineral density

Over 1,100 independent signals have been identified with genome-wide association studies (GWAS) for bone mineral density (BMD), a key risk factor for mortality-increasing fragility fractures; however, the effector gene(s) for most remain unknown. Informed by a variant-to-gene mapping strategy implicating 89 non-coding elements predicted to regulate osteoblast gene expression at BMD GWAS loci, we executed a single-cell CRISPRi screen in human fetal osteoblast 1.19 cells (hFOBs). The BMD relevance of hFOBs was supported by heritability enrichment from cross-cell type stratified LD-score regression involving 98 cell types grouped into 15 tissues. 24 genes showed perturbation in the screen, with four (ARID5B, CC2D1B, EIF4G2, and NCOA3) exhibiting consistent effects upon siRNA knockdown on three measures of osteoblast maturation and mineralization. Lastly, additional heritability enrichments, genetic correlations, and multi-trait fine-mapping revealed that many BMD GWAS signals are pleiotropic and likely mediate their effects via non-bone tissues that warrant attention in future screens.

Variant to gene mapping for carpal tunnel syndrome risk loci implicates skeletal muscle regulatory elements

BACKGROUND

Carpal tunnel syndrome (CTS) is a common disorder caused by compression of the median nerve in the wrist, resulting in pain and numbness throughout the hand and forearm. While multiple behavioural and physiological factors influence CTS risk, a growing body of evidence supports a strong genetic contribution. Recent genome-wide association study (GWAS) efforts have reported 53 independent signals associated with CTS. While GWAS can identify genetic loci conferring risk, it does not determine which cell types drive the genetic aetiology of the trait, which variants are "causal" at a given signal, and which effector genes correspond to these non-coding variants. These obstacles limit interpretation of potential disease mechanisms.

METHODS

We analysed CTS GWAS findings in the context of chromatin conformation between gene promoters and accessible chromatin regions across cellular models of bone, skeletal muscle, adipocytes and neurons. We identified proxy variants in high LD with the lead CTS sentinel SNPs residing in promoter connected open chromatin in the skeletal muscle and bone contexts.

FINDINGS

We detected significant enrichment for heritability in skeletal muscle myotubes, as well as a weaker correlation in human mesenchymal stem cell-derived osteoblasts. In myotubes, our approach implicated 117 genes contacting 60 proxy variants corresponding to 20 of the 53 GWAS signals. In the osteoblast context we implicated 30 genes contacting 24 proxy variants coinciding with 12 signals, of which 19 genes shared. We subsequently prioritized BZW2 as a candidate effector gene in CTS and implicated it as novel gene that perturbs myocyte differentiation in vitro.

INTERPRETATION

Taken together our results suggest that the CTS genetic component influences the size, integrity, and organization of multiple tissues surrounding the carpal tunnel, in particular muscle and bone, to predispose the nerve to being compressed in this disease setting.

FUNDING

This work was supported by NIH Grant UM1 DK126194 (SFAG and WY), R01AG072705 (SFAG & KDH) and the Center for Spatial and Functional Genomics at CHOP (SFAG & ADW). SFAG is supported by the Daniel B. Burke Endowed Chair for Diabetes Research. WY is supported by the Perelman School of Medicine of the University of Pennsylvania.

Deficiency of the Heterogeneous Nuclear Ribonucleoprotein U locus leads to delayed hindbrain neurogenesis

Genetic variants affecting Heterogeneous Nuclear Ribonucleoprotein U (HNRNPU) have been identified in several neurodevelopmental disorders (NDDs). HNRNPU is widely expressed in the human brain and shows the highest postnatal expression in the cerebellum. Recent studies have investigated the role of HNRNPU in cerebral cortical development, but the effects of HNRNPU deficiency on cerebellar development remain unknown. Here, we describe the molecular and cellular outcomes of HNRNPU locus deficiency during in vitro neural differentiation of patient-derived and isogenic neuroepithelial stem cells with a hindbrain profile. We demonstrate that HNRNPU deficiency leads to chromatin remodeling of A/B compartments, and transcriptional rewiring, partly by impacting exon inclusion during mRNA processing. Genomic regions affected by the chromatin restructuring and host genes of exon usage differences show a strong enrichment for genes implicated in epilepsies, intellectual disability, and autism. Lastly, we show that at the cellular level HNRNPU downregulation leads to an increased fraction of neural progenitors in the maturing neuronal population. We conclude that the HNRNPU locus is involved in delayed commitment of neural progenitors to differentiate in cell types with hindbrain profile.

Sex Significantly Impacts the Function of Major Depression-Linked Variants In Vivo

BACKGROUND

Genome-wide association studies have discovered blocks of common variants-likely transcriptional-regulatory-associated with major depressive disorder (MDD), though the functional subset and their biological impacts remain unknown. Likewise, why depression occurs in females more frequently than males is unclear. We therefore tested the hypothesis that risk-associated functional variants interact with sex and produce greater impact in female brains.

METHODS

We developed techniques to directly measure regulatory variant activity and sex interactions using massively parallel reporter assays in the mouse brain in vivo, in a cell type-specific manner, and applied these approaches to measure activity of >1000 variants from >30 MDD loci.

RESULTS

We identified extensive sex-by-allele effects in mature hippocampal neurons, suggesting that sex-differentiated impacts of genetic risk may underlie sex bias in disease. Unbiased informatics approaches indicated that functional MDD variants recurrently disrupt a number of transcription factor binding motifs, including those of sex hormone receptors. We confirmed a role for the latter by performing massively parallel reporter assays in neonatal mice on the day of birth (during a sex-differentiating hormone surge) and hormonally quiescent juveniles.

CONCLUSIONS

Our study provides novel insights into the influence of age, biological sex, and cell type on regulatory variant function and provides a framework for in vivo parallel assays to functionally define interactions between organismal variables such as sex and regulatory variation. Moreover, we experimentally demonstrate that a portion of the sex differences seen in MDD occurrence may be a product of sex-differentiated effects at associated regulatory variants.

Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2

The ch12q13 obesity locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via an influence on cis-regulation within the genomic region. We implicated rs7132908 as a putative causal variant at this locus leveraging a combination of our inhouse 3D genomic data, public domain datasets, and several computational approaches. Using a luciferase reporter assay in human primary astrocytes, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. Motivated by this finding, we went on to generate isogenic human embryonic stem cell lines homozygous for either rs7132908 allele with CRISPR-Cas9 homology-directed repair to assess changes in gene expression due to genotype and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. We observed that the rs7132908 obesity risk allele influenced the expression of FAIM2 along with other genes, decreased the proportion of neurons produced during differentiation, up-regulated cell death gene sets, and conversely down-regulated neuron differentiation gene sets. We have therefore functionally validated rs7132908 as a causal obesity variant which temporally regulates nearby effector genes at the ch12q13 locus and influences neurodevelopment and survival.

Neurodevelopmental disorders, like cancer, are connected to impaired chromatin remodelers, PI3K/mTOR, and PAK1-regulated MAPK

Neurodevelopmental disorders (NDDs) and cancer share proteins, pathways, and mutations. Their clinical symptoms are different. However, individuals with NDDs have higher probabilities of eventually developing cancer. Here, we review the literature and ask how the shared features can lead to different medical conditions and why having an NDD first can increase the chances of malignancy. To explore these vital questions, we focus on dysregulated PI3K/mTOR, a major brain cell growth pathway in differentiation, and MAPK, a critical pathway in proliferation, a hallmark of cancer. Differentiation is governed by chromatin organization, making aberrant chromatin remodelers highly likely agents in NDDs. Dysregulated chromatin organization and accessibility influence the lineage of specific cell brain types at specific embryonic development stages. PAK1, with pivotal roles in brain development and in cancer, also regulates MAPK. We review, clarify, and connect dysregulated pathways with dysregulated proliferation and differentiation in cancer and NDDs and highlight PAK1 role in brain development and MAPK regulation. Exactly how PAK1 activation controls brain development, and why specific chromatin remodeler components, e.g., BAF170 encoded by SMARCC2 in autism, await clarification.

Technologies, strategies, and cautions when deconvoluting genome-wide association signals: FTO in focus

Genome-wide association studies have revealed a plethora of genetic variants that correlate with polygenic conditions. However, causal molecular mechanisms have proven challenging to fully define. Without such information, the associations are not physiologically useful or clinically actionable. By reviewing studies of the FTO locus in the genetic etiology of obesity, we wish to highlight advances in the field fueled by the evolution of technical and analytic strategies in assessing the molecular bases for genetic associations. Particular attention is drawn to extrapolating experimental findings from animal models and cell types to humans, as well as technical aspects used to identify long-range DNA interactions and their biological relevance with regard to the associated trait. A unifying model is proposed by which independent obesogenic pathways regulated by multiple FTO variants and genes are integrated at the primary cilium, a cellular antenna where signaling molecules that control energy balance convene.

Variant-to-gene mapping followed by cross-species genetic screening identifies GPI-anchor biosynthesis as a regulator of sleep

Genome-wide association studies (GWAS) in humans have identified loci robustly associated with several heritable diseases or traits, yet little is known about the functional roles of the underlying causal variants in regulating sleep duration or quality. We applied an ATAC-seq/promoter focused Capture C strategy in human iPSC-derived neural progenitors to carry out a "variant-to-gene" mapping campaign that identified 88 candidate sleep effector genes connected to relevant GWAS signals. To functionally validate the role of the implicated effector genes in sleep regulation, we performed a neuron-specific RNA interference screen in the fruit fly, Drosophila melanogaster, followed by validation in zebrafish. This approach identified a number of genes that regulate sleep including a critical role for glycosylphosphatidylinositol (GPI)-anchor biosynthesis. These results provide the first physical variant-to-gene mapping of human sleep genes followed by a model organism-based prioritization, revealing a conserved role for GPI-anchor biosynthesis in sleep regulation.

Abnormal Chromatin Folding in the Molecular Pathogenesis of Epilepsy and Autism Spectrum Disorder: a Meta-synthesis with Systematic Searching

How DNA is folded and packaged in nucleosomes is an essential regulator of gene expression. Abnormal patterns of chromatin folding are implicated in a wide range of diseases and disorders, including epilepsy and autism spectrum disorder (ASD). These disorders are thought to have a shared pathogenesis involving an imbalance in the number of excitatory-inhibitory neurons formed during neurodevelopment; however, the underlying pathological mechanism behind this imbalance is poorly understood. Studies are increasingly implicating abnormal chromatin folding in neural stem cells as one of the candidate pathological mechanisms, but no review has yet attempted to summarise the knowledge in this field. This meta-synthesis is a systematic search of all the articles on epilepsy, ASD, and chromatin folding. Its two main objectives were to determine to what extent abnormal chromatin folding is implicated in the pathogenesis of epilepsy and ASD, and secondly how abnormal chromatin folding leads to pathological disease processes. This search produced 22 relevant articles, which together strongly implicate abnormal chromatin folding in the pathogenesis of epilepsy and ASD. A range of mutations and chromosomal structural abnormalities lead to this effect, including single nucleotide polymorphisms, copy number variants, translocations and mutations in chromatin modifying. However, knowledge is much more limited into how abnormal chromatin organisation subsequently causes pathological disease processes, not yet showing, for example, whether it leads to abnormal excitation-inhibitory neuron imbalance in human brain organoids.

Disruption of the gene regulatory programme in neurodevelopmental disorders

Cortical development consists of a series of synchronised events, including fate transition of cortical progenitors, neuronal migration, specification and connectivity. It is becoming clear that gene expression programs governing these events rely on the interplay between signalling molecules, transcription factors and epigenetic mechanisms. When genetic or environmental factors disrupt expression of genes involved in important brain development processes, neurodevelopmental disorders can occur. This review aims to highlight how recent advances in technologies have helped uncover and imitate the gene regulatory mechanisms commonly disrupted in neurodevelopmental disorders.

Identifying differential regulatory control of APOE ɛ4 on African versus European haplotypes as potential therapeutic targets

We previously demonstrated that in Alzheimer's disease (AD) patients, European apolipoprotein E (APOE) ε4 carriers express significantly more APOE ε4 in their brains than African AD carriers. We examined single nucleotide polymorphisms near APOE with significant frequency differences between African and European/Japanese APOE ε4 haplotypes that could contribute to this difference in expression through regulation. Two enhancer massively parallel reporter assay (MPRA) approaches were performed, supplemented with single fragment reporter assays. We used Capture C analyses to support interactions with the APOE promoter. Introns within TOMM40 showed increased enhancer activity in the European/Japanese versus African haplotypes in astrocytes and microglia. This region overlaps with APOE promoter interactions as assessed by Capture C analysis. Single variant analyses pinpoints rs2075650/rs157581, and rs59007384 as functionally different on these haplotypes. Identification of the mechanisms for differential regulatory function for APOE expression between African and European/Japanese haplotypes could lead to therapeutic targets for APOE ε4 carriers.

3D chromatin maps of the human pancreas reveal lineage-specific regulatory architecture of T2D risk

Three-dimensional (3D) chromatin organization maps help dissect cell-type-specific gene regulatory programs. Furthermore, 3D chromatin maps contribute to elucidating the pathogenesis of complex genetic diseases by connecting distal regulatory regions and genetic risk variants to their respective target genes. To understand the cell-type-specific regulatory architecture of diabetes risk, we generated transcriptomic and 3D epigenomic profiles of human pancreatic acinar, alpha, and beta cells using single-cell RNA-seq, single-cell ATAC-seq, and high-resolution Hi-C of sorted cells. Comparisons of these profiles revealed differential A/B (open/closed) chromatin compartmentalization, chromatin looping, and transcriptional factor-mediated control of cell-type-specific gene regulatory programs. We identified a total of 4,750 putative causal-variant-to-target-gene pairs at 194 type 2 diabetes GWAS signals using pancreatic 3D chromatin maps. We found that the connections between candidate causal variants and their putative target effector genes are cell-type stratified and emphasize previously underappreciated roles for alpha and acinar cells in diabetes pathogenesis.

Implicating effector genes at COVID-19 GWAS loci using promoter-focused Capture-C in disease-relevant immune cell types

BACKGROUND

SARS-CoV-2 infection results in a broad spectrum of COVID-19 disease, from mild or no symptoms to hospitalization and death. COVID-19 disease severity has been associated with some pre-existing conditions and the magnitude of the adaptive immune response to SARS-CoV-2, and a recent genome-wide association study (GWAS) of the risk of critical illness revealed a significant genetic component. To gain insight into how human genetic variation attenuates or exacerbates disease following SARS-CoV-2 infection, we implicated putatively functional COVID risk variants in the cis-regulatory landscapes of human immune cell types with established roles in disease severity and used high-resolution chromatin conformation capture to map these disease-associated elements to their effector genes.

RESULTS

This functional genomic approach implicates 16 genes involved in viral replication, the interferon response, and inflammation. Several of these genes (PAXBP1, IFNAR2, OAS1, OAS3, TNFAIP8L1, GART) were differentially expressed in immune cells from patients with severe versus moderate COVID-19 disease, and we demonstrate a previously unappreciated role for GART in T cell-dependent antibody-producing B cell differentiation in a human tonsillar organoid model.

CONCLUSIONS

This study offers immunogenetic insight into the basis of COVID-19 disease severity and implicates new targets for therapeutics that limit SARS-CoV-2 infection and its resultant life-threatening inflammation.

Whole Genome DNA Methylation Profiling of D2 Medium Spiny Neurons in Mouse Nucleus Accumbens Using Two Independent Library Preparation Methods

DNA methylation plays essential roles in various cellular processes. Next-generation sequencing has enabled us to study the functional implication of DNA methylation across the whole genome. However, this approach usually requires a substantial amount of genomic DNA, which limits its application to defined cell types within a discrete brain region. Here, we applied two separate protocols, Accel-NGS Methyl-Seq (AM-seq) and Enzymatic Methyl-seq (EM-seq), to profile the methylome of D2 dopamine receptor-expressing medium spiny neurons (D2-MSNs) in mouse nucleus accumbens (NAc). Using 40 ng DNA extracted from FACS-isolated D2-MSNs, we found that both methods yielded comparably high-quality methylome data. Additionally, we identified numerous unmethylated regions (UMRs) as cell type-specific regulatory regions. By comparing the NAc D2-MSN methylome with the published methylomes of mouse prefrontal cortex excitatory neurons and neural progenitor cells (NPCs), we identified numerous differentially methylated CpG and non-CpG regions. Our study not only presents a comparison of these two low-input DNA whole genome methylation profiling protocols, but also provides a resource of DNA methylome of mouse accumbal D2-MSNs, a neuron type that has critical roles in addiction and other neuropsychiatric disorders.

How can same-gene mutations promote both cancer and developmental disorders?

The question of how same-gene mutations can drive both cancer and neurodevelopmental disorders has been puzzling. It has also been puzzling why those with neurodevelopmental disorders have a high risk of cancer. Ras, MEK, PI3K, PTEN, and SHP2 are among the oncogenic proteins that can harbor mutations that encode diseases other than cancer. Understanding why some of their mutations can promote cancer, whereas others promote neurodevelopmental diseases, and why even the same mutations may promote both phenotypes, has important clinical ramifications. Here, we review the literature and address these tantalizing questions. We propose that cell type–specific expression of the mutant protein, and of other proteins in the respective pathway, timing of activation (during embryonic development or sporadic emergence), and the absolute number of molecules that the mutations activate, alone or in combination, are pivotal in determining the pathological phenotypes—cancer and (or) developmental disorders.

Cis-regulatory architecture of human ESC-derived hypothalamic neuron differentiation aids in variant-to-gene mapping of relevant complex traits

The hypothalamus regulates metabolic homeostasis by influencing behavior and endocrine systems. Given its role governing key traits, such as body weight and reproductive timing, understanding the genetic regulation of hypothalamic development and function could yield insights into disease pathogenesis. However, given its inaccessibility, studying human hypothalamic gene regulation has proven challenging. To address this gap, we generate a high-resolution chromatin architecture atlas of an established embryonic stem cell derived hypothalamic-like neuron model across three stages of in vitro differentiation. We profile accessible chromatin and identify physical contacts between gene promoters and putative cis-regulatory elements to characterize global regulatory landscape changes during hypothalamic differentiation. Next, we integrate these data with GWAS loci for various complex traits, identifying multiple candidate effector genes. Our results reveal common target genes for these traits, potentially affecting core developmental pathways. Our atlas will enable future efforts to determine hypothalamic mechanisms influencing disease susceptibility.

Maternal immune activation primes deficiencies in adult hippocampal neurogenesis

Neurogenesis, the process in which new neurons are generated, occurs throughout life in the mammalian hippocampus. Decreased adult hippocampal neurogenesis (AHN) is a common feature across psychiatric disorders, including schizophrenia, depression- and anxiety-related behaviours, and is highly regulated by environmental influences. Epidemiological studies have consistently implicated maternal immune activation (MIA) during neurodevelopment as a risk factor for psychiatric disorders in adulthood. The extent to which the reduction of hippocampal neurogenesis in adulthood may be driven by early life exposures, such as MIA, is however unclear. We therefore reviewed the literature for evidence of the involvement of MIA in disrupting AHN. Consistent with our hypothesis, data from both in vivo murine and in vitro human models of AHN provide evidence for key roles of specific cytokines induced by MIA in the foetal brain in disrupting hippocampal neural progenitor cell proliferation and differentiation early in development. The precise molecular mechanisms however remain unclear. Nonetheless, these data suggest a potential latent vulnerability mechanism, whereby MIA primes dysfunction in the unique hippocampal pool of neural stem/progenitor cells. This renders offspring potentially more susceptible to additional environmental exposures later in life, such as chronic stress, resulting in the unmasking of psychopathology. We highlight the need for studies to test this hypothesis using validated animal models of MIA, but also to test the relevance of such data for human pathology at a molecular basis through the use of patient-derived induced pluripotent stem cells (hiPSC) differentiated into hippocampal progenitor cells.