Next Article in Journal
Retinal Alterations Predict Early Prodromal Signs of Neurodegenerative Disease
Previous Article in Journal
Bioactive Hydrogel Based on Collagen and Hyaluronic Acid Enriched with Freeze-Dried Sheep Placenta for Wound Healing Support
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Placenta-Specific Transcripts Containing Androgen Response Elements Are Altered In Silico by Male Growth Outcomes

1
Early Origins of Adult Health Research Group, Health and Biomedical Innovation, UniSA: Clinical and Health Sciences, University of South Australia, Adelaide, SA 5000, Australia
2
Mater Medical Research Institute, The University of Queensland, Brisbane, QLD 4101, Australia
3
Flinders Health and Medical Research Institute, College of Medicine and Public Health, Flinders University, Adelaide, SA 5042, Australia
4
Department of Obstetrics, Gynaecology and Newborn Health, Mercy Hospital for Women, The University of Melbourne, Heidelberg, VIC 3084, Australia
5
Department of Obstetrics and Gynaecology, Monash University, Melbourne, VIC 3168, Australia
6
The Ritchie Centre, Hudson Institute of Medical Research, Melbourne, VIC 3168, Australia
7
School of Pharmacy and Medical Sciences, Griffith University, Gold Coast, QLD 4222, Australia
8
School of Health, University of the Sunshine Coast, Sunshine Coast, QLD 4556, Australia
9
Women and Babies Research, Faculty of Medicine and Health, The University of Sydney, Camperdown, NSW 2050, Australia
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(3), 1688; https://doi.org/10.3390/ijms25031688
Submission received: 27 November 2023 / Revised: 18 January 2024 / Accepted: 23 January 2024 / Published: 30 January 2024
(This article belongs to the Section Molecular Pathology, Diagnostics, and Therapeutics)

Abstract

:
A birthweight centile (BWC) below the 25th is associated with an elevated risk of adverse perinatal outcomes, particularly among males. This male vulnerability may stem from alterations in placenta-specific androgen signalling, a signalling axis that involves the androgen receptor (AR)-mediated regulation of target genes containing androgen response elements (AREs). In this study, we examined global and ARE-specific transcriptomic signatures in term male placentae (≥37 weeks of gestation) across BWC subcategories (<10th, 10th–30th, >30th) using RNA-seq and gene set enrichment analysis. ARE-containing transcripts in placentae with BWCs below the 10th percentile were upregulated compared to those in the 10th–30th and >30th percentiles, which coincided with the enrichment of gene sets related to hypoxia and the suppression of gene sets associated with mitochondrial function. In the absence of ARE-containing transcripts in silico, <10th and 10th–30th BWC subcategory placentae upregulated gene sets involved in vasculature development, immune function, and cell adhesion when compared to those in the >30th BWC subcategory. Collectively, our in silico findings suggest that changes in the expression of ARE-containing transcripts in male placentae may contribute to impaired placental vasculature and therefore result in reduced fetal growth outcomes.

1. Introduction

Disrupted fetal growth can result in a small for gestational age (SGA; <10th birthweight centile (BWC)) neonate who is at increased risk of morbidity and mortality [1], but the contributing molecular mechanisms remain unclear [2,3,4,5]. While SGA is a known risk factor for adverse neonatal outcomes, recent work identified that the probability of adverse outcomes was significantly increased up to the 25th BWC when compared with neonates born between the 25th and 90th BWC, and that male neonates were more likely to have a severe adverse outcome compared to females [6]. These observations favour previous studies describing the male disadvantage hypothesis [7] and emphasise that fetal growth to less than the 25th BWC contributes to adverse outcomes. However, studies characterising molecular pathways associated with reduced fetal growth primarily focus comparisons between the <10th and >10th BWC groups. Based on initial observations in the placenta [7], Clifton and colleagues hypothesised that males prioritise growth over adaptation in utero in the presence of a stressor. However, this male-specific response impairs their ability to adequately respond to shifts in the intrauterine environment such as hypoxic insults and inflammation, which then places them at a greater risk for adverse perinatal outcomes. Despite this, the exact mechanisms that contribute to this male-specific response is unclear. Previous work from our group indicates that male fetal growth outcomes may be mediated through differences in placental androgen signalling regulated by the differential expression of androgen receptor (AR) variants [8].
Multiple AR protein variants have been characterised in human [8] and sheep [9] placenta that challenge the current understanding of the placenta-specific androgen signalling axis. This work showed that placental AR protein variants may influence the regulation of key AR-mediated transcripts, particularly in males, associated with altered growth outcomes. Specifically, we showed that increased nuclear expression of the N-terminally truncated variant, AR-45, was associated with continued male growth and increased angiogenic downstream gene expression [8]. However, in males with growth <10th BWC, this AR-45-specific adaptation was impaired. Indeed, different AR protein variants can initiate transcription through an interaction with the cis-regulatory element (androgen response element (ARE)) associated with the AR [10]. These 15 base-pair sequences are composed of two palindromic half sites separated by a three-base nucleotide spacer (i.e., AGAACANNNTGTTCT). ARE half and full sites are important initial factors that drive the transcriptional regulation of genes belonging to the AR-mediated transcriptome. Despite this growing understanding of the AR-mediated transcriptome, its function in the placenta is largely uncharacterised. Thus, in the current study, we characterised global and ARE-specific transcriptomic signatures in male placentae from different BWCs (<10th, 10th–30th, and >30th) and performed gene set enrichment analysis (GSEA) on global and ARE-containing transcripts.

2. Results

2.1. Maternal and Neonatal Characteristics

No maternal measures were significantly different between the BWC subgroups; however, gestation, BWC, birthweight and placenta weight were significantly different (Table 1).

2.2. Placental Transcript Expression Variation in Response to Fetal Growth Outcomes

