Glial dysregulation in the human brain in fragile X-associated tremor/ataxia syndrome.

Short trinucleotide expansions at the FMR1 locus are associated with the late-onset condition fragile X-associated tremor/ataxia syndrome (FXTAS), which shows very different clinical and pathological features from fragile X syndrome (associated with longer expansions), with no clear molecular explanation for these marked differences. One prevailing theory posits that the shorter, premutation expansion uniquely causes extreme neurotoxic increases in FMR1 mRNA (i.e., four to eightfold increases), but evidence to support this hypothesis is largely derived from analysis of peripheral blood. We applied single-nucleus RNA sequencing to postmortem frontal cortex and cerebellum from 7 individuals with premutation and matched controls (n = 6) to assess cell type-specific molecular neuropathology. We found only modest upregulation (~1.3-fold) of FMR1 in some glial populations associated with premutation expansions. In premutation cases, we also identified decreased astrocyte proportions in the cortex. Differential expression and gene ontology analysis demonstrated altered neuroregulatory roles of glia. Using network analyses, we identified cell type-specific and region-specific patterns of FMR1 protein target gene dysregulation unique to premutation cases, with notable network dysregulation in the cortical oligodendrocyte lineage. We used pseudotime trajectory analysis to determine how oligodendrocyte development was altered and identified differences in early gene expression in oligodendrocyte trajectories in premutation cases specifically, implicating early cortical glial developmental perturbations. These findings challenge dogma regarding extremely elevated FMR1 increases in FXTAS and implicate glial dysregulation as a critical facet of premutation pathophysiology, representing potential unique therapeutic targets directly derived from the human condition.

The SET oncoprotein promotes estrogen-induced transcription by facilitating establishment of active chromatin.

SET is a multifunctional histone-binding oncoprotein that regulates transcription by an unclear mechanism. Here we show that SET enhances estrogen-dependent transcription. SET knockdown abrogates transcription of estrogen-responsive genes and their enhancer RNAs. In response to 17ß-estradiol (E2), SET binds to the estrogen receptor α (ERα) and is recruited to ERα-bound enhancers and promoters at estrogen response elements (EREs). SET functions as a histone H2 chaperone that dynamically associates with H2A.Z via its acidic C-terminal domain and promotes H2A.Z incorporation, ERα, MLL1, and KDM3A loading and modulates histone methylation at EREs. SET depletion diminishes recruitment of condensin complexes to EREs and impairs E2-dependent enhancer-promoter looping. Thus, SET boosts E2-induced gene expression by establishing an active chromatin structure at ERα-bound enhancers and promoters, which is essential for transcriptional activation.

Single-cell transcriptome profile of mouse skin undergoing antigen-driven allergic inflammation recapitulates findings in atopic dermatitis skin lesions.

BACKGROUND: Allergic skin inflammation elicited in mice by epicutaneous (EC) sensitization with antigen shares characteristics with human atopic dermatitis (AD). OBJECTIVE: We characterized gene expression by single cells in mouse skin undergoing antigen-driven allergic inflammation and compared the results with findings in AD skin lesions. METHODS: Mice were EC sensitized by application of ovalbumin (OVA) or saline to tape-stripped skin. Single-cell RNA sequencing was performed on skin cells 12 days later. Flow cytometry analysis was performed to validate results. RESULTS: Sequencing identified 7 nonhematopoietic and 6 hematopoietic cell subsets in EC-sensitized mouse skin. OVA sensitization resulted in the expansion in the skin of T cells, dendritic cells, macrophages, mast cells/basophils, fibroblasts, and myocytes cell clusters, and in upregulation of TH2 cytokine gene expression in CD4+ T cells and mast cells/basophils. Genes differentially expressed in OVA-sensitized skin included genes important for inflammation in dendritic cells and macrophages, collagen deposition, and leukocyte migration in fibroblasts, chemotaxis in endothelial cells and skin barrier integrity, and differentiation in KCs-findings that recapitulate those in AD skin lesions. Unexpectedly, mast cells/basophils, rather than T cells, were the major source of Il4 and ll13 in OVA-sensitized mouse skin. In addition, our results suggest novel pathways in fibroblast and endothelial cells that may contribute to allergic skin inflammation. CONCLUSION: The gene expression profile of single cells in mouse skin undergoing antigen-driven shares many features with that in AD skin lesions and unveils novel pathways that may be involved in allergic skin inflammation.

