Skip to main content

X chromosome inactivation skewing is common in advanced carotid atherosclerotic lesions in females and predicts secondary peripheral artery events

Abstract

Background and aim

Sex differences in atherosclerosis have been described with female plaques being mostly perceived as stable and fibrous. Sex-specific mechanisms such as mosaic loss of the Y chromosome in men have been linked to cardiovascular health. In women, X-linked mechanisms such as X chromosome inactivation (XCI) skewing is common in several tissues. Yet, information on the role of XCI in female atherosclerotic plaques is lacking. Here, we investigated the presence of XCI skewing in advanced atherosclerotic lesions and its association with cardiovascular risk factors, histological plaque data, and clinical data.

Methods

XCI skewing was quantified in 154 atherosclerotic plaque and 55 blood DNA samples of women included in the Athero-Express study. The skewing status was determined performing the HUMARA assay. Then, we studied the relationship of XCI skewing in female plaque and cardiovascular risk factors using regression models. In addition, we studied if plaque XCI predicted plaque composition, and adverse events during 3-years follow-up using Cox proportional hazard models.

Results

XCI skewing was detected in 76 of 154 (49.4%) plaques and in 27 of 55 (67%) blood samples. None of the clinical risk factors were associated with plaque skewing. Plaque skewing was more often detected in plaques with a plaque hemorrhage (OR [95% CI]: 1.44 [1.06–1.98], P = 0.02). Moreover, skewed plaques were not associated with a higher incidence of composite and major events but were specifically associated with peripheral artery events during a 3-year follow-up period in a multivariate model (HR [95%CI]: 1.46 [1.09–1.97]; P = 0.007).

Conclusions

XCI skewing is common in carotid plaques of females and is predictive for the occurrence of peripheral artery events within 3 years after carotid endarterectomy.

Plain language summary

Sex-differences have been observed in the development of atherosclerosis between men and women. Women tend to have more stable and fibrous plaques compared to men. Sex-specific mechanisms such as mosaic loss of the Y chromosome in men, were associated with cardiovascular health. In women, despite X-linked mechanisms like X chromosome inactivation (XCI) skewing was identified in various tissues. However, its relationship with atherosclerosis has not yet been investigated. In our study, we explored if prevalence of XCI skewing in advanced atherosclerotic lesions related to cardiovascular risk factors, histological plaque data, and clinical information. We found that XCI skewing was present in approximately 50% of human plaques, particularly those with plaque hemorrhage. Interestingly, we did not find any notable relationship between plaque skewing and clinical risk factors. However, we found that XCI was more present in women with peripheral artery events during the 3 years period following carotid endarterectomy. In summary, our findings indicate that XCI skewing is commonly observed in carotid plaques among females and may serve as a predictive factor for the occurrence of peripheral artery events within 3 years after carotid endarterectomy.

Highlights

  • XCI skewing is common in half of the carotid plaques of females and mostly in those plaques with a plaque hemorrhage.

  • None of the clinical risk factors were associated with plaque skewing.

  • Skewed plaques were specifically associated with peripheral artery events during a 3-year follow-up period after carotid endarterectomy.

Introduction

Atherosclerosis is often the underlying pathology for cardiovascular disease, characterized as a chronic and progressive inflammatory condition in which oxidized low-density lipoproteins infiltrate into the vascular wall, followed by inflammatory cells and leading to the formation of atherosclerotic plaques [1]. These plaques present significant phenotypic differences between the sexes [2]. Females have a more fibrous phenotype, more prone to erode while in men a vulnerable rupture-prone phenotype is more often reported [3, 4]. The differences in the clinical presentation of atherosclerosis between males and females are mainly attributed to sex hormones and sex chromosomes [5].

Females have two X chromosomes, one active (Xa) and one inactive (Xi) to ensure that the expression levels of X-linked genes are similar between females and males [6]. This dosage compensation mechanism is known as X-chromosome inactivation (XCI). Briefly, XCI is a process regulated by the expression of XIST on the Xi [7], and it is thought to start as early as the eight-cell stage during human embryonic development [8]. The inactivation of one of the parental X chromosomes seems to randomly occur within each cell during preimplantation phases and it gets mitotically inherited by all the future somatic daughter cells [6, 9]. This random inactivation results in a mosaic of cells within female tissues, where both parental X-linked alleles are expected to have balanced expression (50:50). The inactivation of a particular parental X is translated with an unbalanced expression and known as XCI skewing or nonrandom XCI [9].

The preferential X inactivation (XCI skewing) is the result of multiple mechanisms that can either take place directly during embryonic development (primary skewing) [10] or later in life (acquired or secondary skewing) due to positive selection of cells that after having inactivated a particular parental X, acquire a survival advantage over cells who inactivated the other parental X chromosome [9]. Twin studies have reported that genetic factors contribute to XCI skewing (primary) in blood-derived cells [11, 12], while other studies indicated that most of the XCI skewing levels in human are acquired secondarily [13], associated with ageing [11, 14,15,16,17,18] and frequent in the healthy and diseased female tissues [19]. Besides, females affected by X-linked disease [20,21,22], autoimmune disorders [23, 24] and cancers [25, 26] commonly present XCI skewing.

Hence, given the lack of information on X-linked mechanisms in female atherosclerosis, we investigated the prevalence of XCI skewing in advanced atherosclerotic plaque and blood of female patients who underwent carotid endarterectomy (CEA). Then, we studied if age and cardiovascular risk factors were related to plaque skewing. We also studied if plaque skewing was related to histological plaque characteristics and with an increased risk of having secondary cardiovascular events during 3-years follow-up.

Materials and methods

Patient characteristics

