- Research article
- Open Access
Vitamin D pathway gene polymorphisms, diet, and risk of postmenopausal breast cancer: a nested case-control study
Breast Cancer Researchvolume 9, Article number: R9 (2007)
Vitamin D receptor (VDR) polymorphisms have been inconsistently associated with breast cancer risk. Whether risk is influenced by polymorphisms in other vitamin D metabolism genes and whether calcium or vitamin D intake modifies risk by genotype have not been evaluated.
We conducted a nested case-control study within the Cancer Prevention Study II Nutrition Cohort of associations between breast cancer and four VDR single-nucleotide polymorphisms (SNPs), Bsm1,Apa1,Taq1, and Fok1, a poly(A) microsatellite, and associated haplotypes (baTL and BAtS). We also examined one SNP in the 24-hydroxylase gene (CYP24A1) and two in the vitamin D-binding protein (group-specific component [GC]) gene. Participants completed a questionnaire on diet and medical history at baseline in 1992. This study includes 500 postmenopausal breast cancer cases and 500 controls matched by age, race/ethnicity, and date of blood collection.
Incident breast cancer was not associated with any genotype examined. However, women with the Bsm1 bb SNP who consumed greater than the median intake of total calcium (≥902 mg/day) had lower odds of breast cancer compared to women with the Bb or BB genotype and less than the median calcium intake (odds ratio 0.61, 95% confidence interval 0.38 to 0.96; p interaction = 0.01). Similar interactions were observed for Taq1 (T allele) and the poly(A) (LL) repeat.
We found no overall association between selected vitamin D pathway genes and postmenopausal breast cancer risk. However, certain VDR gene polymorphisms were associated with lower risk in women consuming high levels of calcium, suggesting that dietary factors may modify associations by VDR genotype.
Vitamin D has been associated with a lower risk of several types of cancer, including breast cancer [1, 2]. Its action is mediated through the vitamin D receptor (VDR), a nuclear transcription-regulating factor that signals the synthesis of proteins involved in bone mineral homeostasis and cell-cycle regulation . The receptor is present in most cell types, including normal and neoplastic breast tissue . Clinically, breast tumors with greater VDR expression relapse more slowly after first diagnosis .
The VDR gene contains several polymorphisms. Three single-nucleotide polymorphisms (SNPs) located near the 3' region have been examined in relation to breast cancer . These polymorphisms, identified by their restriction endonuclease cleavage sites (Taq1,Bsm1, and Apa1), include an intronic G→A replacement (Bsm1, rs1544410, G = b), an intronic C→A replacement (Apa1, rs7975232, C = a), and a T→C replacement in exon 9 (Taq1, rs731236, C = t), which result in a synonymous change at codon 352. Although not functional, these SNPs are strongly linked with a poly(A) microsatellite repeat in the 3' untranslated region [5, 6] which may influence VDR mRNA stability. Five [7–11] of nine [7–15] studies reported higher breast cancer risk associated with the Bsm1 bb genotype. Other VDR 3' SNPs or haplotypes have been inconsistently related to breast cancer .
The presence of the Fok1 f allele in the 5' promoter region of the VDR (rs2228570) results in production of a VDR protein that is three amino acids longer and less effective as a transcriptional activator . The ff genotype was associated with a statistically significant 34% higher breast cancer risk in the Nurses' Health Study . Previous smaller studies of Fok1 and breast cancer risk were null [7, 13, 17].
Gene-environment interactions may explain in part the inconsistencies in the association between VDR SNPs and breast cancer across studies. Dietary factors are known to impact the vitamin D endocrine system , and recent studies suggest that diet may also influence autocrine/paracrine vitamin D metabolism . Few studies have examined interactions between dietary factors and VDR polymorphisms in relation to breast cancer risk. Furthermore, polymorphisms in other genes in the vitamin D metabolic pathway have not been examined in relation to breast cancer risk.
We examined several polymorphisms in the VDR gene, their associated haplotypes, and polymorphisms in the vitamin D-binding (group-specific component [GC]) protein and 1,25(OH)2D degradation enzyme (CYP24A1) genes in relation to breast cancer risk in a nested case-control study of postmenopausal women enrolled in a large prospective cohort. Of particular a priori interest was whether calcium and vitamin D intakes modified associations by genotype.
Materials and methods
Women in this analysis were participants in the Cancer Prevention Study (CPS) II Nutrition Cohort, a prospective study of cancer incidence among men and women from 21 states in the US, established by the American Cancer Society in 1992 . At enrollment, 86,404 men and 97,788 women completed a 10-page, self-administered questionnaire on medical, dietary, reproductive, and lifestyle variables. Beginning in 1997, follow-up questionnaires were sent to cohort members every 2 years to update exposure information and to ascertain newly diagnosed tumors.
From June 1998 through June 2001, blood specimens were collected from 21,965 Nutrition Cohort postmenopausal women. Using blood samples, we identified 502 breast cancer cases who were cancer-free in 1992. For each case, we randomly selected one postmenopausal female control who had provided a blood sample and was cancer-free (except non-melanoma skin cancer) at the time of the case's diagnosis . Controls were individually matched to cases on birth date (± 6 months), race/ethnicity (Caucasian, African-American, Hispanic, Asian, other/unknown), and date of blood collection (± 6 months) using risk-set sampling. After excluding 2 cases with no match, 500 cases and 500 controls were available for analysis. All aspects of the CPS II cohort were approved by the Emory University (Atlanta, GA, USA) Institutional Review Board.
Genotyping was performed at Applied Biosystems (ABI, Foster City, CA, USA) by laboratory personnel blinded to case-control status; blinded quality control samples for 10% of participants were inserted to validate genotyping procedures. TaqMan (ABI, Foster City, CA, USA) was used to distinguish alleles of the VDR gene at restriction enzyme sites Fok1, Bsm1, Apa1, and Taq1. A variable-length poly(A) sequence in the 3' untranslated region was assessed using methods previously described  and was classified as short (S, A14-A17) or long (L, A18-A22).
We also genotyped selected SNPs in vitamin D pathway genes that were relatively common and validated at study inception. Using TaqMan, we genotyped a T432G aspartate to glutamate SNP (rs7041) and an A436C threonine to lysine SNP (rs4588), both in exon 11 of the GC gene, which codes for the vitamin D-binding protein (DBP). We also measured a synonymous A→G SNP in exon 4 of the CYP24A1 gene (rs2296241), which encodes 24(OH)ase, which initiates degradation of 1,25(OH)2D.
Concordance for quality control samples was 100%. The overall success rate for genotyping assays was at least 95%, except for Bsm1 (93%). All SNPs were in Hardy-Weinberg equilibrium (p >0.05).
Calcium, vitamin D, and dairy product intakes were measured in 1992 to 1993  using a semi-quantitative 68-item food frequency questionnaire (FFQ), which is a modification of the brief 'Health Habits and History Questionnaire' developed by Block and colleagues . Dietary and total nutrient intakes were estimated using the Diet Analysis System version 3.8a (Division of Cancer Prevention and Control, National Cancer Institute, Bethesda, MD, USA) . Total calcium estimates (mg/day) included contributions from diet, individual calcium supplements (250, 500, 600, or 750 mg per tablet), and multivitamin pills (estimated at 130 mg per multivitamin pill). Vitamin D values were added to the nutrient database using US Department of Agriculture sources . Total intake of vitamin D (IU/day) was estimated from diet (fortified milk, cereals, and fish) and multivitamin usage (400 IU per pill). Nutrient estimates were adjusted for total energy using the residual method .
The FFQ was validated using four random 24-hour recalls collected throughout a 1-year period as the comparison measure among 441 Nutrition Cohort participants . Median energy-adjusted, attenuation-corrected Pearson validity correlations and 95% confidence intervals (CIs) among women were high for dietary calcium (r = 0.66, 95% CI 0.51 to 0.76) and for dairy products (r = 0.63, 95% CI 0.45 to 0.77).
We used a χ2 test to assess whether VDR, GC, and CYP24A1 genotype distributions were in Hardy-Weinberg equilibrium. Haplotypes of the VDR SNPs Bsm1,Apa1, and Taq1 and the poly(A) repeat were estimated using the expectation-maximization algorithm implemented in the TAGSNPS program . To account for matching, conditional logistic regression was used to examine the association between all SNPs and the VDR haplotypes and postmenopausal breast cancer risk. We adjusted for additional risk factors and confounders as noted below. The main analysis combined invasive (n = 393) and in situ (n = 107) breast tumors; sensitivity analyses excluding in situ cases showed similar findings; thus, all cases are presented.
Covariates included education, age at menopause and menarche, parity, age at first live birth, personal history of benign breast disease, family history of breast cancer in mother or sisters, postmenopausal hormone replacement therapy use, and body mass index. Adjustment for alcohol consumption did not change the effect estimates and was not included.
We examined whether dietary or total calcium or vitamin D modified associations by genotype; interaction terms between genotypes and each dietary variable in continuous form were included in the conditional logistic regression models. Total energy intake was included in the models examining dietary interactions. Two-way interactions of the VDR haplotypes, Fok1, GC SNPs and the CYP24A1 SNP were examined. Statistical interaction was evaluated using the Wald test; a p value of less than 0.05 was considered statistically significant.
Descriptive characteristics of participants in this study are presented in Table 1. Cases and controls had a median age of 62 years (range 43 to 75 years) at enrollment in 1992, and 99% were Caucasian. None of the VDR genotypes or polymorphisms in the vitamin D-binding (GC) or CYP24A1 gene was significantly associated with overall postmenopausal breast cancer risk (Table 2). Heterozygote and homozygote variants for SNPs near the 3' end of the VDR gene were combined because of similar odds ratios (ORs). Three VDR haplotypes were identified that explained more than 96% of the variance among these SNPs: baTL (45.2%), BAtS (38.7%), and bATL (12.3%). Haplotype associations were similarly null: baTL OR = 0.89 (95% CI 0.66 to 1.20), BAtS OR = 1.24 (95% CI 0.92 to 1.66), and bATL OR = 0.96 (95% CI 0.69 to 1.32) for one or two copies versus none.
In analyses stratified by dietary factors, total calcium intake (but not vitamin D intake) modified the associations between polymorphisms near the 3' end of the VDR and breast cancer risk (Table 3). Compared to women with less than median total calcium intakes (<902 mg/day) and the BB or Bb Bsm1 genotype (reference group), women with calcium intakes of more than 902 mg/day and the bb genotype had a 39% lower risk of breast cancer (p interaction = 0.01). Similar findings were observed for those with the Taq1 TT and the poly(A) LL tail. These interactions were similar and statistically significant for dietary calcium, and no effect modification was observed with dietary vitamin D or dairy product consumption. Diet did not modify the associations between GC or CYP24A1 gene polymorphisms and breast cancer risk (not shown). We also observed no statistically significant effect modification in analyses stratified by median UV index for the state of residence as a proxy for cutaneous production of vitamin D. In addition, we observed no significant interactions between any two-way combination of VDR haplotypes, Fok1, or CYP24A1 or GC genotypes or in specific combinations of Fok1 and poly(A) genotypes.
None of the genotypes or haplotypes examined from the vitamin D pathway was associated with postmenopausal breast cancer risk in all women. However, women with a calcium intake higher than the median were at lower breast cancer risk if they had the VDR Bsm1 bb, Taq1 TT, or the poly(A) long repeat (LL) genotype. Among women with a calcium intake lower than the median, these VDR genotypes were not associated with lower breast cancer risk.
The lack of association between the Bsm1 bb genotype or associated polymorphisms in the 3' region and breast cancer risk is in contrast to findings in five studies that reported increased risk of breast cancer, mostly among Caucasian women with the bb genotype [7–11]. Four other studies, plus ours, have not replicated this association [12–15]; two of these were in Hispanic  and Taiwanese  women. Other SNPs in the baT haplotype allele likewise have been inconsistently related to breast cancer risk . It is unclear whether chance or underlying differences in the populations studied explain these inconsistencies. It is conceivable that differences in background risk factor profiles (menopausal status, reproductive factors, obesity, or dietary intake) contribute to these variations. Most of these studies included both pre- and postmenopausal cases; in the one study that presented results stratified by menopausal status , results for the BsmI BB versus the bb allele were similar for premenopausal (OR = 0.92) and postmenopausal (OR = 0.94) women.
In our study, dietary calcium intake, but not vitamin D intake, significantly modified the association between VDR SNPs near the 3' end of the gene and breast cancer risk. Both nutrients have direct effects on cell proliferation and differentiation of several cancer cell lines in vitro [28, 29]. In a rodent model, calcium reduces fat-induced cell proliferation by maintaining intracellular calcium concentrations . It is plausible that polymorphisms in the VDR gene could modify the effects of exposures such as vitamin D or calcium on breast cancer risk, although the specific mechanisms in humans are unknown. Although vitamin D intake did not modify risk in our analysis, we did not measure plasma 25(OH)D, an integrated measure of vitamin D status from diet and UVB synthesis. Higher circulating levels of the storage form of vitamin D (25(OH)D) may lower breast cancer risk [9, 31], potentially by providing substrate for breast tissue-specific synthesis of 1,25(OH)2D . In one case-control study in the UK, women with low plasma 25(OH)D levels and the BsmI bb VDR genotype had a 6.8-fold higher risk of breast cancer than individuals with higher 25(OH)D levels and the BB or Bb genotype . Thus, calcium or related nutrients may influence tissue-specific vitamin D metabolism in ways that are as yet unclear, possibly varying by VDR genotype. The interactions we observed may also be due to chance and should be replicated.
Our null findings for the Fok1 ff SNP and breast cancer risk do not replicate the positive association observed in the Nurses' Health Study  but do agree with those of three smaller studies [7, 13, 17], nor did we replicate previously reported potential interactions between different VDR polymorphisms observed by others [8, 15]. Guy and colleagues  reported an approximately two-fold increase in risk among women with the FFLL or FfLL genotypes versus other combinations; among patients with metastases, the risk was three-fold greater. The mix of grade, lymph node involvement, or metastases in cases studied may also contribute to variability in findings across studies.
We observed no association with two non-synonymous polymorphisms in the GC gene (codes for DBP) or one in CYP24A1 and postmenopausal breast cancer. The polymorphic DBP facilitates vitamin D actions, and DBP alleles differ in their affinity for 1,25(OH)2 D . The CYP24A1 gene codes for 24(OH)ase, which degrades 1,25(OH)2D. The potential for the SNPs we examined to alter DBP or CYP24A1 function is possible but has not been studied. Some studies have suggested a relationship of these SNPs with other outcomes, including asthma, glucose tolerance, and bone mineral density [33–35]. We are not aware of any studies that examined these SNPs in relation to breast cancer.
Limitations of this analysis include the use of a single assessment of diet, on average 6 years before diagnosis of breast cancer, which may have led to misclassification of dietary intake. We did not measure blood levels of 25(OH)D to assess vitamin D status. However, our FFQ was validated  and diet was assessed before cancer diagnosis. Larger case numbers would have provided greater statistical power to examine gene-gene and gene-diet interactions; however, this study is among the largest published to date for postmenopausal breast cancer.
Our research does not support an independent association between selected polymorphisms in the VDR, GC, or CYP24A1 genes and postmenopausal breast cancer risk among US Caucasian women. However, calcium intake appeared to modify the association of VDR SNPs with breast cancer risk. More research is needed to determine whether environmental factors, particularly diet, influence cancer risk by vitamin D pathway genotype.
Cancer Prevention Study
vitamin D-binding protein
food frequency questionnaire
vitamin D receptor.
Giovannucci E: The epidemiology of vitamin D and cancer incidence and mortality: a review (United States). Cancer Causes Control. 2005, 16: 83-95. 10.1007/s10552-004-1661-4.
Colston K, Welsh J: Vitamin D and breast cancer. Vitamin D. Edited by: Feldman D, Glorieux FH, Pike JW. 2005, Stanford, CA: Elsevier Academic Press, II: 1663-1677.
Townsend K, Banwell CM, Guy M, Colston KW, Mansi JL, Stewart PM, Campbell MJ, Hewison M: Autocrine metabolism of vitamin D in normal and malignant breast tissue. Clin Cancer Res. 2005, 11: 3579-3586. 10.1158/1078-0432.CCR-04-2359.
Colston KW, Berger U, Coombes RC: Possible role for vitamin D in controlling breast cancer cell proliferation. Lancet. 1989, 1: 188-191. 10.1016/S0140-6736(89)91204-X.
Ingles SA, Haile RW, Henderson BE, Kolonel LN, Nakaichi G, Shi C, Yu MC, Ross RK, Coetzee GA: Strength of linkage disequilibrium between two vitamin D receptor markers in five ethnic groups: implications for association studies. Cancer Epidemiol Biomarkers Prev. 1997, 6: 93-98.
Slattery ML, Yakumo K, Hoffman M, Neuhausen S: Variants of the VDR gene and risk of colon cancer (United States). Cancer Causes Control. 2001, 12: 359-364. 10.1023/A:1011280518278.
Bretherton-Watt D, Given-Wilson R, Mansi JL, Thomas V, Carter N, Colston KW: Vitamin D receptor gene polymorphisms are associated with breast cancer risk in a UK Caucasian population. Br J Cancer. 2001, 85: 171-175. 10.1054/bjoc.2001.1864.
Guy M, Lowe LC, Bretherton-Watt D, Mansi JL, Peckitt C, Bliss J, Wilson RG, Thomas V, Colston KW: Vitamin D receptor gene polymorphisms and breast cancer risk. Clin Cancer Res. 2004, 10: 5472-5481. 10.1158/1078-0432.CCR-04-0206.
Lowe LC, Guy M, Mansi JL, Peckitt C, Bliss J, Wilson RG, Colston KW: Plasma 25-hydroxy vitamin D concentrations, vitamin D receptor genotype and breast cancer risk in a UK Caucasian population. Eur J Cancer. 2005, 41: 1164-1169. 10.1016/j.ejca.2005.01.017.
Yamagata Z, Zhang Y, Asaka A: Association of breast cancer with vitamin D receptor gene polymorphism (abstract). Am J Hum Genet. 1997, 61: A388-
Ruggiero M, Pacini S, Aterini S, Fallai C, Ruggiero C, Pacini P: Vitamin D receptor gene polymorphism is associated with metastatic breast cancer. Oncol Res. 1998, 10: 43-46.
Buyru N, Tezol A, Yosunkaya-Fenerci E, Dalay N: Vitamin D receptor gene polymorphisms in breast cancer. Exp Mol Med. 2003, 35: 550-555.
Ingles SA, Garcia DG, Wang W, Nieters A, Henderson BE, Kolonel LN, Haile RW, Coetzee GA: Vitamin D receptor genotype and breast cancer in Latinas (United States). Cancer Causes Control. 2000, 11: 25-30. 10.1023/A:1008979417618.
Hou MF, Tien YC, Lin GT, Chen CJ, Liu CS, Lin SY, Huang TJ: Association of vitamin D receptor gene polymorphism with sporadic breast cancer in Taiwanese patients. Breast Cancer Res Treat. 2002, 74: 1-7. 10.1023/A:1016048900049.
Chen WY, Bertone-Johnson ER, Hunter DJ, Willett WC, Hankinson SE: Associations between polymorphisms in the vitamin D receptor and breast cancer risk. Cancer Epidemiol Biomarkers Prev. 2005, 14: 2335-2339. 10.1158/1055-9965.EPI-05-0283.
Uitterlinden AG, Fang Y, van Meurs JBJ, Pols HAP: Genetic vitamin D receptor polymorphisms and risk of disease. Vitamin D. Edited by: Feldman D, Glorieux FH, Pike JW. 2005, Stanford, CA: Elsevier Academic Press, II:
Curran JE, Vaughan T, Lea RA, Weinstein SR, Morrison NA, Griffiths LR: Association of a vitamin D receptor polymorphism with sporadic breast cancer development. Int J Cancer. 1999, 83: 723-726. 10.1002/(SICI)1097-0215(19991210)83:6<723::AID-IJC4>3.0.CO;2-3.
Holick M: Vitamin D. Modern Nutrition in Health and Disease. Edited by: Shils ME, Olson JA, Shike M, Ross AC. 1999, Philadelphia, PA: Lippincott Williams & Wilkins, 9
Cross HS, Kallay E, Lechner D, Gerdenitsch W, Adlercreutz H, Armbrecht HJ: Phystoestrogens and vitamin D metabolism: a new concept for the prevention and therapy of colorectal, prostate, and mammary carcinomas. J Nutr. 2004, 134: 1207S-1212S.
Calle EE, Rodriguez C, Jacobs EJ, Almon ML, Chao A, McCullough ML, Feigelson HS, Thun MJ: The American Cancer Society Cancer Prevention Study II Nutrition Cohort: rationale, study design, and baseline characteristics. Cancer. 2002, 94: 2490-2501. 10.1002/cncr.101970.
Rothman KJ, Greenland S: Modern Epidemiology. 1998, Philadelphia, PA: Lippincott Williams & Wilkins, 2
Block G, Hartman AM, Naughton D: A reduced dietary questionnaire: development and validation. Epidemiology. 1990, 1: 58-64. 10.1097/00001648-199001000-00013.
Block G, Coyl L, Smucker R, Harlan LC: Health Habits and History Questionnaire: diet history and other risk factors. [Personal computer system documentation]. 1989, Bethesda, MD: National Cancer Institute Division of Cancer Prevention and Control, National Institutes of Health
US Department of Agriculture: Provisional Table On Vitamin D Content Of Foods. [http://www.nal.usda.gov/fnic/foodcomp/Data/Other/vit_d99.pdf]
Willett W, Stampfer MJ: Total energy intake: implications for epidemiologic analyses. Am J Epidemiol. 1986, 124: 17-27.
Flagg EW, Coates RJ, Calle EE, Potischman N, Thun MJ: Validation of the American Cancer Society Cancer Prevention Study II Nutrition Survey Cohort Food Frequency Questionnaire. Epidemiology. 2000, 11: 462-468. 10.1097/00001648-200007000-00017.
Stram DO, Leigh Pearce C, Bretsky P, Freedman M, Hirschhorn JN, Altshuler D, Kolonel LN, Henderson BE, Thomas DC: Modeling and E-M estimation of haplotype-specific relative risks from genotype data for a case-control study of unrelated individuals. Hum Hered. 2003, 55: 179-190. 10.1159/000073202.
Lipkin M, Newmark HL: Vitamin D, calcium, and prevention of breast cancer: a review. J Am Coll Nutr. 1999, 18: 392S-397S.
Martinez ME, Willett WC: Calcium, vitamin D and colorectal cancer: a review of the epidemiologic evidence. Cancer Epidemiol Biomarkers Prev. 1998, 7: 163-168.
Jacobson EA, James KA, Newmark HL, Carroll KK: Effects of dietary fat, calcium, and vitamin D on growth and mammary tumorigenesis induced by 7,12-dimethylbenz(a)anthracene in female Sprague-Dawley rats. Cancer Res. 1989, 49: 6300-6303.
Bertone-Johnson ER, Chen WY, Holick MF, Hollis BW, Colditz GA, Willett WC, Hankinson SE: Plasma 25-hydroxyvitamin D and 1,25-dihydroxyvitamin D and risk of breast cancer. Cancer Epidemiol Biomarkers Prev. 2005, 14: 1991-1997. 10.1158/1055-9965.EPI-04-0722.
Pani MA, Regulla K, Segni M, Hofmann S, Hufner M, Pasquino AM, Usadel KH, Badenhoop K: A polymorphism within the vitamin D-binding protein gene is associated with Graves' disease but not with Hashimoto's thyroiditis. J Clin Endocrinol Metab. 2002, 87: 2564-2567. 10.1210/jc.87.6.2564.
Ezura Y, Nakajima T, Kajita M, Ishida R, Inoue S, Yoshida H, Suzuki T, Shiraki M, Hosoi T, Orimo H, Emi M, et al: Association of molecular variants, haplotypes, and linkage disequilibrium within the human vitamin D-binding protein (DBP) gene with postmenopausal bone mineral density. J Bone Miner Res. 2003, 18: 1642-1649. 10.1359/jbmr.2003.18.9.1642.
Wjst M, Altmuller J, Faus-Kessler T, Braig C, Bahnweg M, Andre E: Asthma families show transmission disequilibrium of gene variants in the vitamin D metabolism and signalling pathway. Respir Res. 2006, 7: 60-10.1186/1465-9921-7-60.
Hirai M, Suzuki S, Hinokio Y, Hirai A, Chiba M, Akai H, Suzuki C, Toyota T: Variations in vitamin D-binding protein (group-specific component protein) are associated with fasting plasma insulin levels in Japanese with normal glucose tolerance. J Clin Endocrinol Metab. 2000, 85: 1951-1953. 10.1210/jc.85.5.1951.
The authors declare that they have no competing interests.
MLM conceived the study and drafted the manuscript. VLS, HSF, EEC, and RMB made substantial contributions to the conception and design of the study. MJT, CR and EEC made substantial contributions to the acquisition and interpretation of data. WRD conducted the statistical analyses. All authors were involved in revising the manuscript critically for important intellectual content. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.