Chromatin Accessibility Profiles Research Articles

Analysis of potential biomarkers for diabetic kidney disease based on single-cell RNA-sequencing integrated with a single-cell sequencing assay for transposase-accessible chromatin.

Diabetic kidney disease (DKD) is a renal microvascular disease caused by hyperglycemia that involves metabolic remodeling, oxidative stress, inflammation, and other factors. The mechanism is complex and not fully unraveled. We performed an integrated single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) and single-cell RNA-sequencing (scRNA-seq) analyses of kidneys from db/db and db/m mice to identify differential open chromatin regions and gene expression, particularly in genes related to proximal tubular reabsorption and secretion. We identified 9,776 differentially expressed genes (DEGs) and 884 cell type-specific transcription factors (TFs) across 15 cell types. Glucose and lipid transporters, and TFs related to the circadian rhythm in the proximal tubules had significantly higher expression in db/db mice than in db/m mice (P<0.01). Crosstalk between podocytes and tubular cells in the proximal tubules was enhanced, and renal inflammation, oxidative stress, and fibrosis pathways were activated in db/db mice. Western blotting and immunohistochemical staining results showed that Wfdc2 expression in the urine and kidneys of DKD patients was higher than that in non-diabetic kidney disease (NDKD) controls. The revealed landscape of chromatin accessibility and transcriptional profiles in db/db mice provide insights into the pathological mechanism of DKD.

Remodeling of Chromatin Accessibility Regulates the Radiological Responses of NSCLC A549 Cells to High-LET Carbon Ions.

Carbon-ion radiation therapy (CIRT) may offer remarkable advantages in cancer treatment with its unique physical and biological characteristics. However, the underlying epigenetic regulatory mechanisms of cancer response to CIRT remain to be identified. In this study, we showed consistent but different degrees of biological effects induced in NSCLC A549 cells by carbon ions of different LET. The genome-wide chromatin accessibility and transcriptional profiles of carbon ion-treated A549 cells were performed using transposase-accessible chromatin sequencing (ATAC-seq) and RNA-seq, respectively, and further gene regulatory network analysis was performed by integrating the two sets of genomic data. Alterations in chromatin accessibility by carbon ions of different LET predominantly occurred in intron, distal intergenic and promoter regions of differential chromatin accessibility regions. The transcriptional changes were mainly regulated by proximal chromatin accessibility. Notably, CCCTC-binding factor (CTCF) was identified as a key transcription factor in the cellular response to carbon ions. The target genes regulated by CTCF in response to carbon ions were found to be closely associated with the LET of carbon ions, particularly in the regulation of gene transcription within the DNA replication- and metabolism-related signaling pathways. This study provides a regulatory profile of genes involved in key signaling pathways and highlighted key regulatory elements in NSCLC A549 cells during CIRT, which expands our understanding of the epigenetic mechanisms of carbon ion-induced biological effects and reveals an important role for LET in the regulation of changes in chromatin accessibility, although further research is needed.

Integrative single-nucleus multi-omics analysis prioritizes candidate cis and trans regulatory networks and their target genes in Alzheimer’s disease brains

BackgroundThe genetic underpinnings of late-onset Alzheimer’s disease (LOAD) are yet to be fully elucidated. Although numerous LOAD-associated loci have been discovered, the causal variants and their target genes remain largely unknown. Since the brain is composed of heterogenous cell subtypes, it is imperative to study the brain on a cell subtype specific level to explore the biological processes underlying LOAD.MethodsHere, we present the largest parallel single-nucleus (sn) multi-omics study to simultaneously profile gene expression (snRNA-seq) and chromatin accessibility (snATAC-seq) to date, using nuclei from 12 normal and 12 LOAD brains. We identified cell subtype clusters based on gene expression and chromatin accessibility profiles and characterized cell subtype-specific LOAD-associated differentially expressed genes (DEGs), differentially accessible peaks (DAPs) and cis co-accessibility networks (CCANs).ResultsIntegrative analysis defined disease-relevant CCANs in multiple cell subtypes and discovered LOAD-associated cell subtype-specific candidate cis regulatory elements (cCREs), their candidate target genes, and trans-interacting transcription factors (TFs), some of which, including ELK1, JUN, and SMAD4 in excitatory neurons, were also LOAD-DEGs. Finally, we focused on a subset of cell subtype-specific CCANs that overlap known LOAD-GWAS regions and catalogued putative functional SNPs changing the affinities of TF motifs within LOAD-cCREs linked to LOAD-DEGs, including APOE and MYO1E in a specific subtype of microglia and BIN1 in a subpopulation of oligodendrocytes.ConclusionsTo our knowledge, this study represents the most comprehensive systematic interrogation to date of regulatory networks and the impact of genetic variants on gene dysregulation in LOAD at a cell subtype resolution. Our findings reveal crosstalk between epigenetic, genomic, and transcriptomic determinants of LOAD pathogenesis and define catalogues of candidate genes, cCREs, and variants involved in LOAD genetic etiology and the cell subtypes in which they act to exert their pathogenic effects. Overall, these results suggest that cell subtype-specific cis–trans interactions between regulatory elements and TFs, and the genes dysregulated by these networks contribute to the development of LOAD.

Hybridization led to a rewired pluripotency network in the allotetraploid Xenopus laevis.

After fertilization, maternally contributed factors to the egg initiate the transition to pluripotency to give rise to embryonic stem cells, in large part by activating de novo transcription from the embryonic genome. Diverse mechanisms coordinate this transition across animals, suggesting that pervasive regulatory remodeling has shaped the earliest stages of development. Here, we show that maternal homologs of mammalian pluripotency reprogramming factors OCT4 and SOX2 divergently activate the two subgenomes of Xenopus laevis, an allotetraploid that arose from hybridization of two diploid species ~18 million years ago. Although most genes have been retained as two homeologous copies, we find that a majority of them undergo asymmetric activation in the early embryo. Chromatin accessibility profiling and CUT&RUN for modified histones and transcription factor binding reveal extensive differences in predicted enhancer architecture between the subgenomes, which likely arose through genomic disruptions as a consequence of allotetraploidy. However, comparison with diploid X. tropicalis and zebrafish shows broad conservation of embryonic gene expression levels when divergent homeolog contributions are combined, implying strong selection to maintain dosage in the core vertebrate pluripotency transcriptional program, amid genomic instability following hybridization.

Chromatin accessibility profiling of targeted cell populations with laser capture microdissection coupled to ATAC-seq

Spatially resolved omics technologies reveal context-dependent cellular regulatory networks in tissues of interest. Beyond transcriptome analysis, information on epigenetic traits and chromatin accessibility can provide further insights on gene regulation in health and disease. Nevertheless, compared to the enormous advancements in spatial transcriptomics technologies, the field of spatial epigenomics is much younger and still underexplored. In this study, we report laser capture microdissection coupled to ATAC-seq (LCM-ATAC-seq) applied to fresh frozen samples for the spatial characterization of chromatin accessibility. We first demonstrate the efficient use of LCM coupled to in situ tagmentation and evaluate its performance as a function of cell number, microdissected areas, and tissue type. Further, we demonstrate its use for the targeted chromatin accessibility analysis of discrete contiguous or scattered cell populations in tissues via single-nuclei capture based on immunostaining for specific cellular markers.

Chromatin accessibility, not 5mC methylation covaries with partial dosage compensation in crows.

The evolution of genetic sex determination is often accompanied by degradation of the sex-limited chromosome. Male heterogametic systems have evolved convergent, epigenetic mechanisms restoring the resulting imbalance in gene dosage between diploid autosomes (AA) and the hemizygous sex chromosome (X). Female heterogametic systems (AAf Zf, AAm ZZm) tend to only show partial dosage compensation (0.5 < Zf:AAf < 1) and dosage balance (0.5<Zf:ZZm<1). The underlying mechanism remains largely elusive. Here, we quantified gene expression for a total of 15 male and female Eurasian crows (Corvus (corone) spp.) raised under common garden conditions. In addition, we characterized aspects of the regulatory epigenetic landscape quantifying chromatin accessibility (ATAC-seq) and 5mC methylation profiles. Partial dosage balance and compensation was due to female upregulation of Z-linked genes which covaried significantly with increased chromatin accessibility of the female Z chromosome. 5mC methylation was tissue and sex chromosome-specific, but unrelated to dosage. With the exception of the pseudo-autosomal region (PAR), female upregulation of gene expression was evenly spread across the Z chromosome without evidence for regional centers of epigenetic regulation, as has, for example, been suggested for the male hypermethylated region (MHM) in chicken. Our results suggest that partial dosage balance and compensation in female heterogametic systems are tightly linked to chromosome-wide, epigenetic control of the female Z chromosome mediated by differential chromatin accessibility.

Dissection of tumor-intrinsic and tumor-extrinsic features of MiT/TFE translocation renal cell carcinoma via single-cell RNA sequencing

Abstract Background Translocation renal cell carcinoma (tRCC) is a rare and aggressive subtype of renal cell carcinoma driven by oncogenic gene fusions involving MiT/TFE family transcription factors, most commonly TFE3. Currently, there are no molecularly tailored treatments for tRCC, and standard-of-care therapies utilized for other RCC subtypes are typically less effective in tRCC. Emerging data suggest that tRCC is molecularly distinct from more common RCC subtypes. However, an incomplete understanding of both tumor-intrinsic drivers and the tumor microenvironment (TME) features of tRCC presents barriers to developing effective therapeutics for this cancer. Methods We used single-nucleus RNA sequencing (snRNA) to profile ten samples of tRCC tumors and one adjacent normal tissue. Five of the samples also underwent multiome profiling with snATAC-seq while nine of the tumor samples underwent whole genome sequencing and RNA-sequencing. We preprocessed and obtained the raw fasta files using cellranger and extracted RNA counts and chromatin accessibility peaks using the Seurat/Signac algorithm. Scrublet was used to exclude data from droplets containing more than one cell by performing doublet detection and removal on gene-barcode matrices. Cells with fewer than 200 genes detected or more than 5% of counts attributed to mitochondrially-encoded transcripts were removed before across-sample integration. We also excluded genes detected in fewer than three cells across all samples. The combined cohort was analyzed using the MergeData function in Seurat, and cancer cells were selected based on known tRCC gene markers and inferred transcriptional copy number variations estimated via the inferCNV package and verified from WGS data. Non-tumor cell types were determined through manual annotation via known marker genes. We analyzed patterns of gene expression at the single-cell level using the Seurat V4 package and module score functions. Differential expression analysis comparing cells from different clusters or treatment exposure groups was performed using a two-sided Wilcoxon rank-sum test with Bonferroni FDR correction. Results Following quality control and integration, 71,124 single-cells were used for analysis (67,762 from tumor samples) – representing 43,214 tumor cells, 8,352 monocytes, 4,692 T cells, 4,614 endothelial/stromal cells, and normal kidney cells. The normal-adjacent sample represented multiple cell types from the normal kidney, including cells from distal connecting tubule, connecting duct, proximal tubule, thick ascending limb, podocytes, and endothelial cells. Analysis of chromatin accessibility profiles from snATAC-seq and snRNA-seq data of tumor cells was employed to identify unique cell states and a putative cell of origin for tRCCs. Tumor cell analysis also demonstrated intra- and inter-tumoral transcriptional variability in genetically similar tumor subclusters, each with distinct regulon activity. Immune subpopulation analyses revealed higher proportions of resting T cells in treatment-naive tRCC samples and higher proportions of progenitor-exhausted and terminal-exhausted populations in an immune checkpoint inhibitor (ICI) treated sample, indicating immune reprogramming in response to ICI in this mutationally quiet tumor. To further understand the role of, and factors permissive to, anti-tumor T cell responses in tRCC (which is characterized by a low tumor mutational burden) we explored the hypothesis that MiT/TFE fusions may constitute tumor neoantigens. We computationally identified fusion-associated neoantigens corresponding to these tRCC samples, as well as across a larger aggregate dataset of tRCC samples profiled by bulk RNA-Seq. We predicted binding affinities between fusion-associated neoantigens and HLA types and identified multiple peptide candidates with high HLA-binding affinity, which will be validated for their immunogenicity in vitro. Conclusions Overall, our results highlight tumor-intrinsic and tumor-extrinsic determinants of immunogenicity in tRCC and may guide the development of immunotherapeutic strategies with strong mechanistic rationale in tRCC. CDMRP DOD Funding: yes

Regulatory changes associated with the head to trunk developmental transition

