Abstract
Several variants in the TSHR and RET signaling pathways genes have been reported to be related to cancer risk. We hypothesized that polymorphic variants in these genes are associated with the risk of papillary thyroid cancer. A nested case-control study was conducted within the U.S. Radiologic Technologists cohort. Eligible validated papillary thyroid cancer cases (n = 167) and frequency-matched (by sex and birth year) controls (n = 491) donated blood for analysis. There were no statistically significant associations between papillary thyroid cancer and 10 selected polymorphic variants in analyses of men and women combined. A borderline significant increasing risk was found for RET G691S (Ptrend = 0.05) and was especially pronounced among young women. For women under 38 years (the median age at diagnosis), the odds ratios were 2.1 (95% confidence interval, 1.2-3.7) for those heterozygous for the RET G691S polymorphism and 3.7 (95% confidence interval, 1.1-11.8) for those who were homozygous (Ptrend = 0.001). Our data provide limited evidence that TSHR- and RET-related genes are related to papillary thyroid cancer risk. (Cancer Epidemiol Biomarkers Prev 2007;16(1):174–7)
Introduction
Thyroid cancer is the most prevalent endocrine malignancy, and papillary thyroid cancer is the most common histologic type. In the United States, ∼90% of all reported thyroid cancer cases are of papillary histology (1). The age-adjusted incidence is low compared with more common cancers but has reportedly increased during the 1990s in several developed countries (1-4). Thyroid cancer occurs more frequently in women, in ∼3:1 ratio, with a median age at diagnosis of ∼46 years (1). The etiology of thyroid cancers remains largely unknown. Some factors, such as height, weight, diet, benign thyroid diseases, reproductive outcomes, and residing in areas of iodine deficiency or excess, have been reported as possible risk factors for thyroid cancer (5-7), but the only established risk factors are a previous history of benign thyroid adenoma and ionizing radiation exposure during childhood (5, 8). Several genetic factors related to thyroid tumorigenesis have been reported (9), but many studies have been limited by small numbers.
The thyrotropin receptor (TSHR) is a key protein in the control of thyroid function and a major thyroid autoantigen. Mutations of the TSHR gene are associated with benign hyperfunctioning thyroid nodules and adenomas (9) and are possibly related to thyroid cancer (10). Rearranged during transfection (RET) gene products are involved in kinase signaling pathways that control the transcription of genes that regulate cell growth and differentiation. The RET gene is a potent oncogene that is associated with several cancers, including thyroid cancer (9, 11-13). We hypothesized that polymorphic variants in TSHR- and RET-related genes might be associated the risk of papillary thyroid cancer. We tested this hypothesis in a nested case-control study within a large cohort of radiologic technologists.
Materials and Methods
The study population was the U.S. Radiologic Technologists cohort established in the early 1980s to study cancer and other health effects. The details of the U.S. Radiologic Technologists cohort have been reported previously (14). Briefly, the initial cohort, including 146,022 persons (73% women), were radiologic technologists who were certified by the American Registry of Radiologic Technologists for ≥2 years during 1926 to 1982 and resided in the United States. Two postal surveys were administered to this cohort, the first during 1984 to 1989 and the second during 1995 to 1998. The surveys included questions related to work history, health outcomes, and various environmental risk factors. Approximately 110,000 technologists responded to the first and/or second survey.
Thyroid cancer cases (n = 334) were ascertained through self-reports. Living participants reporting thyroid cancer (n = 313) were contacted to obtain consent to acquire medical records to validate the diagnosis. If the participant declined to release medical records (n = 66), or if records verifying the diagnosis could not be obtained (n = 66), the case was excluded from the study. In total, 181 validated papillary thyroid cancer cases were targeted for recruitment, and 167 of them donated blood for analysis. Controls were randomly selected from a pool of controls recruited for ongoing case-control studies of breast, testis, and early-onset cancer within the U.S. Radiologic Technologists cohort. Three controls per case were frequency matched by sex and year of birth (in 5-year strata). In total, 491 control blood samples were included for analysis.
The selection of marker single nucleotide polymorphisms (SNP) was based on review of the published literature and other public databases to cover a part of the TSHR or RET signaling pathways. Genetic variants of interest were EPAC G332S, GFRA1 −193C>G, GFRA3 IVS7+39G>A, RET A45A, RET G691S, RET L769L, RET S686S, TSHR D727E, TSHR IVS1+8651A>G, and TSHR N187N. Samples were genotyped for SNPs using standard Taqman assays, as described previously (15, 16). Primer and probe sequences are available upon request to J.P. Struewing.4
Statistical Analysis
All analyses were done using the Intercooled Stata 9.1 statistical software (Stata Co., College Station, TX). SNP frequencies were assessed for deviation from Hardy-Weinberg equilibrium in the controls. Associations between SNPs and risk of papillary thyroid cancer were assessed using unconditional logistic regression analyses. Odds ratios and 95% confidence intervals were calculated as an estimate of the relative risk. We assumed a codominant relationship as numbers permitted; dominant analyses were done when homozygote variant frequencies were <1%. The analyses were adjusted for sex, age, race, smoking habits, exposure to ionizing radiation (including occupational exposure and radiation therapy to the head and neck), year began working as a radiologic technologist, years worked as a radiologic technologist, previous history of benign thyroid conditions, personal and family history of cancer, and age at diagnosis; none of these factors influenced the results (data not shown). Because of the limited effect of the covariates, only the matching variables (age and sex) were included in the final model.
To test any differences between age at diagnosis and gender, we stratified the data by sex and age using two age groups with the cutoff point at the median age of diagnosis.
Results
Distributions for selected variables are shown for 167 papillary cancer cases and 491 controls in Table 1. There were no differences in age at blood collection and race between cases and controls, but the proportion of men among cases (13%) was higher compared with controls (7%) because of a higher nonresponse among male controls. The median age of papillary thyroid cancer diagnosis was 38 years.
. | Cases (n = 167) . | Controls (n = 491) . | P* . | |||
---|---|---|---|---|---|---|
Sex | ||||||
Men | 22 (13%) | 32 (7%) | 0.01 | |||
Women | 145 (87%) | 459 (93%) | ||||
Age at diagnosis (y) | ||||||
20-29 | 36 (22%) | NA | NA | |||
30-39 | 62 (37%) | NA | ||||
40-69 | 69 (41%) | NA | ||||
Mean age (y) at blood collection (95% confidence interval) | 56 (55-57) | 57 (56-58) | 0.13 | |||
Race | ||||||
White | 163 (98%) | 474 (97%) | 0.50 | |||
Other | 4 (2%) | 17 (3%) | ||||
History of benign thyroid nodules | ||||||
Yes | 36 (23%) | 3 (1%) | <0.01 | |||
No | 124 (77%) | 455 (99%) | ||||
History of other thyroid conditions† | ||||||
Yes | 27 (18%) | 42 (10%) | 0.11 | |||
No | 122 (82%) | 371 (90%) | ||||
Personal history of cancer | ||||||
Yes | 8 (5%) | 27 (6%) | 0.66 | |||
No | 159 (95%) | 448 (94%) | ||||
Family history of cancer | ||||||
Yes | 72 (47%) | 236 (54%) | 0.13 | |||
No | 82 (53%) | 203 (46%) | ||||
Family history of thyroid cancer | ||||||
Yes | 5 (3%) | 7 (1%) | 0.20 | |||
No | 158 (97%) | 468 (99%) | ||||
Radiation therapy to the head and neck | ||||||
Yes | 9 (6%) | 3 (1%) | <0.01 | |||
No | 151 (94%) | 455 (99%) | ||||
Year began working as a radiologic technologist | ||||||
1930-1959 | 42 (26%) | 118 (24%) | 0.61 | |||
1960-1969 | 50 (31%) | 169 (35%) | ||||
1970-1982 | 72 (44%) | 200 (41%) | ||||
Years worked as a radiologic technologist (y) | ||||||
0-5 | 67 (43%) | 148 (30%) | <0.01 | |||
6-10 | 35 (22%) | 197 (40%) | ||||
>10 | 55 (35%) | 146 (30%) | ||||
Mean thyroid radiation dose in mGy (95% confidence interval)‡ | 43 (34-52) | 46 (41-50) | 0.57 |
. | Cases (n = 167) . | Controls (n = 491) . | P* . | |||
---|---|---|---|---|---|---|
Sex | ||||||
Men | 22 (13%) | 32 (7%) | 0.01 | |||
Women | 145 (87%) | 459 (93%) | ||||
Age at diagnosis (y) | ||||||
20-29 | 36 (22%) | NA | NA | |||
30-39 | 62 (37%) | NA | ||||
40-69 | 69 (41%) | NA | ||||
Mean age (y) at blood collection (95% confidence interval) | 56 (55-57) | 57 (56-58) | 0.13 | |||
Race | ||||||
White | 163 (98%) | 474 (97%) | 0.50 | |||
Other | 4 (2%) | 17 (3%) | ||||
History of benign thyroid nodules | ||||||
Yes | 36 (23%) | 3 (1%) | <0.01 | |||
No | 124 (77%) | 455 (99%) | ||||
History of other thyroid conditions† | ||||||
Yes | 27 (18%) | 42 (10%) | 0.11 | |||
No | 122 (82%) | 371 (90%) | ||||
Personal history of cancer | ||||||
Yes | 8 (5%) | 27 (6%) | 0.66 | |||
No | 159 (95%) | 448 (94%) | ||||
Family history of cancer | ||||||
Yes | 72 (47%) | 236 (54%) | 0.13 | |||
No | 82 (53%) | 203 (46%) | ||||
Family history of thyroid cancer | ||||||
Yes | 5 (3%) | 7 (1%) | 0.20 | |||
No | 158 (97%) | 468 (99%) | ||||
Radiation therapy to the head and neck | ||||||
Yes | 9 (6%) | 3 (1%) | <0.01 | |||
No | 151 (94%) | 455 (99%) | ||||
Year began working as a radiologic technologist | ||||||
1930-1959 | 42 (26%) | 118 (24%) | 0.61 | |||
1960-1969 | 50 (31%) | 169 (35%) | ||||
1970-1982 | 72 (44%) | 200 (41%) | ||||
Years worked as a radiologic technologist (y) | ||||||
0-5 | 67 (43%) | 148 (30%) | <0.01 | |||
6-10 | 35 (22%) | 197 (40%) | ||||
>10 | 55 (35%) | 146 (30%) | ||||
Mean thyroid radiation dose in mGy (95% confidence interval)‡ | 43 (34-52) | 46 (41-50) | 0.57 |
NOTE: Data for unknown categories are not shown.
Abbreviation: NA, not available.
Based on χ2 testing or t-test (for means).
Including goiter, hyperthyroidism, hypothyroidism, and thyroiditis.
Reconstructed cumulative dose in mGy truncated 5 years before diagnosis.
Table 2 shows that there were no significant deviations from Hardy-Weinberg equilibrium. The associations between papillary thyroid cancer and polymorphic variants in TSHR- and RET-related genes are shown in Table 2. The results were similar for the univariate analysis and the age- and sex-adjusted analysis. There were no statistically significant associations between papillary thyroid cancer and several polymorphic variants in the overall analysis; however, an association was indicated for RET G691S with a borderline significant trend and potentially also for TSHR N187N, which displayed a borderline significant odds ratio for heterozygous carriers (Ptrend = 0.11).
Gene/dbSNP ID* . | PHWE† . | Genotype* . | Cases (n = 167) . | Controls (n = 491) . | Crude OR (95% CI) . | Adjusted OR‡ (95% CI) . |
---|---|---|---|---|---|---|
EPAC G332S, rs12422983 | 0.07 | GG | 121 | 352 | 1.0 | 1.0 |
GS | 34 | 114 | 0.9 (0.6-1.3) | 0.9 (0.6-1.3) | ||
SS | 4 | 16 | 0.7 (0.2-2.2) | 0.7 (0.2-2.2) | ||
Ptrend = 0.41 | ||||||
GFRA1 −193C>G, rs not assigned | 0.37 | CC | 152 | 439 | 1.0 | 1.0 |
CG and GG | 9 | 38 | 0.7 (0.3-1.5) | 0.7 (0.3-1.5) | ||
GFRA3 IVS7 +39G>A, rs not assigned | 0.68 | GG | 131 | 416 | 1.0 | 1.0 |
GA and AA | 29 | 69 | 1.3 (0.8-2.2) | 1.4 (0.8-2.2) | ||
RET A45A G>A, rs1800858 | 0.90 | GG | 95 | 276 | 1.0 | 1.0 |
GA | 51 | 178 | 0.8 (0.6-1.2) | 0.8 (0.6-1.2) | ||
AA | 10 | 30 | 1.0 (0.5-2.1) | 1.0 (0.5-2.1) | ||
Ptrend = 0.52 | ||||||
RET G691S, rs1799939 | 0.48 | GG | 100 | 344 | 1.0 | 1.0 |
GS | 51 | 128 | 1.4 (0.9-2.0) | 1.4 (0.9-2.0) | ||
SS | 8 | 15 | 1.8 (0.8-4.5) | 1.9 (0.8-4.6) | ||
Ptrend = 0.05§ | ||||||
RET L769L T>G, rs1800861 | 0.45 | TT | 100 | 282 | 1.0 | 1.0 |
TG | 49 | 177 | 0.8 (0.5-1.2) | 0.8 (0.5-1.2) | ||
GG | 9 | 23 | 1.1 (0.5-2.5) | 1.1 (0.5-2.6) | ||
Ptrend = 0.48 | ||||||
RET S836S C>T, rs1800862 | 0.33 | CC | 148 | 442 | 1.0 | 1.0 |
CT | 13 | 42 | 0.9 (0.5-1.8) | 0.9 (0.5-1.7) | ||
TSHR D727E, rs1991517 | 0.74 | DD | 142 | 414 | 1.0 | 1.0 |
DE and EE | 18 | 69 | 0.8 (0.4-1.3) | 0.8 (0.4-1.3) | ||
TSHR IVS1 +8651A>G, rs not assigned | 0.70 | AA | 104 | 318 | 1.0 | 1.0 |
AG | 46 | 149 | 0.9 (0.6-1.4) | 1.0 (0.6-1.4) | ||
GG | 9 | 15 | 1.8 (0.8-4.3) | 1.9 (0.8-4.5) | ||
Ptrend = 0.54 | ||||||
TSHR N187N T>C, rs2075179 | 0.95 | TT | 128 | 351 | 1.0 | 1.0 |
TC | 30 | 117 | 0.7 (0.5-1.1) | 0.7 (0.4-1.0) | ||
CC | 2 | 9 | 0.6 (0.1-2.9) | 0.6 (0.1-2.9) | ||
Ptrend = 0.11 |
Gene/dbSNP ID* . | PHWE† . | Genotype* . | Cases (n = 167) . | Controls (n = 491) . | Crude OR (95% CI) . | Adjusted OR‡ (95% CI) . |
---|---|---|---|---|---|---|
EPAC G332S, rs12422983 | 0.07 | GG | 121 | 352 | 1.0 | 1.0 |
GS | 34 | 114 | 0.9 (0.6-1.3) | 0.9 (0.6-1.3) | ||
SS | 4 | 16 | 0.7 (0.2-2.2) | 0.7 (0.2-2.2) | ||
Ptrend = 0.41 | ||||||
GFRA1 −193C>G, rs not assigned | 0.37 | CC | 152 | 439 | 1.0 | 1.0 |
CG and GG | 9 | 38 | 0.7 (0.3-1.5) | 0.7 (0.3-1.5) | ||
GFRA3 IVS7 +39G>A, rs not assigned | 0.68 | GG | 131 | 416 | 1.0 | 1.0 |
GA and AA | 29 | 69 | 1.3 (0.8-2.2) | 1.4 (0.8-2.2) | ||
RET A45A G>A, rs1800858 | 0.90 | GG | 95 | 276 | 1.0 | 1.0 |
GA | 51 | 178 | 0.8 (0.6-1.2) | 0.8 (0.6-1.2) | ||
AA | 10 | 30 | 1.0 (0.5-2.1) | 1.0 (0.5-2.1) | ||
Ptrend = 0.52 | ||||||
RET G691S, rs1799939 | 0.48 | GG | 100 | 344 | 1.0 | 1.0 |
GS | 51 | 128 | 1.4 (0.9-2.0) | 1.4 (0.9-2.0) | ||
SS | 8 | 15 | 1.8 (0.8-4.5) | 1.9 (0.8-4.6) | ||
Ptrend = 0.05§ | ||||||
RET L769L T>G, rs1800861 | 0.45 | TT | 100 | 282 | 1.0 | 1.0 |
TG | 49 | 177 | 0.8 (0.5-1.2) | 0.8 (0.5-1.2) | ||
GG | 9 | 23 | 1.1 (0.5-2.5) | 1.1 (0.5-2.6) | ||
Ptrend = 0.48 | ||||||
RET S836S C>T, rs1800862 | 0.33 | CC | 148 | 442 | 1.0 | 1.0 |
CT | 13 | 42 | 0.9 (0.5-1.8) | 0.9 (0.5-1.7) | ||
TSHR D727E, rs1991517 | 0.74 | DD | 142 | 414 | 1.0 | 1.0 |
DE and EE | 18 | 69 | 0.8 (0.4-1.3) | 0.8 (0.4-1.3) | ||
TSHR IVS1 +8651A>G, rs not assigned | 0.70 | AA | 104 | 318 | 1.0 | 1.0 |
AG | 46 | 149 | 0.9 (0.6-1.4) | 1.0 (0.6-1.4) | ||
GG | 9 | 15 | 1.8 (0.8-4.3) | 1.9 (0.8-4.5) | ||
Ptrend = 0.54 | ||||||
TSHR N187N T>C, rs2075179 | 0.95 | TT | 128 | 351 | 1.0 | 1.0 |
TC | 30 | 117 | 0.7 (0.5-1.1) | 0.7 (0.4-1.0) | ||
CC | 2 | 9 | 0.6 (0.1-2.9) | 0.6 (0.1-2.9) | ||
Ptrend = 0.11 |
NOTE: Totals for variables are not equal because of missing information.
Abbreviations: HWE, Hardy-Weinberg equilibrium; OR, odds ratio; 95% CI, 95% confidence interval.
Amino acids and their symbols: A, alanine; G, glycine; D, aspartic acid; E, glutamic acid; S, serine. If there was no amino acid change, the “wild-type” nucleotide is underlined.
P for the deviations from expectation under the assumptions of Hardy-Weinberg equilibrium in controls.
Adjusted for sex and age.
P > 0.05.
All gender-specific analyses yielded similar results as for the combined analysis (data not shown). For women, the age-stratified odds ratios for RET G691S among those who were under age 38 years at diagnosis were 2.1 (95% confidence interval, 1.2-3.7) for heterozygous carriers and 3.7 (95% confidence interval, 1.1-11.8) for homozygous carriers (Ptrend = 0.001). The odds ratios for those who were 38 years or older were 1.2 (95% confidence interval, 0.6-2.3) and 0.8 (95% confidence interval, 0.2-4.2), respectively (Ptrend = 0.8). The same age-stratified analysis for men did not indicate any associations between the SNP studied and papillary thyroid cancer risk (data not shown). Stratifying by age at diagnosis did not modify the risks for any other SNPs (data not shown).
Discussion
The results of this study provide limited support for the hypothesis that polymorphic variants in TSHR- or RET-related genes play a role in the etiology of papillary thyroid cancer. We observed a suggestive decreased risk for papillary thyroid cancer with one SNP in the TSHR pathway (N187N), and we detected a positive association for one SNP in the RET pathway (G691S). The results for RET G691S were especially pronounced among younger women. We have limited support for our findings. Ho et al. (12) reported an association between differentiated thyroid cancer and RET S836S, but not with RET G691S, and Lesueur et al. (17) observed an association between papillary thyroid cancer and RET A45A and RET L769L, findings that we were not able to confirm in this study. A previous study found an association with the RET G691S and medullary thyroid cancer (18). This is, however, a histologic type that shares few features with papillary thyroid cancer.
It is likely that there are different etiologic factors for specific histologic types of thyroid cancer, which may be obscured in a study with multiple histologic types. To increase etiologic homogeneity, we included only confirmed papillary thyroid cancers in our study. This study is, to our knowledge, one of the largest analyzing genotypes in association with papillary thyroid cancer risk but is still limited by a relatively small number of cases; thus, the findings need to be carefully interpreted. The small number of men in the study limits the gender comparison, but it is possible that there is a difference in thyroid cancer etiology between men and women. The fact that thyroid cancer is more common among women may reflect differences in etiology (1).
Selection bias is universally a concern in case-control studies. The exclusion of cases that could not be verified by medical records and participant refusals could potentially have introduced selection bias in our study. The magnitude and direction of this potential bias is difficult to quantify; however, previous analyses that evaluated respondent characteristics in this and other study populations suggest the effect would be small (19). Another concern in genetic analysis is multiple statistical comparisons. Our strategy was to limit the number of SNPs to be tested for associations. The majority of SNPs were chosen based on previous literature indicating an association with several cancers, including thyroid cancer (9, 11-13, 17, 18).
Genetic factors often play an important role in disease development in early life rather than later. Our finding of a stronger association among younger women for RET G691S and thyroid cancer supports this idea. The median age at diagnosis for thyroid cancer in the United States is 46 years (1), which is relatively young compared with many other cancers, and suggests that genetic factors may play an important role with respect to risk.
All control participants in this study were also controls in other nested case-control genetic studies within the cohort and were used for efficiency reasons. The controls were randomly chosen from the pool of controls with blood samples available, and it is unlikely that genotype was related to the chance of being selected. We therefore have no reason to believe that the control selection could affect the estimated variables in any of the studies. Furthermore, most of the allele frequencies among the controls were similar to other populations (12).5
In conclusion, our data provide limited support that RET- and TSHR-related genes are associated with papillary thyroid cancer risk. We found a possible association for one SNP in the RET pathway among younger women, but the relatively small sample size suggests caution in interpretation. Thyroid cancer is a rare disease, and it is difficult to obtain biological samples from large numbers of thyroid cancer cases in a defined population, especially for analyses by histologic type. To more definitively assess risks associated with genetic polymorphisms, pooled analyses by combining study populations should be considered.
Grant support: Intramural Research Program of the Division of Cancer Epidemiology and Genetics and Center for Cancer Research, Department of Health and Human Services, National Cancer Institute, NIH grants NO1-CP-51016, NO2-CP-81121, NO2-CP-81005, and NO1-CP-15673.
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.
Acknowledgments
We thank the radiologic technologists who participated in the U.S. Radiologic Technologists Study, Jerry Reid of the American Registry of Radiologic Technologists for continued support of this study, Diane Kampa and Allison Iwan (University of Minnesota) for data collection and study coordination, Laura Bowen (Information Management Services, Inc., Silver Spring, MD) for data management and computing, and Chris McClure (Research Triangle International, Inc., Research Triangle Park, NC) for tracing and data management.