The Athero-Express biobank study is an ongoing cohort study that includes atherosclerotic plaques and blood of patients undergoing either carotid or femoral endarterectomy in 2 large tertiary referral hospitals (University Medical Center Utrecht and St. Antonius Hospital Nieuwegein) in the Netherlands. The patients used for this study have been selected based on the availability of the DNA samples isolated from carotid plaque and matched blood of female donors. Clinical data were obtained from medical files and standardized questionnaires. Age was determined as age at surgery. Current smoking was determined as patient-reported smoking in the past year. Hypertension and hypercholesterolemia were self-reported. Diabetes mellitus was considered present in any of the following cases: use of insulin or oral glucose inhibitors, self-reported diabetes mellitus in the patient questionnaire, or diabetes mellitus extracted from the medical file. A history of coronary artery disease was considered present if the patient had experienced a myocardial infarction or underwent a percutaneous coronary intervention or coronary artery bypass grafting surgery. Peripheral arterial occlusive disease was considered present if the patient either presented with an ankle-brachial index < 0.7, claudication complaints, or underwent percutaneous or surgical intervention for peripheral arterial occlusive disease.

Follow-up was obtained by questionnaires sent to the patients by mail 1, 2, and 3 years post-operatively. Composite events include myocardial infarction (MI), hemorrhagic and ischemic stroke, coronary angioplasty, peripheral intervention, cardiovascular death, coronary bypass, leg amputation, sudden death, fatal aneurysm rupture, other cardiovascular death. Major adverse cardiovascular events (MACE) include MI, hemorrhagic and ischemic stroke, cardiovascular death, sudden death, fatal heart failure fatal aneurysm rupture, other cardiovascular death. Peripheral artery events include percutaneous transluminal angioplasty (PTA), peripheral (re-) intervention of the limb arteries, and leg amputation [27]. All the endpoints were validated using medical records. The medical ethics boards of both hospitals approved of the study, which is conducted in accordance with the declaration of Helsinki, and the subjects gave informed consent.

Sample collection

Blood was obtained before surgery and plasma was subsequently stored at − 80 °C. Plaque specimens were immediately processed after removal during surgery. After identification of the area with the largest plaque burden (culprit lesion), the plaque was cut transversely into segments of 5 mm. The culprit lesion was fixed in 4% formaldehyde and subsequently decalcified and embedded in paraffin. Cross-sections were stained for histological examination. Remaining segments were stored at − 80 °C and used for the measurement of other paraments (e.g. inflammatory cytokines) and isolation of DNA. A detailed description of the samples phenotyping within the Athero-Express study can be found elsewhere [28].

HUMARA assay

The HUMARA assay was performed on the DNA isolated from blood and plaque as following. Working DNA solutions were prepared by diluting 10 × with DNAse free H2O and added, in an equal volume, to a 2 × enzyme reaction mix (2 × NH4 Taq buffer, 3u HpaII (NEB) and 4 mM MgCl2) and incubated at 37 °C of 1 h. then, an equal volume of PCR reaction mix (1 × NH4 Taq buffer, 1.5 mM MgCl2, 0.8 mM dNTP, 1uM Humara primer mix (fw = 5′-ACCGAGGAGCTTTCC AGAAT-3′; rv = 5′-TGGGGAGAACCATCCTCAC-3′) [29] and 1uM GAPDH primer mix (fw = 5′-CGCAGGCCGGATGTGTTC-3′; rv = 5′-ACACACACGCCTCCCCTC-3′), 0.25u Hotstart Taq polymerase (NEB) and 5% DMSO), was added to the restriction mix and a standard PCR program with a 60 °C annealing temperature and 13 cycles was performed. The PCR mix was diluted 2.5 × with MQ.

A second PCR reaction performed by adding 5ul of the diluted PCR mix to 15ul of reaction mix (1 × NH4 Taq buffer, 3 mM MgCl2, 0.8 mM dNTP, 1 µM HUMARA primer mix and 1 µM GAPDH primer mix, 0.25u Hotstart Taq polymerase and 5% DMSO).

A standard PCR program with a 60 °C annealing temperature and 15 cycles was performed. Samples were analyzed on 8% PAGE gel running at 200 V for 100 min and visualized with sybr safe.

Quantification of XCI skewing

The quantification of the XCI skewing levels has been assessed by measuring the band intensities after HUMARA assay. Interpretation of PCR results from each case required the following criteria: (1) successful PCR pre-amplification and amplification of DNA, both undigested and HpaII-digested; and (2) presence of two different androgen receptor alleles (heterozygous). An internal homozygous control, GAPDH containing HpaII digestion sites to exclude experimental artifact of incomplete enzyme digestion was included in each reaction. If GAPDH control showed undigested DNA, the reaction was considered as failed, repeated and if the similar result was obtained, the sample was removed from the analysis. Each reaction gave two intensity percentages (one per allele), the one equal or greater than 50% was selected and used in our data to determine presence of skewed XCI.

Determination of plaque skewing

Dichotomous plaque skewing variable was determined by distribution method using bestNormalize package in R Studio V. 1.1.456. The cut-off used to determine skewed plaque was if greater than 63.9% (Additional file 1: Fig. S1).

Binned plaque skewing variable has been defined by four skewing levels. The plaques having skewing percentage < 60% were defined as non-skewed (No, 37%), ≥ 60% and < 70%, as lowly skewed (Low, 30%), ≥ 70% and < 80% as skewed (Mid, 18%), and ≥ 80% as highly skewed (High, 15%).

Statistical analyses

XCI plaque skewing variables were associated with baseline characteristics. To determine possible confounders χ2 (categorical variables), ANOVA (continuous variables), and Kruskal–Wallis (non-normal continuous variables) tests, where applicable, were used.

Histological plaque characteristic variables (fat content, amount of calcification, collagen, presence of plaque hemorrhage, macrophage, smooth muscle cell, neo-vessel and glycophorin content) were used to assess whether plaque skewing associates with any of those plaque characteristics. We performed Saphiro test on the continuous variables to determine their distribution. All variables having Saphiro test with P ≤ 0.05 were considered non-normal distributed. These non-normal distributed variables (e.g. Glycophorin content) were transformed using BestNormalize package in R Studio to achieve normal distribution and used as outcomes in linear regression models. The transformation was based on the order quantile normalization technique and it was achieved applying the following transformation: g(x) = Φ − 1(rank(x) − 1/2 / length(x)) where x refer to the original data. The continuous variables that cannot be normalized (e.g. neo-vessel content) were binarized using BestNormalize package in R Studio and used as dichotomous outcomes for logistic regression.

The plaque characteristics (e.g. plaque hemorrhage) associating with plaque skewing were used to assess whether classical cardiovascular risk factors influence the association of plaque skewing with plaque characteristics. The significance of the association in the logistic regression models was assessed using Wald test and a P ≤ 0.05 was considered significant.

Kaplan–Meier, for univariate curves, and cox proportional hazard regression, for adjusted curves for all covariates that associated with outcome, methods were used to investigate the association between plaque skewing and secondary cardiovascular events (composite, major adverse cardiovascular and peripheral artery events) during 3-year follow-up. Statistical significance was assessed using Wald test and a P ≤ 0.05. Global statistical significance of the model was assessed using Mantel–Cox (log-rank) test. All statistical analyses were performed using R statistical package v. 4.1.2 within R Studio v. 1.1.456.

Differential expression analysis

Read counts obtained from RNA sequencing data and DESeq2 package in R were used for differential gene expression analysis between groups. Genes meeting the following criteria were considered as differentially expressed genes (DEGs): a P-value < 10–6 and a Log2 fold change ≥ 0.5 (up-) and ≤ − 0.5 (down-regulated). The volcano plot was generated using the “EnhancedVolcano” package in R. The genes that exhibited significant up- and down-regulation relative to the first condition (e.g., First Condition vs Second Condition) are depicted as red dots on the right and left sides, respectively.

Results

Plaque XCI skewing is highly prevalent in female atherosclerotic patients

We selected all female atherosclerotic patients who underwent carotid endarterectomy (CEA) and had DNA stored for XCI determination (n = 154). Of these 154 women, a subset of 55 women had DNA stored from blood samples. Baseline characteristics are displayed in Table 1. In short, these women were 66 years old on average and had a median BMI of 26 at baseline. At the day of the hospitalization, 32% of the patients had a previous history of peripheral artery occlusive disease (PAOD), while 26% had a previous history of coronary artery disease (CAD). Most of the patients suffered from hypertension (88%) and/or hypercholesterolemia (70%) and were on antiplatelet (94%) and/or lipid-lowering (77%) medications. Also, 71 patients (47%) reported to smoke and 19% had diabetes.

Table 1 Baseline characteristics of patients with and without XCI skewing in plaque

XCI skewing was determined in 154 carotid plaque and 55 matched blood DNA samples by performing a HUMARA assay followed by PAGE gel quantification. We used these quantitative data to determine XCI skewing status using a computed cut-off value of 63.9% (Additional file 1: Fig. S1). XCI skewing occurred in 76 of 154 plaques (49.4%), while in blood, it was higher with 67.2% (37/55) (Fig. 1 A). Blood skewing was strongly correlated with plaque skewing (P = 0.0001) (Fig. 1B). Age was not associated with blood (10-year β = − 0.14; P = 0.89, Additional file 1: Fig. S2 A, B) nor with plaque skewing (10-year β = 0.54; P = 0.43, Additional file 1: Fig. S2 C, D). The baseline characteristics showed no significant difference between the populations with skewed and non-skewed plaque (Table 1) or blood (Additional file 1: Table S1). The patients presenting skewed plaques were on average 67 years old and had a median BMI of 25, while the ones presenting with non-skewed plaques were on average 66 years old and had a median BMI of 26. Also, 37% of the skewed group had a history of PAOD and 29% of CAD while in the non-skewed group, 28% had history PAOD and 23% of CAD. Smokers were reported to be 41% in the skewed and 53% in the non-skewed group and the presence of diabetes was 18% and 19%, respectively.

Fig. 1
figure 1

A Bar plot showing the occurrence of skewing in percentage in plaque (49.4%) and blood (67.2%). B Scatter plot showing strong association between plaque and blood skewing (P = 0.0001)

To assess whether the extent of skewing (i.e., how skewed a plaque is) was important in relation to baseline characteristics, we created a categorical binned variable where X-skewing was classified in 4 levels (No, Low, Mid, High). Baseline characteristics of these 4 levels are shown in Additional file 1: Table S2.

XCI skewing in different plaque phenotypes

To investigate whether XCI skewing in plaque was more common in a specific plaque phenotype, we assessed the relation between plaque skewing and histological plaque characteristics. For this, fat content (> 10% and > 40%), amount of calcification, amount of collagen, presence of plaque hemorrhage, macrophage, smooth muscle cell, neo-vessel and glycophorin content within the plaque (Table 2) was studied [28]. Univariate analysis showed a linear association between plaque skewing (Continuous variable) and plaque hemorrhage (OR: 1.44 [1.06–1.98]; P = 0.02), but not when using the dichotomous variable (Table 2). Data on the binned plaque skewing levels are showed in Additional file 1: Table S3. As plaque skewing showed a significant association with plaque hemorrhage, we examined the correlation between blood skewing and plaque hemorrhage using the XCI data from blood of 55 female patients. Our univariate analysis revealed that the confidence intervals were very wide, indicative of low power. Nevertheless, effect sizes were small suggesting that there is no linear association between blood skewing and plaque hemorrhage (OR: 1.10 [0.70–1.75], P = 0.67). Accordingly, the limited sample size renders our analysis underpowered, precluding the derivation of significant conclusions.

Table 2 Association of plaque skewing with plaque characteristics

In addition, in our pursuit to elucidate the potential mechanisms behind XCI skewing in plaques, we used RNA samples from 73 out of the 154 donors available from a previous study [30]. We conducted a differential gene expression analysis, to compare skewed and non-skewed plaques, correcting for plaque hemorrhage (Additional file 1: Fig. S3A). We listed the top 10 differentially expressed genes based on the nominal P-value (Additional file 1: Fig. S3B). However, upon adjusting for multiple testing, we found no statistically significant differences between the skewed and non-skewed plaques.

Risk factors do not influence the relation between plaque skewing and plaque hemorrhage

Although risk factors did not show any significant differences between skewed and non-skewed plaques, we studied whether they modulate the association between plaque skewing and plaque hemorrhage (OR [95% CI]: 1.44 [1.06–1.98]; P = 0.02). We used classical cardiovascular risk factors, such as age, body mass index (BMI), smoking, diabetes mellitus and glomerular filtration rate to correct the model for one risk factor at the time. We found that BMI slightly attenuated the association of plaque skewing with plaque hemorrhage (OR [95% CI]: 1.37 [0.99–1.94]; P = 0.07). Also, smoking marginally increased this association (OR [95% CI]: 1.50 [1.10–2.09]; P = 0.01) (Table 3). Additional file 1: Table S4 shows the data on the binned plaque skewing levels.

Table 3 Effect of classical cardiovascular risk factors on the association between plaque skewing and plaque hemorrhage

Association with cardiovascular endpoints

We investigated if XCI skewing in plaque was associated with the incidence of secondary cardiovascular events in a 3-year follow-up after endarterectomy. For this, we used survival analysis, Kaplan–Meier (Additional file 1: Fig. S4) and cox proportional hazard model, corrected for BMI and current smoking (Fig. 2), to assess the association of plaque skewing with composite (41 events, 27%), major adverse cardiovascular (MACE; 13 events, 9%) and peripheral artery events (28 events, 19%) during 3-years follow-up.

Fig. 2
figure 2

XCI plaque skewing (yes/no), cox proportional hazards models for: A composite event-free survival (P = 0.08); B major event-free survival (P = 0.5); C peripheral artery event-free survival (P = 0.02). Models adjusted for BMI and current smoking. The survival probability is predicted until the occurrence of the last event

Women with skewed plaques had no more MACE (Fig. 2B) compared to women without skewing (HR [95%CI]: 0.89 [0.30–2.65]; P = 0.84), also when correcting for BMI and current smoking (HR [95%CI]: 0.80 [0.25–2.53]; P = 0.70) (Table 4). During 3 years of follow-up, plaque skewing was also not significantly associated with composite endpoints (Table 4; Fig. 2 A). However, a trend was observed in the association of dichotomous plaque skewing with composite events (HR [95%CI]: 1.70 [0.91–3.19]; P = 0.09) which increased slightly in strength when correcting the model for BMI and current smoking (HR [95%CI]: 1.82 [0.96–3.45]; P = 0.07). We found that women with skewed plaques had significantly more peripheral artery events during 3 years follow up (Table 4; Fig. 2C) when using both continuous (HR [95%CI]: 3.14 [1.38–7.18]; P = 0.01) and dichotomous XCI (HR [95%CI]: 1.46 [1.09–1.97]; P = 0.007) in multivariate analyses. Also adjusting for plaque hemorrhage did not change the relation between XCI and outcome (data not shown). We performed similar analysis using the binned plaque skewing levels which are showed in Additional file 1: Table S5, Fig. S5 and Fig. S6.

Table 4 Association of plaque skewing with secondary cardiovascular endpoints during 3-years follow-up

Discussion

In our study, we investigated XCI status of atherosclerotic plaques by examining the DNA methylation status of the polymorphic AR locus (HUMARA assay) and studied the occurrence of XCI skewing in carotid plaques on the DNA of 154 female atherosclerotic patients and its association with cardiovascular risk factors, histological and clinical outcomes.

We report that XCI skewing in plaques was present in half (49.6%) of the population. XCI plaque skewing was not associated with cardiovascular risk factors, but it was strongly associated with the presence of plaque hemorrhage. Plaque hemorrhage has been widely studied in the atherosclerosis field, and known to be a common phenomenon in advanced atherosclerotic plaques [4, 31]. In carotid lesions, it was found associated with an increased risk of secondary cardiovascular events in men but not in women [4]. Also, in our study, adjusting for plaque hemorrhage in the multivariate analysis did not change the relation between XCI and incident peripheral artery events. This may be explained by different mechanisms underlying symptomatic plaques in women and men [32]. In female plaques, cell plasticity seems to be dominant pointing to plaque erosion as a symptomatic mechanism in women as compared to plaque rupture in male plaques. We studied the plaque skewing and plaque hemorrhage further and hypothesized that infiltration of erythrocytes from either the bloodstream or leaky neo-vessels within the plaque [31] might be the driver of the association between plaque skewing and plaque hemorrhage in these female plaques. However, we did not find an association with glycophorin which is an erythrocyte-specific protein used to identify plaque bleeding [33]. A possible explanation for the lack of this association may be found in the quantification method. Indeed, plaque hemorrhage includes an expert assessment of the bleeding in the plaque while glycophorin staining has been measured via a computed process that cannot be subjected to the user experience [33, 34]. Hence, the associations with plaque hemorrhage or glycophorin data can be considered independent because none of them fully represents bleeding in atherosclerotic plaques. We replicated the same association analysis using the binned plaque skewing variable and we confirmed the association between plaque skewing and plaque hemorrhage also finding a dose-depended relationship.

We also found that age was not associated with plaque skewing although several studies showed that age was dose-dependently associated with skewing in general population [35,36,37]. This may be explained by the differences in the population studied, the tissue analyzed, and/or the assay used. We studied the XCI status of end-stage carotid plaques obtained from symptomatic women that have limited age distribution which may explain the lack of association. Also, the lack of healthy individuals in our study limits us to comment of the whether the prevalence of skewing in plaque is disease-related.

One of the possible mechanisms to explain XCI skewing in atherosclerotic lesions is the clonal expansion of cells involved in plaque formation [38]. Studies conducted in mice showed that macrophages undergo polyclonal expansion within atherosclerotic plaques, indicating that the proliferation of macrophages is already present in the plaque [39, 40]. However, in our study, we did not identify differences in the content of macrophages associated with skewed XCI. Potentially this could be due to the differences between mice and human plaques and/or the fact that the plaques in our study are in an advanced stage compared to plaques studied in mice.

Besides macrophages, the clonal expansion of smooth muscle cells was also proposed as a potential mechanism for XCI skewing in human carotid plaques, resulting in a population of monoclonal cells with the same XCI pattern [41, 42]. So, if the proposed hypothesis was true, a skewed XCI pattern would lead to a SMC-rich stable/fibrotic atherosclerotic plaque phenotype. However, this was not observed in our analysis arguing against clonal expansion being the driver of XCI. On the other hand, we observed a relation with peripheral artery event which are characterized by fibrous lesions, known to contain a large number of smooth muscle cells with enhanced proliferation capacities. This proliferative activity of smooth muscle cells seems to differ between vascular beds and in peripheral artery has a role in stabilizing of those lesions [42]. These SMC-mechanisms may be explained by clonal expansion which takes place in peripheral arteries, leading to events affecting limb arteries. In our population, we were unable to identify a higher content of smooth muscle cells associated with skewed XCI. One of the possible explanations for this may be due to the advanced disease stage on the patients in our study. It may be that the XCI skewing in carotid plaques plays a role on the cellular composition of plaques in an earlier phase of the disease and initially affects the plaque to become fibrotic. As this population is female only, and thereby the prevalence of fibrous plaques is already high, we may not find an association due to the lack of variation in smooth muscle cell content.

We also assessed the association between plaque skewing and secondary cardiovascular events to assess the clinical relevance of the presence of skewing. In our study we did not find an association of plaque skewing with either MACE or composite endpoints. Interestingly, in a secondary analysis we found that skewed plaques were associated with a higher occurrence of peripheral artery events during a 3-year follow-up period. Indeed, it is known that plaques found in the ilio-femoral arteries are often more fibrous than carotid plaques [43, 44], this in combination with our results may suggest XCI may be related to fibrotic processes in the periphery rather than the carotid arteries. However, peripheral artery endpoint was not the primary focus of the study and therefore this observation can only be seen as a sub-analysis which merits careful consideration and should be validated in an independent study.

This study had some limitations, the HUMARA assay may be influenced by preferential amplification of AR alleles with shorter repeats [45]. Also, the expression of a single X-linked locus may not reflect the expression status of the entire X-chromosome as there are genes with variable levels of escape from X-inactivation in the healthy population. Another limitation is the lack of healthy control tissues, which were previously reported to present extensive XCI skewing in a population of healthy women [19]. In our study, we did not examine whether XCI skewing is more prevalent in diseased arteries compared to normal ones. This might be challenging to study as atherosclerosis is present in the arteries of women without apparent arterial disease. Nevertheless, we observed a high occurrence of XCI skewing in blood samples (67% of the samples). This aligns with previous reports where the prevalence of XCI skewing is around 70% in individuals older than 55 years [16]. A possible explanation for the observed higher XCI skewing in blood compared to plaques may lie in the cellular composition of these tissues. Plaques consist of a complex mixture of cellular and acellular components, with a prominent presence of mesenchymal cells [46, 47]. Conversely, blood is characterized by a high abundance of leukocytes [48] that commonly undergo clonal proliferation, known to occur more frequently with advancing age [49]. This variation in cell composition likely contributes to the observed differences in XCI skewing prevalence between these tissue types. However, it is important to note that the mechanisms underlying XCI skewing are complex and multifactorial, and further research is needed to understand the relationship between tissue composition and XCI skewing.

Additionally, the major limitation of our study is the lack of power given by the limited number of patients included together with the fact that they only represent a snapshot of end-stage carotid plaques. However, this exploratory study highlights the need for further studies on the role of XCI skewing in atherosclerosis.

In conclusion, XCI skewing is common in carotid plaques of females and is predictive for the occurrence of peripheral artery events within 3 years after carotid endarterectomy.

Perspective and significance

Despite the well-known differences in the behavior of plaques between males and females, with female plaques being perceived as more stable and fibrous [2,3,4], the knowledge on XCI has been largely unexplored. Our study sheds light on the prevalence of X chromosome inactivation (XCI) skewing in female atherosclerotic plaques and its association with cardiovascular risk factors, histological plaque data, and clinical data. We found that XCI skewing was common in carotid plaques of females, with 49% of the plaques exhibiting skewing. None of the clinical risk factors were associated with plaque skewing. Furthermore, we found that skewed plaques were specifically associated with peripheral artery events within 3 years after carotid endarterectomy. These findings highlight the significance of understanding XCI skewing in female atherosclerotic plaques and its potential as a predictor of cardiovascular events. This study sheds light on the importance of sex-specific mechanisms in atherosclerosis.

Availability of data and materials

All data generated or analyzed during this study are included in this article.

Abbreviations

BMI:

Body mass index

CAD:

Coronary artery disease

CEA:

Carotid endarterectomy

HUMARA:

Human androgen receptor assay

MACE:

Major adverse cardiovascular events

PAOD:

Peripheral artery occlusive disease

SMC:

Smooth muscle cell

XCI:

X chromosome inactivation

References

  1. Libby P, Buring JE, Badimon L, et al. Atherosclerosis. Nat Rev Dis Primer. 2019;5(1):1–18. https://0-doi-org.brum.beds.ac.uk/10.1038/s41572-019-0106-z.

    Article  Google Scholar 

  2. Hartman RJG, Huisman SE, den Ruijter HM. Sex differences in cardiovascular epigenetics—a systematic review. Biol Sex Differ. 2018;9(1):19. https://0-doi-org.brum.beds.ac.uk/10.1186/s13293-018-0180-z.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Hellings WE, Pasterkamp G, Verhoeven BAN, et al. Gender-associated differences in plaque phenotype of patients undergoing carotid endarterectomy. J Vasc Surg. 2007;45(2):289–96. https://0-doi-org.brum.beds.ac.uk/10.1016/j.jvs.2006.09.051.

    Article  PubMed  Google Scholar 

  4. Vrijenhoek JEP, Den Ruijter HM, De Borst GJ, et al. Sex is associated with the presence of atherosclerotic plaque hemorrhage and modifies the relation between plaque hemorrhage and cardiovascular outcome. Stroke. 2013;44(12):3318–23. https://0-doi-org.brum.beds.ac.uk/10.1161/STROKEAHA.113.002633.

    Article  PubMed  Google Scholar 

  5. Den Ruijter HM, Haitjema S, Asselbergs FW, Pasterkamp G. Sex matters to the heart: a special issue dedicated to the impact of sex related differences of cardiovascular diseases. Atherosclerosis. 2015;241(1):205–7. https://0-doi-org.brum.beds.ac.uk/10.1016/j.atherosclerosis.2015.05.003.

    Article  CAS  Google Scholar 

  6. Lyon MF. Gene action in the X-chromosome of the mouse (Mus musculus L.). Nature. 1961;190:372–3. https://0-doi-org.brum.beds.ac.uk/10.1038/190372a0.

    Article  CAS  PubMed  Google Scholar 

  7. Brown CJ, Ballabio A, Rupert JL, et al. A gene from the region of the human X inactivation centre is expressed exclusively from the inactive X chromosome. Nature. 1991;349(6304):38–44. https://0-doi-org.brum.beds.ac.uk/10.1038/349038a0.

    Article  CAS  PubMed  Google Scholar 

  8. van den Berg IM, Laven JSE, Stevens M, et al. X chromosome inactivation is initiated in human preimplantation embryos. Am J Hum Genet. 2009;84(6):771–9. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ajhg.2009.05.003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Zito A, Davies MN, Tsai PC, et al. Heritability of skewed X-inactivation in female twins is tissue-specific and associated with age. Nat Commun. 2019;10(1):5339. https://0-doi-org.brum.beds.ac.uk/10.1038/s41467-019-13340-w.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Gribnau J, Barakat TS. X-chromosome inactivation and its implications for human disease. 2017:076950. https://0-doi-org.brum.beds.ac.uk/10.1101/076950.

  11. Kristiansen M, Knudsen GPS, Bathum L, et al. Twin study of genetic and aging effects on X chromosome inactivation. Eur J Hum Genet EJHG. 2005;13(5):599–606. https://0-doi-org.brum.beds.ac.uk/10.1038/sj.ejhg.5201398.

    Article  CAS  PubMed  Google Scholar 

  12. Vickers MA, McLeod E, Spector TD, Wilson IJ. Assessment of mechanism of acquired skewed X inactivation by analysis of twins. Blood. 2001;97(5):1274–81. https://0-doi-org.brum.beds.ac.uk/10.1182/blood.V97.5.1274.

    Article  CAS  PubMed  Google Scholar 

  13. Bolduc V, Chagnon P, Provost S, et al. No evidence that skewing of X chromosome inactivation patterns is transmitted to offspring in humans. J Clin Invest. 2008;118(1):333–41. https://0-doi-org.brum.beds.ac.uk/10.1172/JCI33166.

    Article  CAS  PubMed  Google Scholar 

  14. Christensen K, Kristiansen M, Hagen-Larsen H, et al. X-linked genetic factors regulate hematopoietic stem-cell kinetics in females. Blood. 2000;95(7):2449–51.

    Article  CAS  PubMed  Google Scholar 

  15. Wong CCY, Caspi A, Williams B, Houts R, Craig IW, Mill J. A Longitudinal twin study of skewed X chromosome-inactivation. PLoS ONE. 2011;6(3):e17873. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0017873.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Knudsen GPS, Pedersen J, Klingenberg O, Lygren I, Ørstavik KH. Increased skewing of X chromosome inactivation with age in both blood and buccal cells. Cytogenet Genome Res. 2007;116(1–2):24–8. https://0-doi-org.brum.beds.ac.uk/10.1159/000097414.

    Article  CAS  PubMed  Google Scholar 

  17. Hatakeyama C, Anderson CL, Beever CL, Peñaherrera MS, Brown CJ, Robinson WP. The dynamics of X-inactivation skewing as women age. Clin Genet. 2004;66(4):327–32. https://0-doi-org.brum.beds.ac.uk/10.1111/j.1399-0004.2004.00310.x.

    Article  CAS  PubMed  Google Scholar 

  18. Busque L, Mio R, Mattioli J, et al. Nonrandom X-inactivation patterns in normal females: lyonization ratios vary with age. Blood. 1996;88(1):59–65.

    Article  CAS  PubMed  Google Scholar 

  19. Shvetsova E, Sofronova A, Monajemi R, et al. Skewed X-inactivation is common in the general female population. Eur J Hum Genet. 2019;27(3):455–65. https://0-doi-org.brum.beds.ac.uk/10.1038/s41431-018-0291-3.

    Article  CAS  PubMed  Google Scholar 

  20. Migeon BR, Moser HW, Moser AB, Axelman J, Sillence D, Norum RA. Adrenoleukodystrophy: evidence for X linkage, inactivation, and selection favoring the mutant allele in heterozygous cells. Proc Natl Acad Sci U S A. 1981;78(8):5066–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Plenge RM, Stevenson RA, Lubs HA, Schwartz CE, Willard HF. Skewed X-chromosome inactivation is a common feature of X-linked mental retardation disorders. Am J Hum Genet. 2002;71(1):168–73. https://0-doi-org.brum.beds.ac.uk/10.1086/341123.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Puck JM, Nussbaum RL, Conley ME. Carrier detection in X-linked severe combined immunodeficiency based on patterns of X chromosome inactivation. J Clin Invest. 1987;79(5):1395–400. https://0-doi-org.brum.beds.ac.uk/10.1172/JCI112967.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Ozbalkan Z, Bagişlar S, Kiraz S, et al. Skewed X chromosome inactivation in blood cells of women with scleroderma. Arthritis Rheum. 2005;52(5):1564–70. https://0-doi-org.brum.beds.ac.uk/10.1002/art.21026.

    Article  CAS  PubMed  Google Scholar 

  24. Chabchoub G, Uz E, Maalej A, et al. Analysis of skewed X-chromosome inactivation in females with rheumatoid arthritis and autoimmune thyroid diseases. Arthritis Res Ther. 2009;11(4):R106. https://0-doi-org.brum.beds.ac.uk/10.1186/ar2759.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Kristiansen M, Langerød A, Knudsen GP, Weber BL, Børresen-Dale AL, Orstavik KH. High frequency of skewed X inactivation in young breast cancer patients. J Med Genet. 2002;39(1):30–3. https://0-doi-org.brum.beds.ac.uk/10.1136/jmg.39.1.30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Buller RE, Sood AK, Lallas T, Buekers T, Skilling JS. Association between nonrandom X-chromosome inactivation and BRCA1 mutation in germline DNA of patients with ovarian cancer. J Natl Cancer Inst. 1999;91(4):339–46. https://0-doi-org.brum.beds.ac.uk/10.1093/jnci/91.4.339.

    Article  CAS  PubMed  Google Scholar 

  27. van Haelst STW, Haitjema S, Derksen W, et al. Atherosclerotic plaque characteristics are not associated with future cardiovascular events in patients undergoing iliofemoral endarterectomy. J Vasc Surg. 2018;67(3):809-816.e1. https://0-doi-org.brum.beds.ac.uk/10.1016/j.jvs.2017.07.112.

    Article  PubMed  Google Scholar 

  28. Verhoeven BAN, Velema E, Schoneveld AH, 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(12):1127–33. https://0-doi-org.brum.beds.ac.uk/10.1007/s10564-004-2304-6.

    Article  CAS  PubMed  Google Scholar 

  29. Comertpay S, Pastorino S, Tanji M, et al. Evaluation of clonal origin of malignant mesothelioma. J Transl Med. 2014;12:301. https://0-doi-org.brum.beds.ac.uk/10.1186/s12967-014-0301-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Mokry M, Boltjes A, Slenders L, et al. Transcriptomic-based clustering of human atherosclerotic plaques identifies subgroups with different underlying biology and clinical presentation. Nat Cardiovasc Res. 2022;1(12):1140–55. https://0-doi-org.brum.beds.ac.uk/10.1038/s44161-022-00171-0.

    Article  Google Scholar 

  31. Michel JB, Virmani R, Arbustini E, Pasterkamp G. Intraplaque haemorrhages as the trigger of plaque vulnerability. Eur Heart J. 2011;32(16):1977–85. https://0-doi-org.brum.beds.ac.uk/10.1093/eurheartj/ehr054.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Hartman RJG, Owsiany K, Ma 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(7):713–26. https://0-doi-org.brum.beds.ac.uk/10.1161/CIRCULATIONAHA.120.051231.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Mekke JM, Sakkers TR, Verwer MC, et al. Glycophorin C in atherosclerotic plaque is associated with major adverse cardiovascular events after carotid endarterectomy. Published online July 19, 2021:2021.07.15.21260570. https://0-doi-org.brum.beds.ac.uk/10.1101/2021.07.15.21260570.

  34. Redgrave JNE, Lovett JK, Gallagher PJ, Rothwell PM. Histological assessment of 526 symptomatic carotid plaques in relation to the nature and timing of ischemic symptoms: the Oxford plaque study. Circulation. 2006;113(19):2320–8. https://0-doi-org.brum.beds.ac.uk/10.1161/CIRCULATIONAHA.105.589044.

    Article  CAS  PubMed  Google Scholar 

  35. Amos-Landgraf JM, Cottle A, Plenge RM, et al. X chromosome-inactivation patterns of 1,005 phenotypically unaffected females. Am J Hum Genet. 2006;79(3):493–9. https://0-doi-org.brum.beds.ac.uk/10.1086/507565.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Sharp A, Robinson D, Jacobs P. Age- and tissue-specific variation of X chromosome inactivation ratios in normal women. Hum Genet. 2000;107(4):343–9. https://0-doi-org.brum.beds.ac.uk/10.1007/s004390000382.

    Article  CAS  PubMed  Google Scholar 

  37. Sandovici I, Naumova AK, Leppert M, Linares Y, Sapienza C. A longitudinal study of X-inactivation ratio in human females. Hum Genet. 2004;115(5):387–92. https://0-doi-org.brum.beds.ac.uk/10.1007/s00439-004-1177-8.

    Article  PubMed  Google Scholar 

  38. Misra A, Rehan R, Lin A, Patel S, Fisher EA. Emerging concepts of vascular cell clonal expansion in atherosclerosis. Arterioscler Thromb Vasc Biol. 2022;42(3):e74–84. https://0-doi-org.brum.beds.ac.uk/10.1161/ATVBAHA.121.316093.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Robbins CS, Hilgendorf I, Weber GF, et al. Local proliferation dominates lesional macrophage accumulation in atherosclerosis. Nat Med. 2013;19(9):1166–72. https://0-doi-org.brum.beds.ac.uk/10.1038/nm.3258.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Härdtner C, Kornemann J, Krebs K, et al. Inhibition of macrophage proliferation dominates plaque regression in response to cholesterol lowering. Basic Res Cardiol. 2020;115(6):78. https://0-doi-org.brum.beds.ac.uk/10.1007/s00395-020-00838-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Murry CE, Gipaya CT, Bartosek T, Benditt EP, Schwartz SM. Monoclonality of smooth muscle cells in human atherosclerosis. Am J Pathol. 1997;151(3):697–705.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Espinosa-Diez C, Mandi V, Du M, Liu M, Gomez D. Smooth muscle cells in atherosclerosis: clones but not carbon copies. JVS-Vasc Sci. 2021;2:136–48. https://0-doi-org.brum.beds.ac.uk/10.1016/j.jvssci.2021.02.002.

    Article  PubMed  PubMed Central  Google Scholar 

  43. de Bakker M, Timmerman N, van Koeverden ID, et al. The age- and sex-specific composition of atherosclerotic plaques in vascular surgery patients. Atherosclerosis. 2020;310:1–10. https://0-doi-org.brum.beds.ac.uk/10.1016/j.atherosclerosis.2020.07.016.

    Article  CAS  PubMed  Google Scholar 

  44. Espitia O, Chatelais M, Steenman M, et al. Implication of molecular vascular smooth muscle cell heterogeneity among arterial beds in arterial calcification. PLoS ONE. 2018;13(1):e0191976. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0191976.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Swierczek SI, Piterkova L, Jelinek J, et al. Methylation of AR locus does not always reflect X chromosome inactivation state. Blood. 2012;119(13):e100–9. https://0-doi-org.brum.beds.ac.uk/10.1182/blood-2011-11-390351.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Hellings WE, Peeters W, Moll FL, et al. Composition of carotid atherosclerotic plaque is associated with cardiovascular outcome. Circulation. 2010;121(17):1941–50. https://0-doi-org.brum.beds.ac.uk/10.1161/CIRCULATIONAHA.109.887497.

    Article  PubMed  Google Scholar 

  47. Hellings WE, Moll FL, De Vries JPPM, et al. Atherosclerotic plaque composition and occurrence of restenosis after carotid endarterectomy. JAMA. 2008;299(5):547–54. https://0-doi-org.brum.beds.ac.uk/10.1001/jama.299.5.547.

    Article  CAS  PubMed  Google Scholar 

  48. Blumenreich MS. The white blood cell and differential count. In: Walker HK, Hall WD, Hurst JW, editors. Clinical methods: the history, physical, and laboratory examinations. 3rd ed. Butterworths; 1990.

    Google Scholar 

  49. Jaiswal S, Ebert BL. Clonal hematopoiesis in human aging and disease. Science. 2019;366(6465):eaan4673. https://0-doi-org.brum.beds.ac.uk/10.1126/science.aan4673.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We acknowledge the patients of Athero-Express biobank.

Funding

This study was funded by EU Project ERC Consolidator Grant 866478 (UCARE).

Author information

Authors and Affiliations

Authors

Contributions

MFB wrote the manuscript and performed the data analysis, MFB and MD performed HUMARA assay and quantification analysis, COM and EDB guided for the data analysis, BMM, JMM, GJdB, DPVdK and SWvdL provided datasets, performed clinical work and patients’ inclusions, HdR, MM, EDB and MB conceptualized the project, HdR, MM, EDB reviewed the manuscript.

Corresponding author

Correspondence to Hester M. den Ruijter.

Ethics declarations

Ethics approval and consent to participate

The Medical Ethical Committee of the UMC Utrecht has approved this study and all patients provided written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflicts of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Supplementary Data.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Buono, M.F., Benavente, E.D., Daniels, M. et al. X chromosome inactivation skewing is common in advanced carotid atherosclerotic lesions in females and predicts secondary peripheral artery events. Biol Sex Differ 14, 43 (2023). https://0-doi-org.brum.beds.ac.uk/10.1186/s13293-023-00527-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s13293-023-00527-6

Keywords