Gene Expression of the Endocannabinoid System in Endometrium Through Menstrual Cycle

doi:10.21203/rs.3.rs-1084037/v1

Download PDF

Research Article

Gene Expression of the Endocannabinoid System in Endometrium Through Menstrual Cycle

https://doi.org/10.21203/rs.3.rs-1084037/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Endocannabinoids mediate cellular functions and their activity is controlled by a complex system of enzymes, membrane receptors and transport molecules. Endocannabinoids are present in endometrium, a cyclical regenerative tissue requiring tightly regulated cellular mechanisms for maturation. The objective of this study was to chart the gene expression of key elements involved in the endocannabinoid system across the menstrual cycle. RNA was isolated from endometrial tissue and genome-wide gene expression datasets were generated using RNA-sequencing. An a priori set of 70 genes associated with endocannabinoid system were selected from published literature. Gene expression across the menstrual cycle was analyzed using a moderated t test, corrected for multiple testing with Bonferroni’s method. A total of 40 of the 70 genes were present in >90% of the samples, and significant differential gene expression identified for 29 genes. We identified 4 distinct regulation patterns for synthesizing enzymes, as well as a distinct regulation pattern for degradations and transporting enzymes. This study charts the expression of endometrial endocannabinoid system genes across the menstrual cycle. Altered expression of genes that control endocannabinoid may allow fine control over endocannabinoid concentrations and their influence on cellular function, maturation and differentiation as the endometrium matures through the menstrual cycle.

Physiology

Obstetrics & Gynecology

Endometrium

endocannabinoids

menstrual cycle

phospholipase

FABP (fatty acid-binding protein)

Endocannabinoids are small lipid-based molecules that play critical roles in regulating cellular functions. The most well characterised endocannabinoids are anandamide (AEA) and 2-arachdonoylglycerol (2-AG). Their production, concentration and cellular influence is mediated by a complex interaction of proteins that regulate their synthesis, transport, membrane binding, trafficking and degradation ¹. Collectively endocannabinoids and these proteins are termed the endocannabinoid system (ECS) ². Altered expression of any ECS component can lead to variation in endocannabinoid activity and alter their impact of cellular function.

Endocannabinoids are found in many organs and tissues including the endometrium and influence cell migration ^{3, 4}, proliferation ⁵, survival ⁶ inflammation ^{7, 8} and cellular differentiation ⁹. The endometrium lines the uterus and is comprised of epithelial glandular structures, vascularised stroma and infiltrating immune cells. It is unique in that it undergoes cyclical growth, regeneration and shedding each month. Endometrial maturation during the menstrual cycle is controlled by female sex hormones, but also requires tightly regulated cellular responses to facilitate the cellular and structural changes. The proliferative phase is dominated by estrogen-driven cellular proliferation, vascularisation and immune stimulation while the secretory phase is dominated by progesterone mediated stromal cell decidualisation, vascular remodelling and immune cell modulation ¹⁰. Regeneration of the endometrium following menses is mediated through adult stem and progenitor cells located in the basalis layer of the endometrium that proliferate, mature and differentiate across the menstrual cycle ^{11, 12}. The consistent cyclical shedding and regrowth creates the monthly potential for errors in cell replication that may underlie both permanent, or transient endometrial abnormalities.

Endocannabinoids are abundantly expressed in the endometrium and influence endometrial cellular function. Synthetic endocannabinoid, methandamide induces endometrial stromal cell migration ³, and AEA impairs both cellular proliferation and differentiation of an immortalised endometrial stromal cell line (St-T1b) and human decidual fibroblast from placenta, suggesting a crucial role in decidualisation ¹³. AEA was also shown to influence migration in the immortalized human endometrial epithelial cell line HEC-1B ⁴. Variations in concentrations of endocannabinoids may influence cellular proliferation and differentiation as the menstrual cycle progresses and the endometrium matures. Genes within the ECS are increasingly being recognised as targets to modulate endocannabinoid activity for multiple reproductive diseases including endometriosis, miscarriage, ectopic pregnancy, pre-eclampsia and endometrial cancer ^{14, 15}. How these are altered across the menstrual cycle is not yet clear.

We hypothesised that the ECS plays a role in the cellular function of the endometrium, that they are regulated across the menstrual cycle and that this may have implications for how the endometrium matures. The goal of this study was to define an a priori set of ECS genes from the literature, analyse their mRNA expression in human endometrial tissue and determine variations in their expression across the menstrual cycle.

ECS gene expression in the endometrium

Of the 70 genes selected from a priori list for examination we found only 1 gene with no discernible expression in any endometrial samples (PLA2G2E), 64 were expressed in at least 1 sample with >10 counts, 61 were expressed in 3 or more samples, 40 were expressed in greater than 90% of all samples (Supplementary table 1). Synthesizing enzymes accounted for 45 of the investigated genes, 29 (64.4%) of which were expressed in >90% of samples. Endocannabinoid degrading enzyme genes accounted for 12 of the investigated genes, of which 8 (66.7%) were expressed in >90% of samples. Genes encoding transport proteins accounted for 7 of the investigated genes of which 2 were consistently expressed in >90% of the samples. Finally, 6 genes for endocannabinoid membrane receptors were investigated and only one of these was expressed in >90% of the samples. We focused the subsequent analysis on these 40 genes expressed in >90% of the samples and their expression patterns across the menstrual cycle. The relative mean expression of each gene showed significant variation (Figure 1A) and a heatmap analysis charting gene expression against the menstrual cycle indicated a large degree of this variability could be explained by changes associated with cycle stage (Figure 1B).

ECS gene expression in women with and without endometriosis

Endometriosis is one of the most common endometrial pathologies occurring in 1 in 9 reproductive age women ^{16, 17} with suspected endometriosis the most common indication for surgery in this cohort. To determine whether endometrial pathologies, such as endometriosis, could influence endometrial ESC gene expression, we assessed the differences between women with and without endometriosis. Of the 206 samples collected 143 were confirmed endometriosis cases from visual inspection of the pelvis at laparoscopy and 63 were without evidence of endometriosis. Using the a priori set of selected genes and comparing the expression between women with and without endometriosis, no genes showed evidence of significant difference of expression when adjusted for multiple testing either compared together and taking into account menstrual cycle as a co-variate (Supplementary Table S2), or when comparing each cycle stage individually (Supplementary Table S3-S8). Previous studies on gene expression differences in these samples showed no genome-wide significance differences between women with and without endometriosis ^{18, 19}. In subsequent analysis of expression across the menstrual stage therefore, we assessed all samples together, regardless of endometriosis status.

Gene expression of enzymes responsible for endocannabinoid production

We performed 6 comparisons of gene expression between the 7 consecutive stages of the menstrual cycles, adjusting the p value to account for multiple testing of both multiple genes and menstrual cycle stage tests (Supplementary Table S9-14). We identified a significant difference in 21 genes between the 7 consecutive stages of the menstrual cycle. Using a heatmap analysis we compared the expression of these genes across all stages and grouped those that showed similar levels of expression and pattern of regulation (Figure 2), identifying four distinct patterns of interest.

The first group contained four genes with a high median gene expression that increased consistently across the menstrual cycle from the EP stage and peaked in either the ES (PLA2G4A, PLCL1, PLCB1) or MS stage (PLCB4) followed by a decline in expression in the LS (Figure 3A). Significant increases were observed for PLCL1 between the MP to LP stage (p = 0.004) and the LP to ES stage (p = 8.28 x 10⁻¹⁵), but was subsequently down regulated in the MS stage (p = 0.007). There was a significant increase in PLCB4 from the LP to ES stage (p = 0.007) with a further increase from the ES to MS stage (p = 0.003). A significant difference was recorded in the transition from the ES to MS stage for PLCB1 (p = 1.31 x 10⁻¹¹). A significant decrease was observed for PLA2G4A in the transition from the MS stage to the LS stage (p = 1.31 x 10⁻¹¹).

The second group of genes was characterised by a moderate mean expression with a gradual decrease across the cycle after initially peaking during the EP stage. An exception was observed with PLCB3 that showed a consistent expression across the cycle peaking in the MS stage (Figure 3B). Significant increases in expression were observed during the transition from the M to EP stage for PLCE1 (p = 0.004) and PLA2G12A (p = 0.001) that also showed a significant decrease from the LP to ES (p = 0.007), as well as the ES to MS stage (p = 1.44 x 10⁻⁸). Significant decreases were also observed for GDE1 from the LP to ES stage (p = 7.81 x 10⁻⁸) and for PLCB3 in transition from the MS to LS stage (p = 0.003).

The third group was characterised by medium mean expression and that decreased directly after the menstrual stage and then subsequently increased again after the ES stage (Figure 3C). Significant decreases in PLD1 were observed in the transition from the M to EP stage (p = 0.0005) and again from the LP to ES transition (p = 0.002). A significant increase was observed for PLCG2 (p = 0.005) and PLCD3 (p = 8.01 x10⁻⁵) during the LP to ES transition. PLD1 began to increase in the MS stage (p = 0.0003) continuing into the LS stage (p = 2.8 x 10⁻⁸). Both PLCG2 (p = 0.005) and PLCD3 (p = 1.44 x 10⁻⁸) continued their increase during the early to mid-secretory stage.

The final group was characterised by low mean expression that showed multiple and varied expression changes across the EP to ES stage (Figure 3D). Significant differences for the individual genes were observed between the M to EP stage for DAGLA that was increased (p = 0.003) and PTPN22 that was decreased (p = 0.003). There was minimal regulation across the proliferative stage, however in the transition from the LP to ES stage there was a decrease in GDE4 (p = 5.48 x10⁻⁵), PLCL2 (p = 0.007) and an increase in GDE7 (p = 2.93 x 10⁻¹¹), ABHD4 (p = 0.002), PLCH2 (p = 0.0001) and PLA2G6 (p = 0.007). Significant changes also occurred in the ES to MS stage with a reduction in PLCH2 (p = 2.19 x10⁻¹⁶), DAGL (p = 0.0004), GDE7 (p = 0.0009) and an increase in PLCB2 (p = 0.002) and PTPN22 (p = 0.001). NAPE-PLD decreased from the MS to LS (p = 0.003). Finally, PLCH2 (p = 0.003), PTPN22 (p = 0.0002) and GDE4 (p = 0.003) were all significantly increased in LS compared to the MS stage.

Gene expression encoding endocannabinoid transporters

We identified 7 genes known to encode fatty acid-binding proteins (FABPs) that transport endocannabinoids intracellularly. The largest changes were observed for FABP5 which was significantly upregulated from the LP to ES stage (p = 0.002) (Supplementary Table S12) with a subsequent down regulation from the MS to LS stage (p = 0.01) (Figure 4A) (Supplementary Table S14). FABP3 was significantly upregulated in the transition from the ES to MS stage (p = 0.0005) (Supplementary Table S13).

Gene expression of enzymes responsible for degradation of endocannabinoids

We identified expression of 8 of the 12 genes related to endocannabinoid degradation in endometrial tissue. Of these ALOX5 (p = 0.0005) and COX-2 (p = 0.008) were significantly down regulated and NAAA significantly upregulated in the transition from the M to the EP stage (p = 0.008) (Figure 4B) (Supplementary Table S9). In the transition from the LP to ES stage, FAAH was decreased (p = 0.002) and ABHD12 was increased (p = 0.001) (Figure 5) (Supplementary Table S12). In contrast in the ES to MS stage FAAH was upregulated (p = 0.0003) and ABHD12 downregulated (p = 0.004) (Supplementary Table S13). Other significant changes from the ES to the MS stage included; NAAA (p = 5.04 x10⁻⁷) and COX-2 (p = 1.15 x 10⁻⁶) that were downregulated and ALOX5 (p = 0.0002) that was upregulated.

Gene expression of endocannabinoid receptors

We assessed gene expression of 6 ECS membrane receptors across the menstrual cycle and between endometriosis status. Of these 6 genes we found that only one, TRPV1 was expressed in >90% of the samples. It had a low but consistent expression and was not significantly regulated across the menstrual cycle or influenced by endometriosis status. The expression of all other membrane receptors occurred in <90% of the samples. This suggests either a consistently low expression below read depth, or a binary on/off regulation. Based on the potential for binary regulation we determined the percentage of samples expressing membranes receptor genes in each menstrual stage (Figure 5). For both CNR1 and CNR2 the number of samples with positive expression was between 70-100% for the various stages of the cycle, decreasing through the proliferative stage with a subsequent rise during the secretory stage. The changes however did not reach significance. GPR55 expression occurred in approximately 70% of samples, also decreasing through the proliferative stage reaching a minimum during the ES stage and subsequently increasing during the MS to LS stages. TRPA1 expression was observed only between the MP and MS stage and only in a low number of samples. Statistical analysis revealed there was no significant difference for any endocannabinoid receptor across the different stages of the menstrual cycle. Furthermore, no difference was observed in the percentage of samples expressing any of the membrane receptors when comparing endometriosis status and adjusting for menstrual cycle phase.

Endocannabinoids, naturally synthesized lipids are produced in abundance in the female reproductive tract ²⁰. They are present in normal endometrium and significantly influence cellular maturation and function ²¹. Modulation of the ECS system occurs through enzymes regulating production, transport and degradation of the endocannabinoid signaling molecules. Understanding enzyme expression may provide insights into how endocannabinoid levels are modulated in endometrium. Utilizing genome-wide gene expression data of endometrial samples we performed a comprehensive analysis of ECS gene expression during the menstrual cycle. We identified significant regulation across the menstrual cycle for many genes and identified distinct patterns of regulation that strongly support alteration in endocannabinoid activity across the menstrual cycle that could be both a consequence of, or contribution to the changes that occur during cyclical endometrial maturation.

Enzymes capable of synthesizing endocannabinoid production were the most abundant in the endometrium, and also showed the greatest level of menstrual cycle regulation. The majority of gene expression changes occurred during the transition from the late proliferative to early secretory stage and continued across the secretory stage. The transition from proliferative to secretory stage is characterized by ovulation and an increase in the production of progesterone. Progesterone levels also fluctuate during the secretory stage by increasing up to the mid secretory stage and decreasing during the late secretory stage ¹⁰. Gene expression changes across the cycle suggest an important role of steroid hormones in ECS modulation, in particular progesterone. Previous studies have reported progesterone is involved in the maintenance of endocannabinoid levels of the endometrium ²². Although in contrast, one study reported no correlation between plasma levels of AEA and progesterone in normal cycling women ²³.

The most abundant endocannabinoids in the endometrium, AEA and 2-AG, are predominantly synthesized via NAPE-PLD and DAGL respectively ¹, representing two of the most important enzymes in the ECS. We found both enzymes were downregulated during the transition from menstrual to early proliferative stage and then subsequently upregulated from early to mid-secretory and again from the mid to late secretory stage. This suggests an increased capacity for endocannabinoid production as the menstrual cycle progresses. This increased capacity may not be linked to progesterone as it decreases towards the end of the cycle, a period where endocannabinoid production capacity is still increasing. A previous study reported NAPE-PLD immunoreactivity was increased in the menstrual, early proliferative and late secretory glands with its lowest levels in the early secretory phase ²⁴. In a mouse model progesterone was reported to decrease NAPE-PLD expression ²⁵. The increase in these enzymes and the potential increase in the endocannabinoids may have implications for cellular function and differentiation during the menstrual cycle in humans and rodents.

A number of other enzymes can regulate endocannabinoid production through alternate pathways. One such group includes the many isoforms of phospholipase C (PLC). PLCs selectively catalyze the degradation of phosphatidylinositol 4,5 bisphosphate resulting in soluble inositol-1,4,5-triphosphate (IP3) and membrane delimited 1,2 diacylglycerol (DAG) ²⁶. DAG is an important precursor for 2-AG ²⁷. We observed multiple, but distinct patterns of PLC regulation across the menstrual cycle, particularly during the transition from the proliferative to secretory stage. This regulation would impact DAG expression subsequently influencing 2-AG production. Whether each of these isoforms has similar efficacy in the metabolism of IP3 and the subsequent production of DAG in the endometrium is unclear from this study.

Endocannabinoids continue to exert their activity until degraded. A number of enzymes responsible for the degradation of endocannabinoids, including FAAH, MGLL and ABHD12 were also significantly downregulated during the transition from the late proliferative to early secretory stage. A down regulation of these enzymes would potentiate the ability of endocannabinoids to continue signaling. Of these enzymes FAAH is of particular interest. FAAH is a principal enzyme in the degradation of AEA and previous evidence suggests a relationship to progesterone ^{28, 29}. FAAH has previously been implicated in the survival of endometrial stromal cells from ectopic endometriosis lesions ³⁰ and has received significant interest as a pharmacological target, modulating pain through both peripheral and central mechanisms ³¹. Physiological concentrations of progesterone stimulated FAAH activity in human lymphocytes, decreasing AEA levels ^{28, 32, 33}, while in a mouse model progesterone was reported to decrease uterine FAAH activity ³⁴. Progesterone levels and FAAH expression have been correlated during the menstrual cycle ³⁵, in agreement with the finding that progesterone up-regulates FAAH gene expression ^{32, 33}. Combined with the regulation of the PLCs during this period, this suggests a regulation of endocannabinoid turnover during the transition between menstrual cycle stages that is tightly regulated.

Additional factors, other than their production and inhibition, can influence endocannabinoids concentrations and their activity. Intracellular degradation relies on lipophilic endocannabinoid reaching the intracellular enzymes. One of the most significant changes we observed in this gene set was the regulation of FABP5 across the menstrual cycle. FABP5 is a fatty acid binding protein that transports endocannabinoids through the aqueous cellular environment, preferentially directing it towards FAAH ³⁶. FABP5 was also significantly downregulated during the transition from the late proliferative to the early secretory stage. A reduction in the transport towards the degrading enzyme could potentiate the presence of the endocannabinoids during this period and may contribute to cellular and inflammatory changes ^{10, 37} that enables the regeneration of the endometrium in the early and mid-secretory phase in non-conceptive cycles.

Although we included an extensive list of genes encoding endocannabinoid receptors, including CNR1, CNR2, TRPV1, TRPA1, GPR55, and GPR119 we found that only one of these (TRPV1) was consistently expressed in >90% of the samples. Analysis of gene expression for the remaining receptors as a binary variable found changes across the menstrual cycle. Although we did not find any that reached significance after multiple testing, this may be related to the power of the sample size, or the complex milieu of individual cell types and maturing cell states that exist within the endometrium at any point in time. The results showed non-significant changes in number of samples expressing CNR1, CNR2 and GPR55 during the late proliferative to early secretory period. Previous studies on CB1 expression (the protein encoded by the CNR1 gene) in the epithelial glands found no significant regulation across the menstrual cycle ²⁴, and CB2 expression in stromal cells was similar between the proliferative and the secretory phase ⁶. Previous studies have also reported that progesterone exerted minimal effects on CB1 expression in lymphocytes ^{28, 32, 33}. It is therefore likely that the lack of significant changes in the expression of these receptors across samples and across the menstrual cycle indicates they may not be significantly regulated by reproductive hormones in endometrial tissue. The ability of these genes to be switched on or off may reflect the selective use in patients, depended on local requirements and provide the ability to fine tune the activity of the ECS.

From our study cohort, 143 (69.4%) were diagnosed with endometriosis. Previous studies on the ECS in women with endometriosis have shown conflicting results. One reported no significant difference in NAPE-PLD and FAAH expression in the endometrium of patients compared to women without endometriosis throughout the menstrual cycle ³⁸. Another study reported no difference in CB1 protein expression during the proliferative phase between patients with and without endometriosis ³⁸, whereas lower levels of CNR1 and the protein CB1 have been reported in endometrial tissue from women with endometriosis compared to controls regardless of the cycle phase ⁷. No difference in TRPV1 expression between women with and without endometriosis throughout the menstrual cycle has also been reported ³⁸. In this current study, when accounting for menstrual stage and multiple testing, we found no significant difference in any of the ECS genes investigated between women with and without endometriosis. Our data therefore correlates with previously published data that shows the ECS is not significantly dysregulated in eutopic endometrial tissue from patients with endometriosis.

In summary, we have utilized genome wide gene expression to investigate a specific predefined gene set based on their previously identified involvement in the ECS. Targeted analysis of a defined gene set will reduce the burden of multiple testing in gene expression data of modest size providing the potential to uncover subtle variations. We detected differences in the gene expression across the menstrual cycle that are most pronounced during the transition from the proliferative to secretory stage and reflect a potential to dynamically modulate endocannabinoid concentrations during this period. As the activity of endocannabinoids are rapid and short-lived, we speculate that the altered regulation of expression for many components of the ECS across the menstrual cycle provides a quick acting, fine tuning potential to modulate endocannabinoid expression as needed, particularly during the transition period in which the endometrium undergoes major structural rearrangement. Of particular interest, we noted the large changes in a number of PLC enzymes, and the significant variation in FAAH and FABP5 that work in concert to deactivate endocannabinoid activity and are already receiving considerable attention as drug targets for pain, inflammation and cancer ^{39, 40}. The important ability of endocannabinoid receptors to be regulated may also hold particular interest for the potential of selective targeting of patients and the personalisation of medication.

Sample collection

Samples were collected as described previously (n=206) ¹⁸ and inclusion criteria were European ancestry and within reproductive age (18-49, mean = 31.85). Exclusion criteria were the use of hormonal medication within 3 months prior to surgery, abnormalities in histopathological examination, ambiguous disease status or menstrual cycle stage. Informed consent was obtained from all subjects. A histological assessment was performed on formalin fixed paraffin embedded tissue with each sample assigned to 1 of 7 possible menstrual stages following the Noyes criteria ⁴¹: menstrual (M) (n =14), early proliferative (EP) (n = 5), mid proliferative (MP) (n =72), late proliferative (LP) (n =22), early secretory (ES) (n =31), mid secretory (MS) (n=41) and late secretory (LS) (n=21). If endometriosis was observed during laparoscopic surgery lesions were excised and sent to pathology for histopathological confirmation (n = 143 endometriosis cases and n = 63 non-endometriosis controls). The study was approved by the Human Research Ethics Committee (HREC) of the Royal Women’s Hospital (Projects 11-24 and 16-43 and the Melbourne IVF Project 05-11) and the University of Queensland HREC (2016000746), and all experiments were performed in accordance with relevant guidelines and regulations.

RNA extraction and Sequencing

Total RNA was isolated from endometrial biopsies stored in RNAlater using the Allprep DNA/RNA Mini Kit (Qiagen, USA) as per the manufacturer’s instructions and isolated RNA samples were treated with the Turbo DNA-free kit (Thermo Fischer Scientific, USA). RNA quality was confirmed using the Agilent Bioanalyzer 2100 (Agilent Technologies, USA) with RNA integrity number (RIN) cut off values set at above 8 for inclusion and final RNA concentrations determined by the Nanodrop ND -6000 (Thermo Fisher Scientific, USA). Stranded RNA-sequencing (RNA-seq) libraries were prepared with the Illumina TruSeq Stranded Total RNA Gold protocol incorporating ribosomal depletion (Illumina, USA). The resulting libraries were pooled and sequenced with 75bp paired-end reads on the Illumina HiSeq 4000 to a mean depth of 37,490,673 for 178 samples and with 120bp paired-end reads on an Illumina Hi Seq 2000 (Illumina, USA) for a mean depth of 40,818,062 reads for 28 samples.

Preparation of RNA-seq data

The quality of raw RNA-seq reads were confirmed with FastQC v0.11.7 ⁴² and MultiQC v1.6 ⁴³. Low quality reads, or reads containing HiSeq Illumina adapter sequences were trimmed using Trimmomatic v0.36 ⁴⁴. The resulting trimmed reads were aligned to a reference assembly (Ensembl Homo sapiens GRCh38 release 91) with HISAT2 v2.0.5 and the transcript assembly performed with StringTie v1.3.1 ^{45, 46} and each read mapping to a known transcript counted. Transcript, exon and expression matrices in Fragments Per Kilobase of transcript per million mapped reads (FPKM) were determined with StringTie counts for each individual. Prior to normalisation of RNA-seq counts lowly expressed genes (counts per million (CPM) < 0.22) and expressed in <90% of samples were removed. Gene counts were normalised for composition bias and total raw reads using TMM ^47–49 in edgeR R package v 3.22.3 ⁵⁰. The normalized counts were converted to CPM and log2 transformed (log2-CPM).

Selection of the ECS genes

A priori gene selection was performed through a literature search using the search term ‘endocannabinoid’, and combinations of ‘pathway’, ‘synthesis’, ‘metabolism’, ‘transport’ or ‘endometrium’. All genes relevant to the control of the ECS were catalogued by ENSEMBL gene expression ID and gene expression extracted from the curated database. ECS genes of interest were split into four categories based on their main function in the ECS including (i) synthesizing enzymes, (ii) receptors, (iii) transporters, and (iv) degradation enzymes. In total 70 potential genes known to play a role in the ECS were assessed (Supplementary Table 1).

Gene expression analysis

Differential gene expression was investigated both across the menstrual cycle and in women with and without endometriosis. Only genes with CPM>0.22 in at least 90% of samples were analysed for differential expression. Seven differential gene expressions comparisons were performed across the menstrual cycle: (i) between M and EP, (ii) between EP and MP, (iii) between MP and LP, (iv) between LP and ES, (v) between ES and MS, (vi) between MS and LS and (vii) between MP and MS. Seven differential gene expressions comparisons were performed between women with and without endometriosis: (i) including stage of menstrual cycle as a covariate, (ii) within M, (iii) within MP, (iv) within LP, (v) within ES, (vi) within MS and (vii) within LS. Batch effects (flow-cell and lane) were fitted as covariates in all models. The voom function in limma ⁵¹ was used to transform the normalised counts prior to using the eBayes method to test between groups. P-values were adjusted for the number of genes tested and for the comparison of multiple cycle stage comparisons using Benjamini–Hochberg with a significance threshold of 0.05. The data underlying this article are available on the article and in its online supplementary material.

In the case of membrane receptors where positive gene expression was consistently <90% of the samples, we explored gene expression as a binary function of expressed/not expressed and performed a logistic regression analysis as described previously ¹⁸ to determine if there was as significant variation in the number of samples that show a positive expression of these genes across the menstrual cycle.

Author contributions

KT, AAA, GWM and BM were responsible for the conception and design of the present study. PAWR, SJH, JFD and WTT acquired the clinical and surgical data. PAWR, SJH, JFD and WTT acquired the histological and RNA data. SM was responsible for statistical analysis. All authors drafted the work and revised it critically for important intellectual content. All authors have read and approved the final manuscript.

Tanaka, K. et al. The role of the endocannabinoid system in aetiopathogenesis of endometriosis: A potential therapeutic target. Eur J Obstet Gynecol Reprod Biol, 244, 87–94 (2020).
Lu, H. C. & Mackie, K. An Introduction to the Endogenous Cannabinoid System. Biol Psychiatry, 79 (7), 516–25 (2016).
Gentilini, D. et al. Endocannabinoid system regulates migration of endometrial stromal cells via cannabinoid receptor 1 through the activation of PI3K and ERK1/2 pathways. Fertil Steril, 93 (8), 2588–93 (2010).
McHugh, D., Page, J., Dunn, E. & Bradshaw, H. B. Delta(9) -Tetrahydrocannabinol and N-arachidonyl glycine are full agonists at GPR18 receptors and induce migration in human endometrial HEC-1B cells. Br J Pharmacol, 165 (8), 2414–24 (2012).
Leconte, M. et al. Antiproliferative effects of cannabinoid agonists on deep infiltrating endometriosis. Am J Pathol, 177 (6), 2963–70 (2010).
Bilgic, E. et al. Endocannabinoids modulate apoptosis in endometriosis and adenomyosis. Acta Histochem, 119 (5), 523–32 (2017).
Resuehr, D., Glore, D. R., Taylor, H. S., Bruner-Tran, K. L. & Osteen, K. G. Progesterone-dependent regulation of endometrial cannabinoid receptor type 1 (CB1-R) expression is disrupted in women with endometriosis and in isolated stromal cells exposed to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Fertil Steril, 98 (4), 948–561 (2012).
Iuvone, T. et al. Selective CB2 up-regulation in women affected by endometrial inflammation. J Cell Mol Med, 12 (2), 661–70 (2008).
Galve-Roperh, I. et al. Cannabinoid receptor signaling in progenitor/stem cell proliferation and differentiation. Progress in Lipid Research, 52 (4), 633–50 (2013).
Critchley, H. O. D., Maybin, J. A., Armstrong, G. M. & Williams, A. R. W. Physiology of the Endometrium and Regulation of Menstruation. Physiol Rev, 100 (3), 1149–79 (2020).
Gargett, C. E., Schwab, K. E. & Deane, J. A. Endometrial stem/progenitor cells: the first 10 years. Hum Reprod Update, 22 (2), 137–63 (2016).
Cousins, F. L., Pandoy, R., Jin, S. & Gargett, C. E. The Elusive Endometrial Epithelial Stem/Progenitor Cells. Front Cell Dev Biol, 9, 640319 (2021).
Almada, M. et al. The endocannabinoid anandamide impairs in vitro decidualization of human cells., 152 (4), 351–61 (2016).
Stasiulewicz, A., Znajdek, K., Grudzien, M., Pawinski, T. & Sulkowska, A. J. I. A Guide to Targeting the Endocannabinoid System in Drug Design.Int J Mol Sci. 2020;21(8)
Maia, J., Fonseca, B., Teixeira, N. & Correia-da-Silva, G. The fundamental role of the endocannabinoid system in endometrium and placenta: implications in pathophysiological aspects of uterine and pregnancy disorders. Hum. Reprod. Update, 26 (4), 586–602 (2020).
Giudice, L. C., Kao, L. C. & EndometriosisLancet. 2004; 364(9447):1789–99
Rowlands, I. J. et al. Prevalence and incidence of endometriosis in Australian women: a data linkage cohort study. BJOG, 128 (4), 657–65 (2021).
Fung, J. N. et al. Genetic regulation of disease risk and endometrial gene expression highlights potential target genes for endometriosis and polycystic ovarian syndrome. Sci. Rep, 8 (1), 11424 (2018).
Mortlock, S. et al. Tissue specific regulation of transcription in endometrium and association with disease. Hum. Reprod, 35 (2), 377–93 (2020).
Guindon, J. & Hohmann, A. G. The endocannabinoid system and pain. CNS Neurol Disord Drug Targets, 8 (6), 403–21 (2009).
Di Blasio, A. M., Vignali, M. & Gentilini, D. The endocannabinoid pathway and the female reproductive organs. J Mol Endocrinol, 50 (1), R1–9 (2013).
Karasu, T., Marczylo, T. H., Maccarrone, M. & Konje, J. C. The role of sex steroid hormones, cytokines and the endocannabinoid system in female fertility. Hum Reprod Update, 17 (3), 347–61 (2011).
El-Talatini, M. R., Taylor, A. H. & Konje, J. C. The relationship between plasma levels of the endocannabinoid, anandamide, sex steroids, and gonadotrophins during the menstrual cycle. Fertil Steril, 93 (6), 1989–96 (2010).
Taylor, A. H., Abbas, M. S., Habiba, M. A. & Konje, J. C. Histomorphometric evaluation of cannabinoid receptor and anandamide modulating enzyme expression in the human endometrium through the menstrual cycle. Histochem Cell Biol, 133 (5), 557–65 (2010).
Guo, Y. et al. N-acylphosphatidylethanolamine-hydrolyzing phospholipase D is an important determinant of uterine anandamide levels during implantation. J Biol Chem, 280 (25), 23429–32 (2005).
Rajala, R. V. S. & Anderson, R. E. Focus on molecules: phosphatidylinositol-4,5-bisphosphate (PIP2). Exp Eye Res, 91 (3), 324–5 (2010).
Murataeva, N., Straiker, A. & Mackie, K. Parsing the players: 2-arachidonoylglycerol synthesis and degradation in the CNS. British journal of pharmacology, 171 (6), 1379–91 (2014).
Maccarrone, M. et al. Progesterone up-regulates anandamide hydrolase in human lymphocytes: role of cytokines and implications for fertility. J Immunol, 166 (12), 7183–9 (2001).
Gasperi, V., Fezza, F., Spagnuolo, P., Pasquariello, N. & Maccarrone, M. Further insights into the regulation of human FAAH by progesterone and leptin implications for endogenous levels of anandamide and apoptosis of immune and neuronal cells., 26 (5), 811–7 (2005).
Pavone, M. E. et al. Altered retinoid uptake and action contributes to cell survival in endometriosis. J Clin Endocrinol Metab, 95 (11), E300–9 (2010).
Peng, X. et al. Fatty-acid-binding protein inhibition produces analgesic effects through peripheral and central mechanisms. Mol Pain, 13, 1744806917697007 (2017).
Maccarrone, M., Di Rienzo, M., Finazzi-Agro, A. & Rossi, A. Leptin activates the anandamide hydrolase promoter in human T lymphocytes through STAT3. J Biol Chem, 278 (15), 13318–24 (2003).
Maccarrone, M., Bari, M., Di Rienzo, M., Finazzi-Agro, A. & Rossi, A. Progesterone activates fatty acid amide hydrolase (FAAH) promoter in human T lymphocytes through the transcription factor Ikaros. Evidence for a synergistic effect of leptin. J Biol Chem, 278 (35), 32726–32 (2003).
MacCarrone, M. et al. Down-regulation of anandamide hydrolase in mouse uterus by sex hormones. Eur J Biochem, 267 (10), 2991–7 (2000).
Lazzarin, N. et al. Fluctuations of fatty acid amide hydrolase and anandamide levels during the human ovulatory cycle. Gynecol Endocrinol, 18 (4), 212–8 (2004).
Kaczocha, M., Glaser, S. T. & Deutsch, D. G. Identification of intracellular carriers for the endocannabinoid anandamide. Proc Natl Acad Sci U S A, 106 (15), 6375–80 (2009).
Hall, O. J. & Klein, S. L. Progesterone-based compounds affect immune responses and susceptibility to infections at diverse mucosal sites. Mucosal Immunol, 10 (5), 1097–107 (2017).
Sanchez, A. M. et al. Elevated Systemic Levels of Endocannabinoids and Related Mediators Across the Menstrual Cycle in Women With Endometriosis. Reprod Sci, 23 (8), 1071–9 (2016).
Schlosburg, J. E., Kinsey, S. G. & Lichtman, A. H. Targeting fatty acid amide hydrolase (FAAH) to treat pain and inflammation. Aaps j, 11 (1), 39–44 (2009).
Brunetti, L., Loiodice, F., Piemontese, L., Tortorella, P. & Laghezza, A. New Approaches to Cancer Therapy: Combining Fatty Acid Amide Hydrolase (FAAH) Inhibition with Peroxisome Proliferator-Activated Receptors (PPARs) Activation. Journal of Medicinal Chemistry, 62 (24), 10995–1003 (2019).
Noyes, R. W., Hertig, A. T. & Rock, J. Dating the Endometrial Biopsy. Fertility and Sterility, 1 (1), 3–25 (1950).
Andrews FastQC: a quality control tool for high throughput sequence data. Available online at http://www.bioinformatics.babraham.ac.uk/projects/fastqc. 2010
Ewels, P., Magnusson, M., Lundin, S. & Käller, M. MultiQC: summarize analysis results for multiple tools and samples in a single report., 32 (19), 3047–8 (2016).
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data., 30 (15), 2114–20 (2014).
Pertea, M., Kim, D., Pertea, G., Leek, J. & Salzberg, S. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc, 11 (9), 1650–67 (2016).
Pertea, M. et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol, 33 (3), 290–5 (2015).
Fadista, J. et al. Global genomic and transcriptomic analysis of human pancreatic islets reveals novel genes influencing glucose metabolism. Proceedings of the National Academy of Sciences. 2014;111(38):13924-9
Seo, M. et al. RNA-seq analysis for detecting quantitative trait-associated genes. Sci. Rep, 6, 24375 (2016).
Taneera, J. et al. Identification of novel genes for glucose metabolism based upon expression pattern in human islets and effect on insulin secretion and glycemia. Hum. Mol. Genet, 24 (7), 1945–5 (2015).
Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data., 26 (1), 139–40 (2010).
Law, C. W., Chen, Y., Shi, W. & Smyth, G. K. voom: Precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol, 15 (2), R29 (2014).

No competing interests reported.

Download PDF

Editorial decision: Major revision
30 Mar, 2022
Reviews received at journal
17 Mar, 2022
Reviewers agreed at journal
08 Mar, 2022
Reviews received at journal
28 Dec, 2021
Reviewers agreed at journal
23 Dec, 2021
Reviewers invited by journal
22 Dec, 2021
Editor assigned by journal
22 Dec, 2021
Editor invited by journal
25 Nov, 2021
Submission checks completed at journal
25 Nov, 2021
First submitted to journal
15 Nov, 2021

You are reading this latest preprint version

Gene Expression of the Endocannabinoid System in Endometrium Through Menstrual Cycle

Status:

Version 1

Abstract

Figures

Introduction

Results

ECS gene expression in the endometrium

ECS gene expression in women with and without endometriosis

Gene expression of enzymes responsible for endocannabinoid production

Gene expression encoding endocannabinoid transporters

Gene expression of enzymes responsible for degradation of endocannabinoids

Gene expression of endocannabinoid receptors

Discussion

Methods

Sample collection

RNA extraction and Sequencing

Preparation of RNA-seq data

Selection of the ECS genes

Gene expression analysis

Declarations

Author contributions

References

Additional Declarations

Supplementary Files

Status:

Version 1