BackgroundDevelopment of vertebrate embryos is characterized by early formation of the anterior tissues followed by the sequential extension of the axis at their posterior end to build the trunk and tail structures, first by the activity of the primitive streak and then of the tail bud. Embryological, molecular and genetic data indicate that head and trunk development are significantly different, suggesting that the transition into the trunk formation stage involves major changes in regulatory gene networks.ResultsWe explored those regulatory changes by generating differential interaction networks and chromatin accessibility profiles from the posterior epiblast region of mouse embryos at embryonic day (E)7.5 and E8.5. We observed changes in various cell processes, including several signaling pathways, ubiquitination machinery, ion dynamics and metabolic processes involving lipids that could contribute to the functional switch in the progenitor region of the embryo. We further explored the functional impact of changes observed in Wnt signaling associated processes, revealing a switch in the functional relevance of Wnt molecule palmitoleoylation, essential during gastrulation but becoming differentially required for the control of axial extension and progenitor differentiation processes during trunk formation. We also found substantial changes in chromatin accessibility at the two developmental stages, mostly mapping to intergenic regions and presenting differential footprinting profiles to several key transcription factors, indicating a significant switch in the regulatory elements controlling head or trunk development. Those chromatin changes are largely independent of retinoic acid, despite the key role of this factor in the transition to trunk development. We also tested the functional relevance of potential enhancers identified in the accessibility assays that reproduced the expression profiles of genes involved in the transition. Deletion of these regions by genome editing had limited effect on the expression of those genes, suggesting the existence of redundant enhancers that guarantee robust expression patterns.ConclusionsThis work provides a global view of the regulatory changes controlling the switch into the axial extension phase of vertebrate embryonic development. It also revealed mechanisms by which the cellular context influences the activity of regulatory factors, channeling them to implement one of several possible biological outputs.

The chromatin and single-cell transcriptional landscapes of CD4 T cells in inflammatory bowel disease link risk loci with a proinflammatory Th17 cell population.

The imbalance between Th17 and regulatory T cells in inflammatory bowel diseases (IBD) promotes intestinal epithelial cell damage. In this scenario, T helper cell lineage commitment is accompanied by dynamic changes to the chromatin that facilitate or repress gene expression. Here, we characterized the chromatin landscape and heterogeneity of intestinal and peripheral CD4 T cellsfrom IBD patients using in house ATAC-Seq and single cell RNA-Seq libraries. We show that chromatin accessibility profiles of CD4 T cells from inflamed intestinal biopsies relate to genes associated with a network of inflammatory processes. After integrating the chromatin profiles of tissue-derived CD4 T cells and in-vitro polarized CD4 T cell subpopulations, we found that the chromatin accessibility changes of CD4 T cells were associated with a higher predominance of pathogenic Th17 cells (pTh17 cells) in inflamed biopsies. In addition, IBD risk loci in CD4 T cells were colocalized with accessible chromatin changes near pTh17-related genes, as shown in intronic STAT3 and IL23R regions enriched in areas of active intestinal inflammation. Moreover, single cell RNA-Seq analysis revealed a population of pTh17 cells that co-expresses Th1 and cytotoxic transcriptional programs associated with IBD severity. Altogether, we show that cytotoxic pTh17 cells were specifically associated with IBD genetic variants and linked to intestinal inflammation of IBD patients.

Core conserved transcriptional regulatory networks define the invasive trophoblast cell lineage.

The invasive trophoblast cell lineages in rat and human share crucial responsibilities in establishing the uterine-placental interface of the hemochorial placenta. These observations have led to the rat becoming an especially useful animal model for studying hemochorial placentation. However, our understanding of similarities or differences between regulatory mechanisms governing rat and human invasive trophoblast cell populations is limited. In this study, we generated single-nucleus ATAC-seq data from gestation day 15.5 and 19.5 rat uterine-placental interface tissues, and integrated the data with single-cell RNA-seq data generated at the same stages. We determined the chromatin accessibility profiles of invasive trophoblast, natural killer, macrophage, endothelial and smooth muscle cells, and compared invasive trophoblast chromatin accessibility with extravillous trophoblast cell accessibility. In comparing chromatin accessibility profiles between species, we found similarities in patterns of gene regulation and groups of motifs enriched in accessible regions. Finally, we identified a conserved gene regulatory network in invasive trophoblast cells. Our data, findings and analysis will facilitate future studies investigating regulatory mechanisms essential for the invasive trophoblast cell lineage.

Multimodal single-cell analysis of non-random heteroplasmy distribution in human retinal mitochondrial disease.

