Zfp296 negatively regulates H3K9 methylation in embryonic development as a component of heterochromatin

Specific human gene expression in response to infection is an effective marker for diagnosis of latent and active tuberculosis

RNA sequencing and microarray analysis revealed transcriptional markers expressed in whole blood can differentiate active pulmonary TB (ATB) from other respiratory diseases (ORDs), and latent TB infection (LTBI) from healthy controls (HC). Here we describe a streamlined reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) assay that could be applied at near point-of-care for diagnosing and distinguishing ATB from ORDs and LTBI from HC. A literature review was undertaken to identify the most plausible host-gene markers (HGM) of TB infection. Primers, and dual labelled hydrolysis probes were designed and analytically evaluated for accuracy in an in-vitro model of infection using a lung fibroblast cell-line. Best performing genes were multiplexed into panels of 2-4 targets and taken forward for clinical evaluation. Mycobacteria Growth Indicator Tube and QuantiFERON-TB Gold Plus were used as reference tests for ATB and LTBI respectively. A total of 16 HGM were selected and incorporated into five multiplex RT-qPCR panels. PCR assay efficiency of all evaluated targets was ≥ 90% with a median analytical sensitivity of 292 copies/µl [IQR: 215.0-358.3 copies/µl], and a median limit of quantification of 61.7 copies/µl [IQR: 29.4-176.3 copies/µl]. Clinically, ATB was characterised by higher gene expression than ORDs, while LTBI was associated with lower gene expression than HC, Kruskal-Wallis p < 0.0001. Crucially, BATF2, CD64, GBP5, C1QB, GBP6, DUSP3, and GAS6 exhibited high differentiative ability for ATB from ORDs, LTBI or HC while KLF2, PTPRC, NEMF, ASUN, and ZNF296 independently discriminated LTBI from HC. Our results show that different HGM maybe required for ATB and LTBI differentiation from ORDs or HC respectively and demonstrate the feasibility of host gene-based RT-qPCR to diagnose ATB and LTBI at near point-of-care.

Suppression of host gene expression is associated with latent TB infection: a possible diagnostic biomarker

The World Health Organization End TB strategy aims for a 90% reduction of tuberculosis (TB) incidence by 2035. Systematic testing and treatment of latent TB infection (LTBI) among contacts of active TB patients is recommended as one of the ways to curtail TB incidence. However, there is a shortage of tools to accurately diagnose LTBI. We assessed the appropriateness of whole blood host transcriptomic markers (TM) to diagnose LTBI among household contacts of bacteriologically confirmed index cases compared to HIV negative healthy controls (HC). QuantiFERON-TB Gold Plus Interferon gamma release assay (IGRA) and reverse-transcriptase quantitative PCR were used to determine LTBI and quantify TM expression respectively. Association between TM expression and LTBI was evaluated by logistic regression modelling. A total of 100 participants, 49 TB exposed (TBEx) household contacts and 51 HC, were enrolled. Twenty-five (51%) TBEx individuals tested positive by IGRA, and were denoted as LTBI individuals, and 37 (72.5%) HC were IGRA-negative. Expression of 11 evaluated TM was significantly suppressed among LTBI compared to HC. Out of the 11 TM, ZNF296 and KLF2 expression were strongly associated with LTBI and successfully differentiated LTBI from HC. Paradoxically, 21 (49%) TBEx participants who tested IGRA negative exhibited the same pattern of suppressed TM expression as IGRA positive (LTBI-confirmed individuals). Results suggest that suppression of gene expression underlies LTBI and may be a more sensitive diagnostic biomarker than standard-of-care IGRA.

Analysis and Validation of Critical Signatures and Immune Cell Infiltration Characteristics in Doxorubicin-Induced Cardiotoxicity by Integrating Bioinformatics and Machine Learning

Purpose

Doxorubicin-induced cardiotoxicity (DIC) is a severe side reaction in cancer chemotherapy that greatly impacts the well-being of cancer patients. Currently, there is still an insufficiency of effective and reliable biomarkers in the field of clinical practice for the early detection of DIC. This study aimed to determine and validate the potential diagnostic and predictive values of critical signatures in DIC.

Methods

