Black African populations are more genetically diverse than others, but genetic variants have been studied primarily in European populations. The present study examined the association of four single nucleotide polymorphisms (SNPs) of the fibroblast growth factor receptor 2, associated with breast cancer in non-African populations, with breast cancer in Black, southern African women. Genomic DNA was extracted from whole blood samples of 1,001 patients with breast cancer and 1,006 controls (without breast cancer), and the rs2981582, rs35054928, rs2981578, and rs11200014 polymorphisms were analyzed using allele-specific Kompetitive allele-specific PCR™, and the χ2 or Fisher's exact tests were used to compare the genotype frequencies. There was no association between those SNPs and breast cancer in the studied cohort, although an association was identified between the C/C homozygote genotype for rs2981578 and invasive lobular carcinoma. These results show that genetic biomarkers of breast cancer risk in European populations are not necessarily associated with risk in sub-Saharan African populations. African populations are more heterogenous than other populations, and the information from this population can help focus genetic risks of cancer in this understudied population.
FGFR2SNPAfricanbreast cancerheterogenousriskNational Institutes of Health of the National Cancer Institute01-CA192627P30-CA136696University of the Witwatersrand/South African Medical Research Council Common Epithelial Cancer Research Center Grant; a South African National Research Foundation105646Cancer Association of South Africa; the National Research FoundationNRF 87935Capex from the University of the WitwatersrandAORTIC/NCI BIG CAT 2This study was funded by the National Institutes of Health of the National Cancer Institute (grant nos. 01-CA192627 and P30-CA136696); the University of the Witwatersrand/South African Medical Research Council Common Epithelial Cancer Research Center Grant; a South African National Research Foundation grant award (grant no. 105646); the Cancer Association of South Africa; the National Research Foundation (grant no. NRF 87935); the minor Capex from the University of the Witwatersrand; and an AORTIC/NCI BIG CAT 2 Research grant.Introduction
African populations are more genetically diverse than Asian or European populations (1,2), but few genomic studies have been conducted in African populations (3). Only 2.4% of genome-wide association studies (GWAS) have included Africans or African Americans (2), and very little information is available regarding cancer genomics in African populations (3,4). A meta-analysis by Rotimi et al (3) found that between January 1990 and December 2019, only 0.329% of cancer-related publications globally focused on African populations, and only 0.016% were related to cancer genetics or genomics in Africa. Breast cancer is the most common type of cancer in women worldwide, including African women. However, it is much less common in African than in European populations. The 2018 GLOBOCAN report showed an estimated age-standardized incidence rate (ASIR) of 37.9/100,000 (compared with 113.2/100,000 in Belgium) and a lifetime risk (LR) of 1 in 25 for women under 75 years (5). Similarly, the South African national cancer registry reported an ASIR of 20.4/100,000 and an LR of 1 in 47 among black South African women (6).
Since the 1990s, the breast cancer 1 (BRCA1) (7) and BRCA2 (8) genes have been associated with hereditary breast cancer. Other rare, but highly penetrant genes include phosphatase and tensin homolog (PTEN), tumor protein P53 (TP53), epithelial cadherin (CDH1), and serine/threonine kinase (9,10). Moderate penetrance genes include checkpoint kinase 2 (CHEK2), BRCA1 interacting helicase 1, ataxia-telangiectasia mutated (ATM), or partner and localizer of BRCA2 (PALB2) (8,10). Studies have been performed in South Africa to examine some of these genes, particularly examining the effects of BRCA1 and BRCA2 in various ethnic populations. In the self-identified black population (with a sample size of 165), Eygelaar et al (11) found 1.2% BRCA1, 0.6% BRCA2; 0.6% ATM, 0.6% CHEK2, and 0.6% PALB deleterious variants associated with breast cancer. Similarly, in 78 black patients, Francies et al (12) found 3.8% BRCA1 and 3.8% BRCA2 pathogenic mutations, but no deleterious mutations in PALB2 or CHEK2 in this group. Van der Merwe et al (13,14) identified larger rearrangements of the BRCA1 and BRCA2 genes that are specific to the black South African population. Deleterious mutations in high and medium penetrance genes do not explain the vast majority of breast cancers in the black population. Over the past 15 years, GWAS has led to the detection of over 200 loci associated with breast cancer (15–18). Among the top hits for these loci are variants in the fibroblast growth factor receptor 2 (FGFR2) gene, a low penetrance gene. The FGFRs are receptor tyrosine kinases involved in signaling pathways that catalyze a variety of biological processes, including cell growth, survival, differentiation, angiogenesis, tumorigenesis (19), and epithelial-to-mesenchymal transition (20). Variants in intron 2 of FGFR2 have been found to be highly associated with breast cancer (15,17,18,21–23). Although most single nucleotide polymorphisms (SNPs) have small effects individually, polygenic models indicate that an accumulation of small mutations may increase the risk of cancer (24). Included in the top hits of FGFR2 SNPs associated with breast cancer are rs2981582, rs35054928, and rs2981578 in women of European ancestry, and rs11200014 in African American women.
The FGFR2 SNP that is most commonly associated with breast cancer in women of European ancestry, rs2981582 (17), was recently also associated with an increased risk of breast cancer in Saudi Arabian women (25), and with luminal A breast cancer in Han Chinese women (26) and Korean women (27). The expression of the minor allele was associated with early-onset breast cancer in Indonesian women (28). However, a study of women from Argentina and Uruguay did not find an association between rs2981582 and breast cancer, possibly because the populations of those countries include subpopulations of varied ethnicity (admixed populations) (29). Likewise, among postmenopausal Turkish women, rs2981582 was not associated with breast cancer (30). Admixed populations, such as that of Turkey (31), may have different allele frequencies than European populations and, given similar sample sizes, less power to detect associations.
The risk allele for rs35054928 appears to bind the transcription factor, E2F1 (32). Both rs35054928 and rs2981578 are reported to be part of a response element, a sequence within the promoter of a gene that regulates transcription. The presence of the risk allele for rs2981578 for example, substantially increases the binding to FOXA1 in MCF7 epithelial, hormone receptor-positive cells, increasing chromatin accessibility and allowing access to transcriptional repressors such as Yin Yang 1 (YY1), SIN3A, and histone deacetylase (HDAC) (33). The DNA binding protein YY1, co-repressor, SIN3A, and histone-modifying HDAC form a complex that can inhibit promoter activity (34). The SNPs rs35054928 and rs2981578 are located next to an organic cation transporter (OCT)-binding site. The risk allele of rs2981578 also creates a potential binding site for runt-related transcription factor 2 (RUNX2) (20,35). Runx2 is primarily known for its role in osteoclast development, but it is also a regulator of mammary development and breast cancer (36). The rs2981578 variants have equal affinity for OCT1, but the high-risk allele has a much higher affinity for RUNX2 (33,35), possibly because the SNP sites differ in histone acetylation (37). OCT1 promotes cell proliferation in estrogen receptor (ER) positive breast cancer cells (38).
The SNP rs11200014 has been associated with breast cancer in African American women (18,39). African Americans generally have admixed African and European ancestry, and a small proportion of Native American ancestry (1,40), and their African ancestry is primarily from West or central West Africa (41). South African populations are genetically different from West and Central Africans, and differ even more from African Americans.
Thus, in the present study, the association between FGFR2 and breast cancer was explored, and their association with hormone receptor subtypes of breast cancer in an urban South African Black female population was assessed.
Materials and methodsStudy population
The median age [interquartile range (IQR)] of the cases was 53 years (44–64 years), and 51 years (40–62 years) in the control participants. Participants in the present study were drawn from the South African Breast Cancer and HIV Outcome (SABCHO) study, a cohort of breast cancer patients diagnosed and treated at five hospitals in Gauteng and KwaZulu Natal, South Africa (42). For this study, women diagnosed with breast cancer at the Charlotte Maxeke Johannesburg Academic Hospital (CMJAH) Surgical Breast Unit or the Batho Pele Breast Unit of the Chris Hani Baragwanath Academic Hospital (CHBAH) were selected. The CMJAH Surgical Breast Unit is located in central Johannesburg and identifies ~250 new breast cancer cases each year. The Batho Pele Breast Unit serves patients from Soweto and surrounding areas and diagnoses about 350 patients with breast cancer yearly (42). As controls, women seen in the CMJAH breast unit or the Batho Pele Breast unit who were found not to have breast cancer, and patients undergoing routine assessment at other clinics at CMJAH not related to cancer were used. Eligible cases were self-identified black, southern African women >18 years of age with histologically confirmed invasive breast cancer; exclusion criteria were patients with ductal carcinoma in situ or lobular carcinoma in situ. Eligible controls were self-identified black, southern African women >18 years of age with no history of breast or ovarian cancer; and neither pregnant nor breastfeeding. All participants had a sample of peripheral blood drawn (2–8 ml) and collected into EDTA vacutainer tubes (Becton, Dickson, and Company) between October 2014 to March 2020. Ethics clearance for this study was granted from the Human Research Ethics Committee (Medical) of the University of the Witwatersrand (Ethics numbers M140980, M161116). Written permission was granted by the CEOs of both CHBAH and CMJAH for the study.
Genomic DNA extraction, SNP selection, and analysis
DNA was extracted from whole blood using a modified salting out method. Briefly, whole blood was lysed with 320 mM sucrose in ice cold buffer (10 mM Tris-Cl, 5 mM MgCl2, 1% Triton X; 1 part blood: 4 parts lysis buffer), and was centrifuged at 900 × g for 10 min at 4°°C. The pellet was resuspended in fresh lysis buffer, recentrifuged and the resulting pellet was resuspended in 3 ml T20E5 (20 mM Tris-HCl, 5 mM EDTA), 200 µl 10% SDS and 495 µl proteinase K solution (2 mg/ml proteinase K; 1% SDS, 2mM EDTA) and incubated overnight at 42–50°C. Subsequently, 1 ml saturated NaCl (40%) was added to the solution, incubated on ice for 5 min and centrifuged at 900 × g for 30 min at 4°C. The supernatant was transferred to a clean tube where 20 ml absolute ethanol was added, which caused the DNA to precipitate out of solution. The DNA could then be spooled and transferred a clean 1.5-ml microcentrifuge tube. The DNA was airdried and dissolved in low TE buffer (10 mM Tris HCl; 0.1 mM EDTA, pH 8.3) and diluted to a final concentration of 25 ng/µl. This was based on the method of Miller et al (43). A Nanodrop 2000™ spectrophotometer (Thermo Fisher Scientific, Inc.) was used to determine DNA concentrations and the A260/280 ratios; a ratio between 1.7 and 2.0 indicated adequate DNA purity for this genotype analysis.
For this confirmatory candidate gene study, 4 SNPs were selected for genotyping; specifically, rs2981582, rs35054928, rs2981578, and rs11200014, which are located in intron 2 of the FGFR2 gene. Power for this study was assessed using the University of Michigan School of Public Health Genetic Association Study Power Calculator (https://csg.sph.umich.edu/abecasis/cats/gas_power_calculator/index.html). A sample size of 1,000 cases and 1,000 controls was chosen along with a significance level of 0.0125 (α=0.05/4). Assuming a dominant pattern of inheritance, odds ratios of 1.5 (for rs2981582) and 1.4 (for rs35054928, rs2981578, and rs11200014) could be detected with 80% power; and reflect odds ratios reported in the literature (23,44,45). Table I shows the minor allele frequencies (MAFs) for each SNP and compares them with the African and global allele frequencies in other studies. Information regarding allele frequencies from the African and global populations was obtained from the 1000 genomes project (46). Investigations of ancestry informative markers on a similar cohort sourced from this region showed limited evidence of population substructure (47).
The FGFR2 SNP polymorphisms were genotyped using Kompetitive allele-specific PCR (KASP™) technology at LGC Genomics Ltd. This trademarked method has 2 allele-specific forward primers and a common reverse primer. The forward primers each have a unique tail sequence that corresponds with a fluorescent resonant energy transfer cassette; one is labeled with FAM™ dye and the other with HEX™ dye. One allele binds the forward primer with FAM™ and the second allele binds the HEX™ labeled forward primer. During PCR, the allele-specific forward primer binds the DNA template, and subsequent PCR rounds generate the complement, which unquenches the fluorescent tag. If the genotype at a given SNP is homozygous, only one of two possible fluorescent signals will be generated, while a heterozygous genotype will generate a mixed fluorescent signal (48). The primers for each SNP are listed in Table II. For each SNP, the deviation of genotype frequencies from the Hardy-Weinberg equilibrium (HWE) in controls was determined using a χ2 test (Table I).
Classification of tumors
Histopathological characteristics, including histological diagnosis, tumor subtype and grade, and immunolocalization of estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and Ki67, were determined by histopathologists at the National Health Laboratory Services (NHLS) at CHBAH and CMJAH as part of patient standard of care. Pathological reports were produced for the clinical care of the patients, and selected data from these reports were included in the database for the present study. Data on these reports included tumor type, IHC, stage, in situ component. Immunostaining was performed on the benchmark XT automatic platform (Roche Diagnostics). Immunohistochemistry was performed according to the College of American Pathologists guidelines (49).
The ER/PR status was determined using Allred scoring as described previously (50) and scored as follows: 0–2, negative; and 3–8, positive. HER2 was regarded as positive if the score was 3+; negative when it was 0 or 1+ and equivocal when it was 2+; equivocal HER2 results were confirmed as positive by in situ hybridization as described previously (50). Specimens in which ≤14% of cells expressed Ki67 were categorized as having low expression, as per the St Gallen 2011 guidelines which form the basis of the current South African Guidelines (51,52). Immunostaining was performed according to the protocol in Bancroft's theory and practice of histological techniques (53) on an accredited diagnostic instrument as part of routine standard of care and according to the College of American Pathologists guidelines (49), as aforementioned. Based on IHC subtyping, the breast tumors were further categorized as: A-like or luminal A (ER and/or PR positive, HER2 negative, Ki67 ≤14%); B-like or luminal B (ER and/or PR positive and HER2 negative with Ki67% >14%); B/HER2-like or luminal B/HER2+ (ER and/or PR positive, HER2 positive with any Ki67); HER2 positive subtype (ER and PR negative, HER2 positive); and triple-negative breast cancer (TNBC; ER, PR, and HER2 negative) (51,54).
Statistical analysis
Continuous variables were assessed for normality using the Shapiro-Wilkes test. Normally distributed variables are presented as the mean ± SD, and non-normally distributed data as the median and IQR. Normally distributed continuous variables were compared using a Student's t-test, while non-normally distributed variables were compared using a Mann-Whitney U test. Categorical variables are presented as frequencies and percentages and were compared using a Pearson's χ2 test or Fisher's exact test if frequencies of >20% of cells were <5. All statistical analyses were performed using STATA version 14.2 (StataCorp LP). P<0.05 was considered to indicate a statistically significant difference.
Using the genotyping data provided by LGC genomics, both genotype and allele frequencies for each SNP of interest were calculated. The frequencies of each genotype (homozygous major, heterozygous, and homozygous minor) were compared between cases and controls for the SNPs of interest using a Pearson's χ2 test or Fisher's exact test as appropriate. Genotype frequencies were also examined under four different models of genetic disease risk (55). Results were analyzed using dominant, recessive, multiplicative, additive, and homozygous models. A Bonferroni correction for multiple comparison testing was applied. A value of P=0.008 (0.05/6) was considered to indicate a statistically significant difference.
Results
Overall, the 1,001 women with breast cancer were significantly older than the 1,006 women without breast cancer; the medians [IQRs] were 53 [44–64] and 51 [40–62], respectively (Table III). Most patients with cancer (57.4%) had been diagnosed with Stage III or IV cancer at presentation. The majority of tumors (82.5%) were invasive ductal carcinomas; 75.2% were ER-positive (76.7%), 63.8% were PR-positive, 26.9% were HER2-positive, irrespective of hormone receptor status, and 14.4% were TNBC (Tables III and IV).
The present study showed all 4 SNPs were in HWE (Table I). The MAFs of the SNPs in this study were more similar to those of the African sample group than those of the global population, as expected. For example, the global MAF of rs2981578 (C/T) was T=0.372, while it was T=0.078 in the African population and T=0.061 in the present study cohort (Table I).
Cases and controls did not differ in allele or genotype frequencies of the four FGFR2 SNPs. The odds ratios (ORs) (95% confidence intervals) were: rs2981582, OR=1.10 (0.98–1.25) and P=0.156; rs35054928, OR=0.95 (0.84–1.40) and P=0.604; rs2981578, OR=1.07 (0.82–1.40) and P=0.637; and rs11200014, OR=1.04 (0.87–1.25) and P=0.666 (Table V).
Although in rs11200014, the recessive genotype (AA) was associated with an increased risk of HER2-positive breast cancer (0.038) compared with HER2-negative tumors (0.014) (Table VI), the associations were not statistically significant after Bonferroni correction.
The different immunohistochemical subtypes were not associated with the FGFR2 SNPs (Table VII). In addition, ductal carcinoma, lobular carcinoma, and other breast cancer types were not associated with the FGFR2 SNPs, except for invasive lobular cancers exclusively, which had a protective allele in rs2981578 (P=0.016), whereas the more aggressive invasive ductal cancers had the risk allele (T) (Table VIII).
Immunohistochemical images of the various subtypes are shown in Fig. S1, Fig. S2, Fig. S3, Fig. S4, Fig. S5, Fig. S6. Each figure has a representative hematoxylin and eosin-stained image, and an IHC-stained image for ER, PR, HER2 and Ki67 (51). Fig. S1 shows an example of an ILC positive for hormone receptors (HRs) but negative for HER2 and with a Ki67 of 10%. Fig. S2, Fig. S3, Fig. S4, Fig. S5, Fig. S6 show IDCs. Fig. S2 is indicative of an A-like IHC, with positive HR status, no HER2 localization and low Ki67 expression. B-like IHC (Fig. S3) shows positive HR, negative HER2 and high Ki67 expression. B/HER2 (Fig. S4) is a sample that is HR-positive and HER2-positive, and HER2-like (Fig. S5) is HR-negative and HER2-positive. Fig. S6 shows a triple negative cancer, where there is no expression of ER, PR or HER2. Images were selected by a senior pathologist at the NHLS.
Discussion
There is very little data on the genetics of cancer in sub-Saharan Africa in general, and here no association between FGFR2 variants with breast cancer within the Saharan or sub-Saharan populations was shown. However, FGFR2 variants have been studied in African American populations (23). In the present study, the potential association of FGFR2 intronic SNPs with breast cancer in black southern African women was assessed. FGFR2 belongs to a tyrosine kinase receptor family that catalyzes multiple processes, including pro-survival signals, anti-apoptotic signals, cell proliferation, and cell migration (19). In GWAS studies, the rs2981582 SNP was strongly associated with breast cancer risk (17,18), particularly with ER-positive cancers [P=6×10−7]. The present candidate gene replication study however, found no association between breast cancer and rs2981582, similar to the results of Udler et al (23) who also found no significant association. In further analysis, ER, PR, HER2, and immunohistochemical tumor types luminal A, luminal B, HER2 enriched, or TNBC, were also not associated with rs2981582.
Intron 2 of FGFR2 contains putative transcription factor binding sites. Meyer et al (32) showed that the risk-associated allele (C) of rs2981578 preferentially bound FOXA1, and was able to recruit ERα to this site, and that rs35054928 preferentially bound to E2F1. E2F1 is important in the regulation of the proliferative response of breast cancer cells to estrogen. The expression of E2F1 increases with more advanced stages of breast cancer (56). Breast cancer was not associated with rs35054928 in the present study. However, rs2981578 was associated with ILC. Classic ILCs are typically of low histological grade, express ER and PR, and rarely show HER2 protein overexpression or amplification. In the present study, all ILCs were homozygous for the major allele, C/C. A defining feature of ILC is a lack of CDH1 protein expression. Ciriello et al (57) found ILC had mutations in PTEN, TBX3, and FOXA1. These suggest that ILC has mutations in FOXA1 leading to increased ER recruitment as well as increased FOXA1 binding to FGFR2. These results, however, must be taken cautiously, as only 14 ILC samples out of 989 breast cancer cases had data on rs2981578. Nevertheless, this association may be interesting to investigate in a larger cohort of patients with ILC cases.
Udler et al (23) showed a significant association between rs2981578 and ER-positive cases compared with the controls in African American women, whereas no significant associations were found in the black southern African women in the present study. Similarly, Barnholtz-Sloan et al (39) found rs11200014 to be associated with breast cancer in African American women, but no association was found between rs11200014 and breast cancer in the present study. This difference may be caused by the genetic heterogeneity of the two populations; several African American women have both European and African ancestry, and their African ancestry is predominantly West or West-Central African (1). Indeed, multiple studies have shown that African Americans are most closely related to the Yoruba or Esan groups of Nigeria (58–60), or to groups from Sierra Leone (61), which are West African states, and to a lesser extent to the people from the Gambia, also a West African state, or people from Central Africa (58,59). African Americans are less related to people from East African states such as Kenya, and to groups from Southern Africa such as the Xhosa or the San, and people from Northern Africa (41,58,61). Most African Americans are in populations that form a continuum from Europeans to West Africans (40,41,60,62).
Sub-Saharan Africans, including western, eastern and southern Africans as well as African Americans, are genetically diverse groups. Indeed, the human hereditary and health in Africa (1,2,63) consortium and malaria genomic epidemiology network (64) have shown that while there is genetic transfer between different African groups, there are distinct geographic and genomic groups. In sub-Saharan Africa, the movement of people speaking Niger-Congo languages, seems to have been from Nigeria (West Africa), through central Africa to Zambia (East Africa). From Zambia, there was a movement of people north and east (to present-day Uganda, Kenya, and Ethiopia), which makes up the Eastern Africa group. From Zambia, there was movement south to Botswana, South Africa, Namibia, and Eswatini, which make up the southern African group. These groups are genetically different from each other. The southern African group has also interacted with the Khoe and San groups, which are as distinct from people who speak Niger-Congo languages as are people from Europe (63). Additionally, the FGFR2 intron 2 block is in strong linkage disequilibrium (LD) among European populations while this LD is weaker among the African populations and thus the selected SNPs in this replication study may not be in LD with the causal variant.
Among the cases of the present study, the homozygote recessive genotype AA for rs11200014 was more prevalent in HER2-positive than in HER2-negative tumors; however, given the Bonferroni correction, the association was not statistically significant. Fernández-Noguiera et al (65) found that activation of FGFR2 increased resistance to HER2 therapy. Conversely, when FGFR2 was inactivated, HER2 activity decreased, and therapy against resistant HER2 breast cancer cells improved. Hanker et al (66) suggested that resistance to HER2 therapy may be caused by a change from an ER/HER2 signaling pathway to an FGFR2 signaling pathway. In ~25% of breast cancer cases in the South African population, patients are positive for HER2 expression (67,68).
In conclusion, GWAS studies in other populations have highlighted intron 2 of FGFR2 as a region of interest in breast cancer. In the present study, the rs2981582, rs35054928, rs2981578, and rs11200014 SNPs were investigated in samples from black South African women but found no significant association with breast cancer. Thus, it is surmised that the difference between black southern African women and African American women is caused by the genetic diversity between southern Africans and west Africans, as well as the historical influence of European ancestry in the African American population. A limitation of this study is that the scope did not allow for interrogation of environmental factors that could cause epigenetic or germline mutations, and hence affect breast cancer susceptibility in this population. This study was used to investigate the low penetrance gene, FGFR2, that was highlighted by GWAS, as has been done in other geographic regions. As such, high and medium penetrance genes such as BRCA1/2, CHEK2, or PALB2, were not interrogated in the black South African population. These studies will be performed in the future. The black South African population may have weak LD with the causal allele, and the true causal variant may not yet be defined. Some interesting findings, albeit with low numbers, are that rs2981578 is associated with invasive lobular cancer, possibly through the FOXA1 pathway, and that the recessive homozygote of rs11200014 is associated with HER2-positive breast cancer.
Supplementary Material
Supporting Data
Acknowledgements
The authors would like to thank Mr. Eric Liebenberg (National Health Laboratory Services, Johannesburg, South Africa) for his help with the presentation of the immunohistochemistry micrographs; and Ms. Confidence Makgoro (Internal Medicine, Johannesburg, South Africa); as well as Mr. Victor Shandukani, Ms. Nontlanta Mkwanazi, Ms. Mokgadi Mawela, Ms. Thandi Mtyapi, Ms. Sihle Sibiya, Ms. Yvonne Chaka and Ms. Olebogeng Mokgadi (Strengthening Oncology Research Unit, University of the Witwatersrand, Johannesburg, South Africa), for help recruiting patients and control participants.
Availability of data and materials
The datasets used and/or analyzed during the present study are available from the corresponding author on reasonable request.
Authors' contributions
TDP, CD and RABD designed the study with input from CGM and MH on GWAS, SNP selection and power calculations. TDP, CD, BPP, EJvdB, MJ, OAA, HC, SN, AIN, JSJ, PR and RABD collected and checked patient data. TDP, CD, MJ, OA and RABD collected and checked control participant data. TDP, CD, TNA and RABD performed the experiments. TDP, CD and RABD collected and cross referenced the SNP data to patients and controls. TDP, TNA, CD and RABD analyzed the data. TDP, CD and RABD confirm the authenticity of all the raw data. TDP, CD, TNA, CGM, MH, and RABD wrote the manuscript. All authors edited the manuscript. All authors read and approved the final manuscript.
Ethics approval and participant consent
The protocol used in the present study regarding human participants complied with the guidelines described in the Declaration of Helsinki. The present study was approved by the Human Research Ethics Committee (Medical) of the University of the Witwatersrand (approval no. M140980, M161116). All participants provided signed informed consent prior to being enrolled in the study.
Patient consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
AbbreviationsASIR
age standardized incidence rate
ATM
ataxia-telangiectasia mutated
BRCA1
breast cancer 1
CDH1
epithelial cadherin 1
CHBAH
Chris Hani Baragwanath Academic Hospital
CHEK2
checkpoint kinase 2
CMJAH
Charlotte Maxeke Johannesburg Academic Hospital
ER
estrogen receptor
FGFR2
fibroblast growth factor receptor 2
GWAS
genome-wide association studies
HDAC
histone deacetylase
HER2
human epidermal growth factor receptor 2
HWE
Hardy-Weinberg equilibrium
ILC
invasive lobular carcinoma
KASP
Kompetitive allele-specific PCR
LD
linkage disequilibrium
LR
lifetime risk
MAF
minor allele frequency
NHLS
National Health Laboratory Services
PALB2
partner and localizer of BRCA2
PR
progesterone receptor
PTEN
phosphatase and tensin homolog
RUNX2
runt-related transcription factor 2
SABCHO
South African Breast Cancer and HIV Outcome
SNP
single nucleotide polymorphism
TNBC
triple-negative breast cancer
TP53
tumor protein P53
YY1
Yin Yang 1
ReferencesGurdasaniDCarstensenTTekola-AyeleFPaganiLTachmazidouIHatzikotoulasKKarthikeyanSIlesLPollardMOChoudhuryAThe African Genome Variation Project shapes medical genetics in AfricaNature517327332201510.1038/nature1399725470054GurdasaniDBarrosoIZegginiESandhuMSGenomics of disease risk in globally diverse populationsNat Rev Genet20520535201910.1038/s41576-019-0153-z31235872RotimiSORotimiOASalhiaBA Review of cancer genetics and genomics studies in AfricaFront Oncol10606400202110.3389/fonc.2020.60640033659210HayatMChenWCBrandenburgJTBabb de VilliersCRamsayMMathewCGGenetic susceptibility to breast cancer in Sub-Saharan African populationsJCO Glob Oncol714621471202110.1200/GO.21.0008934623906SharmaRBreast cancer burden in Africa: Evidence from GLOBOCAN 2018J Public Health (Oxf)43763771202110.1093/pubmed/fdaa09932657321South African National Cancer Registry. National Cancer RegistryCancer in South Africa 2019https://www.nicd.ac.za/wp-content/uploads/2021/12/NCR_Path_2019_Full_Report_8dec2021.pdfMay132022HallJMLeeMKNewmanBMorrowJEAndersonLAHueyBKingMCLinkage of Early-onset familial breast cancer to chromosome 17q21Science25016841689199010.1126/science.22704822270482WoosterRNeuhausenSLMangionJQuirkYFordDCollinsNNguyenKSealSTranTAverillDLocalization of a breast cancer susceptibility gene, BRCA2, to chromosome 13q12-13Science26520882090199410.1126/science.80912318091231ShiovitzSKordeLAGenetics of breast cancer: A topic in evolutionAnn Oncology2612911299201510.1093/annonc/mdv02225605744MavaddatNAntoniouACEastonDFGarcia-ClosasMGenetic susceptibility to breast cancerMol Oncol4174191201010.1016/j.molonc.2010.04.01120542480EygelaarDvan RensburgEJJoubertFGermline sequence variants contributing to cancer susceptibility in South African breast cancer patients of African ancestrySci Rep12802202210.1038/s41598-022-04791-135039564FranciesFZWainsteinTDe LeeneerKCairnsAMurdochMNietzSCubaschHPoppeBVan MaerkenTCrombezBBRCA1, BRCA2 and PALB2 mutations and CHEK2 c.1100delC in different South African ethnic groups diagnosed with premenopausal and/or triple negative breast cancerBMC Cancer15912201510.1186/s12885-015-1913-626577449Van der MerweNCCombrinkHMENtaitaKSOosthuizenJPrevalence of clinically relevant germline BRCA variants in a large unselected South African breast and ovarian cancer cohort: A public sector experienceFront Genet13834265202210.3389/fgene.2022.83426535464868Van Der MerweNCOosthuizenJTheronMChongGFoulkesWDThe contribution of large genomic rearrangements in BRCA1 and BRCA2 to South African familial breast cancerBMC Cancer20391202010.1186/s12885-020-06917-y32375709FachalLAschardHBeesleyJBarnesDRAllenJKarSPooleyKADennisJMichailidouKTurmanCFine-mapping of 150 breast cancer risk regions identifies 191 likely target genesNat Genet525673202010.1038/s41588-019-0537-131911677ZhangHAhearnTULecarpentierJBarnesDBeesleyJQiGJiangXO'MaraTAZhaoNBollaMKGenome-wide association study identifies 32 novel breast cancer susceptibility loci from overall and subtype-specific analysesNat Genet52572581202010.1038/s41588-020-0609-232424353EastonDFPooleyKADunningAMPharoahPDPBallingerDGStruewingJPMorrisonJFieldHLubenRWarehamNGenome-wide association study identifies novel breast cancer susceptibility lociNature44710871093200710.1038/nature0588717529967HunterDJKraftPJacobsKBCoxDGYeagerMHankinsonSEWangZWelchRHutchinsonAWangJA genome-wide association study identifies alleles in FGFR2 associated with risk of sporadic postmenopausal breast cancerNat Genet39870874200710.1038/ng207517529973ZhouWYZhengHDuXLYangJLCharacterization of FGFR signaling pathway as therapeutic targets for sarcoma patientsCancer Biol Med13260268201610.20892/j.issn.2095-3941.2015.010227458533KatohYKatohMFGFR2-related pathogenesis and FGFR2-targeted therapeutics (Review)Int J Mol Med23307311200910.3892/ijmm_0000013219212647AgarwalDPinedaSMichailidouKHerranzJPitaGMorenoLTAlonsoMRDennisJWangQBollaMKFGF receptor genes and breast cancer susceptibility: Results from the Breast Cancer Association ConsortiumBr J Cancer11010881100201410.1038/bjc.2013.76924548884MountjoyESchmidtEMCarmonaMSchwartzentruberJPeatGMirandaAFumisLHayhurstJBunielloAKarimMAAn open approach to systematically prioritize causal variants and genes at all published human GWAS trait-associated lociNat Genet5315271533202110.1038/s41588-021-00945-534711957UdlerMSMeyerKBPooleyKAKarlinsEStruewingJPZhangJDoodyDRMacArthurSTyrerJPharoahPDFGFR2 variants and breast cancer risk: Fine-scale mapping using African American studies and analysis of chromatin conformationHum Mol Genet1816921703200910.1093/hmg/ddp07819223389MavaddatNMichailidouKDennisJLushMFachalLLeeATyrerJPChenTHWangQBollaMKPolygenic risk scores for prediction of breast cancer and breast cancer subtypesAm J Hum Genet1042134201910.1016/j.ajhg.2018.11.00230554720AlRaddadiRIRAlamriRJNShebliWTYFallatahEIYAlhujailyASMohamedHSAlotibiMKHFibroblast growth factor receptor 2 gene (FGFR2) rs2981582T/C polymorphism and susceptibility to breast cancer in Saudi womenSaudi J Biol Sci2861126115202110.1016/j.sjbs.2021.07.00534759736LiangHYangXChenLLiHZhuASunMWangHLiMHeterogeneity of breast cancer associations with common genetic variants in FGFR2 according to the Intrinsic Subtypes in Southern Han Chinese WomenBiomed Res Int2015626948201510.1155/2015/62694826421298HanWWooJHYuJHLeeMJMoonHGKangDNohDYCommon genetic variants associated with breast cancer in korean women and differential susceptibility according to intrinsic subtypeCancer Epidemiol Biomarkers Prev20793798201110.1158/1055-9965.EPI-10-128221415360PurnomosariDRaharjoCKalimASHerviastutiRYushanMFajarRAHarrisKKWahyonoAP21 Ser31Arg and FGFR2 rs2981582 polymorphisms as risk factors for early onset of breast cancer in Yogyakarta, IndonesiaAsian Pac J Cancer Prev2033053309201910.31557/APJCP.2019.20.11.330531759353BrignoniLCappettaMColistroVSansMArtagaveytiaNBonillaCBertoniBGenomic diversity in sporadic breast cancer in a Latin American populationGenes (Basel)111272202010.3390/genes1111127233126731ÖzgözAMutlu İçduyguFYükseltürkAŞamlIHHekİmler ÖztürkKBaşkanZLow-penetrance susceptibility variants and postmenopausal oestrogen receptor positive breast cancerJ Genet9915202010.1007/s12041-019-1174-232366738HodoǧlugilUMahleyRWTurkish population structure and genetic ancestry reveal relatedness among eurasian populationsAnn Hum Genet76128141201210.1111/j.1469-1809.2011.00701.x22332727MeyerKBO'ReillyMMichailidouKCarleburSEdwardsSLFrenchJDPrathalinghamRDennisJBollaMKWangQFine-scale mapping of the FGFR2 breast cancer risk locus: Putative functional variants differentially bind FOXA1 and E2F1Am J Hum Genet9310461060201310.1016/j.ajhg.2013.10.02624290378CampbellTMCastroMAADe SantiagoIFletcherMNCHalimSPrathalingamRPonderBAJMeyerKBFGFR2 risk SNPs confer breast cancer risk by augmenting oestrogen responsivenessCarcinogenesis37741750201610.1093/carcin/bgw06527236187LuPHankelILHostagerBSSwartzendruberJAFriedmanADBrentonJLRothmanPBColganJDThe developmental regulator protein Gon4l associates with protein YY1, co-repressor Sin3a, and histone deacetylase 1 and mediates transcriptional repressionJ Biol Chem2861831118319201110.1074/jbc.M110.13360321454521MeyerKBMaiaATO'ReillyMTeschendorffAEChinSFCaldasCPonderBAAllele-specific up-regulation of FGFR2 increases susceptibility to breast cancerPLoS Biol6e108200810.1371/journal.pbio.006010818462018OwensTWRogersRLBestSALedgerAMooneyAMFergusonAShorePSwarbrickAOrmandyCJSimpsonPTRunx2 is a novel regulator of mammary epithelial cell fate in development and breast cancerCancer Res7452775286201410.1158/0008-5472.CAN-14-005325056120ZhuXAsaSLEzzatSHistone-acetylated control of fibroblast growth factor receptor 2 intron 2 polymorphisms and isoform splicing in breast cancerMol Endocrinol2313971405200910.1210/me.2009-007119497954OguraTAzumaKSatoJKinowakiKTakayamaKITakeiwaTKawabataHInoueSOCT1 is a poor prognostic factor for breast cancer patients and promotes cell proliferation via inducing NCAPHInt J Mol Sci2211505202110.3390/ijms22211150534768935Barnholtz-SloanJSShettyPBGuanXNyanteSJLuoJBrennanDJMillikanRCFGFR2 and other loci identified in genome-wide association studies are associated with breast cancer in African-American and younger womenCarcinogenesis3114171423201010.1093/carcin/bgq12820554749BrycKDurandEYMacphersonJMReichDMountainJLThe genetic ancestry of African Americans, Latinos, and European Americans across the United StatesAm J Hum Genet963753201510.1016/j.ajhg.2014.11.01025529636ZakhariaFBasuAAbsherDAssimesTLGoASHlatkyMAIribarrenCKnowlesJWLiJNarasimhanBCharacterizing the admixed African ancestry of African AmericansGenome Biol10R141200910.1186/gb-2009-10-12-r14120025784CubaschHRuffPJoffeMNorrisSChirwaTNietzSSharmaVDuarteRBuccimazzaIČačalaSSouth African breast cancer and HIV outcomes study: Methods and baseline assessmentJ Glob Oncol3114124201710.1200/JGO.2015.00267528706996MillerSADykesDDPoleskyHFA simple salting out procedure for extracting DNA from human nucleated cellsNucleic Acids Res161215198810.1093/nar/16.3.12153344216ZhangYZengXLiuPHongRLuHJiHLuLLiYAssociation between FGFR2 (rs2981582, rs2420946 and rs2981578) polymorphism and breast cancer susceptibility: A meta-analysisOncotarget834543470201710.18632/oncotarget.1383927966449ZhangYLuHJiHLuLLiuPHongRLiYAssociation between rs11200014, rs2981579, and rs1219648 polymorphism and breast cancer susceptibility: A meta-analysisMedicine (Baltimore)96913201710.1097/MD.00000000000092461000 Genomes Project ConsortiumAutonABrooksLDDurbinRMGarrisonEPKangHMKorbelJOMarchiniJLMcCarthySMcVeanGAAbecasisGRA global reference for human genetic variationNature5266874201510.1038/nature1539326432245KasembeliANDuarteRRamsayMMosianePDickensCDix-PeekTLimouSSezginENelsonGWFogoABAPOL1 risk variants are strongly associated with HIV-associated nephropathy in black South AfricansJ Am Soc Nephrol2628822890201510.1681/ASN.201405046925788523HeCHolmeJAnthonyJSNP genotyping: The KASP assayMethods Mol Biol11457586201410.1007/978-1-4939-0446-4_724816661FitzgibbonsPLConnollyJLBoseSChenYYde BacaMEEdgertonMHayesDFHilKALesterSCSimpsonJFProtocol for the examination of resection specimens from patients with invasive carcinoma of the breastCollege of American PathologistsNorthfield (IL)2020FitzgibbonsPLDillonDAAlsabehRBermanMAHayesDFHicksDGHughesKSNofech-MozesSTemplate for reporting results of biomarker testing of specimens from patients with carcinoma of the breastArch Pathol Lab Med138595601201410.5858/arpa.2013-0566-CP24236805GoldhirschAWoodWCCoatesASGelberRDThurlimannBSennHJStrategies for subtypes-dealing with the diversity of breast cancer: Highlights of the St Gallen international expert consensus on the primary therapy of early breast cancer 2011Ann Oncol2217361747201110.1093/annonc/mdr30421709140Department of Health Republic of South AfricaClinical Guidelines for Breast Cancer Control and Management, 1–123 2018Available fromhttps://cansa.org.za/files/2019/08/DOH-Breast-Cancer-Guidelines-Final.pdfAugust312022SuvarnaSKLaytonCBancroftJDBancroft's Theory and Practice of Histological Techniques8th editionElsevier2019Dix-PeekTPhakathiBPvan den BergEJDickensCAugustineTNCubaschHNeugutAIJacobsonJSJoffeMRuffPDuarteRABDiscordance between PAM50 intrinsic subtyping and immunohistochemistry in South African women with breast cancerBreast Cancer Res Treat199112202310.1007/s10549-023-06886-336867282LewisCMGenetic association studies: Design, analysis and interpretationBrief Bioinform3146153200210.1093/bib/3.2.14612139434StenderJDFrasorJKommBChangKCNKrausWLKatzenellenbogenBSEstrogen-regulated gene networks in human breast cancer cells: Involvement of E2F1 in the regulation of cell proliferationMol Endocrinol2121122123200710.1210/me.2006-047417550982CirielloGGatzaMLBeckAHWilkersonMDRhieSKPastoreAZhangHMcLellanMYauCKandothCComprehensive molecular portraits of invasive lobular breast cancerCell163506519201510.1016/j.cell.2015.09.03326451490TishkoffSAReedFAFriedlaenderFREhretCRanciaroAFromentAHirboJBAwomoyiAABodoJMDoumboOThe genetic structure and history of Africans and African AmericansScience32410351044200910.1126/science.117225719407144MartinARGignouxCRWaltersRKWojcikGLNealeBMGravelSDalyMJBustamanteCDKennyEEHuman demographic history impacts genetic risk prediction across diverse populationsAm J Hum Genet100635649201710.1016/j.ajhg.2017.03.00428366442DaiCLVazifehMMYeangCHTachetRWellsRSVilarMGDalyMJRattiCMartinARPopulation histories of the United States revealed through fine-scale migration and haplotype analysisAm J Hum Genet106371388202010.1016/j.ajhg.2020.02.00232142644ZimmermanKDSchurrTGChenWMNayakUMychaleckyjJCQuetQMoultrieLHDiversJKeeneKLKamenDLGenetic landscape of Gullah African AmericansAm J Phys Anthropol175905919202110.1002/ajpa.2433334008864SuchestonLEBensenJTXuZSinghPKPreusLMohlerJLSuLJFonthamETRuizBSmithGJTaylorJAGenetic ancestry, self-reported race and ethnicity in African Americans and European Americans in the PCaP cohortPLoS One7e30950201210.1371/journal.pone.003095022479307ChoudhuryAAronSBotiguéLRSenguptaDBothaGBensellakTWellsGKumuthiniJShrinerDFakimYJHigh-depth African genomes inform human migration and healthNature586741748202010.1038/s41586-020-2859-733116287BusbyGBBandGSi LeQJallowMBougamaEManganoVDAmenga-EtegoLNEnimilAApinjohTNdilaCMAdmixture into and within sub-Saharan AfricaElife5e15266201610.7554/eLife.1526627324836Fernández-NogueiraPMancinoMFusterGLópez-PlanaAJaureguiPAlmendroVEnreigEMenéndezSRojoFNoguera-CastellsATumor-associated fibroblasts promote HER2-targeted therapy resistance through FGFR2 activationClin Cancer Res2614321448202010.1158/1078-0432.CCR-19-035331699826HankerABGarrettJTEstradaMVMoorePDEricssonPGKochJPLangleyESinghSKimPSFramptonGMHER2-overexpressing breast cancers amplify FGFR signaling upon acquisition of resistance to dual therapeutic blockade of HER2Clin Cancer Res2343234334201710.1158/1078-0432.CCR-16-228728381415PhakathiBCubaschHNietzSDickensCDix-PeekTJoffeMNeugutAIJacobsonJDuarteRRuffPClinico-pathological characteristics among South African women with breast cancer receiving anti-retroviral therapy for HIVBreast43123129201910.1016/j.breast.2018.12.00530550925AchilonuOJSinghENimakoGEijkemansRMJCMusengeERule-based information extraction from free-text pathology reports reveals trends in South African female breast cancer molecular subtypes and Ki67 expressionBiomed Res Int20226157861202210.1155/2022/615786135355821MayAHazelhurstSLiYNorrisSAGovindNTiklyMHonCJohnsonKJHartmannNStaedtlerFRamsayMGenetic diversity in black South Africans from SowetoBMC Genomics14644201310.1186/1471-2164-14-64424059264
Allele frequency distribution in South African, African, and global populations.
Single nucleotide polymorphism
Position
Allele
MAF Present study
MAF aSoweto
MAF bAfrican
MAF bGlobal
Hardy-Weinberg Equilibrium P-value
rs2981582
Chr10:121592803
G/A
A=0.465
A=0.452
A=0.495
A=0.404
0.09
rs35054928
Chr10:12150918
C-/CC
C-=0.426
Not available
C-=0.331
CC=0.491
0.11
rs2981578
Chr10:121580797
C/T
T=0.061
T=0.106
T=0.078
T=0.372
0.39
rs11200014
Chr10:121575416
G/A
A=0.135
A=0.112
A=0.190
A=0.313
0.50
MAF, minor allele frequencies.
Allele frequencies looking at genetic diversity in the black population from Soweto, South Africa (69);
Frequencies available from the 1000 genomes project (46). The African group consisted of African Caribbean in Barbados; African Ancestry in Southwest USA; Esan in Nigeria; Yoruba in Ibadan, Nigeria; Gambian in Western division, the Gambia; Mende in Sierra Leone; and Luhya in Webuye, Kenya.
Sequences of the FGFR2 primers for detection of SNPs.
SNP
Primer ID
Allele
Primer sequence, 5′-3′
rs2981582
2981582_A
A
GGCACCAGGTGGACTCTCCA-FAM
2981582_G
G
GCACCAGGTGGACTCTCCG-HEX
2981582_Common
TAAAACGGCAGATCCCAGCACTCAT
rs35054928
35054928_CC
CC
CTCAGAAGGGCTGTGCGCC-FAM
35054928_C
C-
TCTCAGAAGGGCTGTGCGCG-HEX
35054928_Common
GCCCTGTCCCAGAAAGCCTACAT
rs2981578
2981578_T
T
TAACCTTTCTTCCCTGCTCCAAACT-FAM
2981578_C
C
CCTTTCTTCCCTGCTCCAAACC-HEX
2981578_Common
GTTTTCTTGAAGCTTTTACCTCTATGCAAA
rs11200014
11200014_A
A
CTCCAAAAAAAGATGCACAGAGGGAAGA-FAM
11200014_G
G
CAAAAAAAGATGCACAGAGGGAAGG-HEX
11200014_Common
ACACGTGTTGGGACCAGAGAGAAAA
SNP, single nucleotide polymorphism.
Characteristics of the study participants.
Characteristic
Case
Control
P-value
Age, n (%)
n=1,001
n=1,006
Overall age in years, median (IQR)
53 (44–64)
51 (40–62)
<0.001a
18–39 years
142 (14.2%)
241 (24.0%)
40–49 years
259 (25.9%)
232 (23.1%)
50–59 years
239 (23.9%)
223 (22.2%)
≥60+ years
358 (35.8%)
308 (30.7%)
Missing
3 (0.3%)
2 (0.2%)
Cancer stage, n (%)
Stage I
38 (3.8%)
Stage II
378 (37.8%)
Stage III
439 (43.9%)
Stage IV
139 (13.9%)
Stage unavailable
7 (0.7%)
Tumor type, n (%)
Invasive ductal carcinoma
826 (82.5%)
Invasive lobular carcinoma
14 (1.4%)
Otherb, n (%)
71 (7.1%)
Tumor type unavailable
90 (9.0%)
IHC phenotype, n (%)
A-like (ER and/or PR+/HER2-/Ki67≤14%)
108 (10.8%)
B-like (ER and/or PR+/HER2-/Ki67>14%)
447 (44.7%)
B/HER2-like (ER and/or PR+/HER2+)
211 (21.1%)
HER2 (ER-/PR-/HER2+)
58 (5.8%)
TNBC (ER-/PR-/HER2-)
144 (14.4%)
IHC unavailable
10 (1.0%)
IHC equivocalc
23 (2.3%)
P<0.001.
Other tumor types included apocrine, cribriform, medullary, mesenchymal, metaplastic, mucinous, neuroendocrine, papillary, pleomorphic, squamous, tubular, anaplastic, micropapillary carcinomas.
HER2 had a value of 2+, but no fluorescence in situ hybridization was performed to confirm, or the Ki67 values were unavailable or ambiguous. ER, estrogen receptor; PR, progesterone receptor, HER2, human epidermal growth factor receptor 2; IQR, interquartile range; Ki67, proliferation marker; A-like, luminal-A; B-like, luminal B/HER2-; B/HER2-like, luminal B/HER2+; TNBC, triple-negative breast cancer.
ER, PR and HER2 expression in the breast tumors.
Receptor status
Positive (%)
Negative (%)
Equivocala (%)
Unknown (%)
ER status
753 (75.2)
237 (23.7)
11 (1.1)
PR status
639 (63.8)
347 (34.7)
15 (1.5)
HER2 status
269 (26.9)
704 (70.3)
18 (18.0)
10 (1.0)
HER2 had a value 2+, but no fluorescence in situ hybridization was performed to confirm, or the Ki67 values were unavailable or ambiguous.
Genotypic frequencies of FGFR2 variants in the case and control populations.
Allele
Case, n (freq.)
Control, n (freq.)
P-value
Odds ratioa (95% confidence interval)
rs2981582
984
996
0.347b
GG
286 (0.291)
311 (0.312)
1 (Ref)
GA
458 (0.465)
467 (0.469)
1.07 (0.87–1.31)
AA
240 (0.244)
218 (0.219)
1.22 (0.96–1.57)
P-trendd
0.156
1.10 (0.98–1.25)
GA+AA vs. GG
698 (0.709)
685 (0.688)
0.295b
1.11 (0.92–1.35)
AA vs. GG+GA
240 (0.244)
218 (0.219)
0.187b
1.18 (0.95–1.45)
AA vs. GG
240 (0.456)
218 (0.412)
0.148b
1.23 (0.96–1.57)
G
1030 (0.523)
1089 (0.547)
1 (Ref)
A
938 (0.477)
903 (0.453)
0.141b
1.11 (0.98–1.26)
rs35054928
989
996
0.848b
CC/CC
340 (0.344)
335 (0.336)
1 (Ref)
CC/C-
463 (0.468)
464 (0.466)
0.96 (0.78–1.17)
C-/C-
186 (0.188)
197 (0.198)
0.91 (0.70–1.17)
P-trendd
0.604
0.95 (0.84–1.08)
CC/C-+C-/C- vs. CC/CC)
649 (0.656)
661 (0.664)
0.727b
0.94 (0.78–1.13)
C-/C- vs. CC/CC+CC/C-)
186 (0.188)
197 (0.198)
0.583b
0.93 (0.74–1.17)
C-/C- vs. CC/CC
186 (0.354)
197 (0.370)
0.572b
0.90 (0.69–1.16)
CC
1143 (0.578)
1134 (0.569)
C-
835 (0.422)
858 (0.431)
0.585b
0.95 (0.84–1.08)
rs2981578
989
995
0.813c
CC
867 (0.877)
879 (0.883)
1 (Ref)
CT
119 (0.120)
114 (0.115)
1.06 (0.80–1.40)
TT
3 (0.003)
2 (0.002)
1.52 (0.25–9.15)
P-trendd
0.637
1.07 (0.82–1.40)
CT+TT vs. CC
122 (0.123)
116 (0.117)
0.642b
1.07 (0.81–1.40)
TT vs. CC+CT
3 (0.003)
2 (0.002)
0.686c
1.50 (0.25–9.08)
TT vs. CC
3 (0.0034)
2 (0.0023)
0.685c
1.52 (0.25–9.14)
C
1853 (0.937)
1872 (0.941)
1 (Ref)
T
125 (0.063)
118 (0.059)
0.609b
1.07 (0.82–1.39)
rs11200014
988
994
0.903b
GG
737 (0.746)
750 (0.755)
1 (Ref)
GA
231 (0.234)
224 (0.225)
1.06 (0.85–1.31)
AA
20 (0.020)
20 (0.020)
1.01 (0.54–1.90)
P-trendd
0.666
1.04 (0.87–1.25)
GA+AA vs. GG
251 (0.254)
244 (0.246)
0.659b
1.05 (0.86–1.29)
AA vs. GG+GA
20 (0.020)
20 (0.020)
0.985b
1.00 (0.53–1.87)
AA vs. GG
20 (0.026)
20 (0.026)
0.956b
1.01 (0.54–1.90)
G
1705 (0.863)
1724 (0.867)
1 (Ref)
A
271 (0.137)
264 (0.133)
0.689b
1.04 (0.87–1.25)
Adjusted for age.
Pearson χ2,
Fisher's exact test.
Non-parametric test for trend.
Association between genotypic frequencies of FGFR2 variants and immunohistochemical subtyping of patients with breast cancer.