There were 26 transcripts differentially expressed between >30th and <10th BWC males, and 34 transcripts were differentially expressed between 10th–30th and <10th BWC males (Figure 1A). Three transcripts (PNCK, MXI1, and ARNT2) were commonly differentially expressed in the >30th vs. <10th BWC and 10th–30th vs. <10th BWC comparisons (Figure 1B). XIST was commonly differentially expressed in the 10th–30th vs. <10th BWC and >30th vs. 10th–30th BWC comparisons. Twenty-three transcripts were unique to the >30th vs. <10th BWC comparison (1 upregulated, 22 downregulated), and 30 were unique to the 10th–30th vs. <10th BWC comparison (2 upregulated, 28 downregulated; Figure 1B). A complete list of differentially expressed transcripts between the BWC subcategories can be found in Table 2.

2.3. Differentially Expressed Transcripts Contain ARE Full and Half Sites

Using a comprehensive and refined ARE dataset previously published [10], differentially expressed transcripts from the current study were integrated to identify whether they contained ARE full (AGAACANNNTGTTCT) and/or half sites, and whether these transcripts were previously reported to be positively or negatively regulated by these ARE sites. Of the 34 transcripts that were differentially expressed in the <10th vs. 10th–30th BWC comparison, 16 contained ARE full sites and 17 contained ARE half sites, two of which have previously been reported to be regulated by both ARE full and half sites (HK2 and SLCO5A1; Table 2). In the <10th vs. >30th BWC comparison, 17 transcripts contained ARE full sites and 21 contained ARE half sites (Table 2); however, none of these have been previously reported to be regulated by either ARE full or half sites.

2.4. Gene Set Enrichment Analysis

We performed GSEA on global gene lists, as well as gene lists of transcripts with or without ARE full sites (ARE-positive and ARE-negative, respectively). For global GSEA (see Supplementary Table S1), terms upregulated in <10th BWC relative to >30th BWC were primarily involved in hypoxia and vasculature development. Comparatively, terms downregulated in <10th BWC relative to >30th BWC were primarily involved in mitochondrial respiration. Consistent with the <10th BWC relative to >30th BWC comparison, terms involved in hypoxia were upregulated in the <10th BWC relative to the 10th–30th BWC comparison. In >30th BWC placentae, terms involved in protein localisation to phagophore assembly sites and COPII-coated vesicle cargo loading were upregulated when compared with 10th–30th BWC placentae. No terms were downregulated in male >30th BWC vs. 10th–30th BWC.
Given our primary research question focused on understanding placental molecular signatures in males that may be regulated by androgen signalling, we performed subsequent GSEA on ARE-positive and ARE-negative gene lists and compared any findings to male-specific global GSEA terms (Supplementary Table S1). Interestingly, in the ARE-negative GSEA, there were 106 unique terms upregulated in the <10th BWC vs. >30th BWC comparison; these terms were primarily involved in cell adhesion and vasculature development (Figure 2A). Similarly, in the ARE-negative GSEA, 460 unique terms were upregulated in the 10th–30th BWC vs. >30th BWC comparison; these terms were primarily involved in immune function and vasculature development (Figure 2C). Twenty-two terms were commonly upregulated in smaller male placentae (<10th and 10th–30th BWC) when compared with >30th BWC placentae, including myogenesis, extracellular matrix structure, and integrin-mediated signalling.

3. Discussion