Evaluation of copy number variants for genetic hearing loss: a review of current approaches and recent findings.

Structural variation includes a change in copy number, orientation, or location of a part of the genome. Copy number variants (CNVs) are a common cause of genetic hearing loss, comprising nearly 20% of diagnosed cases. While large deletions involving the gene STRC are the most common pathogenic CNVs, a significant proportion of known hearing loss genes also contain pathogenic CNVs. In this review, we provide an overview of currently used methods for detection of CNVs in genes known to cause hearing loss including molecular techniques such as multiplex ligation probe amplification (MLPA) and digital droplet polymerase chain reaction (ddPCR), array-CGH and single-nucleotide polymorphism (SNP) arrays, as well as techniques for detection of CNVs using next-generation sequencing data analysis including targeted gene panel, exome, and genome sequencing data. In addition, in this review, we compile published data on pathogenic hearing loss CNVs to provide an up-to-date overview. We show that CNVs have been identified in 29 different non-syndromic hearing loss genes. An understanding of the contribution of CNVs to genetic hearing loss is critical to the current diagnosis of hearing loss and is crucial for future gene therapies. Thus, evaluation for CNVs is required in any modern pipeline for genetic diagnosis of hearing loss.

Congenital X-linked neutropenia with myelodysplasia and somatic tetraploidy due to a germline mutation in SEPT6.

Septins play key roles in mammalian cell division and cytokinesis but have not previously been implicated in a germline human disorder. A male infant with severe neutropenia and progressive dysmyelopoiesis with tetraploid myeloid precursors was identified. No known genetic etiologies for neutropenia or bone marrow failure were found. However, next-generation sequencing of germline samples from the patient revealed a novel, de novo germline stop-loss mutation in the X-linked gene SEPT6 that resulted in reduced SEPT6 staining in bone marrow granulocyte precursors and megakaryocytes. Patient skin fibroblast-derived induced pluripotent stem cells (iPSCs) produced reduced myeloid colonies, particularly of the granulocyte lineage. CRISPR/Cas9 knock-in of the patient's mutation or complete knock-out of SEPT6 was not tolerated in non-patient-derived iPSCs or human myeloid cell lines, but SEPT6 knock-out was successful in an erythroid cell line and resulting clones revealed a propensity to multinucleation. In silico analysis predicts that the mutated protein hinders the dimerization of SEPT6 coiled-coils in both parallel and antiparallel arrangements, which could in turn impair filament formation. These data demonstrate a critical role for SEPT6 in chromosomal segregation in myeloid progenitors that can account for the unusual predisposition to aneuploidy and dysmyelopoiesis.

Mechanisms underlying genetic susceptibility to multisystem inflammatory syndrome in children (MIS-C).

BACKGROUND: Multisystem inflammatory syndrome in children (MIS-C) is a pediatric complication of severe acute respiratory syndrome coronavirus 2 infection that is characterized by multiorgan inflammation and frequently by cardiovascular dysfunction. It occurs predominantly in otherwise healthy children. We previously reported haploinsufficiency of suppressor of cytokine signaling 1 (SOCS1), a negative regulator of type I and II interferons, as a genetic risk factor for MIS-C. OBJECTIVES: We aimed to identify additional genetic mechanisms underlying susceptibility to severe acute respiratory syndrome coronavirus 2-associated MIS-C. METHODS: In a single-center, prospective cohort study, whole exome sequencing was performed on patients with MIS-C. The impact of candidate variants was tested by using patients' PBMCs obtained at least 7 months after recovery. RESULTS: We enrolled 18 patients with MIS-C (median age = 8 years; interquartile range = 5-12.25 years), of whom 89% had no conditions other than obesity. In 2 boys with no significant infection history, we identified and validated hemizygous deleterious defects in XIAP, encoding X-linked inhibitor of apoptosis, and CYBB, encoding cytochrome b-245, beta subunit. Including the previously reported SOCS1 haploinsufficiency, a genetic diagnosis was identified in 3 of 18 patients (17%). In contrast to patients with mild COVID-19, patients with defects in SOCS1, XIAP, or CYBB exhibit an inflammatory immune cell transcriptome with enrichment of differentially expressed genes in pathways downstream of IL-18, oncostatin M, and nuclear factor κB, even after recovery. CONCLUSIONS: Although inflammatory disorders are rare in the general population, our cohort of patients with MIS-C was enriched for monogenic susceptibility to inflammation. Our results support the use of next-generation sequencing in previously healthy children who develop MIS-C.

