Next Article in Journal
Characteristics of Serum Exosomes after Burn Injury and Dermal Fibroblast Regulation by Exosomes In Vitro
Next Article in Special Issue
Long-Term Maintenance of Viable Human Endometrial Epithelial Cells to Analyze Estrogen and Progestin Effects
Previous Article in Journal
Dopamine Transporter Deficiency Syndrome (DTDS): Expanding the Clinical Phenotype and Precision Medicine Approaches
Previous Article in Special Issue
Nociceptin/Orphanin FQ Opioid Peptide-Receptor Expression in the Endometriosis-Associated Nerve Fibers—Possible Treatment Option?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Analysis of Transcription Start Sites and Enhancers in Endometrial Stromal Cells and Differences Associated with Endometriosis

1
The Institute for Molecular Bioscience, The University of Queensland, Brisbane, QLD 4072, Australia
2
The Genome Innovation Hub, The University of Queensland, Brisbane, QLD 4072, Australia
3
Department of Gynecology and Gynecological Oncology, Inselspital, Bern University Hospital, University of Bern, 3010 Berne, Switzerland
*
Author to whom correspondence should be addressed.
Cells 2023, 12(13), 1736; https://doi.org/10.3390/cells12131736
Submission received: 9 May 2023 / Revised: 20 June 2023 / Accepted: 21 June 2023 / Published: 28 June 2023
(This article belongs to the Special Issue Molecular Advances and New Therapeutic Approaches in Endometriosis)

Abstract

:
Identifying tissue-specific molecular signatures of active regulatory elements is critical to understanding gene regulatory mechanisms. In this study, transcription start sites (TSS) and enhancers were identified using Cap analysis of gene expression (CAGE) across endometrial stromal cell (ESC) samples obtained from women with (n = 4) and without endometriosis (n = 4). ESC TSSs and enhancers were compared to those reported in other tissue and cell types in FANTOM5 and were integrated with RNA-seq and ATAC-seq data from the same samples for regulatory activity and network analyses. CAGE tag count differences between women with and without endometriosis were statistically tested and tags within close proximity to genetic variants associated with endometriosis risk were identified. Over 90% of tag clusters mapping to promoters were observed in cells and tissues in FANTOM5. However, some potential cell-type-specific promoters and enhancers were also observed. Regions of open chromatin identified using ATAC-seq provided further evidence of the active transcriptional regions identified by CAGE. Despite the small sample number, there was evidence of differences associated with endometriosis at 210 consensus clusters, including IGFBP5, CALD1 and OXTR. ESC TSSs were also located within loci associated with endometriosis risk from genome-wide association studies. This study provides novel evidence of transcriptional differences in endometrial stromal cells associated with endometriosis and provides a valuable cell-type specific resource of active TSSs and enhancers in endometrial stromal cells.

Graphical Abstract

1. Introduction

Human endometrium, the inner lining of the uterus, is a vitally important reproductive tissue, essential for fertility and associated with many reproductive disorders [1]. The endometrium consists of distinct cell types: epithelial cells (luminal and glandular) and the supporting mesenchymal cells (stromal cells). The tissue undergoes cyclical changes in cellular composition and gene expression, largely driven by hormonal regulation [2,3,4]. One associated disease is endometriosis, a common gynaecological disease defined by the presence of endometrial glands and stroma outside the uterus. The commonly accepted theory for the origin of cells in lesions is viable endometrial cells being transported in retrograde menstruation [5]. A major approach to discovering the fundamental mechanisms and causes of endometriosis is through identification and understanding of changes in transcriptional regulation associated with disease and the functions of the genetic variants responsible for the genetic component of the disease. Although genome-wide association study (GWAS) results [6] represent a significant breakthrough towards understanding genetic risk factors for endometriosis, the cellular mechanisms through which these variants drive the disease remain largely unknown.
Understanding gene expression in endometriosis [3] and endometrial cancer, especially in the context of fertility-sparing treatment [7], requires analysis of the transcriptional landscape of several uterine cell types. This approach may help identify disease processes that are unique to certain cell types but would otherwise be averaged out in bulk tissue datasets. One study [8] reported that the development of endometriosis may be influenced by an epigenetic switch that regulates the expression of the GATA isoform in stromal cells. Gene transcription is controlled by functional interactions between promoters and enhancers. Identifying active regulatory regions in the genome is critical for understanding gene regulation and assessing the impact of genetic variation on phenotype. The Cap analysis of gene expression (CAGE) method is a high-throughput, transcriptome-wide and unbiased (towards polyA-tailed RNA) approach, well established at mapping the transcription start sites (TSSs) by capturing the 5′ ends of the transcribed and capped mRNAs [9]. CAGE employs a cap-trapping method and captures 5′ ends of cDNAs, which yields short sequences (tags) of 20 nucleotides that can be mapped back to the reference genome to infer the exact position of the TSSs of captured RNAs. The number of CAGE tags supporting each TSS reflects the relative frequency of its usage and can be used as a measure of expression from that specific TSS, reflecting the activities of a promoter or an enhancer [10]. The functional annotation of the mammalian genome (FANTOM) project [11] generated CAGE expression measurements across 573 primary cell types and tissues.
To date, there are no publicly available CAGE datasets for endometrium/endometrial stromal cells. In this study, for the first time, CAGE data were produced from endometrial stromal cells of endometriosis patients and controls. We built a genome-wide TSS map for stromal cells and we added a deep ATAC-seq dataset, as well as a total RNA-seq dataset, to the valuable CAGE resource, forming the most comprehensive genomics regulatory database of the endometrial stromal cells. We performed integrative analysis across existing epigenomic datasets to reveal molecular signatures of non-coding stromal cell elements and differences identified in endometriosis patients.

2. Materials and Methods

2.1. Samples

Endometriotic biopsies were collected via soft curette (Pipelle de Cornier, Laboratorie CCD, France) from eight women undergoing laparoscopic surgery for suspected endometriosis. Samples were stored in Complete IMDM media (10% foetal calf serum, 1% antibiotics/antimycotics (ThermoFisher Scientific, Waltham, MA, USA)) supplemented with 10% Dimethyl sulfoxide (DMSO) (Thermo Fischer Scientific, Waltham, MA, USA) and stored at −80 °C using the slow-freezing method in a Bicell vessel. The pelvic cavity of each patient was subsequently examined and any suspected endometriotic lesions were removed and disease stage reported according to the revised American Fertility Society staging system (rAFS) [12]. The location of any lesion identified was noted using the ENZIAN classification [13] (Table 1) and subsequently categorized as either superficial peritoneal, ovarian or deep-infiltrating endometriosis (DIE) in the base of the pelvis.
The endometrial stromal cells (ESCs) were isolated and purified according to our published procedure [14,15]. Cells were cultured/passaged as per standard methods with trypsin/EDTA when approximately 80% confluent. Cells were counted at each passage, population doubling was determined and only cells that remained within a log phase of growth were included in the analysis.

2.2. CAGE Sequencing

CAGE data were generated from eight endometriosis stromal cell lines, four of which were derived from women with endometriosis (cases) and four of them from women without endometriosis (controls). For a comprehensive analysis of transcription start sites, we performed two different CAGE protocols, each with its own advantages and disadvantages, as described below.
Commercial protocol: CAGE libraries were generated via the cap-trapping method using a commercially available kit from DNAFORM. The resulting libraries were PCR-free, full-length (un-fragmented) Illumina-compatible libraries containing both polyA and non-polyA transcripts. The entire process from RNA to final library took eight days and required 5 ug RNA/sample with RIN > 7. The kit contained most of the necessary components for the entire workflow.
In-house protocol: CAGE libraries were also generated via the cap-trapping method using an optimized protocol designed and validated in house. This protocol was developed to increase throughput, decrease processing time and offer the flexibility to run samples with standard molecular reagents. The in-house protocol shares many common steps with the commercially available kit from DNAFORM, but it is shorter (taking about five days) and the protocol has a PCR step to generate higher concentration libraries. The resulting libraries were full-length (un-fragmented), Illumina-compatible libraries containing both polyA and non-polyA transcripts that had been minimally amplified by PCR. The entire process required 5 ug RNA/sample with RIN > 7.
Sequencing, read processing and alignment: Sequencing was performed on an Illumina NextSeq500 with 5% PhiX spike-in using the following configurations: Read1-76 bp, Index1-6 bp. Loading concentrations for all the samples were 1.6 pM using a modified denaturation and dilution protocol for low-concentration libraries. Fastq files were generated from the sequencing BCL files using fastqc and raw reads were filtered according to quality using multiqc followed by trimming the reads to remove the ‘N’ base at the 5′ and 3′ ends using MORAI20140528/bin/trimBaseN package. The resulting reads were then aligned to the human genome hg38 assembly https://ftp.ensembl.org/pub/release-100/fasta/homo_sapiens/dna/Homo_sapiens.GRCh38.dna.primary_assembly.fa (accessed on 25 April 2020) using BWA v0.7.15. The uniquely mapped CAGE tags with a minimal mapping quality of 10 were used in this study. The BWA SAM-formatted alignments were converted to BAM format using samtools. The BAM files were then converted to bed files and BigWig files using bedtools2 v 2.24.0 and samtools.

2.3. Identification of CAGE-Defined Transcription Start Sites (CTSS) and Tag Clusters (TC)

