Integration of a multi-omics stem cell differentiation dataset using a dynamical model.

Stem cell differentiation is a highly dynamic process involving pervasive changes in gene expression. The large majority of existing studies has characterized differentiation at the level of individual molecular profiles, such as the transcriptome or the proteome. To obtain a more comprehensive view, we measured protein, mRNA and microRNA abundance during retinoic acid-driven differentiation of mouse embryonic stem cells. We found that mRNA and protein abundance are typically only weakly correlated across time. To understand this finding, we developed a hierarchical dynamical model that allowed us to integrate all data sets. This model was able to explain mRNA-protein discordance for most genes and identified instances of potential microRNA-mediated regulation. Overexpression or depletion of microRNAs identified by the model, followed by RNA sequencing and protein quantification, were used to follow up on the predictions of the model. Overall, our study shows how multi-omics integration by a dynamical model could be used to nominate candidate regulators.

Mpox outbreak in the Netherlands, 2022: public health response, characteristics of the first 1,000 cases and protection of the first-generation smallpox vaccine.

In early May 2022, a global outbreak of mpox started among persons without travel history to regions known to be enzootic for monkeypox virus (MPXV). On 8 August 2022, the Netherlands reported its 1,000th mpox case, representing a cumulative incidence of 55 per million population, one of the highest cumulative incidences worldwide. We describe characteristics of the first 1,000 mpox cases in the Netherlands, reported between 20 May and 8 August 2022, within the context of the public health response. These cases were predominantly men who have sex with men aged 31-45 years. The vast majority of infections were acquired through sexual contact with casual partners in private or recreational settings including LGBTQIA+ venues in the Netherlands. This indicates that, although some larger upsurges occurred from point-source and/or travel-related events, the outbreak was mainly characterised by sustained transmission within the Netherlands. In addition, we estimated the protective effect of first-generation smallpox vaccine against moderate/severe mpox and found a vaccine effectiveness of 58% (95% CI: 17-78%), suggesting moderate protection against moderate/severe mpox symptoms on top of any possible protection by this vaccine against MPXV infection and disease. Communication with and supporting the at-risk population in following mitigation measures remains essential.

A gastruloid model of the interaction between embryonic and extra-embryonic cell types.

Stem-cell derived in vitro systems, such as organoids or embryoids, hold great potential for modeling in vivo development. Full control over their initial composition, scalability, and easily measurable dynamics make those systems useful for studying specific developmental processes in isolation. Here we report the formation of gastruloids consisting of mouse embryonic stem cells (mESCs) and extraembryonic endoderm (XEN) cells. These XEN-enhanced gastruloids (XEGs) exhibit the formation of neural epithelia, which are absent in gastruloids derived from mESCs only. By single-cell RNA-seq, imaging, and differentiation experiments, we demonstrate the neural characteristics of the epithelial tissue. We further show that the mESCs induce the differentiation of the XEN cells to a visceral endoderm-like state. Finally, we demonstrate that local inhibition of WNT signaling and production of a basement membrane by the XEN cells underlie the formation of the neuroepithelial tissue. In summary, we establish XEGs to explore heterotypic cellular interactions and their developmental consequences in vitro.

Dynamic Enhancer DNA Methylation as Basis for Transcriptional and Cellular Heterogeneity of ESCs.

Variable levels of DNA methylation have been reported at tissue-specific differential methylation regions (DMRs) overlapping enhancers, including super-enhancers (SEs) associated with key cell identity genes, but the mechanisms responsible for this intriguing behavior are not well understood. We used allele-specific reporters at the endogenous Sox2 and Mir290 SEs in embryonic stem cells and found that the allelic DNA methylation state is dynamically switching, resulting in cell-to-cell heterogeneity. Dynamic DNA methylation is driven by the balance between DNA methyltransferases and transcription factor binding on one side and co-regulated with the Mediator complex recruitment and H3K27ac level changes at regulatory elements on the other side. DNA methylation at the Sox2 and the Mir290 SEs is independently regulated and has distinct consequences on the cellular differentiation state. Dynamic allele-specific DNA methylation at the two SEs was also seen at different stages in preimplantation embryos, revealing that methylation heterogeneity occurs in vivo.

Single-cell transcriptomics reveals gene expression dynamics of human fetal kidney development.

The current understanding of mammalian kidney development is largely based on mouse models. Recent landmark studies revealed pervasive differences in renal embryogenesis between mouse and human. The scarcity of detailed gene expression data in humans therefore hampers a thorough understanding of human kidney development and the possible developmental origin of kidney diseases. In this paper, we present a single-cell transcriptomics study of the human fetal kidney. We identified 22 cell types and a host of marker genes. Comparison of samples from different developmental ages revealed continuous gene expression changes in podocytes. To demonstrate the usefulness of our data set, we explored the heterogeneity of the nephrogenic niche, localized podocyte precursors, and confirmed disease-associated marker genes. With close to 18,000 renal cells from five different developmental ages, this study provides a rich resource for the elucidation of human kidney development, easily accessible through an interactive web application.

The miR-96 and RARγ signaling axis governs androgen signaling and prostate cancer progression.

Expression levels of retinoic acid receptor gamma (NR1B3/RARG, encodes RARγ) are commonly reduced in prostate cancer (PCa). Therefore, we sought to establish the cellular and gene regulatory consequences of reduced RARγ expression, and determine RARγ regulatory mechanisms. RARG shRNA approaches in non-malignant (RWPE-1 and HPr1-AR) and malignant (LNCaP) prostate models revealed that reducing RARγ levels, rather than adding exogenous retinoid ligand, had the greatest impact on prostate cell viability and gene expression. ChIP-Seq defined the RARγ cistrome, which was significantly enriched at active enhancers associated with AR binding sites. Reflecting a significant genomic role for RARγ to regulate androgen signaling, RARγ knockdown in HPr1-AR cells significantly regulated the magnitude of the AR transcriptome. RARγ downregulation was explained by increased miR-96 in PCa cell and mouse models, and TCGA PCa cohorts. Biochemical approaches confirmed that miR-96 directly regulated RARγ expression and function. Capture of the miR-96 targetome by biotin-miR-96 identified that RARγ and a number of RARγ interacting co-factors including TACC1 were all targeted by miR-96, and expression of these genes were prominently altered, positively and negatively, in the TCGA-PRAD cohort. Differential gene expression analyses between tumors in the TCGA-PRAD cohort with lower quartile expression levels of RARG and TACC1 and upper quartile miR-96, compared to the reverse, identified a gene network including several RARγ target genes (e.g., SOX15) that significantly associated with worse disease-free survival (hazard ratio 2.23, 95% CI 1.58 to 2.88, p = 0.015). In summary, miR-96 targets a RARγ network to govern AR signaling, PCa progression and disease outcome.

Allele-specific repression of Sox2 through the long non-coding RNA Sox2ot.

The transcription factor Sox2 controls the fate of pluripotent stem cells and neural stem cells. This gatekeeper function requires well-regulated Sox2 levels. We postulated that Sox2 regulation is partially controlled by the Sox2 overlapping long non-coding RNA (lncRNA) gene Sox2ot. Here we show that the RNA levels of Sox2ot and Sox2 are inversely correlated during neural differentiation of mouse embryonic stem cells (ESCs). Through allele-specific enhanced transcription of Sox2ot in mouse Sox2eGFP knockin ESCs we demonstrate that increased Sox2ot transcriptional activity reduces Sox2 RNA levels in an allele-specific manner. Enhanced Sox2ot transcription, yielding lower Sox2 RNA levels, correlates with a decreased chromatin interaction of the upstream regulatory sequence of Sox2 and the ESC-specific Sox2 super enhancer. Our study indicates that, in addition to previously reported in trans mechanisms, Sox2ot can regulate Sox2 by an allele-specific mechanism, in particular during development.

Integration of VDR genome wide binding and GWAS genetic variation data reveals co-occurrence of VDR and NF-κB binding that is linked to immune phenotypes.

BACKGROUND: The nuclear hormone receptor superfamily acts as a genomic sensor of diverse signals. Their actions are often intertwined with other transcription factors. Nuclear hormone receptors are targets for many therapeutic drugs, and include the vitamin D receptor (VDR). VDR signaling is pleotropic, being implicated in calcaemic function, antibacterial actions, growth control, immunomodulation and anti-cancer actions. Specifically, we hypothesized that the biologically significant relationships between the VDR transcriptome and phenotype-associated biology could be discovered by integrating the known VDR transcription factor binding sites and all published trait- and disease-associated SNPs. By integrating VDR genome-wide binding data (ChIP-seq) with the National Human Genome Research Institute (NHGRI) GWAS catalog of SNPs we would see where and which target gene interactions and pathways are impacted by inherited genetic variation in VDR binding sites, indicating which of VDR's multiple functions are most biologically significant. RESULTS: To examine how genetic variation impacts VDR function we overlapped 23,409 VDR genomic binding peaks from six VDR ChIP-seq datasets with 191,482 SNPs, derived from GWAS-significant SNPs (Lead SNPs) and their correlated variants (r 2 > 0.8) from HapMap3 and the 1000 genomes project. In total, 574 SNPs (71 Lead and 503 SNPs in linkage disequilibrium with Lead SNPs) were present at VDR binding loci and associated with 211 phenotypes. For each phenotype a hypergeometric test was used to determine if SNPs were enriched at VDR binding sites. Bonferroni correction for multiple testing across the 211 phenotypes yielded 42 SNPs that were either disease- or phenotype-associated with seven predominately immune related including self-reported allergy; esophageal cancer was the only cancer phenotype. Motif analyses revealed that only two of these 42 SNPs reside within a canonical VDR binding site (DR3 motif), and that 1/3 of the 42 SNPs significantly impacted binding and gene regulation by other transcription factors, including NF-κB. This suggests a plausible link for the potential cross-talk between VDR and NF-κB. CONCLUSIONS: These analyses showed that VDR peaks are enriched for SNPs associated with immune phenotypes suggesting that VDR immunomodulatory functions are amongst its most important actions. The enrichment of genetic variation in non-DR3 motifs suggests a significant role for the VDR to bind in multimeric complexes containing other transcription factors that are the primary DNA binding component. Our work provides a framework for the combination of ChIP-seq and GWAS findings to provide insight into the underlying phenotype-associated biology of a given transcription factor.

Integrative genomic analysis in K562 chronic myelogenous leukemia cells reveals that proximal NCOR1 binding positively regulates genes that govern erythroid differentiation and Imatinib sensitivity.

To define the functions of NCOR1 we developed an integrative analysis that combined ENCODE and NCI-60 data, followed by in vitro validation. NCOR1 and H3K9me3 ChIP-Seq, FAIRE-seq and DNA CpG methylation interactions were related to gene expression using bootstrapping approaches. Most NCOR1 combinations (24/44) were associated with significantly elevated level expression of protein coding genes and only very few combinations related to gene repression. DAVID's biological process annotation revealed that elevated gene expression was uniquely associated with acetylation and ETS binding. A matrix of gene and drug interactions built on NCI-60 data identified that Imatinib significantly targeted the NCOR1 governed transcriptome. Stable knockdown of NCOR1 in K562 cells slowed growth and significantly repressed genes associated with NCOR1 cistrome, again, with the GO terms acetylation and ETS binding, and significantly dampened sensitivity to Imatinib-induced erythroid differentiation. Mining public microarray data revealed that NCOR1-targeted genes were significantly enriched in Imatinib response gene signatures in cell lines and chronic myelogenous leukemia (CML) patients. These approaches integrated cistrome, transcriptome and drug sensitivity relationships to reveal that NCOR1 function is surprisingly most associated with elevated gene expression, and that these targets, both in CML cell lines and patients, associate with sensitivity to Imatinib.