Extracellular vesicles from mouse trophoblast cells: Effects on neural progenitor cells and potential participants in the placenta-brain axis†

The fetal brain of the mouse is thought to be dependent upon the placenta as a source of serotonin (5-hydroxytryptamine; 5-HT) and other factors. How factors reach the developing brain remains uncertain but are postulated here to be part of the cargo carried by placental extracellular vesicles (EV). We have analyzed the protein, catecholamine, and small RNA content of EV from mouse trophoblast stem cells (TSC) and TSC differentiated into parietal trophoblast giant cells (pTGC), potential primary purveyors of 5-HT. Current studies examined how exposure of mouse neural progenitor cells (NPC) to EV from either TSC or pTGC affect their transcriptome profiles. The EV from trophoblast cells contained relatively high amounts of 5-HT, as well as dopamine and norepinephrine, but there were no significant differences between EV derived from pTGC and from TSC. Content of miRNA and small nucleolar (sno)RNA, however, did differ according to EV source, and snoRNA were upregulated in EV from pTGC. The primary inferred targets of the microRNA (miRNA) from both pTGC and TSC were mRNA enriched in the fetal brain. NPC readily internalized EV, leading to changes in their transcriptome profiles. Transcripts regulated were mainly ones enriched in neural tissues. The transcripts in EV-treated NPC that demonstrated a likely complementarity with miRNA in EV were mainly up- rather than downregulated, with functions linked to neuronal processes. Our results are consistent with placenta-derived EV providing direct support for fetal brain development and being an integral part of the placenta-brain axis.

Conditionally mutant animal model for investigating the invasive trophoblast cell lineage

Placental development involves coordinated expansion and differentiation of trophoblast cell lineages possessing specialized functions. Among the differentiated trophoblast cell lineages are invasive trophoblast cells, which exit the placenta and invade the uterus, where they restructure the uterine parenchyma and facilitate remodeling of uterine spiral arteries. The rat exhibits deep intrauterine trophoblast cell invasion, a feature shared with human placentation, and is also amenable to gene manipulation using genome-editing techniques. In this investigation, we generated a conditional rat model targeting the invasive trophoblast cell lineage. Prolactin family 7, subfamily b, member 1 (Prl7b1) is uniquely and abundantly expressed in the rat invasive trophoblast cell lineage. Disruption of Prl7b1 did not adversely affect placental development. We demonstrated that the Prl7b1 locus could be effectively used to drive the expression of Cre recombinase in invasive trophoblast cells. Our rat model represents a new tool for investigating candidate genes contributing to the regulation of invasive trophoblast cells and their roles in trophoblast-guided uterine spiral artery remodeling.

CONDITIONALLY MUTANT ANIMAL MODEL FOR INVESTIGATING THE INVASIVE TROPHOBLAST CELL LINEAGE

Placental development involves coordinated expansion and differentiation of trophoblast cell lineages possessing specialized functions. Among the differentiated trophoblast cell lineages are invasive trophoblast cells, which exit the placenta and invade into the uterus where they restructure the uterine parenchyma and facilitate remodeling of uterine spiral arteries. The rat exhibits deep intrauterine trophoblast cell invasion, a feature shared with human placentation, and is also amenable to gene manipulation using genome editing techniques. In this investigation, we generated a conditional rat model targeting the invasive trophoblast cell lineage. Prolactin family 7, subfamily b, member 1 ( Prl7b1 ) is uniquely and abundantly expressed in the rat invasive trophoblast cell lineage. Disruption of Prl7b1 did not adversely affect placental development. We demonstrated that the Prl7b1 locus could be effectively used to drive the expression of Cre recombinase in invasive trophoblast cells. Our rat model represents a new tool for investigating candidate genes contributing to the regulation of invasive trophoblast cells and their contributions to trophoblast-guided uterine spiral artery remodeling.

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.

Unsupervised contrastive peak caller for ATAC-seq

The assay for transposase-accessible chromatin with sequencing (ATAC-seq) is a common assay to identify chromatin accessible regions by using a Tn5 transposase that can access, cut, and ligate adapters to DNA fragments for subsequent amplification and sequencing. These sequenced regions are quantified and tested for enrichment in a process referred to as "peak calling." Most unsupervised peak calling methods are based on simple statistical models and suffer from elevated false positive rates. Newly developed supervised deep learning methods can be successful, but they rely on high quality labeled data for training, which can be difficult to obtain. Moreover, though biological replicates are recognized to be important, there are no established approaches for using replicates in the deep learning tools, and the approaches available for traditional methods either cannot be applied to ATAC-seq, where control samples may be unavailable, or are post hoc and do not capitalize on potentially complex, but reproducible signal in the read enrichment data. Here, we propose a novel peak caller that uses unsupervised contrastive learning to extract shared signals from multiple replicates. Raw coverage data are encoded to obtain low-dimensional embeddings and optimized to minimize a contrastive loss over biological replicates. These embeddings are passed to another contrastive loss for learning and predicting peaks and decoded to denoised data under an autoencoder loss. We compared our replicative contrastive learner (RCL) method with other existing methods on ATAC-seq data, using annotations from ChromHMM genomic labels and transcription factor ChIP-seq as noisy truth. RCL consistently achieved the best performance.

CORE CONSERVED TRANSCRIPTIONAL REGULATORY NETWORKS DEFINE THE INVASIVE TROPHOBLAST CELL LINEAGE

The invasive trophoblast cell lineage 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 to study 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 (sn) ATAC-seq data from gestation day (gd) 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 to extravillous trophoblast (EVT) 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.

Identifying novel regulators of placental development using time-series transcriptome data

The placenta serves as a connection between the mother and the fetus during pregnancy, providing the fetus with oxygen, nutrients, and growth hormones. However, the regulatory mechanisms and dynamic gene interaction networks underlying early placental development are understudied. Here, we generated RNA-sequencing data from mouse fetal placenta at embryonic days 7.5, 8.5, and 9.5 to identify genes with timepoint-specific expression, then inferred gene interaction networks to analyze highly connected network modules. We determined that timepoint-specific gene network modules were associated with distinct developmental processes, and with similar expression profiles to specific human placental cell populations. From each module, we identified hub genes and their direct neighboring genes, which were predicted to govern placental functions. We confirmed that four novel candidate regulators identified through our analyses regulate cell migration in the HTR-8/SVneo cell line. Overall, we predicted several novel regulators of placental development expressed in specific placental cell types using network analysis of bulk RNA-sequencing data. Our findings and analysis approaches will be valuable for future studies investigating the transcriptional landscape of early development.

Unsupervised Contrastive Peak Caller for ATAC-seq

The assay for transposase-accessible chromatin with sequencing (ATAC-seq) is a common assay to identify chromatin accessible regions by using a Tn5 transposase that can access, cut, and ligate adapters to DNA fragments for subsequent amplification and sequencing. These sequenced regions are quantified and tested for enrichment in a process referred to as "peak calling". Most unsupervised peak calling methods are based on simple statistical models and suffer from elevated false positive rates. Newly developed supervised deep learning methods can be successful, but they rely on high quality labeled data for training, which can be difficult to obtain. Moreover, though biological replicates are recognized to be important, there are no established approaches for using replicates in the deep learning tools, and the approaches available for traditional methods either cannot be applied to ATAC-seq, where control samples may be unavailable, or are post-hoc and do not capitalize on potentially complex, but reproducible signal in the read enrichment data. Here, we propose a novel peak caller that uses unsupervised contrastive learning to extract shared signals from multiple replicates. Raw coverage data are encoded to obtain low-dimensional embeddings and optimized to minimize a contrastive loss over biological replicates. These embeddings are passed to another contrastive loss for learning and predicting peaks and decoded to denoised data under an autoencoder loss. We compared our Replicative Contrastive Learner (RCL) method with other existing methods on ATAC-seq data, using annotations from ChromHMM genome and transcription factor ChIP-seq as noisy truth. RCL consistently achieved the best performance.

Transcriptional regulation of Satb1 in mouse trophoblast stem cells

SATB homeobox proteins are important regulators of developmental gene expression. Among the stem cell lineages that emerge during early embryonic development, trophoblast stem (TS) cells exhibit robust SATB expression. Both SATB1 and SATB2 act to maintain the trophoblast stem-state. However, the molecular mechanisms that regulate TS-specific Satb expression are not yet known. We identified Satb1 variant 2 as the predominant transcript in trophoblasts. Histone marks, and RNA polymerase II occupancy in TS cells indicated an active state of the promoter. A novel cis-regulatory region with active histone marks was identified ∼21 kbp upstream of the variant 2 promoter. CRISPR/Cas9 mediated disruption of this sequence decreased Satb1 expression in TS cells and chromosome conformation capture analysis confirmed looping of this distant regulatory region into the proximal promoter. Scanning position weight matrices across the enhancer predicted two ELF5 binding sites in close proximity to SATB1 sites, which were confirmed by chromatin immunoprecipitation. Knockdown of ELF5 downregulated Satb1 expression in TS cells and overexpression of ELF5 increased the enhancer-reporter activity. Interestingly, ELF5 interacts with SATB1 in TS cells, and the enhancer activity was upregulated following SATB overexpression. Our findings indicate that trophoblast-specific Satb1 expression is regulated by long-range chromatin looping of an enhancer that interacts with ELF5 and SATB proteins.

The product of BMP-directed differentiation protocols for human primed pluripotent stem cells is placental trophoblast and not amnion

The observation that trophoblast (TB) can be generated from primed pluripotent stem cells (PSCs) by exposure to bone morphogenetic protein-4 (BMP4) when FGF2 and ACTIVIN signaling is minimized has recently been challenged with the suggestion that the procedure instead produces amnion. Here, by analyzing transcriptome data from multiple sources, including bulk and single-cell data, we show that the BMP4 procedure generates bona fide TB with similarities to both placental villous TB and TB generated from TB stem cells. The analyses also suggest that the transcriptomic signatures between embryonic amnion and different forms of TB have commonalities. Our data provide justification for the continued use of TB derived from PSCs as a model for investigating placental development.

•

Cells differentiated by using BAP protocols resemble TB more than embryonic amnion

•

Deviation from the standard BAP protocol results in less differentiated TB

•

Single-cell/nucleus RNA-seq analysis identifies two syncytiotrophoblast populations

The impact of bisphenol a (BPA) on the placenta

Bisphenol A (BPA) is an endocrine-disrupting chemical (EDC) that is used in a wide-variety of plastic and common house-hold items. Therefore, there is potential continual exposure to this compound. BPA exposure has been linked to certain placenta-associated obstetric complications such as preeclampsia, fetal growth restriction, miscarriage, and preterm birth. However, how BPA exposure results in these disorders remains uncertain. Hence, we have herein summarized the reported impact of BPA on the morphology and metabolic state of the placenta and have proposed mechanisms by which BPA affects placentation, potentially leading to obstetric complications. Current findings suggest that BPA induces pathological changes in the placenta and disrupts its metabolic activities. Based on exposure concentrations, BPA can elicit apoptotic or anti-apoptotic signals in the trophoblasts; and can exaggerate trophoblast fusion while inhibiting trophoblast migration and invasion to affect pregnancy. Accordingly, the usage of BPA products by pregnant women should be minimized and less harmful alternative chemicals should be explored and employed where possible.

Mapping cis-regulatory elements in the midgestation mouse placenta

The placenta is a temporary organ that provides the developing fetus with nutrients, oxygen, and protection in utero. Defects in its development, which may be caused by misregulated gene expression, can lead to devastating outcomes for the mother and fetus. In mouse, placental defects during midgestation commonly lead to embryonic lethality. However, the regulatory mechanisms controlling expression of genes during this period have not been thoroughly investigated. Therefore, we generated and analyzed ChIP-seq data for multiple histone modifications known to mark cis-regulatory regions. We annotated active and poised promoters and enhancers, as well as regions generally associated with repressed gene expression. We found that poised promoters were associated with neuronal development genes, while active promoters were largely associated with housekeeping genes. Active and poised enhancers were associated with placental development genes, though only active enhancers were associated with genes that have placenta-specific expression. Motif analysis within active enhancers identified a large network of transcription factors, including those that have not been previously studied in the placenta and are candidates for future studies. The data generated and genomic regions annotated provide researchers with a foundation for future studies, aimed at understanding how specific genes in the midgestation mouse placenta are regulated.

Spatial Transcriptomics analysis of uterine gene expression in enhancer of Zeste homolog 2 (Ezh2) conditional knockout mice

Histone proteins undergo various modifications that alter chromatin structure, including addition of methyl groups. Enhancer of homolog 2 (EZH2), is a histone methyltransferase that methylates lysine residue 27, and thereby, suppresses gene expression. EZH2 plays integral role in the uterus and other reproductive organs. We have previously shown that conditional deletion of uterine EZH2 results in increased proliferation of luminal and glandular epithelial cells, and RNAseq analyses reveal several uterine transcriptomic changes in Ezh2 conditional (c) knockout (KO) mice that can affect estrogen signaling pathways. To pinpoint the origin of such gene expression changes, we used the recently developed spatial transcriptomics (ST) method with the hypotheses that Ezh2cKO mice would predominantly demonstrate changes in epithelial cells and/or ablation of this gene would disrupt normal epithelial/stromal gene expression patterns. Uteri were collected from ovariectomized adult WT and Ezh2cKO mice and analyzed by ST. Asb4, Cxcl14, Dio2, and Igfbp5 were increased, Sult1d1, Mt3, and Lcn2 were reduced in Ezh2cKO uterine epithelium vs. WT epithelium. For Ezh2cKO uterine stroma, differentially expressed key hub genes included Cald1, Fbln1, Myh11, Acta2, and Tagln. Conditional loss of uterine Ezh2 also appears to shift the balance of gene expression profiles in epithelial vs. stromal tissue toward uterine epithelial cell and gland development and proliferation, consistent with uterine gland hyperplasia in these mice. Current findings provide further insight into how EZH2 may selectively affect uterine epithelial and stromal compartments. Additionally, these transcriptome data might provide the mechanistic understanding and valuable biomarkers for human endometrial disorders with epigenetic underpinnings.

Single Nucleus RNA Sequence (snRNAseq) Analysis of the Spectrum of Trophoblast Lineages Generated From Human Pluripotent Stem Cells in vitro

One model to study the emergence of the human trophoblast (TB) has been the exposure of pluripotent stem cells to bone morphogenetic protein 4 (BMP4) in presence of inhibitors of ACTIVIN/TGFB; A83-01 and FGF2; PD173074 (BAP), which generates a mixture of cytotrophoblast, syncytiotrophoblast, and cells with similarities to extravillous trophoblast. Here, H1 human embryonic stem cells were BAP-exposed under two O2 conditions (20% and 5%, respectively). At day 8, single nuclei RNA sequencing was used for transcriptomics analysis, thereby allowing profiling of fragile syncytial structures as well as the more resilient mononucleated cells. Following cluster analysis, two major groupings, one comprised of five (2,4,6,7,8) and the second of three (1,3,5) clusters were evident, all of which displayed recognized TB markers. Of these, two (2 and 3) weakly resembled extravillous trophoblast, two (5 and 6) strongly carried the hallmark transcripts of syncytiotrophoblast, while the remaining five were likely different kinds of mononucleated cytotrophoblast. We suggest that the two populations of nuclei within syncytiotrophoblast may have arisen from fusion events involving two distinct species of precursor cells. The number of differentially expressed genes between O2 conditions varied among the clusters, and the number of genes upregulated in cells cultured under 5% O2 was highest in syncytiotrophoblast cluster 6. In summary, the BAP model reveals an unexpectedly complex picture of trophoblast lineage emergence that will need to be resolved further in time-course studies.

Transcription Factor PLAGL1 Is Associated with Angiogenic Gene Expression in the Placenta

During pregnancy, the placenta is important for transporting nutrients and waste between the maternal and fetal blood supply, secreting hormones, and serving as a protective barrier. To better understand placental development, we must understand how placental gene expression is regulated. We used RNA-seq data and ChIP-seq data for the enhancer associated mark, H3k27ac, to study gene regulation in the mouse placenta at embryonic day (e) 9.5, when the placenta is developing a complex network of blood vessels. We identified several upregulated transcription factors with enriched binding sites in e9.5-specific enhancers. The most enriched transcription factor, PLAGL1 had a predicted motif in 233 regions that were significantly associated with vasculature development and response to insulin stimulus genes. We then performed several experiments using mouse placenta and a human trophoblast cell line to understand the role of PLAGL1 in placental development. In the mouse placenta, Plagl1 is expressed in endothelial cells of the labyrinth layer and is differentially expressed in placentas from mice with gestational diabetes compared to placentas from control mice in a sex-specific manner. In human trophoblast cells, siRNA knockdown significantly decreased expression of genes associated with placental vasculature development terms. In a tube assay, decreased PLAGL1 expression led to reduced cord formation. These results suggest that Plagl1 regulates overlapping gene networks in placental trophoblast and endothelial cells, and may play a critical role in placental development in normal and complicated pregnancies.

Maternal oxycodone treatment causes pathophysiological changes in the mouse placenta

INTRODUCTION

Pregnant women are increasingly being prescribed and abusing opioid drugs. As the primary communication organ between mother and conceptus, the placenta may be vulnerable to opioid effects but also holds the key to better understanding how these drugs affect long-term offspring health. We hypothesized that maternal treatment with oxycodone (OXY), the primary opioid at the center of the current crisis, deleteriously affects placental structure and gene expression patterns.

METHODS

Female mice were treated daily with 5 mg OXY/kg or saline solution (Control, CTL) for two weeks prior to breeding and until placenta were collected at embryonic age 12.5. A portion of the placenta was fixed for histology, and the remainder was frozen for RNA isolation followed by RNAseq.

RESULTS

Maternal OXY treatment reduced parietal trophoblast giant cell (pTGC) area and decreased the maternal blood vessel area within the labyrinth region. OXY exposure affected placental gene expression profiles in a sex dependent manner with female placenta showing up-regulation of many placental enriched genes, including Ceacam11, Ceacam14, Ceacam12, Ceacam13, Prl7b1, Prl2b1, Ctsq, and Tpbpa. In contrast, placenta of OXY exposed males had alteration of many ribosomal proteins. Weighted correlation network analysis revealed that in OXY female vs. CTL female comparison, select modules correlated with OXY-induced placental histological changes. Such associations were lacking in the male OXY vs. CTL male comparison.

DISCUSSION

Results suggest OXY exposure alters placental histology. In response to OXY exposure, female placenta responds by upregulating placental enriched transcripts that are either unchanged or downregulated in male placenta. Such changes may shield female offspring from developmental origins of health and disease-based diseases.

TissueEnrich: Tissue-specific gene enrichment analysis

SUMMARY

RNA-Seq data analysis results in lists of genes that may have a similar function, based on differential gene expression analysis or co-expression network analysis. While tools have been developed to identify biological processes that are enriched in the genes sets, there remains a need for tools that identify enrichment of tissue-specific genes. Therefore, we developed TissueEnrich, a tool that calculates tissue-specific gene enrichment in an input gene set. We demonstrated that TissueEnrich can assign tissue identities to single cell clusters and differentiated embryonic stem cells.

AVAILABILITY AND IMPLEMENTATION

The TissueEnrich web application is freely available at http://tissueenrich.gdcb.iastate.edu/. The R package is available through Bioconductor at https://bioconductor.org/packages/TissueEnrich. Both the web application and R package are for non-profit academic use under the MIT license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Bisphenol A and bisphenol S disruptions of the mouse placenta and potential effects on the placenta-brain axis

Placental trophoblast cells are potentially at risk from circulating endocrine-disrupting chemicals, such as bisphenol A (BPA). To understand how BPA and the reputedly more inert bisphenol S (BPS) affect the placenta, C57BL6J mouse dams were fed 200 μg/kg body weight BPA or BPS daily for 2 wk and then bred. They continued to receive these chemicals until embryonic day 12.5, whereupon placental samples were collected and compared with unexposed controls. BPA and BPS altered the expression of an identical set of 13 genes. Both exposures led to a decrease in the area occupied by spongiotrophoblast relative to trophoblast giant cells (GCs) within the junctional zone, markedly reduced placental serotonin (5-HT) concentrations, and lowered 5-HT GC immunoreactivity. Concentrations of dopamine and 5-hydroxyindoleacetic acid, the main metabolite of serotonin, were increased. GC dopamine immunoreactivity was increased in BPA- and BPS-exposed placentas. A strong positive correlation between 5-HT+ GCs and reductions in spongiotrophoblast to GC area suggests that this neurotransmitter is essential for maintaining cells within the junctional zone. In contrast, a negative correlation existed between dopamine+ GCs and reductions in spongiotrophoblast to GC area ratio. These outcomes lead to the following conclusions. First, BPS exposure causes almost identical placental effects as BPA. Second, a major target of BPA/BPS is either spongiotrophoblast or GCs within the junctional zone. Third, imbalances in neurotransmitter-positive GCs and an observed decrease in docosahexaenoic acid and estradiol, also occurring in response to BPA/BPS exposure, likely affect the placental-brain axis of the developing mouse fetus.

Comparative Analysis of the Transcriptome and Proteome during Mouse Placental Development

The condition of the placenta is a determinant of the short- and long-term health of the mother and the fetus. However, critical processes occurring in early placental development, such as trophoblast invasion and establishment of placental metabolism, remain poorly understood. To gain a better understanding of the genes involved in regulating these processes, we utilized a multiomics approach, incorporating transcriptome, proteome, and phosphoproteome data generated from mouse placental tissue collected at two critical developmental time points. We found that incorporating information from both the transcriptome and proteome identifies genes associated with time point-specific biological processes, unlike using the proteome alone. We further inferred genes upregulated on the basis of the proteome data but not the transcriptome data at each time point, leading us to identify 27 genes that we predict to have a role in trophoblast migration or placental metabolism. Finally, using the phosphoproteome data set, we discovered novel phosphosites that may play crucial roles in the regulation of placental transcription factors. By generating the largest proteome and phosphoproteome data sets in the developing placenta, and integrating transcriptome analysis, we uncovered novel aspects of placental gene regulation.

Combined analysis of dissimilar promoter accessibility and gene expression profiles identifies tissue-specific genes and actively repressed networks

BACKGROUND

The assay for transposase-accessible chromatin (ATAC-seq) is a powerful method to examine chromatin accessibility. While many studies have reported a positive correlation between gene expression and promoter accessibility, few have investigated the genes that deviate from this trend. In this study, we aimed to understand the relationship between gene expression and promoter accessibility in multiple cell types while also identifying gene regulatory networks in the placenta, an understudied organ that is critical for a successful pregnancy.

RESULTS

We started by assaying the open chromatin landscape in the mid-gestation placenta, when the fetal vasculature has started developing. After incorporating transcriptomic data generated in the placenta at the same time point, we grouped genes based on their expression levels and ATAC-seq promoter coverage. We found that the genes with the strongest correlation (high expression and high coverage) are likely involved in housekeeping functions, whereas tissue-specific genes were highly expressed and had only medium-low coverage. We also predicted that genes with medium-low expression and high promoter coverage were actively repressed. Within this group, we extracted a protein-protein interaction network enriched for neuronal functions, likely preventing the cells from adopting a neuronal fate. We further confirmed that a repressive histone mark is bound to the promoters of genes in this network. Finally, we ran our pipeline using ATAC-seq and RNA-seq data generated in ten additional cell types. We again found that genes with the strongest correlation are enriched for housekeeping functions and that genes with medium-low promoter coverage and high expression are more likely to be tissue-specific. These results demonstrate that only two data types, both of which require relatively low starting material to generate and are becoming more commonly available, can be integrated to understand multiple aspects of gene regulation.

CONCLUSIONS

Within the placenta, we identified an active placenta-specific gene network as well as a repressed neuronal network. Beyond the placenta, we demonstrate that ATAC-seq data and RNA-seq data can be integrated to identify tissue-specific genes and actively repressed gene networks in multiple cell types.

Deciphering transcriptional regulation in human embryonic stem cells specified towards a trophoblast fate

Differentiated human embryonic stem cells (hESC) continue to provide a model for studying early trophoblast cells (TB), but many questions have been raised regarding their true identity. Therefore, we carried out a global and unbiased analysis on previously published transcriptomic profiles for hESC differentiated to TB by means of bone morphogenetic protein-4 and inhibitors of activin A and fibroblast growth factor-2 signaling (BAP treatment). Our results confirm that BAP treated hESC (ESCd) lack a mesoderm signature and are a subtype of placental cells unlike those present at term. ESCd display a high level of expression of genes implicated in migration and invasion compared to commonly used, immortalized TB cell lines and primary cells from term placenta. Co-expression network analysis also identified gene modules involved in cell migration and adhesion, processes that are likely critical during the beginning stages of placentation. Finally, protein-protein interaction analysis predicted several additional genes that may play important roles in early stages of placental development. Together, our analyses provide novel insights into the transcriptional programs that are active in ESCd.

IFPA meeting 2016 workshop report I: Genomic communication, bioinformatics, trophoblast biology and transport systems

Workshops are an important part of the IFPA annual meeting as they allow for discussion of specialized topics. At IFPA meeting 2016 there were twelve themed workshops, four of which are summarized in this report. These workshops covered innovative technologies applied to new and traditional areas of placental research: 1) genomic communication; 2) bioinformatics; 3) trophoblast biology and pathology; 4) placental transport systems.

Changes in the enhancer landscape during early placental development uncover a trophoblast invasion gene-enhancer network

INTRODUCTION

Trophoblast invasion establishes adequate blood flow between mother and fetus in early placental development. However, little is known about the cis-regulatory mechanisms underlying this important process. We aimed to identify enhancer elements that are active during trophoblast invasion, and build a trophoblast invasion gene-enhancer network.

METHODS

We carried out ChIP-Seq for an enhancer-associated mark (H3k27Ac) at two time points during early placental development in mouse. One time point when invasion is at its peak (e7.5) and another time point shortly afterwards (e9.5). We use computational analysis to identify putative enhancers, as well as the transcription factor binding sites within them, that are specific to the time point of trophoblast invasion.

RESULTS

We compared read profiles at e7.5 and e9.5 to identify 1,977 e7.5-specific enhancers. Within a subset of e7.5-specific enhancers, we discovered a cell migration associated regulatory code, consisting of three transcription factor motifs: AP1, Ets, and Tcfap2. To validate differential expression of the transcription factors that bind these motifs, we performed RNA-Seq in the same context. Finally, we integrated these data with publicly available protein-protein interaction data and constructed a trophoblast invasion gene-enhancer network.

DISCUSSION

The data we generated and analysis we carried out improves our understanding of the regulatory mechanisms of trophoblast invasion, by suggesting a transcriptional code exists in the enhancers of cell migration genes. Furthermore, the network we constructed highlights novel candidate genes that may be critical for trophoblast invasion.

Automated discovery of tissue-targeting enhancers and transcription factors from binding motif and gene function data

Identifying enhancers regulating gene expression remains an important and challenging task. While recent sequencing-based methods provide epigenomic characteristics that correlate well with enhancer activity, it remains onerous to comprehensively identify all enhancers across development. Here we introduce a computational framework to identify tissue-specific enhancers evolving under purifying selection. First, we incorporate high-confidence binding site predictions with target gene functional enrichment analysis to identify transcription factors (TFs) likely functioning in a particular context. We then search the genome for clusters of binding sites for these TFs, overcoming previous constraints associated with biased manual curation of TFs or enhancers. Applying our method to the placenta, we find 33 known and implicate 17 novel TFs in placental function, and discover 2,216 putative placenta enhancers. Using luciferase reporter assays, 31/36 (86%) tested candidates drive activity in placental cells. Our predictions agree well with recent epigenomic data in human and mouse, yet over half our loci, including 7/8 (87%) tested regions, are novel. Finally, we establish that our method is generalizable by applying it to 5 additional tissues: heart, pancreas, blood vessel, bone marrow, and liver.

Enhancers are distal gene regulatory elements that can activate tissue- and time-point specific gene expression. Identification of active enhancers is challenging, and is the subject of intense investigation. We developed an automated computational framework to predict transcription factors (TFs) and enhancers that target a tissue of interest by combining two growing resources: TF binding motifs and target gene function annotations. We applied our framework to the placenta, and confirmed our enhancer predictions are more active in placental cell types than others. To demonstrate generalizability, we applied our approach to 5 additional tissues. The combination of experimental sampling with computational prediction approaches will aid in the identification of those enhancers that are most likely active in a particular tissue, as well as the characterization of groups of TFs associated with these enhancers.

The enhancer landscape during early neocortical development reveals patterns of dense regulation and co-option

Genetic studies have identified a core set of transcription factors and target genes that control the development of the neocortex, the region of the human brain responsible for higher cognition. The specific regulatory interactions between these factors, many key upstream and downstream genes, and the enhancers that mediate all these interactions remain mostly uncharacterized. We perform p300 ChIP-seq to identify over 6,600 candidate enhancers active in the dorsal cerebral wall of embryonic day 14.5 (E14.5) mice. Over 95% of the peaks we measure are conserved to human. Eight of ten (80%) candidates tested using mouse transgenesis drive activity in restricted laminar patterns within the neocortex. GREAT based computational analysis reveals highly significant correlation with genes expressed at E14.5 in key areas for neocortex development, and allows the grouping of enhancers by known biological functions and pathways for further studies. We find that multiple genes are flanked by dozens of candidate enhancers each, including well-known key neocortical genes as well as suspected and novel genes. Nearly a quarter of our candidate enhancers are conserved well beyond mammals. Human and zebrafish regions orthologous to our candidate enhancers are shown to most often function in other aspects of central nervous system development. Finally, we find strong evidence that specific interspersed repeat families have contributed potentially key developmental enhancers via co-option. Our analysis expands the methodologies available for extracting the richness of information found in genome-wide functional maps.

Foxa2 and H2A.Z mediate nucleosome depletion during embryonic stem cell differentiation

Nucleosome occupancy is fundamental for establishing chromatin architecture. However, little is known about the relationship between nucleosome dynamics and initial cell lineage specification. Here, we determine the mechanisms that control global nucleosome dynamics during embryonic stem (ES) cell differentiation into endoderm. Both nucleosome depletion and de novo occupation occur during the differentiation process, with higher overall nucleosome density after differentiation. The variant histone H2A.Z and the winged helix transcription factor Foxa2 both act to regulate nucleosome depletion and gene activation, thus promoting ES cell differentiation, whereas DNA methylation promotes nucleosome occupation and suppresses gene expression. Nucleosome depletion during ES cell differentiation is dependent on Nap1l1-coupled SWI/SNF and INO80 chromatin remodeling complexes. Thus, both epigenetic and genetic regulators cooperate to control nucleosome dynamics during ES cell fate decisions.

Abstract 5112: Foxa1 and Foxa2 are essential for gender dimorphism in liver cancer

Hepatocellular carcinoma (HCC) is sexually dimorphic in both rodents and humans, with significantly higher incidence in males, an effect that is dependent on sex hormones. The molecular mechanisms by which estrogens prevent and androgens promote liver cancer remain unclear. Here, we discover that sexually-dimorphic HCC is completely reversed in Foxa1- and Foxa2-deficient mice after diethylnitrosamine-induced hepatocarcinogenesis. Co-regulation of target genes by Foxa1/a2 and either the estrogen receptor (ERalpha) or the androgen receptor (AR) was increased during hepatocarcinogenesis in normal female or male mice, respectively, but was lost in Foxa1/2-deficient mice. Thus, both estrogen-dependent resistance to and androgen-mediated facilitation of HCC depend on Foxa1/2. Strikingly, single nucleotide polymorphisms at FOXA2 binding sites reduce binding of both FOXA2 and ERalpha to their targets in human liver, and correlate with HCC development in women. Thus, Foxa factors and their targets are central for the sexual dimorphism of HCC. (Cell, accepted in press)

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 5112. doi:1538-7445.AM2012-5112

Species-specific strategies underlying conserved functions of metabolic transcription factors

The winged helix protein FOXA2 and the nuclear receptor peroxisome proliferator-activated receptor-γ (PPARγ) are highly conserved, regionally expressed transcription factors (TFs) that regulate networks of genes controlling complex metabolic functions. Cistrome analysis for Foxa2 in mouse liver and PPARγ in mouse adipocytes has previously produced consensus-binding sites that are nearly identical to those used by the corresponding TFs in human cells. We report here that, despite the conservation of the canonical binding motif, the great majority of binding regions for FOXA2 in human liver and for PPARγ in human adipocytes are not in the orthologous locations corresponding to the mouse genome, and vice versa. Of note, TF binding can be absent in one species despite sequence conservation, including motifs that do support binding in the other species, demonstrating a major limitation of in silico binding site prediction. Whereas only approximately 10% of binding sites are conserved, gene-centric analysis reveals that about 50% of genes with nearby TF occupancy are shared across species for both hepatic FOXA2 and adipocyte PPARγ. Remarkably, for both TFs, many of the shared genes function in tissue-specific metabolic pathways, whereas species-unique genes fail to show enrichment for these pathways. Nonetheless, the species-unique genes, like the shared genes, showed the expected transcriptional regulation by the TFs in loss-of-function experiments. Thus, species-specific strategies underlie the biological functions of metabolic TFs that are highly conserved across mammalian species. Analysis of factor binding in multiple species may be necessary to distinguish apparent species-unique noise and reveal functionally relevant information.

Extracting transcription factor targets from ChIP-Seq data

ChIP-Seq technology, which combines chromatin immunoprecipitation (ChIP) with massively parallel sequencing, is rapidly replacing ChIP-on-chip for the genome-wide identification of transcription factor binding events. Identifying bound regions from the large number of sequence tags produced by ChIP-Seq is a challenging task. Here, we present GLITR (GLobal Identifier of Target Regions), which accurately identifies enriched regions in target data by calculating a fold-change based on random samples of control (input chromatin) data. GLITR uses a classification method to identify regions in ChIP data that have a peak height and fold-change which do not resemble regions in an input sample. We compare GLITR to several recent methods and show that GLITR has improved sensitivity for identifying bound regions closely matching the consensus sequence of a given transcription factor, and can detect bona fide transcription factor targets missed by other programs. We also use GLITR to address the issue of sequencing depth, and show that sequencing biological replicates identifies far more binding regions than re-sequencing the same sample.

Foxa1 and Foxa2 regulate bile duct development in mice

The forkhead box proteins A1 and A2 (Foxa1 and Foxa2) are transcription factors with critical roles in establishing the developmental competence of the foregut endoderm and in initiating liver specification. Using conditional gene ablation during a later phase of liver development, we show here that deletion of both Foxa1 and Foxa2 (Foxa1/2) in the embryonic liver caused hyperplasia of the biliary tree. Abnormal bile duct formation in Foxa1/2-deficient liver was due, at least in part, to activation of IL-6 expression, a proliferative signal for cholangiocytes. The glucocorticoid receptor is a negative regulator of IL-6 transcription; in the absence of Foxa1/2, the glucocorticoid receptor failed to bind to the IL-6 promoter, causing enhanced IL-6 expression. Thus, after liver specification, Foxa1/2 are required for normal bile duct development through prevention of excess cholangiocyte proliferation. Our data suggest that Foxa1/2 function as terminators of bile duct expansion in the adult liver through inhibition of IL-6 expression.

Cis-regulatory modules in the mammalian liver: composition depends on strength of Foxa2 consensus site

Foxa2 is a critical transcription factor that controls liver development and plays an important role in hepatic gluconeogensis in adult mice. Here, we use genome-wide location analysis for Foxa2 to identify its targets in the adult liver. We then show by computational analyses that Foxa2 containing cis-regulatory modules are not constructed from a random assortment of binding sites for other transcription factors expressed in the liver, but rather that their composition depends on the strength of the Foxa2 consensus site present. Genes containing a cis-regulatory module with a medium or weak Foxa2 consensus site are much more liver-specific than the genes with a strong consensus site. We not only provide a better understanding of the mechanisms of Foxa2 regulation but also introduce a novel method for identification of different cis-regulatory modules involving a single factor.