The vast number of small CAGE tag sequences demands considerable data processing, and statistical approaches are required for biological inferences and discoveries. The follow-up data analysis was based primarily on the annotation of CAGE sequence tags to known transcripts/genes. CAGE data processing, data formats and analysis methodologies are not yet standardised, and the scientific community is still working to develop them. CAGE tags/sequences are aligned to genome sequences using simple computational processes (called BLAST) and then counted, which provides the frequency of RNA expression. CAGE expression is measured specifically for each transcription start site (TSS).
TSSs frequently cluster into tag clusters (TC) that represent transcriptionally active regions [9]. CTSS (CAGE-detected transcription start site) is a cluster of CAGE tags that share the same nucleotide position at their 5′-ends which are defined to reflect a high-resolution map of each individual TSS/TC (tag clusters), spanning a genomic area indicating a group of CTSS. TCs describe a larger area, serving as units for identifying potential core promoters. The shapes of TCs change in a regular manner. In mammalian cells, TCs are often smaller than 10 basepairs (Sharp TSSs) or hundreds of basepairs wide (Broad TSSs) (Figure 1).
The uniquely aligned reads from the bed files were provided as input for the CAGEr pipeline in R [16]. The 5′ end of each read was considered as the single-nucleotide resolution CTSS, and counts were generated for each position to represent the number of unique tags starting from that position. These tag counts were normalised to one million reads (TPM), and CTSS that were within a 20 bp distance were merged to generate TCs. Each TC was considered as an active transcription start site, as a functional unit of a promoter or as an enhancer. Prior to clustering, we filtered out low-fidelity TSSs, i.e., the ones supported by less than 2 normalised tag counts in all of the samples and included only singleton CTSS which had a normalised signal above three. The TCs from all the samples that were within 100 bp of each other were clustered together into a single set of non-overlapping consensus clusters to define the promoter regions.

2.4. Correlation of Tag Counts within and between Protocols

We performed Pearson’s correlation for normalised CAGE tag counts per TSS between the samples separately for commercial and in-house data. For robustness, only TSSs with at least five tag counts were considered for calculating the correlation coefficient. The ggcorrplot () function in R was used to plot correlation coefficients between all possible pairs of samples.
We obtained the mean consensus cluster tag counts for each promoter/gene from all the samples in each of the commercial and in-house datasets and performed Pearson’s correlation to calculate correlation coefficients between same-sample pairs obtained from the two protocols.

2.5. Merging Data from In-House and Commercial Protocols

To increase TC detection sensitivity, we combined CAGE data from the commercial and in-house protocols. This approach has been widely applied in integrating CAGE datasets. Fastq files from the same sample obtained from two protocols were concatenated. Alignment, mapping, CTSS and TC identification were repeated for the merged files using the method described above.

2.6. Promoter Shapes

The CAGEr pipeline was used to evaluate endometrial stromal cell promoters predicted from CAGE clusters. Broad or narrow promoter types were determined by calculating the inter-quantile range of promoter CAGE clusters by assessing the base pair distance between 10% and 90% of a promoter’s total signal. The shapes of the promoter, as determined by the depth of the mapped reads and the width of the clusters, are a unique feature associated with transcriptional regulation activities.

2.7. Enhancer Identification

Enhancers are distal regulatory elements that regulate the transcription of their target genes [17,18]. RNApolII binds to enhancer regions and eRNAs are transcribed bidirectionally from active enhancers [19]. Similar to mRNAs, eRNAs are transcribed by RNApolII and capped in the process, which can be captured by the CAGE method. Enhancers are transcribed in a cell-type-specific manner [20]. Active enhancers can be detected by CAGE signals at the two ends of the enhancers [21]. CAGEfightR [22] was used to analyse the enhancers from our CAGE data. The CTSS coordinates and counts from the CAGEr pipeline were input as bigwig files into the CAGEfightR program. CTSS appearing in only a single sample were removed. The remaining CTSS were normalised to tags per million (TPM) and summed across the samples to yield a pooled CTSS signal. CTSS within 20 bp distance were merged to obtain unidirectional clusters and were required to have 0.3 TPM in at least two samples. Bidirectional clusters were required to have a balance score ≥ 0.95 and to be bidirectional in at least a single sample with at least ≥ 2 counts in at least one sample. Enhancers identified were then annotated to genomic locations (TxDb.Hsapiens.UCSC.hg38.knownGene: Annotation package for TxDb object(s). R package version 4.1.3). Enhancers annotated to introns and intergenic regions were used in subsequent analyses.

2.8. Interaction of Transcription Start Sites and Enhancers

TSS-enhancer candidates were identified using the findLinks function from the InteractionSet R package by looking for very closely spaced TSSs and enhancers that had a highly correlated expression within 50 bp distance. This analysis linked promoters and enhancers which were likely to regulate similar genes/pathways.

2.9. Differential Promoter and Enhancer Expression

Differential expression analyses were performed using the Bioconductor package DESeq2 [23] in R to build a statistical model and perform contrasts of endometriosis cases versus controls for consensus cluster tag counts and enhancer counts. Fold changes were obtained along with their associated p-value. The Benjamini–Hochberg method was used to control the false-discovery rate (FDR) by adjusting p-values for multiple testing correction. A promoter or enhancer was defined as significantly differentially expressed if it had a Benjamini–Hochberg-adjusted (FDR) p-adjusted < 0.05.

2.10. RNA Extraction and Sequencing

The same eight cell lines were subjected to RNA sequencing. The AllprepRNA Mini kit was used to extract total RNA from endometrial stromal cells that were actively growing. Before collecting the mixture in 1.5 mL Eppendorf tubes, Qiagen lysis buffer was added to the cells directly. RNA was then isolated according to the manufacturer′s instructions (Qiagen, Redwood City, CA, USA). RNA was treated with a Turbo DNA-free kit (Thermofisher Scientific, USA) and the RNA concentrations and integrity of each sample were evaluated using the NanoDRopBD-6000 instrument (ThermoFisher Scientific, Waltham, MA, USA). For RNA sequencing, only samples with an RNA integrity number of >8 were used. Stranded RNA-seq libraries were generated using the Illumina Tru-Seq Stranded Total RNA Gold technique, including ribosomal depletion, per the manufacturer′s instructions (Illumina, San Diego, CA, USA). The Illumina HISeq 4000 was used to sequence libraries that were pooled. Low-quality sequencing reads and contaminated HiSeq Illuminata adaptor sequences were trimmed using Trimmomatic v0.36 [24].

2.11. Correlation of CAGE-Seq Data with RNA-Seq Data

The average tag counts annotated to the promoter of each gene were correlated to count estimates from RNA-seq data in the same samples to evaluate concordance in gene expression intensities between RNA-seq and CAGE. Pearson’s correlation was used to estimate the correlation between the log-transformed count data of CAGE and RNA-seq.

2.12. ATAC-Seq Using the Omni-ATAC Protocol from Actively Growing Endometrial Stromal Cells

Sample preparation: ATAC-seq was performed on the same eight cell lines. Cells growing in the tissue culture were Trypsinized and the medium was then washed out. The cells were resuspended in cold PBS and counted. Cell pellets of 50,000 cells were resuspended in 50 µL of ice-cold ATAC-Resuspension Buffer (0.1% NP40, 0.1% Tween-20 and 0.01% Digitonin) by pipetting up and down three times. This cell lysis reaction was incubated on ice for 3 min. After lysis, 1 mL of ATAC-seq RSB containing 0.1% Tween-20 (without NP40 or digitonin) was added to wash the cells, and the tubes were inverted three times to mix. Cells were transported to a pre-chilled (4 °C) fixed-angle centrifuge to pellet the nuclei at 600× g for 10 min. The supernatant was removed with two pipetting steps (aspirate down to 100 µL with a p1000 pipette and remove final 100 µL with a p200 pipette) and nuclei were resuspended in 50 μL of transposition mix by pipetting up and down six times. The transposition reaction mixture was then incubated at 37 °C for 30 min in a thermomixer with 1000 r.p.m mixing. The transposed DNA was pre-amplified and purified by AMPure DNA magnetic beads. Fragmented DNA was eluted in 20.5 μL RSB containing 10 mM Tris, pH = 7.4. The quality of the libraries was assessed using Bioanalyzer High-Sensitivity DNA Analysis kit (Agilent, Santa Barbara, CA, USA).
Sequencing, data processing and peak calling: The libraries were pooled and sequenced to a mean depth of 80–100 million reads per sample on NovaSeq 50 bp paired end SP run. The quality of raw reads was checked using FastQC v0.11.7 [25] and MultiQC v1.6 [26]. Low-quality reads and adapter sequences were trimmed from fastQ files using Trimmomatic v0.36 [24]. Reads were mapped to the GRCh38 human reference genome using Bowtie2 v2.3.4.3 [27]. “MarkDuplicates” and “REMOVE_DUPLICATES” functions from Picard (v2.23.4) were then used to remove duplicate reads. Reads were filtered for non-unique alignment (MAPQ > 30) and low-quality mapping (-F 1804) using samtools (v1.1). Using alignmentSieve v3.5, reads were shifted to adjust for the inserted adapter sequences. Finally, peaks were called using MACS2 v2.2.7.1 [28].

2.13. Overlap of ESC CAGE Tag Clusters with ESC ATAC-Seq Peaks

CAGE-detected tag clusters were overlapped with ATAC-seq peaks generated in the same samples using BEDtools intersect [29]. Our approach is a stringent method to define regulatory genomic regions, where both CAGE and ATAC-seq signals are present.

2.14. Overlap of ESC Consensus Clusters with FANTOM5 and ENSEMBL Databases

Human transcription start site CAGE peaks data generated from 573 tissues and cell lines were obtained from FANTOM5 https://fantom.gsc.riken.jp/5/datafiles/reprocessed/hg38_latest/extra/CAGE_peaks_expression/hg38_liftover+new_CAGE_peaks_phase1and2.bed.gz (accessed on 30 June 2021). ENSEMBL TSS annotations for the human genome (ENSEMBL annotation version hg38) were obtained from ENSEMBL-Biomart. CAGE-detected consensus clusters were overlapped with human TSS from FANTOM5 and ENSEMBL TSS using the intersect BEDtools function -u for unique overlap.

2.15. Overlap of ESC Consensus Clusters with Endometriosis GWAS Signals

Endometriosis GWA meta-analysis summary statistics were available from Sapkota et al. 2017. CAGE-detected consensus clusters were overlapped with the SNPs reaching nominal significance (p < 1 × 10−5) in the GWAS using -u (unique) function in BEDtools. To account for the promoter region, genomic windows of 2 kb were extended on either side of the consensus clusters and overlapped with the GWAS signals using the intersect function in BEDtools. By mapping GWAS signals to promoter and enhancer regions, as well as genes regulated by these regulatory regions, we expected to find potential mechanisms of significant SNPs that were not within the protein-coding regions.

2.16. Pathway Analysis

To gain greater biological insight on differentially expressed consensus clusters and consensus clusters that were not annotated to FANTOM5 data, we used clusterProfiler [30], which has the ability to analyse and visualize data for enrichment analysis. Gene lists included those identified from the differential consensus cluster expression analysis and those annotated to the potential ESC specific consensus clusters not present in FANTOM5. Benjamini–Hochberg-corrected p-values of <0.05 were considered significant enrichment.

3. Results

3.1. Comparing Quality of the Commercial and In-House CAGE Sequencing Protocols

Raw tag counts were normalised to quantify the expression from each individual TSS and to enable comparison between multiple samples. Tag counts ranged from 10–20 million for the commercial protocol and 5–9 million for the in-house protocol (Table 2). The CAGEr package was used to count the number of CAGE-detected transcription start sites (CTSS normalised to the sequencing depth) in each sample processed using the in-house and commercial protocols. It was observed that the number of CTSS counts obtained from in-house protocol samples (mean = 95,000) were more than the commercial protocol samples (mean = 75,000) (Figure 2a). The relationship between CAGE tag counts in different samples of each protocol was assessed using Pearson’s correlation to calculate correlation coefficients (r2) between all pairs of samples. Samples processed using the commercial protocol were highly correlated, r2 ≥ 0.8, as were those processed using the in-house analysis, r2 ≥ 0.8.

3.2. Genomic Distribution of CAGE-Defined Transcription Start Sites (CTSS)

To investigate the genomic distribution of CAGE tags, CTSS were annotated to the human hg38 reference genome https://ftp.ensembl.org/pub/release-100/fasta/homo_sapiens/dna/Homo_sapiens.GRCh38.dna.primary_assembly.fa (accessed on 25 April 2020). Distribution of CTSS counts across all the samples and between the two protocols were found to have similar genomic distribution (Figure 2b). As anticipated, the majority (>80%) of CTSS were annotated to promoter regions, followed by an average of 10% annotated to exons, 5% annotated to introns and 2% intergenic.

3.3. Identification of ESC Promoter Elements

TSSs in close proximity give rise to functionally equivalent sets of transcripts and are the elements of the same promoter/enhancer. As such, to understand the transcriptional activity of a promoter, CTSS were clustered into tag clusters (TCs). An average of 8313 tag clusters were identified using the commercial protocol and 7770 tag clusters from the in-house protocol. To compare genome-wide transcriptional activity across the samples and to perform expression profiling, tag clusters were aggregated from all the samples into a single set of non-overlapping consensus clusters. Consensus clusters observed from in-house and commercial protocols were 6588 and 7136, respectively. These consensus clusters were annotated to genes. The relationship between consensus cluster tag counts obtained from different samples in each protocol was assessed using Pearson’s correlation to calculate correlation coefficients between all pairs of samples. Consensus cluster tag counts from samples processed using the commercial protocol were highly correlated (r2 ≥ 0.8) with counts from samples processed using in-house protocol.
Despite the use of two protocols in generating the CAGE data from ESCs, results of downstream analyses were consistent between both the protocols. Considering the cost and time, our inhouse protocol generated equivalent results to the commercial protocol but in a shorter time and for less cost.

3.4. Concordance between Sequencing Technologies

To determine if expression of TSSs, as estimated by the CAGE data, is consistent with overall gene expression in ESCs, tag counts were correlated for each gene for the single consensus clusters (annotated to genes) with RNA-seq data generated from the same eight ESC lines. It was observed that consensus cluster counts were significantly (p = 2.2 × 10−16) positively correlated (cor = 0.334) with RNA-seq data annotated to the same genes.

3.5. Transcriptional and Regulatory Elements Are Supported by ATAC-Seq Data

Tag clusters and enhancers identified from CAGE were overlapped with the chromatin accessibility data profiled using ATAC-seq. It was observed that almost 90% of the tag clusters and enhancers fall in open chromatin regions (Table 3).

3.6. Differential Expression of Promoters between Endometriosis Cases and Controls

Given the consistency between the two CAGE protocols, and to maximise read depth for further analysis, the data from the two methods were combined and re-analysed. An average of 8351 tag clusters (Table 4) and 7125 consensus clusters were generated from the merged data and the genomic distribution of CTSS counts was consistent with that observed for the commercial and in-house data separately. As expected, we observed two distinct populations of promoters, defined as sharp (interquartile range < 10 bp) and broad promoters (interquartile range > 10 bp) (Figure 3).
Differential expression analysis resulted in 210 consensus clusters which were significantly differentially expressed between endometriosis cases and controls following correction for multiple testing (FDR-adjusted p-value < 0.05) (Table 5 and Figure 4a). After Bonferroni correction for multiple testing (p < 7.01 × 10−6), 53 out of 210 consensus clusters remained significant (Table 5 and Figure 4a). Eight consensus clusters annotated to the gene insulin-like growth factor binding protein-5 (IGFBP5) and two annotated to Caldesmon 1 (CALD1) had significantly lower expression in cases compared to controls. Consensus clusters associated with genes for oxytocin receptor (OXTR) (Figure 4c) and pregnancy-specific beta-1-glycoprotein 4 (PSG4) also had lower expression in endometriosis cases in comparison with controls.
To gain further understanding into the biological mechanisms underlying the consensus clusters, pathway analysis was performed using the 210 significantly differentially expressed consensus clusters using gene ontology databases (Figure 4b). This revealed significantly enriched processes involved in focal adhesion (p = 3.8 × 10−7) and cell substrate adhesion (p = 4.22 × 10−7).
Annotations of consensus clusters to the genes IGFBP5, CALD1 and OXTR, and their overlap with open chromatin and ensemble regulatory regions, are shown in Figure 5. The majority of clusters annotated to IGFBP5 were located in the 3′ and 5′ UTR regions overlapping regions of open chromatin and predicted promoter and enhancer regions (Figure 5a). IGFBP5 clusters differentially expressed between cases and controls were located in promoters and enhancers. Clusters annotated to the CALD1 gene were observed in the 5′UTR, promoter and open chromatin regions near exonic regions. The differentially expressed clusters between cases and controls were found to overlap with enhancers and promoters. Only one cluster was discovered to be differently expressed between cases and controls in the 5′ UTR region of the OXTR gene, which overlapped with the promoter region.

3.7. Predicted Enhancers and Differential Expression Analysis between Endometriosis Cases and Controls

Using the CAGEfightR package, 1561 enhancers were identified. Of the 1561 predicted enhancers, 1081 and 480 enhancers were annotated to intron and intergenic regions, respectively. Following interaction analyses, 1176 predicted enhancer–TSS interactions were found where expression of a TSS and nearby enhancer was significantly correlated; however, none of the 1176 interactions remained significant following correction for multiple testing (FDR < 0.05).
Differential expression analysis resulted in two enhancers which were significantly differentially expressed between endometriosis cases and controls following correction for multiple testing (FDR adjusted p-value < 0.05) located in the intergenic region nearby to the genes LINC02547 and HTR1D.

3.8. ESC-Specific Promoters and Enhancers

To determine if consensus clusters identified in ESC were consistent with previously reported TSSs and enhancers, consensus clusters identified in this study were compared with human TSS data from FANTOM5. The FANTOM5 dataset consisted of 1,048,125 TSSs from 573 tissues and cell types. It was observed that 96.5% (6880) of the CAGE consensus clusters overlapped with the FANTOM5 TSSs. The remaining 244 consensus clusters that did not overlap with the FANTOM5 data were annotated to 194 genes, some of which were found to have multiple promoters; for example, there were eight novel promoters annotated to COL4A2, six to MALAT1, five to MYH9 and two to IGFBP5 genes. Of the 244 consensus clusters, 31 were found to be annotated to unknown promoters. Pathway analysis for the 194 genes revealed significantly enriched processes involved in focal adhesion (p = 1.34 × 10−27) and cell substrate adhesion (p = 1.34 × 10−27).
Enhancers identified in the ESC were compared with previously reported enhancers from the FANTOM5 to determine if results were consistent. The FANTOM5 dataset consisted of 63,285 enhancers. It was found that 60% of the ESC enhancers overlapped with FANTOM5 enhancers.

3.9. Comparison of ESC CAGE Data with ENSEMBL TSS Annotation