Maternal Iron Deficiency Modulates Placental Transcriptome and Proteome in Mid-Gestation of Mouse Pregnancy.

BACKGROUND: Maternal iron deficiency (ID) is associated with poor pregnancy and fetal outcomes. The effect is thought to be mediated by the placenta but there is no comprehensive assessment of placental responses to maternal ID. Additionally, whether the influence of maternal ID on the placenta differs by fetal sex is unknown. OBJECTIVES: To identify gene and protein signatures of ID mouse placentas at mid-gestation. A secondary objective was to profile the expression of iron genes in mouse placentas across gestation. METHODS: We used a real-time PCR-based array to determine the mRNA expression of all known iron genes in mouse placentas at embryonic day (E) 12.5, E14.5, E16.5, and E19.5 (n = 3 placentas/time point). To determine the effect of maternal ID, we performed RNA sequencing and proteomics in male and female placentas from ID and iron-adequate mice at E12.5 (n = 8 dams/diet). RESULTS: In female placentas, 6 genes, including transferrin receptor (Tfrc) and solute carrier family 11 member 2, were significantly changed by maternal ID. An additional 154 genes were altered in male ID placentas. A proteomic analysis quantified 7662 proteins in the placenta. Proteins translated from iron-responsive element (IRE)-containing mRNA were altered in abundance; ferritin and ferroportin 1 decreased, while TFRC increased in ID placentas. Less than 4% of the significantly altered genes in ID placentas occurred both at the transcriptional and translational levels. CONCLUSIONS: Our data demonstrate that the impact of maternal ID on placental gene expression in mice is limited in scope and magnitude at mid-gestation. We provide strong evidence for IRE-based transcriptional and translational coordination of iron gene expression in the mouse placenta. Finally, we discover sexually dimorphic effects of maternal ID on placental gene expression, with more genes and pathways altered in male compared with female mouse placentas.

Pioneer factors govern super-enhancer dynamics in stem cell plasticity and lineage choice.

Adult stem cells occur in niches that balance self-renewal with lineage selection and progression during tissue homeostasis. Following injury, culture or transplantation, stem cells outside their niche often display fate flexibility. Here we show that super-enhancers underlie the identity, lineage commitment and plasticity of adult stem cells in vivo. Using hair follicle as a model, we map the global chromatin domains of hair follicle stem cells and their committed progenitors in their native microenvironments. We show that super-enhancers and their dense clusters ('epicentres') of transcription factor binding sites undergo remodelling upon lineage progression. New fate is acquired by decommissioning old and establishing new super-enhancers and/or epicentres, an auto-regulatory process that abates one master regulator subset while enhancing another. We further show that when outside their niche, either in vitro or in wound-repair, hair follicle stem cells dynamically remodel super-enhancers in response to changes in their microenvironment. Intriguingly, some key super-enhancers shift epicentres, enabling their genes to remain active and maintain a transitional state in an ever-changing transcriptional landscape. Finally, we identify SOX9 as a crucial chromatin rheostat of hair follicle stem cell super-enhancers, and provide functional evidence that super-enhancers are dynamic, dense transcription-factor-binding platforms which are acutely sensitive to pioneer master regulators whose levels define not only spatial and temporal features of lineage-status but also stemness, plasticity in transitional states and differentiation.

An H3K36 methylation-engaging Tudor motif of polycomb-like proteins mediates PRC2 complex targeting.

Polycomb repressive complex 2 (PRC2) regulates pluripotency, differentiation, and tumorigenesis through catalysis of histone H3 lysine 27 trimethylation (H3K27me3) on chromatin. However, the mechanisms that underlie PRC2 recruitment and spreading on chromatin remain unclear. Here we report that histone H3 lysine 36 trimethylation (H3K36me3) binding activity is harbored in the Tudor motifs of PRC2-associated polycomb-like (PCL) proteins PHF1/PCL1 and PHF19/PCL3. Ectopically expressed PHF1 induced Tudor-dependent stabilization of PRC2 complexes on bulk chromatin and mediated spreading of PRC2 and H3K27me3 into H3K36me3-containing chromatin regions. In murine pluripotent stem cells, we identified coexistence of H3K36me3, H3K27me3, and PHF19/PCL3 at a subset of poised developmental genes and demonstrated that PHF19/PCL3 Tudor function is required for optimal H3K27me3 and repression of these loci. Collectively, our data suggest that PCL recognition of H3K36me3 promotes intrusion of PRC2 complexes into active chromatin regions to promote gene silencing and modulate the chromatin landscape during development.

Genetic diagnoses in epilepsy: The impact of dynamic exome analysis in a pediatric cohort.

OBJECTIVE: We evaluated the yield of systematic analysis and/or reanalysis of whole exome sequencing (WES) data from a cohort of well-phenotyped pediatric patients with epilepsy and suspected but previously undetermined genetic etiology. METHODS: We identified and phenotyped 125 participants with pediatric epilepsy. Etiology was unexplained at the time of enrollment despite clinical testing, which included chromosomal microarray (57 patients), epilepsy gene panel (n = 48), both (n = 28), or WES (n = 8). Clinical epilepsy diagnoses included developmental and epileptic encephalopathy (DEE), febrile infection-related epilepsy syndrome, Rasmussen encephalitis, and other focal and generalized epilepsies. We analyzed WES data and compared the yield in participants with and without prior clinical genetic testing. RESULTS: Overall, we identified pathogenic or likely pathogenic variants in 40% (50/125) of our study participants. Nine patients with DEE had genetic variants in recently published genes that had not been recognized as epilepsy-related at the time of clinical testing (FGF12, GABBR1, GABBR2, ITPA, KAT6A, PTPN23, RHOBTB2, SATB2), and eight patients had genetic variants in candidate epilepsy genes (CAMTA1, FAT3, GABRA6, HUWE1, PTCHD1). Ninety participants had concomitant or subsequent clinical genetic testing, which was ultimately explanatory for 26% (23/90). Of the 67 participants whose molecular diagnoses were "unsolved" through clinical genetic testing, we identified pathogenic or likely pathogenic variants in 17 (25%). SIGNIFICANCE: Our data argue for early consideration of WES with iterative reanalysis for patients with epilepsy, particularly those with DEE or epilepsy with intellectual disability. Rigorous analysis of WES data of well-phenotyped patients with epilepsy leads to a broader understanding of gene-specific phenotypic spectra as well as candidate disease gene identification. We illustrate the dynamic nature of genetic diagnosis over time, with analysis and in some cases reanalysis of exome data leading to the identification of disease-associated variants among participants with previously nondiagnostic results from a variety of clinical testing strategies.

Selective inhibition of EZH2 and EZH1 enzymatic activity by a small molecule suppresses MLL-rearranged leukemia.

Enhancer of zeste homolog 2 (EZH2) and related EZH1 control gene expression and promote tumorigenesis via methylating histone H3 at lysine 27 (H3K27). These methyltransferases are ideal therapeutic targets due to their frequent hyperactive mutations and overexpression found in cancer, including hematopoietic malignancies. Here, we characterized a set of small molecules that allow pharmacologic manipulation of EZH2 and EZH1, which include UNC1999, a selective inhibitor of both enzymes, and UNC2400, an inactive analog compound useful for assessment of off-target effect. UNC1999 suppresses global H3K27 trimethylation/dimethylation (H3K27me3/2) and inhibits growth of mixed lineage leukemia (MLL)-rearranged leukemia cells. UNC1999-induced transcriptome alterations overlap those following knockdown of embryonic ectoderm development, a common cofactor of EZH2 and EZH1, demonstrating UNC1999's on-target inhibition. Mechanistically, UNC1999 preferentially affects distal regulatory elements such as enhancers, leading to derepression of polycomb targets including Cdkn2a. Gene derepression correlates with a decrease in H3K27me3 and concurrent gain in H3K27 acetylation. UNC2400 does not induce such effects. Oral administration of UNC1999 prolongs survival of a well-defined murine leukemia model bearing MLL-AF9. Collectively, our study provides the detailed profiling for a set of chemicals to manipulate EZH2 and EZH1 and establishes specific enzymatic inhibition of polycomb repressive complex 2 (PRC2)-EZH2 and PRC2-EZH1 by small-molecule compounds as a novel therapeutics for MLL-rearranged leukemia.

Significant expansion of the REST/NRSF cistrome in human versus mouse embryonic stem cells: potential implications for neural development.

Recent studies have employed cross-species comparisons of transcription factor binding, reporting significant regulatory network 'rewiring' between species. Here, we address how a transcriptional repressor targets and regulates neural genes differentially between human and mouse embryonic stem cells (ESCs). We find that the transcription factor, Repressor Element 1 Silencing Transcription factor (REST; also called neuron restrictive silencer factor) binds to a core group of ∼1200 syntenic genomic regions in both species, with these conserved sites highly enriched with co-factors, selective histone modifications and DNA hypomethylation. Genes with conserved REST binding are enriched with neural functions and more likely to be upregulated upon REST depletion. Interestingly, we identified twice as many REST peaks in human ESCs compared to mouse ESCs. Human REST cistrome expansion involves additional peaks in genes targeted by REST in both species and human-specific gene targets. Genes with expanded REST occupancy in humans are enriched for learning or memory functions. Analysis of neurological disorder associated genes reveals that Amyotrophic Lateral Sclerosis and oxidative stress genes are particularly enriched with human-specific REST binding. Overall, our results demonstrate that there is substantial rewiring of human and mouse REST cistromes, and that REST may have human-specific roles in brain development and functions.

Identification of in vivo DNA-binding mechanisms of Pax6 and reconstruction of Pax6-dependent gene regulatory networks during forebrain and lens development.

The transcription factor Pax6 is comprised of the paired domain (PD) and homeodomain (HD). In the developing forebrain, Pax6 is expressed in ventricular zone precursor cells and in specific subpopulations of neurons; absence of Pax6 results in disrupted cell proliferation and cell fate specification. Pax6 also regulates the entire lens developmental program. To reconstruct Pax6-dependent gene regulatory networks (GRNs), ChIP-seq studies were performed using forebrain and lens chromatin from mice. A total of 3514 (forebrain) and 3723 (lens) Pax6-containing peaks were identified, with ∼70% of them found in both tissues and thereafter called 'common' peaks. Analysis of Pax6-bound peaks identified motifs that closely resemble Pax6-PD, Pax6-PD/HD and Pax6-HD established binding sequences. Mapping of H3K4me1, H3K4me3, H3K27ac, H3K27me3 and RNA polymerase II revealed distinct types of tissue-specific enhancers bound by Pax6. Pax6 directly regulates cortical neurogenesis through activation (e.g. Dmrta1 and Ngn2) and repression (e.g. Ascl1, Fezf2, and Gsx2) of transcription factors. In lens, Pax6 directly regulates cell cycle exit via components of FGF (Fgfr2, Prox1 and Ccnd1) and Wnt (Dkk3, Wnt7a, Lrp6, Bcl9l, and Ccnd1) signaling pathways. Collectively, these studies provide genome-wide analysis of Pax6-dependent GRNs in lens and forebrain and establish novel roles of Pax6 in organogenesis.

Chromatin features, RNA polymerase II and the comparative expression of lens genes encoding crystallins, transcription factors, and autophagy mediators.

PURPOSE: Gene expression correlates with local chromatin structure. Our studies have mapped histone post-translational modifications, RNA polymerase II (pol II), and transcription factor Pax6 in lens chromatin. These data represent the first genome-wide insights into the relationship between lens chromatin structure and lens transcriptomes and serve as an excellent source for additional data analysis and refinement. The principal lens proteins, the crystallins, are encoded by predominantly expressed mRNAs; however, the regulatory mechanisms underlying their high expression in the lens remain poorly understood. METHODS: The formaldehyde-assisted identification of regulatory regions (FAIRE-Seq) was employed to analyze newborn lens chromatin. ChIP-seq and RNA-seq data published earlier (GSE66961) have been used to assist in FAIRE-seq data interpretation. RNA transcriptomes from murine lens epithelium, lens fibers, erythrocytes, forebrain, liver, neurons, and pancreas were compared to establish the gene expression levels of the most abundant mRNAs versus median gene expression across other differentiated cells. RESULTS: Normalized RNA expression data from multiple tissues show that crystallins rank among the most highly expressed genes in mammalian cells. These findings correlate with the extremely high abundance of pol II all across the crystallin loci, including crystallin genes clustered on chromosomes 1 and 5, as well as within regions of "open" chromatin, as identified by FAIRE-seq. The expression levels of mRNAs encoding DNA-binding transcription factors (e.g., Foxe3, Hsf4, Maf, Pax6, Prox1, Sox1, and Tfap2a) revealed that their transcripts form "clusters" of abundant mRNAs in either lens fibers or lens epithelium. The expression of three autophagy regulatory mRNAs, encoding Tfeb, FoxO1, and Hif1α, was found within a group of lens preferentially expressed transcription factors compared to the E12.5 forebrain. CONCLUSIONS: This study reveals novel features of lens chromatin, including the remarkably high abundance of pol II at the crystallin loci that exhibit features of "open" chromatin. Hsf4 ranks among the most abundant fiber cell-preferred DNA-binding transcription factors. Notable transcripts, including Atf4, Ctcf, E2F4, Hey1, Hmgb1, Mycn, RXRß, Smad4, Sp1, and Taf1 (transcription factors) and Ctsd, Gabarapl1, and Park7 (autophagy regulators) have been identified with high levels of expression in lens fibers, which suggests specific roles in lens fiber cell terminal differentiation.

