Skip to main content

An investigation of genetic polymorphisms in heparan sulfate proteoglycan core proteins and key modification enzymes in an Australian Caucasian multiple sclerosis population

Abstract

Multiple sclerosis (MS) is a chronic inflammatory demyelinating disease affecting the central nervous system in young adults. Heparan sulfate proteoglycans (HSPGs) are ubiquitous to the cell surface and the extracellular matrix. HSPG biosynthesis is a complex process involving enzymatic attachment of heparan sulfate (HS) chains to a core protein. HS side chains mediate specific ligand and growth factor interactions directing cellular processes including cell adhesion, migration and differentiation. Two main families of HSPGs exist, the syndecans (SDC1-4) and glypicans (GPC1-6). The SDCs are transmembrane proteins, while the GPC family are GPI linked to the cell surface. SDC1 has well-documented interactions with numerous signalling pathways. Genome-wide association studies (GWAS) have identified regions of the genome associated with MS including a region on chromosome 13 containing GPC5 and GPC6. International studies have revealed significant associations between this region and disease development. The exostosin-1 (EXT1) and sulfatase-1 (SULF1) are key enzymes contributing to the generation of HS chains. EXT1, with documented tumour suppressor properties, is involved in the initiation and polymerisation of the growing HS chain. SULF1 removes 6-O-sulfate groups from HS chains, affecting protein-ligand interactions and subsequent downstream signalling with HS modification potentially having significant effects on MS progression. In this study, we identified significant associations between single nucleotide polymorphisms in SDC1, GPC5 and GPC6 and MS in an Australian Caucasian case-control population. Further significant associations in these genes were identified when the population was stratified by sex and disease subtype. No association was found for EXT1 or SULF1.

Introduction

Multiple sclerosis (MS) is one of the most common neurological diseases affecting young adults in Western society [1]. It is an inflammatory, demyelinating, neurodegenerative disease affecting the central nervous system (CNS) destroying myelin, oligodendrocytes, axons and neurons [2, 3]. This disease is predominantly found in Caucasians with an onset between 18 and 40 years [4]. Onset of MS varies from early childhood to adult life, affecting more than 25,000 people in Australia in 2017 [5] and 2.3 million people worldwide in 2013 [6]. Symptoms include weakness in one or more limbs, visual disturbances and balance problems. As MS progresses, fatigue, bladder and bowel dysfunction, numbness, tremor, spasticity and/or depression may occur [7]. Individuals of Northern European ancestry, including those from Scandinavia, Iceland, the British Isles and North America, exhibit increased risk compared with African Americans (0.005%) [8, 9]. There is a female predominance of about 3:1 [10].

MS is a complex genetic disease characterised by polygenic heritable components and multifaceted gene-environment interactions and factors [11]. Several studies have investigated the interaction between genes and environmental factors [12] with underlying genetic factors implicated in determining familial clustering and individual susceptibility [13]. MS can be categorised into four subtypes: relapsing remitting MS (55%; RRMS); secondary progressive MS (30%; SPMS); primary progressive MS (10%; PPMS) and progressive relapsing MS (5%; PRMS) [9, 14]. The severity and frequency of attacks as well as the reversibility of damage suffered during these attacks varies with disease state [9, 14]. Females are twice as likely to develop RRMS as males [14, 15] and due to the progressive axonal loss, 50% of RRMS cases convert to the late phase SPMS within 8 to 20 years [9]. The conversion to more severe forms of disease results in progressive disability and increasing irreversible damage with fewer remitting stages and no symptomatic relief [7, 9, 13, 16,17,18].

Neurodegeneration is the major cause of permanent neurological disability in MS patients [13]. In the CNS, neurons are composed of an oligodendrocyte axon surrounded by myelin, a lipid-rich insulating sheath [19] that allows rapid conduction of nerve impulses in the form of an action potential [4]. Degradation of the myelin sheath, and hence axonal damage, resulting from deregulation of the immune system results in partial or complete blockage of CNS signal pathways [9, 13]. The formation of scar tissue (“sclerosis” or demyelinated plaques) in multiple locations within the CNS [9, 13] produces symptoms common to MS including sensory disturbances [9, 16]. The mechanisms by which damage occurs in different subtypes is thought to differ; however, the outcome remains the same—primary demyelination and secondary axonal degeneration [14].

