ORIGINAL RESEARCH article
Sec. Sex and Gender in Cardiovascular Medicine
Volume 8 - 2021 | https://doi.org/10.3389/fcvm.2021.693351
Plasma Testosterone Levels and Atherosclerotic Plaque Gene Expression in Men With Advanced Atherosclerosis
- 1Laboratory of Experimental Cardiology, University Medical Center Utrecht, Utrecht University, Utrecht, Netherlands
- 2Central Diagnostic Laboratory, University Medical Center Utrecht, Utrecht University, Utrecht, Netherlands
- 3Division of Surgical Specialties, Department of Vascular Surgery, University Medical Centre Utrecht, Utrecht University, Utrecht, Netherlands
- 4Julius Center for Health Sciences and Primary Care, University Medical Center Utrecht, Utrecht University, Utrecht, Netherlands
Aims: Low plasma testosterone levels have been shown to predict worse outcome in men with severe atherosclerotic disease. We hypothesized that a low plasma testosterone level affects disease risk through changes in gene expression in atherosclerotic plaques. Therefore, we studied plasma testosterone levels in relation to gene expression levels in atherosclerotic plaque tissue of men with advanced atherosclerotic disease.
Methods: Plasma testosterone levels were measured in 203 men undergoing carotid endarterectomy. The corresponding atherosclerotic plaque tissue was used for RNA sequencing. First, we assessed how often the androgen receptor gene was expressed in the plaque. Second, correlations between plasma testosterone levels and pre-selected testosterone-sensitive genes were assessed. Finally, differences within the RNA expression profile of the plaque as a whole, characterized into gene regulatory networks and at individual gene level were assessed in relation to testosterone levels.
Results: Testosterone plasma levels were low with a median of 11.6 nmol/L (IQR: 8.6–13.8). RNA-seq of the plaque resulted in reliable expression data for 18,850 genes to be analyzed. Within the RNA seq data, the androgen-receptor gene was expressed in 189 out of 203 (93%) atherosclerotic plaques of men undergoing carotid endarterectomy. The androgen receptor gene expression was not associated with testosterone plasma levels. There were no significant differences in gene expression of atherosclerotic plaques between different endogenous testosterone levels. This remained true for known testosterone-sensitive genes, the complete transcriptomic profile, male-specific gene co-expression modules as well as for individual genes.
Conclusion: In men with severe atherosclerotic disease the androgen receptor is highly expressed in plaque tissue. However, plasma testosterone levels were neither associated with pre-selected testosterone sensitive genes, gene expression profiles nor gene regulatory networks in late-stage atherosclerotic plaques. The effect of testosterone on gene expression of the late-stage atherosclerotic plaque appears limited, suggesting that alternate mechanisms explain its effect on clinical outcomes.
Testosterone is the most important sex hormone in men. The decline in testosterone levels with age in men has been linked to increased prevalence of (coronary) artery disease and cardiovascular events (1). In men, higher plasma testosterone levels have shown to protect against cardiovascular disease (2, 3). Therefore, supplementation of testosterone in men with low testosterone levels was considered to be beneficial. However, treatment of low testosterone levels in men have resulted in ambiguous results (4), stressing the need to unravel the underlying mechanisms responsible for the detrimental effect of low endogenous testosterone in men.
Atherosclerosis underlies the majority of cardiovascular disease (5). The composition of the plaque is an important indicator of cardiovascular events (6). Animal studies have shown a beneficial effect of testosterone on atherosclerotic plaque formation, suspectedly via its anti-inflammatory properties (7, 8). In humans, even though the relationship between endogenous testosterone levels and cardiovascular disease has extensively been studied (3), the role of testosterone in the development and composition of the plaque is still unknown.
Testosterone is the main endogenous androgen (9) and is involved in many physiological processes, among which the function of the vascular endothelium. For example, by binding to the nuclear androgen receptor, testosterone can alter downstream gene expression resulting in vasodilatation. To exert its regulatory effect, testosterone needs to bind to the androgen receptor (10), this receptor is nuclear and the predominant mediator of the effects of endogenous testosterone. For example, androgen receptor ablation in the vascular smooth muscle cell inhibits vascular calcification (11). The androgen receptor regulates downstream effects of testosterone as a transcription factor and it is expressed in many tissues, including the vascular endothelium. However, besides genomic actions of testosterone on the vasculature, it is thought to also have non-genomic vasoactive effects which might be exerted through interaction with proteins, receptors, or ion channels nested in the plasma membrane (9).
As we previously showed that plasma testosterone was a strong predictor of a secondary stroke in our cohort of men undergoing carotid endarterectomy (12), we hypothesized that the mechanism by which low testosterone levels affect cardiovascular risk is through the atherosclerotic plaque. To test this hypothesis, we studied if plasma testosterone associates with plaque tissue gene expression. We used several approaches including analyses of complete transcriptomic profiles as well as gene co-expression networks as these are able to capture relatively subtle changes in overall expression within tissues (13). We also studied if pre-selected testosterone-sensitive genes and receptors were associated with plasma testosterone levels.
We analyzed transcriptional activity by RNA-sequencing of plaques from men included in the Athero-Express Biobank Study (14) of whom testosterone level and RNA-sequencing data were available (n = 203). Details of the Athero-Express study protocol have been described previously (12, 15). In short, patients undergoing endarterectomy of the carotid artery in two Dutch tertiary referral centers between 2002 and 2015 were included in this study. Study procedures comprise of a baseline blood withdrawal, an extensive questionnaire filled in by the participants verified against medical records, and collection of carotid arterial plaque material during surgery. All patients provided written informed consent before surgery, the study was approved by the Local Medical Ethical Committee and conducted according to the Declaration of Helsinki (16).
Sex Steroid Measurements
Serum testosterone was measured by immunoassay on an ARCHITECT ci8200 system (Abbott Diagnostics, Abbott laboratories, USA) (12).
As described previously (12, 14, 17), the atherosclerotic plaque was processed directly after surgery and (immune-) histochemical staining was routinely performed on the culprit lesion (segment with the highest plaque burden) for identification of macrophages (CD 68), calcification (haematoxylin-eosin (HE)), smooth muscle cells (alfa actin), collagen [Picro Sirius red (PSR)], intra-plaque hemorrhage (HE, Elastin von Gieson staining), vessel density (CD34) and fat (PSR, HE).
Bulk RNA Sequencing
As the culprit lesion is used for plaque histology following the standardized Athero-Express protocol (14), the adjacent plaque segments were used for RNA sequencing. To measure bulk RNA expression in the plaques total RNA was isolated according to the manufacturers protocol after processing of the plaque segments using ceramic beads and tissue homogenizer (Precellys, Bertin intruments, Montigny-le-Bretonneux) with use of TriPure (Sigma Aldrich). After precipitating RNA in the aqueous phase with propanolol, RNA was washed with 75% ethanol and either used immediately after an additional washing step with 75% ethanol or stored in 75% ethanol for later use. Subsequently, library preparation was performed as described before (18–20). Ethanol was removed and the pellet air-dried. Then, primer mix (5 ng primer per reaction) was added to initiate primer annealing at 65 degrees Celsius for 5 min. Subsequently, reverse transcription (RT) was executed. Subsequent reverse-transcription RT reaction; first strand reaction for 1 h at 42°C, heat inactivated for 10 min at 70°C, second strand reaction for 2 h at 16°C, and then put on ice until proceeding to sample pooling.
This initial RT reaction used the following primer design: an anchored polyT, a unique 6 bp barcode, a unique molecular identifier (UMI) of 6 bp, the 5' Illumina adapter and a T7 promoter, as described (18–20). Each sample now contained its own unique barcode making it possible to pool together cDNA samples at 7 samples per pool. Complementary DNA (cDNA) was cleaned using AMPure XP beads (Beckman Coulter), washed with 80% ethanol, and resuspended in water before proceeding to the in vitro transcription (IVT) reaction (AM1334; Thermo-Fisher) incubated at 37°C for 13 h. Next, Exo-SAP (Affymetrix, Thermo-Fisher) was used to remove primers, upon which amplified RNA (aRNA) was fragmented, cleaned with RNAClean XP (Beckman-Coulter), washed with 70% ethanol, air-dried, and resuspended in water. RNA yield and quality in the suspension were checked by Bioanalyzer (Agilent) after removal of the beads with use of a magnetic stand. By performing an RT reaction using SuperScript II reverse transcriptase (Invitrogen/Thermo-Fisher) according to the protocol of the manufacturer cDNA library construction was initiated. Next, PCR amplification was performed as described previously (18–20). PCR products were cleaned twice using AMPure XP beads (Beckman Coulter). Qubit fluorometric quantification (Thermo-Fisher) and Bioanalyzer (Agilent) were used to checked Library cDNA yield and quality. Illumina Nextseq500 platform was used to sequence the libraries; paired end, 2 × 75 bp. After sequencing, retrieved fastq files were de-barcoded and split into forward and reverse reads. From there, the reads were mapped using Burrows-Wheel aligner (BWA37) version 0.7.17–r1188, calling “bwa aln” with settings -B 6 -q 0 -n 0.00 -k 2 -l 200 -t 6 for R1 and -B 0 -q 0 -n 0.04 -k 2 -l 200 -t 6 for R2, “bwa sampe” with settings -n 100 -N 100, and a cDNA reference (assembly hg19, Ensembl release 84). Read and UMI counts were acquired from SAM files with use of custom perl code and collected into count matrices. Further analyses were performed using R (21) version 3.6.2 and later and its IDE Rstudio (22) version 1.2 and later. Genes were annotated with Ensembl ID's, basic quality control was performed [filtering out samples with low gene numbers (<10,000 genes) and read counts (<18,000 reads)].
Counts, metadata and clinical data were combined into a SummarizedExperiment (23) object. Counts were pre-filtered, normalized and transformed making use of the variance stabilization transformation function (vst) in DESeq2 (24), resulting in transformed data of n = 203 on a log2-scale, normalized for library size (24). These were used to visualize gene expression, to assess correlations between gene expression and androgen receptor expression and plasma testosterone levels, to construct correlation heatmaps, correlation scatter plots, and Manhattan plots, furthermore we used these to calculate differentially expressed genes between high and low testosterone tertiles. Correlation estimates and p-values were calculated using Spearman's rank correlation. Correlations between the expression levels of genes and testosterone level were explored using scatterplots. Heatmaps were drawn using Pheatmap from the Pheatmap package (25), applying hierarchical clustering based on correlation estimates with standard settings: complete linkage and Euclidean distance.
Gene Regulatory Networks Analysis
The software package WGCNA (13) was used to generate modules of co-expression genes on the 486 available men RNAseq samples. After excluding all the ribosomal genes and included only the protein-coding genes with annotated HGCN names, a set of 12,765 protein-coding genes which passed quality control (average >1 count per sample) was used for module generation. The raw read counts were corrected for UMI sampling [corrected_count = −096*(ln (1–(raw_count/4096)))], normalized by sample sequencing depth and log-transformed. A signed network was constructed using the robust “bicor” correlation measure. To determine the exponent used for the adjacency matrix construction, soft thresholding analysis was performed with the WGCNA package for powers ranging from 2 to 30. The cut-off for assuming scale-free topology was set at an R-squared of 0.8, while having a median connectivity of lower than 100.The chosen lowest complying power was 26. The network was constructed by first generating an adjacency matrix, which was transformed into a topological overlay matrix (13). Modules were detected by clustering the average distance of the dissimilarity matrix (defined as 1- topological overlay matrix) and cutting the subsequent dendrogram by using the cutreeDynamic function (deepSplit = 2, minClusterSize = 20). Module eigengenes were calculated by taking the first principal component of gene expression in that module. Module eigengene values were selected for the 203 samples with testosterone measurements available. The module eigengenes were correlated to clinical traits by Pearson correlation, with a Student asymptomatic p-value test for significance.
Data and Scripts Availability
Data and scripts are available upon reasonable request.
We analyzed data of 203 men with a mean age of 68 years old. Most men presented with a transient ischemic attack (42%) or stroke (23%). Median testosterone level was 11.6 nmol/L (IQR: 8.6–13.8). The participants of the study population were allocated to tertiles based on testosterone levels: high (median 16.5, IQR: 14.9–18.9 nmol/L), medium (11.6, IQR: 10.6–12.7 nmol/L), and low (7.6, IQR: 8.6–13.8 nmol/L). Patient characteristics of our study population stratified by testosterone tertiles are presented in Table 1. Age, body mass index and renal function did not differ between the tertiles. Hypercholesterolemia seemed more prevalent in the low and medium testosterone level group (73 and 72%, respectively) as compared to the group with the highest testosterone levels (60%), yet this difference was not statistically significant. Hypertension was most prevalent in individuals with the lowest testosterone levels [83% in the low, 69% in the medium and 59% in the highest testosterone tertile, respectively (p < 0.05)]. In all groups, a transient ischemic attack was the most common clinical presentation and there was a balanced distribution of plaque phenotypes based on histology with similar prevalence of atheromatous, fibro-atheromatous, and fibrous plaques (17).
Table 1. Baseline characteristics by testosterone tertiles of men undergoing carotid endarterectomy.
Androgen Receptor Expression Within the Atherosclerotic Plaque
To exert its function, testosterone needs to bind to the androgen receptor (AR) (10) which mediates almost all the known genomic effects of testosterone. To assess the level of AR expression within the atherosclerotic plaque we compared its expression to the mean expression of 1,000 randomly selected genes (Figure 1A), the mean expression of the AR gene was 2.2, SD: 0.7 vs. 1.1, SD: 0.1 for the average of 1,000 genes (log transformed counts). We found that the AR gene is expressed within the atherosclerotic plaques of 93% (189 of 203 patients) in our male cohort (count > 1). AR expression levels were not significantly correlated to testosterone levels in serum of the same patients (Figure 1B) (Spearman rank correlation (rho): 0.08, unadjusted p-value = 0.25). In addition, we explored the expression of other nuclear receptors in relationship to the AR gene, several receptors presented with comparable expression levels with NR6A, HNF4A, and THRB presenting the highest expression in this cohort (Supplementary Figure 1).
Figure 1. Expression of the Androgen Receptor (AR) gene and its association with testosterone levels in men undergoing carotid endarterectomy. (A) AR gene expression levels as compared with the mean expression level of 1,000 randomly selected genes. (B) Scatterplot of testosterone level (rank-normalized) and AR gene expression in plaque.
Testosterone-Sensitive Genes and Plasma Testosterone Levels
A list of testosterone-sensitive genes was selected based on literature (26). In addition, an unbiased list of genes was selected based on significant correlation (after Bonferroni adjustment) of their gene expression with expression of the AR gene (Supplementary Figure 2, Supplementary Table 1), 12 out of these 20 genes (60%) were identified as being targets of the AR gene as a transcription factor (27). The combination of these two lists of candidate testosterone-sensitive genes was used, and an association of their expression with testosterone levels within the plaque was further evaluated (Table 2). Correlation analysis of the testosterone-sensitive genes' expression (Table 2) and continuous testosterone levels showed no significant correlations (Figure 2). Clustering of individuals using hierarchical clustering and based on sample distances restricted to the testosterone-sensitive gene subset did not show an overlap with the established testosterone-level tertiles (Supplementary Figure 3).
Figure 2. Scatterplots for testosterone levels and gene expression of selected hormone-dependent genes based on literature and gene co-expression with androgen receptor (AR) gene in the atherosclerotic plaque.
Plaque Transcriptional States Do Not Correlate With Testosterone Levels
Hierarchical clustering based on sample distances did not highlight a group of patients who are similar with respect to their specific testosterone levels (Figure 3A). In order to identify individual genes associated with testosterone levels, correlations between gene expression from each gene with testosterone levels were tested. No gene was statistically significantly correlated with testosterone (after Bonferroni correction). The top 15 genes nominally correlated with testosterone levels are labeled in a Manhattan plot in Figure 3B and presented in Supplementary Table 2. We also performed differential expressed analysis between high and low testosterone levels. Yet, there were no differentially expressed genes (data not shown).
Figure 3. Association of testosterone levels with transcriptome signatures in male plaque tissue. (A) Heatmap based on sample-to-sample distances colored by continuous testosterone levels. (B) Manhattan plot of associations between testosterone levels and individual gene expression.
No Evidence of an Effect of Testosterone Levels on Male-Specific Gene Co-expression Modules
Gene co-expression modules were constructed following previously described methods (28, 29) (Supplementary Figure 4), identifying 14 modules (Figure 4A) in the bulk plaque RNA-seq data. Eigen genes were calculated for each individual module, representing a proxy for the combined expression of the genes that belong to each module. Using this composite measure of module expression, we correlated each module with the rank-normalized testosterone levels in the 203 patients with data available. Two modules (purple and turquoise) were identified with nominally significant negative correlations (p-value < 0.1) (Figure 4B). However, no significant associations remained after adjusting the p-values using Holm method for false discovery rate (30). To understand the biological relevance of such modules, gene enrichment analysis was performed and the modules that were nominally associated were pointing to endoplasmic reticulum processes (Figure 4C).
Figure 4. Association analysis of gene co-expression networks with testosterone levels. (A) Hierarchical cluster tree of 12,756 genes, height (y-axis) is determined using the average gene linkage from the topological overlap measure, which incorporates information about both co-expression and connectivity in the network. A total of 14 modules were assigned based on the dynamic hybrid branch cutting method and are shown along the x-axis in different colors. Genes in the same module have higher co-expression and present higher connection in the network. (B) Correlation of module eigengenes (a composite measure of the combined expression of the genes from each module) with testosterone levels and age of the participants, p-values for the correlation are in parenthesis. (C) Gene enrichment GO terms obtained for the purple module, nominally associated with testosterone levels (no significant enrichment was found for turquoise module).
By using statistical and system biology approaches, we show that gene expression of the plaque tissue of men with severe atherosclerosis is unrelated to plasma testosterone levels, despite a high relative expression of the androgen receptor.
Testosterone level itself was not correlated with the overall transcriptional state of the atherosclerotic plaques in our study and did not significantly influence gene expression at individual gene level. Furthermore, dissecting the gene expression profiles into gene co-expression modules, which are known to reveal subtle differences in tissue gene expression (29) also showed no effect of testosterone on gene regulatory networks. Also, known testosterone-sensitive (26) genes as well as genes that were correlated to the androgen receptor expression were not correlated to plasma testosterone levels.
Several factors might explain our unexpected results. First of all, levels in our study were overall low as our population consists mainly of older men. It may well be that testosterone does affect plaque expression, and thereby affects atherosclerosis, but rather at an earlier age when testosterone levels are higher. Despite the low levels, we did observe a strong relation between sex-hormone levels and secondary outcomes (12) in this cohort. Low testosterone to estradiol predicted future major cardiovascular events (hazard ratio (HR) 1.67; 95% CI: 1.02–2.76, an effect even stronger in obese men (HR 2.42; 95% CI: 1.09–5.38), This coincided with an unfavorable inflammatory pattern shown by elevated hsCRP levels.
It may be that the systemic inflammation we observed is an important factor in secondary outcome in these men. Indeed, it is known that stroke patients with a high inflammatory profile have a poor outcome (31–33), possibly by time-dependent recruitment and activation of inflammatory cells at the site of the injury. Indeed, testosterone is known to affect the immune system (34). Interestingly, a high prevalence of venous thrombo-embolic complications has been documented in Klinefelter Syndrome (35), characterized by genetically low testosterone and often accompanied by abdominal adiposity and metabolic syndrome. It could well be that the relation between low testosterone and secondary outcome is mediated by changes in vascular hemostasis and thrombosis. Testosterone is known to affect platelets (36), but an effect on thrombosis has not been found (37).
From a disease progression point of view, testosterone could be more important in the development and progression of atherosclerotic disease in earlier stage of the disease, instead of in older patients with advanced atherosclerosis as we have studied here. Animal studies indeed point toward a beneficial effect of testosterone during the development of the atherosclerotic plaque (7, 8), but longitudinal studies in humans are inconclusive as excellently reviewed elsewhere (38). The premature interruption of the Testosterone in Older Men With Mobility Limitations (TOM) study raised serious concerns on the potential increased risk of CVD upon testosterone supplementation (39).
Limitations of the Study
Testosterone levels were available in 203 men and our study may not be sufficiently powered to assess associations of each individual gene with plasma testosterone levels. However, a predefined testosterone responsive list of genes also did not reveal any association.
In literature, we searched for candidate genes that are influenced by the androgen receptor. Expression levels of these genes were associated with plaque androgen receptor expression and testosterone levels. Publicly available data was scarce and limited to expression levels obtained in endothelial progenitor cells. Although the endothelial cell is considered to be a central player in the effects of testosterone on the vessel wall, we cannot rule out the possibility that co-expression of relevant genes can be found in other cell clusters present in atherosclerotic lesions.
In our study, the expression of the androgen receptor does show a correlation with downstream gene expression, hinting toward an association between the testosterone pathway and down-stream gene expression, but was independent of plasma testosterone levels. As dihydrotestosterone (DHT) is the most active androgen, for future studies, it might be interesting to assess the association between different DHT levels and the plaque transcriptome, instead of total testosterone alone. Unfortunately, we did not measure DHT in this study.
Lastly, testosterone levels depend on aromatase-activity, as testosterone is converted to estradiol (1). This aromatization is dependent on white adipose tissue, and thereby obese subjects often have lower testosterone levels. Thus, after conversion to estradiol, testosterone can theoretically also have an effect on plaque expression indirectly via estrogen-receptors, a link that needs to be studied further.
In men with severe atherosclerotic disease the androgen receptor is highly expressed in plaque tissue. However, plasma testosterone levels were neither associated with gene expression profiles nor with gene regulatory networks or with pre-selected testosterone sensitive genes in late-stage atherosclerotic plaques. The effect of testosterone on gene expression of the late-stage atherosclerotic plaque appears limited, suggesting that alternative mechanisms explain its effect on clinical outcomes.
Data Availability Statement
The datasets generated for this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
The studies involving human participants were reviewed and approved by University Medical Center Utrecht Medical Research Ethical Committee. The patients/participants provided their written informed consent to participate in this study.
FG undertook the writing of the paper and initiated the analysis plan. EDB guided the writing of the paper and analyzed the data. AB analyzed the data and made substantial improvements to the paper. NT, FW, RH and GP made substantial improvements to the paper. NO-M, GP and HR supervised the study and made substantial improvements to the paper. All authors have read and given final approval of the submitted manuscript.
FG is supported by funding from CVON 2017-22 ARGUS. HD is supported by funding from The DCVA (2020B004—IMPRESS), an ERC consolidator Grant (866478—UCARE).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2021.693351/full#supplementary-material
1. Kloner RA, Carson C, Dobs A, Kopecky S, Mohler ER. Testosterone and cardiovascular disease. J Am Coll Cardiol. (2016) 67:545–57. doi: 10.1016/j.jacc.2015.12.005
2. Mathur P, Ostadal B, Romeo F, Mehta JL. Gender-related differences in atherosclerosis. Cardiovasc Drugs Ther. (2015) 29:319–27. doi: 10.1007/s10557-015-6596-3
3. Moreau KL, Babcock MC, Hildreth KL. Sex differences in vascular aging in response to testosterone. Biol Sex Differ. (2020) 11:18. doi: 10.1186/s13293-020-00294-8
4. Elagizi A, Köhler TS, Lavie CJ. Testosterone and cardiovascular health. Mayo Clin Proc. (2018) 93:83–100. doi: 10.1016/j.mayocp.2017.11.006
5. Lusis AJ. Atherosclerosis. Nature. (2000) 407:233–41. doi: 10.1038/35025203
6. Corti R, Hutter R, Badimon JJ, Fuster V. Evolving concepts in the triad of atherosclerosis, inflammation and thrombosis. J Thromb Thrombolysis. (2004) 17:35–44. doi: 10.1023/B:THRO.0000036027.39353.70
7. Malkin CJ, Pugh PJ, Jones RD, Jones TH, Channer KS. Testosterone as a protective factor against atherosclerosis-immunomodulation and influence upon plaque development and stability. J Endocrinol. (2003) 178:373–80. doi: 10.1677/joe.0.1780373
8. Li SJ, Li XY, Li Y. Regulation of atherosclerotic plaque growth and stability by testosterone and its receptor via influence of inflammatory reaction. Vascul Pharmacol. (2008) 49:14–8. doi: 10.1016/j.vph.2008.03.004
9. Lorigo M, Mariana M, Lemos MC, Cairrao E. Vascular mechanisms of testosterone: the non-genomic point of view. J Steroid Biochemistry and Molecular Biology. (2020) 196:105496. doi: 10.1016/j.jsbmb.2019.105496
10. Matsumoto T, Sakari M, Okada M, Yokoyama A, Takahashi S, Kouzmenko A, et al. The androgen receptor in health and disease. Annu Rev Physiol. (2013) 75:201–24. doi: 10.1146/annurev-physiol-030212-183656
11. Zhu D, Hadoke PWF, Wu J, Vesey AT, Lerman DA, Dweck MR, et al. Ablation of the androgen receptor from vascular smooth muscle cells demonstrates a role for testosterone in vascular calcification. Sci Rep. (2016) 6:1–9. doi: 10.1038/srep24807
12. Van Koeverden ID, De Bakker M, Haitjema S, Van Der Laan SW, De Vries JPPM, Hoefer IE, et al. Testosterone to oestradiol ratio reflects systemic and plaque inflammation and predicts future cardiovascular events in men with severe atherosclerosis. Cardiovasc Res. (2019) 115:453–62. doi: 10.1093/cvr/cvy188
13. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformat. (2008) 9:559. doi: 10.1186/1471-2105-9-559
14. Verhoeven BAN, Velema E, Schoneveld AH, Paul J, De Vries PM, De Bruin P, et al. Athero-Express: differential atherosclerotic plaque expression of mRNA and protein in relation to cardiovascular events and patient characteristics. Rationale and design. Eur J Epidemiol. (2004) 19:1127–33. doi: 10.1007/s10564-004-2304-6
15. Van Der Laan SW, Foroughi Asl H, van den Borne P, van Setten J, van der Perk MEM, van de Weg SM, et al. Variants in ALOX5, ALOX5AP and LTA4H are not associated with atherosclerotic plaque phenotypes: the athero-express genomics study. Atherosclerosis. (2015) 239:528–38. doi: 10.1016/j.atherosclerosis.2015.01.018
16. Adopted by the 18th WMA General Assembly, Helsinki, Finland, June 1964 amended by the 29th WMA General Assembly, Tokyo, Japan, October 1975, 35th WMA General Assembly, Venice, Italy, October 1983, 41st WMA General Assembly HK, September (1989), 48th WMA General Assembly, Somerset West, the 52nd WMA General Assembly, Edinburgh, Scotland O 2000, 53rd WMA General Assembly, Washington DC, USA, October 2002, et al. World medical association declaration of Helsinki. Bull World Health Organ. (2001) 79:373–4. doi: 10.1001/jama.2013.281053
17. de Bakker M, Timmerman N, van Koeverden ID, de Kleijn DPV, de Borst GJ, Pasterkamp G, et al. The age- and sex-specific composition of atherosclerotic plaques in vascular surgery patients. Atherosclerosis. (2020) 310:1–10. doi: 10.1016/j.atherosclerosis.2020.07.016
18. Hashimshony T, Senderovich N, Avital G, Klochendler A, de Leeuw Y, Anavy L, et al. CEL-Seq2: sensitive highly-multiplexed single-cell RNA-Seq. Genome Biol. (2016) 17:77. doi: 10.1186/s13059-016-0938-8
19. Hashimshony T, Wagner F, Sher N, Yanai I. CEL-seq: single-cell RNA-seq by multiplexed linear amplification. Cell Rep. (2012) 2:666–73. doi: 10.1016/j.celrep.2012.08.003
20. Ferraz MAMM, Rho HS, Hemerich D, Henning HHW, van Tol HTA, Hölker M, et al. An oviduct-on-a-chip provides an enhanced in vitro environment for zygote genome reprogramming. Nat Commun. (2018) 9:4934. doi: 10.1038/s41467-018-07119-8
21. R Core Team 2020. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing (2020). Available online at: http://www.eea.europa.eu/data-and-maps/indicators/oxygen-consuming-substances-in-rivers/r-development-core-team-2006 (accessed December, 2019).
23. Morgan M, Obenchain V, Hester J, Pagès H. SummarizedExperiment: SummarizedExperiment container. R package version 1.20.0. (2021). Available online at: https://bioconductor.org/packages/SummarizedExperiment
24. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. (2014) 15:550. doi: 10.1186/s13059-014-0550-8
25. Kolde R. Package ‘Pheatmap’: Pretty Heat Map. (2019). p. 1–8.
26. Ye Y, Li X, Zhang Y, Shen Z, Yang J. Androgen modulates functions of endothelial progenitor cells through activated Egr1 signaling. Stem Cells Int. (2016) 2016:7057894. doi: 10.1155/2016/7057894
27. Rouillard AD, Gundersen GW, Fernandez NF, Wang Z, Monteiro CD, McDermott MG, et al. The harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins. Database. (2016) 2016:baw100. doi: 10.1093/database/baw100
28. Hartman RJG, Korporaal SJA, Mokry M, de Jager SCA, Meeuwsen JAL, van der Laan SW, et al. Platelet RNA modules point to coronary calcification in asymptomatic women with former preeclampsia. Atherosclerosis. (2019) 291:114–21. doi: 10.1016/j.atherosclerosis.2019.10.009
29. Hartman RJG, Owsiany K, Ma L, Koplev S, Hao K, Slenders L, et al. Sex-stratified gene regulatory networks reveal female key driver genes of atherosclerosis involved in smooth muscle cell phenotype switching. Circulation. (2021) 143:713–26. doi: 10.1161/CIRCULATIONAHA.120.051231
30. Holm S. A simple sequentially rejective multiple test procedure. Scand J Stat. (1979) 6:65–70.
31. Whiteley W, Jackson C, Lewis S, Lowe G, Rumley A, Sandercock P, et al. Association of circulating inflammatory markers with recurrent vascular events after stroke: a prospective cohort study. Stroke. (2011) 42:10–6. doi: 10.1161/STROKEAHA.110.588954
32. Whiteley W, Jackson C, Lewis S, Lowe G, Rumley A, Sandercock P, et al. Inflammatory markers and poor outcome after stroke: a prospective cohort study and systematic review of interleukin-6. PLoS Med. (2009) 6:e1000145. doi: 10.1371/journal.pmed.1000145
33. Zhang XG, Xue J, Yang WH, Xu XS, Sun HX, Hu L, et al. Inflammatory markers as independent predictors for stroke outcomes. Brain Behav. (2021) 11:1–9. doi: 10.1002/brb3.1922
34. Mohamad NV, Wong SK, Wan Hasan WN, Jolly JJ, Nur-Farhana MF, Ima-Nirwana S, et al. The relationship between circulating testosterone and inflammatory cytokines in men. Aging Male. (2019) 22:129–40. doi: 10.1080/13685538.2018.1482487
35. Zöller B, Ji J, Sundquist J, Sundquist K. High risk of venous thromboembolism in klinefelter syndrome. J Am Heart Assoc. (2016) 5:1–6. doi: 10.1161/JAHA.116.003567
36. Spranger M, Aspey BS, Harrison MJG. Sex difference in antithrombotic effect of aspirin. Stroke. (1989) 20:34–7. doi: 10.1161/01.STR.20.1.34
37. Mygind ND, Pena A, Mide Michelsen M, Ali Qayyum A, Frestad D, Emil Christensen T, et al. Myocardial first pass perfusion assessed by cardiac magnetic resonance and coronary microvascular dysfunction in women with angina and no obstructive coronary artery disease. Scand J Clin Lab Invest. (2019) 79:238–46. doi: 10.1080/00365513.2019.1587670
38. Gencer B, Bonomi M, Adorni MP, Sirtori CR, Mach F, Ruscica M. Cardiovascular risk and testosterone - from subclinical atherosclerosis to lipoprotein function to heart failure. Rev Endocr Metab Disord. (2021) 22:257–74. doi: 10.1007/s11154-021-09628-2
39. Basaria S, Coviello AD, Travison TG, Storer TW, Farwell WR, Jette AM, et al. Adverse events associated with testosterone administration. N Engl J Med. (2010) 363:109–22. doi: 10.1056/NEJMoa1000485
Keywords: testosterone, atherosclerosis, RNA-expression, transcriptome (RNA-seq), men
Citation: Groepenhoff F, Diez Benavente E, Boltjes A, Timmerman N, Waissi F, Hartman RJG, Onland-Moret NC, Pasterkamp G and Den Ruijter H (2021) Plasma Testosterone Levels and Atherosclerotic Plaque Gene Expression in Men With Advanced Atherosclerosis. Front. Cardiovasc. Med. 8:693351. doi: 10.3389/fcvm.2021.693351
Received: 10 April 2021; Accepted: 19 May 2021;
Published: 14 June 2021.
Edited by:Susana Novella, University of Valencia, Spain
Reviewed by:Daniel Perez-Cremades, Brigham and Women's Hospital and Harvard Medical School, United States
Nabila Bouatia-Naji, Institut National de la Santé et de la Recherche Médicale (INSERM), France
Copyright © 2021 Groepenhoff, Diez Benavente, Boltjes, Timmerman, Waissi, Hartman, Onland-Moret, Pasterkamp and Den Ruijter. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Floor Groepenhoff, firstname.lastname@example.org
†These authors have contributed equally to this work