In this study, we explored the existence of a transcriptional network co-regulated by E2F7 and HIF1α, as we show that expression of E2F7, like HIF1α, is induced in hypoxia, and because of the previously reported ability of E2F7 to interact with HIF1α. Our genome-wide analysis uncovers a transcriptional network that is directly controlled by HIF1α and E2F7, and demonstrates both stimulatory and repressive functions of the HIF1α -E2F7 complex. Among this network we reveal Neuropilin 1 (NRP1) as a HIF1α-E2F7 repressed gene. By performing in vitro and in vivo reporter assays we demonstrate that the HIF1α-E2F7 mediated NRP1 repression depends on a 41 base pairs 'E2F-binding site hub', providing a molecular mechanism for a previously unanticipated role for HIF1α in transcriptional repression. To explore the biological significance of this regulation we performed in situ hybridizations and observed enhanced nrp1a expression in spinal motorneurons (MN) of zebrafish embryos, upon morpholino-inhibition of e2f7/8 or hif1α Consistent with the chemo-repellent role of nrp1a, morpholino-inhibition of e2f7/8 or hif1α caused MN truncations, which was rescued in TALEN-induced nrp1a(hu10012) mutants, and phenocopied in e2f7/8 mutant zebrafish. Therefore, we conclude that repression of NRP1 by the HIF1α-E2F7 complex regulates MN axon guidance in vivo.

Routine methods for assaying steady-state mRNA levels such as RNA-seq and micro-arrays are commonly used as readouts to study the role of transcription factors (TFs) in gene expression regulation. However, cellular RNA levels do not solely depend on activity of TFs and subsequent transcription by RNA polymerase II (Pol II), but are also affected by RNA turnover rate. Here, we demonstrate that integrated analysis of genome-wide TF occupancy, Pol II binding and steady-state RNA levels provide important insights in gene regulatory mechanisms. Pol II occupancy, as detected by Pol II ChIP-seq, was found to correlate better with TF occupancy compared to steady-state RNA levels and is thus a more precise readout for the primary transcriptional mechanisms that are triggered by signal transduction. Furthermore, analysis of differential Pol II occupancy and RNA-seq levels identified genes with high Pol II occupancy and relatively low RNA levels and vice versa. These categories are strongly enriched for genes from different functional classes. Our results demonstrate a complementary value in Pol II chip-seq and RNA-seq approaches for better understanding of gene expression regulation.

Although the fetal immune system is considered tolerogenic, preterm infants can suffer from severe intestinal inflammation, including necrotizing enterocolitis (NEC). Here, we demonstrate that human fetal intestines predominantly contain tumor necrosis factor-α (TNF-α)+CD4+CD69+ T effector memory (Tem) cells. Single-cell RNA sequencing of fetal intestinal CD4+ T cells showed a T helper 1 phenotype and expression of genes mediating epithelial growth and cell cycling. Organoid co-cultures revealed a dose-dependent, TNF-α-mediated effect of fetal intestinal CD4+ T cells on intestinal stem cell (ISC) development, in which low T cell numbers supported epithelial development, whereas high numbers abrogated ISC proliferation. CD4+ Tem cell frequencies were higher in inflamed intestines from preterm infants with NEC than in healthy infant intestines and showed enhanced TNF signaling. These findings reveal a distinct population of TNF-α-producing CD4+ T cells that promote mucosal development in fetal intestines but can also mediate inflammation upon preterm birth.

Worldwide over 5 million children have been conceived using assisted reproductive technology, and research has concentrated on increasing the likelihood of ongoing pregnancy. However, studies using animal models have indicated undesirable effects of in vitro embryo culture on offspring development and health. In vivo, the oviduct hosts a period in which the early embryo undergoes complete reprogramming of its (epi)genome in preparation for the reacquisition of (epi)genetic marks. We designed an oviduct-on-a-chip platform to better investigate the mechanisms related to (epi)genetic reprogramming and the degree to which they differ between in vitro and in vivo embryos. The device supports more physiological (in vivo-like) zygote genetic reprogramming than conventional IVF. This approach will be instrumental in identifying and investigating factors critical to fertilization and pre-implantation development, which could improve the quality and (epi)genetic integrity of IVF zygotes with likely relevance for early embryonic and later fetal development.

Gastro-intestinal helminth infections trigger the release of interleukin-33 (IL-33), which induces type-2 helper T cells (Th2 cells) at the site of infection to produce IL-13, thereby contributing to host resistance in a T cell receptor (TCR)-independent manner. Here, we show that, as a prerequisite for IL-33-induced IL-13 secretion, Th2 cells required the expression of the epidermal growth factor receptor (EGFR) and of its ligand, amphiregulin, for the formation of a signaling complex between T1/ST2 (the IL-33R) and EGFR. This shared signaling complex allowed IL-33 to induce the EGFR-mediated activation of the MAP-kinase signaling pathway and consequently the expression of IL-13. Lack of EGFR expression on T cells abrogated IL-13 expression in infected tissues and impaired host resistance. EGFR expression on Th2 cells was TCR-signaling dependent, and therefore, our data reveal a mechanism by which antigen presentation controls the innate effector function of Th2 cells at the site of inflammation.

The ability to accurately detect DNA copy number variation in both a sensitive and quantitative manner is important in many research areas. However, genome-wide DNA copy number analyses are complicated by variations in detection signal.

H3K27ac histone acetylome changes contribute to the phenotypic response in heart diseases, particularly in end-stage heart failure. However, such epigenetic alterations have not been systematically investigated in remodeled non-failing human hearts. Therefore, valuable insight into cardiac dysfunction in early remodeling is lacking. This study aimed to reveal the acetylation changes of chromatin regions in response to myocardial remodeling and their correlations to transcriptional changes of neighboring genes.

A missense variant of the sushi, von Willebrand factor type A, EGF and pentraxin domain containing protein 1 (SVEP1) is genome-wide significantly associated with coronary artery disease. The mechanisms how SVEP1 impacts atherosclerosis are not known. We found endothelial cells (EC) and vascular smooth muscle cells to represent the major cellular source of SVEP1 in plaques. Plaques were larger in atherosclerosis-prone Svep1 haploinsufficient (ApoE-/-Svep1+/-) compared to Svep1 wild-type mice (ApoE-/-Svep1+/+) and ApoE-/-Svep1+/- mice displayed elevated plaque neutrophil, Ly6Chigh monocyte, and macrophage numbers. We assessed how leukocytes accumulated more inside plaques in ApoE-/-Svep1+/- mice and found enhanced leukocyte recruitment from blood into plaques. In vitro, we examined how SVEP1 deficiency promotes leukocyte recruitment and found elevated expression of the leukocyte attractant chemokine (C-X-C motif) ligand 1 (CXCL1) in EC after incubation with missense compared to wild-type SVEP1. Increasing wild-type SVEP1 levels silenced endothelial CXCL1 release. In line, plasma Cxcl1 levels were elevated in ApoE-/-Svep1+/- mice. Our studies reveal an atheroprotective role of SVEP1. Deficiency of wild-type Svep1 increased endothelial CXCL1 expression leading to enhanced recruitment of proinflammatory leukocytes from blood to plaque. Consequently, elevated vascular inflammation resulted in enhanced plaque progression in Svep1 deficiency.

The myriad of available hepatocyte in vitro models provides researchers the possibility to select hepatocyte-like cells (HLCs) for specific research goals. However, direct comparison of hepatocyte models is currently challenging. We systematically searched the literature and compared different HLCs, but reported functions were limited to a small subset of hepatic functions. To enable a more comprehensive comparison, we developed an algorithm to compare transcriptomic data across studies that tested HLCs derived from hepatocytes, biliary cells, fibroblasts, and pluripotent stem cells, alongside primary human hepatocytes (PHHs). This revealed that no HLC covered the complete hepatic transcriptome, highlighting the importance of HLC selection. HLCs derived from hepatocytes had the highest transcriptional resemblance to PHHs regardless of the protocol, whereas the quality of fibroblasts and PSC derived HLCs varied depending on the protocol used. Finally, we developed and validated a web application (HLCompR) enabling comparison for specific pathways and addition of new HLCs. In conclusion, our comprehensive transcriptomic comparison of HLCs allows selection of HLCs for specific research questions and can guide improvements in culturing conditions.

Microvascular homeostasis is strictly regulated, requiring close interaction between endothelial cells and pericytes. Here, we aimed to improve our understanding of how microvascular crosstalk affects pericytes. Human-derived pericytes, cultured in absence, or presence of human endothelial cells, were studied by RNA sequencing. Compared with mono-cultured pericytes, a total of 6704 genes were differentially expressed in co-cultured pericytes. Direct endothelial contact induced transcriptome profiles associated with pericyte maturation, suppression of extracellular matrix production, proliferation, and morphological adaptation. In vitro studies confirmed enhanced pericyte proliferation mediated by endothelial-derived PDGFB and pericyte-derived HB-EGF and FGF2. Endothelial-induced PLXNA2 and ACTR3 upregulation also triggered pericyte morphological adaptation. Pathway analysis predicted a key role for TGFβ signaling in endothelial-induced pericyte differentiation, whereas the effect of signaling via gap- and adherens junctions was limited. We demonstrate that endothelial cells have a major impact on the transcriptional profile of pericytes, regulating endothelial-induced maturation, proliferation, and suppression of ECM production.

Tregs are crucial for maintaining maternal immunotolerance against the semiallogeneic fetus. We investigated the elusive transcriptional profile and functional adaptation of human uterine Tregs (uTregs) during pregnancy. Uterine biopsies, from placental bed (materno-fetal interface) and incision site (control) and blood were obtained from women with uncomplicated pregnancies undergoing cesarean section. Tregs and CD4+ non-Tregs were isolated for transcriptomic profiling by Cel-Seq2. Results were validated on protein and single cell levels by flow cytometry. Placental bed uTregs showed elevated expression of Treg signature markers, including FOXP3, CTLA-4, and TIGIT. Their transcriptional profile was indicative of late-stage effector Treg differentiation and chronic activation, with increased expression of immune checkpoints GITR, TNFR2, OX-40, and 4-1BB; genes associated with suppressive capacity (HAVCR2, IL10, LAYN, and PDCD1); and transcription factors MAF, PRDM1, BATF, and VDR. uTregs mirrored non-Treg Th1 polarization and tissue residency. The particular transcriptional signature of placental bed uTregs overlapped strongly with that of tumor-infiltrating Tregs and was remarkably pronounced at the placental bed compared with uterine control site. In conclusion, human uTregs acquire a differentiated effector Treg profile similar to tumor-infiltrating Tregs, specifically at the materno-fetal interface. This introduces the concept of site-specific transcriptional adaptation of Tregs within 1 organ.

Sex differences in endothelial cell (EC) biology may reflect intrinsic differences driven by chromosomes or sex steroid exposure and gender differences accumulated over life. We analysed EC gene expression data from boy-girl twins at birth and in non-twin adults to detect sex differences at different stages of life, and show that 14-25% of the EC transcriptome is sex-biased. By combining data from both stages of life, we identified sex differences that are present at birth and maintained throughout life, and those that are acquired over life. Promisingly, we found that genes that present with an acquired sex difference in ECs are more likely to be targets of sex steroids. Annotating both gene sets with data from multiple genome-wide association studies (GWAS) revealed that genes with an intrinsic sex difference in ECs are enriched for coronary artery disease GWAS hits. This study underscores the need for treating sex as a biological variable.

Gene knockdown techniques are widely used to examine the function of specific genes or proteins. While a variety of techniques are available, a technique commonly used on mammalian oocytes is mRNA knockdown by microinjection of small interfering RNA (siRNA), with non-specific siRNA injection used as a technical control. Here, we investigate whether and how the microinjection procedure itself affects the transcriptome of bovine oocytes. Injection of non-specific siRNA resulted in differential expression of 119 transcripts, of which 76 were down-regulated. Gene ontology analysis revealed that the differentially regulated genes were enriched in the biological processes of ATP synthesis, molecular transport and regulation of protein polyubiquitination. This study establishes a background effect of the microinjection procedure that should be borne in mind by those using microinjection to manipulate gene expression in oocytes.

For many autoimmune diseases, the underlying mechanism is still unknown. In order to get more insight into the etiology of autoimmune diseases, we recently published a study were we performed epigenetic profiling and RNA sequencing on CD4(+)CD45RO(+) T cells derived from the site of inflammation of Juvenile Idiopathic Arthritis (JIA) patients and compared this with healthy controls [1]. In this "Data in Brief", we focus on the analysis of our RNA sequencing data reported in this study, of which the raw and processed files can be found in GEO under GSE71595. We provide a detailed description of the downstream analysis, quality controls, and different analysis methods or techniques that validate the results obtained with RNA-sequencing.