We obtained high-throughput sequencing data from the GEO database and performed data analysis and visualization using R software, GO, KEGG and Cytoscape. Machine learning methods and weighted gene coexpression network (WGCNA) were used to identify key genes for diagnostic model construction. Receiver operating characteristic (ROC) analysis and a nomogram were used to assess their diagnostic values. A multiregulatory network was built to reveal the possible regulatory relationships of critical signatures. Cell-type identification by estimating relative subsets of RNA transcript (CIBERSORT) analysis was used to investigate differential immune cell infiltration. Additionally, a cell and animal model were constructed to investigate the relationship between the identified genes and DIC.

Results

Among the 3713 differentially expressed genes, three key genes (CSGALNACT1, ZNF296 and FANCB) were identified. A nomogram and ROC curves based on three key genes showed excellent diagnostic predictive performance. The regulatory network analysis showed that the TFs CREB1, EP300, FLI1, FOXA1, MAX, and MAZ modulated three key genes. An analysis of immune cell infiltration indicated that many immune cells (activated NK cells, M0 macrophages, activated dendritic cells and neutrophils) might be related to the progression of DIC. Furthermore, there may be various degrees of correlation between the three critical signatures and immune cells. RT‒qPCR demonstrated that the mRNA expression of CSGALNACT1 and ZNF296 was significantly upregulated, while FANCB was significantly downregulated in DOX-treated cardiomyocytes in vitro and in vivo.

Conclusion

Our study suggested that the differential expression of CSGALNACT1, ZNF296 and FANCB is associated with cardiotoxicity and is also involved in immune cell infiltration in DIC. They might be potential biomarkers for the early occurrence of DIC.

Acute irradiation induces a senescence-like chromatin structure in mammalian oocytes

The mechanisms leading to changes in mesoscale chromatin organization during cellular aging are unknown. Here, we used transcriptional activator-like effectors, RNA-seq and superresolution analysis to determine the effects of genotoxic stress on oocyte chromatin structure. Major satellites are organized into tightly packed globular structures that coalesce into chromocenters and dynamically associate with the nucleolus. Acute irradiation significantly enhanced chromocenter mobility in transcriptionally inactive oocytes. In transcriptionally active oocytes, irradiation induced a striking unfolding of satellite chromatin fibers and enhanced the expression of transcripts required for protection from oxidative stress (Fermt1, Smg1), recovery from DNA damage (Tlk2, Rad54l) and regulation of heterochromatin assembly (Zfp296, Ski-oncogene). Non-irradiated, senescent oocytes exhibit not only high chromocenter mobility and satellite distension but also a high frequency of extra chromosomal satellite DNA. Notably, analysis of biological aging using an oocyte-specific RNA clock revealed cellular communication, posttranslational protein modifications, chromatin and histone dynamics as the top cellular processes that are dysregulated in both senescent and irradiated oocytes. Our results indicate that unfolding of heterochromatin fibers following acute genotoxic stress or cellular aging induced the formation of distended satellites and that abnormal chromatin structure together with increased chromocenter mobility leads to chromosome instability in senescent oocytes.

Zfp296 knockout enhances chromatin accessibility and induces a unique state of pluripotency in embryonic stem cells

The Zfp296 gene encodes a zinc finger-type protein. Its expression is high in mouse embryonic stem cells (ESCs) but rapidly decreases following differentiation. Zfp296-knockout (KO) ESCs grew as flat colonies, which were reverted to rounded colonies by exogenous expression of Zfp296. KO ESCs could not form teratomas when transplanted into mice but could efficiently contribute to germline-competent chimeric mice following blastocyst injection. Transcriptome analysis revealed that Zfp296 deficiency up- and down-regulates a distinct group of genes, among which Dppa3, Otx2, and Pou3f1 were markedly downregulated. Chromatin immunoprecipitation sequencing demonstrated that ZFP296 binding is predominantly seen in the vicinity of the transcription start sites (TSSs) of a number of genes, and ZFP296 was suggested to negatively regulate transcription. Consistently, chromatin accessibility assay clearly showed that ZFP296 binding reduces the accessibility of the TSS regions of target genes. Zfp296-KO ESCs showed increased histone H3K9 di- and trimethylation. Co-immunoprecipitation analyses revealed interaction of ZFP296 with G9a and GLP. These results show that ZFP296 plays essential roles in maintaining the global epigenetic state of ESCs through multiple mechanisms including activation of Dppa3, attenuation of chromatin accessibility, and repression of H3K9 methylation, but that Zfp296-KO ESCs retain a unique state of pluripotency while lacking the teratoma-forming ability.

The Novel Role of Zfp296 in Mammalian Embryonic Genome Activation as an H3K9me3 Modulator

The changes in epigenetic modifications during early embryonic development significantly impact mammalian embryonic genome activation (EGA) and are species-conserved to some degree. Here, we reanalyzed the published RNA-Seq of human, mouse, and goat early embryos and found that Zfp296 (zinc finger protein 296) expression was higher at the EGA stage than at the oocyte stage in all three species (adjusted p-value < 0.05 |log2(foldchange)| ≥ 1). Subsequently, we found that Zfp296 was conserved across human, mouse, goat, sheep, pig, and bovine embryos. In addition, we identified that ZFP296 interacts with the epigenetic regulators KDM5B, SMARCA4, DNMT1, DNMT3B, HP1β, and UHRF1. The Cys₂-His₂(C2H2) zinc finger domain TYPE2 TYPE3 domains of ZFP296 co-regulated the modification level of the trimethylation of lysine 9 on the histone H3 protein subunit (H3K9me3). According to ChIP-seq analysis, ZFP296 was also enriched in Trim28, Suv39h1, Setdb1, Kdm4a, and Ehmt2 in the mESC genome. Then, knockdown of the expression of Zfp296 at the late zygote of the mouse led to the early developmental arrest of the mouse embryos and failure resulting from a decrease in H3K9me3. Together, our results reveal that Zfp296 is an H3K9me3 modulator which is essential to the embryonic genome activation of mouse embryos.

Genome activation in equine in vitro-produced embryos

Embryonic genome activation is a critical event in embryo development, in which the transcriptional program of the embryo is initiated. The timing and regulation of this process are species-specific. In vitro embryo production is becoming an important clinical and research tool in the horse; however, very little is known about genome activation in this species. The objective of this work was to identify the timing of genome activation, and the transcriptional networks involved, in in vitro-produced horse embryos. RNA-Seq was performed on oocytes and embryos at eight stages of development (MII, zygote, 2-cell, 4-cell, 8-cell, 16-cell, morula, blastocyst; n = 6 per stage, 2 from each of 3 mares). Transcription of seven genes was initiated at the 2-cell stage. The first substantial increase in gene expression occurred at the 4-cell stage (minor activation), followed by massive gene upregulation and downregulation at the 8-cell stage (major activation). An increase in intronic nucleotides, indicative of transcription initiation, was also observed at the 4-cell stage. Co-expression network analyses identified groups of genes that appeared to be regulated by common mechanisms. Investigation of hub genes and binding motifs enriched in the promoters of co-expressed genes implicated several transcription factors. This work represents, to the best of our knowledge, the first genomic evaluation of embryonic genome activation in horse embryos.

The Nucleosome Remodelling and Deacetylation complex suppresses transcriptional noise during lineage commitment

Multiprotein chromatin remodelling complexes show remarkable conservation of function amongst metazoans, even though components present in invertebrates are often found as multiple paralogous proteins in vertebrate complexes. In some cases, these paralogues specify distinct biochemical and/or functional activities in vertebrate cells. Here, we set out to define the biochemical and functional diversity encoded by one such group of proteins within the mammalian Nucleosome Remodelling and Deacetylation (NuRD) complex: Mta1, Mta2 and Mta3. We find that, in contrast to what has been described in somatic cells, MTA proteins are not mutually exclusive within embryonic stem (ES) cell NuRD and, despite subtle differences in chromatin binding and biochemical interactions, serve largely redundant functions. ES cells lacking all three MTA proteins exhibit complete NuRD loss of function and are viable, allowing us to identify a previously unreported function for NuRD in reducing transcriptional noise, which is essential for maintaining a proper differentiation trajectory during early stages of lineage commitment.

Single cell transcriptome analysis of human, marmoset and mouse embryos reveals common and divergent features of preimplantation development

The mouse embryo is the canonical model for mammalian preimplantation development. Recent advances in single cell profiling allow detailed analysis of embryogenesis in other eutherian species, including human, to distinguish conserved from divergent regulatory programs and signalling pathways in the rodent paradigm. Here, we identify and compare transcriptional features of human, marmoset and mouse embryos by single cell RNA-seq. Zygotic genome activation correlates with the presence of polycomb repressive complexes in all three species, while ribosome biogenesis emerges as a predominant attribute in primate embryos, supporting prolonged translation of maternally deposited RNAs. We find that transposable element expression signatures are species, stage and lineage specific. The pluripotency network in the primate epiblast lacks certain regulators that are operative in mouse, but encompasses WNT components and genes associated with trophoblast specification. Sequential activation of GATA6, SOX17 and GATA4 markers of primitive endoderm identity is conserved in primates. Unexpectedly, OTX2 is also associated with primitive endoderm specification in human and non-human primate blastocysts. Our cross-species analysis demarcates both conserved and primate-specific features of preimplantation development, and underscores the molecular adaptability of early mammalian embryogenesis.

Highlighted Article: Analysis of stage-matched, single-cell gene expression data from three mammalian species reveals conserved and primate-specific regulation of early embryogenesis and lineage specification.

NuRD-interacting protein ZFP296 regulates genome-wide NuRD localization and differentiation of mouse embryonic stem cells

The nucleosome remodeling and deacetylase (NuRD) complex plays an important role in gene expression regulation, stem cell self-renewal, and lineage commitment. However, little is known about the dynamics of NuRD during cellular differentiation. Here, we study these dynamics using genome-wide profiling and quantitative interaction proteomics in mouse embryonic stem cells (ESCs) and neural progenitor cells (NPCs). We find that the genomic targets of NuRD are highly dynamic during differentiation, with most binding occurring at cell-type specific promoters and enhancers. We identify ZFP296 as an ESC-specific NuRD interactor that also interacts with the SIN3A complex. ChIP-sequencing in Zfp296 knockout (KO) ESCs reveals decreased NuRD binding both genome-wide and at ZFP296 binding sites, although this has little effect on the transcriptome. Nevertheless, Zfp296 KO ESCs exhibit delayed induction of lineage-specific markers upon differentiation to embryoid bodies. In summary, we identify an ESC-specific NuRD-interacting protein which regulates genome-wide NuRD binding and cellular differentiation.

The NuRD complex plays an important role in regulating lineage commitment and cell fate during early embryonic development. Here the authors present an integrative analysis of MBD3/NuRD composition and binding in mouse embryonic stem cells and neural progenitor cells, providing a molecular basis for genome-wide NuRD localization

Tracing the transitions from pluripotency to germ cell fate with CRISPR screening

Early mammalian development entails transit through naive pluripotency towards post-implantation epiblast, which subsequently gives rise to primordial germ cells (PGC), the founding germline population. To investigate these cell fate transitions, we developed a compound-reporter to track cellular identity in a model of PGC specification (PGC-like cells; PGCLC), and coupled it with genome-wide CRISPR screening. We identify key genes both for exit from pluripotency and for acquisition of PGC fate, and characterise a central role for the transcription regulators Nr5a2 and Zfp296 in germline ontogeny. Abrogation of these genes results in widespread activation (Nr5a2^−/−) or inhibition (Zfp296^−/−) of WNT pathway factors in PGCLC. This leads to aberrant upregulation of the somatic programme or failure to activate germline genes, respectively, and consequently loss of germ cell identity. Our study places Zfp296 and Nr5a2 as key components of an expanded PGC gene regulatory network, and outlines a transferable strategy for identifying critical regulators of complex cell fate decisions.

Primordial Germ Cell-Like Cells (PGCLCs) are an in vitro model for primordial germ cell development. Here, the authors couple a novel compound reporter with CRISPR screening to identify key genes for exit from pluripotency and acquisition of PGCLC fate; specifically identifying Nr5a2 and Zfp296.

Integrated analysis of single-cell embryo data yields a unified transcriptome signature for the human pre-implantation epiblast

Single-cell profiling techniques create opportunities to delineate cell fate progression in mammalian development. Recent studies have provided transcriptome data from human pre-implantation embryos, in total comprising nearly 2000 individual cells. Interpretation of these data is confounded by biological factors, such as variable embryo staging and cell-type ambiguity, as well as technical challenges in the collective analysis of datasets produced with different sample preparation and sequencing protocols. Here, we address these issues to assemble a complete gene expression time course spanning human pre-implantation embryogenesis. We identify key transcriptional features over developmental time and elucidate lineage-specific regulatory networks. We resolve post-hoc cell-type assignment in the blastocyst, and define robust transcriptional prototypes that capture epiblast and primitive endoderm lineages. Examination of human pluripotent stem cell transcriptomes in this framework identifies culture conditions that sustain a naïve state pertaining to the inner cell mass. Our approach thus clarifies understanding both of lineage segregation in the early human embryo and of in vitro stem cell identity, and provides an analytical resource for comparative molecular embryology.