CAGE-detected consensus clusters were evaluated against ENSEMBL TSS annotations for the human genome (ENSEMBL annotation version hg38). Out of 7125 consensus clusters that were within the regions, only 2948 matched perfectly with the ENSEMBL TSS annotations. After extending the regions by 100 bp either side, an average of 4293 tags overlapped.

3.10. Evidence of Transcriptional Elements in Endometriosis Risk Regions

Consensus clusters and enhancers from endometrial stromal cells were overlapped with endometriosis risk SNPs. None of the consensus clusters or enhancers mapped to known GWAS signals. By extending the genomic window around the consensus clusters by 2 kb to account for the promoter region, consensus clusters were found to overlap 11 endometriosis-associated SNPs across 6 genomic regions (Table 6), including gene regions near FGD6, ARL14EP and CDC42 (Figure 6).

4. Discussion

Gene transcription is controlled by functional interactions between promoters and enhancers [31]. Identifying tissue-specific molecular signatures of active regulatory elements is critical to understanding gene regulatory mechanisms. A comprehensive understanding of factors affecting the regulation of transcription in the endometrium and endometrial cell types is important for the accurate analysis and interpretation of differential gene expression data from endometrium across biological and disease contexts because the endometrium is the likely source of cells initiating endometriotic lesions. In this study, for the first time, CAGE data were generated from endometrial stromal cells of endometriosis patients and controls to map genomic regulatory regions, including promoters and enhancers in endometrial stromal cells and the location of endometriosis risk loci within genomic regulatory regions.
Profiling TSSs using CAGE in eight endometrial stromal cells identified a consensus set of ~7000 reproducible TCs (or promoters) and 1500 enhancers, consistent with numbers of previously identified TCs from human pancreatic islets [32]. ESC TCs were enriched to occur in stromal cell promoter states and ATAC-seq peaks, which reflects the expected chromatin landscape at regions where transcription initiation occurs. Theoretically, it is understood that gene expression data obtained from different technologies should be consistent if they measure the same outcomes. CAGE-profiled gene expression measured in this study was comparable with the RNA-seq data from the same samples, as evidenced by a significant positive correlation. However, as expected, the genes with low expression were poorly quantified with CAGE profiling, and this probably reflects the modest sequencing depth in these samples [33]. Using CAGE, the FANTOM consortium has mapped TCs across multiple tissue and cell types from mice and humans [11]. Comparison of endometrial stromal cell TCs with those identified from FANTOM5 data showed a high degree of concordance and revealed that 6.3% of TCs were potentially stromal-cell-specific.
Despite the small sample size, there was evidence of endometriosis case–control differences at 210 consensus clusters. Of the 210 consensus clusters showing differential expression, 53 promoters passed Bonferroni correction, including genes (IGFBP5, OXTR CALD1) previously associated with endometriosis [15,34,35]. Of particular interest is the differential expression of several TCs for IGFBP5. IGFBP5 is one of a family of binding proteins involved in the regulation of insulin-like growth factor (IGF) signalling. Of the 15 consensus TCs observed for this gene, 13 TCs showed FDR-significant differential expression between endometriosis cases and controls. It was observed that 5 out of 284 TCs which were not annotated to FANTOM5 data belonged to IGFBP5. IGFBP5 showed consistently low expression in endometrial stromal cells obtained from endometriosis patients. IGFBP5 is selectively expressed in the proliferative phase of human endometrium, indicating that IGFBP5 may play a role in promoting endometrial cell proliferation [36]. Previous studies suggest that IGFBP5 is influenced by oestrogen and is consistently overexpressed in endometriotic tissue [37]. Differential expression of IGFBP5 in bulk endometrial tissue was not observed [3,4]. However, a recent study [38] found that IGFBP5 was significantly downregulated in a sub-population of stromal fibroblasts that was associated with endometriosis cases. Discrepancies in the level of expression of IGFBP5 between eutopic and ectopic endometrial tissue and in endometrial cell types is consistent with reports that IGFBP5 is expressed in a wide range of cell types but exhibits a broad range of biological functions depending on the context [38]. Studies have shown that IGFBP5 can both positively and negatively regulate IGF signalling in different cell types and tissues and similarly it has been shown to both inhibit and promote cell survival, proliferation and migration [38]. The potential for dysregulation of IGFBP5 to influence cell survival and proliferation highlight this gene as a good candidate for endometriosis.
Twelve consensus clusters were identified for CALD1 in endometrial stromal cells, six consensus clusters showed significantly low expression in women with endometriosis compared to women without endometriosis. The CALD1 gene encodes the protein caldesmon found in the cytoskeleton and plays an important role in smooth muscle contraction and relaxation, as well as cellular functions, such as proliferation, adhesion and cell motility [39,40]. CALD1 isoforms have also been associated with tumour malignancy in several cancers [41,42,43,44]. Previous studies have shown increased expression of the CALD1 gene in endometriotic tissue obtained from women with endometriosis; however, they also report a significantly lower expression of the caldesmon protein in eutopic endometrium of women with endometriosis, which is consistent with the lower levels of expression observed in this study [40,45,46]. It has been hypothesised that low levels of caldesmon protein in eutopic endometrium may increase the motility and invasiveness of the endometrium via the inability of the protein to bind to actin and inhibit the activity of actomyosin ATPase [40]. Alternatively, caldesmon activity is also negatively regulated by calmodulin (CALM2), which has increased expression in ectopic lesions alongside CALD1. Increased expression of CALM2 may also prevent caldesmon from binding to actin, allowing lesions to have a greater invasive potential [40]. The potential role of CALD1 in endometriosis pathogenesis in both the context of eutopic and ectopic endometrial cell types warrants further functional investigation.
This study also supports an association between the variation in expression of OXTR and endometriosis. One of the consensus clusters for OXTR was significantly different between women with and without endometriosis. Oxytocin and its receptor appear to play an important role in the regulation of contractions in the uterus [47]. OXTR is expressed in the smooth muscle cells and epithelial cells in the uterus and endometriotic lesions [48]. We show that OXTR is also expressed in stromal cells. Studies have shown variable expression of OXTR across the menstrual cycle [3,35,49] and in response to endometrial pathologies, including endometriosis, adenomyosis and recurrent implantation failure [35,50,51].
CAGE identified 1561 elements showing enhancer activity located in intron and intergenic regions. Transcribed enhancers were identified across a wide variety of human cells and tissues from the analysis of the FANTOM5 dataset [21]. In this study, 60% of the identified enhancers overlapped with the FANTOM5-identified enhancers. Potential tissue/cell-type-specific enhancers identified in this study may provide valuable insights into regulatory elements that play a distinct role in endometrial structure, function and disease. Ninety-three percent of the CAGE-defined enhancers were enriched in the open chromatin regions, providing further evidence for the presence of potential cell-type-specific enhancers.
Endometriosis is driven by both genetic and environmental factors. A significant breakthrough in genetic research on endometriosis came through GWAS, where 27 genomic risk regions have been identified [6]. However, many implicated variants are classified as non-coding and their significance can be realized only when the associated DNA sequence variants in the tissues relevant to endometriosis are determined. There was no evidence that risk variants associated with endometriosis risk were located within tag clusters or enhancers in endometrial stromal cells; however, there is evidence that functional sites regulating the activity of promoters are concentrated within a few thousand bases of the TSS [52]. Analysis of SNPs with regulatory effects on the expression of neighbouring genes shows that most variants with a significant effect are located within 2 kb of a TSS [52]. Considering the length of these regulatory regions, 11 GWAS SNPs were found to overlap within the 2 kb genomic window of the identified consensus clusters (Table 5). A few of the important gene regions identified here were near FGD6, ARL14EP and CDC42. In relation to endometriosis, several studies have reported association of endometriosis risk with increased expression of VEZT and FGD6 [3,4]. Both VEZT and FGD6 play an important role in plasma membrane, cell adhesion and cytoskeletal remodelling, which are important in lesion formation [53,54]. A recent study suggests the evidence of bidirectional promoters for VEZT and FGD6 [4]. The RNA-seq data revealed that both FGD6 and VEZT were expressed in endometrial stromal cells, although VEZT was expressed at a higher level than FGD6. From the CAGE data, consensus clusters were only detected for VEZT in this investigation. The absence of FGD6 expression may be attributed to insufficient sequencing depth.
Active TSSs in ESCs in the FSHB endometriosis risk locus were also identified. The region upstream of the promoter of FSHB (follicle-stimulating hormone subunit B), regulating FSH (follicle-stimulating hormone) concentrations, has been associated with risk of endometriosis and several other gynaecological traits and diseases [55]. ARL14EP is located on the same chromosomal locus as FSHB. ARL14EP is expressed in many tissue types, with relatively high levels in the ovary, testis and uterus, and plays a role in controlling the export of major histocompatibility class II molecules along the actin cytoskeleton [56,57]. FSHB is only expressed in the pituitary gland and is the beta subunit for FSH, a gonadotropin with key regulatory roles in reproductive function [58].
Endometriosis risk alleles in the chr1 region have been associated with increased expression of CDC42 [59], and active TSSs in ESCs were also located in the CDC42 endometriosis risk locus. CDC42 is a member of the Rho family of small GTPases and is thought to regulate a variety of cellular processes, including cell cycle progression, cell polarity, cytoskeletal reorganization and transcription [60]. GnRH-activated CDC42 regulates FSH and LH in response to pulsatile GnRH [61], there is also evidence that CDC42 regulates ovarian reserve, follicle activation and granulosa cell function [62], suggesting that altered CDC42 regulation may be implicated in fertility issues associated with endometriosis.
This study has a number of strengths, the most notable being the generation of the first CAGE dataset generated from ESC, providing cell-type specific data on gene regulatory mechanisms. The ability to compare this dataset with other OMIC datasets generated from the same samples (ATAC-seq, RNA-seq, genotypes) and other publicly available datasets (e.g., FANTOM5) provides additional evidence for the validation of gene regulatory mechanisms and facilitates the identification of 244 TSSs that are not observed in the FANTOM5 datasets. These TSSs may be specific to endometrial stromal cells and has led to the identification of some active TSSs with distinct transcriptional regulation patterns linked to endometriosis status. Differential regulation of these genes was not previously detected in bulk endometrial tissue. Further studies of differential regulation of IGFBP5, CALD1 and OXTR in specific endometrial cell types with disease status will be required to validate the reported results. The power of this dataset is, however, limited by the small sample size and limited sequencing depth. Individual results have not been validated and will need confirmation from alternative methods such as RT-PCR, Western blotting or immunohistochemistry.
Characterising active TSSs in endometrial stromal cells revealed differences in transcriptional regulation associated with endometriosis that have not been observed in bulk endometrial tissue. This study highlights the importance of mapping cell-type-specific transcriptional events that may contribute to disease processes and provides a unique resource to investigate transcriptional regulation in a cell type important in many reproductive traits and conditions.

Author Contributions

Conceptualization, S.M. (Sushma Marla), S.M. (Sally Mortlock) and G.W.M.; methodology, M.D.M., B.M., S.Y., S.A. and J.C.; software, S.M. (Sushma Marla), S.M. (Sally Mortlock) and Q.N.; validation, S.M. (Sushma Marla), S.M. (Sally Mortlock) and G.W.M.; formal analysis, S.M. (Sushma Marla) and S.M. (Sally Mortlock); investigation, S.M. (Sushma Marla), S.M. (Sally Mortlock) and G.W.M.; resources, G.W.M.; data curation, S.M. (Sushma Marla), S.M. (Sally Mortlock) and Q.N.; writing—original draft preparation, S.M. (Sushma Marla); writing—review and editing, S.M. (Sally Mortlock) and G.W.M.; supervision, S.M. (Sally Mortlock) and G.W.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by NHMRC fellowships GNT1078399, GNT1177194 and GNT2008928.

Institutional Review Board Statement

Tissue sample collection was approved by the Cantonal Ethics Commission Bern (149/03). Experimental procedures were approved by the Cantonal Ethics Commission Bern (2019-01146) and the University of Queensland Human Research Ethics Committee (2016001723) (2019/HE002744).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data underlying this article will be shared on reasonable request to the corresponding author.

Acknowledgments

We thank the Genome Innovation Hub at the University of Queensland, for support for this project and Jun Xu for his advice and assistance with data analysis.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Evans, J.; Salamonsen, L.A.; Winship, A.; Menkhorst, E.; Nie, G.; Gargett, C.E.; Dimitriadis, E. Fertile ground: Human endometrial programming and lessons in health and disease. Nat. Rev. Endocrinol. 2016, 12, 654–667. [Google Scholar] [CrossRef] [PubMed]
  2. Spencer, T.E.; Hayashi, K.; Hu, J.; Carpenter, K.D. Comparative Developmental Biology of the Mammalian Uterus. In Current Topics in Developmental Biology; Academic Press: Cambridge, MA, USA, 2005; Volume 68, pp. 85–122. [Google Scholar] [CrossRef]
  3. Fung, J.N.; Mortlock, S.; Girling, J.E.; Holdsworth-Carson, S.J.; Teh, W.T.; Zhu, Z.; Lukowski, S.W.; McKinnon, B.D.; McRae, A.; Yang, J.; et al. Genetic regulation of disease risk and endometrial gene expression highlights potential target genes for endometriosis and polycystic ovarian syndrome. Sci. Rep. 2018, 8, 11424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Mortlock, S.; Kendarsari, I.R.; Fung, J.N.; Gibson, G.; Yang, F.; Restuadi, R.; Girling, E.J.; Holdsworth-Carson, S.J.; Teh, W.T.; Lukowski, S.W.; et al. Tissue specific regulation of transcription in endometrium and association with disease. Hum. Reprod. 2020, 35, 377–393. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Giudice, L.C.; Kao, L.C. Endometriosis. Lancet 2004, 364, 1789–1799. [Google Scholar] [CrossRef]
  6. Nilufer, R.; Banasik, K.; Paraskevi, C.; Danning, R.; Genevieve, G.; Ayush, G.; MacGregor, S.; Mortlock, S.; Sapkota, Y.; Schork, J.A.; et al. Large-scale genome-wide association meta-analysis of endometriosis reveals 13 novel loci and genetically-associated comorbidity with other pain conditions. bioRxiv 2018, 406967. [Google Scholar]
  7. Tanos, P.; Dimitriou, S.; Gullo, G.; Tanos, V. Biomolecular and Genetic Prognostic Factors That Can Facilitate Fertility-Sparing Treatment (FST) Decision Making in Early Stage Endometrial Cancer (ES-EC): A Systematic Review. Int. J. Mol. Sci. 2022, 23, 2653. [Google Scholar] [CrossRef]
  8. Dyson, M.; Roqueiro, D.; Monsivais, D.; Ercan, C.M.; Pavone, M.E.; Brooks, D.C.; Kakinuma, T.; Ono, M.; Jafari, N.; Dai, Y.; et al. Genome-Wide DNA Methylation Analysis Predicts an Epigenetic Switch for GATA Factor Expression in Endometriosis. PLoS Genet. 2014, 10, e1004158. [Google Scholar] [CrossRef] [Green Version]
  9. Shiraki, T.; Kondo, S.; Katayama, S.; Waki, K.; Kasukawa, T.; Kawaji, H.; Kodzius, R.; Watahiki, A.; Nakamura, M.; Arakawa, T.; et al. Cap analysis gene expression for high-throughput analysis of transcriptional starting point and identification of promoter usage. Proc. Natl. Acad. Sci. USA 2003, 100, 15776–15781. [Google Scholar] [CrossRef] [Green Version]
  10. Kodzius, R.; Kojima, M.; Nishiyori, H.; Nakamura, M.; Fukuda, S.; Tagami, M.; Sasaki, D.; Imamura, K.; Kai, C.; Harbers, M.; et al. CAGE: Cap analysis of gene expression. Nat. Methods 2006, 3, 211–222. [Google Scholar] [CrossRef]
  11. Forrest, A.R.R.; Kawaji, H.; Rehli, M.; Baillie, K.J.; de Hoon, M.J.L.; Haberle, V.; Lassmann, T.; Kulakovskiy, I.V.; Lizio, M.; Itoh, M.; et al. A promoter-level mammalian expression atlas. Nature 2014, 507, 462–470. [Google Scholar]
  12. American Society for Reproductive Medicine. Revised American Society for Reproductive Medicine classification of endometriosis: 1996. Fertil. Steril. 1997, 67, 817–821. [Google Scholar] [CrossRef] [PubMed]
  13. Burla, L.; Scheiner, D.; Samartzis, E.P.; Seidel, S.; Eberhard, M.; Fink, D.; Boss, A.; Imesch, P. The ENZIAN score as a preoperative MRI-based classification instrument for deep infiltrating endometriosis. Arch. Gynecol. Obstet. 2019, 300, 109–116. [Google Scholar] [CrossRef] [PubMed]
  14. McKinnon, B.D.; Evers, J.; Bersinger, N.A.; Mueller, M. Induction of the Neurokinin 1 Receptor by TNFα in Endometriotic Tissue Provides the Potential for Neurogenic Control Over Endometriotic Lesion Growth. J. Clin. Endocrinol. Metab. 2013, 98, 2469–2477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. McKinnon, B.D.; Lukowski, S.W.; Mortlock, S.; Crawford, J.; Atluri, S.; Subramaniam, S.; Johnston, R.L.; Nirgianakis, K.; Tanaka, K.; Amoako, A.; et al. Altered differentiation of endometrial mesenchymal stromal fibroblasts is associated with endometriosis susceptibility. Commun. Biol. 2022, 5, 600. [Google Scholar] [CrossRef]
  16. Haberle, V.; Forrest, A.; Hayashizaki, Y.; Carninci, P.; Lenhard, B. CAGEr: Precise TSS data retrieval and high-resolution promoterome mining for integrative analyses. Nucleic Acids Res. 2015, 43, e51. [Google Scholar] [CrossRef] [Green Version]
  17. Ong, C.; Corces, V.G. Enhancers: Emerging roles in cell fate specification. EMBO Rep. 2012, 13, 423–430. [Google Scholar] [CrossRef] [Green Version]
  18. Murakawa, Y.; Yoshihara, M.; Kawaji, H.; Nishikawa, M.; Zayed, H.; Suzuki, H.; FANTOM Consortium; Hayashizaki, Y. Enhanced Identification of Transcriptional Enhancers Provides Mechanistic Insights into Diseases. Trends Genet. 2016, 32, 76–88. [Google Scholar] [CrossRef]
  19. Kim, T.-K.; Hemberg, M.; Gray, J.M.; Costa, A.M.; Bear, D.M.; Wu, J.; Harmin, D.A.; Laptewicz, M.; Barbara-Haley, K.; Kuersten, S.; et al. Widespread transcription at neuronal activity-regulated enhancers. Nature 2010, 465, 182–187. [Google Scholar] [CrossRef] [Green Version]
  20. Heinz, S.; Romanoski, C.E.; Benner, C.; Glass, C.K. The selection and function of cell type-specific enhancers. Nat. Rev. Mol. Cell Biol. 2015, 16, 144–154. [Google Scholar] [CrossRef] [Green Version]
  21. Andersson, R.; Gebhard, C.; Miguel-Escalada, I.; Hoof, I.; Bornholdt, J.; Boyd, M.; Chen, Y.; Zhao, X.; Schmidl, C.; Suzuki, T.; et al. An atlas of active enhancers across human cell types and tissues. Nature 2014, 507, 455–461. [Google Scholar] [CrossRef] [Green Version]
  22. Thodberg, M.; Thieffry, A.; Vitting-Seerup, K.; Andersson, R.; Sandelin, A. CAGEfightR: Analysis of 5′-end data using R/Bioconductor. BMC Bioinform. 2019, 20, 487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  25. Andrews, S. FastQC: A quality control tool for high throughput sequence data. In Babraham Bioinformatics; Babraham Institute: Cambridge, UK, 2010. [Google Scholar]
  26. Ewels, P.; Magnusson, M.; Lundin, S.; Käller, M. MultiQC: Summarize analysis results for multiple tools and samples in a single report. Bioinformatics 2016, 32, 3047–3048. [Google Scholar] [CrossRef] [Green Version]
  27. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [Green Version]
  28. Zhang, Y.; Liu, T.; Meyer, C.A.; Eeckhoute, J.; Johnson, D.S.; Bernstein, B.E.; Nusbaum, C.; Myers, R.M.; Brown, M.; Li, W.; et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9, R137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Quinlan, A.R.; Hall, I.M. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 2010, 26, 841–842. [Google Scholar] [CrossRef] [Green Version]
  30. Yu, G.; Wang, L.-G.; Han, Y.; He, Q.-Y. clusterProfiler: An R Package for Comparing Biological Themes Among Gene Clusters. OMICS J. Integr. Biol. 2012, 16, 284–287. [Google Scholar] [CrossRef]
  31. Mikhaylichenko, O.; Bondarenko, V.; Harnett, D.; Schor, I.E.; Males, M.; Viales, R.R.; Furlong, E.E.M. The degree of enhancer or promoter activity is reflected by the levels and directionality of eRNA transcription. Genes Dev. 2018, 32, 42–57. [Google Scholar] [CrossRef] [Green Version]
  32. Varshney, A.; Kyono, Y.; Elangovan, V.R.; Wang, C.; Erdos, M.R.; Narisu, N.; Albanus, R.D.; Orchard, P.; Stitzel, M.L.; Collins, F.S.; et al. A Transcription Start Site Map in Human Pancreatic Islets Reveals Functional Regulatory Signatures. Diabetes 2021, 70, 1581–1591. [Google Scholar] [CrossRef]
  33. Kawaji, H.; Lizio, M.; Itoh, M.; Kanamori-Katayama, M.; Kaiho, A.; Nishiyori-Sueki, H.; Shin, J.W.; Kojima-Ishiyama, M.; Kawano, M.; Murata, M.; et al. Comparison of CAGE and RNA-seq transcriptome profiling using clonally amplified and single-molecule next-generation sequencing. Genome Res. 2014, 24, 708–717. [Google Scholar] [CrossRef] [Green Version]
  34. Meola, J.; Hidalgo, G.D.S.; e Silva, J.C.R.; Silva, L.E.C.M.; Paz, C.C.P.; Ferriani, R.A. Caldesmon: New Insights for Diagnosing Endometriosis1. Biol. Reprod. 2013, 88, 122. [Google Scholar] [CrossRef]
  35. Huang, M.; Li, X.; Guo, P.; Yu, Z.; Xu, Y.; Wei, Z. The abnormal expression of oxytocin receptors in the uterine junctional zone in women with endometriosis. Reprod. Biol. Endocrinol. 2017, 15, 1–10. [Google Scholar] [CrossRef] [Green Version]
  36. Zhou, J.; Dsupin, B.A.; Guidice, L.C.; Bondy, C.A. Insulin-like growth factor system gene expression in human endometrium during the menstrual cycle. J. Clin. Endocrinol. Metab. 1994, 79, 1723–1734. [Google Scholar] [PubMed] [Green Version]
  37. Hu, W.-P.; Tay, S.K.; Zhao, Y. Endometriosis-Specific Genes Identified by Real-Time Reverse Transcription-Polymerase Chain Reaction Expression Profiling of Endometriosis Versus Autologous Uterine Endometrium. J. Clin. Endocrinol. Metab. 2006, 91, 228–238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Duan, C.; Allard, J.B. Insulin-Like Growth Factor Binding Protein-5 in Physiology and Disease. Front. Endocrinol. 2020, 11, 100. [Google Scholar] [CrossRef] [Green Version]
  39. Gusev, N.B. Some properties of caldesmon and calponin and the participation of these proteins in regulation of smooth muscle contraction and cytoskeleton formation. Biochemistry 2001, 66, 1112–1121. [Google Scholar] [CrossRef]
  40. Meola, J.; Rosa e Silva, J.C.; Dentillo, D.B.; da Silva, W.A., Jr.; Veiga-Castelli, L.C.; Bernardes, L.A.; Ferriani, R.A.; de Paz, C.C.; Giuliatti, S.; Martelli, L. Differentially expressed genes in eutopic and ectopic endometrium of women with endometriosis. Fertil. Steril. 2010, 93, 1750–1773. [Google Scholar] [CrossRef]
  41. Kim, K.-H.; Yeo, S.-G.; Kim, W.K.; Kim, D.Y.; Yeo, H.Y.; Hong, J.P.; Chang, H.J.; Park, J.W.; Kim, S.Y.; Kim, B.C.; et al. Up-regulated expression of l-caldesmon associated with malignancy of colorectal cancer. BMC Cancer 2012, 12, 601. [Google Scholar] [CrossRef] [Green Version]
  42. Chang, K.-P.; Wang, C.-L.A.; Kao, H.-K.; Liang, Y.; Bs, S.-C.L.; Bs, L.-L.H.; Hseuh, C.; Hsieh, Y.-J.; Chien, K.-Y.; Chang, Y.-S.; et al. Overexpression of caldesmon is associated with lymph node metastasis and poorer prognosis in patients with oral cavity squamous cell carcinoma. Cancer 2013, 119, 4003–4011. [Google Scholar] [CrossRef]
  43. Lee, M.-S.; Lee, J.; Kim, J.H.; Kim, W.T.; Kim, W.-J.; Ahn, H.; Park, J. Overexpression of caldesmon is associated with tumor progression in patients with primary non-muscle-invasive bladder cancer. Oncotarget 2015, 6, 40370–40384. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Lian, H.; Wang, A.; Shen, Y.; Wang, Q.; Zhou, Z.; Zhang, R.; Li, K.; Liu, C.; Jia, H. Identification of novel alternative splicing isoform biomarkers and their association with overall survival in colorectal cancer. BMC Gastroenterol. 2020, 20, 171. [Google Scholar] [CrossRef] [PubMed]
  45. Vouk, K.; Šmuc, T.; Guggenberger, C.; Ribič-Pucelj, M.; Šinkovec, J.; Husen, B.; Thole, H.; Houba, P.; Thaete, C.; Adamski, J.; et al. Novel estrogen-related genes and potential biomarkers of ovarian endometriosis identified by differential expression analysis. J. Steroid Biochem. Mol. Biol. 2011, 125, 231–242. [Google Scholar] [CrossRef]
  46. Mettler, L.; Salmassi, A.; Schollmeyer, T.; Schmutzler, A.G.; Püngel, F.; Jonat, W. Comparison of c-DNA microarray analysis of gene expression between eutopic endometrium and ectopic endometrium (endometriosis). J. Assist. Reprod. Genet. 2007, 24, 249–258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Shojo, H.; Kaneko, Y. Characterization and Expression of Oxytocin and the Oxytocin Receptor. Mol. Genet. Metab. 2000, 71, 552–558. [Google Scholar] [CrossRef]
  48. Mechsner, S.; Bartley, J.; Loddenkemper, C.; Salomon, D.S.; Starzinski-Powitz, A.; Ebert, A.D. Oxytocin receptor expression in smooth muscle cells of peritoneal endometriotic lesions and ovarian endometriotic cysts. Fertil. Steril. 2005, 83 (Suppl. 4), 1220–1231. [Google Scholar] [CrossRef] [PubMed]
  49. Fuchs, A.-R.; Fuchs, F.; Soloff, M.S. Oxytocin Receptors in Nonpregnant Human Uterus. J. Clin. Endocrinol. Metab. 1985, 60, 37–41. [Google Scholar] [CrossRef]
  50. Guo, S.-W.; Mao, X.; Ma, Q.; Liu, X. Dysmenorrhea and its severity are associated with increased uterine contractility and overexpression of oxytocin receptor (OTR) in women with symptomatic adenomyosis. Fertil. Steril. 2013, 99, 231–240. [Google Scholar] [CrossRef] [PubMed]
  51. Liu, C.; Li, Y.; Li, L.; Shui, S.; Yang, L.; Sui, C.; Zhang, H. Aberrant expression of oxytocin receptor in endometrium and decidua in women who have experienced recurrent implantation failure. F&S Sci. 2020, 1, 183–187. [Google Scholar] [CrossRef]
  52. Stranger, E.B.; Nica, A.C.; Forrest, M.S.; Dimas, A.; Bird, C.P.; Beazley, C.; Ingle, E.C.; Dunning, M.; Flicek, P.; Koller, D.; et al. Population genomics of human gene expression. Nat. Genet. 2007, 39, 1217–1224. [Google Scholar] [CrossRef] [Green Version]
  53. Guo, X.; Jing, C.; Li, L.; Zhang, L.; Shi, Y.; Wang, J.; Liu, J.; Li, C. Down-regulation of VEZT gene expression in human gastric cancer involves promoter methylation and miR-43c. Biochem. Biophys. Res. Commun. 2011, 404, 622–627. [Google Scholar] [CrossRef] [PubMed]
  54. Holdsworth-Carson, S.J.; Fung, J.N.; Luong, H.T.; Sapkota, Y.; Bowdler, L.M.; Wallace, L.; Teh, W.T.; Powell, J.E.; Girling, J.E.; Healey, M.; et al. Endometrial vezatin and its association with endometriosis risk. Hum. Reprod. 2016, 31, 999–1013. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. McGrath, I.M.; Mortlock, S.; Montgomery, G.W. Genetic Regulation of Physiological Reproductive Lifespan and Female Fertility. Int. J. Mol. Sci. 2021, 22, 2556. [Google Scholar] [CrossRef] [PubMed]
  56. GTEx Consortium. Erratum: Genetic effects on gene expression across human tissues. Nature 2018, 553, 530. [Google Scholar] [CrossRef] [Green Version]
  57. Paul, P.; Hoorn, T.V.D.; Jongsma, M.L.; Bakker, M.J.; Hengeveld, R.; Janssen, L.; Cresswell, P.; Egan, D.A.; van Ham, M.; Brinke, A.T.; et al. A Genome-wide Multidimensional RNAi Screen Reveals Pathways Controlling MHC Class II Antigen Presentation. Cell 2011, 145, 268–283. [Google Scholar] [CrossRef] [Green Version]
  58. Trevisan, C.M.; de Oliveira, R.; Christofolini, D.M.; Barbosa, C.P.; Bianco, B. Effects of a Polymorphism in the Promoter Region of the Follicle-Stimulating Hormone Subunit Beta (FSHB) Gene on Female Reproductive Outcomes. Genet. Test. Mol. Biomarkers 2018, 23, 39–44. [Google Scholar] [CrossRef]
  59. Powell, J.E.; Fung, J.N.; Shakhbazov, K.; Sapkota, Y.; Cloonan, N.; Hemani, G.; Hillman, K.M.; Kaufmann, S.; Luong, H.T.; Bowdler, L.; et al. Endometriosis risk alleles at 1p36.12 act through inverse regulation of CDC42 and LINC00339. Hum. Mol. Genet. 2016, 25, 5046–5058. [Google Scholar]
  60. Stengel, K.; Zheng, Y. Cdc42 in oncogenic transformation, invasion, and tumorigenesis. Cell. Signal. 2011, 23, 1415–1423. [Google Scholar] [CrossRef] [Green Version]
  61. Stamatiades, G.A.; Kaiser, U.B. Gonadotropin regulation by pulsatile GnRH: Signaling and gene expression. Mol. Cell. Endocrinol. 2018, 463, 131–141. [Google Scholar] [CrossRef]
  62. Mei, Q.; Li, H.; Liu, Y.; Wang, X.; Xiang, W. Advances in the study of CDC42 in the female reproductive system. J. Cell. Mol. Med. 2022, 26, 16–24. [Google Scholar] [CrossRef]
Figure 1. An overview of the CAGE bioinformatics analysis workflow.
Figure 1. An overview of the CAGE bioinformatics analysis workflow.
Cells 12 01736 g001
Figure 2. CAGE-identified transcription start sites (CTSS) and genomic distribution of CAGE tags of endometrial stromal cells obtained from commercial and in-house protocols. (a) CTSS count comparisons between commercial (orange) and in-house (blue) protocols of CAGE data. X-axis represents sample numbers of endometrial stromal cell lines and Y-axis represents number of CTSS counts; (b) genomic distribution of CAGE tags of endometrial stromal cells obtained from commercial and in-house protocols.
Figure 2. CAGE-identified transcription start sites (CTSS) and genomic distribution of CAGE tags of endometrial stromal cells obtained from commercial and in-house protocols. (a) CTSS count comparisons between commercial (orange) and in-house (blue) protocols of CAGE data. X-axis represents sample numbers of endometrial stromal cell lines and Y-axis represents number of CTSS counts; (b) genomic distribution of CAGE tags of endometrial stromal cells obtained from commercial and in-house protocols.
Cells 12 01736 g002
Figure 3. Endometrial stromal cells show two types of promoters: those that initiate transcription over a narrow region (<10 bp, sharp promoter) or over a broad range (10–100 bp).
Figure 3. Endometrial stromal cells show two types of promoters: those that initiate transcription over a narrow region (<10 bp, sharp promoter) or over a broad range (10–100 bp).
Cells 12 01736 g003
Figure 4. (a) Volcano plot for differential expression analysis of consensus clusters/promoters (from merged) between endometriosis controls and cases; (b) pathway analysis for 210 FDR significant consensus clusters; (c) top 14 differentially expressed consensus clusters in cases and controls obtained from the merged data.
Figure 4. (a) Volcano plot for differential expression analysis of consensus clusters/promoters (from merged) between endometriosis controls and cases; (b) pathway analysis for 210 FDR significant consensus clusters; (c) top 14 differentially expressed consensus clusters in cases and controls obtained from the merged data.
Cells 12 01736 g004
Figure 5. Maps of the chromosomal positions and gene context with consensus CAGE clusters annotated to the genes IGFBP5 (a), CALD1 (b) and OXTR (c), highlighting the overlap of cage clusters, open chromatin and Ensembl regulatory features. CAGE clusters are denoted by peaks on the CAGE clusters track. CAGE clusters in green were significantly differentially expressed following multiple testing correction (FDR p < 0.05) between endometriosis cases and controls, while clusters in purple were not annotated in FANTOM5. ATAC-seq signals are displayed as red peaks, and the corresponding narrow peak call is indicated by the red bars underneath the peaks. In the Ensembl regulatory features track, the following features are depicted as rectangular blocks with promoters (green), enhancers (brown), CTCF binding sites (blue) and open chromatin (pink). Where multiple features are clustered and would only be distinguished by enlarging the scale, these regions are denoted by dark green rectangular blocks.
Figure 5. Maps of the chromosomal positions and gene context with consensus CAGE clusters annotated to the genes IGFBP5 (a), CALD1 (b) and OXTR (c), highlighting the overlap of cage clusters, open chromatin and Ensembl regulatory features. CAGE clusters are denoted by peaks on the CAGE clusters track. CAGE clusters in green were significantly differentially expressed following multiple testing correction (FDR p < 0.05) between endometriosis cases and controls, while clusters in purple were not annotated in FANTOM5. ATAC-seq signals are displayed as red peaks, and the corresponding narrow peak call is indicated by the red bars underneath the peaks. In the Ensembl regulatory features track, the following features are depicted as rectangular blocks with promoters (green), enhancers (brown), CTCF binding sites (blue) and open chromatin (pink). Where multiple features are clustered and would only be distinguished by enlarging the scale, these regions are denoted by dark green rectangular blocks.
Cells 12 01736 g005
Figure 6. Maps of the chromosomal positions and gene context with consensus CAGE clusters and overlap with endometriosis risk SNPs in genomic regions for the CDC42 region (a), FGD6/VEZT region (b) and ARL14EP region (c). Location of SNPs significantly associated with endometriosis from GWAS studies are shown on the GWAS SNPs track with the corresponding SNP IDs. The promoter region of each gene is denoted by the grey-shaded area on the figures and demonstrates the overlap between CAGE clusters (peaks on the CAGE cluster track), open chromatin (red peaks on the ATAC-seq track with corresponding narrow peak call shown by the dark red bar underneath) and Ensembl promoter features shown in green on the Ensembl regulatory features track.
Figure 6. Maps of the chromosomal positions and gene context with consensus CAGE clusters and overlap with endometriosis risk SNPs in genomic regions for the CDC42 region (a), FGD6/VEZT region (b) and ARL14EP region (c). Location of SNPs significantly associated with endometriosis from GWAS studies are shown on the GWAS SNPs track with the corresponding SNP IDs. The promoter region of each gene is denoted by the grey-shaded area on the figures and demonstrates the overlap between CAGE clusters (peaks on the CAGE cluster track), open chromatin (red peaks on the ATAC-seq track with corresponding narrow peak call shown by the dark red bar underneath) and Ensembl promoter features shown in green on the Ensembl regulatory features track.
Cells 12 01736 g006
Table 1. Information on clinical parameters for participants included in the CAGE analysis.
Table 1. Information on clinical parameters for participants included in the CAGE analysis.
Sample IDCase/ControlrAFS %ENZIAN &AgeBMISmoking #Infertility *
Sample 1CaseI0432500
Sample 2Control005124.300
Sample 3Control004124.201a
Sample 4CaseIA0B2C04620.700
Sample 5CaseI03622.911a
Sample 6CaseII2B3121.300
Sample 7Control0025.728.7n.d.0
Sample 8Control0021.717.200
% rAFS score for endometriosis disease stage for cases were given according to the American Fertility Society scoring system (12). & ENZIAN scores provide a detailed description of the location and size or severity of the lesions and were given according to the Enzian classification (13). # Smoking status; 1 = Yes and n.d.= no data. * Infertility status; 1a = primary and 1b = secondary.
Table 2. Number of CAGE tag counts obtained in each sample (women with endometriosis = cases; women without endometriosis = controls) from commercial and in-house protocols.
Table 2. Number of CAGE tag counts obtained in each sample (women with endometriosis = cases; women without endometriosis = controls) from commercial and in-house protocols.
Sample IDCase/ControlCommercialIn-House
Sample 1Case13,878,6687,424,199
Sample 2Control11,976,9556,373,486
Sample 3Control10,337,5866,380,118
Sample 4Case13,202,9816,579,239
Sample 5Case14,025,0095,394,546
Sample 6Case14,100,3807,617,671
Sample 7Control16,107,9179,768,981
Sample 8Control20,367,5987,825,180
Table 3. Overlap of CAGE tag clusters with ATAC-seq data of endometrial stromal cells.
Table 3. Overlap of CAGE tag clusters with ATAC-seq data of endometrial stromal cells.
Sample_IDCase/ControlNumber of Tag ClustersATAC-Seq PeaksNo. of Tag Clusters Overlapping with ATAC-Seq
Sample 1Case83431145836688
Sample 2Control83371380936933
Sample 3Control83271256156869
Sample 4Case83461169576860
Sample 5Case83191372157005
Sample 6Case83461215186913
Sample 7Control83891221496834
Sample 8Control84031222146802
Table 4. Number of CAGE tag clusters identified in each sample after merging commercial and in-house analysis data.
Table 4. Number of CAGE tag clusters identified in each sample after merging commercial and in-house analysis data.
Sample IDCase/ControlTag Clusters
Sample 1Case8343
Sample 2Control8337
Sample 3Control8327
Sample 4Case8346
Sample 5Case8319
Sample 6Case8346
Sample 7Control8389
Sample 8Control8403
Table 5. List of 53 consensus promoters differentially expressed in women without endometriosis and with endometriosis following Bonferroni correction for multiple testing (p < 7.02 × 10−6).
Table 5. List of 53 consensus promoters differentially expressed in women without endometriosis and with endometriosis following Bonferroni correction for multiple testing (p < 7.02 × 10−6).
ChrStartEndStrandTPM #AnnotationGenesLog2FCp-ValueFDR *
chr221669483821669491513.00promoterIGFBP53.053.00 × 10−162.13 × 10−12
chr38769610876962826.97promoterOXTR2.652.24 × 10−157.93 × 10−12
chr22166955472166955591476.3promoterIGFBP52.595.23 × 10−141.24 × 10−10
chr514790653614790654110.87promoterC5orf463.621.29 × 10−122.29 × 10−9
chr14101964571101964575+44.23promoterDYNC1H11.006.75 × 10−129.57 × 10−9
chr59877366198773701+41.68promoterRGMB1.351.23 × 10−101.45 × 10−7
chr194320563343205648108.38promoterPSG43.741.89 × 10−101.92 × 10−7
chr418672666518672672828.46promoterFAT11.582.54 × 10−102.25 × 10−7
chr221667645021667670980.13exonIGFBP52.473.49 × 10−102.75 × 10−7
chr411436491111436495127.86unknown 2.378.15 × 10−105.78 × 10−7
chr221667432121667445220.13exonIGFBP52.761.67 × 10−91.08 × 10−6
chr1218345334218345346+76.83promoterTGFB22.443.02 × 10−91.78 × 10−6
chr7134711476134711515+34.95unknown 1.453.99 × 10−92.17 × 10−6
chr97829712278297161+36.46promoterPSAT1−1.946.68 × 10−93.38 × 10−6
chr3188212669188212713+37.34promoterLPP1.8581.19 × 10−85.62 × 10−6
chr7134646835134646861+29.33promoterBPGM0.951.91 × 10−88.47 × 10−6
chr22166950592166951408.79promoterIGFBP52.653.88 × 10−81.62 × 10−5
chr49445190194451973+109.58promoterPDLIM50.835.87 × 10−82.31 × 10−5
chr62602157226021654+23.35promoter −1.436.45 × 10−82.40 × 10−5
chr54067991440679919+5.58promoterPTGER4−1.227.21 × 10−82.55 × 10−5
chr1630656313065641+6.75promoterIL322.339.46 × 10−83.19 × 10−5
chr9729530397295307330.88promoterALDH1A1−2.262.21 × 10−77.13 × 10−5
chr61319495561319495694.73exonCCN21.852.71 × 10−78.33 × 10−5
chr178026083280260882+19.59promoterAC124319.11.022.94 × 10−78.69 × 10−5
chr12563158755631603986.96promoterAC073896.1;
CNPY2
-0.443.73 × 10−70.00010
chr18459672674596733021.44promoterEPG50.8384.59 × 10−70.00012
chr22166756812166756832.52exonIGFBP52.854.70 × 10−70.00012
chr156304268063042756+659.91promoterTPM10.864.96 × 10−70.00012
chr22166761282166761313.19exonIGFBP53.595.28 × 10−70.00012
chr7134867737134867770+7.57exonCALD11.496.20 × 10−70.00014
chr5178204530178204537+114.22promoterHNRNPAB−0.757.73 × 10−70.00017
chr93839267138392799+44.95promoterALDH1B11.9549.92 × 10−70.00021
chr955104985510556+16.50promoterPDCD1LG21.1491.25 × 10−60.00026
chr3156674590156674634+27.80promoterTIPARP0.931.29 × 10−60.00026
chr82340411823404156732.82promoterLOXL2;
ENTPD4
0.8231.43 × 10−60.00028
chr12297839102978392217.55promoterTMTC1−1.761.55 × 10−60.00029
chr841309471413094743.24promoterSFRP1−2.641.53 × 10−60.00029
chr5139293674139293754+75.81promoterMATR3-0.531.61 × 10−60.00029
chr5843843808438448319.82promoterEDIL31.621.89 × 10−60.00034
chr205019082950190835+51.41promoterCEBPB−0.8552.31 × 10−60.00040
chr116254674962546845101.41promoterAHNAK1.0312.43 × 10−60.00041
chr7134928752134928863+21.03exonCALD10.992.80 × 10−60.00047
chr11117204261117204391+42.70exonTAGLN0.932.98 × 10−60.00047
chr5141969105141969137+6.875promoterRNF141.662.95 × 10−60.00047
chr10329578843295798058.70promoterITGB11.1923.16 × 10−60.00049
chr206369664663696657+16.99promoterRTEL1-TNFRSF6B;
TNFRSF6B
1.4053.96 × 10−60.00060
chr1109687817109687847+4.17promoterGSTM2;
GSTM1
5.4764.47 × 10−60.00067
chr167135872371358731+93.20promoterCALB22.514.96 × 10−60.00073
chr23023170930231716+7.77promoterLBH1.245.44 × 10−60.00078
chr221669535721669537052.35promoterIGFBP52.435.64 × 10−60.00079
chr22163720532163720786.62promoterMARCHF41.085.83 × 10−60.00080
chr17800492078004954+30.40promoterDNAJB4;
GIPC2
1.315.90 × 10−60.00080
chr9116153791116153813+44.29214promoterPAPPA1.596.49 × 10−60.00086
# TPM—tags per million; * FDR—false discovery rate.
Table 6. Overlap of endometriosis GWAS SNPs within 2 kb genomic window from the identified CAGE tag clusters.
Table 6. Overlap of endometriosis GWAS SNPs within 2 kb genomic window from the identified CAGE tag clusters.
ChrBPSNPp-ValueAlleleGene Region
chr122051787rs7251584.88 × 10−16aCDC42
chr122052387rs37544964.99 × 10−16aCDC42
chr2215433073rs12502448.93 × 10−8cFN1
chr7137345599rs1613354.67 × 10−6tnext to PTN
chr1130322210rs38584295.45 × 10−8tXR_931152.2
chr1130323044rs40715585.62 × 10−8tnear ARL14EP
chr1130323178rs40715595.60 × 10−8tARL14EP
chr1295216444rs73108338.07 × 10−9aFGD6
chr1295217365rs65386173.60 × 10−6tFGD6
chr1295217409rs65386186.68 × 10−9tFGD6
chr1744898335rs356531929.15 × 10−6aEFTUD2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Marla, S.; Mortlock, S.; Yoon, S.; Crawford, J.; Andersen, S.; Mueller, M.D.; McKinnon, B.; Nguyen, Q.; Montgomery, G.W. Global Analysis of Transcription Start Sites and Enhancers in Endometrial Stromal Cells and Differences Associated with Endometriosis. Cells 2023, 12, 1736. https://doi.org/10.3390/cells12131736

AMA Style

Marla S, Mortlock S, Yoon S, Crawford J, Andersen S, Mueller MD, McKinnon B, Nguyen Q, Montgomery GW. Global Analysis of Transcription Start Sites and Enhancers in Endometrial Stromal Cells and Differences Associated with Endometriosis. Cells. 2023; 12(13):1736. https://doi.org/10.3390/cells12131736

Chicago/Turabian Style

Marla, Sushma, Sally Mortlock, Sohye Yoon, Joanna Crawford, Stacey Andersen, Michael D. Mueller, Brett McKinnon, Quan Nguyen, and Grant W. Montgomery. 2023. "Global Analysis of Transcription Start Sites and Enhancers in Endometrial Stromal Cells and Differences Associated with Endometriosis" Cells 12, no. 13: 1736. https://doi.org/10.3390/cells12131736

APA Style

Marla, S., Mortlock, S., Yoon, S., Crawford, J., Andersen, S., Mueller, M. D., McKinnon, B., Nguyen, Q., & Montgomery, G. W. (2023). Global Analysis of Transcription Start Sites and Enhancers in Endometrial Stromal Cells and Differences Associated with Endometriosis. Cells, 12(13), 1736. https://doi.org/10.3390/cells12131736

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop