Integrated RNA-seq and DNase-seq analyses identify phenotype-specific BMP4 signaling in breast cancer

Cellular zinc status alters chromatin accessibility and binding of p53 to DNA

Zn2+ is an essential metal required by approximately 850 human transcription factors. How these proteins acquire their essential Zn2+ cofactor and whether they are sensitive to changes in the labile Zn2+ pool in cells remain open questions. Using ATAC-seq to profile regions of accessible chromatin coupled with transcription factor enrichment analysis, we examined how increases and decreases in the labile zinc pool affect chromatin accessibility and transcription factor enrichment. We found 685 transcription factor motifs were differentially enriched, corresponding to 507 unique transcription factors. The pattern of perturbation and the types of transcription factors were notably different at promoters versus intergenic regions, with zinc-finger transcription factors strongly enriched in intergenic regions in elevated Zn2+ To test whether ATAC-seq and transcription factor enrichment analysis predictions correlate with changes in transcription factor binding, we used ChIP-qPCR to profile six p53 binding sites. We found that for five of the six targets, p53 binding correlates with the local accessibility determined by ATAC-seq. These results demonstrate that changes in labile zinc alter chromatin accessibility and transcription factor binding to DNA.

Enhancer selectivity across cell types delineates three functionally distinct enhancer-promoter regulation patterns

BACKGROUND

Multiple enhancers co-regulating the same gene is prevalent and plays a crucial role during development and disease. However, how multiple enhancers coordinate the same gene expression across various cell types remains largely unexplored at genome scale.

RESULTS

We develop a computational approach that enables the quantitative assessment of enhancer specificity and selectivity across diverse cell types, leveraging enhancer-promoter (E-P) interactions data. We observe two well-known gene regulation patterns controlled by enhancer clusters, which regulate the same gene either in a limited number of cell types (Specific pattern, Spe) or in the majority of cell types (Conserved pattern, Con), both of which are enriched for super-enhancers (SEs). We identify a previously overlooked pattern (Variable pattern, Var) that multiple enhancers link to the same gene, but rarely coexist in the same cell type. These three patterns control the genes associating with distinct biological function and exhibit unique epigenetic features. Specifically, we discover a subset of Var patterns contains Shared enhancers with stable enhancer-promoter interactions in the majority of cell types, which might contribute to maintaining gene expression by recruiting abundant CTCF.

CONCLUSIONS

Together, our findings reveal three distinct E-P regulation patterns across different cell types, providing insights into deciphering the complexity of gene transcriptional regulation.

Chromatin Profiles Are Prognostic of Clinical Response to Bortezomib-Containing Chemotherapy in Pediatric Acute Myeloid Leukemia: Results from the COG AAML1031 Trial

The addition of the proteasome inhibitor bortezomib to standard chemotherapy did not improve survival in pediatric acute myeloid leukemia (AML) when all patients were analyzed as a group in the Children's Oncology Group phase 3 trial AAML1031 (NCT01371981). Proteasome inhibition influences the chromatin landscape and proteostasis, and we hypothesized that baseline proteomic analysis of histone- and chromatin-modifying enzymes (HMEs) would identify AML subgroups that benefitted from bortezomib addition. A proteomic profile of 483 patients treated with AAML1031 chemotherapy was generated using a reverse-phase protein array. A relatively high expression of 16 HME was associated with lower EFS and higher 3-year relapse risk after AML standard treatment compared to low expressions (52% vs. 29%, p = 0.005). The high-HME profile correlated with more transposase-accessible chromatin, as demonstrated via ATAC-sequencing, and the bortezomib addition improved the 3-year overall survival compared with standard therapy (62% vs. 75%, p = 0.033). These data suggest that there are pediatric AML populations that respond well to bortezomib-containing chemotherapy.

Chromatin Accessibility Profiling and Data Analysis Using ATAC-seq in Nothobranchius furzeri

The state of genome-wide chromatin accessibility in cells, tissues, or organisms can be investigated with a technique called assay for transposase-accessible chromatin using sequencing (ATAC-seq). ATAC-seq is a powerful approach for profiling the epigenomic landscape of cells using very low input materials. Analysis of chromatin accessibility data allows for prediction of gene expression and identification of regulatory elements such as potential enhancers and specific transcription-factor binding sites. Here, we describe an optimized ATAC-seq protocol for the preparation of isolated nuclei and subsequent next-generation sequencing from whole embryos and tissues of the African turquoise killifish (Nothobranchius furzeri). Importantly, we provide an overview of a pipeline for processing and analyzing ATAC-seq data from killifish.

Cellular zinc status alters chromatin accessibility and binding of transcription factor p53 to genomic sites

Zinc (Zn2+) is an essential metal required by approximately 2500 proteins. Nearly half of these proteins act on DNA, including > 850 human transcription factors, polymerases, DNA damage response factors, and proteins involved in chromatin architecture. How these proteins acquire their essential Zn2+ cofactor and whether they are sensitive to changes in the labile Zn2+ pool in cells remain open questions. Here, we examine how changes in the labile Zn2+ pool affect chromatin accessibility and transcription factor binding to DNA. We observed both increases and decreases in accessibility in different chromatin regions via ATAC-seq upon treating MCF10A cells with elevated Zn2+ or the Zn2+-specific chelator tris(2-pyridylmethyl)amine (TPA). Transcription factor enrichment analysis was used to correlate changes in chromatin accessibility with transcription factor motifs, revealing 477 transcription factor motifs that were differentially enriched upon Zn2+ perturbation. 186 of these transcription factor motifs were enriched in Zn2+ and depleted in TPA, and the majority correspond to Zn2+ finger transcription factors. We selected TP53 as a candidate to examine how changes in motif enrichment correlate with changes in transcription factor occupancy by ChIP-qPCR. Using publicly available ChIP-seq and nascent transcription datasets, we narrowed the 50,000+ ATAC-seq peaks to 2164 TP53 targets and subsequently selected 6 high-probability TP53 binding sites for testing. ChIP-qPCR revealed that for 5 of the 6 targets, TP53 binding correlates with the local accessibility determined by ATAC-seq. These results demonstrate that changes in labile zinc directly alter chromatin accessibility and transcription factor binding to DNA.

Single nucleus multiomics identifies ZEB1 and MAFB as candidate regulators of Alzheimer's disease-specific cis-regulatory elements

Cell type-specific transcriptional differences between brain tissues from donors with Alzheimer's disease (AD) and unaffected controls have been well documented, but few studies have rigorously interrogated the regulatory mechanisms responsible for these alterations. We performed single nucleus multiomics (snRNA-seq plus snATAC-seq) on 105,332 nuclei isolated from cortical tissues from 7 AD and 8 unaffected donors to identify candidate cis-regulatory elements (CREs) involved in AD-associated transcriptional changes. We detected 319,861 significant correlations, or links, between gene expression and cell type-specific transposase accessible regions enriched for active CREs. Among these, 40,831 were unique to AD tissues. Validation experiments confirmed the activity of many regions, including several candidate regulators of APP expression. We identified ZEB1 and MAFB as candidate transcription factors playing important roles in AD-specific gene regulation in neurons and microglia, respectively. Microglia links were globally enriched for heritability of AD risk and previously identified active regulatory regions.

Towards the Improvement of Ornamental Attributes in Chrysanthemum: Recent Progress in Biotechnological Advances

Incessant development and introduction of novel cultivars with improved floral attributes are vital in the dynamic ornamental industry. Chrysanthemum (Chrysanthemum morifolium) is a highly favored ornamental plant, ranking second globally in the cut flower trade, after rose. Development of new chrysanthemum cultivars with improved and innovative modifications in ornamental attributes, including floral color, shape, plant architecture, flowering time, enhanced shelf life, and biotic and abiotic stress tolerance, is a major goal in chrysanthemum breeding. Despite being an economically important ornamental plant, the application of conventional and molecular breeding approaches to various key traits of chrysanthemum is hindered owing to its genomic complexity, heterozygosity, and limited gene pool availability. Although classical breeding of chrysanthemum has resulted in the development of several hundreds of cultivars with various morphological variations, the genetic and transcriptional control of various important ornamental traits remains unclear. The coveted blue colored flowers of chrysanthemums cannot be achieved through conventional breeding and mutation breeding due to technical limitations. However, blue-hued flower has been developed by genetic engineering, and transgenic molecular breeding has been successfully employed, leading to substantial progress in improving various traits. The recent availability of whole-genome sequences of chrysanthemum offers a platform to extensively employ MAS to identify a large number of markers for QTL mapping, and GWAS to dissect the genetic control of complex traits. The combination of NGS, multi-omic platforms, and genome editing technologies has provided a tremendous scope to decipher the molecular and regulatory mechanisms. However, the application and integration of these technologies remain inadequate for chrysanthemum. This review, therefore, details the significance of floral attributes, describes the efforts of recent advancements, and highlights the possibilities for future application towards the improvement of crucial ornamental traits in the globally popular chrysanthemum plant.

Changes in chromatin accessibility are not concordant with transcriptional changes for single-factor perturbations

A major goal in the field of transcriptional regulation is the mapping of changes in the binding of transcription factors to the resultant changes in gene expression. Recently, methods for measuring chromatin accessibility have enabled us to measure changes in accessibility across the genome, which are thought to correspond to transcription factor‐binding events. In concert with RNA‐sequencing, these data in principle enable such mappings; however, few studies have looked at their concordance over short‐duration treatments with specific perturbations. Here, we used tandem, bulk ATAC‐seq, and RNA‐seq measurements from MCF‐7 breast carcinoma cells to systematically evaluate the concordance between changes in accessibility and changes in expression in response to retinoic acid and TGF‐β. We found two classes of genes whose expression showed a significant change: those that showed some changes in the accessibility of nearby chromatin, and those that showed virtually no change despite strong changes in expression. The peaks associated with genes in the former group had lower baseline accessibility prior to exposure to signal. Focusing the analysis specifically on peaks with motifs for transcription factors associated with retinoic acid and TGF‐β signaling did not reduce the lack of correspondence. Analysis of paired chromatin accessibility and gene expression data from distinct paths along the hematopoietic differentiation trajectory showed a much stronger correspondence, suggesting that the multifactorial biological processes associated with differentiation may lead to changes in chromatin accessibility that reflect rather than driving altered transcriptional status. Together, these results show many gene expression changes can happen independently of changes in the accessibility of local chromatin in the context of a single‐factor perturbation.

Increased Methyl-CpG-Binding Domain Protein 2 Promotes Cigarette Smoke-Induced Pulmonary Hypertension

Pulmonary hypertension (PH) is a chronic vascular proliferative disorder. While cigarette smoke (CS) plays a vital part in PH related to chronic obstructive pulmonary disease (COPD). Methyl-CpG-Binding Domain Protein 2 (MBD2) has been linked to multiple proliferative diseases. However, the specific mechanisms of MBD2 in CS-induced PH remain to be elucidated. Herein, the differential expression of MBD2 was tested between the controls and the PH patients’ pulmonary arteries, CS-exposed rat models’ pulmonary arteries, and primary human pulmonary artery smooth muscle cells (HPASMCs) following cigarette smoke extract (CSE) stimulation. As a result, PH patients and CS-induced rats and HPASMCs showed an increase in MBD2 protein expression compared with the controls. Then, MBD2 silencing was used to investigate the function of MBD2 on CSE-induced HPASMCs’ proliferation, migration, and cell cycle progression. As a consequence, CSE could induce HPASMCs’ increased proliferation and migration, and cell cycle transition, which were suppressed by MBD2 interference. Furthermore, RNA-seq, ChIP-qPCR, and MassARRAY were conducted to find out the downstream mechanisms of MBD2 for CS-induced pulmonary vascular remodeling. Subsequently, RNA-seq revealed MBD2 might affect the transcription of BMP2 gene, which furtherly altered the expression of BMP2 protein. ChIP-qPCR demonstrated MBD2 could bind BMP2’s promotor. MassARRAY indicated that MBD2 itself could not directly affect DNA methylation. In sum, our results indicate that increased MBD2 expression promotes CS-induced pulmonary vascular remodeling. The fundamental mechanisms may be that MBD2 can bind BMP2’s promoter and downregulate its expression. Thus, MBD2 may promote the occurrence of the CS-induced PH.

The Key Network of mRNAs and miRNAs Regulated by HIF1A in Hypoxic Hepatocellular Carcinoma Cells

Purpose: Hypoxia plays an essential role in the progression of hepatocellular carcinoma (HCC), whereas hypoxia inducible factor-1 (HIF-1) is the key transcription factor allowing HCC to survive hypoxia. The aim of this study was to define the essential mRNAs and miRNAs regulated by HIF1A and dissect their functions, interactions, and tumor-infiltrating immune cells in HCC.

Methods: A human HCC cell line HepG2 was used as a cell model of HCC. The CRISPR/Cas9 system was used to knock out HIF1A in HepG2 cells, and RNA sequencing was utilized to characterize differentially expressed mRNAs and miRNAs in the HIF1A-knockout HepG2 cells; the identified candidates were then analyzed by GO annotation and KEGG pathway enrichment to study their function and establish a PPI network. Quantitative (q) PCR was used to verify if there were significant differences in the expression of mRNAs, and the association of the selected mRNAs expression with immune cell infiltration levels was further analyzed using The Cancer Genome Atlas (TCGA) pan-cancer data.

Results: Using RNA-sequencing, we discovered that there were 1535 mRNAs differentially expressed (adjusted p < 0.05, |fold change|>1.5) in the HIF1A-knockout HepG2 cells, among which there were 644 mRNAs upregulated and 891 mRNAs downregulated. GO annotation and KEGG pathway enrichment showed that these mRNAs were involved in glycolysis/gluconeogenesis, PI3K-Akt signaling pathways, and HIF-1 signaling pathways. In addition, we found that there were 309 miRNAs differentially expressed (adjusted p < 0.05, |fold change|>1.5) in the HIF1A-knockout HepG2 cells, of which there were 213 miRNAs upregulated and 96 miRNAs downregulated. Our further analyses uncovered that these miRNA putative targets were involved in the hippo signaling pathway, axon guidance, and tight junction. Moreover, the construction and analysis of the PPI network showed that OASL, IL6, and TAF1 were recognized as hub genes with the highest connectivity degrees. Importantly, in the HIF1A-knockout HepG2 cells, our qRT-PCR data confirmed the selected mRNA changes revealed by RNA-sequencing, and with TCGA pan-cancer data, we revealed that the expressional levels of these three genes, LUM, SCOC, and CCL2, were associated with immune cell infiltration levels.

Conclusion: The identified potential key network of mRNAs and miRNAs regulated by HIF1A in the HCC cells suggests a key role of HIF1A in the tumorigenesis of HCC.

ATAC-Seq Reveals the Landscape of Open Chromatin and cis-Regulatory Elements in the Phytophthora sojae Genome

Nucleosome-free open chromatin often harbors transcription factor (TF)-binding sites that are associated with active cis-regulatory elements. However, analysis of open chromatin regions has rarely been applied to oomycete or fungal plant pathogens. In this study, we performed the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) to identify open chromatin and cis-regulatory elements in Phytophthora sojae at the mycelial stage. We identified 10,389 peaks representing nucleosome-free regions (NFRs). The peaks were enriched in gene-promoter regions and associated with 40% of P. sojae genes; transcription levels were higher for genes with multiple peaks than genes with a single peak and were higher for genes with a single peak than genes without peak. Chromatin accessibility was positively correlated with gene transcription level. Through motif discovery based on NFR peaks in core promoter regions, 25 candidate cis-regulatory motifs with evidence of TF-binding footprints were identified. These motifs exhibited various preferences for location in the promoter region and associations with the transcription level of their target genes, which included some putative pathogenicity-related genes. As the first study revealing the landscape of open chromatin and the correlation between chromatin accessibility and gene transcription level in oomycetes, the results provide a technical reference and data resources for future studies on the regulatory mechanisms of gene transcription.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Differential analysis of chromatin accessibility and gene expression profiles identifies cis-regulatory elements in rat adipose and muscle

Chromatin accessibility is a key factor influencing gene expression. We optimized the Omni-ATAC-seq protocol and used it together with RNA-seq to investigate cis-regulatory elements in rat white adipose and skeletal muscle, two tissues with contrasting metabolic functions. While promoter accessibility correlated with RNA expression, integration of the two datasets identified tissue-specific differentially accessible regions (DARs) that predominantly localized in intergenic and intron regions. DARs were mapped to differentially expressed (DE) genes enriched in distinct biological processes in each tissue. Randomly selected DE genes were validated by qPCR. Top enriched motifs in DARs predicted binding sites for transcription factors (TFs) showing tissue-specific up-regulation. The correlation between differential chromatin accessibility at a given TF binding motif and differential expression of target genes further supported the functional relevance of that motif. Our study identified cis-regulatory regions that likely play a major role in the regulation of tissue-specific gene expression in adipose and muscle.

A long noncoding RNA, LncMyoD, modulates chromatin accessibility to regulate muscle stem cell myogenic lineage progression

Epigenetic regulations control the accessibility of transcription factors to their target regions. Modulation of chromatin accessibility determines which transcripts to be expressed and therefore, defines cell identity. Chromatin modulation during cell fate determination involves a complex regulatory network, yet the comprehensive view remains to be explored. Here, we provide a global view of chromatin accessibility during muscle stem cell activation. We identified a long noncoding RNA (lncRNA), LncMyoD, which regulates lineage determination and progression through modulating chromatin accessibility. Functional analysis showed that loss of LncMyoD strongly impairs reprogramming of fibroblasts into myogenic lineage and causes defects in muscle stem cell differentiation. Our findings provide an epigenetic mechanism for the regulation of muscle stem cell myogenic lineage progression by an lncRNA.

Epigenetics regulation plays a critical role in determining cell identity by controlling the accessibility of lineage-specific regulatory regions. In muscle stem cells, epigenetic mechanisms of how chromatin accessibility is modulated during cell fate determination are not fully understood. Here, we identified a long noncoding RNA, LncMyoD, that functions as a chromatin modulator for myogenic lineage determination and progression. The depletion of LncMyoD in muscle stem cells led to the down-regulation of myogenic genes and defects in myogenic differentiation. LncMyoD exclusively binds with MyoD and not with other myogenic regulatory factors and promotes transactivation of target genes. The mechanistic study revealed that loss of LncMyoD prevents the establishment of a permissive chromatin environment at myogenic E-box–containing regions, therefore restricting the binding of MyoD. Furthermore, the depletion of LncMyoD strongly impairs the reprogramming of fibroblasts into the myogenic lineage. Taken together, our study shows that LncMyoD associates with MyoD and promotes myogenic gene expression through modulating MyoD accessibility to chromatin, thereby regulating myogenic lineage determination and progression.

Genomic Rewiring of SOX2 Chromatin Interaction Network during Differentiation of ESCs to Postmitotic Neurons

Cellular differentiation requires dramatic changes in chromatin organization, transcriptional regulation, and protein production. To understand the regulatory connections between these processes, we generated proteomic, transcriptomic, and chromatin accessibility data during differentiation of mouse embryonic stem cells (ESCs) into postmitotic neurons and found extensive associations between different molecular layers within and across differentiation time points. We observed that SOX2, as a regulator of pluripotency and neuronal genes, redistributes from pluripotency enhancers to neuronal promoters during differentiation, likely driven by changes in its protein interaction network. We identified ATRX as a major SOX2 partner in neurons, whose co-localization correlated with an increase in active enhancer marks and increased expression of nearby genes, which we experimentally confirmed for three loci. Collectively, our data provide key insights into the regulatory transformation of SOX2 during neuronal differentiation, and we highlight the significance of multi-omic approaches in understanding gene regulation in complex systems.

•

Complex interplay of RNA, protein, and chromatin during neuronal differentiation

•

Multi-omic profiling reveals divergent roles of SOX2 in stem cells and neurons

•

SOX2 on-chromatin interaction network changes from pluripotent to neuronal factors

•

ATRX interacts with SOX2 in neurons and co-binds highly expressed neuronal genes

Cis-regulatory basis of sister cell type divergence in the vertebrate retina

Multicellular organisms evolved via repeated functional divergence of transcriptionally related sister cell types, but the mechanisms underlying sister cell type divergence are not well understood. Here, we study a canonical pair of sister cell types, retinal photoreceptors and bipolar cells, to identify the key cis-regulatory features that distinguish them. By comparing open chromatin maps and transcriptomic profiles, we found that while photoreceptor and bipolar cells have divergent transcriptomes, they share remarkably similar cis-regulatory grammars, marked by enrichment of K50 homeodomain binding sites. However, cell class-specific enhancers are distinguished by enrichment of E-box motifs in bipolar cells, and Q50 homeodomain motifs in photoreceptors. We show that converting K50 motifs to Q50 motifs represses reporter expression in bipolar cells, while photoreceptor expression is maintained. These findings suggest that partitioning of Q50 motifs within cell type-specific cis-regulatory elements was a critical step in the evolutionary divergence of the bipolar transcriptome from that of photoreceptors.

miR-26a promotes hepatocellular carcinoma invasion and metastasis by inhibiting PTEN and inhibits cell growth by repressing EZH2

A previous study revealed that therapeutic miR-26a delivery suppresses tumorigenesis in a murine liver cancer model, whereas we found that forced miR-26a expression increased hepatocellular carcinoma (HCC) cell migration and invasion, which prompted us to characterize the causes and mechanisms underlying enhanced invasion due to ectopic miR-26a expression. Gain-of-function and loss-of-function experiments demonstrated that miR-26a promoted migration and invasion of BEL-7402 and HepG2 cells in vitro and positively modulated matrix metalloproteinase (MMP)-1, MMP-2, MMP-9, and MMP-10 expression. In addition, exogenous miR-26a expression significantly enhanced the metastatic ability of HepG2 cells in vivo. miR-26a negatively regulated in vitro proliferation of HCC cells, and miR-26a overexpression suppressed HepG2 cell tumor growth in nude mice. Further studies revealed that miR-26a inhibited cell growth by repressing the methyltransferase EZH2 and promoted cell migration and invasion by inhibiting the phosphatase PTEN. Furthermore, PTEN expression negatively correlated with miR-26a expression in HCC specimens from patients with and without metastasis. Thus, our findings suggest for the first time that miR-26a promotes invasion/metastasis by inhibiting PTEN and inhibits cell proliferation by repressing EZH2 in HCC. More importantly, our data also suggest caution if miR-26a is used as a target for cancer therapy in the future.

Unveiling the gene regulatory landscape in diseases through the identification of DNase I-hypersensitive sites

DNase I-hypersensitive sites (DHSs) serve key roles in the regulation of gene transcription as markers of cis-regulatory elements (CREs). Recent advances in next-generation sequencing have enabled the genome-wide location and annotation of DHSs in a variety of cells. Numerous studies have confirmed that DHSs are involved in several processes in cell fate decision and development. DHSs have also been indicated in cancer and inherited diseases as driver distal regulatory elements. Here, the definition of DHSs is reviewed, in addition to high-throughput methods of DHS identification. Furthermore, the function of DHSs in gene expression is probed. The roles of DHSs in disease occurrence are also reviewed and discussed. Concomitant advances in the identification of essential roles of DHSs will assist in disclosing the underlying molecular mechanisms, supplementing gene transcription and enlarging the molecular basis of DHS-related bioprocesses, phenotypes, distinct traits and diseases.

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.