Comparison of REST cistromes across human cell types reveals common and context-specific functions.

Recent studies have shown that the transcriptional functions of REST are much broader than repressing neuronal genes in non-neuronal systems. Whether REST occupies similar chromatin regions in different cell types and how it interacts with other transcriptional regulators to execute its functions in a context-dependent manner has not been adequately investigated. We have applied ChIP-seq analysis to identify the REST cistrome in human CD4+ T cells and compared it with published data from 15 other cell types. We found that REST cistromes were distinct among cell types, with REST binding to several tumor suppressors specifically in cancer cells, whereas 7% of the REST peaks in non-neuronal cells were ubiquitously called and <25% were identified for ≥ 5 cell types. Nevertheless, using a quantitative metric directly comparing raw ChIP-seq signals, we found the majority (∼80%) was shared by ≥ 2 cell types. Integration with RNA-seq data showed that REST binding was generally correlated with low gene expression. Close examination revealed that multiple contexts were correlated with reduced expression of REST targets, e.g., the presence of a cognate RE1 motif and cellular specificity of REST binding. These contexts were shown to play a role in differential corepressor recruitment. Furthermore, transcriptional outcome was highly influenced by REST cofactors, e.g., SIN3 and EZH2 co-occupancy marked higher and lower expression of REST targets, respectively. Unexpectedly, the REST cistrome in differentiated neurons exhibited unique features not observed in non-neuronal cells, e.g., the lack of RE1 motifs and an association with active gene expression. Finally, our analysis demonstrated how REST could differentially regulate a transcription network constituted of miRNAs, REST complex and neuronal factors. Overall, our findings of contexts playing critical roles in REST occupancy and regulatory outcome provide insights into the molecular interactions underlying REST's diverse functions, and point to novel roles of REST in differentiated neurons.

Zebrafish and cellular models of SELENON-Related Myopathy exhibit novel embryonic and metabolic phenotypes.

SELENON-Related Myopathy (SELENON-RM) is a rare congenital myopathy caused by mutations of the SELENON gene characterized by axial muscle weakness and progressive respiratory insufficiency. Muscle histopathology commonly includes multiminicores or a dystrophic pattern but is often non-specific. The SELENON gene encodes selenoprotein N (SelN), a selenocysteine-containing redox enzyme located in the endo/sarcoplasmic reticulum membrane where it colocalizes with mitochondria-associated membranes. However, the molecular mechanism(s) by which SelN deficiency causes SELENON-RM are undetermined. A hurdle is the lack of cellular and animal models that show assayable phenotypes. Here we report deep-phenotyping of SelN-deficient zebrafish and muscle cells. SelN-deficient zebrafish exhibit changes in embryonic muscle function and swimming activity in larvae. Analysis of single cell RNAseq data in a zebrafish embryo-atlas revealed coexpression between selenon and genes involved in glutathione redox pathway. SelN-deficient zebrafish and mouse myoblasts exhibit changes in glutathione and redox homeostasis, suggesting a direct relationship with SelN function. We report changes in metabolic function abnormalities in SelN-null myotubes when compared to WT. These results suggest that SelN has functional roles during zebrafish early development and myoblast metabolism.

Cell-type-specific effects of autism-associated chromosome 15q11.2-13.1 duplications in human brain.