Variants within the high copy number mitochondrial genome (mtDNA) can disrupt organelle function and lead to severe multi-system disease. The wide range of manifestations observed in mitochondrial disease patients results from varying fractions of abnormal mtDNA molecules in different cells and tissues, a phenomenon termed heteroplasmy. However, the landscape of heteroplasmy across cell types within tissues and its influence on phenotype expression in affected patients remains largely unexplored. Here, we identify non-random distribution of a pathogenic mtDNA variant across a complex tissue using single-cell RNA sequencing, mitochondrial single-cell ATAC sequencing, and multimodal single-cell sequencing. We profile the transcriptome, chromatin accessibility state, and heteroplasmy in cells from the eyes of a patient with mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) and healthy control donors. Utilizing the retina as a model for complex multi-lineage tissues, we found that the proportion of the pathogenic m.3243A>G allele was neither evenly nor randomly distributed across diverse cell types. All neuroectoderm-derived neural cells exhibited a high percentage of the mutant variant. However, a subset of mesoderm-derived lineage, namely the vasculature of the choroid, was near homoplasmic for the wildtype allele. Gene expression and chromatin accessibility profiles of cell types with high and low proportions of m.3243A>G implicate mTOR signaling in the cellular response to heteroplasmy. We further found by multimodal single-cell sequencing of retinal pigment epithelial cells that a high proportion of the pathogenic mtDNA variant was associated with transcriptionally and morphologically abnormal cells. Together, these findings show the non-random nature of mitochondrial variant partitioning in human mitochondrial disease and underscore its implications for mitochondrial disease pathogenesis and treatment.

Single-cell chromatin accessibility profiling of acute myeloid leukemia reveals heterogeneous lineage composition upon therapy-resistance

Acute myeloid leukemia (AML) is a heterogeneous disease characterized by high rate of therapy resistance. Since the cell of origin can impact response to therapy, it is crucial to understand the lineage composition of AML cells at time of therapy resistance. Here we leverage single-cell chromatin accessibility profiling of 22 AML bone marrow aspirates from eight patients at time of therapy resistance and following subsequent therapy to characterize their lineage landscape. Our findings reveal a complex lineage architecture of therapy-resistant AML cells that are primed for stem and progenitor lineages and spanning quiescent, activated and late stem cell/progenitor states. Remarkably, therapy-resistant AML cells are also composed of cells primed for differentiated myeloid, erythroid and even lymphoid lineages. The heterogeneous lineage composition persists following subsequent therapy, with early progenitor-driven features marking unfavorable prognosis in The Cancer Genome Atlas AML cohort. Pseudotime analysis further confirms the vast degree of heterogeneity driven by the dynamic changes in chromatin accessibility. Our findings suggest that therapy-resistant AML cells are characterized not only by stem and progenitor states, but also by a continuum of differentiated cellular lineages. The heterogeneity in lineages likely contributes to their therapy resistance by harboring different degrees of lineage-specific susceptibilities to therapy.

SCENIC+: single-cell multiomic inference of enhancers and gene regulatory networks

Joint profiling of chromatin accessibility and gene expression in individual cells provides an opportunity to decipher enhancer-driven gene regulatory networks (GRNs). Here we present a method for the inference of enhancer-driven GRNs, called SCENIC+. SCENIC+ predicts genomic enhancers along with candidate upstream transcription factors (TFs) and links these enhancers to candidate target genes. To improve both recall and precision of TF identification, we curated and clustered a motif collection with more than 30,000 motifs. We benchmarked SCENIC+ on diverse datasets from different species, including human peripheral blood mononuclear cells, ENCODE cell lines, melanoma cell states and Drosophila retinal development. Next, we exploit SCENIC+ predictions to study conserved TFs, enhancers and GRNs between human and mouse cell types in the cerebral cortex. Finally, we use SCENIC+ to study the dynamics of gene regulation along differentiation trajectories and the effect of TF perturbations on cell state. SCENIC+ is available at scenicplus.readthedocs.io.

Developmental dynamics of the single nucleus regulatory landscape of pig hippocampus.

