ABSTRACT
Glioblastoma multiforme (GBM) is one of the most lethal brain tumors, characterized by profound heterogeneity. While single-cell transcriptomic studies have revealed extensive intra-tumor heterogeneity, shed light on intra-tumor diversity, spatial intricacies remain largely unexplored. Leveraging clinical GBM specimens, this study employs spatial transcriptomics technology to delve into gene expression heterogeneity. Our investigation unveils a significant enrichment of tissue stem cell signature in regions bordering necrosis and the peritumoral area, positively correlated with the mesenchymal subtype signature. Moreover, upregulated genes in these regions are linked with extracellular matrix (ECM)-receptor interaction, proteoglycans, as well as vascular endothelial growth factor (VEGF) and angiopoietin-Tie (ANGPT) signaling pathways. In contrast, signatures related to glycogen metabolism and oxidative phosphorylation show no relevance to pathological zoning, whereas creatine metabolism signature is notably exclusive to vascular-enriched areas. These spatial profiles not only offer valuable references but also pave the way for future in-depth functional and mechanistic investigations into GBM progression.
ABSTRACT
Genome-wide association studies of brain imaging phenotypes are mainly performed in European populations, but other populations are severely under-represented. Here, we conducted Chinese-alone and cross-ancestry genome-wide association studies of 3,414 brain imaging phenotypes in 7,058 Chinese Han and 33,224 white British participants. We identified 38 new associations in Chinese-alone analyses and 486 additional new associations in cross-ancestry meta-analyses at P < 1.46 × 10-11 for discovery and P < 0.05 for replication. We pooled significant autosomal associations identified by single- or cross-ancestry analyses into 6,443 independent associations, which showed uneven distribution in the genome and the phenotype subgroups. We further divided them into 44 associations with different effect sizes and 3,557 associations with similar effect sizes between ancestries. Loci of these associations were shared with 15 brain-related non-imaging traits including cognition and neuropsychiatric disorders. Our results provide a valuable catalog of genetic associations for brain imaging phenotypes in more diverse populations.
Subject(s)
Brain , East Asian People , Neuroimaging , White People , Adult , Female , Humans , Male , Asian People/genetics , Brain/diagnostic imaging , Genome-Wide Association Study , Magnetic Resonance Imaging , Phenotype , Polymorphism, Single Nucleotide , White People/genetics , East Asian People/genetics , United Kingdom , ChinaABSTRACT
Background: Sarcoma is a highly heterogeneous malignancy with a poor prognosis. Although chemotherapy and targeted therapy have improved the prognosis to some extent, the efficacy remains unsatisfactory in some patients. The efficacy and safety of immunotherapy in sarcoma need further evaluation. Methods: We conducted a two-center study of sarcoma patients receiving PD-1 immunotherapy at Tianjin Medical University Cancer Institute and Hospital and Henan Provincial Cancer Hospital. The treatment regimens included PD-1 inhibitor monotherapy and combination therapy based on PD-1 inhibitors. The observed primary endpoints were median progression-free survival (mPFS) and median overall survival (mOS). Survival curves were compared using the Kaplan-Meier method. Results: A total of 43 patients were included from the two centers. The median follow-up time for all patients was 13 months (range, 1-48 months). In the group of 37 patients with advanced or unresectable sarcoma, the mPFS was 6 months (95%CI: 5-12 months), and the mOS was 16 months (95%CI: 10-28 months). The ORR was 10.8% (4/37), and the DCR was 18.9% (7/37). Subgroup analysis showed no significant differences in mPFS (p=0.11) and mOS (p=0.88) between patients with PD-L1 negative/positive expression. There were also no significant differences in mPFS (p=0.13) or mOS (p=0.72) between PD-1 inhibitor monotherapy and combination therapy. Additionally, there were no significant differences in mPFS (p=0.52) or mOS (p=0.49) between osteogenic sarcoma and soft tissue sarcoma. Furthermore, the results showed no significant differences in mPFS (p=0.66) or mOS (p=0.96) between PD-1 inhibitors combined with targeted therapy and PD-1 inhibitors combined with AI chemotherapy. Among the 6 patients receiving adjuvant therapy after surgery, the mPFS was 15 months (95%CI: 6-NA months), and the mOS was not reached. In terms of safety, most adverse events were mild (grade 1-2) and manageable. The most severe grade 4 adverse events were bone marrow suppression, which occurred in 4 patients but resolved after treatment. There was also one case of a grade 4 adverse event related to hypertension. Conclusion: Immunotherapy is an effective treatment modality for sarcoma with manageable safety. Further inclusion of more patients or prospective clinical trials is needed to validate these findings.
Subject(s)
Bone Neoplasms , Sarcoma , Soft Tissue Neoplasms , Humans , Immune Checkpoint Inhibitors/adverse effects , Prospective Studies , Sarcoma/drug therapy , Immunotherapy/adverse effectsABSTRACT
Disruption of circadian rhythm during pregnancy produces adverse health outcomes in offspring; however, the role of maternal circadian rhythms in the immune system of infants and their susceptibility to inflammation remains poorly understood. Here we show that disruption of circadian rhythms in pregnant mice profoundly aggravates the severity of neonatal inflammatory disorders in both male and female offspring, such as necrotizing enterocolitis and sepsis. The diminished maternal production of docosahexaenoic acid (DHA) and the impaired immunosuppressive function of neonatal myeloid-derived suppressor cells (MDSCs) contribute to this phenomenon. Mechanistically, DHA enhances the immunosuppressive function of MDSCs via PPARγ-mediated mitochondrial oxidative phosphorylation. Transfer of MDSCs or perinatal supplementation of DHA relieves neonatal inflammation induced by maternal rhythm disruption. These observations collectively demonstrate a previously unrecognized role of maternal circadian rhythms in the control of neonatal inflammation via metabolic reprograming of myeloid cells.
Subject(s)
Animals, Newborn , Circadian Rhythm , Inflammation , Myeloid Cells , Animals , Female , Mice , Inflammation/metabolism , Pregnancy , Myeloid Cells/metabolism , Male , Docosahexaenoic Acids/metabolism , Docosahexaenoic Acids/pharmacology , Myeloid-Derived Suppressor Cells/metabolism , Mice, Inbred C57BLABSTRACT
The co-occurrence and familial clustering of neurodevelopmental disorders and immune disorders suggest shared genetic risk factors. Based on genome-wide association summary statistics from five neurodevelopmental disorders and four immune disorders, we conducted genome-wide, local genetic correlation and polygenic overlap analysis. We further performed a cross-trait GWAS meta-analysis. Pleotropic loci shared between the two categories of diseases were mapped to candidate genes using multiple algorithms and approaches. Significant genetic correlations were observed between neurodevelopmental disorders and immune disorders, including both positive and negative correlations. Neurodevelopmental disorders exhibited higher polygenicity compared to immune disorders. Around 50%-90% of genetic variants of the immune disorders were shared with neurodevelopmental disorders. The cross-trait meta-analysis revealed 154 genome-wide significant loci, including 8 novel pleiotropic loci. Significant associations were observed for 30 loci with both types of diseases. Pathway analysis on the candidate genes at these loci revealed common pathways shared by the two types of diseases, including neural signaling, inflammatory response, and PI3K-Akt signaling pathway. In addition, 26 of the 30 lead SNPs were associated with blood cell traits. Neurodevelopmental disorders exhibit complex polygenic architecture, with a subset of individuals being at a heightened genetic risk for both neurodevelopmental and immune disorders. The identification of pleiotropic loci has important implications for exploring opportunities for drug repurposing, enabling more accurate patient stratification, and advancing genomics-informed precision in the medical field of neurodevelopmental disorders.
Subject(s)
Genetic Predisposition to Disease , Genome-Wide Association Study , Immune System Diseases , Multifactorial Inheritance , Neurodevelopmental Disorders , Polymorphism, Single Nucleotide , Humans , Neurodevelopmental Disorders/genetics , Immune System Diseases/genetics , Genetic Predisposition to Disease/genetics , Polymorphism, Single Nucleotide/genetics , Multifactorial Inheritance/geneticsABSTRACT
Mapping of expression quantitative trait loci (eQTLs) and other molecular QTLs can help characterize the modes of action of disease-associated genetic variants. However, current eQTL databases present data from bulk RNA-seq approaches, which cannot shed light on the cell type- and environment-specific regulation of disease-associated genetic variants. Here, we introduce our Single-cell eQTL Interactive Database which collects single-cell eQTL (sc-eQTL) datasets and provides online visualization of sc-eQTLs across different cell types in a user-friendly manner. Although sc-eQTL mapping is still in its early stage, our database curates the most comprehensive summary statistics of sc-eQTLs published to date. sc-eQTL studies have revolutionized our understanding of gene regulation in specific cellular contexts, and we anticipate that our database will further accelerate the research of functional genomics. Database URL: http://www.sqraolab.com/scqtl.
Subject(s)
Gene Expression Regulation , Quantitative Trait Loci , Humans , Quantitative Trait Loci/genetics , RNA-Seq , Genomics , Genome-Wide Association Study , Polymorphism, Single NucleotideABSTRACT
BACKGROUND: CLEC16A intron 19 has been identified as a candidate locus for common variable immunodeficiency (CVID). OBJECTIVES: This study sought to elucidate the molecular mechanism by which variants at the CLEC16A intronic locus may contribute to the pathogenesis of CVID. METHODS: The investigators performed fine-mapping of the CLEC16A locus in a CVID cohort, then deleted the candidate functional SNP in T-cell lines by the CRISPR-Cas9 technique and conducted RNA-sequencing to identify target gene(s). The interactions between the CLEC16A locus and its target genes were identified using circular chromosome conformation capture. The transcription factor complexes mediating the chromatin interactions were determined by proteomic approach. The molecular pathways regulated by the CLEC16A locus were examined by RNA-sequencing and reverse phase protein array. RESULTS: This study showed that the CLEC16A locus is an enhancer regulating expression of multiple target genes including a distant gene ATF7IP2 through chromatin interactions. Distinct transcription factor complexes mediate the chromatin interactions in an allele-specific manner. Disruption of the CLEC16A locus affects the AKT signaling pathway, as well as the molecular response of CD4+ T cells to immune stimulation. CONCLUSIONS: Through multiomics and targeted experimental approaches, this study elucidated the underlying target genes and signaling pathways involved in the genetic association of CLEC16A with CVID, and highlighted plausible molecular targets for developing novel therapeutics.
Subject(s)
Common Variable Immunodeficiency , Introns , Lectins, C-Type , Monosaccharide Transport Proteins , Humans , Lectins, C-Type/genetics , Introns/genetics , Monosaccharide Transport Proteins/genetics , Common Variable Immunodeficiency/genetics , Common Variable Immunodeficiency/immunology , Polymorphism, Single Nucleotide , Gene Expression Regulation , Female , Male , Signal Transduction/genetics , CD4-Positive T-Lymphocytes/immunology , AdultABSTRACT
BACKGROUND: The visual outcome of open globe injury (OGI)-no light perception (NLP) eyes is unpredictable traditionally. This study aimed to develop a model to predict the visual outcomes of vitrectomy surgery in OGI-NLP eyes using a machine learning algorithm and to provide an interpretable system for the prediction results. METHODS: Clinical data of 459 OGI-NLP eyes were retrospectively collected from 19 medical centres across China to establish a training data set for developing a model, called 'VisionGo', which can predict the visual outcome of the patients involved and compare with the Ocular Trauma Score (OTS). Another 72 cases were retrospectively collected and used for human-machine comparison, and an additional 27 cases were prospectively collected for real-world validation of the model. The SHapley Additive exPlanations method was applied to analyse feature contribution to the model. An online platform was built for real-world application. RESULTS: The area under the receiver operating characteristic curve (AUC) of VisionGo was 0.75 and 0.90 in previtrectomy and intravitrectomy application scenarios, which was much higher than the OTS (AUC=0.49). VisionGo showed better performance than ophthalmologists in both previtrectomy and intravitrectomy application scenarios (AUC=0.73 vs 0.57 and 0.87 vs 0.64). In real-world validation, VisionGo achieved an AUC of 0.60 and 0.91 in previtrectomy and intravitrectomy application scenarios. Feature contribution analysis indicated that wound length-related indicators, vitreous status and retina-related indicators contributed highly to visual outcomes. CONCLUSIONS: VisionGo has achieved an accurate and reliable prediction in visual outcome after vitrectomy for OGI-NLP eyes.
Subject(s)
Eye Injuries, Penetrating , Eye Injuries , Humans , Retrospective Studies , Visual Acuity , Retina , Vitrectomy , Prognosis , Eye Injuries, Penetrating/diagnosis , Eye Injuries, Penetrating/surgeryABSTRACT
Telomere length (TL) shortening is a pivotal indicator of biological aging and is associated with many human diseases. The genetic determinates of human TL have been widely investigated, however, most existing studies were conducted based on adult tissues which are heavily influenced by lifetime exposure. Based on the analyses of terminal restriction fragment (TRF) length of telomere, individual genotypes, and gene expressions on 166 healthy placental tissues, we systematically interrogate TL-modulated genes and their potential functions. We discover that the TL in the placenta is comparatively longer than in other adult tissues, but exhibiting an intra-tissue homogeneity. Trans-ancestral TL genome-wide association studies (GWASs) on 644,553 individuals identify 20 newly discovered genetic associations and provide increased polygenic determination of human TL. Next, we integrate the powerful TL GWAS with placental expression quantitative trait locus (eQTL) mapping to prioritize 23 likely causal genes, among which 4 are functionally validated, including MMUT, RRM1, KIAA1429, and YWHAZ. Finally, modeling transcriptomic signatures and TRF-based TL improve the prediction performance of human TL. This study deepens our understanding of causal genes and transcriptomic determinants of human TL, promoting the mechanistic research on fine-grained TL regulation.
Subject(s)
Genome-Wide Association Study , Placenta , Adult , Humans , Female , Pregnancy , Placenta/metabolism , Telomere Shortening , Telomere/genetics , Gene Expression ProfilingABSTRACT
Microbiota have an important function in shaping and priming neonatal immunity, although the cellular and molecular mechanisms underlying these effects remain obscure. Here we report that prenatal antibiotic exposure causes significant elevation of group 2 innate lymphoid cells (ILC2s) in neonatal lungs, in both cell numbers and functionality. Downregulation of type 1 interferon signaling in ILC2s due to diminished production of microbiota-derived butyrate represents the underlying mechanism. Mice lacking butyrate receptor GPR41 (Gpr41-/-) or type 1 interferon receptor IFNAR1 (Ifnar1-/-) recapitulate the phenotype of neonatal ILC2s upon maternal antibiotic exposure. Furthermore, prenatal antibiotic exposure induces epigenetic changes in ILC2s and has a long-lasting deteriorative effect on allergic airway inflammation in adult offspring. Prenatal supplementation of butyrate ameliorates airway inflammation in adult mice born to antibiotic-exposed dams. These observations demonstrate an essential role for the microbiota in the control of type 2 innate immunity at the neonatal stage, which suggests a therapeutic window for treating asthma in early life.
Subject(s)
Anti-Bacterial Agents , Immunity, Innate , Interferon Type I , Lymphocytes , Animals , Mice , Butyrates , Cytokines , Down-Regulation , Inflammation , Lung , Anti-Bacterial Agents/adverse effects , Anti-Bacterial Agents/pharmacology , Maternal Exposure , Interferon Type I/drug effects , Interferon Type I/metabolismABSTRACT
BACKGROUND: The high mutation rate throughout the entire melanoma genome presents a major challenge in stratifying true driver events from the background mutations. Numerous recurrent non-coding alterations, such as those in enhancers, can shape tumor evolution, thereby emphasizing the importance in systematically deciphering enhancer disruptions in melanoma. RESULTS: Here, we leveraged 297 melanoma whole-genome sequencing samples to prioritize highly recurrent regions. By performing a genome-scale CRISPR interference (CRISPRi) screen on highly recurrent region-associated enhancers in melanoma cells, we identified 66 significant hits which could have tumor-suppressive roles. These functional enhancers show unique mutational patterns independent of classical significantly mutated genes in melanoma. Target gene analysis for the essential enhancers reveal many known and hidden mechanisms underlying melanoma growth. Utilizing extensive functional validation experiments, we demonstrate that a super enhancer element could modulate melanoma cell proliferation by targeting MEF2A, and another distal enhancer is able to sustain PTEN tumor-suppressive potential via long-range interactions. CONCLUSIONS: Our study establishes a catalogue of crucial enhancers and their target genes in melanoma growth and progression, and illuminates the identification of novel mechanisms of dysregulation for melanoma driver genes and new therapeutic targeting strategies.
Subject(s)
Enhancer Elements, Genetic , Melanoma , Humans , Melanoma/genetics , Melanoma/pathology , MutationABSTRACT
Unscheduled R-loops are a major source of replication stress and DNA damage. R-loop-induced replication defects are sensed and suppressed by ATR kinase, whereas it is not known whether R-loop itself is actively involved in ATR activation and, if so, how this is achieved. Here, we report that the nuclear form of RNA-editing enzyme ADAR1 promotes ATR activation and resolves genome-wide R-loops, a process that requires its double-stranded RNA-binding domains. Mechanistically, ADAR1 interacts with TOPBP1 and facilitates its loading on perturbed replication forks by enhancing the association of TOPBP1 with RAD9 of the 9-1-1 complex. When replication is inhibited, DNA-RNA hybrid competes with TOPBP1 for ADAR1 binding to promote the translocation of ADAR1 from damaged fork to accumulate at R-loop region. There, ADAR1 recruits RNA helicases DHX9 and DDX21 to unwind R-loops, simultaneously allowing TOPBP1 to stimulate ATR more efficiently. Collectively, we propose that the tempo-spatially regulated assembly of ADAR1-nucleated protein complexes link R-loop clearance and ATR activation, while R-loops crosstalk with blocked replication forks by transposing ADAR1 to finetune ATR activity and safeguard the genome.
Subject(s)
DNA-Binding Proteins , R-Loop Structures , Ataxia Telangiectasia Mutated Proteins/genetics , Ataxia Telangiectasia Mutated Proteins/metabolism , Cell Cycle Proteins/metabolism , DNA Replication , DNA-Binding Proteins/genetics , RNA/genetics , Humans , Animals , MiceABSTRACT
CCCTC-binding factor (CTCF) is a transcription regulator with a complex role in gene regulation. The recognition and effects of CTCF on DNA sequences, chromosome barriers, and enhancer blocking are not well understood. Existing computational tools struggle to assess the regulatory potential of CTCF-binding sites and their impact on chromatin loop formation. Here we have developed a deep-learning model, DeepAnchor, to accurately characterize CTCF binding using high-resolution genomic/epigenomic features. This has revealed distinct chromatin and sequence patterns for CTCF-mediated insulation and looping. An optimized implementation of a previous loop model based on DeepAnchor score excels in predicting CTCF-anchored loops. We have established a compendium of CTCF-anchored loops across 52 human tissue/cell types, and this suggests that genomic disruption of these loops could be a general mechanism of disease pathogenesis. These computational models and resources can help investigate how CTCF-mediated cis-regulatory elements shape context-specific gene regulation in cell development and disease progression.
ABSTRACT
The hippocampus is critical for memory and cognition and neuropsychiatric disorders, and its subfields differ in architecture and function. Genome-wide association studies on hippocampal and subfield volumes are mainly conducted in European populations; however, other ancestral populations are under-represented. Here we conduct cross-ancestry genome-wide association meta-analyses in 65,791 individuals for hippocampal volume and 38,977 for subfield volumes, including 7,009 individuals of East Asian ancestry. We identify 339 variant-trait associations at P < 1.13 × 10-9 for 44 hippocampal traits, including 23 new associations. Common genetic variants have similar effects on hippocampal traits across ancestries, although ancestry-specific associations exist. Cross-ancestry analysis improves the fine-mapping precision and the prediction performance of polygenic scores in under-represented populations. These genetic variants are enriched for Wnt signaling and neuron differentiation and affect cognition, emotion and neuropsychiatric disorders. These findings may provide insight into the genetic architectures of hippocampal and subfield volumes.
Subject(s)
Genome-Wide Association Study , Magnetic Resonance Imaging , Humans , Hippocampus/diagnostic imaging , CognitionABSTRACT
A previous genome-wide association study (GWAS) revealed an association of the noncoding SNP rs1663689 with susceptibility to lung cancer in the Chinese population. However, the underlying mechanism is unknown. In this study, using allele-specific 4C-seq in heterozygous lung cancer cells combined with epigenetic information from CRISPR/Cas9-edited cell lines, we show that the rs1663689 C/C variant represses the expression of ADGRG6, a gene located on a separate chromosome, through an interchromosomal interaction of the rs1663689 bearing region with the ADGRG6 promoter. This reduces downstream cAMP-PKA signaling and subsequently tumor growth both in vitro and in xenograft models. Using patient-derived organoids, we show that rs1663689 T/T-but not C/C-bearing lung tumors are sensitive to the PKA inhibitor H89, potentially informing therapeutic strategies. Our study identifies a genetic variant-mediated interchromosomal interaction underlying ADGRG6 regulation and suggests that targeting the cAMP-PKA signaling pathway may be beneficial in lung cancer patients bearing the homozygous risk genotype at rs1663689.
Subject(s)
Genome-Wide Association Study , Lung Neoplasms , Humans , Lung Neoplasms/genetics , Lung , Receptors, G-Protein-Coupled/genetics , Gene Expression RegulationABSTRACT
To date, genome-wide association studies (GWAS) have revealed over 200 genetic risk loci associated with prostate cancer; yet, true disease-causing variants remain elusive. Identification of causal variants and their targets from association signals is complicated by high linkage disequilibrium and limited availability of functional genomics data for specific tissue/cell types. Here, we integrated statistical fine-mapping and functional annotation from prostate-specific epigenomic profiles, 3D genome features, and quantitative trait loci data to distinguish causal variants from associations and identify target genes. Our fine-mapping analysis yielded 3,395 likely causal variants, and multiscale functional annotation linked them to 487 target genes. We prioritized rs10486567 as a genome-wide top-ranked SNP and predicted HOTTIP as its target. Deletion of the rs10486567-associated enhancer in prostate cancer cells decreased their capacity for invasive migration. HOTTIP overexpression in enhancer-KO cell lines rescued defective invasive migration. Furthermore, we found that rs10486567 regulates HOTTIP through allele-specific long-range chromatin interaction.
ABSTRACT
Whole -genome sequencing projects of millions of subjects contain enormous genotypes, entailing a huge memory burden and time for computation. Here, we present GBC, a toolkit for rapidly compressing large-scale genotypes into highly addressable byte-encoding blocks under an optimized parallel framework. We demonstrate that GBC is up to 1000 times faster than state-of-the-art methods to access and manage compressed large-scale genotypes while maintaining a competitive compression ratio. We also showed that conventional analysis would be substantially sped up if built on GBC to access genotypes of a large population. GBC's data structure and algorithms are valuable for accelerating large-scale genomic research.
Subject(s)
Data Compression , Software , Humans , Algorithms , Genotype , Data Compression/methods , Genomics/methodsABSTRACT
Age-related macular degeneration (AMD), cataract, and glaucoma are leading causes of blindness worldwide. Previous genome-wide association studies (GWASs) have revealed a variety of susceptible loci associated with age-related ocular disorders, yet the genetic pleiotropy and causal genes across these diseases remain poorly understood. By leveraging large-scale genetic and observational data from ocular disease GWASs and UK Biobank (UKBB), we found significant pairwise genetic correlations and consistent epidemiological associations among these ocular disorders. Cross-disease meta-analysis uncovered seven pleiotropic loci, three of which were replicated in an additional cohort. Integration of variants in pleiotropic loci and multiple single-cell omics data identified that Müller cells and astrocytes were likely trait-related cell types underlying ocular comorbidity. In addition, we comprehensively integrated eye-specific gene expression quantitative loci (eQTLs), epigenomic profiling, and 3D genome data to prioritize causal pleiotropic genes. We found that pleiotropic genes were essential in nerve development and eye pigmentation, and targetable by aflibercept and pilocarpine for the treatment of AMD and glaucoma. These findings will not only facilitate the mechanistic research of ocular comorbidities but also benefit the therapeutic optimization of age-related ocular diseases.
Subject(s)
Glaucoma , Macular Degeneration , Humans , Genetic Pleiotropy , Genome-Wide Association Study , Genetic Predisposition to Disease , Macular Degeneration/genetics , Glaucoma/genetics , Polymorphism, Single NucleotideABSTRACT
Genetic sharing is extensively observed for autoimmune diseases, but the causal variants and their underlying molecular mechanisms remain largely unknown. Through systematic investigation of autoimmune disease pleiotropic loci, we found most of these shared genetic effects are transmitted from regulatory code. We used an evidence-based strategy to functionally prioritize causal pleiotropic variants and identify their target genes. A top-ranked pleiotropic variant, rs4728142, yielded many lines of evidence as being causal. Mechanistically, the rs4728142-containing region interacts with the IRF5 alternative promoter in an allele-specific manner and orchestrates its upstream enhancer to regulate IRF5 alternative promoter usage through chromatin looping. A putative structural regulator, ZBTB3, mediates the allele-specific loop to promote IRF5-short transcript expression at the rs4728142 risk allele, resulting in IRF5 overactivation and M1 macrophage polarization. Together, our findings establish a causal mechanism between the regulatory variant and fine-scale molecular phenotype underlying the dysfunction of pleiotropic genes in human autoimmunity.
Subject(s)
Autoimmune Diseases , DNA-Binding Proteins , Interferon Regulatory Factors , Humans , Alleles , Autoimmunity , Chromatin , Polymorphism, Single Nucleotide , Promoter Regions, GeneticABSTRACT
Polycomb group (PcG) proteins are critical chromatin regulators for cell fate control. The mono-ubiquitylation on histone H2AK119 (H2AK119ub1) is one of the well-recognized mechanisms for Polycomb repressive complex 1 (PRC1)-mediated transcription repression. Unexpectedly, the specific H2AK119 deubiquitylation complex composed by additional sex comb-like proteins and BAP1 has also been genetically characterized as Polycomb repressive deubiquitnase (PR-DUB) for unclear reasons. However, it remains a mystery whether and how PR-DUB deficiency affects chromatin states and cell fates through impaired PcG silencing. Here through a careful epigenomic analysis, we demonstrate that a bulk of H2AK119ub1 is diffusely distributed away from promoter regions and their enrichment is positively correlated with PRC1 occupancy. Upon deletion of Asxl2 in mouse embryonic stem cells (ESCs), a pervasive gain of H2AK119ub1 is coincident with increased PRC1 sampling at chromatin. Accordingly, PRC1 is significantly lost from a subset of highly occupied promoters, leading to impaired silencing of associated genes before and after lineage differentiation of Asxl2-null ESCs. Therefore, our study highlights the importance of genome-wide H2AK119ub1 restriction by PR-DUB in safeguarding robust PRC1 deposition and its roles in developmental regulation.