
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
1 Centre for Molecular, Environmental, Genetic, and Analytic Epidemiology and 2 Genetic Epidemiology Laboratory, University of Melbourne, Melbourne, Victoria, Australia; 3 Genetic Epidemiology Laboratory, Queensland Institute of Medical Research, Brisbane, Queensland, Australia; and 4 Garvan Institute of Medical Research, Cancer Research Program, and 5 Department of Medicine, St. Vincent's Hospital Clinical School, University of New South Wales, Sydney, New South Wales, Australia
Requests for reprints: John Hopper, Centre for Molecular, Environmental, Genetic and Analytic Epidemiology, University of Melbourne, 2/723 Swanston Street, Carlton, Victoria 3053, Australia. Phone: 61-0383440697; Fax: 61-0393495815. E-mail: j.hopper{at}unimelb.edu.au
| Abstract |
|---|
|
|
|---|
Methods: We studied 457 women from 207 sisterhoods (104 monozygotic twins, 182 dizygotic twins, and 171 singletons). Percentage mammographic density (PMD) as well as dense area and nondense area were calculated using a computer-assisted method. We measured six single nucleotide polymorphisms from six candidate genes (COMT, HSD3B1, IGFBP3, HER2, XPD, and XRCC3). Associations between genotypes and mammographic measures were tested (a) cross-sectionally using a multivariate normal model fitted using FISHER that allowed separate correlations for monozygotic, dizygotic, and nontwin pairs and (b) within sister pairs using paired t tests.
Results: Cross-sectionally, each additional copy of the HSD3B1 Asn367Thr variant allele was associated with lower PMD (3.47% per allele; SE = 1.65; P = 0.035). Within-pair regression estimates confirmed this association. There was no evidence for an association between the mammographic density measures and any of the other variants studied.
Conclusion: We have replicated an association between a variant in the HSD3B1 gene and PMD, which suggests that HSD3B1 may be genetic determinant of mammographic density. (Cancer Epidemiol Biomarkers Prev 2007;16(7):147984)
| Introduction |
|---|
|
|
|---|
Very little is known about the genetic determinants of these mammographic density measures. A handful of studies have examined associations between genetic variants and PMD (4-8). Four of these have considered polymorphisms in the genes that code for enzymes that produce or metabolize the sex hormone estrogen. Of these, for only two polymorphisms has their been nominally significant evidence of an association with PMD, one in each of the COMT and HSD3B1 genes (4, 6, 8). Another pathway that has been investigated has been insulin-like growth factor-I, particularly the gene that controls levels of its principle binding protein, IGFBP3 (7). Other potential candidates that have not yet been investigated for their association with mammographic density include another growth factor gene, human epidermal growth factor receptor 2 (HER2), and DNA repair genes, XPD and XRCC3, all of which have a plausible link with breast cancer risk.
In this study, we have measured six polymorphisms selected from six candidate genes: COMT, HSD3B1, IGFBP3, HER2 (ERBB2), XPD (ERCC2), and XRCC3. We have examined potential associations with the three measures of mammographic densityPMD, the absolute area of dense tissue, and the absolute area of nondense tissue. We have used a twins and sisters design and have estimated both cross-sectional and within-sibling pair associations.
| Materials and Methods |
|---|
|
|
|---|
At the time of sample section, 956 subjects from 445 families had provided a blood sample; however, >60% of these subjects were monozygotic twins. Monozygotic twins were included only if an additional sister was participating. A total of 457 subjects from 207 sisterhoods (104 monozygotic twins, 182 dizygotic twins, and 171 singletons) provided a blood sample, completed a questionnaire, and had a (prediagnosis) mammogram measured (see below).
Measurement of Mammographic Density
Mammograms were obtained with the written consent of all subjects. Films were retrieved from state BreastScreen Services (80%), private clinics (12%), and participants who kept their films at home (8%). The craniocaudal views were digitized using a Lumysis 85 scanner in the Australian Mammographic Density Research Facility at the University of Melbourne. For most subjects, the right breast was measured by one reader (J.S.) using a computer-assisted thresholding technique (9). For women with a previous breast cancer in the right breast, the prediagnosis mammogram of the left breast was used, along with the left breast of the subject's sister(s). The interactive thresholding technique involves an operator selecting a gray value as a threshold to separate the image of the breast from the background, which determines breast size. A second threshold is then selected to identify the edge(s) of the mammographically dense tissue. The computer then records the number of pixels in the digitized image that lie within the defined areas. The result is an absolute measure of the total area of the breast and an absolute measure of dense area, respectively, which when subtracted gives the area of nondense tissue and when expressed as a ratio gives a percentage of mammographic density.
The films were randomized first by family into reading sets of
100, ensuring that all twins and/or sisters of the same family were measured in the same set. Within a set, the films were randomized by subject, so that twins from the original study who contributed two mammograms were viewed one after the other, but the order was unknown to the reader. This procedure has been shown in previous work to be the best method for randomization and viewing of multiple mammograms (10). The reader was blinded to all other identifying information. A 10% random sample of repeats was included in each set and between every fourth set to assess repeatability.
Data Collection
A questionnaire was administered by telephone to each participating subject and included demographic information, weight, height, smoking history, alcohol consumption, reproductive history, cessation of menstruation, use of oral contraceptives and hormone replacement therapy (HRT), and family history of cancer. A woman was defined as postmenopausal if she had a hysterectomy, both ovaries removed, or radiation; was not on HRT at the time of the mammogram and had not menstruated 12 months prior; or was on HRT at the time of the mammogram and had not menstruated 12 months prior and was not menstruating before commencing HRT. All subjects who did not fit the postmenopausal criteria were considered premenopausal.
For twins, zygosity was determined by a standard question that describes the differences between identical and nonidentical pairs. For pairs whose answers contradicted each other, or who were unsure, zygosity was determined using additional questions and methods of classifying responses that have been shown to give 95% agreement with zygosity based on blood typing in middle-aged adults (11-13).
Blood Collection
Participants were posted a blood collection kit containing tubes, a shipping container, and detailed instructions. Collection of a 27 mL blood sample was arranged with a contracted pathology laboratory. Blood samples were couriered to the Genetic Epidemiology Laboratory at the University of Melbourne within 48 h of collection. In the event that participants living within the Melbourne metropolitan area were unable to travel to a designated pathology laboratory, a trained phlebotomist visited the subject(s) at their home.
DNA Extraction and Storage and Genotyping
Blood samples were processed to generate Guthrie cards (stored at room temperature), whole blood, plasma, non lymphocytes (stored at -70°C) and monocytes stored in liquid nitrogen. DNA was extracted from whole blood aliquots using a Corbett Robotics DNA Xtractor and placed in a 96-well plate format. The XRCC3 single nucleotide polymorphism (SNP; rs861539) was genotyped by Taqman allelic discrimination. Genotyping of the COMT SNP (rs17295216) was done by direct sequencing. The other polymorphisms were genotyped using homogenous MassEXTEND matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (HSD3B1 rs1047303, IGFBP3 rs2854744, HER2 rs1801200, and XPD rs17359310).
Taqman Allelic Discrimination
The XRCC3 rs861539 SNP was genotyped using a custom-designed Taqman genotyping assay (Applied Biosystems) according to the instructions of the manufacturer and using fluorescent allele-specific probes. Briefly, 5 ng of template DNA, 2.5 µL 2x Taqman Universal PCR Master Mix, and 0.125 µL 20x SNP Genotyping Assay Mix (including primers and fluorescently labeled probes, available on request) were added to a total reaction volume of 5 µL. PCR cycling was done using an ABI Prism 7900HT Sequence Detection System under the following conditions: 95°C for 10 min followed by 40 cycles of 92°C for 15 s and 60°C for 1 min. Allelic discrimination (via end-point fluorescence) was done in 384-well format using an ABI 7900HT Sequence Detection System and the allelic discrimination analysis module (Applied Biosystems).
Sequencing
The COMT rs17295216 (or rs4680) SNP was genotyped via direct sequencing after amplification using the following primers: forward, 5'-CTGTGGCTACTCAGCTGTGC-3' and reverse, 5'-GCATGCACACCTTGTCCTT-3' (using standard PCR conditions and annealing temperature of 64°C). Direct sequencing was achieved using the forward primer, ABI Prism BigDye Termination kit, and the ABI Prism 3100 Genetic Analyzer (Applied Biosystems).
Homogenous MassEXTEND Matrix-Assisted Laser Desorption/Ionization Time-of-Flight Mass Spectrometry
PCR and homogenous MassEXTEND extension primer sequences for the four SNPs were designed using Sequenom RealSNP7 and multiplexed into three assays (two 9-plex and one 4-plex, other SNPs assessed in these assays do not form the scope of this study). A 5 µL PCR contained 2.5 ng (4-plex assay) or 5 ng (9-plex assay) of template DNA, 10x Qiagen HotStar Taq PCR buffer,8 25 mmol/L MgCl2, 25 mmol/L deoxynucleotide triphosphates, 200 nmol/L of each PCR primer, and 0.1 unit Qiagen HotStar Taq polymerase. PCR cycling conditions included 95°C for 15 min, 45 cycles of 95°C for 20 s, 56°C for 30 s, and 72°C for 1 min followed by a final extension step of 3 min at 72°C. After exonuclease cleanup, the MassEXTEND reaction was done using 1x the appropriate ACT termination mix, 600 nmol/L of each extension primer, and 0.054 unit (5-plex) or 1.25 units (9-plex) of Thermo Sequenase (Sequenom). The cycling conditions were as follows: 94°C for 2 min, 60 (5-plex) or 100 (9-plex) cycles of 94°C for 5 s, 52°C for 5 s, and 72°C for 5 s. Samples were spotted onto SpectroCHIP Bioarrays using the MassARRAY nanodispenser (Sequenom) and run on a MassARRAY Compact System (Sequenom). Results were analyzed using Sequenom MassARRAY RT software.
Statistical Analysis
Deviation of observed genotype frequencies from those expected under Hardy-Weinberg equilibrium was assessed by the usual
2 test with one degree of freedom. Polymorphisms from the candidate genes were coded as 0, 1, and 2 for homozygous wild-type, heterozygaus, and homozygous variant genotypes, respectively. Dense and nondense area measurements were log transformed.
Cross-Sectional Analysis. The cross-sectional mean was estimated for each genotype while concurrently adjusting for the potential confounders age, weight, height, number of live births, years of HRT use, pack-years of smoking, menopausal status, and years between time of mammogram and time of interview. The variance was assumed to be a constant and the covariance between sister pairs was allowed to differ according to whether they were both members of a monozygotic pair, both members of a dizygotic pair, or nontwin sisters. The effects of genotype were first parameterized in terms of the mean for the homozygous wild-type genotype and the differences between that mean and the means for each of the other two genotypes. They were also parameterized in terms of an additive genetic model, which assumed that the difference in mean between each genotype was the same, represented by a regression coefficient ß. Variables and SEs were estimated under maximum likelihood theory on the assumption that, for each sisterhood, the (adjusted) measures followed a multivariate normal distribution. Models were fitted using the statistical software package FISHER (14).
Within-Pair Analysis. The associations between within-pair differences in mammographic measures and within-pair differences in genotypes were estimated using linear regression. The mean within-pair difference was parameterized first depending on whether sisters in the same pair differed from each other by one or two alleles, and these means were compared with zero. It was also parameterized under the assumption that the mean within-pair difference in mammographic measure was twice that in pairs who differed by two alleles than in pairs who differed by one allele and represented by a regression coefficient ßW. This was achieved by regressing within-pair difference in mammographic measure against within-pair difference in number of alleles and fitting a straight line through the origin. Note that pairs who did not differ in genotype did not contribute to the within-pair analysis of that genotype. All within-pair linear regression models were adjusted for the within-pair differences in other measured covariates. Models were fitted using PROC MIXED in Statistical Analysis System (version 8.2; SAS Institute), with a family indicator term included to take into account that some families contributed multiple pairs.
| Results |
|---|
|
|
|---|
|
|
|
Within-Pair Analysis
Table 4
shows the results of the within-pair regression analyses. Estimates replicate the weak association between log dense area and the HSD3B1 genotype found in the cross-sectional analysis. Sisters with one copy of the variant allele had 0.26 cm2 less log dense area than those with none (SE = 0.14; P = 0.06).
|
The magnitude of the within-pair regression coefficient between PMD and the HSD3B1 genotype was stronger than that from the cross-sectional estimates given in Table 3 (ß = 5.85; SE = 2.07; Ptrend = 0.005). There was no nominally significant evidence that any of the polymorphisms in the COMT, IGFBP3, XPD, and XRCC3 genes was associated with any of the mammographic measures (all Ptrend > 0.2).
| Discussion |
|---|
|
|
|---|
The association we found between the HSD3B1 gene and PMD is similar to that previously reported by Haiman et al. (4). They found that Caucasian women (with breast cancer) homozygous for the Thr allele of HSD3B1 had 5.1 percentage points less PMD than those homozygous for the Asn allele. Interestingly, they found this effect to be in the opposite direction for African-American women (absolute difference +19.7% Thr/Thr versus Asn/Asn; P = 0.02), although there were only four subjects with the Thr/Thr genotype. More than 95% of our study sample indicated their ethnicity as European, of which only seven women indicated that they also came from North African descent. The other 5% of women were of Pacific Island descent. Little is known about the function of this variant and its relationship to steroid hormone levels, let alone breast cancer risk, but our apparent replication of its association with PMD merits further investigation.
Four studies have investigated the relationship between the COMT Val158Met polymorphism and mammographic density, particularly for premenopausal women. Two found little evidence of an association between PMD and the COMT gene (4, 5). However, in the 2002 study, all of the subjects had breast cancer, and therefore the distribution of the genotypes may not have been representative of the general population. In the 2003 study, the premenopausal subgroup had only 94 subjects and may not be sufficiently powered to rule out an association. The two remaining studies found that decreased PMD was associated with the low-activity homozygous genotype, Met/Met, in premenopausal women (6, 8). However, the direction of these associations is opposite to expectation, as COMT is involved in the detoxification of harmful estrogen metabolites thought to increase breast cancer risk. It was therefore hypothesized that the low activity of the COMT gene results in increased mammographic density and increased breast cancer risk (8). We found little evidence of an association between the COMT gene and PMD in the cross-sectional analysis (including separate examination of the premenopausal subgroup) or the within-pair analysis. We pooled the mean estimates across all five studies for each COMT genotype, weighted by the inverse of the variance of each estimate. The overall pooled estimates of mean PMD for each of the genotypes were 40.8 (SE = 2.1), 40.3 (SE = 3.2), and 37.4 (SE = 4.0) for Val/Val, Val/Met, and Met/Met genotypes, respectively, showing little evidence of an association.
To date, the commonly reported Val655Ile polymorphism in the HER2 gene has not previously been investigated for an association with mammographic density. We found, at best, only weak evidence of a positive association between nondense area and the variant allele. One of the activities of the HER2 gene is the control of cell amplification and it may play a role in the development of breast cancer, as overexpression of this gene is found in >20% of human breast tumors (15). This would not explain why the variant allele would be associated with increased nondense or fatty tissue.
Insulin-like growth factor-I is a growth hormone that is produced in the breast stroma and promotes the proliferation of normal and malignant breast epithelial cells. Its principal binding protein, IGFBP3, slows breast cell proliferation and promotes apoptosis. Lai et al. (7) found a positive association of the A allele of the single base pair polymorphism at position 202 relative to the CAP site in the IGFBP3 gene (A-202C) with serum levels of IGFBP3 and with PMD in premenopausal women. Whereas they found no association between serum levels of IGFBP3 and PMD, three other studies reported a negative association between levels of IGFBP3 and PMD in premenopausal women (16-18). This is consistent with the biological theory that IGFBP3 slows breast cell proliferation and promotes apoptosis, inferring that increased levels of IGFBP3 should be protective against breast cancer. We found no evidence that the IGFBP3 A-202C polymorphism, which has recently been shown to be associated with both levels of IGFBP3 and breast cancer risk (19), was associated with any of the mammographic measures.
Two coding DNA repair gene polymorphisms were genotyped because there is a plausible link with breast cancer risk. There was no evidence that the Lys751Gln polymorphism from the XPD gene or the Thr241Met polymorphism from the XRCC3 gene was associated with any of the mammographic measures.
A potential limitation of this study is its modest sample size so that small real effects would not have been detectable. For example, for PMD, the SE of the linear effect term (ß) ranged from
1.5 to <3 across the cross-sectional and within-pair analysis, so only effects of about 4 to 7 units per allele or more would have been detectable with 80% power. Similar post hoc power calculations can be calculated for the other mammographic measures, and the size of these compared with the means of these measures presented in Table 1. Therefore, in general, we cannot exclude effects of <10% of the mean per allele.
In conclusion, we have replicated an association between the HSD3B1 gene and PMD. The discovery of genetic determinants of mammographic density is of considerable importance as it is a strong and highly heritable risk factor for breast cancer and may identify genetic variants associated with the disease or, at least, may lead to the discovery of new genetic pathways of relevance to breast cancer control. Clearly, much larger studies of mammographic density using candidate genes or genome-wide scans are needed.
| Acknowledgments |
|---|
| Footnotes |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Received 2/ 4/07; revised 4/29/07; accepted 5/10/07.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Kataoka, A. Antoniou, R. Warren, J. Leyland, J. Brown, T. Audley, and D. Easton Genetic Models for the Familial Aggregation of Mammographic Breast Density Cancer Epidemiol. Biomarkers Prev., April 1, 2009; 18(4): 1277 - 1284. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |