Next Article in Journal
Synthesis, Structure, and Antiproliferative Action of 2-Pyridyl Urea-Based Cu(II) Complexes
Next Article in Special Issue
SARS-CoV-2 Infection, Sex-Related Differences, and a Possible Personalized Treatment Approach with Valproic Acid: A Review
Previous Article in Journal
Extracellular Vesicles—New Players in Cell-to-Cell Communication in Gestational Diabetes Mellitus
Previous Article in Special Issue
Sex, Allergic Diseases and Omalizumab
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Epigenetic Effects of Gender-Affirming Hormone Treatment: A Pilot Study of the ESR2 Promoter’s Methylation in AFAB People

1
Laboratory of Seminology—Sperm Bank “Loredana Gandini”, Department of Experimental Medicine, University of Rome “Sapienza”, 00161 Rome, Italy
2
Department of Psychological Health and Territorial Sciences, School of Medicine and Health Sciences, “G. d’Annunzio” University of Chieti-Pescara, 66100 Chieti, Italy
3
Unit of Molecular Genetics, Center for Advanced Studies and Technology (CAST), “G. d’Annunzio” University of Chieti-Pescara, 66100 Chieti, Italy
4
Department of Experimental and Clinical Biomedical Sciences “Mario Serio”, University of Florence, 50139 Florence, Italy
*
Author to whom correspondence should be addressed.
Biomedicines 2022, 10(2), 459; https://doi.org/10.3390/biomedicines10020459
Submission received: 15 January 2022 / Revised: 9 February 2022 / Accepted: 14 February 2022 / Published: 16 February 2022
(This article belongs to the Special Issue Gender Medicine and Pharmacology)

Abstract

:
Virilization of gender-incongruent subjects to whom were assigned the female gender at birth (AFAB) is achieved through testosterone administration. Inter-individual differences in the timing and acquisition of phenotypic characteristics, even if the same hormone preparations and regimens are used, are frequently observed. Polymorphisms of sex hormone receptors and methylation of their gene promoters, as well of several imprinted genes as H19, may underlie the differential response to treatment. Thus, the aim of this study was to examine the possible relationship between the CpG methylation profile of the estrogen receptor 2 gene (ESR2) and H19 promoters and their influence on phenotype modifications in a cohort of AFAB people at baseline (T0) and after 6 mo (T6) and 12 mo (T12) of testosterone therapy (testosterone enanthate, 250 mg i.m. every 28 d). A total of 13 AFAB subjects (mean age 29.3 ± 12.6) were recruited. The percentage of methylation of the ESR2 promoter significantly increased at T6 (adj. p = 0.001) and T12 (adj. p = 0.05), while no difference was detected for H19 (p = 0.237). Methylation levels were not associated with androgen receptor (AR)/estrogen receptor beta (ERβ) polymorphisms nor hormone levels at baseline and after six months of treatment. On the other hand, total testosterone level and patient age resulted in being significantly associated with ESR2 methylation after twelve months of treatment. Finally, the difference in ESR2 promoter methylation between T6 and baseline was significantly associated with the number of CA repeats of the ERβ receptor, adjusted vs. all considered variables (R2 = 0.62, adj. R2 = 0.35). No associations were found with CAG repeats of the AR, age, and estradiol and testosterone levels. Despite the small sample size, we can hypothesize that treatment with exogenous testosterone can modify the ESR2 methylation pattern. Our data also indicated that epigenetic changes may be regulated, suggesting that the modulation of estrogen signaling is relevant shortly after the beginning of the treatment up to T6, with no further significant modification at T12. Furthermore, estrogen receptor methylation appears to be associated with the age of the subjects and exogenous testosterone administration, representing a marker of androgenic treatment. Nonetheless, it will be necessary to increase the number of subjects to evaluate how epigenetic regulation might play a relevant role in the modulation of phenotypical changes after testosterone treatment.

1. Introduction

Gender-affirming hormonal therapy (GAHT) is critical for phenotypical and psychological transition in adults with gender incongruence (GI), which is an enduring and marked incongruence between the individual’s desired gender and the assigned sex at birth. While this definition includes a broad spectrum of transgender people (stretching between binary and non-binary subjects), the current scientific literature generally refers to assigned females at birth (AFAB) and assigned males at birth (AMAB) subjects [1].
Transgender people are often considered a psychologically vulnerable group compared to the general population, partly due to distress deriving from GI itself and partly to a non-inclusive cultural background, which might also include difficulties in accessing general healthcare and GAHT [1,2]. GAHT, in particular, is a paramount tool to alleviate the incongruence between biological and experienced gender and, therefore, any psychological burden arising because of it. Moreover, physical changes induced by a hormone gender-affirming treatment are associated with an improvement of the subjects’ quality of life and self-esteem, reducing anxiety, depression, and social distress [3].
This is achieved through the reduction of endogenous sex hormones and their replacement with exogenous hormones, allowing the progressive substitution of the incongruent patient’s secondary sexual characteristics with the desired ones. In virilizing therapy for incongruent subjects who have been assigned the female gender at birth (AFAB people), the treatment of choice relies on one of the many testosterone preparations available, following the principles of hormone replacement therapy in hypogonadal patients [4]. Cessation of menses, clitoridomegaly, increased muscle/fat mass ratio, facial and body hair growth, voice deepening, and other desired features usually appear within months to years. Nonetheless, it has been observed that inter-individual differences in the timing and acquisition of phenotypic characteristics are frequent, even in subjects who undergo the same hormone preparation and regimen. There are several explanations for this phenomenon, many of which are linked to inter-individual differences in hormone signaling due to the polymorphisms of sex hormone receptors [5]. It is likely that hormone receptor activity might be modulated by both congenital and acquired factors. Unfortunately, literature data on genetic and epigenetic contributions in this setting are lacking. Nonetheless, gene promoters are capable of being activated or inhibited by physiological changes or by the action of exogenous compounds with estrogenic activity [6]. The most frequently characterized epigenetic modification is DNA methylation, which involves inherent and acquired gene transcription changes, independent of the DNA sequence.
Scant data are available for the androgen receptor (AR), and a few authors have investigated the occurrence of modifications of the estrogen receptor 1 gene (ESR1) promoter methylation after GAHT [7,8]. The functional biological meaning of the induced methylation pattern is still unknown, but it is possible to speculate that AR and ESR1 genes’ methylation could play a role in modulating the induction of desired phenotypical changes, as well as being a putative marker of hormonal treatment. On the other hand, no data are available on the estrogen receptor 2 (ESR2) gene.
Moreover, H19 is a paternal imprinted gene whose hypomethylation is associated with the Silver Russell syndrome and other anomalies (pre-natal and post-natal growth deficit, insulin resistance, polycystic ovarian syndrome, and precocious pubarche) [9]. Nonetheless, DNA methylation of H19 has also been associated with pubic hair onset and genital or breast development in adolescent boys and girls, respectively [10]. Thus, the H19 promoter methylation could be potentially related to a differential development of secondary sex characteristics after GAHT.
Thus, the main purpose of this study was to examine the possible relationship between the CpG methylation profile of ESR2 and H19 promoters and their possible influence on acquired changes in a cohort of AFAB people at baseline (T0) and after 6 mo (T6) and 12 mo (T12) of testosterone therapy.

2. Materials and Methods

We prospectively recruited consecutive AFAB people referring to the Gender Incongruence Ambulatory of the Department of Experimental Medicine (AOU Policlinico Umberto I—“Sapienza” University of Rome) before administration of GAHT. The diagnosis of gender dysphoria was confirmed by mental health specialists (Servizio di Adeguamento tra Identità Fisica e Identità Psichica—SAIFIP, AO San Camillo Forlanini), according to the DSM-V criteria. All patients provided their informed consent before the start of the treatment and the study was approved by Ethics Committee “Sapienza” (Ref. 6554; Protocol 1057/2021; 30 November 2021). Each patient underwent physical examination, Ferriman–Gallwey scoring [11], and hormonal blood testing before and after the administration of testosterone enantate 250 mg/mL i.m. once every 28 d. Additional hematological and biochemical data were retrieved from patient’s medical records.

2.1. Hormone Profile

Each recruited subject provided a peripheral blood sample (8 a.m.) after overnight fasting. Serum follicle-stimulating hormone (FSH), luteinizing hormone (LH), and testosterone were quantified by Chemiluminescent Microparticle ImmunoAssay (CMIA, Architect System; Abbott Laboratories, Abbott Park, IL, USA). Detection limits, intra- and inter-assay coefficients of variation, as well as normal ranges from our laboratory were previously described [12].

2.2. DNA Extraction and AR and ER Polymorphisms

DNA was extracted from peripheral blood leukocytes using the Wizard Genomic DNA Purification Kit (Promega, Madison, WI, USA). Extracted DNA was quantified by NanoDrop ND-2000 (Thermo Fisher Scientific, Waltham, MA, USA). For the androgen receptor, located on the X chromosome, the analysis of the percentage of the inactivation of the alleles was performed as proposed by Zitzmann et al. (2004) [13]. Briefly, methylation-sensitive restriction enzyme HpaII (Promega, Madison, WI, USA) was used to digest non-methylated DNA, hence allowing the selective PCR amplification of methylated (inactive) DNA regions. Genomic DNA (1 µg) was digested for 6 h at 37 °C with 20 U of HpaII, according to the manufacturer’s instructions, followed by incubation at 65 °C for 30 min to inactivate the restriction enzyme. Amplification for both digested and non-digested samples was performed for a volume of 25 μL (0.5 ng gDNA, 0.8 μM of each primer, 12.5 μL of Ampli Taq Gold 360 Master Mix—Applied Biosystems, Carlsbad, CA, USA).
The percentage of X chromosome inactivation was then calculated using Gene Mapper Analysis software v. 4.1 (Applied Biosystem, Carlsbad, CA, USA). Variables were defined as: digested allele 1 (A), digested allele 2 (B), non-digested allele 1 (C), and non-digested allele 2 (D). Normalization of the results from the two independent reactions and compensation for iniquitous allele amplifications due to confounding factors other than methylation were achieved through a correction factor [13]:
Allele 1 inactivation: (A/C)/(A/C + B/D)
Allele 2 inactivation: (B/D)/(A/C + B/D)
The extreme values of 0 and 1 correspond to no inactivation and complete inactivation, respectively. A value of 0.5 was interpreted as a random inactivation. The sum of inactivation percentages of alleles 1 and 2 must be equal to 1 (Equation (1) + Equation (2) = 1). To determine the length of the polymorphic fragments in the androgen receptor (AR) and estrogen receptor beta (ERβ) genes (CAG and CA repeats, respectively), DNA samples were analyzed by linear sequencing (3500 Genetic Analyzer, Applied Biosystems) using primers flanking the triplets’ repeat regions. The forward primer for CA and CAG was fluorescently labelled with FAM at 5′ to enable the fragment to be seen during electrophoresis (fragment analysis) [14]. In Supplementary Material Table S1 shows the primers used for the CA and CAG analyses. Raw data from the capillary electrophoresis were analyzed by Gene Mapper Analysis (Applied Biosystems).

2.3. Epigenetic Analysis

One-hundred nanograms of DNA was bisulfite-converted using the EpiTect Plus DNA Bisulfite Kit (Qiagen). PCR amplification of the ESR2 and H19 promoters was performed as previously described by different authors [15,16] Bisulfite-converted DNA served as the template for the polymerase chain reaction (PCR) followed by pyrosequencing (PSQ). The PCR mix included PyroMark PCR Master Mix, 2×, and CoralLoad Concentrate, 10× (Qiagen), 0.2 µM of each primer, 1 µL of converted DNA, and nuclease-free water to a final volume of 25 µL. Primers used for DNA methylation analysis and PCR cycling conditions are shown in Table 1. The pyrosequencing reaction was run on a PyroMark Q96ID (Qiagen), and CpGs methylation analysis was conducted by the PyroMark CpG software (Qiagen). The methylation for each amplicon was calculated as the median of the methylation status of each analyzed CpG.

2.4. Statistical Analysis

Continuous variables are presented as the mean ± SD or the median and IQR, depending on the shape of the distribution curve evaluated by the Kolmogorov–Smirnov test. Categorical variables are presented as counts and percentages. Differences in methylation pattern across samples were calculated by Friedman’s test. A two-tailed p-value of 0.05 was considered significant. Associations between methylation and the other evaluated parameters were evaluated through linear regressions and generalized linear models. All computations were carried out with the Statistical Package for the Social Sciences (SPSS) 25.0 (SPSS Inc., Chicago, IL, USA).

3. Results

A total of 13 AFAB subjects (mean age 29.3 ± 12.6, median 21 (IQR 19.5–43.5) y) provided a blood sample for epigenetic analysis. The baseline characteristics of the caseload and data on AR/ERβ receptors are shown in Table 2 and Table 3, respectively. No significant variations of the investigated serum hormones, with the obvious exception of total testosterone, were detected after GAHT (Figure 1). Figures S1 and S2 show the variations of the other available clinical parameters during the 12 mo after GAHT. All available clinical and biochemical variables were within normal ranges both at baseline (T0) and during follow-up (T6 and T12). Cessation of menses occurred within the second testosterone vial (around two months). As a measure of successful virilization, we monitored the Ferriman–Gallwey score at baseline and during treatment, detecting a significant increase of the scores at T6 and T12 (p < 0.001) (Figure 2).
Post-treatment data from this pilot study showed that the percentage of methylation of ESR2 increased significantly at T6 (adj. p = 0.001) and T12 (adj. p = 0.05) (Figure 3), while no significant difference was found in the H19 promoter methylation (p = 0.237).
Linear regressions were performed to evaluate the associations of the ESR2 methylation variation with available genetic and hormone parameters. Baseline methylation levels were not associated with either AR/ERβ genetic polymorphisms or hormone levels (Table 4). Likewise, no significant associations with ESR2 methylation were found after six months of testosterone treatment (Table 4). On the other hand, total testosterone and the patient’s age resulted in being significantly associated with ESR2 methylation after twelve months of GAHT (Table 4). Moreover, although we could not find an association between ESR2 promoter methylation and the Ferriman–Gallwey score at any time point of the study, the increased virilization of subjects was strongly associated with time from the beginning of treatment (β = 0.918; p < 0.001; adj. R2 = 0.84) and blood testosterone levels (β = 0.774; p < 0.001; adj. R2 = 0.55).
Finally, we used linear models to attempt to evaluate predictors of the difference in methylation after starting testosterone treatment. We could detect that the difference in ESR2 promoter methylation between T6 and baseline was significantly associated with the number of CA repeats of the ERβ receptor, adjusted vs. all considered variables (R2 = 0.62, adj. R2 = 0.35) (Table 5). On the other hand, only age was significantly associated with the ESR2 promoter methylation difference between T12 and baseline (R2 = 0.73, adj. R2 = 0.53) (Table 5).

4. Discussion

Gender dysphoria is a relatively rare condition, where the strict link between hormone receptor signaling and exogenous hormone administration is likely to affect clinical outcomes and phenotypical changes. However, the few genetic and epigenetic studies in this setting have focused only on their possible contribution on the development of gender identity/incongruence. A recent metanalysis from D’Andrea et al. (2020) confirmed the presence of a significantly higher number of CAG repeats of AR in transgender women compared to cisgender controls [5]. Caution must be taken in interpreting this result as a possible etiological genetic contribution, as it must be remarked that it is still unknown to what extent functional variants of hormone receptor genes might be implicated in gender incongruence. Besides, it is likely that hormonal signaling may be modulated by both congenital and acquired epigenetic factors. Currently, the most characterized epigenetic modification is DNA methylation.
Several articles support the effect of methylation on imprinted genes with growth, growth-related hormone concentrations, adiposity, and birth weight [17,18]. Despite the biological plausibility, we did not observe any evidence supporting an association between H19 promoter methylation and response to GAHT.
Hormone receptor (AR, ERα, and ERβ) promoters’ methylation is involved in the development of many diseases, such as atherosclerosis, endometriosis, and cancer [15,19]. DNA methylation is mainly associated with the function of three enzymes belonging to the DNA methyltransferase family (DNMT1, DNMT3A, and DNMT3B), which maintain or induce de novo methylation of a specific DNA sequence [20]. Fernandez et al. (2020) recently demonstrated different methylation patterns of the ESR1 gene promoter in cisgender and transgender subjects [7]. Although the authors identified specific methylation patterns in both assigned male at birth (AMAB) and AFAB people before GAHT, the paucity of data does not allow solid conclusions. In fact, it remains to be ascertained whether epigenetic changes in gender incongruence are causative or simply represent an epiphenomenon of other environmental agents.
From a clinical point of view, the main impact from functional variants in hormone receptor signaling is expected to be found both in the timing and amplitude of the phenotypical changes after GAHT. In physiological puberty, hormonal activation is expected to be modulated by the genetic background and the epigenome [21,22]. Although human studies are lacking, environmental influences in the form of the so-called endocrine-disrupting chemicals (EDCs) are also known to impact physiological puberty, possibly affecting the development of sex hormone-dependent organs (external genitalia and other secondary sex characteristics) through interference on both steroidal and non-steroidal hormone receptors either directly or through epigenetic changes [23,24]. Thus, differences in the timing and acquisition of secondary sexual characteristics may depend on differences in epigenomes induced by GAHT and environmental factors. Therefore, epigenetic changes might ultimately affect the response to treatment with sex steroids, even though this hypothesis has not been thoughtfully investigated in the literature. Sader et al. (2005) were among the first to observe a significant downregulation of AR mRNA expression in blood after both testosterone and estradiol treatment in transgender subjects (six AFAB and six AMAB, respectively) [25]. More recently, Aranda et al. (2017) reported the effect of exogenous sex steroids on the methylation and expression patterns of hormone receptors [8]. The authors observed a downregulation of AR mRNA in the blood of transgender men, as well as an increase in the methylation of AR and ESR1 in transgender women and men, respectively, after 12 mo of treatment. Remarkably, these changes were correlated with several anthropometric (waist circumference), metabolic (HDL, hematocrit), and hormonal (estradiol, testosterone) parameters. This result is particularly important, as it suggests how a fine regulation of the hormone signaling may take place under GAHT. The small sample size of this study does not allow generalizing this observation. Finally, the previously cited study from Fernandez et al. (2020) also analyzed the methylation pattern of ESR1 after 6 mo of hormonal treatment [7]. The authors found that only transgender men after testosterone treatment had a statistically significant increase of methylation. Similarly, our data showed an increase in methylation of ESR2 after testosterone treatment of AFAB subjects within a similar time frame (6 mo). Furthermore, these changes appeared to remain constant up to 12 mo. Therefore, it is possible to hypothesize that treatment with exogenous testosterone can modify the pattern of methylation and, possibly, gene expression. However, our data based on a small study population need to be validated on a larger number of subjects in order to establish the role of epigenetic regulation in the modulation of phenotypical changes after GAHT. Nonetheless, our data suggest that epigenetic changes may be strictly regulated by testosterone. In fact, the modulation of estrogen signaling seems relevant within the first six months of treatment, when methylation might work in tandem with CA polymorphisms, as we detected through linear models between baseline and T6, maintaining itself for more prolonged treatments. Age might be involved as it seems to be constantly and positively associated with methylation of the ESR2 promoter during GATH. In fact, it is known that human ageing is associated with genetic, epigenetic, and environmental factors, which may be, in turn, associated with the development of diseases [26]. Age-dependent modifications in DNA methylations have been already described [27,28], and at least some of these changes might reflect specific environmental exposure of the individual rather than the contribution of genetic factors [29]. Another relevant aspect of this preliminary data stands in the confirmation that estrogen receptor methylation is associated with the administration of exogenous testosterone, while other clinical variables (such as BMI, baseline hormone levels, etc.) seem to play only a marginal role in the induction of epigenetic changes of the ESR2 promoter. If confirmed on a larger caseload and longer follow-up, we may hypothesize exploiting the modulation of androgenic and estrogenic signaling for GAHT tailoring in transgender men and women. Another potential application of these data could be to use the methylation pattern as a marker of androgenic treatment. Being exogenous androgens are the main determinant of this epigenetic marker at T12, it is possible to speculate that this change may also be present following other forms of exogenous testosterone administration, including doping. Epigenetic modifications have been investigated in the setting of sports medicine, hypothesizing an association between epigenome and physical performance traits [30]. Histone modifications (acetylation, methylation) seem to play a role in exercise adaptation and in particular in the expression of musculoskeletal protein expression [31], although their effects on health and physical performance need to be fully elucidated. As our data indicate an influence of exogenous testosterone administration on the epigenetic markers of ESR2, the investigation of these effects in the setting of illegal hormone administration (testosterone, but also other doping agents such as growth hormone, etc.) may be warranted.
In conclusion, data from this pilot study confirm the safety and efficacy of testosterone treatment in AFAB people. ESR2 promoter methylation is increased after testosterone treatment and remains constant up to one year of treatment. The associations we detected with age and testosterone levels suggest that epigenetic modifications might be modulated from these parameters during treatment. Instead, we could not demonstrate a role for H19 promoter methylation, a parentally imprinted gene. Future studies are needed in order to explore the sexually differentiated mechanisms behind DNA methylation considering epigenetic regulation in a larger panel of genes that may directly or indirectly influence hormone treatment-related response.
Nonetheless, it is likely that epigenetically modulated hormone signaling may not only affect the timing of phenotype changes, but also may represent a molecular marker of hormone utilization/abuse in both transgender and cisgender people [32]. Current data will assist in the development of studies with larger caseloads to fully elucidate the role of epigenetic modifications in the pharmacological treatment of transgender people.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biomedicines10020459/s1, Figure S1: Boxplots of the hemocromocytometric parameters before and during testosterone treatment; Figure S2: Boxplots of the glycolipid parameters before and during testosterone treatment; Table S1: Primers used and amplification protocol for fragment analysis of AR and ERβ polymorphisms.

Author Contributions

Conceptualization, F.P. and D.P.; data curation, F.P., G.S., S.D.C., F.C. and V.G.; formal analysis, F.P.; investigation, F.K. and F.R.; resources, F.P., G.S., S.D.C., F.F., F.C., D.P. and V.G.; writing—original draft, F.P. and D.P.; writing—review and editing, C.K., A.L., L.S., V.G. and F.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Italian Ministry of University and Research (Grant Number MIUR-PRIN 2017FC4BS9).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee “Sapienza” (Ref. 6554; Protocol 1057/2021).

Informed Consent Statement

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

Data Availability Statement

The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fisher, A.D.; Senofonte, G.; Cocchetti, C.; Guercio, G.; Lingiardi, V.; Meriggiola, M.C.; Mosconi, M.; Motta, G.; Ristori, J.; Speranza, A.M.; et al. SIGIS–SIAMS–SIE position statement of gender affirming hormonal treatment in transgender and non-binary people. J. Endocrinol. Investig. 2021, 1–17. [Google Scholar] [CrossRef]
  2. Mirabella, M.; Giovanardi, G.; Fortunato, A.; Senofonte, G.; Lombardo, F.; Lingiardi, V.; Speranza, A.M. The Body I Live in. Perceptions and Meanings of Body Dissatisfaction in Young Transgender Adults: A Qualitative Study. J. Clin. Med. 2020, 9, 3733. [Google Scholar] [CrossRef] [PubMed]
  3. Nguyen, H.B.; Chavez, A.M.; Lipner, E.; Hantsoo, L.; Kornfield, S.L.; Davies, R.D.; Epperson, C.N. Gender-Affirming Hormone Use in Transgender Individuals: Impact on Behavioral Health and Cognition. Curr. Psychiatry Rep. 2018, 20, 110. [Google Scholar] [CrossRef] [PubMed]
  4. Hembree, W.C.; Cohen-Kettenis, P.T.; Gooren, L.; Hannema, S.E.; Meyer, W.J.; Murad, M.H.; Rosenthal, S.M.; Safer, J.D.; Tangpricha, V.; T’Sjoen, G.G. Endocrine Treatment of Gender-Dysphoric/Gender-Incongruent Persons: An Endocrine Society* Clinical Practice Guideline. J. Clin. Endocrinol. Metab. 2017, 102, 3869–3903. [Google Scholar] [CrossRef]
  5. D’Andrea, S.; Pallotti, F.; Senofonte, G.; Castellini, C.; Paoli, D.; Lombardo, F.; Lenzi, A.; Francavilla, S.; Francavilla, F.; Barbonetti, A. Polymorphic Cytosine-Adenine-Guanine Repeat Length of Androgen Receptor Gene and Gender Incongruence in Trans Women: A Systematic Review and Meta-Analysis of Case-Control Studies. J. Sex. Med. 2020, 17, 543–550. [Google Scholar] [CrossRef] [PubMed]
  6. Almstrup, K.; Frederiksen, H.; Andersson, A.-M.; Juul, A. Levels of endocrine-disrupting chemicals are associated with changes in the peri-pubertal epigenome. Endocr. Connect. 2020, 9, 845–857. [Google Scholar] [CrossRef]
  7. Fernández, R.; Ramírez, K.; Gómez-Gil, E.; Cortés-Cortés, J.; Mora, M.; Aranda, G.; Zayas, E.D.; Esteva, I.; Almaraz, M.C.; Guillamon, A.; et al. Gender-Affirming Hormone Therapy Modifies the CpG Methylation Pattern of the ESR1 Gene Promoter After Six Months of Treatment in Transmen. J. Sex. Med. 2020, 17, 1795–1806. [Google Scholar] [CrossRef]
  8. Aranda, G.; Fernandez-Rebollo, E.; Pradas-Juni, M.; Hanzu, F.A.; Kalko, S.; Halperin, I.; Mora, M. Effects of sex steroids on the pattern of methylation and expression of the promoter region of estrogen and androgen receptors in people with gender dysphoria under cross-sex hormone treatment. J. Steroid Biochem. Mol. Biol. 2017, 172, 20–28. [Google Scholar] [CrossRef]
  9. Murphy, R.; Ibáñez, L.; Hattersley, A.; Tost, J. IGF2/H19 hypomethylation in a patient with very low birthweight, preocious pubarche and insulin resistance. BMC Med Genet. 2012, 13, 42. [Google Scholar] [CrossRef] [Green Version]
  10. Wu, Y.; Peterson, K.E.; Sánchez, B.N.; Dolinoy, D.C.; Mercado-Garcia, A.; Rojo, M.T.; Goodrich, J.M. Association of blood leukocyte DNA methylation at LINE-1 and growth-related candidate genes with pubertal onset and progression. Epigenetics 2018, 13, 1222–1233. [Google Scholar] [CrossRef] [Green Version]
  11. Yildiz, B.O.; Bolour, S.; Woods, K.; Moore, A.; Azziz, R. Visually scoring hirsutism. Hum. Reprod. Updat. 2009, 16, 51–64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Pallotti, F.; Senofonte, G.; Pelloni, M.; Cargnelutti, F.; Carlini, T.; Radicioni, A.F.; Rossi, A.; Lenzi, A.; Paoli, D.; Lombardo, F. Androgenetic alopecia: Effects of oral finasteride on hormone profile, reproduction and sexual function. Endocrinology 2020, 68, 688–694. [Google Scholar] [CrossRef] [PubMed]
  13. Zitzmann, M.; Depenbusch, M.; Gromoll, J.; Nieschlag, E. X-Chromosome Inactivation Patterns and Androgen Receptor Functionality Influence Phenotype and Social Characteristics as Well as Pharmacogenetics of Testosterone Therapy in Klinefelter Patients. J. Clin. Endocrinol. Metab. 2004, 89, 6208–6217. [Google Scholar] [CrossRef]
  14. Grassetti, D.; Giannandrea, F.; Paoli, D.; Masciandaro, P.; Figura, V.; Carlini, T.; Rizzo, F.; Lombardo, F.; Lenzi, A.; Gandini, L. Androgen receptor polymorphisms and testicular cancer risk. Andrology 2015, 3, 27–33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Spitzwieser, M.; Entfellner, E.; Werner, B.; Pulverer, W.; Pfeiler, G.; Hacker, S.; Cichna-Markl, M. Hypermethylation of CDKN2A exon 2 in tumor, tumor-adjacent and tumor-distant tissues from breast cancer patients. BMC Cancer 2017, 17, 260. [Google Scholar] [CrossRef] [Green Version]
  16. Montrose, L.; Padmanabhan, V.; Goodrich, J.M.; Domino, S.E.; Treadwell, M.C.; Meeker, J.D.; Watkins, D.J.; Dolinoy, D.C. Maternal levels of endocrine disrupting chemicals in the first trimester of pregnancy are associated with infant cord blood DNA methylation. Epigenetics 2018, 13, 301–309. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Huang, R.-C.; Galati, J.C.; Burrows, S.; Beilin, L.J.; Li, X.; Pennell, C.E.; Van Eekelen, J.; Mori, T.A.; Adams, L.A.; Craig, J.M. DNA methylation of the IGF2/H19 imprinting control region and adiposity distribution in young adults. Clin. Epigenetics 2012, 4, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Deodati, A.; Inzaghi, E.; Liguori, A.; Puglianiello, A.; Germani, D.; Brufani, C.; Fintini, D.; Cappa, M.; Barbetti, F.; Cianfarani, S. IGF2Methylation Is Associated with Lipid Profile in Obese Children. Horm. Res. Paediatr. 2013, 79, 361–367. [Google Scholar] [CrossRef]
  19. Holtzman, L.; Gersbach, C.A. Editing the Epigenome: Reshaping the Genomic Landscape. Annu. Rev. Genom. Hum. Genet. 2018, 19, 43–71. [Google Scholar] [CrossRef]
  20. Chen, H.; Malentacchi, F.; Fambrini, M.; Harrath, A.H.; Huang, H.; Petraglia, F. Epigenetics of Estrogen and Progesterone Receptors in Endometriosis. Reprod. Sci. 2020, 27, 1967–1974. [Google Scholar] [CrossRef]
  21. Lomniczi, A.; Wright, H.; Ojeda, S.R. Epigenetic regulation of female puberty. Front. Neuroendocr. 2015, 36, 90–107. [Google Scholar] [CrossRef] [PubMed]
  22. Shalev, D.; Melamed, P. The role of the hypothalamus and pituitary epigenomes in central activation of the reproductive axis at puberty. Mol. Cell. Endocrinol. 2020, 518, 111031. [Google Scholar] [CrossRef] [PubMed]
  23. Cargnelutti, F.; Di Nisio, A.; Pallotti, F.; Sabovic, I.; Spaziani, M.; Tarsitano, M.G.; Paoli, D.; Foresta, C. Effects of endocrine disruptors on fetal testis development, male puberty, and transition age. Endocrinology 2021, 72, 358–374. [Google Scholar] [CrossRef] [PubMed]
  24. Lopez-Rodriguez, D.; Franssen, D.; Bakker, J.; Lomniczi, A.; Parent, A.-S. Cellular and molecular features of EDC exposure: Consequences for the GnRH network. Nat. Rev. Endocrinol. 2021, 17, 83–96. [Google Scholar] [CrossRef] [PubMed]
  25. Sader, M.A.; McGrath, K.; Hill, M.D.; Bradstock, K.F.; Jimenez, M.; Handelsman, D.J.; Celermajer, D.S.; Death, A.K. Androgen receptor gene expression in leucocytes is hormonally regulated: Implications for gender differences in disease pathogenesis. Clin. Endocrinol. 2005, 62, 56–63. [Google Scholar] [CrossRef]
  26. Ben-Avraham, D.; Muzumdar, R.H.; Atzmon, G. Epigenetic genome-wide association methylation in aging and longevity. Epigenomics 2012, 4, 503–509. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Bocklandt, S.; Lin, W.; Sehl, M.E.; Sánchez, F.J.; Sinsheimer, J.S.; Horvath, S.; Vilain, E. Epigenetic Predictor of Age. PLoS ONE 2011, 6, e14821. [Google Scholar] [CrossRef] [Green Version]
  28. Florath, I.; Butterbach, K.; Müller, H.; Bewerunge-Hudler, M.; Brenner, H. Cross-sectional and longitudinal changes in DNA methylation with age: An epigenome-wide analysis revealing over 60 novel age-associated CpG sites. Hum. Mol. Genet. 2014, 23, 1186–1201. [Google Scholar] [CrossRef]
  29. Tan, Q.; Heijmans, B.T.; Hjelmborg, J.V.B.; Soerensen, M.; Christensen, K.; Christiansen, L. Epigenetic drift in the aging genome: A ten-year follow-up in an elderly twin cohort. Int. J. Epidemiol. 2016, 45, 1146–1158. [Google Scholar] [CrossRef] [Green Version]
  30. Ehlert, T.; Simon, P.; Moser, D.A. Epigenetics in Sports. Sports Med. 2013, 43, 93–110. [Google Scholar] [CrossRef]
  31. McGee, S.L.; Fairlie, E.; Garnham, A.P.; Hargreaves, M. Exercise-induced histone modifications in human skeletal muscle. J. Physiol. 2009, 587, 5951–5958. [Google Scholar] [CrossRef]
  32. Schwarzenbach, H. Impact of Physical Activity and Doping on Epigenetic Gene Regulation. Drug Test. Anal. 2011, 3, 682–687. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Hormone profile of the recruited subjects at baseline and after 6 mo (T6) and 12 mo (T12) of GAHT. Friedman test. p-values are Bonferroni adjusted.
