1. Sex is associated with acidic mucin staining in the proximal colon of wild Mus musculus.
Following capture of animals and excision of proximal colonic tissue, Alcian blue combined with Periodic acid, Schiff’s (AB+PAS) staining was carried out. This allowed visualisation of the goblet cells, and examination of the acidic nature of mucins stored within secretory granules. Alcian blue is a positively charged dye molecule which binds to and stains negatively charged sialic acid and sulphate residues on mucin glycans. PAS stains neutral glycan structures, but can also stain non-acetylated sialic acids in mucin glycans (Table 1). Importantly for this study, AB-positive but PAS-negative staining goblet cells contain predominantly acidic glycan structures (Figure 1A).
The proportion of PAS-negative staining goblet cells within the colons of 268 wild mice was highly variable, ranging from 0% to 77%, while 31.5% of mice (48 individuals) possessed no PAS-negative goblet cells (Figure 1B). Following this, model averaging using general linear models was undertaken to determine which environmental parameters were most associated with an animal’s possession of a population of PAS-negative goblet cells within the colon (parameters fully listed in Suppl. Figure 1A). Sex was reported to be the most significant predictor of PAS-negative goblet cells in the wild mouse population (full model output in Suppl. Figure 1B).
When comparing the cohort as a whole, the median percentage of goblet cells staining PAS-negative was similar between males, at 19.1% (inter quartile range (IQR): 34.4), and females, at 18.0% (IQR: 22.5) (Figure 1C). However, when the whole cohort was stratified by whether or not an individual possessed any PAS-negative goblet cells, a significantly greater proportion of female mice were seen to possess PAS-negative colonic goblet cells than male mice, a finding consistent across both 2018 (Figure 1H) and 2019 (Figure 1I). This suggests acidic glycan modifications are more widespread in the colonic goblet cells of wild female mice compared to male mice.
To further understand the differences in mucin-glycosylation between male and female mice, High-Iron diamine (HID) staining was performed to determine if the variation in the percentage of PAS-negative goblet cells was due to alterations in the sulphate or sialic acid content of mucin glycans. The HID procedure allows the visualisation of sulphate structures specifically, and when counterstained with Alcian blue, indicates sulphomucin-containing goblet cells (staining black), and goblet cells containing exclusively sialomucins (light blue staining) (Figure 1D).
The percentage of HID-positive goblet cells/crypt varied between individual mice in the wild population, with the median percentage of HID-positive goblet cells per crypt being 22.6% (IQR: 23.8) (Figure 1E). The proportion of PAS-negative and HID-positive goblet cells correlated significantly, but with a weak coefficient (R = 0.25) (Figure 1G). This suggests that although increased mucin sulphation may contribute to an increased prevalence of PAS-negative goblet cells, it may not be the major determinant. From this we conclude that increased sialylation of mucin glycans is the major change in glycosylation giving rise to increased population of PAS-negative staining goblet cells.
2. In laboratory C57BL/6 mice, sulphation of colonic mucins is associated with the estrus cycle of female mice, but colonic mucin sialylation underpins the difference between males and females.
Following the observation that sex associates with the proportion of sialomucin-containing goblet cells in wild mice, this was validated in laboratory-housed C57BL/6 mice. Goblet cells in the proximal colons of naive male and female C57BL/6 mice were examined using AB+PAS and HID staining. Differences between male and female laboratory mice were more pronounced than among wild mice, though variation within both sexes was reduced. Male mice presented with no PAS-negative stained goblet cells, whereas in female mice a median of 36.9% (IQR: 10.5) of goblet cells stained PAS-negative (Figure 2A). The extent of mucin sulphation as assessed by HID staining appeared similar among male (median = 3.4% ) and female mice (median = 4.6%) (Figure 2B). This supports our observation from wild mice that mucin sialylation underpins the differences in PAS-negative goblet cell incidence observed between male and female mice.
Next, we explored the role of systemic sex hormone signalling in relation to the observed sex-dependent differences of the colonic goblet cell staining among wild and laboratory mice. Female C57BL/6 mice were staged via histology of vaginal tissue [18] (Suppl. Figure 2), and after staging colonic tissue histology performed to characterise the acidic nature of the mucins within goblet cells. Quantifying the percentage of PAS- negative goblet cells did not reveal any differences between estrus cycles stages (Figure 2C). However, HID staining revealed a cyclical relationship between the extent of mucin sulphation within colonic goblet cells with respect to the estrus cycle stage of female mice, ranging from a median of 3.3% in the estrus stage, to 11.2% in the diestrus stage (Figure 2D).
3. An estrogenic diet is able to induce changes in the sulphation of goblet cell mucins.
In the wild mouse study, sunflower seeds were used to bait traps, resulting in a short (<12hr) availability of a novel, phytoestrogen-rich diet. Sunflower seeds contain phytoestrogens which can influence sex hormone signalling in mammals [19]. In order to assess the influence of sunflower seeds upon colonic mucin glycosylation, male C57BL/6 were fed on sunflower seeds for a period of 12 hours, 36 hours or 2 weeks.
No change in the percentage of PAS-negative goblet cells per crypt was seen at 12 hours, although the median percentage of HID-positive goblet cells increased from 1.5% to 14.5% at this timepoint. Further changes in mucin glycosylation were apparent by 2 weeks, with male mice displaying populations of PAS-negative goblet cells as high as 20% in some cases (Figure 2E). Additionally, an increase in mucin sulphation was observed with 100% of goblet cells containing HID-positive mucins in 3 of the 4 mice at the 2 week time point (Figure 2F).
4. Esr1 expression associates with the expression of sialyltransferase enzymes.
Transcriptomic analysis was carried out to identify what differential signalling events may be occurring in mice with varying degrees of mucin sialylation in the proximal colon. Twenty-four wild mice were selected for bulk RNA-sequencing analysis of colonic tissue, based upon the results of AB+PAS staining in the colon, which reflects changes in sialic acid and sulphate species within the mucin glycans. Mice were selected with the aim of encompassing a normal distribution of staining ‘profiles’ (Figure 3A). Sex, month and location were all controlled for within the 24 mice. Using IPA pathway analysis of the top 1000 most variable genes with respect to colour score, an index of AB+PAS staining. Esr1 activation was identified as the strongest predicted regulator of changes in gene expression (Figure 3B).
RNA sequencing confirmed that Esr1 was expressed in the colonic tissue of mice, while Esr2 was not, as observed in previous publications [19], [20] (Figure 3C). The expression of other sex-hormone receptors (Esr2, Pgr, Ar) was examined in the RNA sequencing dataset, and were not detected in proximal colonic tissue. Given that Esr1 was linked to PAS-negative goblet cell staining via RNA sequencing analysis, qPCR analysis of proximal colon tissue in the 2019 cohort was conducted to determine if mice possessing greater relative proportions of colonic sialomucins had increased expression of Esr1 (Suppl. Figure 3). However, levels of Esr1 expression had no association with PAS-negative goblet cell staining (Figure 3D). Female wild mice also exhibited no significant difference in the relative expression of Esr1 compared to male wild mice (Figure 3E).
Sialyltransferase enzymes catalyse the transfer of sialic acid from CMP-sialic acid to glycan structures [21]. In the context of mucin O-glycans, common substrates for the sialyltransferase enzymes are N-acetylgalactosamine (typically favoured by the St6GalNAc family), and galactose (favoured by the St3Gal family) [22]. The sialyltransferase genes St3gal1, St3gal2, St3gal3, St6galNAc1, St6galNAc2 and St6galNAc4 are all known to be expressed in the mouse colon, therefore the expression of those sialyltransferases was examined via qPCR. Interestingly, all sialyltransferase enzymes examined demonstrated a strong positive association with the expression of Esr1 (Suppl. Figure 4).
Finally, given the multi-variate environment within which the wild mice live and the many host-intrinsic and extrinsic variables at play, we asked whether certain variables could explain the variation in sialyltransferase enzyme expression via a redundancy analysis. The choice of variables used was informed by the ecological factors found to be most influential on goblet cell staining from general linear models (Suppl. Figure 1F) as well as Esr 1 expression, diet and age. Esr1 expression was the most highly significant factor (p = 0.001) in describing changes in the expression of sialyltransferases and associated with a general increase in the expression of all the sialyltransferase enzymes quantified (Figure 4, full output in Suppl. Figure 4). The age of mice was also approaching significance (p=0.07) (Suppl. Figure 4 D), and was particularly associated with the expression of St6galNAc1 (Figure 4).