Despite growing evidence that adverse outcomes for the developing male fetus are more likely to occur up to the 25th BWC when compared with those born within the 25th to 90th BWC [6], the underlying molecular mechanisms driving these outcomes remain poorly understood. Herein, we show for the first time that male placenta-specific transcriptomic signatures are altered by distinct fetal growth outcomes. Importantly, GSEA identified that in the absence of transcripts containing ARE full sites, there was increased enrichment in gene sets involved in cell adhesion, vasculature development, and immune function in the <10th and 10th–30th BWC, when compared with the >30th BWC. In line with previous work by our team and others [8,11,12,13,14,15,16], these in silico findings suggest that the increased risk of poor outcomes associated with reduced male fetal growth may be caused, in part, by AR-mediated placental dysfunction that results in perturbations to vasculature development, cell adhesion, and immune function.
Although several of the differentially expressed transcripts in the <10th BWC males contained AREs, only HK2 and SLCO5A1 have previously been reported to be positively or negatively regulated by ARE full sites, respectively. The observed increase in HK2 expression in the <10th when compared with >30th BWC in male placentae is indicative of increased activity of placental androgen signalling. HK2 encodes the hexokinase 2 protein, is involved in the metabolism of glucose [17], and has also been reported to promote the proliferation of certain cancer cells [18,19,20]. The overexpression of AR results in increased HK2 mRNA and protein expression [21]. HK2 expression can also be activated by HIF1α [22,23], which suggests that increased HK2 expression in male placentae <10th BWC is the result of either androgen- or hypoxia-mediated transactivation. Indeed, there is growing evidence that androgen signalling can result in placental hypoxia via a reduction in uterine blood flow [24,25]. Likewise, hypoxia can increase androgen concentrations and enhance AR activity in endocrine tissue, including the prostate [26,27,28,29,30,31]. Thus, the observed increase in HK2 expression in male placentae <10th BWC may be directly regulated by either androgen- or hypoxia-mediated signalling. It is also possible that the enhanced activation of androgen-mediated signalling via a hypoxic in utero environment drives HK2 expression in male placentae <10th BWC. Evidently, further studies are needed to untangle androgen- and hypoxia-dependent signalling pathways in the placenta and define how they contribute to male fetal growth outcomes.
SLCO5A1 expression was also increased in male placentae <10th BWC when compared with male >30th placentae. Whereas HK2 is positively regulated by ARE full sites, SLCO5A1 is reported to be negatively regulated by ARE full sites [10]. SLCO5A1 encodes the solute carrier organic anion transporter family member 5A1 protein. SLCO5A1 not only functions to transport amphipathic molecules across the plasma membrane, but may also function to mediate pathways involved in differentiation and migration [32]. In normal and cancerous prostate tissue, SLCO5A1 expression is significantly reduced in response to androgen deprivation treatment [33], which may be the result of reactivated AR signalling following such treatment. Like HK2, SLCO5A1 expression can be altered by hypoxic insults [34]. Indeed, given that previously published studies report the negative regulation of SLCO5A1 via AR activation, it is likely that alternative pathways such as those regulated by hypoxia are involved in its increased expression in the placentae of small males. Collectively, the observed increase in HK2 and SLCO5A1 expression in male placentae <10th BWC may indicate their direct regulation by either androgen- or hypoxia-mediated signalling. It is important to note that the interpretation of these data is based on comprehensive studies of the AR cistrome in prostate cancer. Given that tissue-specific responses to androgens are reported [35], future studies would benefit from defining the placenta-specific AR cistrome; however, this is beyond the scope of the current study.
Our global GSEA revealed term enrichments consistent with other studies demonstrating that impaired oxygen supply to the fetoplacental unit can result in fetoplacental hypoxia and reduced fetal growth [36]. In <10th BWC males compared with larger males, we identified that upregulated terms were involved in the response to hypoxia, whereas terms that were significantly downregulated were involved in mitochondrial respiration. This finding is supported by previous studies that demonstrate that male-specific alterations to mitochondrial respiration in response to prenatal hypoxia may impair mitochondrial function [37,38]. Our current data suggest that males are less adaptable to intrauterine environments that contribute to reduced fetal growth, as evidenced by disrupted mitochondrial function. This impairment may contribute to an increase in the prevalence of intrauterine morbidity and mortality associated with reduced fetal growth in males [39]. Although androgens are increasingly being recognised as important regulators of mitochondrial function [40], recent work demonstrated that androgen exposure in the trophoblast cell line ACH-3P resulted in mitochondrial dysfunction and enhanced reactive oxygen species (ROS) production [41]. Despite this recent work demonstrating androgen-mediated mitochondrial dysfunction in trophoblast cells, it cannot be concluded that the observed changes to mitochondrial processes in the current study are due to altered androgen signalling in utero. Nonetheless, aberrant ROS production can result in oxidative stress and subsequent placental insufficiency [42], and thus, reduced fetal growth outcomes. Therefore, it is to be expected that an enrichment of terms associated with these processes was reported in <10th BWC males. Indeed, comprehensive studies in the placenta are needed to understand whether androgen-dependent mechanisms contribute to the regulation of mitochondrial function in utero and determine what this means for sex-specific fetal growth outcomes.
Intriguingly, our in silico analysis demonstrated that the absence of ARE-containing transcript results in an enrichment of gene sets involved in vasculature development, cell adhesion, and immune function in small-male placentae (i.e., <10th and 10th–30th BWC) when compared with the placentae males with normal growth (>30th BWC). These findings support growing evidence that functional AR signalling contributes to placental vasculature impairment, and we postulate that this may be due to changes in the expression of AR protein variants. The placentae of pregnancies associated with elevated androgen concentrations during pregnancy have reduced volume and weight and enhanced placental inflammatory pathways and lesions [15,43] when compared with the placentae of healthy controls. In support of these observations, animal studies have demonstrated that elevated androgen concentrations reduce uterine artery blood flow, spiral artery elongation, and placental oxygenation, which, in turn, results in reduced fetal growth [12,24,44]. Our seminal findings of multiple AR protein variants in the human placenta demonstrated that increased androgen concentrations in vitro reduced the nuclear expression of AR-45, whereas the nuclear expression of other measured AR protein variants remained unchanged [8]. AR-45 functions antagonistically to the full-length AR protein variant, is associated with increased angiogenic downstream target gene expression, and may therefore modulate the AR-mediated transcriptomic signature to enhance placental vascular function [45,46]. Thus, impaired placental AR-45 function may enable the dysregulated activity of placenta-specific androgen signalling and exacerbate placental dysfunction and fetal growth perturbations in the presence of chronic intrauterine stressors. Although we did not measure AR-45 protein expression or circulating maternal androgen concentrations in the current cohort due to sample availability, future studies would benefit from examining the associations between androgen ligand concentrations (e.g., testosterone or dihydrotestosterone), AR protein variant expression, and ARE-containing transcript abundance, to identify how any observed interactions vary amongst fetal growth outcomes.
Despite our work highlighting distinct changes to placenta-specific transcriptomic signatures in males from different BWCs, and that changes to gene set enrichment were associated with the presence or absence of ARE-containing transcripts, this study is not without limitations. Indeed, the decision to analyse male placentae only captures the ‘male’ side of an important story of sex differences in placentology. However, the decision to analyse male samples only was hypothesis-driven and based on our previously published works that show that male, but not female, placentae modulate androgen signalling pathways in response to a maternal stressor or reduced fetal growth outcomes [8], or ex vivo inflammatory challenges [11]. Indeed, these studies support our seminal findings that female and male placentae respond differently to the same intrauterine environment [7], and we expect that the molecular adaptations that occur in females across different BWCs would be distinct to those observed in males. Nonetheless, future studies would benefit from interrogating the transcriptomic profiles and molecular regulatory systems that occur in the placentae of small females; however, this is beyond the scope of our current study. Another important consideration is that we used publicly available prostate cancer data [10] to inform us of transcripts that may be regulated by the AR in the placenta, as determined by the presence of an ARE full site. Despite this, while AREs are global motifs, the direction of regulation may differ between tissues. Indeed, similarities between cancer and placenta biology exist [47]; it may therefore be inferred that comparable responses to androgens, as measured by the expression of ARE containing transcripts, between the two also exist. However, given that differences in AR protein variant profiles in placenta and prostate cancer have been observed [8], it is likely that tissue-specific responses to androgens, and thus, changes in the direction of the regulation of AR-mediated transcripts occur. Future studies that apply AR ChIP-seq coupled with RNA-seq studies are therefore needed to define the placenta-specific AR-mediated pathways altered in males of different BWCs; however, this is beyond the scope of the current study. Regardless, our in silico work does provide the rationale to facilitate such future studies.

4. Materials and Methods

4.1. Sample Cohort