Early genetic studies associated the major histone compatibility (MHC) complex with MS and identified alleles of the human leukocyte antigen (HLA) region on chromosome 6q21 (6q21-23) common to MS [20]. Between 15 and 60% of the genetic aetiology of MS has been attributed to this region [21]. A more recent strategy employed to identify genetic markers of MS is genome-wide association studies (GWAS). These studies examine single nucleotide polymorphisms (SNPs) across the genome and provide additional genetic information about common complex disorders including identification of polymorphisms associated with diseases such as cancer and MS ([22], 2011). In 2008, a GWAS identified non-HLA genes involved in MS [11]. Most of the genes identified were in immunological pathways (eg. interleukin-7 receptor [11, 23,24,25,26,27]; with confirmation of their impact gained through gene expression studies.

To date, the complex processes and factors that lead to demyelination, axonal damage, progressive neurological symptoms and the loss of immune homeostasis remain poorly understood. However, during nervous system development, the growing axons are thought to be guided towards their target by a complex that may include heparan sulfate proteoglycans (HSPGs) [28].

HSPGs are found in both vertebrates and invertebrates [29, 30] and are ubiquitous to the cell surface and the extracellular matrix (ECM). HSPGs belong to a large family of proteoglycans (PG) that are composed of one or more variably sulfated glycosaminoglycan chains (GAGs) attached to a core protein. Heparan sulfate (HS) chains undergo a complex, highly coordinated, sequential biosynthesis process resulting in a disaccharide chain consisting of N- and O-sulfated residues (Habuchi, [31]). This process involves a number of enzymes responsible for the initiation and subsequent modification of the growing GAG chain, resulting in a chain of variable length and sulfation pattern (Habuchi, [31]). Two of these enzymes are exostosin-1 (EXT1) and sulfatase-1 (SULF1). EXT1 encodes a glycotransferase protein that catalyses polymerisation of the GAG chain [32, 33]. Complete abrogation of EXT1 results in embryonic lethality during gastrulation due to the absence of heparan sulfate [34, 35].

HSPGs interact and bind a variety of growth factors (eg. FGF-2, TGF-β), chemokines and morphogenes [36] to mediate downstream signalling and subsequent cellular processes. The ligand-binding properties of HS are associated with chain length (initiated by EXT1 and others), but more particularly, specific sulfation patterns conferred through the action of enzymes such as SULF1. SULF1 removes 6-O-sulfate groups, altering binding sites for signalling molecules (eg. FGF2) [37] resulting in altered cell function. Mutations of SULF1 can result in a loss of function, potentially leading to lethality [38] while a double knockout of SULF1 results in neonatal lethality, similar to that seen in EXT1 knockouts [39] although probably through a different mechanism.

When bound to the basal membrane, HSPGs mediate cell-cell and cell-ECM interactions [40,41,42,43] with the HS chains acting as receptors or co-receptors triggering cell responses such as proliferation, adhesion, migration and signalling [44, 45]. Two major membrane bound families of HSPG core proteins are the four transmembrane syndecans (SDC; SDC1-4 [46];) and the six glycosyl-phosphatidylinositol (GPI) anchored glypicans (GPC; GPC1-6 [44, 47];). While SDC1 can carry other GAG chains, it predominantly carries HS chains [48] and is a known binding partner of FGF-2 [49]. GPCs carry only HS chains with attachment points close to the cell surface [50]. They are known to have a role in cell proliferation and differentiation and interact with the Wnt signalling pathway. In addition, GPCs have also been suggested to have a role in the inflammatory response [51].

One promising non-immune region, not previously identified in MS genome linkage screens, was 13q31 [52]. One GWAS analysed over 500,000 SNPs in 978 MS cases and 883 control samples and identified the HSPG GPC5 as a novel candidate gene for increased risk in MS [23]. A follow-up study genotyped an additional 974 MS affected individuals with similar results obtained for the GPC5 region [23]. A subsequent study investigating the GPC5 and GPC6 genes in susceptibility to MS determined their strongest association to be a specific GPC5 SNP (rs9523762) [51].

With the genetics of MS not completely understood, genes found to have a modest effect may provide potential new markers of therapy. Previous studies have identified GPCs in active MS plaques [53] and have associated these genes with disease. In contrast, little research has been conducted to elucidate the genetic involvement of SDCs, or HS chain modification enzymes, in MS development and progression. However, there is sufficient evidence to suggest involvement of these genes through their participation in cellular signalling. Indeed, FGF-2, a known signalling partner of the SDCs, has a key role in the regulation of differentiation and proliferation and therefore may have a role in remyelination [49].

To further investigate the role of HSPGs in MS susceptibility, this study set out to replicate previously identified GWAS SNP associations (GPC5, GPC6, 13q31-32 [52]) and to extend our knowledge of the potential involvement of a number of related genes in an Australian Caucasian case-control population. We investigated polymorphisms within HS chain modification enzymes, EXT1 and SULF1, as well as those within HSPG core proteins SDC1, GPC5 and GPC6. Disruption to the critical enzymes responsible for the diverse functions of the HS side chains, or the HSPG core proteins themselves, could have an adverse effect on the progression of MS.

Methods

Population DNA

This study utilised a Caucasian (Northern European descent) case-control population. The population consisted of 205 multiple sclerosis (MS) patient samples and 194 age (+/− 5 years) and sex matched control samples. The case samples were then further divided into three disease states: relapsing-remitting (RR), secondary progressive (SP) and primary progressive (PP) MS. The case group consisted of 160 females and 45 males. A summary of the numbers of males and females in each group can be found in Table 1. Genomic DNA was extracted from peripheral blood using a standard salting-out procedure described previously ([54]; Griffith Ethics Approval: 1300000488; QUT Ethics Approval: 1400000502).

Table 1 Population demographics

DNA was analysed for quantity and quality using the Nanodrop (Thermo Scientific, Australia). Where necessary samples were cleaned using standard ethanol precipitation protocols and reanalysed for quantity and quality. Samples were diluted to a working concentration of 20 ng/μL. Aliquots were stored at −20 °C.

Primers

Primers for EXT1, SULF1 and SDC1 were designed using NCBI primer blast (NCBI). Primers for GPC5 and GPC6 were designed using PyroMark Assay Design Software v2.0 (Qiagen, Australia) for use in pyrosequencing and subsequent analysis. For GPC6, rs9524260, a sequencing primer was included in the design process. Primer sequences can be found in Table 2. Primers were synthesised by Integrated DNA Technologies (IDT), USA.

Table 2 Assay details and SNP information including RFLP fragment sizes where appropriate

Polymerase chain reaction

For EXT1 and SULF1, 40 ng of DNA was amplified with 1× PCR buffer, 100 nM each forward and reverse primers (IDT, USA), 200 μM dNTPs (NEB, Australia), 1.75 mM MgCl2, 0.5 U GoTaq Flexi DNA Polymerase (Promega, Australia) in a 15-μL reaction. Optimal reaction conditions for SDC1 amplified 40 ng of DNA with 1× PCR buffer, 200 nM each forward and reverse primers, 200 μM dNTPs, 1.75 mM MgCl2, 0.5 U GoTaq Flexi DNA polymerase in a 15-μL reaction.

Annealing temperatures (Ta) for individual assays can be found in Table 2. Cycling conditions for these three SNPs were as follows: an initial denaturation step at 95 °C for 10 min was followed by 35 cycles of 95 °C for 30 s, annealing at Ta for 30 s and extension at 72 °C for 30 s. This was followed by a final extension step at 72 °C for 5 min.

PCR products for all SNPs were analysed on a 2% agarose gel in 1× TAE with ethidium bromide at 90 V for at least 30 min to confirm amplification of single PCR products of the correct size. A 100 bp ladder was used alongside the samples for sizing purposes. Fragments were visualised using UV light.

Restriction fragment length polymorphism

Genotypes for EXT1, SULF1 and SDC1 were determined by RFLP. EXT1 PCR product (7 μL, approximately 1 μg DNA) was digested with 0.5 U Cac8I enzyme, 1× NEB reaction buffer 4 in a total volume of 15 μL. Digestion occurred overnight at 37 °C followed by an inactivation step at 65 °C for 20 min. Digest products were then analysed on a 3% agarose gel in 1× TAE with ethidium bromide at 90 V for 45 min. DNA fragments were visualised under UV light. A 100 bp ladder was used alongside digest products for sizing purposes. The enzyme Cac8I recognises the sequence GCNNGC. The presence of the wild-type G allele allows for digestion producing fragments of 281 bp and 115 bp. The homozygous variant (AA) remained uncut with a single band at 397 bp.

SULF1 and SDC1 PCR products (7 μL, approximately 1 μg DNA) were digested with 1 U PspGI enzyme, 1× NEB reaction buffer 4 in a 15-μL reaction. Samples were digested at 75 °C for 16 h. No inactivation step was required. SULF1 digest products were analysed on a 3% agarose gel as for EXT1. SDC1 digest products were assessed on a 4% agarose gel in 1× TAE with ethidium bromide at 75 V for 1 h along with a 100 bp ladder for sizing purposes. Fragments were visualised under UV light. The enzyme PspGI recognises the sequence CCWGG. The digestion of SULF1 PCR products produces bands at 212 bp and 78 bp in the presence of the wild-type C allele while the homozygous variant remained uncut with a single band visible at 291 bp. When PCR product from SDC1 was digested, bands were produced at 188 bp, 37 bp and 17 bp regardless of the allele present. In the presence of the wild-type G allele, bands are also produced at 53 bp and 43 bp, while the homozygous variant (CC) contains a single extra band at 96 bp.

High resolution melt

High resolution melt (HRM) analyses were performed on the Rotor Gene™ 6000 (Corbett, Australia), Rotor-Gene® Q (Qiagen, Australia) or with the 7900 HT Fast Real-Time PCR System (Applied Biosystems, Australia). All DNA case-control samples were assayed in duplicate. For each SNP, positive control samples were identified from samples not included in the analysis population. Where possible, positive controls representing each of the three genotypes were identified for each SNP.

HRM reaction conditions amplified 40 ng of DNA with 1× reaction buffer, 1.75 mM MgCl2, 100 nm each of forward and reverse primers, 100 nM dNTPs, 50 μM Syto®9 and 0.5 U GoTaq Hot Start DNA polymerase in a 15-μL reaction. Cycling conditions on the ABI 7900 HT system included an initial 2-min hold at 50 °C followed by denaturation and HotStart polymerase activation step of 95 °C for 10 min. This was followed by 40 cycles of 95 °C for 15 s and Ta for 1 min. Finally, products were melted between 60 °C and 90 °C to produce the melt curves required to differentiate between genotypes. The Rotor Gene 6000 system utilised the following cycling conditions: denaturation and HotStart polymerase activation at 95 °C for 10 min followed by 40 cycles of 95 °C for 5 s, Ta for 10s, and a final melt between 70 °C and 90 °C. The ABI 7900 HT system was used to genotype rs7333912 and rs9523787 (GPC5) while the Rotor Gene 6000 system was used for genotyping rs10492503 (GPC5) and rs17267815 (GPC6).

Pyrosequencing

Pyrosequencing was performed as per the manufacturer’s instructions on the Pyromark Q24 (Qiagen, Australia). Briefly, 15 μL of PCR product was combined with 2 μL of Streptavidin Sepharose beads (GE Healthcare, Australia), 40 μL of binding buffer (10 mM Tris-HCl, 2 M NaCl, 1 mM EDTA, 0.1% Tween™ 20, pH 7.6) in a reaction volume of 80 μL. Amplicons were denatured (denaturation buffer: 0.2 M NaOH) and washed with washing buffer (10 mM Tris-Acetate) and 70% ethanol. The clean biotin-labelled amplicons were transferred to a sequencing plate containing 0.3 μM sequencing primer diluted in annealing buffer (20 mM Tris-Acetate, 5 mM MgAc2) and denatured for 2 min at 80 °C. Samples were loaded into the PyroMark Q24 chamber for analysis. The reagent cartridge was loaded with appropriate volumes of dNTPs and Enzyme and Substrate mix (Qiagen) as determined by proprietary Q24 equipment software.

Sequencing validation

To verify the results obtained from the HRM genotyping represent the three different genotypes, positive controls and examples of each genotype from the population samples were sequenced. This was performed by Sanger sequencing using the BigDye® Terminator (BDT) v3.1 Cycle sequencing kit (Thermo Fisher Scientific, Australia). Briefly, PCR product was cleaned with ExoSAP-IT® (Affimetrix) following manufacturers protocols. PCR product concentration was estimated and adjusted for optimal sequencing conditions. Following the BDT reaction, the samples were then cleaned using a standard ethanol precipitation method, dried and resuspended in water. Forward and reverse reactions for each sample were prepared. Samples were then loaded into a 96-well plate followed by separation on a four capillary 3130 genetic analyser system (Thermo Fisher Scientific, Australia). The results were later analysed with AB Sequencing Analysis Software v5.3.

Statistical analysis

Genotype and allele frequencies were calculated using a standard counting method. Populations were analysed by Hardy-Weinberg equilibrium (HWE) and chi-square tests. A significance level of p < 0.05 was used. Where allele or genotype analysis was significant, the population was stratified by either disease type or sex and reanalysed. Global p values were adjusted for multiple testing using the Benjimini-Hocherg and Bonferroni methods. Corrections were conducted in R v3.3.0 and Rstudio v0.99.896. Corrected p values are presented in results tables.

For the GPC5 and GPC6 SNPs within the chromosome 13 risk region, linkage disequilibrium (LD) analysis was preformed using Haploview v4.2.

Results

We examined eight single nucleotide polymorphisms (SNPs) within HSPG initiation and modification enzyme (EXT1 and SULF1) and HSGP core protein (SDC1, GPC5 and GPC6) genes in an Australian case-control population for their role in MS susceptibility. Genotype and allele frequencies were also compared to results from the HapMap CEU reference population. The chromosomal region of 13q31-32, where GPC5 and GPC6 are located, has been previously identified in a GWAS as showing a significant association with genetic susceptibility of multiple sclerosis (MS). We aimed to investigate the role of these genes as well as the potential role of enzymes responsible for the modification of the HS chains attached to these and other core proteins. This study investigated eight SNPs in five genes. Differences in final population numbers successfully genotyped for each SNP analysis are due to variation in assay efficiency as well as DNA availability.

Modification enzymes

EXT1 and SULF1 enzymes initiate and modify HS side chains attached to core proteins. The action of these enzymes determines the final length and sulfation pattern of the side chain and mediate signalling interactions (eg. growth factors). They are critical for HSPG synthesis and any mutation could result in an autosomal dominant disorder [33].

EXT1, rs11546829

No significant association was found between the rs11546829 SNP in EXT1 and MS. Both case and control populations followed HWE and allele and genotype frequencies matched the HapMap CEU reference population. When the population was stratified by disease type and further analysed, there was no association found with disease state. Results are summarised in Table 3.

Table 3 Genotype and allele frequencies of EXT-1 SNP (rs11546829) within the case-control MS population which is further subdivided into disease states. Corrected P values using the Benjimini-Hochberg (PBH) and Bonferroni (PBon) methods are presented below the uncorrected P value

SULF1, rs2623047

No significant association was found between rs2623047 in SULF1 and MS. HWE was observed in both case and control populations. Allele and genotype frequencies matched the HapMap CEU reference population. Further analysis of the stratified populations revealed no significant association with disease type. Results are summarised in Table 4.

Table 4 Genotype and allele frequencies in the case-control MS population for the SULF-1 SNP (rs262347) further subdivided into disease states. Corrected p values using the Benjimini-Hochberg (PBH) and Bonferroni (PBon) methods are presented below the uncorrected P value

Core proteins

Syndecan-1

SDC1 is known to have a role in various cancers, including breast cancer [56] through its role in cell adhesion, migration and proliferation. Through its interaction with FGF-2, SDC1 has also been proposed to have a role in remyelination [49]. Evidence suggests that SDC1 may be associated with MS through enhanced expression of TGF-β in MS lesions that may lead to increased expression of SDC1 [55]. A link between SDCs and the innate immune response has also been postulated as these HSPGs have been observed in injured tissues, regulating the accompanying inflammatory response [48], suggesting a link to the inflammatory response seen in MS patients.

SDC1, rs1131351

A positive association between the rs1131351 SNP in SDC1 and MS was identified. Both case and control populations followed HWE. When compared to the HapMap CEU, reference population allele and genotype frequencies were similar. Allele and genotype frequencies showed significant differences between the case and control populations. When the population was examined by disease state, significant association for the SDC1 SNP was observed at the allelic level for RRMS and PPMS for the whole population. There was no significant association with SPMS. Results are summarised in Table 5.

Table 5 Results by disease state for SDC1 SNP (rs1131351) and of the MS case-control population. Corrected P values using the Benjimini-Hochberg (PBH) and Bonferroni (PBon) methods are presented below the uncorrected P value

The population was further analysed by sex, where SDC1 showed a further significant association in the female MS population (Table 5). The female population was further stratified by disease state. A significant association was observed in females with PPMS at both the genotype and allelic level and female RRMS cases at the allelic level only with no observed significance with genotypes. The female SPMS group demonstrated no significant association at either allele or genotype level and MS. Results are summarised in Table 6.

Table 6 Female results by disease state for SDC1 SNP (rs1131351)

Glypican-5

The three GPC5 SNPs (rs7333912, rs10492503 and rs9523787) investigated in this study have previously been significantly associated with MS in Caucasian European populations [51, 57]. Another GPC5 SNP (rs9523762), not reported here, was found to be positive in one study [23] while another study identified moderate LD between it and rs9523787 but did not find it to be individually significant [51]. Indeed, in this current study, we also investigated this SNP but have not presented the results as analysis revealed significant deviation from HWE. As this Australian population consists of Caucasian ancestors, we performed an associative study with these SNPs in an Australian Caucasian population to see if the association could be replicated in this cohort. Differences between our results and previous studies may be due to the more mixed heritage of the Australian Caucasian population compared with the purer northern European Caucasian populations previously examined.

GPC5, rs7333912

The GPC5 SNP, rs7333912, is an intergenic SNP, located at 13q31-32 approximately 150,000 bp upstream of the GPC5 gene on chromosome 13. This variation is a C/G polymorphism with no homozygous variants (GG) observed either in our population or in the HapMap CEU reference population. The case and control populations were found to be in HWE with no significant association identified between this SNP and MS. When the population was stratified and analysed by disease state, no significant association was found with disease type. Results are summarised in Table 7.

Table 7 Results for GPC5, rs7333912. Corrected P values using the Benjimini-Hochberg (PBH) and Bonferroni (PBon) methods are presented below the uncorrected P value

GPC5, rs10492503

GPC5-rs10492503 is located in the same intron as rs9523787 (intron7-8) within the 13q31-32 MS susceptibility locus that has previously been associated with MS [52]. In our study, both case and control populations observed HWE. A significant association was identified between this SNP and disease at both the genotype and allele level. When further analysed by disease state, further significant associations were determined in the SPMS group at the allelic level and RRMS case group at the both the genotypic and allelic level. Stratification by sex determined a significant association between this SNP and the female case group at both genotype and allele level. Results are summarised in Table 8. The female case group was further analysed by disease state. A significant association was found at both genotypic and allelic level between the SNP and female SPMS and RRMS disease subtypes. These results are summarised in Table 9.

Table 8 Results for GPC5, rs10492503. Corrected P values using the Benjimini-Hochberg (PBH) and Bonferroni (PBon) methods are presented below the uncorrected P value
Table 9 Female results by disease state for GPC5, rs10492503. Significance for the PP case subgroup is suggestive only as cell counts fell below the minimum required for chi-squared testing (n < 5)

GPC5, rs9523787

This variation is located toward the 3’ end of the gene in intron 7-8, the same intron that contains rs10492503 [52]. GPC5-rs9523787 had previously been associated with MS and to be in modest LD with a variation located close by, rs9523762 [51]. We were unable to replicate either of these findings. Both case and control populations followed HWE and we found no significant association between this variation and the MS population. Populations stratified by disease type also showed no significant association. Results are summarised in Table 10.

Table 10 Results for GPC5, rs9523787. Corrected P values using the Benjimini-Hochberg (PBH) and Bonferroni (PBon) methods are presented below the uncorrected P value. Significance measures are suggestive only as cell counts fell below the minimum required to perform chi-square analysis in the disease subgroups (n < 5)

Glypican-6

GPC6 has previously been implicated in MS [51]. In that study, 22 SNPs were analysed with only three showing a significant association with disease (GPC5-rs7333912, GPC6-rs17267815, GPC6-rs12876985). Only GPC5-rs7333912 and GPC6-rs17267815 were included in this study. The second GPC6 SNP investigated in this study was previously reported to be associated with primary sclerosing cholanigitis (PSC) [58]. GPC6 is located in the chromosome region of 13q31 neighbouring the risk region identified in genome wide screens [52].

GPC6, rs17267815

Of 22 SNPs analysed by Lorentzen and colleagues, this SNP (GPC6-rs17267815) showed the greatest significance [51] with MS. In our study, no significant association was identified in the total case group versus controls. However, when the cases were analysed by disease state, we found a significant association between the RRMS case group and the SNP at both the genotype and allelic level. Results are summarised in Table 11. There was no significance when stratified by sex. However, when the population was stratified by both disease and sex, an association was found at both the genotype and allelic level for the male RRMS subgroup. This is, however, only suggestive as the sample numbers for the homozygous variant genotype fell below the minimum required for reliable chi-square testing (n < 5). These results are summarised in Additional file 1: Table S1.

Table 11 Results for GPC6, rs17267815. Corrected P values using the Benjimini-Hochberg (PBH) and Bonferroni (PBon) methods are presented below the uncorrected P value

GPC6, rs9524260

Located within the 13q31 risk region, rs9524260 is an intronic SNP. Both case and control populations were found to be in HWE with no significant association found between this variation and MS. There was also no significant association when the population was stratified by disease type or sex. The results are summarised in Table 12.

Table 12 Results for GPC6, rs9524260. Corrected P values using the Benjimini-Hochberg (PBH) and Bonferroni (PBon) methods are presented below the uncorrected P value

LD analysis of GPC5 and GPC6 SNPs

GPC5 and GPC6 markers were analysed for LD using Haploview v4.2. No LD was observed in this population. The highest D’ value determined was between rs10492503 and rs9523787 in GPC5 with D’ = 0.15. The LD plot can be seen in Fig. 1.

Fig. 1
figure 1

LD Plot from GPC5/GPC6 haplotype analysis. The figure shows there is no LD between these SNPs in an Australian Caucasian MS population. This analysis was unable to replicate the positive associations and LD found in previous studies

Discussion

Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system (CNS) with a presumed autoimmune origin, triggered by genetic and environmental risk factors. The aetiology of MS is unknown, and the pathology is not well understood. In addition to those SNPs identified as significant in a 2009 GWAS [23], we investigated two SNPs in enzymes responsible for the initiation and modification of the side chain characteristic of HSPGs (EXT1, SULF1) and another HSPG core protein, SDC1. No significant difference was observed in our MS population in relation to the EXT1 and SULF1 SNPs examined. Further analysis revealed no association with disease subtype. However, in this study, we did identify significant associations with SDC1, GPC5 and GPC6 polymorphisms. Overall, analysis by disease subtype maintained this significance, as did analysis by sex.

While this study may suggest that no association exists between EXT1 and MS, previous studies have indicated a strong expression of EXT1 in the developing brain [32]. Additionally, it has been suggested that EXT1 correlates with the sites of active neuron generation [32]. Prenatally EXT1 has been localised in the neuroepithelial cells, which surround the lateral ventricles, cerebral cortex and hippocampus. However, in the postnatal stage, EXT1 is expressed in the cerebellum, which may correlate with the symptoms seen in MS such as ataxia [32]. In a murine model, complete abrogation of EXT1 results in embryonic lethality during gastrulation due to the absence of HS [34]. EXT1 alone is able to polymerise GAG chains for attachment to PG core proteins [59, 60]; however, both EXT1 and EXT2 are required for in vivo HS chain elongation [61].

SULF1 has been suggested to have a role in the brain; however, deficiencies in SULF1 have been associated with developmental abnormalities such as decreased body mass and subtle kidney and bone defects [34]. SULF1 has also been linked to tumour suppressor functions as it has been reported to be downregulated in some cancers [37].

The process of HS chain biosynthesis requires the action of enzymes such as SULF1 to generate complex sulfation patterns through the addition and removal of sulfation sites. Successful binding of growth factors to GAG chains for signalling pathway activation requires specific sulfation patterns on these side chains. The interaction between HS and FGF-2 is a well-documented example where HS requires 2-O sulfation to be able to bind to FGF-2. Without this binding, cellular proliferation and differentiation are inhibited. Although these two SNPs in EXT1 and SULF1 showed no significance in our population of moderate size, further investigation should be undertaken with a larger cohort before excluding the possibility of their involvement in MS susceptibility.

In the first of the HSPG core proteins examined (SDC1), we found a significant association between the SNP, rs1131351 and MS. This association revealed a stronger link between the SNP and females suffering from early-onset forms of the disease (PPMS, RRMS). SPMS occurs 8–20 years after RRMS onset [14] with the negative association seen here with this disease state suggestive of a role for SDC1 in the initiation of disease. Females with PPMS and the minor allele of SDC1 have more than double the risk (OR = 2.24) of developing MS than controls. In patients suffering from RRMS, this increase in risk is approximately 1.5 times (OR = 1.57). This could be due to the fact that PPMS seems to be more aggressive during onset when compared with RRMS. Even though they are both classified as onset stages of the disease, RRMS can progress to SPMS, with reversible damage occurring in this stage, while PPMS damage is irreversible and the symptoms are generally more detrimental (reviewed in [62, 63].

Active MS lesions are characterised by an influx of inflammatory cells and a decrease of chondroitin sulfate proteoglycans (CSPGs) [55]. Furthermore, white matter-associated PGs have been known to accumulate in macrophages, suggesting that CSPGs are phagocytosed with myelin or their breakdown products [55]. SDC1 contains ser-gly sequences that may serve as an attachment site for chondroitin sulfate (CS) [48] while also carrying HS chains. A mutation in SDC1 may contribute to activation of the macrophages causing phagocytosis, subsequently leading to a reduction in SDC1 in MS patients. In addition, TGF-β along with FGF-2 have been linked to enhanced expression of SDC1 [48]. Enhanced expression of TGF-β has been observed in MS lesions causing matrix deposition by the promotion of transcription genes and suppression of degrading enzymes [55]. FGF-2 has been associated with the survival, proliferation and migration of oligodendrocyte precursors leading to the promotion of remyelination [64]. This contradicts the mechanism of neurodegeneration seen in MS patients; however, FGF-2 could be a survival mechanism established to reverse the damage particularly in relapsing and remitting MS patients, through its binding partners other than SDC1.

In addition to the SDC1-FGF-2/TGF-β signalling mechanisms, TNF-α has been demonstrated to decrease SDC1 expression in cultured endothelial cells [48]. TNF-α has been shown to be involved in the inflammatory response [65] and could be involved in the process mimicking the early stages of MS where the breakdown of the blood-brain barrier allows inflammatory cells to cross into the brain and contribute to demyelination and axonal damage [55].

In this study, we aimed to replicate and build on results from a number of previous GWAS and replication studies using an Australian case-control population. These earlier results implicated GPC5 and GPC6 SNPs in MS. Our analysis of three GPC5 and two GPC6 SNPs also identified significant associations between these genes and MS susceptibility. Comparisons between these previous studies and our current study are summarised in Table 13. GPC5-rs10492503 showed a significant association in the total disease population. When analysed further, we found significant associations with two disease states (SPMS and RRMS) and in the total female population and the female SPMS and RRMS subtypes. GPC6-rs17267815 showed a minor significant association within the RRMS subtype only. Further analysis suggested that this association was due to the male RRMS subgroup; however, due to low sample numbers, once the population was stratified, significance values are suggestive only.

Table 13 Comparison of significance obtained in this study compared to previous GWAS. P values from GWAS presented as from the original paper. Pun = uncorrected P value; PC = corrected p value. Baranzini et al. presented their significance as adjusted log P values

We identified no LD between the SNPs studied within the previously identified 13q31-32 risk region containing both these genes, nor could we replicate the moderate LD identified previously in GPC5 [51]. All five GPC5 and GPC6 SNPs investigated in this study had previously been identified as significant in large-scale case/control GWAS and replication studies in Norwegian and Spanish populations [23, 51, 52, 57] with varying and often contradictory levels of significance. SNPs reaching significance in one study were not found to be significant in another [23, 51, 53]. Analysis by the disease state of some of these populations determined significant associations with the RRMS subtype [66]. Indeed, in our population, when a significant association was observed in these genes, it was often significant in the RRMS sub-population. This may be due to the mixture of the populations as patients from pure Northern European ancestry have a higher risk of developing MS [8]. While our Australian population is of Caucasian descent, it is not necessarily of purely Northern European origin, potentially explaining some differences between results and levels of significance identified in these studies. In addition, while our results are not strongly significant on their own, they replicate previous studies and support and strengthen the evidence for the involvement of GPC5 and GPC6 in the development and progression of MS.

Many HSPGs and their associated enzymes have been associated with the disease, with both SDC1 and SDC4 showing strong involvement with breast cancer [46, 67, 68]. As yet, the physiological functions of the GPCs, in both normal and pathological conditions, remain poorly understood. However, data here and in other studies suggest an important function for these proteins in cell growth and regulation of division. Celie and colleagues suggest that HSPGs are involved in the inflammatory response and have a regulatory role in leukocyte extravasation [69], a condition synonymous with MS. Other GPCs have been shown to play roles in diseases such as hepatocellular carcinoma (GPC3 [70];) and Simpson-Golabi-Behmel syndrome (GPC3/GPC4 [71];). While the function of GPC5 remains poorly understood, particularly in MS, different GPC polymorphisms have been reported to increase the risk of lung cancer in non-smokers [72] while decreasing the risk of cancer in MS patients, with this reduced cancer risk stemming from the specific GPC gene [73]. In addition, the gene region 13q31-32 containing both GPC5 and GPC6 has also previously been associated with the increased risk of primary sclerosing cholangitis (PSC), a chronic liver disease where a strong association has been identified between the SNP GPC6-rs9524260 and disease [58].

Due to the interaction of the GPCs with several growth factors, chemokines and ECM proteins, there may also be an effect on neural growth and repair [74]. The results of a study by Cenit et al. [57] not only supported a significant association of GPC5-rs10492503 with MS but also indicated approximately twice the risk of developing the disease in an individual who has one or more copies of the variant allele [57]. GPC5 has been reported to play an important role during the process of cell division and growth regulation. It is predominantly expressed in foetal tissues, including the brain, lung, liver and kidney. However, it has an exclusive expression in adult tissue in the CNS and in its neurons [71, 75]. This suggests a possible and plausible role for this gene in controlling various neurotropic factors and maintenance of neural function. In our study, we found a significant association of this GPC5 variation with the early-onset form of the disease (RRMS) and also the severe form (SPMS), a progression of the disease characterised by irreversible damage suggesting a role for GPC5 in the progression of MS. GPC5 plays an important role in brain patterning, synapse formation, axon regeneration and guidance. Its expression in the developing brain and the adult CNS (the origin of MS) also support a role for this gene in different disease states.

In PPMS, most of the myelin degradation occurs in the cerebrum and cerebellar cortex of the CNS [1]. Dysfunction of GPC5 could affect cell proliferation and tissue growth. With the cells no longer able to interact with the various positively charged growth factors, this would affect brain patterning, synapse formation and an interruption in axon regeneration. This suggests that abnormal GPC5 may play a role in triggering MS and the subsequent disability experienced by sufferers.

Further evidence supporting this hypothesis is data demonstrating that HSPGs have been identified in the active lesions of MS, where they are thought to be involved in the sequestering of pro-inflammatory chemokines [55]. GPC5 expression and interaction with various growth factors and chemokines likely affect growth and repair of neurons, also influencing the guidance of axons and synapse formation [28, 55, 76]. Indeed, another member of the GPC family—GPC1, has been shown to be required for Schwann cell myelination [77]. With documented involvement of other GPCs, it is plausible allelic variants of GPC5 may affect neuronal repair, axon guidance and new synaptic formation.

The embryonic expression of GPC6 is detected in the ovary, liver and kidneys, while in the adult, it is detected only in the ovary and intestine [78]. More recent emerging evidence for GPC6 may indicate a role for the gene in neural diseases with origins in the CNS, the location of the MS-associated lesions. Overall, the functional role for GPC6 is poorly understood, but this study provides some evidence of a potential role for GPC6 in MS.

Conclusion

From this study, we have determined a significant association with the rs1131351 SNP in SDC1, specifically in females suffering from either primary progressive or relapsing-remitting forms of MS. Involvement of SDC1 in the initiation of MS has been suggested through its involvement in the inflammatory response and growth factor interactions. Levels of specific growth factors may vary during MS onset which could be due to dysfunction of HSPGs brought about by their inability to appropriately traffic/sequester growth factors. The specific mechanism of GPC5 and GPC6 involvement in MS has yet to be elucidated. However, a number of genetic studies, this one included, have provided evidence suggesting a role for these genes in the progression of the disease. Evidence already exists for these genes in other diseases utilising similar mechanisms of action. Significant results obtained in this study have been summarised in Fig. 2. This schematic highlights the involvement of specific SNPs in specific disease states as well as whether the SNP is associated with either the male or female subpopulation. Results from this study are by no means conclusive, but they add to the growing body of evidence indicating the involvement of these PGs in the initiation and progression of neurodegenerative diseases. Specifically, this study supports and strengthens evidence suggesting a role for HSPG core proteins, both syndecans and glypicans, in the development and progression of MS.

Fig. 2
figure 2

Schematic of disease progression highlighting the involvement of HSPG core protein SNPs. Significant associations between specific HSPG core protein SNPs and disease subtypes are represented by coloured stars as can be seen in the legend. SDC1: yellow; GPC5: green; GPC6: purple. These SNPs are significantly associated with population stratification by gender as represented by coloured people outlines. Female: pink; Male: blue

Availability of data and materials

All data relevant for this study are included within this manuscript, any further information may be made available on request.

References

  1. Kutzelnigg A, Lucchinetti CF, Stadelmann C, Bruck W, Rauschka H, Bergmann M, Schmidbauer M, Parisi JE, Lassmann H. Cortical demyelination and diffuse white matter injury in multiple sclerosis. Brain. 2005;128(Pt 11):2705–12.

    Article  PubMed  Google Scholar 

  2. Noseworthy JH, Lucchinetti C, Rodriguez M, Weinshenker BG. Multiple sclerosis. N Engl J Med. 2000;343:938–52.

  3. Liguori M, Marrosu MG, Pugliatti M, Giuliani F, De Robertis F, Cocco E, Zimatore GB, Livrea P, Trojano M. Age at onset in multiple sclerosis. Neurol Sci. 2000;21:825–9.

  4. Tajouri L, Fernandez F, Tajouri S, Detriche G, Szvetko A, Colson N, Csurhes P, Pender MP, Griffiths LR. Allelic variation investigation of the estrogen receptor within an Australian multiple sclerosis population. J Neurol Sci. 2007;252(1):9–12.

    Article  CAS  PubMed  Google Scholar 

  5. Ahmad H, Palmer A, Campbell JA, van der Mei I, Taylor B. Health economic impact of multiple sclerosis in Australia in 2017. Hobart: Menzies Institute for Medial Research; 2018.

    Google Scholar 

  6. Browne P, Chandraratna D, Angood C, Tremlett H, Baker C, Taylor BV, Thompson AJ. Atlas of multiple sclerosis 2013: a growing global problem with widespread inequity. Neurology. 2014;83(11):1022–4.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Australia, M. S. S. o. 2017. "Symptoms." 2019, from https://www.msaustralia.org.au/about-ms/symptoms.

  8. Kurtzke JF, Beebe GW, Norman JE Jr. Epidemiology of multiple sclerosis in U.S. veterans: 1. Race, sex, and geographic distribution. Neurology. 1979;29(9 Pt 1):1228–35.

    Article  CAS  PubMed  Google Scholar 

  9. Zuvich RL, McCauley JL, Pericak-Vance MA, Haines JL. Genetics and pathogenesis of multiple sclerosis. Semin Immunol. 2009;21(6):328–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Gooch CL, Pracht E, Borenstein AR. The burden of neurological disease in the United States: a summary report and call to action. Ann Neurol. 2017;81(4):479–84.

    Article  PubMed  Google Scholar 

  11. Oksenberg JR, Baranzini SE, Sawcer S, Hauser SL. The genetics of multiple sclerosis: SNPs to pathways to pathogenesis. Nat Rev Genet. 2008;9(7):516–26.

    Article  CAS  PubMed  Google Scholar 

  12. Willer CJ, Dyment DA, Sadovnick AD, Rothwell PM, Murray TJ, Ebers GC, G. Canadian Collaborative Study. Timing of birth and risk of multiple sclerosis: population based study. BMJ. 2005;330(7483):120.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Compston A, Coles A. Multiple sclerosis. Lancet. 2008;372(9648):1502–17.

    Article  CAS  PubMed  Google Scholar 

  14. Trapp BD, Nave KA. Multiple sclerosis: an immune or neurodegenerative disorder? Annu Rev Neurosci. 2008;31:247–69.

    Article  CAS  PubMed  Google Scholar 

  15. Rosati G. The prevalence of multiple sclerosis in the world: an update. Neurol Sci. 2001;22(2):117–39.

    Article  CAS  PubMed  Google Scholar 

  16. Hauser SL, Oksenberg JR. The neurobiology of multiple sclerosis: genes, inflammation, and neurodegeneration. Neuron. 2006;52(1):61–76.

    Article  CAS  PubMed  Google Scholar 

  17. Huijbregts SC, Kalkers NF, de Sonneville LM, de Groot V, Polman CH. Cognitive impairment and decline in different MS subtypes. J Neurol Sci. 2006;245(1-2):187–94.

    Article  PubMed  Google Scholar 

  18. Huijbregts SC, Kalkers NF, de Sonneville LM, de Groot V, Reuling IE, Polman CH. Differences in cognitive impairment of relapsing remitting, secondary, and primary progressive MS. Neurology. 2004;63(2):335–9.

    Article  CAS  PubMed  Google Scholar 

  19. Brady ST, Witt AS, Kirkpatrick LL, de Waegh SM, Readhead C, Tu PH, Lee VM. Formation of compact myelin is required for maturation of the axonal cytoskeleton. J Neurosci. 1999;19(17):7278–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Olerup O, Hillert J. HLA class II-associated genetic susceptibility in multiple sclerosis: a critical evaluation. Tissue Antigens. 1991;38(1):1–15.

    Article  CAS  PubMed  Google Scholar 

  21. Haines JL, Terwedow HA, Burgess K, Pericak-Vance MA, Rimmler JB, Martin ER, Oksenberg JR, Lincoln R, Zhang DY, Banatao DR, Gatto N, Goodkin DE, Hauser SL. Linkage of the MHC to familial multiple sclerosis suggests genetic heterogeneity. The Multiple Sclerosis Genetics Group. Hum Mol Genet. 1998;7(8):1229–34.

    Article  CAS  PubMed  Google Scholar 

  22. Oksenberg JR, Baranzini SE. Multiple sclerosis genetics--is the glass half full, or half empty? Nat Rev Neurol. 2010;6(8):429–37.

    Article  PubMed  Google Scholar 

  23. Baranzini SE, Wang J, Gibson RA, Galwey N, Naegelin Y, Barkhof F, Radue EW, Lindberg RL, Uitdehaag BM, Johnson MR, Angelakopoulou A, Hall L, Richardson JC, Prinjha RK, Gass A, Geurts JJ, Kragt J, Sombekke M, Vrenken H, Qualley P, Lincoln RR, Gomez R, Caillier SJ, George MF, Mousavi H, Guerrero R, Okuda DT, Cree BA, Green AJ, Waubant E, Goodin DS, Pelletier D, Matthews PM, Hauser SL, Kappos L, Polman CH, Oksenberg JR. Genome-wide association analysis of susceptibility and clinical phenotype in multiple sclerosis. Hum Mol Genet. 2009;18(4):767–78.

    Article  CAS  PubMed  Google Scholar 

  24. Gregory SG, Schmidt S, Seth P, Oksenberg JR, Hart J, Prokop A, Caillier SJ, Ban M, Goris A, Barcellos LF, Lincoln R, McCauley JL, Sawcer SJ, Compston DA, Dubois B, Hauser SL, Garcia-Blanco MA, Pericak-Vance MA, Haines JL, G. Multiple Sclerosis Genetics. Interleukin 7 receptor alpha chain (IL7R) shows allelic and functional association with multiple sclerosis. Nat Genet. 2007;39(9):1083–91.

    Article  CAS  PubMed  Google Scholar 

  25. International Multiple Sclerosis Genetics, Hafler CDA, Compston A, Sawcer S, Lander ES, Daly MJ, De Jager PL, de Bakker PI, Gabriel SB, Mirel DB, Ivinson AJ, Pericak-Vance MA, Gregory SG, Rioux JD, McCauley JL, Haines JL, Barcellos LF, Cree B, Oksenberg JR, Hauser SL. Risk alleles for multiple sclerosis identified by a genomewide study. N Engl J Med. 2007;357(9):851–62.

    Article  Google Scholar 

  26. Lundmark F, Duvefelt K, Iacobaeus E, Kockum I, Wallstrom E, Khademi M, Oturai A, Ryder LP, Saarela J, Harbo HF, Celius EG, Salter H, Olsson T, Hillert J. Variation in interleukin 7 receptor alpha chain (IL7R) influences risk of multiple sclerosis. Nat Genet. 2007;39(9):1108–13.

    Article  CAS  PubMed  Google Scholar 

  27. Xue HH, Kovanen PE, Pise-Masison CA, Berg M, Radovich MF, Brady JN, Leonard WJ. IL-2 negatively regulates IL-7 receptor alpha chain expression in activated T lymphocytes. Proc Natl Acad Sci U S A. 2002;99(21):13759–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Lee JS, Chien CB. When sugars guide axons: insights from heparan sulphate proteoglycan mutants. Nat Rev Genet. 2004;5(12):923–35.

    Article  CAS  PubMed  Google Scholar 

  29. Cassaro CM, Dietrich CP. Distribution of sulfated mucopolysaccharides in invertebrates. J Biol Chem. 1977;252(7):2254–61.

    CAS  PubMed  Google Scholar 

  30. Nader HB, Chavante SF, dos-Santos EA, Oliveira TW, de-Paiva JF, Jeronimo SM, Medeiros GF, de-Abreu LR, Leite EL, de-Sousa-Filho JF, Castro RA, Toma L, Tersariol IL, Porcionatto MA, Dietrich CP. Heparan sulfates and heparins: similar compounds performing the same functions in vertebrates and invertebrates? Braz J Med Biol Res. 1999;32(5):529–38.

    Article  CAS  PubMed  Google Scholar 

  31. Habuchi H, Habuchi O, Kimata K. Sulfation pattern in glycosaminoglycan: does it have a code? Glycoconj J. 2004;21(1-2):47–52.

    Article  CAS  PubMed  Google Scholar 

  32. Inatani M, Yamaguchi Y. Gene expression of EXT1 and EXT2 during mouse brain development. Brain Res Dev Brain Res. 2003;141(1-2):129–36.

    Article  CAS  PubMed  Google Scholar 

  33. Jennes I, Zuntini M, Mees K, Palagani A, Pedrini E, De Cock G, Fransen E, Vanden Berghe W, Sangiorgi L, Wuyts W. Identification and functional characterization of the human EXT1 promoter region. Gene. 2012;492(1):148–59.

    Article  CAS  PubMed  Google Scholar 

  34. Holst CR, Bou-Reslan H, Gore BB, Wong K, Grant D, Chalasani S, Carano RA, Frantz GD, Tessier-Lavigne M, Bolon B, French DM, Ashkenazi A. Secreted sulfatases Sulf1 and Sulf2 have overlapping yet essential roles in mouse neonatal survival. PLoS One. 2007;2(6):e575.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Lin X, Wei G, Shi Z, Dryer L, Esko JD, Wells DE, Matzuk MM. Disruption of gastrulation and heparan sulfate biosynthesis in EXT1-deficient mice. Dev Biol. 2000;224(2):299–311.

    Article  CAS  PubMed  Google Scholar 

  36. Dreyfuss JL, Regatieri CV, Jarrouge TR, Cavalheiro RP, Sampaio LO, Nader HB. Heparan sulfate proteoglycans: structure, protein interactions and cell signaling. An Acad Bras Cienc. 2009;81(3):409–29.

    Article  CAS  PubMed  Google Scholar 

  37. Han CH, Huang YJ, Lu KH, Liu Z, Mills GB, Wei Q, Wang LE. Polymorphisms in the SULF1 gene are associated with early age of onset and survival of ovarian cancer. J Exp Clin Cancer Res. 2011;30:5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Sahota AP, Dhoot GK. A novel SULF1 splice variant inhibits Wnt signalling but enhances angiogenesis by opposing SULF1 activity. Exp Cell Res. 2009;315(16):2752–64.

    Article  CAS  PubMed  Google Scholar 

  39. Nagamine S, Tamba M, Ishimine H, Araki K, Shiomi K, Okada T, Ohto T, Kunita S, Takahashi S, Wismans RG, van Kuppevelt TH, Masu M, Keino-Masu K. Organ-specific sulfation patterns of heparan sulfate generated by extracellular sulfatases Sulf1 and Sulf2 in mice. J Biol Chem. 2012;287(12):9579–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Haupt LM, Griffiths LR. Heparan Sulfate Proteoglycans, Tumour Progression and the Cancer Stem Cell Niche. Curr Cancer Ther Rev. 2009;5(4):256–60.

    Article  CAS  Google Scholar 

  41. Haupt LM, Murali S, Mun FK, Teplyuk N, Mei LF, Stein GS, van Wijnen AJ, Nurcombe V, Cool SM. The heparan sulfate proteoglycan (HSPG) glypican-3 mediates commitment of MC3T3-E1 cells toward osteogenesis. J Cell Physiol. 2009;220(3):780–91.

    Article  CAS  PubMed  Google Scholar 

  42. Lindahl U, Kusche-Gullberg M, Kjellen L. Regulated diversity of heparan sulfate. J Biol Chem. 1998;273(39):24979–82.

    Article  CAS  PubMed  Google Scholar 

  43. Tumova S, Hatch BA, Law DJ, Bame KJ. Basic fibroblast growth factor does not prevent heparan sulphate proteoglycan catabolism in intact cells, but it alters the distribution of the glycosaminoglycan degradation products. Biochem J. 1999;337(Pt 3):471–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Filmus J, Shi W, Wong ZM, Wong MJ. Identification of a new membrane-bound heparan sulphate proteoglycan. Biochem J. 1995;311(Pt 2):561–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Malavaki CJ, Theocharis AD, Lamari FN, Kanakis I, Tsegenidis T, Tzanakakis GN, Karamanos NK. Heparan sulfate: biological significance, tools for biochemical analysis and structural characterization. Biomed Chromatogr. 2011;25(1-2):11–20.

    Article  CAS  PubMed  Google Scholar 

  46. Tkachenko E, Rhodes JM, Simons M. Syndecans: new kids on the signaling block. Circ Res. 2005;96(5):488–500.

    Article  CAS  PubMed  Google Scholar 

  47. David G. Integral membrane heparan sulfate proteoglycans. FASEB J. 1993;7(11):1023–30.

    Article  CAS  PubMed  Google Scholar 

  48. Bernfield M, Gotte M, Park PW, Reizes O, Fitzgerald ML, Lincecum J, Zako M. Functions of cell surface heparan sulfate proteoglycans. Annu Rev Biochem. 1999;68:729–77.

    Article  CAS  PubMed  Google Scholar 

  49. Ling L, Murali S, Dombrowski C, Haupt LM, Stein GS, van Wijnen AJ, Nurcombe V, Cool SM. Sulfated glycosaminoglycans mediate the effects of FGF2 on the osteogenic potential of rat calvarial osteoprogenitor cells. J Cell Physiol. 2006;209(3):811–25.

    Article  CAS  PubMed  Google Scholar 

  50. Filmus J, Capurro M, Rast J. Glypicans. Genome Biol. 2008;9(5):224.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  51. Lorentzen AR, Melum E, Ellinghaus E, Smestad C, Mero IL, Aarseth JH, Myhr KM, Celius EG, Lie BA, Karlsen TH, Franke A, Harbo HF. Association to the Glypican-5 gene in multiple sclerosis. J Neuroimmunol. 2010;226(1-2):194–7.

    Article  CAS  PubMed  Google Scholar 

  52. Comabella M, Craig DW, Camina-Tato M, Morcillo C, Lopez C, Navarro A, Rio J, Biomarker MSSG, Montalban X, Martin R. Identification of a novel risk locus for multiple sclerosis at 13q31.3 by a pooled genome-wide scan of 500,000 single nucleotide polymorphisms. PLoS One. 2008;3(10):e3490.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Cavanillas ML, Fernandez O, Comabella M, Alcina A, Fedetz M, Izquierdo G, Lucas M, Cenit MC, Arroyo R, Vandenbroeck K, Alloza I, Garcia-Barcina M, Antiguedad A, Leyva L, Gomez CL, Olascoaga J, Otaegui D, Blanco Y, Saiz A, Montalban X, Matesanz F, Urcelay E. Replication of top markers of a genome-wide association study in multiple sclerosis in Spain. Genes Immun. 2011;12(2):110–5.

    Article  CAS  PubMed  Google Scholar 

  54. Tajouri L, Mellick AS, Tourtellotte A, Nagra RM, Griffiths LR. An examination of MS candidate genes identified as differentially regulated in multiple sclerosis plaque tissue, using absolute and comparative real-time Q-PCR analysis. Brain Res Brain Res Protoc. 2005;15(2):79–91.

    Article  CAS  PubMed  Google Scholar 

  55. van Horssen J, Bo L, Dijkstra CD, de Vries HE. Extensive extracellular matrix depositions in active multiple sclerosis lesions. Neurobiol Dis. 2006;24(3):484–91.

    Article  PubMed  CAS  Google Scholar 

  56. Okolicsanyi RK, van Wijnen AJ, Cool SM, Stein GS, Griffiths LR, Haupt LM. "Heparan sulfate proteoglycans and human breast cancer epithelial cell tumorigenicity." J Cell Biochem. 2014;115(5):967–76.

  57. Cenit MD, Blanco-Kelly F, de las Heras V, Bartolome M, de la Concha EG, Urcelay E, Arroyo R, Martinez A. Glypican 5 is an interferon-beta response gene: a replication study. Mult Scler. 2009;15(8):913–7.

    Article  CAS  PubMed  Google Scholar 

  58. Karlsen TH, Franke A, Melum E, Kaser A, Hov JR, Balschun T, Lie BA, Bergquist A, Schramm C, Weismuller TJ, Gotthardt D, Rust C, Philipp EE, Fritz T, Henckaerts L, Weersma RK, Stokkers P, Ponsioen CY, Wijmenga C, Sterneck M, Nothnagel M, Hampe J, Teufel A, Runz H, Rosenstiel P, Stiehl A, Vermeire S, Beuers U, Manns MP, Schrumpf E, Boberg KM, Schreiber S. Genome-wide association analysis in primary sclerosing cholangitis. Gastroenterology. 2010;138(3):1102–11.

    Article  PubMed  Google Scholar 

  59. Busse M, Kusche-Gullberg M. In vitro polymerization of heparan sulfate backbone by the EXT proteins. J Biol Chem. 2003;278(42):41333–7.

    Article  CAS  PubMed  Google Scholar 

  60. Kim BT, Kitagawa H, Tanaka J, Tamura J, Sugahara K. In vitro heparan sulfate polymerization: crucial roles of core protein moieties of primer substrates in addition to the EXT1-EXT2 interaction. J Biol Chem. 2003;278(43):41618–23.

    Article  CAS  PubMed  Google Scholar 

  61. Busse-Wicher M, Wicher KB, Kusche-Gullberg M. The exostosin family: proteins with many functions. Matrix Biol. 2014;35:25–33.

    Article  CAS  PubMed  Google Scholar 

  62. Dutta R, Trapp BD. Relapsing and progressive forms of multiple sclerosis: insights from pathology. Curr Opin Neurol. 2014;27(3):271–8.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Goldenberg MM. Multiple sclerosis review. P T. 2012;37(3):175–84.

    PubMed  PubMed Central  Google Scholar 

  64. van Horssen J, Dijkstra CD, de Vries HE. The extracellular matrix in multiple sclerosis pathology. J Neurochem. 2007;103(4):1293–301.

    Article  PubMed  CAS  Google Scholar 

  65. Titelbaum DS, Degenhardt A, Kinkel RP. Anti-tumor necrosis factor alpha-associated multiple sclerosis. AJNR Am J Neuroradiol. 2005;26(6):1548–50.

    PubMed  PubMed Central  Google Scholar 

  66. Poliseno L, Salmena L, Zhang J, Carver B, Haveman WJ, Pandolfi PP. A coding-independent function of gene and pseudogene mRNAs regulates tumour biology. Nature. 2010;465(7301):1033–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Lendorf ME, Manon-Jensen T, Kronqvist P, Multhaupt HA, Couchman JR. Syndecan-1 and syndecan-4 are independent indicators in breast carcinoma. J Histochem Cytochem. 2011;59(6):615–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Okolicsanyi RK, van Wijnen AJ, Cool SM, Stein GS, Griffiths LR, Haupt LM. Heparan sulfate proteoglycans and human breast cancer epithelial cell tumorigenicity. J Cell Biochem. 2014;115(5):967–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Celie JW, Beelen RH, van den Born J. Heparan sulfate proteoglycans in extravasation: assisting leukocyte guidance. Front Biosci (Landmark Ed). 2009;14:4932–49.

    Article  CAS  Google Scholar 

  70. Capurro MI, Xiang YY, Lobe C, Filmus J. Glypican-3 promotes the growth of hepatocellular carcinoma by stimulating canonical Wnt signaling. Cancer Res. 2005;65(14):6245–54.

    Article  CAS  PubMed  Google Scholar 

  71. Veugelers M, De Cat B, Ceulemans H, Bruystens AM, Coomans C, Durr J, Vermeesch J, Marynen P, David G. Glypican-6, a new member of the glypican family of cell surface heparan sulfate proteoglycans. Jbb Biol Chem. 1999;274(38):26968–77.

    Article  CAS  Google Scholar 

  72. Li Y, Sheu CC, Ye Y, de Andrade M, Wang L, Chang SC, Aubry MC, Aakre JA, Allen MS, Chen F, Cunningham JM, Deschamps C, Jiang R, Lin J, Marks RS, Pankratz VS, Su L, Li Y, Sun Z, Tang H, Vasmatzis G, Harris CC, Spitz MR, Jen J, Wang R, Zhang ZF, Christiani DC, Wu X, Yang P. Genetic variants and risk of lung cancer in never smokers: a genome-wide association study. Lancet Oncol. 2010;11(4):321–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Handel AE, Ramagopalan SV. GPC5 and lung cancer in multiple sclerosis. Lancet Oncol. 2010;11(8):714.

    Article  PubMed  Google Scholar 

  74. Byun E, Caillier SJ, Montalban X, Villoslada P, Fernandez O, Brassat D, Comabella M, Wang J, Barcellos LF, Baranzini SE, Oksenberg JR. Genome-wide pharmacogenomic analysis of the response to interferon beta therapy in multiple sclerosis. Arch Neurol. 2008;65(3):337–44.

    Article  PubMed  Google Scholar 

  75. Saunders S, Paine-Saunders S, Lander AD. Expression of the cell surface proteoglycan glypican-5 is developmentally regulated in kidney, limb, and brain. Dev Biol. 1997;190(1):78–93.

    Article  CAS  PubMed  Google Scholar 

  76. Van Vactor D, Wall DP, Johnson KG. Heparan sulfate proteoglycans and the emergence of neuronal connectivity. Curr Opin Neurobiol. 2006;16(1):40–51.

    Article  PubMed  CAS  Google Scholar 

  77. Chernousov MA, Rothblum K, Stahl RC, Evans A, Prentiss L, Carey DJ. Glypican-1 and alpha4(V) collagen are required for Schwann cell myelination. J Neurosci. 2006;26(2):508–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Fransson LA. Glypicans. Int J Biochem Cell Biol. 2003;35(2):125–9.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by infrastructure purchased with Australian Government EIF Super Science Funds as part of the Therapeutic Innovation Australia—Queensland Node project. RKO was supported by an Australian Postgraduate Award and the authors wish to acknowledge the support of the Estate of the late Clem Jones AO (LMH, LRG).

Author information

Authors and Affiliations

Authors

Contributions

Conception and design of the manuscript were completed by RKO and LMH. Design and completion of genotyping assays were conducted by RKO, JB and CM. Analysis of HRM melt curves, RFLP fragment separation and pyrosequencing traces were conducted by RKO, JB and CM. Writing and editing of the manuscript were completed by RKO with substantial editing contributed by LMH and LRG. All authors approved the submission of the manuscript to Human Genetics.

Corresponding authors

Correspondence to Lyn R. Griffiths or Larisa M. Haupt.

Ethics declarations

Ethics approval and consent to participate

Ethical approval for these studies is through Queensland University of Technology (QUT) Human Research Ethics Committee HREC (Approval Number 1400000502) and previously through Griffith University Human Research Ethics Committee (1300000488).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests regarding the information in this manuscript.

Supplementary Table 1. Male results by disease state for GPC6, rs17267815

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: Table S1.

Male results by disease state for GPC6, rs17267815.

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

Okolicsanyi, R.K., Bluhm, J., Miller, C. et al. An investigation of genetic polymorphisms in heparan sulfate proteoglycan core proteins and key modification enzymes in an Australian Caucasian multiple sclerosis population. Hum Genomics 14, 18 (2020). https://doi.org/10.1186/s40246-020-00264-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40246-020-00264-6

Keywords