The Int7G24A variant of transforming growth factor-beta receptor type I is a risk factor for colorectal cancer in the male Spanish population: a case-control study
- Adela Castillejo†1,
- Trinidad Mata-Balaguer†1,
- Carla Guarinos1,
- María-Isabel Castillejo1,
- Ana Martínez-Cantó1,
- Víctor-Manuel Barberá1, 2,
- Paola Montenegro1,
- Enrique Ochoa3,
- Rafael Lázaro4,
- Carmen Guillén-Ponce2,
- Alfredo Carrato1 and
- José-Luís Soto1, 2Email author
© Castillejo et al; licensee BioMed Central Ltd. 2009
Received: 10 June 2009
Accepted: 20 November 2009
Published: 20 November 2009
The Int7G24A variant of transforming growth factor-beta receptor type I (TGFBR1) has been shown to increase the risk for kidney, ovarian, bladder, lung and breast cancers. Its role in colorectal cancer (CRC) has not been established. The aims of this study were to assess the association of TGFBR1*Int7G24A variant with CRC occurrence, patient age, gender, tumour location and stage.
We performed a case-control study with 504 cases of sporadic CRC; and 504 non-cancerous age, gender and ethnically matched controls. Genotyping analysis was performed using allelic discrimination assay by real time PCR.
The Int7G24A variant was associated with increased CRC incidence in an additive model of inheritance (P for trend = 0.005). No significant differences were found between Int7G24A genotypes and tumour location or stage. Interestingly, the association of the Int7G24A variant with CRC risk was significant in men (odds ratio 4.10 with 95% confidence intervals 1.41-11.85 for homozygous individuals; P for trend = 0.00023), but not in women. We also observed an increase in susceptibility to CRC for individuals aged less than 70 years.
Our data suggest that the Int7G24A variant represents a risk factor for CRC in the male Spanish population.
Transforming growth factor beta (TGF-β) is one of the most potent inhibitors of proliferation in epithelial, neuronal and hematopoietic cells. The effects of an activated TGF-β signalling pathway include inhibition of cell cycle progression, promotion of terminal differentiation and activation of cell death .
TGF-β receptor type I (encoded by TGFBR1) is a mediator of TGF-β growth inhibitory signals. Extracellular TGF-β binds first to the type II receptor, which then dimerizes and activates the type I receptor, sending the signal to the nucleus via SMAD proteins . Abnormalities of this signalling pathway are almost universal in cancer cells through a variety of mechanisms . Among them are overexpression of the ligands, downregulation of receptors, and point mutations and deletions in the genes coding for proteins involved in the pathway . Different TGFBR1 allelic variants have been associated with susceptibility to develop different types of tumours. Thus, the TGFBR1*6A allele, located at exon 1 (rs11466445) has been shown to be a susceptibility allele for colorectal cancer (CRC), breast and ovarian cancers, based on a meta-analysis of many case-control studies . In contrast, one recent report found a lack of an association between this variant and CRC risk [5, 6]. The role of this variant in CRC susceptibility remains controversial.
The TGFBR1_Int7G24A polymorphism (rs334354; NM_004612.2:c.1255+24G>A) has also been associated with cancer risk, but the function of this variant has yet to be discovered. Individuals carrying the Int7G24A variant have an increased risk of developing non-small cell lung cancer, renal cell carcinomas, transitional cell carcinomas of the bladder, and breast cancer [7–10]. To our knowledge, there is only one published report on the association between this genetic variant and CRC . In that study, there was no evidence of an association between this polymorphism and colorectal cancer risk in familial cases.
To test whether the Int7G24A variant might be associated with the incidence of CRC in the Spanish population, we carried out a case-control study with 504 cases and 504 matched controls for age, gender and ethnicity. We found a significant association of the Int7G24A variant with the incidence of CRC. It also showed significant associations with patient gender (male) and age, but not with tumour location or stage. These data support the hypothesis that Int7G24A variant may be a low penetrance tumour susceptibility allele that predisposes to CRC mainly in men in the Spanish population.
A total of 504 sporadic cases of CRC and 504 controls from the Elche University Hospital and Castellon Provincial Hospital tissue banks were analysed in this study. Written consent was obtained from every patient to be included in the respective tissue banks. Patients diagnosed with a familial cancer syndrome were excluded. Thus, patients with Amsterdam II criteria for hereditary non-polyposis colorectal cancer were excluded, as were patients with Bethesda criteria and who were negative for microsatellite instability, who were studied in our genetic counselling unit.
This was a hospital-based case-control study. Controls were selected from the same hospitals with no personal history of cancer and selected according to diagnoses unrelated to the disease of interest. They were matched for age, gender and race/ethnicity with the study patients. The study was approved by the Ethics Committees from the Elche University and Castellon Provincial hospitals.
The median age for patients with CRC at diagnosis was 70 years (range 23-92); for controls it was 72 years (range 23-98). Gender distribution for the patients with CRC ('cases') was 59.4% men and 40.6% women; for controls, it was 53.3% and 46.7%, respectively. In 33.3% of the cases, tumours were in the proximal colon and 66.7% were located at the distal colon and rectum; 45% were classified as low stage (I and II) and 55% as high stage (III and IV). The numbers of cases and controls for whom information was unavailable for stratification were 64 and one, respectively.
DNA and RNA isolation and cDNA synthesis
Control DNA samples were obtained from peripheral blood samples. DNA samples from cases were isolated after mechanical homogenization of non-cancerous frozen colorectal tissue. DNA isolation was performed using an EZ1 BioRobot (Qiagen, Valencia, CA, USA) according to the manufacturer's instructions. RNA was isolated from normal colon tissue of selected cases in the same manner. Synthesis of cDNA was performed using random primers and Reverse Transcription reagents (Applied Biosystems, Foster City, CA, USA).
Genotyping the Int7G24Apolymorphism (rs334354)
A total of 1008 individuals (504 cases and 504 controls) were analysed for this polymorphism. Single nucleotide polymorphism (SNP) analysis was performed using a real-time polymerase chain reaction (PCR) TaqMan assay for allelic discrimination (TaqMan SNP Genotyping Assay; ID: C_1413390_20; Applied Biosystems) with minor modifications. All PCR reactions contained 20 ng DNA, 10 μL TaqMan Universal Master Mix (Applied Biosystems), 0.5 μL of the assay and water for a final volume of 20 μL. The appropriate negative controls were also run in parallel. Real time PCR was performed on an ABI Prism 7300 Sequence Detection System (Applied Biosystems) using the following conditions: 95°C for 10 min and 40 cycles of amplification (95°C for 15 s and 60°C for 1 min). Software in the system was used to quantify the fluorescent signals from VIC- or FAM-labelled probes. Two representative cases for each genotype were sequenced directly to confirm the genotyping results. Sequencing reactions were performed using BigDye Terminator (version 3.1) Cycle Sequencing kits and run in an ABI Prism 3100-Avant Genetic Analyser (Applied Biosystems).
Analysis of the lengths of Int7G24 and Int7A24 variant transcripts
Four cases previously genotyped as G/G and four genotyped as A/A were analysed to compare the length of their TGFBR1 transcripts. This experiment was performed using PCR and capillary electrophoresis. Primer sequences designated IVS7cDNA-F (labelled with 6-FAM) were 5'-AAATTGCTCGACGATGTTCC-3' for exon 7, and IVS7cDNA-R, 5'-CTCTGCCATCTGTTTGGGAT-3' for exon 8. The expected size of a normal PCR product was 151 bp. Twenty-five nanograms of DNA were used in a total PCR reaction volume of 25 μL AmpliTaq Gold PCR Master Mix 2 × (Applied Biosystems) and 20 pmol of each primer. After an initial 9 min at 95°C, PCR was run for 40 cycles of 95°C for 20 s; 55°C for 20 s and 72°C for 20 s; followed by 10 min final extension at 72°C. Diluted (1/10) PCR products along with Genescan 500 ROX size standards in HiDi formamide (Applied Biosystems,) were resolved by capillary electrophoresis (ABI 3100 Avant, Applied Biosystems) using POP6 polymer. GeneScan Software was used for the analysis of PCR fragment lengths (Applied Biosystems).
We first examined the Hardy-Weinberg equilibrium of allelic distribution separately for cases and controls and then compared the allele frequencies between cases and controls. We evaluated the association of Int7G24A genotypes with CRC using multivariate unconditional logistic regression models assuming dominant (A/A and G/A, versus G/G), additive (G/G and G/A, versus A/A) or recessive (A/A, versus G/G and G/A) modes of inheritance, respectively. In our basic models, we adjusted for age and gender. We explored the potential effects of modification by age in an analysis stratified by the cases' median age at diagnosis (≤70 or >70 years). A χ2 test was used to evaluate the differences in Int7G24A variant carrier frequencies between tumour and control groups for the dominant and the recessive models and to analyse any association between the SNP and clinical and pathological factors. Armitage's trend test was used to calculate P for trend in the additive model of inheritance. All P values were two-sided and P < 0.05 was considered significant. Results are expressed as the odds ratio (OR) and 95% confidence interval (CI).
Genotypes of the Int7G24A polymorphism in cases and controls and their association with CRC.
(n = 504)
(n = 504)
OR [95% CI]
G/A and A/A
G/G and G/A
≤ 70 yrs
(n = 213)
≤ 70 yrs
(n = 219)
G/G, number (%)
G/A, number (%)
A/A, number (%)
(n = 227)
(n = 284)
G/G, number (%)
G/A, number (%)
A/A, number (%)
(n = 295)
(n = 226)
G/G, number (%)
G/A, number (%)
A/A, number (%)
(n = 202)
(n = 277)
G/G, number (%)
G/A, number (%)
A/A, number (%)
The results after gender adjustment were clear-cut, with no CRC risk association found for women and a strong statistically significant association for the presence of the Int7G24A variant and CRC in men (P for trend = 0.00023). The ORs for heterozygous and homozygous individuals for the Int7G24A variant were 1.74 (95% CI: 1.19-2.54) and 4.10 (95% CI: 1.41-11.85), respectively (Table 1). No significant associations were found between different genotypes and tumour location or stage.
Lengths of Int7G24 and Int7A24transcripts
We performed mRNA analysis from homozygous G/G and A/A individuals to test if the Int7G24A variant might be related to an abnormal mRNA splicing that could explain the function of this allele as a risk factor for CRC. However, the mRNA fragments obtained always had the expected size for normal mRNA splicing independently of the genotype.
The TGFβ pathway plays an important role in tumour development and progression . Recently, there has been interest in the contribution of genetic variation in this pathway to cancer susceptibility. Several polymorphic variants in genes involved in this pathway have been studied . The TGFBR1_Int7G24A variant has been previously associated with cancer risk, but the role of this variant has yet to be discovered. Zhang et al  reported that homozygous but not heterozygous carriers of this variant have an increased risk of developing non-small cell lung cancer. Chen et al  reported the association of this variant with the incidence of renal cell carcinomas (OR: 2.2; 95% CI: 1.22-3.96) and transitional cell carcinoma of the bladder (OR: 2.45; 95% CI: 1.8-3.16). In a meta-analysis of these three studies , the pooled OR was 1.76 (95% CI: 1.33-2.34; P < 0.0001). Another case-control study  found an association with breast cancer (OR: 2.42; 95% CI: 1.55-3.79), suggesting that the Int7G24A variant could also represent a risk factor for invasive and metastatic breast cancers (OR: 2.61; 95% CI: 1.65-4.11). Skoglund-Lundin et al.  did not detect a significant association between CRC risk and the Int7G24A variant when they analysed data from cases with a familial history of colorectal cancer. The expectedly low-penetrance effect of this putative susceptibility allele would be very difficult to detect in cases with high-penetrance alleles. Moreover, the authors did not stratify their data according to sex.
An exclusion criterion for the present study was the diagnosis of any familial cancer syndrome predisposing to CRC, such as hereditary non-polyposis colorectal cancer, familial adenomatous polyposis, Cowden disease and Peutz-Jeghers syndrome. We were interested to evaluate the effect of this SNP for sporadic CRC susceptibility and tumour progression. The potential effect of the SNP on familial syndromes as modifier of the penetrance, or the phenotype expression should be approached independently in further specific studies.
Control individuals were selected from those patients treated at an emergency department with diverse non-cancerous pathologies and with no personal history of cancer. The probability that controls suffered from diseases putatively caused by the same risk factors as the patients with CRC is extremely low.
It is well known that one of the main limitations and sources of confounding the results in association studies comes from population stratification and inappropriate sample sizes. A priori, our study with 504 sporadic cases and 504 controls allowed us to detect an OR of 1.5 for G/A heterozygous individuals (considering a frequency of 0.3) and an OR of 2.0 for A/A homozygous individuals (considering a frequency of 0.025) with 80% power (two-sided test, alpha level 5%). Estimated genotype frequencies were obtained as preliminary results in a subset of our control population and were similar to the frequencies registered at the SNP database for European populations http://www.ncbi.nlm.nih.gov/projects/SNP/. The results presented here show an increase of CRC susceptibility for those individuals carrying the Int7G24A variant allele with a dosage effect (additive model of inheritance). Adjusting the OR by age, we only observed an association of Int7G24A variant carriers with CRC risk in subjects aged less than 70 years (Table 1). Thus, the genetic risk effect is more evident in younger individuals. Ageing implies a huge amount of acquired predisposing carcinogenic factors that could overshadow the mild effect of the susceptibility alleles.
When we analysed the data adjusting by gender, we found that the Int7G24A variant was strongly associated with CRC risk in men but not women. We do not know why. Evidence is accumulating supporting gender-related differences in the development of CRC, but to our knowledge, few reports have been published regarding sex-specific relationship between SNPs and CRC. Slattery et al  reported that the oestrogen receptor beta CA repeat polymorphism was associated with an increase relative risk of colon cancer in women but not in men. On the contrary, increasing number of CAG repeats in the androgen receptor was directly associated with colon cancer among men, but not in women. Both oestrogen and androgen receptors are present in colorectal tissue and it has been suggested that these might be important in regulating the CRC risk associated with these hormones. Activation of the TGFβ signal transduction system is also subject to hormonal regulation . This striking connection might suggest the involvement of sex steroid hormones acting through the TGFβ pathway in the aetiology of CRC.
The apolipoprotein E (ApoE) epsilon 2/3 polymorphism is another example of gender-specific modulation of CRC risk and prognosis where men have a highly significant association but there is no association in women . Bae et al  found a gender-specific association between polymorphisms of the gene for vascular endothelial growth factor (VEGF 936 C>T) and CRC in the Korean population. T allele-bearing genotypes significantly increased the risk for CRC in women but not in men. Moreover, there are gender-related survival differences in association with EGFR (gene for epidermal growth factor receptor) polymorphisms in patients with metastatic colon cancers . Again, oestrogens and androgens may be behind all these associations because of their functional links with ApoE , VEGF  and EGFR .
It is possible that the positive associations shown here were not based on real effects of this polymorphism, but rather reflect unknown differences in population ancestry between the case and control groups . We consider the probability of false-positive inference attributable to population stratification to be small, because the cases and control individuals were recruited from an ethnically homogeneous population with no indication of a significant amount of recent genetic admixture.
Given our results and those in the literature, it seems reasonable to speculate that there might be a gender-dependent allelic architecture for CRC risk associated with sex steroids. This would introduce a new level of complexity in the 'common disease-common variant' hypothesis . Our association studies between this polymorphism and clinicopathological factors did not show any statistically significant results. The lack of difference in genotype distribution between male and female patients with CRC might be attributed to the contribution and compensation of other gender-related susceptibility polymorphisms, where the risks affect women preferentially. To date, how this intronic SNP affects oncogenesis remains to be elucidated. It has been proposed that this SNP creates an alternative splice site within intron 7. In vitro analysis with breast and ovarian cell lines from carriers of the Int7G24A variant have shown the retention of the seventh intron in the mRNA up to the site of the Int7G24A variant . We wanted to study if this alternative splicing also occurred in the normal colorectal tissue from patients with CRC. However, we concluded that factors other than alternative splicing might account for the association of this polymorphism with cancer, although we cannot rule out the possibility that Int7G24A variant is a marker representing a TGFBR1 haplotype. In this regard, it is interesting to note the likely 'bystander susceptibility effect' for the TGFBR1*6A allele because of linkage disequilibrium with the unknown causative mutation of the allele-specific expression (ASE) in the TGFBR1 gene . In our series, no linkage disequilibrium was found between the TGFBR1*6A and Int7G24A variant (unpublished results). No information is available regarding the association between variant and ASE. Further studies are required to characterize the molecular mechanisms by which variant is involved in gender-specific susceptibility to CRC.
The results of the present study strongly suggest that the TGFBR1_Int7G24A allele does confer an increased risk of colorectal cancer in the male Spanish population. The growing list of gender-dependent association of allelic variants and the risk of colorectal cancer might suggest the involvement of sex steroid hormones in the aetiology of CRC.
This action has been performed within the cooperation framework established by the Transversal Cancer Action approved by the Council of Ministers on October 11th 2007 in accordance with the agreement between The Carlos III Health Institute (ISCIII), which is an autonomous entity currently belonging to the Ministry of Science and Innovation, and the Biomedical Research Foundation from the Hospital of Elche. Research supported in part by grants from the Generalitat Valenciana in Spain (AP106/06) and the Biomedical Research Foundation from the Hospital of Elche, Spain (FIBElx-02/2007). T.M-B is recipient of a fellowship from the Spanish Society of Medical Oncology (SEOM).
- Siegel PM, Massagué J: Cytostatic and apoptotic actions of TGF-beta in homeostasis and cancer. Nat Rev Cancer. 2003, 3 (11): 807-821. 10.1038/nrc1208.View ArticlePubMed
- Bian Y, Kaklamani V, Reich J, Pasche B: TGF-beta signaling alterations in cancer. Cancer Treat Res. 2003, 115: 73-94. full_text.View ArticlePubMed
- Massagué J, Blain SW, Lo RS: TGFbeta signaling in growth control, cancer, and heritable disorders. Cell. 2000, 103 (2): 295-309. 10.1016/S0092-8674(00)00121-5.View ArticlePubMed
- Pasche B, Kaklamani V, Hou N, Young T, Rademaker A, Peterlongo P, Ellis N, Offit K, Caldes T, Reiss M, Zheng T: TGFBR1*6A and cancer: a meta-analysis of 12 case-control studies. J Clin Oncol. 2004, 22 (4): 756-758. 10.1200/JCO.2004.99.271.View ArticlePubMed
- Skoglund J, Song B, Dalén J, Dedorson S, Edler D, Hjern F, Holm J, Lenander C, Lindforss U, Lundqvist N, Olivecrona H, Olsson L, Påhlman L, Rutegård J, Smedh K, Törnqvist A, Houlston RS, Lindblom A: Lack of an association between the TGFBR1*6A variant and colorectal cancer risk. Clin Cancer Res. 2007, 13 (12): 3748-3752. 10.1158/1078-0432.CCR-06-2865.View ArticlePubMed
- Castillejo A, Mata-Balaguer T, Montenegro P, Ochoa E, Lázaro R, Martínez-Cantó A, Castillejo MI, Guarinos C, Barberá VM, Guillén-Ponce C, Carrato A, Soto JL: The TGFBR1*6Aallele is not associated with susceptibility to colorectal cancer in a Spanish population: a case-control study. BMC Cancer.
- Zhang HT, Fei QY, Chen F, Qi QY, Zou W, Wang JC, Zhang RM, Tao SH, Chen XF, Luo ZW: Mutational analysis of the transforming growth factor beta receptor type I gene in primary non-small cell lung cancer. Lung Cancer. 2003, 40 (3): 281-287.View ArticlePubMed
- Zhang HT: Int7G24A variant of the TGFBR1 gene and cancer risk: a meta-analysis of three case-control studies. Lung Cancer. 2005, 49 (3): 419-420. 10.1016/j.lungcan.2005.03.038.View ArticlePubMed
- Chen T, Jackson C, Costello B, Singer N, Colligan B, Douglass L, Pemberton J, Deddens J, Graff JR, Carter JH: An intronic variant of the TGFBR1 gene is associated with carcinomas of the kidney and bladder. Int J Cancer. 2004, 112 (3): 420-425. 10.1002/ijc.20419.View ArticlePubMed
- Chen T, Jackson CR, Link A, Markey MP, Colligan BM, Douglass LE, Pemberton JO, Deddens JA, Graff JR, Carter JH: Int7G24A variant of transforming growth factor-beta receptor type I is associated with invasive breast cancer. Clin Cancer Res. 2006, 12 (2): 392-397. 10.1158/1078-0432.CCR-05-1518.View ArticlePubMed
- Skoglund Lundin J, Vandrovcova J, Song B, Zhou X, Zelada-Hedman M, Werelius B, Houlston RS, Lindblom A: TGFBR1 variants TGFBR1(*)6A and Int7G24A are not associated with an increased familial colorectal cancer risk. Br J Cancer. 2009, 100 (10): 1674-1679. 10.1038/sj.bjc.6605054.PubMed CentralView ArticlePubMed
- Loktionov A: Common gene polymorphisms, cancer progression and prognosis. Cancer Lett. 2004, 208 (1): 1-33. 10.1016/j.canlet.2004.02.009.View ArticlePubMed
- Slattery ML, Sweeney C, Murtaugh M, Ma KN, Wolff RK, Potter JD, Caan BJ, Samowitz W: Associations between ERalpha, ERbeta, and AR genotypes and colon and rectal cancer. Cancer Epidemiol Biomarkers Prev. 2005, 14 (12): 2936-2942. 10.1158/1055-9965.EPI-05-0514.View ArticlePubMed
- Buck MB, Knabbe C: TGF-beta signaling in breast cancer. Ann N Y Acad Sci. 2006, 1089: 119-126. 10.1196/annals.1386.024.View ArticlePubMed
- Watson MA, Gay L, Stebbings WS, Speakman CT, Bingham SA, Loktionov A: Apolipoprotein E gene polymorphism and colorectal cancer: gender-specific modulation of risk and prognosis. Clin Sci (Lond). 2003, 104 (5): 537-545. 10.1042/CS20020329.View Article
- Bae SJ, Kim JW, Kang H, Hwang SG, Oh D, Kim NK: Gender-specific association between polymorphism of vascular endothelial growth factor (VEGF 936 C>T) gene and colon cancer in Korea. Anticancer Res. 2008, 28 (2B): 1271-1276.PubMed
- Press OA, Zhang W, Gordon MA, Yang D, Lurje G, Iqbal S, El-Khoueiry A, Lenz HJ: Gender-related survival differences associated with EGFR polymorphisms in metastatic colon cancer. Cancer Res. 2008, 68 (8): 3037-3042. 10.1158/0008-5472.CAN-07-2718.View ArticlePubMed
- Wang JM, Irwin RW, Brinton RD: Activation of estrogen receptor alpha increases and estrogen receptor beta decreases apolipoprotein E expression in hippocampus in vitro and in vivo. Proc Natl Acad Sci USA. 2006, 103 (45): 16983-16988. 10.1073/pnas.0608128103.PubMed CentralView ArticlePubMed
- Ruohola JK, Valve EM, Karkkainen MJ, Joukov V, Alitalo K, Härkönen PL: Vascular endothelial growth factors are differentially regulated by steroid hormones and antiestrogens in breast cancer cells. Mol Cell Endocrinol. 1999, 149 (1-2): 29-40. 10.1016/S0303-7207(99)00003-9.View ArticlePubMed
- Conzen SD: Minireview: nuclear receptors and breast cancer. Mol Endocrinol. 2008, 22 (10): 2215-2228. 10.1210/me.2007-0421.PubMed CentralView ArticlePubMed
- Marchini J, Cardon LR, Phillips MS, Donnelly P: The effects of human population structure on large genetic association studies. Nat Genet. 2004, 36 (5): 512-517. 10.1038/ng1337.View ArticlePubMed
- Pritchard JK, Cox NJ: The allelic architecture of human disease genes: common disease-common variant...or not?. Hum Mol Genet. 2002, 11 (20): 2417-2423. 10.1093/hmg/11.20.2417.View ArticlePubMed
- Markey MP, Graff J, Carter J, Chen T: Alternative splicing of a cancer-associated intronic variant of TGFBR1 in cancer cell lines. Proc Amer Assoc Cancer Research. 2006, 47: 2475-
- Valle L, Serena-Acedo T, Liyanarachchi S, Hampel H, Comeras I, Li Z, Zeng Q, Zhang HT, Pennison MJ, Sadim M, Pasche B, Tanner SM, de la Chapelle A: Germline allele-specific expression of TGFBR1 confers an increased risk of colorectal cancer. Science. 2008, 321 (5894): 1361-1365. 10.1126/science.1159397.PubMed CentralView ArticlePubMed
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/9/406/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.