Figure 1. Hormone profile of the recruited subjects at baseline and after 6 mo (T6) and 12 mo (T12) of GAHT. Friedman test. p-values are Bonferroni adjusted.
Biomedicines 10 00459 g001
Figure 2. Ferriman–Gallwey score at baseline and after gender assignment hormone therapy (T6 and T12). Friedman’s test. p-values are Bonferroni adjusted.
Figure 2. Ferriman–Gallwey score at baseline and after gender assignment hormone therapy (T6 and T12). Friedman’s test. p-values are Bonferroni adjusted.
Biomedicines 10 00459 g002
Figure 3. ESR2 methylation percentage at baseline and after gender assignment hormone therapy (T6 and T12). Friedman’s test. p-values are Bonferroni adjusted.
Figure 3. ESR2 methylation percentage at baseline and after gender assignment hormone therapy (T6 and T12). Friedman’s test. p-values are Bonferroni adjusted.
Biomedicines 10 00459 g003
Table 1. DNA methylation primer sequences and PCR cycling conditions.
Table 1. DNA methylation primer sequences and PCR cycling conditions.
Primer/Sequence to AnalyzeESR2H19
Forward PCR Primer5′-GGAGGTTGAGAGAAATAATTGTTTTTTGA- 3′5′-TTTGTTGATTTTATTAAGGGAG-3′
Reverse PCR Primer5′-[Biotin]-ATAAACACACCCACCTTACCTTCTCTA-3′5′-[Biotin]-CTATAAATAAACCCCAACCAAAC-3′
Sequencing Primer5′-GAAATAATTGTTTTTTGAAATTTG-3′5′-GTGTGGAATTAGAAGT-3′
Sequence to AnalyzeTAGGGYGAAGAGTAGGYGGYGAGYGTTGGGTYGGGGAGGGATTATTYGAGTTGYGAYGGGTTTTGGGGTTGYGGGGTAGGTYGYGYGGYGGTAGTGTAGGTTTATATATTATAGTT
Annealing Temperature (°C)61 °C for 30 s51 °C for 30 s
Table 2. Baseline characteristics of the enrolled AMAB subjects. Continuous data are shown as the mean ± standard deviations and the median (in brackets).
Table 2. Baseline characteristics of the enrolled AMAB subjects. Continuous data are shown as the mean ± standard deviations and the median (in brackets).
Age (Years)29.3 ± 12.6 (21)
BMI26.6 ± 1.4 (27)
CAG repeats expression (AR)23.1 ± 1.9 (23.5)
CA repeats expression (ERβ)22.1 ± 1.4 (21.5)
Ferriman–Gallwey score6.4 ± 1.3 (7.0)
Creatinine (mg/dL)0.6 ± 0.1 (0.6)
Red blood cells (×106/mL)4.6 ± 0.5 (4.6)
Hemoglobin (g/dL)12.7 ± 1.2 (12.6)
Hematocrit (%)39.0 ± 3.1 (39.9)
White blood cells (×103/mL)7.1 ± 1.7 (7.1)
Platelets (×103/mL)243.8 ± 80.3 (252.0)
Glycemia (mg/dL)85.8 ± 11.8 (85.0)
HbA1c (%)5.1 ± 0.4 (5.1)
AST (UI/L)23.5 ± 19.8 (19)
ALT (UI/L)22.1 ± 10.8 (18.0)
γGT (UI/L)19.5 ± 12.9 (15.0)
Total cholesterol (mg/dL)153.8 ± 27.2 (150.0)
HDL (mg/dL)56.3 ± 11.6 (58.0)
LDL (mg/dL)80.5 ± 23.1 (85.0)
Triglycerides (mg/dL)83.8 ± 38.3 (71.0)
FSH (mUI/mL)10.5 ± 13.5 (5.4)
LH (mUI/mL)7.8 ± 4.6 (6.1)
Prolactin (ng/dL)13.2 ± 4.7 (12.2)
17β estradiol (pg/mL)84.0 ± 86.3 (56.1)
Total testosterone (nmol/L)1.2 ± 0.5 (1.3)
Table 3. Polymorphisms of the androgen receptor (CAG repeats) and estrogen receptor beta (CA repeats). For androgen receptor, percentage of X inactivation is provided.
Table 3. Polymorphisms of the androgen receptor (CAG repeats) and estrogen receptor beta (CA repeats). For androgen receptor, percentage of X inactivation is provided.
PatientCAG Repeats (Allele 1)Percentage of Inactivation (Allele 1)CAG Repeats (Allele 2)Percentage of Inactivation (Allele 2)CA Repeats (Allele 1)CA Repeats (Allele 2)
#12441.02559.02424
#22362.12537.92024
#31929.62670.41823
#41750.82449.32020
#52035.82464.22023
#62051.32148.72323
#72456.72643.31825
#82156.62443.42124
#92551.42648.62121
#102258.32541.71726
#1120/20/2424
#122445.42654.62424
#132435.82664.22121
Table 4. Coefficients from regression models predicting baseline ESR2 methylation.
Table 4. Coefficients from regression models predicting baseline ESR2 methylation.
BCI 95%Betap
Coefficients from regression models predicting baseline ESR2 methylation.
CA repeats−0.529−1.302–0.243−0.4140.533
CAG repeats−0.175−0.762–0.413−0.1930.527
Total testosterone1.143−0.852–3.1390.3550.233
Estradiol−0.007−0.019–0.006−0.3340.164
Age0.026−0.103–0.1550.1850.653
Coefficients from regression models predicting T6 ESR2 methylation.
CA repeats0.437−0.428–1.3020.3180.290
CAG repeats−0.111−0.751–0.529−0.1140.711
Total testosterone0.030−0.296–0.3570.0620.841
Estradiol0.009−0.060–0.0770.0830.788
Age−0.034−0.130–0.062−0.2300.450
Coefficients from regression models predicting T12 ESR2 methylation.
CA repeats−0.110−1.062–0.841−0.0770.803
CAG repeats0.110−0.560–0.7800.1080.725
Total testosterone0.2800.086–0.4740.6910.009
Estradiol−0.009−0.028–0.011−0.2760.361
Age0.1060.031–0.1820.6820.010
Table 5. Linear regression model coefficients predicting the increase in ESR2 promoter methylation between baseline and T6.
Table 5. Linear regression model coefficients predicting the increase in ESR2 promoter methylation between baseline and T6.
BCI 95%BetapPartial Eta Squared
Linear regression model coefficients predicting the increase in ESR2 promoter methylation between baseline and T6.
CA repeats1.4000.357–2.4440.8220.0160.590
CAG repeats−0.481−1.385–0.423−0.4000.2490.185
Total testosterone0.445−0.061–0.8500.7280.0760.382
Estradiol0.079−0.016–0.0730.6090.0900.356
Age−0.067−0.183–0.050−0.3610.2200.206
Linear regression model coefficients predicting the increase in ESR2 promoter methylation between baseline and T12.
BCI 95%BetapPartial Eta Squared
CA repeats0.498−0.270–1.2660.3040.1690.251
CAG repeats0.285−0.421–0.9910.2460.3710.115
Total testosterone0.169−0.179–0.5180.3670.2890.159
Estradiol0.014−0.007–0.0360.4070.1520.270
Age0.1090.010–0.2180.6150.0480.449
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pallotti, F.; Senofonte, G.; Konstantinidou, F.; Di Chiano, S.; Faja, F.; Rizzo, F.; Cargnelutti, F.; Krausz, C.; Paoli, D.; Lenzi, A.; et al. Epigenetic Effects of Gender-Affirming Hormone Treatment: A Pilot Study of the ESR2 Promoter’s Methylation in AFAB People. Biomedicines 2022, 10, 459. https://doi.org/10.3390/biomedicines10020459

AMA Style

Pallotti F, Senofonte G, Konstantinidou F, Di Chiano S, Faja F, Rizzo F, Cargnelutti F, Krausz C, Paoli D, Lenzi A, et al. Epigenetic Effects of Gender-Affirming Hormone Treatment: A Pilot Study of the ESR2 Promoter’s Methylation in AFAB People. Biomedicines. 2022; 10(2):459. https://doi.org/10.3390/biomedicines10020459

Chicago/Turabian Style

Pallotti, Francesco, Giulia Senofonte, Fani Konstantinidou, Silvia Di Chiano, Fabiana Faja, Flavio Rizzo, Francesco Cargnelutti, Csilla Krausz, Donatella Paoli, Andrea Lenzi, and et al. 2022. "Epigenetic Effects of Gender-Affirming Hormone Treatment: A Pilot Study of the ESR2 Promoter’s Methylation in AFAB People" Biomedicines 10, no. 2: 459. https://doi.org/10.3390/biomedicines10020459

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