Ethical approval was provided from the Mater Human Research Ethics Committee (53157). Placentae were derived from five biobank sites across Australia based in Melbourne, Brisbane, and Sydney. Each biobank provided term placentae from all BWC categories, and samples were collected and processed following the CoLab protocol as previously described [48]. Samples were matched for maternal BMI and age. For this study, only male placentae with a <10th BWC outcome (n = 16) or >10th BWC outcome (n = 21) were used. Placentae from >10th BWC outcomes were stratified into two subcategories: 10th–30th BWC (n = 12) and >30th BWC (n = 9). The decision to use 10th–30th BWC samples was informed by unpublished gene array data from our group.

4.2. RNA Extraction and RNA Sequencing

RNA extraction was completed using a QIAGEN RNeasy Mini Kit (QIAGEN, Hilden, Germany) as per the manufacturer’s protocol. Total RNA concentration was recorded using NanoDrop 100 spectrophotometer. The RNA integrity numbers (RIN) for samples used in the initial RNA-Seq analysis were measured using an Agilent RNA 6000 Nano Kit (Agilent, Santa Clara, CA, USA), as per the manufacturer’s protocol. RNA-Seq was performed on an Illumina NovaSeq 6000 (Illumina, San Diego, CA, USA),100 base-pair single-end reads, at the Australian Genomic Research Facility (AGRF), Brisbane.

4.3. Data Processing and Differential Gene Expression Analysis

Per-base sequence quality analysis of raw, paired-end sequencing reads was performed using the Rsubread v2.6.4 [49] package. RNA-seq libraries were aligned to the human reference genome (GRCh38.p13) using the Rsubread package. To quantify gene expression, we used featureCounts from the Rsubread package. The resultant count matrix contained expression counts for 43,689 genes for the 76 RNA-seq libraries.
Differential gene expression analysis was performed using R Statistical Software (v4.0.5) using the packages edgeR v3.40.0 and DESeq2 v1.38.0 [50,51]. We retained 23,129 genes with >10 cpm in at least 7 of the 76 RNA-seq libraries. We calculated and applied TMM-normalisation factors and dispersion to account for differences in library size. Differential gene expression was performed using DESeq2. Genes were considered differentially expressed if the Benjamini–Hochberg adjusted p-value associated with the estimated log2 fold change (log2FC) was below 0.05. Data visualisation was prepared using the ggplot v3.4.2 R package [52].

4.4. Gene Set Testing

GSEA v4.1.0 [53,54] was used to perform pre-ranked gene set analysis. For each comparison, genes were ranked based on calculating the sign of the log2FC multiplied by the negative log10 FDR. The gene sets used for the analysis were derived from the Molecular Signature Database (Human MSigDB v2023.2Hs) [53,55] and include the Hallmark (h.all.v2022.1.Hs.symbols.gmt), Curated (c2.all.v2022.1.Hs.symbols.gmt), Gene Ontology (c5.all.v2022.1.Hs.symbols.gmt), and Immunologic signatures (c7.all.v2022.1.Hs.symbols.gmt). Briefly, gene sets were filtered using the following parameters: a minimum gene size of 15 and a maximum gene size of 2000 with a weighted enrichment statistic. Gene features with no symbol match were omitted. For each comparison, ranked genes were compared to the gene sets using a GSEA normalised enrichment score (NES) and an adjusted p-value. The NES is a normalised enrichment score that measures the degree of enrichment for a gene set in each gene expression dataset. Gene sets were considered significant if they had an adjusted p-value of <0.05.

5. Conclusions

The current body of work demonstrates distinct differences in global transcriptomic signatures that may contribute to altered placental function and fetal growth outcomes in males. Indeed, the data presented herein may suggest activated AR signalling in the placentae of males with reduced birthweight (<10th or 10th–30th BWC). We show that the absence of ARE-containing transcripts results in the enrichment of terms involved in vasculature development, cell adhesion, and immune function in silico. These findings are supported by previous in vitro and in vivo studies that show that the activation of AR signalling pathways in the placenta can impair its function and result in reduced fetal growth outcomes. Therefore, comprehensive characterisation studies of the AR-mediated signalling axis are needed to validate the current in silico work, and to identify novel targets to modulate this pathway that result in improved placental function and fetal growth outcomes, thereby reducing adverse perinatal outcomes, particularly for males.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25031688/s1.

Author Contributions

Conceptualisation, A.S.M. and V.L.C.; methodology, all authors; validation, A.S.M. and M.S.; formal analysis, A.S.M. and M.S.; investigation, all authors; resources, J.L.M., C.T.R., M.L., S.J.E., O.H., A.P., S.A.M., V.F. and V.L.C.; data curation, all authors; writing—original draft preparation, A.S.M.; writing—review and editing, all authors; visualisation, A.S.M. and M.S.; supervision, A.S.M. and V.L.C.; project administration, A.S.M. and V.L.C.; funding acquisition, J.L.M., C.T.R., V.F. and V.L.C. All authors have read and agreed to the published version of the manuscript.

Funding

This initiative falls within the work program of the Stillbirth Centre of Research Excellence. Core funding to support the Stillbirth Centre of Research Excellence is provided by the National Health and Medical Research Council. Funding for this project was also provided by Mater Research and the Mater Foundation. C.T.R. was supported by a National Health and Medical Research Council (NHMRC) Investigator Grant (GNT1174971) and a Matthew Flinders Fellowship from Flinders University.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Mater Human Research Ethics Committee (53157).

Informed Consent Statement

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

Data Availability Statement

The data generated and analysed during this study are available from the corresponding author on reasonable request.

Acknowledgments

The authors acknowledge the support of the Stillbirth Centre of Research Excellence, the research midwives involved with initial participant recruitment and collection, all laboratory members involved with biospecimen processing, and the patients for donating their placentae.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Sharma, D.; Shastri, S.; Sharma, P. Intrauterine Growth Restriction: Antenatal and Postnatal Aspects. Clin. Med. Insights Pediatr. 2016, 10, 67–83. [Google Scholar] [CrossRef]
  2. Gaccioli, F.; Lager, S. Placental Nutrient Transport and Intrauterine Growth Restriction. Front. Physiol. 2016, 7, 40. [Google Scholar] [CrossRef]
  3. Krishna, U.; Bhalerao, S. Placental insufficiency and fetal growth restriction. J. Obstet. Gynaecol. Ind. 2011, 61, 505–511. [Google Scholar] [CrossRef] [PubMed]
  4. Burton, G.J.; Fowden, A.L. The placenta: A multifaceted, transient organ. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2015, 370, 20140066. [Google Scholar] [CrossRef]
  5. Dimasuay, K.G.; Boeuf, P.; Powell, T.L.; Jansson, T. Placental Responses to Changes in the Maternal Environment Determine Fetal Growth. Front. Physiol. 2016, 7, 12. [Google Scholar] [CrossRef]
  6. Yu, J.; Flatley, C.; Greer, R.M.; Kumar, S. Birth-weight centiles and the risk of serious adverse neonatal outcomes at term. J. Perinat. Med. 2018, 46, 1048–1056. [Google Scholar] [CrossRef]
  7. Clifton, V.L. Review: Sex and the human placenta: Mediating differential strategies of fetal growth and survival. Placenta 2010, 31, S33–S39. [Google Scholar] [CrossRef]
  8. Meakin, A.S.; Saif, Z.; Tuck, A.R.; Clifton, V.L. Human placental androgen receptor variants: Potential regulators of male fetal growth. Placenta 2019, 80, 18–26. [Google Scholar] [CrossRef] [PubMed]
  9. Meakin, A.S.; Morrison, J.L.; Bradshaw, E.L.; Holman, S.L.; Saif, Z.; Gatford, K.L.; Wallace, M.J.; Bischof, R.J.; Moss, T.J.; Clifton, V.L. Identification of placental androgen receptor isoforms in a sheep model of maternal allergic asthma. Placenta 2021, 104, 232–235. [Google Scholar] [CrossRef] [PubMed]
  10. Wilson, S.; Qi, J.; Filipp, F.V. Refinement of the androgen response element based on ChIP-Seq in androgen-insensitive and androgen-responsive prostate cancer cell lines. Sci. Rep. 2016, 6, 32611. [Google Scholar] [CrossRef]
  11. Meakin, A.S.; Gough, M.; Saif, Z.; Clifton, V.L. An ex vivo approach to understanding sex-specific differences in placental androgen signalling in the presence and absence of inflammation. Placenta 2022, 120, 49–58. [Google Scholar] [CrossRef] [PubMed]
  12. Beckett, E.M.; Astapova, O.; Steckler, T.L.; Veiga-Lopez, A.; Padmanabhan, V. Developmental programing: Impact of testosterone on placental differentiation. Reproduction 2014, 148, 199–209. [Google Scholar] [CrossRef] [PubMed]
  13. Maliqueo, M.; Echiburú, B.; Crisosto, N. Sex Steroids Modulate Uterine-Placental Vasculature: Implications for Obstetrics and Neonatal Outcomes. Front. Physiol. 2016, 7, 152. [Google Scholar] [CrossRef] [PubMed]
  14. Sun, M.; Maliqueo, M.; Benrick, A.; Johansson, J.; Shao, R.; Hou, L.; Jansson, T.; Wu, X.; Stener-Victorin, E. Maternal androgen excess reduces placental and fetal weights, increases placental steroidogenesis, and leads to long-term health effects in their female offspring. Am. J. Physiol.-Endocrinol. Metab. 2012, 303, E1373–E1385. [Google Scholar] [CrossRef] [PubMed]
  15. Sun, M.; Sun, B.; Qiao, S.; Feng, X.; Li, Y.; Zhang, S.; Lin, Y.; Hou, L. Elevated maternal androgen is associated with dysfunctional placenta and lipid disorder in newborns of mothers with polycystic ovary syndrome. Fertil. Steril. 2020, 113, 1275–1285.e2. [Google Scholar] [CrossRef]
  16. Carlsen, S.M.; Jacobsen, G.; Romundstad, P. Maternal testosterone levels during pregnancy are associated with offspring size at birth. Eur. J. Endocrinol. 2006, 155, 365–370. [Google Scholar] [CrossRef]
  17. Roberts, D.J.; Miyamoto, S. Hexokinase II integrates energy metabolism and cellular protection: Akting on mitochondria and TORCing to autophagy. Cell Death Differ. 2015, 22, 248–257. [Google Scholar] [CrossRef]
  18. Dai, W.; Meng, X.; Mo, S.; Xiang, W.; Xu, Y.; Zhang, L.; Wang, R.; Li, Q.; Cai, G. FOXE1 represses cell proliferation and Warburg effect by inhibiting HK2 in colorectal cancer. Cell Commun. Signal. 2020, 18, 7. [Google Scholar] [CrossRef]
  19. Chen, J.; Zhang, S.; Li, Y.; Tang, Z.; Kong, W. Hexokinase 2 overexpression promotes the proliferation and survival of laryngeal squamous cell carcinoma. Tumour Biol. 2014, 35, 3743–3753. [Google Scholar] [CrossRef]
  20. Liu, W.; Li, W.; Liu, H.; Yu, X. Xanthohumol inhibits colorectal cancer cells via downregulation of Hexokinases II-mediated glycolysis. Int. J. Biol. Sci. 2019, 15, 2497–2508. [Google Scholar] [CrossRef]
  21. Moon, J.S.; Jin, W.J.; Kwak, J.H.; Kim, H.J.; Yun, M.J.; Kim, J.W.; Park, S.W.; Kim, K.S. Androgen stimulates glycolysis for de novo lipid synthesis by increasing the activities of hexokinase 2 and 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 2 in prostate cancer cells. Biochem. J. 2011, 433, 225–233. [Google Scholar] [CrossRef]
  22. Riddle, S.R.; Ahmad, A.; Ahmad, S.; Deeb, S.S.; Malkki, M.; Schneider, B.K.; Allen, C.B.; White, C.W. Hypoxia induces hexokinase II gene expression in human lung cell line A549. Am. J. Physiol. Lung Cell. Mol. Physiol. 2000, 278, L407–L416. [Google Scholar] [CrossRef]
  23. Meng, X.; Grötsch, B.; Luo, Y.; Knaup, K.X.; Wiesener, M.S.; Chen, X.-X.; Jantsch, J.; Fillatreau, S.; Schett, G.; Bozec, A. Hypoxia-inducible factor-1α is a critical transcription factor for IL-10-producing B cells in autoimmune disease. Nat. Commun. 2018, 9, 251. [Google Scholar] [CrossRef]
  24. Gopalakrishnan, K.; Mishra, J.S.; Chinnathambi, V.; Vincent, K.L.; Patrikeev, I.; Motamedi, M.; Saade, G.R.; Hankins, G.D.; Sathishkumar, K. Elevated testosterone reduces uterine blood flow, spiral artery elongation and placental oxygenation in pregnant rats. Hypertension 2016, 67, 630–639. [Google Scholar] [CrossRef] [PubMed]
  25. Sathishkumar, K.; Elkins, R.; Chinnathambi, V.; Gao, H.; Hankins, G.D.; Yallampalli, C. Prenatal testosterone-induced fetal growth restriction is associated with down-regulation of rat placental amino acid transport. Reprod. Biol. Endocrinol. 2011, 9, 110. [Google Scholar] [CrossRef] [PubMed]
  26. Hermans, R.H.; Longo, L.D.; Mcgivern, R.F. Decreased postnatal testosterone and corticosterone concentrations in rats following acute intermittent prenatal hypoxia without alterations in adult male sex behavior. Neurotoxicol. Teratol. 1994, 16, 201–206. [Google Scholar] [CrossRef] [PubMed]
  27. Viana, A., Jr.; Daflon, A.C.; Couto, A.; Neves, D.; De Araujo-Melo, M.H.; Capasso, R. Nocturnal Hypoxemia is Associated With Low Testosterone Levels in Overweight Males and Older Men With Normal Weight. J. Clin. Sleep. Med. 2017, 13, 1395–1401. [Google Scholar] [CrossRef] [PubMed]
  28. Fernandez, E.V.; Reece, K.M.; Ley, A.M.; Troutman, S.M.; Sissung, T.M.; Price, D.K.; Chau, C.H.; Figg, W.D. Dual targeting of the androgen receptor and hypoxia-inducible factor 1α pathways synergistically inhibits castration-resistant prostate cancer cells. Mol. Pharmacol. 2015, 87, 1006–1012. [Google Scholar] [CrossRef] [PubMed]
  29. Mitani, T.; Harada, N.; Nakano, Y.; Inui, H.; Yamaji, R. Coordinated Action of Hypoxia-inducible Factor-1α and β-Catenin in Androgen Receptor Signaling. J. Biol. Chem. 2012, 287, 33594–33606. [Google Scholar] [CrossRef] [PubMed]
  30. Mitani, T.; Yamaji, R.; Higashimura, Y.; Harada, N.; Nakano, Y.; Inui, H. Hypoxia enhances transcriptional activity of androgen receptor through hypoxia-inducible factor-1α in a low androgen environment. J. Steroid Biochem. Mol. Biol. 2011, 123, 58–64. [Google Scholar] [CrossRef] [PubMed]
  31. Park, S.-Y.; Kim, Y.-J.; Gao, A.C.; Mohler, J.L.; Onate, S.A.; Hidalgo, A.A.; Ip, C.; Park, E.-M.; Yoon, S.Y.; Park, Y.-M. Hypoxia Increases Androgen Receptor Activity in Prostate Cancer Cells. Cancer Res. 2006, 66, 5121–5129. [Google Scholar] [CrossRef] [PubMed]
  32. Sebastian, K.; Detro-Dassen, S.; Rinis, N.; Fahrenkamp, D.; Müller-Newen, G.; Merk, H.F.; Schmalzing, G.; Zwadlo-Klarwasser, G.; Baron, J.M. Characterization of SLCO5A1/OATP5A1, a solute carrier transport protein with non-classical function. PLoS ONE 2013, 8, e83257. [Google Scholar] [CrossRef] [PubMed]
  33. Alsinnawi, M.; Zhang, A.; Bianchi-Frias, D.; Burns, J.; Cho, E.; Zhang, X.; Sowalsky, A.; Ye, H.; Slee, A.E.; True, L.; et al. Association of prostate cancer SLCO gene expression with Gleason grade and alterations following androgen deprivation therapy. Prostate Cancer Prostatic Dis. 2019, 22, 560–568. [Google Scholar] [CrossRef] [PubMed]
  34. Wu, J.B.; Shao, C.; Li, X.; Shi, C.; Li, Q.; Hu, P.; Chen, Y.-T.; Dou, X.; Sahu, D.; Li, W.; et al. Near-infrared fluorescence imaging of cancer mediated by tumor hypoxia and HIF1α/OATPs signaling axis. Biomaterials 2014, 35, 8175–8185. [Google Scholar] [CrossRef] [PubMed]
  35. De Gendt, K.; Verhoeven, G. Tissue- and cell-specific functions of the androgen receptor revealed through conditional knockout models in mice. Mol. Cell. Endocrinol. 2012, 352, 13–25. [Google Scholar] [CrossRef] [PubMed]
  36. Flouri, D.; Darby, J.R.T.; Holman, S.L.; Cho, S.K.S.; Dimasi, C.G.; Perumal, S.R.; Ourselin, S.; Aughwane, R.; Mufti, N.; Macgowan, C.K.; et al. Placental MRI Predicts Fetal Oxygenation and Growth Rates in Sheep and Human Pregnancy. Adv. Sci. 2022, 9, 2203738. [Google Scholar] [CrossRef] [PubMed]
  37. Napso, T.; Lean, S.C.; Lu, M.; Mort, E.J.; Desforges, M.; Moghimi, A.; Bartels, B.; El-Bacha, T.; Fowden, A.L.; Camm, E.J.; et al. Diet-induced maternal obesity impacts feto-placental growth and induces sex-specific alterations in placental morphology, mitochondrial bioenergetics, dynamics, lipid metabolism and oxidative stress in mice. Acta Physiol. 2022, 234, e13795. [Google Scholar] [CrossRef] [PubMed]
  38. Song, H.; Telugu, B.P.; Thompson, L.P. Sexual dimorphism of mitochondrial function in the hypoxic guinea pig placenta. Biol. Reprod. 2019, 100, 208–216. [Google Scholar] [CrossRef]
  39. Meakin, A.S.; Cuffe, J.S.M.; Darby, J.R.T.; Morrison, J.L.; Clifton, V.L. Let’s Talk about Placental Sex, Baby: Understanding Mechanisms That Drive Female- and Male-Specific Fetal Growth and Developmental Outcomes. Int. J. Mol. Sci. 2021, 22, 6386. [Google Scholar] [CrossRef]
  40. Ahmad, I.; Newell-Fugate, A.E. Role of androgens and androgen receptor in control of mitochondrial function. Am. J. Physiol.-Cell Physiol. 2022, 323, C835–C846. [Google Scholar] [CrossRef]
  41. Parsons, A.M.; Rajendran, R.R.; Whitcomb, L.A.; Bouma, G.J.; Chicco, A.J. Characterization of trophoblast mitochondrial function and responses to testosterone treatment in ACH-3P cells. Placenta 2023, 137, 70–77. [Google Scholar] [CrossRef] [PubMed]
  42. Holland, O.; Nitert, M.D.; Gallo, L.A.; Vejzovic, M.; Fisher, J.J.; Perkins, A.V. Review: Placental mitochondrial function and structure in gestational disorders. Placenta 2017, 54, 2–9. [Google Scholar] [CrossRef] [PubMed]
  43. Koster, M.P.; de Wilde, M.A.; Veltman-Verhulst, S.M.; Houben, M.L.; Nikkels, P.G.; van Rijn, B.B.; Fauser, B.C. Placental characteristics in women with polycystic ovary syndrome. Human. Reprod. 2015, 30, 2829–2837. [Google Scholar] [CrossRef] [PubMed]
  44. Kelley, A.S.; Puttabyatappa, M.; Ciarelli, J.N.; Zeng, L.; Smith, Y.R.; Lieberman, R.; Pennathur, S.; Padmanabhan, V. Prenatal Testosterone Excess Disrupts Placental Function in a Sheep Model of Polycystic Ovary Syndrome. Endocrinology 2019, 160, 2663–2672. [Google Scholar] [CrossRef] [PubMed]
  45. Lim, W.F.; Forouhan, M.; Roberts, T.C.; Dabney, J.; Ellerington, R.; Speciale, A.A.; Manzano, R.; Lieto, M.; Sangha, G.; Banerjee, S.; et al. Gene therapy with AR isoform 2 rescues spinal and bulbar muscular atrophy phenotype by modulating AR transcriptional activity. Sci. Adv. 2021, 7, eabi6896. [Google Scholar] [CrossRef]
  46. Ahrens-Fath, I.; Politz, O.; Geserick, C.; Haendler, B. Androgen receptor function is modulated by the tissue-specific AR45 variant. FEBS J. 2005, 272, 74–84. [Google Scholar]
  47. Costanzo, V.; Bardelli, A.; Siena, S.; Abrignani, S. Exploring the links between cancer and placenta development. Open Biol. 2018, 8, 180081. [Google Scholar] [CrossRef]
  48. Burton, G.; Sebire, N.; Myatt, L.; Tannetta, D.; Wang, Y.-L.; Sadovsky, Y.; Staff, A.; Redman, C. Optimising sample collection for placental research. Placenta 2014, 35, 9–22. [Google Scholar] [CrossRef]
  49. Liao, Y.; Smyth, G.K.; Shi, W. The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads. Nucleic Acids Res. 2019, 47, e47. [Google Scholar] [CrossRef]
  50. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [PubMed]
  51. 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]
  52. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016. [Google Scholar]
  53. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [PubMed]
  54. Mootha, V.K.; Lindgren, C.M.; Eriksson, K.-F.; Subramanian, A.; Sihag, S.; Lehar, J.; Puigserver, P.; Carlsson, E.; Ridderstråle, M.; Laurila, E.; et al. PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat. Genet. 2003, 34, 267–273. [Google Scholar] [CrossRef] [PubMed]
  55. Liberzon, A.; Birger, C.; Thorvaldsdóttir, H.; Ghandi, M.; Mesirov, J.P.; Tamayo, P. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 2015, 1, 417–425. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Differentially expressed transcripts from males of different BWC outcomes. (A) Volcano plot visualisation of differentially expressed transcripts from the >30th BWC vs. <10th BWC and the 10th–30th vs. <10th BWC comparisons. Red data points represent transcripts that were significantly upregulated and blue data points represent transcripts that were significantly downregulated in >30th or 10th–30th BWC males compared with <10th BWC males; grey data points represent non-significant transcripts. (B) Venn diagram representing differentially expressed placental transcripts ( = upregulated; = downregulated) between males with >30th, 10th–30th, and <10th BWC outcomes. An adjusted p-value of 0.05 was used to identify transcripts that were differentially expressed between the BWC subcategories of male placentae.