Recurrent copy number variation represents one of the most well-established genetic drivers in neurodevelopmental disorders, including autism spectrum disorder (ASD). Duplication of 15q11.2-13.1 (dup15q) is a well-described neurodevelopmental syndrome that increases the risk of ASD by over 40-fold. However, the effects of this duplication on gene expression and chromatin accessibility in specific cell types in the human brain remain unknown. To identify the cell-type-specific transcriptional and epigenetic effects of dup15q in the human frontal cortex we conducted single-nucleus RNA-sequencing and multi-omic sequencing on dup15q cases (n=6) as well as non-dup15q ASD (n=7) and neurotypical controls (n=7). Cell-type-specific differential expression analysis identified significantly regulated genes, critical biological pathways, and differentially accessible genomic regions. Although there was overall increased gene expression across the duplicated genomic region, cellular identity represented an important factor mediating gene expression changes. Neuronal subtypes, showed greater upregulation of gene expression across a critical region within the duplication as compared to other cell types. Genes within the duplicated region that had high baseline expression in control individuals showed only modest changes in dup15q, regardless of cell type. Of note, dup15q and ASD had largely distinct signatures of chromatin accessibility, but shared the majority of transcriptional regulatory motifs, suggesting convergent biological pathways. However, the transcriptional binding factor motifs implicated in each condition implicated distinct biological mechanisms; neuronal JUN/FOS networks in ASD vs. an inflammatory transcriptional network in dup15q microglia. This work provides a cell-type-specific analysis of how dup15q changes gene expression and chromatin accessibility in the human brain and finds evidence of marked cell-type-specific effects of this genetic driver. These findings have implications for guiding therapeutic development in dup15q syndrome, as well as understanding the functional effects CNVs more broadly in neurodevelopmental disorders.

Inflammation drives pathogenesis of early intestinal failure-associated liver disease.

Patients with intestinal failure who receive long-term parenteral nutrition (PN) often develop intestinal failure-associated liver disease (IFALD). Although there are identified risk factors, the early pathogenesis is poorly understood and treatment options are limited. Here, we perform a transcriptomic analysis of liver tissue in a large animal IFALD model to generate mechanistic insights and identify therapeutic targets. Preterm Yorkshire piglets were provided PN or bottle-fed with sow-milk replacer for 14 days. Compared to bottle-fed controls, piglets receiving PN developed biochemical cholestasis by day of life 15 (total bilirubin 0.2 vs. 2.9 mg/dL, P = 0.01). RNA-Seq of liver tissue was performed. Ingenuity Pathway Analysis identified 747 differentially expressed genes (343 upregulated and 404 downregulated) with an adjusted P < 0.05 and a fold-change of > |1|. Enriched canonical pathways were identified, demonstrating broad activation of inflammatory pathways and inhibition of cell cycle progression. Potential therapeutics including infliximab, glucocorticoids, statins, and obeticholic acid were identified as predicted upstream master regulators that may reverse the PN-induced gene dysregulation. The early driver of IFALD in neonates may be inflammation with an immature liver; identified therapeutics that target the inflammatory response in the liver should be investigated as potential treatments.

Genome-wide analysis of spina bifida risk variants in a case-control study from Bangladesh.

BACKGROUND: Human studies of genetic risk factors for neural tube defects, severe birth defects associated with long-term health consequences in surviving children, have predominantly been restricted to a subset of candidate genes in specific biological pathways including folate metabolism. METHODS: In this study, we investigated the association of genetic variants spanning the genome with risk of spina bifida (i.e., myelomeningocele and meningocele) in a subset of families enrolled from December 2016 through December 2022 in a case-control study in Bangladesh, a population often underrepresented in genetic studies. Saliva DNA samples were analyzed using the Illumina Global Screening Array. We performed genetic association analyses to compare allele frequencies between 112 case and 121 control children, 272 mothers, and 128 trios. RESULTS: In the transmission disequilibrium test analyses with trios only, we identified three novel exonic spina bifida risk loci, including rs140199800 (SULT1C2, p = 1.9 × 10-7), rs45580033 (ASB2, p = 4.2 × 10-10), and rs75426652 (LHPP, p = 7.2 × 10-14), after adjusting for multiple hypothesis testing. Association analyses comparing cases and controls, as well as models that included their mothers, did not identify genome-wide significant variants. CONCLUSIONS: This study identified three novel single nucleotide polymorphisms involved in biological pathways not previously associated with neural tube defects. The study warrants replication in larger groups to validate findings and to inform targeted prevention strategies.