The hippocampus is a brain region associated with memory, learning and spatial navigation, its aging-related dysfunction is a common sign of Alzheimer's disease. Pig is a good model for human neurodegenerative disease, but our understanding of the regulatory program of the pig hippocampus and its cross-species conservation in humans remains limited. Here, we profiled chromatin accessibility in 33,409 high-quality nuclei and gene expression in 8,122 high-quality nuclei of the pig hippocampus at four postnatal stages. We identified 510,908 accessible chromatin regions (ACRs) in 12 major cell types, among which progenitor cells such as neuroblasts and oligodendrocyte progenitor cells showed a dynamic decrease from early to later developmental stages. We revealed significant enrichment of transposable elements in cell type-specific ACRs, particularly in neuroblasts. We identified oligodendrocytes as the most prominent cell type with the greatest number of genes that showed significant changes during the development. We identified ACRs and key transcription factors underlying the trajectory of neurogenesis (such as POU3F3 and EGR1) and oligodendrocyte differentiation (RXRA and FOXO6). We examined 27 Alzheimer's disease-related genes in our data and found that 15 showed cell type-specific activity (TREM2, RIN3 and CLU), and 15 genes displayed age-associated dynamic activity (BIN1, RABEP1 and APOE). We intersected our data with human genome-wide association study results to detect neurological disease-associated cell types. The present study provides a single nucleus-accessible chromatin landscape of the pig hippocampus at different developmental stages and is helpful for the exploration of pigs as a biomedical model in human neurodegenerative diseases.

RNA-Seq and ATAC-Seq analyses reveal a global transcriptional and chromatin accessibility profiling of γδ T17 differentiation from mouse spleen

IL-17A-producing γδ T cells (γδ T17) are known to play important roles in various autoimmune diseases. However, the molecular mechanisms of γδ T17 differentiation and their functions have not been clarified yet. Here, we sorted IL-17A+ Vγ4, IL-17A- Vγ4, and Vγ1 subsets from mouse spleen by in vitro priming of γδ T17 cells and investigated their differentially expressed genes (DEGs) and differentially accessible regions (DARs) using RNA-seq and ATAC-seq, respectively. Our results showed that DEGs-1 (upregulated genes: 677 and downregulated genes: 821) and DEGs-2 (upregulated genes: 1188 and downregulated genes: 1252) were most closely related to the function and differentiation of peripheral γδ T17. We identified key modules and MCODEs involved in the control of IL-17A+ Vγ4, IL-17A- Vγ4, and Vγ1 subsets using the WGCNA and Metascape analysis. Furthermore, 26 key transcription factors were enriched in three subsets, which contributed to deciphering the potential molecular mechanism driving γδ T17 differentiation. Simultaneously, we conducted chromatin accessibility profiling under γδ T17 differentiation by ATAC-seq. The top six candidate genes were screened for γδ T17 differentiation and function by integrating RNA-seq and ATAC-seq analysis, and the results were further confirmed using RT-qPCR, flow cytometry, and western blot. In addition, the association analysis of candidate genes with the RNA-seq database of psoriasis was performed to elucidate the functional relationship. Our findings provided a novel insight into understanding the molecular mechanisms of γδ T17 differentiation and function and may improve to the development of therapeutic approaches or drugs targeting γδ T17 for autoimmune diseases.

Profiling of Chromatin Accessibility in Pigs across Multiple Tissues and Developmental Stages

The study of chromatin accessibility across tissues and developmental stages is essential for elucidating the transcriptional regulation of various phenotypes and biological processes. However, the chromatin accessibility profiles of multiple tissues in newborn pigs and across porcine liver development remain poorly investigated. Here, we used ATAC-seq and rRNA-depleted RNA-seq to profile open chromatin maps and transcriptional features of heart, kidney, liver, lung, skeletal muscle, and spleen in newborn pigs and porcine liver tissue in the suckling and adult stages, respectively. Specifically, by analyzing a union set of protein-coding genes (PCGs) and two types of transcripts (lncRNAs and TUCPs), we obtained a comprehensive annotation of consensus ATAC-seq peaks for each tissue and developmental stage. As expected, the PCGs with tissue-specific accessible promoters had active transcription and were relevant to tissue-specific functions. In addition, other non-coding tissue-specific peaks were involved in both physical activity and the morphogenesis of neonatal tissues. We also characterized stage-specific peaks and observed a close association between dynamic chromatin accessibility and hepatic function transition during liver postnatal development. Overall, this study expands our current understanding of epigenetic regulation in mammalian tissues and organ development, which can benefit both economic trait improvement and improve the biomedical usage of pigs.

Epigenetically regulated RNA-binding proteins signify malaria hypnozoite dormancy

Dormancy enables relapsing malaria parasites, such as Plasmodium vivax and cynomolgi, to survive unfavorable conditions. It is enabled by hypnozoites, parasites remaining quiescent inside hepatocytes before reactivating and establishing blood-stage infection. We integrate omics approaches to explore gene-regulatory mechanisms underlying hypnozoite dormancy. Genome-wide profiling of activating and repressing histone marks identifies a few genes that get silenced by heterochromatin during hepatic infection of relapsing parasites. By combining single-cell transcriptomics, chromatin accessibility profiling, and fluorescent in situ RNA hybridization, we show that these genes are expressed in hypnozoites and that their silencing precedes parasite development. Intriguingly, these hypnozoite-specific genes mainly encode proteins with RNA-binding domains. We hence hypothesize that these likely repressive RNA-binding proteins keep hypnozoites in a developmentally competent but dormant state and that heterochromatin-mediated silencing of the corresponding genes aids reactivation. Exploring the regulation and exact function of these proteins hence could provide clues for targeted reactivation and killing of these latent pathogens.

Single-cell chromatin accessibility profiling reveals a self-renewing muscle satellite cell state.

A balance between self-renewal and differentiation is critical for the regenerative capacity of tissue-resident stem cells. In skeletal muscle, successful regeneration requires the orchestrated activation, proliferation, and differentiation of muscle satellite cells (MuSCs) that are normally quiescent. A subset of MuSCs undergoes self-renewal to replenish the stem cell pool, but the features that identify and define self-renewing MuSCs remain to be elucidated. Here, through single-cell chromatin accessibility analysis, we reveal the self-renewal versus differentiation trajectories of MuSCs over the course of regeneration in vivo. We identify Betaglycan as a unique marker of self-renewing MuSCs that can be purified and efficiently contributes to regeneration after transplantation. We also show that SMAD4 and downstream genes are genetically required for self-renewal in vivo by restricting differentiation. Our study unveils the identity and mechanisms of self-renewing MuSCs, while providing a key resource for comprehensive analysis of muscle regeneration.

Single-cell chromatin accessibility profiling of cell-state-specific gene regulatory programs during mouse organogenesis.

In mammals, early organogenesis begins soon after gastrulation, accompanied by specification of various type of progenitor/precusor cells. In order to reveal dynamic chromatin landscape of precursor cells and decipher the underlying molecular mechanism driving early mouse organogenesis, we performed single-cell ATAC-seq of E8.5-E10.5 mouse embryos. We profiled a total of 101,599 single cells and identified 41 specific cell types at these stages. Besides, by performing integrated analysis of scATAC-seq and public scRNA-seq data, we identified the critical cis-regulatory elements and key transcription factors which drving development of spinal cord and somitogenesis. Furthermore, we intersected accessible peaks with human diseases/traits-related loci and found potential clinical associated single nucleotide variants (SNPs). Overall, our work provides a fundamental source for understanding cell fate determination and revealing the underlying mechanism during postimplantation embryonic development, and expand our knowledge of pathology for human developmental malformations.

High-throughput identification of regulatory elements and functional assays to uncover susceptibility genes for nasopharyngeal carcinoma

Large-scale genetic association studies have identified multiple susceptibility loci for nasopharyngeal carcinoma (NPC), but the underlying biological mechanisms remain to be explored. To gain insights into the genetic etiology of NPC, we conducted a follow-up study encompassing 6,907 cases and 10,472 controls and identified two additional NPC susceptibility loci, 9q22.33 (rs1867277; OR= 0.74, 95% CI= 0.68-0.81, p= 3.08×10-11) and 17q12 (rs226241; OR= 1.42, 95% CI= 1.26-1.60, p= 1.62×10-8). The two additional loci, together with two previously reported genome-wide significant loci, 5p15.33 and 9p21.3, were investigated by high-throughput sequencing for chromatin accessibility, histone modification, and promoter capture Hi-C (PCHi-C) profiling. Using luciferase reporter assays and CRISPR interference (CRISPRi) to validate the functional profiling, we identified PHF2 at locus 9q22.33 as a susceptibility gene. PHF2 encodes a histone demethylase and acts as a tumor suppressor. The risk alleles of the functional SNPs reduced the expression of the target gene PHF2 by inhibiting the enhancer activity of its long-range (4.3 Mb) cis-regulatory element, which promoted proliferation of NPC cells. In addition, we identified CDKN2B-AS1 as a susceptibility gene at locus 9p21.3, and the NPC risk allele of the functional SNP rs2069418 promoted the expression of CDKN2B-AS1 by increasing its enhancer activity. The overexpression of CDKN2B-AS1 facilitated proliferation of NPC cells. In summary, we identified functional SNPs and NPC susceptibility genes, which provides additional explanations for the genetic association signals and helps to uncover the underlying genetic etiology of NPC development.