- 1Centre of Excellence for Nutrition, North-West University, Potchefstroom, South Africa
- 2Medical Research Council Unit for Hypertension and Cardiovascular Disease, North-West University, Potchefstroom, South Africa
- 3Africa Unit for Transdisciplinary Health Research (AUTHeR), North-West University, Potchefstroom, South Africa
The role of 25-hydroxyvitamin D [25(OH)D] in reducing the risk of cardiovascular disease (CVD) has been recognized, but the mechanisms involved are unclear. Researchers have discovered a link between vitamin D and fibrinogen. Until now, data on the relationship between vitamin D and the γ' splice variant of fibrinogen and fibrin clot characteristics remain unexplored. In this study, 25(OH)D, total and γ' fibrinogen, as well as turbidimetrically determined plasma clot properties, were quantified, and fibrinogen and FXIII SNPs were genotyped in 660 Black, apparently healthy South African women. Alarmingly, 16 and 45% of the women presented with deficient and insufficient 25(OH)D, respectively. Total fibrinogen and maximum absorbance (as a measure of clot density) correlated inversely, whereas γ' fibrinogen correlated positively with 25(OH)D. γ' fibrinogen increased whereas maximum absorbance decreased over the deficient, insufficient, and sufficient 25(OH)D categories before and after adjustment for confounders. 25(OH)D modulated the association of the SNPs regarding fibrinogen concentration and clot structure/properties, but did not stand after correction for false discovery rate. Because only weak relationships were detected, the clinical significance of the findings are questionable and remain to be determined. However, we recommend vitamin D fortification and supplementation to reduce the high prevalence of this micronutrient deficiency and possibly to improve fibrinogen and plasma clot structure if the relationships are indeed clinically significant. There is a need for large cohort studies to demonstrate the relationship between vitamin D and cardiovascular and inflammatory risk factors as well as to uncover the molecular mechanisms responsible.
Introduction
Although cardiovascular (CVD) and other non-communicable diseases are major causes of mortality (1, 2), these conditions also increase the risk of dying of communicable diseases such as coronavirus disease (COVID-19) (3–5). A link between deficiency of vitamin D—a group of fat-soluble secosteroids with the most significant being vitamin D3 (cholecalciferol) and D2 (ergocalciferol), well-known for their role in bone health—and CVD, has been proposed (6–8). However, the underlying mechanisms remain to be clarified. One route that has been recognized is that vitamin D affords protection against circulatory diseases by reducing inflammation (9–12). Less-studied pathways include improving the blood lipid profile (13), downregulation of the renin-angiotensin system improving blood pressure (14, 15), enhancement of insulin secretion and insulin sensitivity (16, 17), protection against uncontrolled angiogenesis (18) and influencing blood coagulation (19, 20).
As vitamin D can be acquired through diet and cholecalciferol synthesized in a chemical reaction when the skin is exposed to sunlight, circulating concentrations of 25-hydroxyvitamin D-−25(OH)D, the sum of vitamin D2 and D3 and the precursor of the active form 1,25-dihydroxyvitamin D3 (1,25(OH)D)—is the hallmark of assessing vitamin D status. The relationships between 25(OH)D concentrations and inflammatory markers, including the glycoprotein fibrinogen, have been investigated in numerous studies, with most sampling small or specific patient groups (19–24). Data for the association between vitamin D and fibrinogen in generally healthy populations is sparse, and evidence from continental Africans, with darker skins that hinder vitamin D production, is limited with some reports on African Americans (12, 25). Whereas some researchers have detected an inverse relationship between vitamin D and fibrinogen (19, 20, 24, 26, 27), information regarding the fibrinogen isoform that stems from a splice variant of the γ-chain messenger RNA resulting from an alternative polyadenylation signal in intron 9, i.e., fibrinogen γ', and clot properties is lacking. Investigating the relationship between fibrinogen and vitamin D is crucial because fibrinogen is the dominant coagulation protein in blood by mass and the precursor of fibrin, the main constituent of blood clots (28–30). Studying fibrinogen γ' is vital because of its influence on clot properties (31–33). Factor XIII (FXIII) is another blood coagulation factor that influences fibrin clot structure and functions by cross-linking fibrin fibers to prevent premature lysis of the clot to allow wound healing (34, 35). Furthermore, total and γ' fibrinogen concentrations, as well as clot structure and characteristics, may be influenced by single nucleotide polymorphisms (SNPs) within the fibrinogen α, β and γ chain genes (FGA, FGB and FGG) as well as those of FXIII (36). Therefore, fibrinogen and FXIII SNPs might interact with 25(OH)D to modulate total or γ' fibrinogen concentration as well as the kinetics of clot formation, structural clot properties and fibrinolysis and warrant consideration.
Thus, we first investigated the associations of 25(OH)D determined in a group of ostensibly healthy Black women with total and γ' fibrinogen and fibrin clot properties. Second, we determined how these associations are influenced by specific genetic variations within the FGA, FGB, FGG and FXIII genes. To provide possible explanations for the observed relationships we considered the inflammatory marker C-reactive protein (CRP) and glycated hemoglobin A1c (HbA1c) which is a long-term indicator of glucose control. By controlling for CRP and HbA1c we explored whether the underlying mechanisms were related to inflammation or to vitamin D's effect on glucose control (16, 17). This research is pressing because unraveling probable mechanisms for the observed relationship between vitamin D and CVD could inform the design of targeted approaches to prevent CVD in those most vulnerable.
Materials and Methods
Research Design and Ethical Considerations
The Prospective Urban and Rural Epidemiology (PURE) study investigated the development of chronic lifestyle diseases, especially CVDs, by stratifying populations at different levels of urbanization (37). For this cross-sectional investigation, we used a subset of 660 randomly selected women derived from the 1,263 healthy women older than 30 years who participated in the North West site of the South African arm of the PURE study at baseline. Goodwill permission was granted by mayors, community leaders and tribal chiefs before researchers accessed the communities included for recruitment. Eligible individuals were well-informed regarding all aspects of the research and, after ample time for reflection, were asked to sign an informed consent form indicating their agreement to join in the study. Final enrolment and measurements for baseline were scheduled between August and November 2005, which is late winter to late spring in the Southern Hemisphere. Due to monetary constraints vitamin D was not measured in the follow-up blood samples. Persons with elevated tympanic temperatures (above 38°C) were excluded. Additional exclusion criteria were using chronic medication (apart from blood pressure medicine), having any known lifestyle disease, being pregnant or lactating, or having a known infection, such as tuberculosis-causing agents and/or human immunodeficiency virus (HIV). However, participants diagnosed with HIV during recruitment were included in the study and, because total fibrinogen, %γ' fibrinogen, lag time, slope, maximum absorbance, clot lysis time (CLT) and 25(OH)D did not differ between infected or non-infected, were also included in the statistical analysis.
Ethical approval of the study protocol was obtained for the larger study as well as this nested study from the Health Research Ethics Committee of the Faculty of Health Sciences, North-West University (NWU–HREC, ethics number: 04M10; NWU-00016-10-A1 and NWU-00034-17-A1-02), honoring the guidelines laid down in the Declaration of Helsinki of 1975 and revised in 2000.
Blood Pressure Measurements
Blood pressure—systolic and diastolic—was measured twice using the OMRON HEM-757 automatic apparatus (Omron Healthcare, Kyoto, Japan) with the cuff on the left arm after 5 min of sitting in a calm environment. The second measurement was recorded.
Anthropometric Measurements
While participants were lightly clothed and their arms hung freely at their sides, their body weights were recorded twice. Heights were measured twice using stadiometers with volunteers' heads in the Frankfort plane, bodies fully extended while inhaling. The means of the participants' weight and height measurements were used, to compute body mass index. Waist circumference was determined at the narrowest point between the 10th rib and iliac crest. Hip circumference was assessed at the level where the protrusion of the buttocks was greatest.
Blood Collection, Sampling, and Storage
Blood samples were collected between 07:00 and 11:00 from fasting participants. Sodium citrate-treated (3.8%) whole blood was centrifuged at 2,000 x g for 15 min to yield a buffy coat for DNA extraction and plasma samples to quantify coagulation markers. Ethylenediamine tetraacetic acid-treated whole blood was collected for HbA1c determination. Coagulated blood was collected to yield serum for lipids, CRP, and γ-glutamyl transferase (GGT) analysis. Aliquots were frozen on dry ice pellets, stored in the field at −18°C and then, after 2–4 days, at −80°C until analyzed.
Biochemical Analysis
Lipograms [including total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C) and triglycerides], high sensitivity CRP and GGT concentrations were assayed using a Sequential Multiple Analyzer (Konelab 20i, Thermo Fisher Scientific Oy, Vantaa, Finland). Low-density lipoprotein cholesterol (LDL-C) was calculated using the Friedewald formula for triglyceride concentrations below 400 mg/dL. HbA1 concentrations were determined using the D-10 Hemoglobin testing system (Bio-Rad Laboratories, Hercules, CA, USA). Serum total 25(OH)D concentrations were quantified using a Roche Elecsys 2010 COBAS system (functional sensitivity: 10.0 nmol/L; Roche Diagnostics, Indianapolis, IN, USA). Whole blood was used for a rapid first response HIV card test 1–2.0 (Transnational Technologies Inc. PMC Medical, Nani Daman, India), and the outcome was confirmed with a Pareeshak test (BHAT Bio-tech, Bangalore, India).
Total fibrinogen concentrations were established with the adapted Clauss method (Multifibrin U-test, BCS coagulation analyzer, Dade Behring, Deerfield, IL, USA). Fibrinogen γ' was quantified with an ELISA using a 2.G2.H9 mouse monoclonal coating antibody against the human γ' sequence (Santa Cruz Biotechnology, Santa Cruz, USA) for antigen capture and a goat polyclonal antibody against human fibrinogen (Antibody 7,539, Abcam for development, Cambridge, USA). Turbidimetric analysis (A405 nm) was used (38) with modified tissue factor and tissue-plasminogen activator (tPA) concentrations to attain plasma fibrinolytic potential reported as CLTs of around 60 to 100 min. Final concentrations were tissue factor (125x diluted – an estimated final concentration of 59 pM; Dade Innovin, Siemens Healthcare Diagnostics Inc., Marburg, Germany), CaCl2 (17 mmol/L), tPA (100 ng/mL; Actilyse, Boehringer, Ingelheim, Germany) and phospholipid vesicles (10 μmol/L; Rossix, Mölndal, Sweden). Kinetics of clot formation (with lag time presenting the time required for fibrin fibers to attain the appropriate length to allow lateral aggregation as well as activation of the coagulation cascade and slope representing the rate of lateral aggregation), and structural clot properties (maximum absorbance representing clot density) were computed from the turbidity curves. The coefficient of variation for all assays was <10%.
Factors Pertaining to Lifestyle
Participants responded to various validated interviewer-administered questionnaires in a language of their choice to obtain demographic and lifestyle information.
Genetic Analyses
For details on the genetic analyses, please refer to Cronjé, Nienaber-Rousseau (39). In short, genomic DNA was isolated from buffy coat using the FlexiGene™ kits (QIAGEN Inc., Valencia, CA, USA). SNPs selected from the literature [FGB-rs7439150, rs1800789 (1420G/A), rs1800791 (−854G/A), rs1800790 (−455G/A), rs1800788 (249C/T), rs1800787 (−148C/T), rs4220, rs4463047, FGA-rs6050, rs2070011 (2224G/A), and FGG rs2066865 and rs1049636 (9340T/C) as well as FXIII His95Arg A/G (rs6003) and Val34Leu, C/A (rs5985)], and two novel FGB SNPs in the promoter region (rs2227385 and rs2227388) identified via sequencing a subpopulation were genotyped. The Thermo Fischer Scientific® Taqman based assay, Illumina® VeraCode GoldenGate assay technology using a BeadXpress® platform and competitive allele-specific polymerase chain reactions (KASP) were used for genotyping.
Statistical Analyses
Statistica® version 13.3 (TIBCO Software Inc., Tulsa, OK, USA) and R version 3.5.0 (40) were used for analysis. Normally distributed data are expressed as mean [95% confidence intervals (CI)], and not-normally distributed data as median (lower and upper quartiles). Differences between HIV and urbanization status in terms of 25(OH)D and coagulation factors were determined by Mann–Whitney U tests and among tobacco use subgroups, using Kruskal–Wallis analysis of variance (ANOVA). Spearman correlations were performed to test for statistical dependence between variables and identify confounders and/or co-variates for inclusion in subsequent analysis (age, LDL-C, body mass index (BMI), HbA1c, GGT, and CRP). Spearman partial correlations were conducted adjusting for these co-variates and additionally for fibrinogen when the outcome variables were clot properties to determine whether the relationship between vitamin D and clot properties was influenced by fibrinogen concentration. The following categories were created based on vitamin D status <50 nmol/L (deficient), 50–75 nmol/L (insufficient) and >75 nmol/L (sufficient) (41, 42). General linear models were used to determine differences in fibrinogen and clot properties among the vitamin D status categories.
To investigate whether vitamin D has interactive effects with the fibrinogen and FXIII gene polymorphisms in predicting total and %γ' fibrinogen concentration and clot properties, general linear models (adjusting for age, LDL-C, BMI, HbA1c, GGT, and CRP) with interactions were adopted. Since the frequency of carrier homozygote variant carrier status was low, we group the genotypes using a dominant genetic mode of action, that is, combining the heterozygotes with those that were homozygous for the variant alleles (0/1/1). Only interactions with a p-value < 0.05 are presented; however, for multiple testing and false discovery rates, the statistical threshold for significance was p < 0.005 [1/(16*3)(0.25)] when the false discovery rate was set at 25% (43), the number of SNPs investigated was 16, the vitamin D one, whereas total and γ' fibrinogen concentrations and fibrin clot properties were regarded as three due to co-linearity. To describe the interaction, Spearman partial correlations were used, adjusting for age, LDL-C, BMI, HbA1c, GGT, and CRP and additionally for fibrinogen, when the interaction was concerning clot properties.
Results
Characteristics of the Study Population in Terms of Vitamin D Status, Fibrinogen, and Clot Properties
In Table 1, we summarize the demographic and clinical characteristics of the women constituting our sample. The median 25(OH)D concentration was 69.5 nmol/L; 16% of the women presented with vitamin D deficiency, 45% with insufficiency and 39% with adequate status.
As can be seen from the table above, 25(OH)D differed between the rural and urban women, with the former presenting with higher concentrations. Total and %γ' fibrinogen, time to clot formation indicated by lag time and the rate of lateral aggregation as indicated by slope also differed. Total fibrinogen and slope were higher in the rural than urban participants, whereas %γ' fibrinogen and lag time were higher in the urban than rural living participants. Clot density as indicated by maximum absorption and CLT did not differ between the urbanization subgroups. No difference in 25(OH)D concentration, fibrinogen, % γ' fibrinogen, lag time, slope, maximum absorbance, or CLT was detected between HIV-seropositive and -negative women.
Association of 25(OH)D Concentrations With Circulating Fibrinogen and Clot Properties
In Table 2, we present the correlations between 25(OH)D and fibrinogen as well as clot properties. Noteworthy observations were total fibrinogen correlating inversely, whereas γ' fibrinogen correlated positively with 25 (OH)D and a negative correlation between maximum absorbance and 25(OH)D. However, after adjusting for variables including CRP (r = 0.42; p < 0.001 with fibrinogen) as well as HbA1c (r = 0.20; p < 0.001 with fibrinogen), only the correlations between γ' fibrinogen and lag time remained and, after further adjustment with fibrinogen, only lag time remained worth mentioning.
Table 2. Correlations of 25(OH)D with total fibrinogen and γ' fibrinogen concentrations and plasma clot properties.
In Table 3, we present the differences in fibrinogen and clot properties among the 25(OH)D status categories. Fibrinogen γ' tended to increase over the 25(OH)D categories (the first represented those with deficiency, the second those with insufficiency and the third those with sufficient vitamin D) before and after adjusting for possible confounders. Maximum absorbance decreased over increasing 25(OH)D categories, with a significant difference detected between those with vitamin D insufficiency and those with sufficient status before and after adjusting for potential confounders, but not after additional adjustment for total fibrinogen. Neither fasting glucose nor HbA1c differed over the vitamin D status categories (p > 0.05).
Table 3. 25(OH)D status categories in relation to total and γ' fibrinogen as well as plasma clot properties.
Gene−25(OH)D Interactions in Relation to Fibrinogen and Plasma Clot Properties
Details of the individual fibrinogen (FGA, FGB and FGG) and FXIII genes SNPs and their associations with total and %γ' fibrinogen concentration and clot properties have been described elsewhere (39). Here, we focus on the conjoined associations between 25(OH)D and genetics of fibrinogen and FXIII with total fibrinogen, %γ' fibrinogen and clot properties. In Table 4, we report only interactions at p < 0.05 in models wherein we controlled for other variables.
Table 4. FGB and FXIII gene polymorphisms' interactions with 25(OH)D in relation to plasma fibrin clot properties.
In those carrying the FGB wildtype allele at rs1800787, a stronger correlation was detected between 25(OH)D and total fibrinogen than those harboring the variant allele, resulting in an interaction. 25(OH)D did not correlate with lag time in those with homozygote carrier status for the wildtype allele. In contrast, when harboring the variant allele at rs5985 on the FXIII gene, a positive association was detected that resulted in an interaction. However, after additional adjustment for fibrinogen, the gene-diet interaction disappeared. An FGB gene (rs2227388)–vitamin D interaction was observed concerning maximum absorbance that became stronger after additional adjustment for fibrinogen. In women homozygous for the wildtype allele at rs2227388, no association between 25(OH)D and maximum absorbance was evident. Contrastingly, in those with one or two variant alleles, an inverse correlation was detected. 25(OH)D interacted with FGA: rs2227385 in relation to CLT before and after additional adjustment for fibrinogen. However, none of the interactions reached the significance required for adjusting for a multiple testing false discovery rate.
Discussion
The role of vitamin D in reducing CVD risk has been recognized (8). In particular, among the favorable pleiotropic effects of vitamin D on CVD, researchers have explored the relationship between vitamin D and total fibrinogen concentrations in small patient populations with large scale epidemiological studies remaining relatively limited. In this study, we determined the relationships between total fibrinogen and—to our knowledge for the first time—the γ' splice variant of fibrinogen, and clot properties and 25(OH)D concentrations in generally healthy, Black South African women, among whom vitamin D deficiency was common. We observed an inverse, albeit weak, relationship between 25(OH)D and fibrinogen as well as clot density as indicated by maximum absorbance and a positive association between 25(OH)D and γ' fibrinogen. This was supported by observing differences in maximum absorbance and γ' fibrinogen among vitamin D status groups as well. Partial correlations with adjustment for, among others, CRP, an inflammatory marker, and HbA1c, which is a long-term indicator of glucose control, indicated that vitamin D only correlated with γ' fibrinogen and lag time. Even though several correlations were statistically significant, possibly due to a loss in statistical power, less differences among fibrinogen and clot properties were observed when the study population was subdivided according to their vitamin D status. Several gene-diet interactions were detected at the p < 0.05 level, but did not stand after correction for false discovery rate.
Here we also report an inverse association between 25(OH)D and fibrinogen before adjusting for other variables. Our study and the results obtained are in line with those of Parikh et al. (25) obtained in 701 Black and White adolescents living in the southeastern region of America; of Mellenthin et al. (19), who investigated 2,723 men and women from a general European population; of Lo Gullo et al. (26), who sampled 27 rheumatoid arthritis patients and 41 healthy controls; of Blondon et al. (12) who studied 6,554 multiethnic American adults free of CVD; and of de Oliveira et al. (20), who investigated 5,870 English adults. Mellenthin et al. (19), Blondon et al. (12) and de Oliveira et al. (20) are, to our knowledge, the only investigations similar to ours, that are, of observational nature involving adult participants from general populations. The observational study of Mellenthin et al. (17) did not include models adjusting for co-variates. Oliveira (18) reported that the relationship between fibrinogen and vitamin D remained whereas Blondon et al. (12) found that it did not after adjusting for age and other demographic variables. Evidence for a relationship between vitamin D and fibrinogen from intervention trials predominantly involved small groups of ill patients, with conflicting results. Even though Akin et al. (21), Fornari et al. (27) and Barnes et al. (44) reported measuring both vitamin D and fibrinogen concentration, they did not report on any association between the variables. An intervention trial conducted in Helsinki, Finland, administered 25(OH)D supplements for 6 months to 218 long-term care patients who were chronically bedridden, but failed to alter fibrinogen concentrations (22). A double-blind, randomized, controlled trial conducted over a 12-week period in 90 type 2 diabetic men and women was used to investigate the effects of vitamin D either with or without extra calcium on fibrinogen (23). Compared to the baseline values, but not among the subgroups, fibrinogen decreased in both those consuming a vitamin D-fortified Persian yogurt drink (or doogh, containing 500 IU vitamin D3 and 150 mg calcium/250 mL) and calcium plus vitamin D3-fortified doogh (containing 500 IU vitamin D3 and 250 mg calcium/250 mL). The reduction was seen in both groups that also received calcium; future research is therefore needed to tease out the effects of vitamin D supplementation per se. Soare et al. (24) failed to detect a difference in fibrinogen between vitamin D-deficient (25(OH)D >20 ng/mL) and -sufficient groups (≥20 ng/mL) in 62 patients with inflammatory bowel disease. Asymptomatic or mildly symptomatic severe acute respiratory syndrome coronavirus 2 RNA positive vitamin D deficient participants were randomized to receive cholecalciferol for a week to correct the deficiency (n = 16) or a placebo (n = 24) (45). Vitamin D supplementation decreased fibrinogen levels, but did not affect other inflammatory markers (45). To summarize, both observational and experimental studies' results allude to an inverse relationship between fibrinogen and vitamin D which is in line with our findings.
Regarding γ' fibrinogen and clot properties, our study is the first to allude to a relationship with vitamin D. Here we expressed γ' as a percentage of total fibrinogen concentrations. In individuals with varying total fibrinogen concentration, the γ' fibrinogen concentration will likely also vary as its concentration is dependent on the total amount of fibrinogen produced. When expressing it as % of total fibrinogen, the potential confounding effect of differing total fibrinogen concentrations is excluded. Previously we reported our observations that while γ' was high when total fibrinogen concentration was high; a decreased γ'/total fibrinogen ratio was observed at the highest total fibrinogen concentrations (32). This finding explains the seemingly conflicting relationships i.e., the positive relationship detected between vitamin D and γ' and the negative relationship between vitamin D and total fibrinogen.
Mechanisms by which vitamin D may influence fibrinogen and clot properties are not known and should be explored in future research to determine the biological relevance of our observations. Possible mechanisms for the inverse association between vitamin D and fibrinogen and clot density might be through inflammation and changes in the glycemic milieu. Thus, regarding inflammation, fibrinogen is itself a downstream inflammatory molecule and therefore directly affected by the suppressor effect of vitamin D on the production of cytokines. In vitro evidence demonstrated that the active metabolite of vitamin D, 1, 25(OH)D, and its analogs could couple to the vitamin D receptor and subsequently downregulate several inflammatory cytokines, including tumor necrosis factor-α, and interleukin 1 and 6 (46–49). This theory is partly supported by our results, because after adjustment for the inflammatory marker CRP (which correlated strongly with fibrinogen), the relationship between fibrinogen and vitamin D disappeared. Moreover, vitamin D's effects on insulin and fasting blood glucose could improve plasma clot properties. Heightened insulin sensitivity (16, 17) could change the glycation of the activated fibrinogen molecules, thereby improving clot structure (making it less dense and therefore less resistant to degradation by plasmin) (50, 51). However, we did not detect any difference in HbA1c or fasting glucose among the vitamin D status groups. Furthermore, the associations between fibrinogen and clot properties remained after adjustment for age and HbA1c alone (results not shown), possibly because the participants did not have diabetes. Therefore, the HbA1c concentrations were normal for most in our study population, so adjusting for it had little effect. To explore the latter mechanism in future, one would probably need to sample diabetic patients and not apparently healthy individuals with normal glycation levels. Whatever the mechanism, decreased plasma total and γ' fibrinogen and improved plasma clot properties (i.e., here vitamin D was associated with less dense plasma clot configurations) may have implications for preventing CVD.
Currently, there is no proof of a direct pathway by which 25(OH)D or its metabolized product, 1, 25(OH)D, influences the expression of the fibrinogen genes, FGA, FGB, and FGG. It seems more plausible that fibrinogen gene expression is influenced by cytokines (36); because they are downregulated by vitamin D (46–49), this might modulate the genes. In terms of FXIII's genetic determinants and possible modulators, apart from the work from our research group, we could not find any other publications. Even though we detected some vitamin D–gene interactions, we question their significance because they did not remain after correction for false discovery rate. Therefore, the value of our findings in terms of nutrigenetics using a candidate gene approach remains to be clarified in studies with more power or through the use of genome wide association studies investigated interactions with vitamin D.
Some limitations that could be addressed by future research must be acknowledged. As women are more likely to develop skeletal disorders associated with low 25(OH)D status than men and the study budget being constrained, we only included women. Our study's power to detect interactions with SNPs was low because we had access to a relatively small subset. Its strengths include that we investigated an ostensibly healthy population of Black South Africans known to present with high fibrinogen and low vitamin D owing to skin pigmentation, thereby increasing our chances of detecting an association should it exist. While total fibrinogen concentrations here reflect functional fibrinogen, the γ' fibrinogen levels are antigen based. Even though the absolute percentage may be slightly different if one measured total fibrinogen with an antigen assay, all samples were analyzed with the same methods, thus allowing comparison between individuals with low or high γ' fibrinogen. We, furthermore, report % γ' fibriniogen for correlations and relationships since the ratio is thought to be more physiologically relevant than the absolute concentration. Currently, our report is the only one of its kind to have considered the γ' fibrinogen molecule and plasma clot properties.
Conclusions
Low vitamin D (7, 52, 53) and elevated total (54–62) and γ' fibrinogen (63–67) as well as dense blood clots resistant to lysis (68, 69) are independently associated with unfavorable health outcomes, including CVD. The current study expanded on existing knowledge of the putative role of vitamin D in terms of fibrinogen and the functionalities of fibrinogen in terms of fibrin clot properties, many of which have not been studied before with respect to vitamin D. In the present cross-sectional study, we exposed an inverse association between 25(OH)D and fibrinogen concentrations as well as clot density. Moreover, γ' fibrinogen correlated positively with 25(OH)D. Because the association between fibrinogen and clot properties failed to remain statistically significant after adjustment with variables including CRP, we believe that plausible mechanisms that could be explored in future are through the anti-inflammatory properties of vitamin D. γ' fibrinogen increased whereas clot density decreased over the 25(OH)D categories before and after adjustment for confounders. Because the observed relationships were weak, the clinical significance of the findings should be determined in future epidemiological and intervention trials. The results reported here could inform future intervention trials in which responsible exposure to sunlight, the inclusion of vitamin D-rich foods and/or whether supplementation or fortification with vitamin D analogs are explored in relation to blood coagulation might ultimately offer therapeutic options for those with a high risk of developing heart disease. The observed relationship between vitamin D and CVD could ultimately reveal promising mechanisms—that might include a link to fibrinogen and plasma clot configurations—which could be modified by lifestyle or supplementation and fortification.
Data Availability Statement
The data analyzed in this study is subject to the following licenses/restrictions: Raw data were generated at the North-West University. Obtained data supporting the findings of this study are available from the PI of the study after obtaining the necessary ethical approval. Due to Ethical restrictions and the POPI act data are not publicly available. Requests to access these datasets should be directed to lanthe.kruger@nwu.ac.za.
Ethics Statement
The studies involving human participants were reviewed and approved by Health Research Ethics Committee of the Faculty of Health Sciences, North-West University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
Conceptualization: CN-R. Methodology: PR, CN-R, ZL-L, IK, and MP. Statistical analysis and writing—original draft preparation: PR and CN-R. Resources and funding acquisition: IK and MP. Data curation and project administration: IK. Writing—review and editing: ZL-L, IK, and MP. All authors contributed to the article and approved the submitted version.
Funding
This research was funded by SANPAD [South Africa–The Netherlands Research Program on Alternatives in Development (08/15)], the South African National Research Foundation (NRF GUN No. 2069139 and FA2006040700010), self-initiated research grants from the South African Medical Research Council, North-West University, Potchefstroom, South Africa, the Population Health Research Institute, Ontario, Canada and the Academy of Medical Sciences UK (Newton Fund Advanced Fellowship Grant [AMS-NAF1-Pieters to MP].
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
PURE Africa: Annamarie Kruger (posthumous), Minrie Greeff, Poncho G. Molaudzi, Vasti Kruger, Juanita Brits, and Paula Jardim working in the Africa Unit for Transdisciplinary Health Research (AUTHeR), Faculty of Health Sciences, North-West University, Potchefstroom, South Africa, as well as all the field workers over the years. PURE International: Salim Yusuf and the PURE project office staff at the Population Health Research Institute (PHRI), Hamilton Health Sciences, and McMaster University, Ontario, Canada. Analyses and visualizations used in this paper were produced with the Giovanni online data system, developed and maintained by the NASA GES DISC.
References
1. MEMBERS WG, Roger VL, Go AS, Lloyd-Jones DM, Benjamin EJ, Berry JD, et al. Heart disease and stroke statistics−2012 update: a report from the American Heart Association. Circulation. (2012) 125:e2. doi: 10.1161/CIR.0b013e31823ac046
2. Bennett JE, Stevens GA, Mathers CD, Bonita R, Rehm J, Kruk ME, et al. NCD countdown 2030: worldwide trends in non-communicable disease mortality and progress towards sustainable development goal target 3.4. Lancet. (2018) 392:1072–88. doi: 10.1016/S0140-6736(18)31992-5
3. Mehra MR, Desai SS, Kuy S, Henry TD, Patel AN. Cardiovascular disease, drug therapy, and mortality in COVID-19. N Engl J Med. (2020) 382:e102. doi: 10.1056/NEJMoa2007621
4. Ruan Q, Yang K, Wang W, Jiang L, Song J. Clinical predictors of mortality due to COVID-19 based on an analysis of data of 150 patients from Wuhan, China. Intensive Care Med. (2020) 46:846–8. doi: 10.1007/s00134-020-05991-x
5. Shi S, Qin M, Shen B, Cai Y, Liu T, Yang F, et al. Association of cardiac injury with mortality in hospitalized patients with COVID-19 in Wuhan, China. JAMA Cardiol. (2020) 5:802–10. doi: 10.1001/jamacardio.2020.0950
6. Oberoi D, Mehrotra V, Rawat A. “Vitamin D” as a profile marker for cardiovascular diseases. Ann Card Anaesth. (2019) 22:47. doi: 10.4103/aca.ACA_66_18
7. Chowdhury R, Kunutsor S, Vitezova A, Oliver-Williams C, Chowdhury S, Kiefte-de-Jong JC, et al. Vitamin D and risk of cause specific death: systematic review and meta-analysis of observational cohort and randomised intervention studies. BMJ. (2014) 348:g1903. doi: 10.1136/bmj.g1903
8. Wang L, Song Y, Manson JE, Pilz S, März W, Michaëlsson K, et al. Circulating 25-hydroxy-vitamin D and risk of cardiovascular disease: a meta-analysis of prospective studies. Circ Cardiovasc Qual Outcomes. (2012) 5:819–29. doi: 10.1161/CIRCOUTCOMES.112.967604
9. Myburgh PH, Towers GW, Kruger IM, Nienaber-Rousseau C. CRP genotypes predict increased risk to co-present with low vitamin D and elevated CRP in a group of healthy black South African women. Int J Environ Res Public Health. (2018) 15:111. doi: 10.3390/ijerph15010111
10. Liefaard MC, Ligthart S, Vitezova A, Hofman A, Uitterlinden AG., Kiefte-de Jong JC, et al. Vitamin D and C-reactive protein: a Mendelian randomization study. PloS ONE. (2015) 10:e0131740. doi: 10.1371/journal.pone.0131740
11. Carlson N, Mah R, Aburto M, Peters MJ, Dupper MV, Chen LH. Hypovitaminosis D correction and high-sensitivity C-reactive protein levels in hypertensive adults. Perm J. (2013) 17:19. doi: 10.7812/TPP/13.963
12. Blondon M, Cushman M, Jenny N, Michos ED, Smith NL, Kestenbaum B, et al. Associations of serum 25-hydroxyvitamin D with hemostatic and inflammatory biomarkers in the multi-ethnic study of atherosclerosis. J Clin Endocrinol Metab. (2016) 101:2348–57. doi: 10.1210/jc.2016-1368
13. Schnatz PF, Jiang X, Vila-Wright S, Aragaki AK, Nudy M, O'Sullivan DM, et al. Calcium/vitamin D supplementation, serum 25-hydroxyvitamin D concentrations, and cholesterol profiles in the women's health initiative Calcium/vitamin D randomized trial. Menopause. (2014) 21:823–33. doi: 10.1097/GME.0000000000000188
14. Li YC, Qiao G, Uskokovic M, Xiang W, Zheng W, Kong J. Vitamin D: a negative endocrine regulator of the renin–angiotensin system and blood pressure. J Steroid Biochem Mol Biol. (2004) 89:387–92. doi: 10.1016/j.jsbmb.2004.03.004
15. Govender D, Damjanovic L, Gaza CA, Meyer V. Vitamin D decreases silencer methylation to downregulate renin gene expression. Gene. (2021) 786:145623. doi: 10.1016/j.gene.2021.145623
16. Norman AW, Frankel J, Heldt AM, Grodsky GM. Vitamin D deficiency inhibits pancreatic secretion of insulin. Science. (1980) 209:823–5. doi: 10.1126/science.6250216
17. Rammos G, Tseke P, Ziakka S. Vitamin D, the renin-angiotensin system, and insulin resistance. Int Urol Nephrol. (2008) 40:419–26. doi: 10.1007/s11255-007-9244-4
18. Carthy E, Yamashita W, Hsu A, Ooi B. 1, 25-Dihydroxyvitamin D3 and rat vascular smooth muscle cell growth. Hypertension. (1989) 13:954–9. doi: 10.1161/01.HYP.13.6.954
19. Mellenthin L, Wallaschofski H, Grotevendt A, Völzke H, Nauck M, Hannemann A. Association between serum vitamin D concentrations and inflammatory markers in the general adult population. Metabolism. (2014) 63:1056–62. doi: 10.1016/j.metabol.2014.05.002
20. de Oliveira C, Biddulph JP, Hirani V, Schneider IJC. Vitamin D and inflammatory markers: cross-sectional analyses using data from the English Longitudinal Study of Ageing (ELSA). J Nutr Sci. (2017) 6:e1. doi: 10.1017/jns.2016.37
21. Akin F, Ayça B, Köse N, Duran M, Sari M, Uysal OK, et al. Serum vitamin D levels are independently associated with severity of coronary artery disease. J Invest Med. (2012) 60:869–73. doi: 10.2310/JIM.0b013e31825457cb
22. Bjorkman M, Sorva A, Tilvis R. C-reactive protein and fibrinogen of bedridden older patients in a six-month vitamin D supplementation trial. J Nutr Health Aging. (2009) 13:435–9. doi: 10.1007/s12603-009-0080-3
23. Neyestani TR, Nikooyeh B, Alavi-Majd H, Shariatzadeh N, Kalayi A, Tayebinejad N, et al. Improvement of vitamin D status via daily intake of fortified yogurt drink either with or without extra calcium ameliorates systemic inflammatory biomarkers, including adipokines, in the subjects with type 2 diabetes. J Clin Endocrinol Metab. (2012) 97:2005–11. doi: 10.1210/jc.2011-3465
24. Soare I, Sirbu A, Diculescu M, Mateescu BR, Tieranu C, Cima L, . editors. Vitamin D deficiency and inflammation in IBD patients. Endocrine Abstracts. (2021) 77:34. doi: 10.1530/endoabs.77.P34
25. Parikh S., Guo D-h, Pollock NK, Petty K, Bhagatwala J, Gutin B, et al. Circulating 25-hydroxyvitamin D concentrations are correlated with cardiometabolic risk among American black and white adolescents living in a year-round sunny climate. Diabetes Care. (2012) 35:1133–8. doi: 10.2337/dc11-1944
26. Lo Gullo A, Mandraffino G, Bagnato G, Aragona CO, Imbalzano E, D'Ascola A, et al. Vitamin D status in rheumatoid arthritis: inflammation, arterial stiffness and circulating progenitor cell number. PLoS ONE. (2015) 10:e0134602. doi: 10.1371/journal.pone.0134602
27. Fornari R, Francomano D, Greco E, Marocco C, Lubrano C, Wannenes F, et al. Lean mass in obese adult subjects correlates with higher levels of vitamin D, insulin sensitivity and lower inflammation. J Endocrinol Invest. (2015) 38:367–72. doi: 10.1007/s40618-014-0189-z
28. De Moerloose P, Boehlen F, Neerman-Arbez M. Fibrinogen and the risk of thrombosis. Semin Thromb Hemost. (2010) 36:7–17. doi: 10.1055/s-0030-1248720
29. Blombäck B. Fibrinogen and fibrin-proteins with complex roles in hemostasis and thrombosis. Thromb Res. (1996) 83:1–75. doi: 10.1016/0049-3848(96)00111-9
30. Herrick S, Blanc-Brude O, Gray A, Laurent G. Fibrinogen. Int J Biochem Cell Biol. (1999) 31:741–6. doi: 10.1016/S1357-2725(99)00032-1
31. Macrae FL, Swieringa F, Heemskerk JW, Ariens RA. High fibrinogen γ' levels in patient plasma increase clot formation at arterial and venous shear. Blood Adv. (2021) 5:3468–77. doi: 10.1182/bloodadvances.2020003346
32. Pieters M, Kotze RC, Jerling JC, Kruger A, Ariëns RA. Evidence that fibrinogen γ′ regulates plasma clot structure and lysis and relationship to cardiovascular risk factors in black Africans. Blood. (2013) 121:3254–60. doi: 10.1182/blood-2012-12-471482
33. Cooper AV, Standeven KF, Ariëns RA. Fibrinogen gamma-chain splice variant γ′ alters fibrin formation and structure. Blood. (2003) 102:535–40. doi: 10.1182/blood-2002-10-3150
34. Ariëns RA, Philippou H, Nagaswami C, Weisel JW, Lane DA, Grant PJ. The factor XIII V34L polymorphism accelerates thrombin activation of factor XIII and affects cross-linked fibrin structure. Blood. (2000) 96:988–95. doi: 10.1182/blood.V96.3.988.015k57_988_995
35. Lim BC, Ariëns RA, Carter AM, Weisel JW, Grant PJ. Genetic regulation of fibrin structure and function: complex gene-environment interactions may modulate vascular risk. Lancet. (2003) 361:1424–31. doi: 10.1016/S0140-6736(03)13135-2
36. Cronje HT, Nienaber-Rousseau C, Zandberg L, De Lange Z, Green FR, Pieters M. Fibrinogen and clot-related phenotypes determined by fibrinogen polymorphisms: Independent and IL-6-interactive associations. PLoS ONE. (2017) 12:e0187712. doi: 10.1371/journal.pone.0187712
37. Teo K, Chow CK, Vaz M, Rangarajan S, Yusuf S. The Prospective Urban Rural Epidemiology (PURE) study: examining the impact of societal influences on chronic noncommunicable diseases in low-, middle-, and high-income countries. Am Heart J. (2009) 158:1–7.e1. doi: 10.1016/j.ahj.2009.04.019
38. Lisman T, de Groot PG, Meijers JC, Rosendaal FR. Reduced plasma fibrinolytic potential is a risk factor for venous thrombosis. Blood. (2005) 105:1102–5. doi: 10.1182/blood-2004-08-3253
39. Cronjé HT, Nienaber-Rousseau C, Zandberg L, Chikowore T, de Lange Z, van Zyl T, et al. Candidate gene analysis of the fibrinogen phenotype reveals the importance of polygenic co-regulation. Matrix Biol. (2017) 60–61:16–26. doi: 10.1016/j.matbio.2016.10.005
40. R Core Team. R: A Language Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria (2018). Available online at https://www.R-project.org/
41. Holick MF, Binkley NC, Bischoff-Ferrari HA, Gordon CM, Hanley DA, Heaney RP, et al. Evaluation, treatment, and prevention of vitamin D deficiency: an endocrine society clinical practice guideline. J Clin Endocrinol Metab. (2011) 96:1911–30. doi: 10.1210/jc.2011-0385
42. Bischoff-Ferrari HA, Giovannucci E, Willett WC, Dietrich T, Dawson-Hughes B. Estimation of optimal serum concentrations of 25-hydroxyvitamin D for multiple health outcomes. Am J Clin Nutr. (2006) 84:18–28. doi: 10.1093/ajcn/84.1.18
43. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B. (1995) 57:289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
44. Barnes MS, Horigan G, Cashman KD, Hill TR, Forsythe LK, Lucey AJ, et al. Maintenance of wintertime vitamin D status with cholecalciferol supplementation is not associated with alterations in serum cytokine concentrations among apparently healthy younger or older adults. J Nutr. (2011) 141:476–81. doi: 10.3945/jn.110.131516
45. Rastogi A, Bhansali A, Khare N, Suri V, Yaddanapudi N, Sachdeva N, et al. Short term, high-dose vitamin D supplementation for COVID-19 disease: a randomised, placebo-controlled, study (SHADE study). Postgrad Med J. (2022) 98:87–90. doi: 10.1136/postgradmedj-2020-139065
46. Nagpal S, Na S, Rathnachalam R. Non-calcemic actions of vitamin D receptor ligands. Endocr Rev. (2005) 26:662–87. doi: 10.1210/er.2004-0002
47. Khoo AL, Chai LY, Koenen HJ, Sweep FC, Joosten I, Netea MG, et al. Regulation of cytokine responses by seasonality of vitamin D status in healthy individuals. Clin Exp Immunol. (2011) 164:72–9. doi: 10.1111/j.1365-2249.2010.04315.x
48. Panichi V, De Pietro S, Andreini B, Bianchi AM, Migliori M, Taccola D, et al. Calcitriol modulates in vivo and in vitro cytokine production: a role for intracellular calcium. Kidney Int. (1998) 54:1463–9. doi: 10.1046/j.1523-1755.1998.00152.x
49. Dickie LJ, Church LD, Coulthard LR, Mathews RJ, Emery P, McDermott MF. Vitamin D3 down-regulates intracellular toll-like receptor 9 expression and toll-like receptor 9-induced IL-6 production in human monocytes. Rheumatology. (2010) 49:1466–71. doi: 10.1093/rheumatology/keq124
50. Carr ME. Diabetes mellitus: a hypercoagulable state. J Diabetes Complications. (2001) 15:44–54. doi: 10.1016/S1056-8727(00)00132-X
51. Pieters M, Covic N, van der Westhuizen FH, Nagaswami C, Baras Y, Jerling JC, et al. Glycaemic control improves fibrin network characteristics in type 2 diabetes–a purified fibrinogen model. Thromb Haemost. (2008) 99:691–700. doi: 10.1160/TH07-11-0699
52. Brøndum-Jacobsen P, Benn M, Jensen GB, Nordestgaard BG. 25-hydroxyvitamin d levels and risk of ischemic heart disease, myocardial infarction, and early death: population-based study and meta-analyses of 18 and 17 studies. Arterioscler Thromb Vasc Biol. (2012) 32:2794–802. doi: 10.1161/ATVBAHA.112.248039
53. Brøndum-Jacobsen P, Nordestgaard BG, Schnohr P, Benn M. 25-hydroxyvitamin D and symptomatic ischemic stroke: an original study and meta-analysis. Ann Neurol. (2013) 73:38–47. doi: 10.1002/ana.23738
54. Koenig W. Fibrin (ogen) in cardiovascular disease: an update. Thromb Haemost. (2003) 89:601–9. doi: 10.1055/s-0037-1613566
55. Danesh J, Collins R, Appleby P, Peto R. Association of fibrinogen, C-reactive protein, albumin, or leukocyte count with coronary heart disease: meta-analyses of prospective studies. JAMA. (1998) 279:1477–82. doi: 10.1001/jama.279.18.1477
56. Danesh J, Lewington S, Thompson SG, Lowe G, Collins R, Kostis J, et al. Plasma fibrinogen level and the risk of major cardiovascular diseases and nonvascular mortality: an individual participant meta-analysis. JAMA. (2005) 294:1799–809. doi: 10.1001/jama.294.14.1799
57. Meade T, Brozovic M, Chakrabarti R, Haines A, Imeson J, Mellows S, et al. Haemostatic function and ischaemic heart disease: principal results of the Northwick Park heart study. Lancet. (1986) 328:533–7. doi: 10.1016/S0140-6736(86)90111-X
58. Wilhelmsen L, Svärdsudd K, Korsan-Bengtsen K, Larsson B, Welin L, Tibblin G. Fibrinogen as a risk factor for stroke and myocardial infarction. N Engl J Med. (1984) 311:501–5. doi: 10.1056/NEJM198408233110804
59. Ernst E. The role of fibrinogen as a cardiovascular risk factor. Atherosclerosis. (1993) 100:1–12. doi: 10.1016/0021-9150(93)90062-Y
60. Maresca G, Di Blasio A, Marchioli R, Di Minno G. Measuring plasma fibrinogen to predict stroke and myocardial infarction: an update. Arterioscler Thromb Vasc Biol. (1999) 19:1368–77. doi: 10.1161/01.ATV.19.6.1368
61. Stone M, Thorp J. Plasma fibrinogen—a major coronary risk factor. J R Coll Gen Pract. (1985) 35:565–9.
62. Kannel WB, D'Agostino RB, Belanger AJ. Fibrinogen, cigarette smoking, and risk of cardiovascular disease: insights from the Framingham Study. Am Heart J. (1987) 113:1006–10. doi: 10.1016/0002-8703(87)90063-9
63. Lovely RS, Falls LA, Al-Mondhiry HA, Chambers CE, Sexton GJ Ni H, et al. Association of γA/γ'fibrinogen levels and coronary artery disease. Thromb Haemost. (2002) 88:26–31. doi: 10.1055/s-0037-1613148
64. Mannila MN, Lovely RS, Kazmierczak S, Eriksson P, Samnegård A, Farrell D, et al. Elevated plasma fibrinogen γ′ concentration is associated with myocardial infarction: effects of variation in fibrinogen genes and environmental factors. J Thromb Haemost. (2007) 5:766–73. doi: 10.1111/j.1538-7836.2007.02406.x
65. Cheung EY, de Willige SU, Vos HL, Leebeek FW, Dippel DW, Bertina RM, et al. Fibrinogen γ′ in ischemic stroke: a case-control study. Stroke. (2008) 39:1033–5. doi: 10.1161/STROKEAHA.107.495499
66. Cheung EY, Vos HL, Kruip MJ, den Hertog HM, Jukema JW, de Maat MP. Elevated fibrinogen γ′ ratio is associated with cardiovascular diseases and acute phase reaction but not with clinical outcome. Blood. (2009) 114:4603–4. doi: 10.1182/blood-2009-08-236240
67. Lovely RS, Yang Q, Massaro JM, Wang J, D'Agostino Sr RB, O'Donnell CJ, et al. Assessment of genetic determinants of the association of γ′ fibrinogen in relation to cardiovascular disease. Arterioscler Thromb Vasc Biol. (2011) 31:2345–52. doi: 10.1161/ATVBAHA.111.232710
68. Ariëns RAS. Fibrin(ogen) and thrombotic disease. J Thromb Haemost. (2013) 11:294–305. doi: 10.1111/jth.12229
69. Undas A, Szułdrzynski K, Stepien E, Zalewski J, Godlewski J, Tracz W, et al. Reduced clot permeability and susceptibility to lysis in patients with acute coronary syndrome: effects of inflammation and oxidative stress. Atherosclerosis. (2008) 196:551–7. doi: 10.1016/j.atherosclerosis.2007.05.028
Keywords: cholecalciferol, ergocalciferol, fiber density, fibrinolysis, lag time, lateral aggregation, nutrigenetics, turbidity
Citation: Rautenbach PH, Nienaber-Rousseau C, de Lange-Loots Z, Kruger IM and Pieters M (2022) Associations Between 25-Hydroxyvitamin D and Total and γ' Fibrinogen and Plasma Clot Properties and Gene Interactions in a Group of Healthy Black South African Women. Front. Cardiovasc. Med. 9:868542. doi: 10.3389/fcvm.2022.868542
Received: 02 February 2022; Accepted: 23 June 2022;
Published: 12 July 2022.
Edited by:
Egidio Imbalzano, University of Messina, ItalyReviewed by:
Zsuzsa Bagoly, University of Debrecen, HungaryMagdolna Nagy, Maastricht University, Netherlands
Antonino Catalano, University of Messina, Italy
Copyright © 2022 Rautenbach, Nienaber-Rousseau, de Lange-Loots, Kruger and Pieters. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Cornelie Nienaber-Rousseau, cornelie.nienaber@nwu.ac.za
†These authors have contributed equally to this work and share first authorship