Figure 1. Differentially expressed transcripts from males of different BWC outcomes. (A) Volcano plot visualisation of differentially expressed transcripts from the >30th BWC vs. <10th BWC and the 10th–30th vs. <10th BWC comparisons. Red data points represent transcripts that were significantly upregulated and blue data points represent transcripts that were significantly downregulated in >30th or 10th–30th BWC males compared with <10th BWC males; grey data points represent non-significant transcripts. (B) Venn diagram representing differentially expressed placental transcripts ( = upregulated; = downregulated) between males with >30th, 10th–30th, and <10th BWC outcomes. An adjusted p-value of 0.05 was used to identify transcripts that were differentially expressed between the BWC subcategories of male placentae.
Ijms 25 01688 g001
Figure 2. Venn diagram visualisation of upregulated and downregulated GSEA terms in (A) <10th vs. >30th, (B) <10th vs. 10th–30th, and (C) 10th–30th vs. >30th. ARE positive = transcripts containing only AREs; ARE negative = transcripts without AREs.
Figure 2. Venn diagram visualisation of upregulated and downregulated GSEA terms in (A) <10th vs. >30th, (B) <10th vs. 10th–30th, and (C) 10th–30th vs. >30th. ARE positive = transcripts containing only AREs; ARE negative = transcripts without AREs.
Ijms 25 01688 g002
Table 1. Maternal and neonatal characteristics.
Table 1. Maternal and neonatal characteristics.
BWC
<10th10th–30th>30th
n = 16n = 12n = 9p
Maternal age (yrs)32.0 (29.5–34.5)33.5 (30.0–35.0)32.0 (29.0–33.0)0.638
BMI20.9 (19.3–21.6)26.5 (21.9–33.2)21.6 (21.1–31.5)0.097
Parity1.0 (0.0–1.0)2.0 (1.0–2.0)1.0 (1.0–2.0)0.106
Gestation (days)262.5 (260.5–267.0) a277.5 (268.0–283.5) b274.0 (267.0–278.0) b0.004
BWC (%)2.5 (2.5–5.2) a20.0 (17.6–25.1) b91.1 (79.2–95.6) b0.000
Birthweight (g)2468.0 (2232.5–2542.5) a3225.0 (3035.0–3375.0) b3900.0 (3760.0–4196.0) b0.000
Placenta weight (g)520.0 (380.0–576.4) a550.0 (510.0–592.0) ab729.5 (600.0–890.0) b0.000
% Asthma12.516.711.1-
% GDM6.300-
% Pre-eclampsia6.300-
Bolded values represent a significant effect from the test variable; different superscripts represent a significant pairwise comparison post hoc test between them. BMI = body mass index; BWC = birthweight centile; GDM = gestational diabetes mellitus.
Table 2. Differentially expressed transcripts between male BWC subcategories contain ARE full and half sites.
Table 2. Differentially expressed transcripts between male BWC subcategories contain ARE full and half sites.
BWCTranscriptlog2 Fold Change
(Relative to <10th)
Adjusted p-ValueARE Full SiteARE Half Site
10th–30thLEP−3.503.49 × 10−3
LOC105370135−3.435.79 × 10−3××
SLCO5A1 1−3.412.91 × 10−2
HK2P1−3.013.49 × 10−3××
HK2 1−2.701.84 × 10−3
AOX3P-AOX2P−2.604.50 × 10−2××
FSTL3−2.402.91 × 10−2×
PNCK−2.011.27 × 10−2××
ARHGAP45−1.749.99 × 10−3××
PDZD7−1.711.30 × 10−2××
ARNT2−1.672.23 × 10−2×
SPAG4−1.492.90 × 10−3
FLT1−1.372.91 × 10−2
CORO2A−1.374.61 × 10−2
HILPDA−1.373.49 × 10−3
PLIN2−1.242.23 × 10−2××
LOC100506358−1.231.30 × 10−2××
SFXN3−1.203.30 × 10−3××
TPBG−1.182.76 × 10−3
LOC102724660−1.141.27 × 10−2××
NDRG1−1.122.91 × 10−2
MGAT3−1.031.30 × 10−2
LIMD1−1.011.27 × 10−2
RH13−1.012.91 × 10−2××
NCOR2−0.984.61 × 10−2
LOC107987470−0.933.51 × 10−2××
MXI1−0.903.69 × 10−2
SYDE1−0.874.61 × 10−2××
EPS8L1−0.822.19 × 10−2
DGKZ−0.663.55 × 10−2×
OBSL1−0.631.84 × 10−3×
LOC1079842441.154.19 × 10−2××
LOC1053788422.022.91 × 10−2××
XIST3.539.76 × 10−9×
>30thPTGS2−2.781.35 × 10−2
BTNL9−2.384.92 × 10−2×
PNCK−2.103.29 × 10−2××
PROK1−1.904.92 × 10−2
IGF2−1.653.29 × 10−2
ACTG2−1.614.92 × 10−2
ARNT2−1.614.92 × 10−2×
PFKP−1.604.92 × 10−2
AFAP1−1.462.67 × 10−2
MRC2−1.414.55 × 10−2
TUBB2A−1.404.92 × 10−2×
NUAK1−1.383.39 × 10−2
P2RX1−1.353.45 × 10−2××
NMUR1−1.264.92 × 10−2×
SYNPO2−1.214.92 × 10−2
SALL1−1.174.92 × 10−2
MXI1−0.993.46 × 10−2
PRKD3−0.944.92 × 10−2×
SEMA3G−0.944.92 × 10−2×
RUSC2−0.884.92 × 10−2
FGD1−0.854.92 × 10−2
S100A16−0.654.92 × 10−2
PIP4K2A−0.564.92 × 10−2
MYH9−0.414.92 × 10−2
DARS1−0.373.39 × 10−2××
LOC1122680560.794.92 × 10−2××
1 = transcripts regulated by ARE half and/or full sites.
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

Meakin, A.S.; Smith, M.; Morrison, J.L.; Roberts, C.T.; Lappas, M.; Ellery, S.J.; Holland, O.; Perkins, A.; McCracken, S.A.; Flenady, V.; et al. Placenta-Specific Transcripts Containing Androgen Response Elements Are Altered In Silico by Male Growth Outcomes. Int. J. Mol. Sci. 2024, 25, 1688. https://doi.org/10.3390/ijms25031688

AMA Style

Meakin AS, Smith M, Morrison JL, Roberts CT, Lappas M, Ellery SJ, Holland O, Perkins A, McCracken SA, Flenady V, et al. Placenta-Specific Transcripts Containing Androgen Response Elements Are Altered In Silico by Male Growth Outcomes. International Journal of Molecular Sciences. 2024; 25(3):1688. https://doi.org/10.3390/ijms25031688

Chicago/Turabian Style

Meakin, Ashley S., Melanie Smith, Janna L. Morrison, Claire T. Roberts, Martha Lappas, Stacey J. Ellery, Olivia Holland, Anthony Perkins, Sharon A. McCracken, Vicki Flenady, and et al. 2024. "Placenta-Specific Transcripts Containing Androgen Response Elements Are Altered In Silico by Male Growth Outcomes" International Journal of Molecular Sciences 25, no. 3: 1688. https://doi.org/10.3390/ijms25031688

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