Overlap between non-coding DNA regulatory sequences and common variant associations can help to identify specific cell and tissue types that are relevant for particular diseases. In a systematic manner, we analyzed variants from 94 genome-wide association studies (reporting at least 12 loci at p<5x10-8) by projecting them onto 466 epigenetic datasets (characterizing DNase I hypersensitive sites; DHSs) derived from various adult and fetal tissue samples and cell lines including many biological replicates. We were able to confirm many expected associations, such as the involvement of specific immune cell types in immune-related diseases and tissue types in diseases that affect specific organs, for example, inflammatory bowel disease and coronary artery disease. Other notable associations include adrenal glands in coronary artery disease, the immune system in Alzheimer's disease, and the kidney for bone marrow density. The association signals for some GWAS (for example, myopia or age at menarche) did not show a clear pattern with any of the cell or tissue types studied. In general, the identified variants from GWAS tend to be located outside coding regions. Altogether, we have performed an extensive characterization of GWAS signals in relation to cell and tissue-specific DHSs, demonstrating a key role for regulatory mechanisms in common diseases and complex traits.

Fetal growth restriction (FGR) is associated with an increased susceptibility for various noncommunicable diseases in adulthood, including cardiovascular and renal disease. During FGR, reduced uteroplacental blood flow, oxygen and nutrient supply to the fetus are hypothesized to detrimentally influence cardiovascular and renal programming. This study examined whether developmental programming profiles, especially related to the cardiovascular and renal system, differ in human umbilical vein endothelial cells (HUVECs) collected from pregnancies complicated by placental insufficiency-induced FGR compared to normal growth pregnancies. Our approach, involving transcriptomic profiling by RNA-sequencing and gene set enrichment analysis focused on cardiovascular and renal gene sets and targeted DNA methylation assays, contributes to the identification of targets underlying long-term cardiovascular and renal diseases.

Treg cells are critical regulators of immune homeostasis, and environment-driven Treg cell differentiation into effector (e)Treg cells is crucial for optimal functioning. However, human Treg cell programming in inflammation is unclear. Here, we combine transcriptional and epigenetic profiling to identify a human eTreg cell signature. Inflammation-derived functional Treg cells have a transcriptional profile characterized by upregulation of both a core Treg cell (FOXP3, CTLA4, TIGIT) and effector program (GITR, BLIMP-1, BATF). We identify a specific human eTreg cell signature that includes the vitamin D receptor (VDR) as a predicted regulator in eTreg cell differentiation. H3K27ac/H3K4me1 occupancy indicates an altered (super-)enhancer landscape, including enrichment of the VDR and BATF binding motifs. The Treg cell profile has striking overlap with tumor-infiltrating Treg cells. Our data demonstrate that human inflammation-derived Treg cells acquire a conserved and specific eTreg cell profile guided by epigenetic changes, and fine-tuned by environment-specific adaptations.

The atherosclerotic plaque microenvironment is highly complex, and selective agents that modulate plaque stability are not yet available. We sought to develop a scRNA-seq analysis workflow to investigate this environment and uncover potential therapeutic approaches. We designed a user-friendly, reproducible workflow that will be applicable to other disease-specific scRNA-seq datasets.

How the local inflammatory environment regulates epigenetic changes in the context of inflammatory arthritis remains unclear. Here we assessed the transcriptional and active enhancer profile of monocytes derived from the inflamed joints of JIA patients, a model well-suited for studying inflammatory arthritis.

Women presenting with coronary artery disease (CAD) more often present with fibrous atherosclerotic plaques, which are currently understudied. Phenotypically modulated smooth muscle cells (SMCs) contribute to atherosclerosis in women. How these phenotypically modulated SMCs shape female versus male plaques is unknown. Here, we show sex-stratified gene regulatory networks (GRNs) from human carotid atherosclerotic tissue. Prioritization of these networks identified two main SMC GRNs in late-stage atherosclerosis. Single-cell RNA-sequencing mapped these GRNs to two SMC phenotypes: a phenotypically modulated myofibroblast-like SMC network and a contractile SMC network. The myofibroblast-like GRN was mostly expressed in plaques that were vulnerable in females. Finally, mice orthologs of the female myofibroblast-like genes showed retained expression in advanced plaques from female mice but were downregulated in male mice during atherosclerosis progression. Female atherosclerosis is driven by GRNs that promote a fibrous vulnerable plaque rich in myofibroblast-like SMCs.