Diet quality, common genetic polymorphisms, and bladder cancer risk in a New England population-based study.

PURPOSE: We examined the interaction between common genetic bladder cancer variants, diet quality, and bladder cancer risk in a population-based case-control study conducted in New England. METHODS: At the time of enrollment, 806 bladder cancer cases and 974 controls provided a DNA sample and completed a diet history questionnaire. Diet quality was assessed using the 2010 Alternate Healthy Eating Index (AHEI-2010) score. Single nucleotide polymorphisms (SNPs) reported in genome-wide association studies to be associated with bladder cancer risk were combined into a polygenic risk score and also examined individually for interaction with the AHEI-2010. Adjusted odds ratios (OR) and 95% confidence intervals (CI) were calculated using logistic regression. RESULTS: A 1-standard deviation increase in polygenic risk score was associated with higher bladder cancer risk (OR, 1.34; 95% CI 1.21-1.49). Adherence to the AHEI-2010 was not associated with bladder cancer risk (OR, 0.99; 95% CI 0.98-1.00) and the polygenic risk score did not appear to modify the association between the AHEI-2010 and bladder cancer risk. In single-SNP analyses, rs8102137 (bladder cancer risk allele, C) modified the association between the AHEI-2010 total score and bladder cancer risk, with the strongest evidence for the AHEI-2010 long chain fat guideline (OR for TT, 0.92; 95% CI 0.87-0.98; OR for CT, 1.02; 95% CI 0.96-1.08; OR for CC, 1.03; 95% CI 0.93-1.14; p for interaction, 0.02). CONCLUSIONS: In conclusion, rs8102137 near the cyclin E1 gene ( CCNE1 ) may be involved in gene-diet interactions for bladder cancer risk.

Bladder cancer risk associated with family history of cancer.

Twin studies suggest a familial aggregation of bladder cancer, but elements of this increased familial risk of bladder cancer are not well understood. To characterize familial risk of bladder cancer, we examined the relationship between family history of bladder and other types of cancer among first-degree relatives and risk of bladder cancer in 1193 bladder cancer cases and 1418 controls in a large population-based case-control study. Multivariate logistic regression was used to estimate odds ratios (ORs) and 95% confidence intervals (CIs) for the association between family history of bladder cancer (defined as at least one first-degree family member with bladder cancer or a cancer of any other site). We also evaluated cancer aggregation of specific sites in family members. Participants with a first-degree relative with bladder cancer had nearly double the risk of bladder cancer (OR = 1.8, 95% CI 1.2-2.9) as those without a family history of bladder cancer. Risk was increased for having a sibling with bladder cancer (OR = 2.6, 95% CI 1.3-5.3) compared to no siblings with cancer. Bladder cancer risk was elevated when participants reported a first-degree relative with a history of female genital cancer (OR = 1.5, 95% CI 1.1-2.1), melanoma (OR = 1.9, 95% CI 1.02-3.6), and tobacco-associated cancer (OR = 1.3, 95% CI 1.06-1.6). These findings add to evidence of a familial predisposition to bladder cancer. Clarification of the aggregation of bladder cancer in families and with other cancer sites will be of interest as many loci and common polymorphisms related to bladder cancer have yet to be identified in large genomic studies.

Occupational insecticide exposure and risk of non-Hodgkin lymphoma: A pooled case-control study from the InterLymph Consortium.

Evidence for the human health effects of pesticides is needed to inform risk assessment. We studied the relationship between occupational insecticide use and risk of non-Hodgkin lymphoma (NHL) by pooling data from nine case-control studies participating in the InterLymph Consortium, including 7909 cases and 8644 controls from North America, the European Union and Australia. Insecticide use was coded using self-report or expert assessment, for insecticide groups (eg, organophosphates, pyrethroids) and active ingredients (eg, malathion, permethrin). Associations with insecticides were estimated using logistic regression to produce odds ratios (ORs) and 95% confidence intervals (CI) for all NHL and NHL subtypes, with adjustment for study site, demographic factors and use of other pesticides. Occupational insecticide use, overall, was not associated with risk of NHL. Use of organophosphate insecticides was associated with increased risk of all NHL and the subtype follicular lymphoma, and an association was found with diazinon, in particular (ever use: OR = 2.05, 95%CI: 1.24-3.37). The carbamate insecticide, carbaryl, was associated with risk of all NHL, and the strongest associations were found with T-cell NHL for ever-use (OR = 2.44, 95%CI: 1.13-5.28) and longer duration (>8 years vs never: OR = 2.90, 95%CI: 1.02-8.25). There was no association of NHL with other broad groups of insecticides, including organochlorines and pyrethroids, and some inverse associations were estimated in relation to historical DDT use. Our findings contribute to the totality of evidence available to help inform risk decisions by public health and regulatory agencies of importance given continued, widespread use of organophosphate and carbamate insecticides.

Ingested Nitrate and Nitrite and Bladder Cancer in Northern New England.

BACKGROUND: N-nitroso compounds are hypothesized human bladder carcinogens. We investigated ingestion of N-nitroso compound precursors nitrate and nitrite from drinking water and diet and bladder cancer in the New England Bladder Cancer Study. METHODS: Using historical nitrate measurements for public water supplies and measured and modeled values for private wells, as well as self-reported water intake, we estimated average nitrate concentrations (mg/L NO3-N) and average daily nitrate intake (mg/d) from 1970 to diagnosis/reference date (987 cases and 1,180 controls). We estimated overall and source-specific dietary nitrate and nitrite intakes using a food frequency questionnaire (1,037 cases and 1,225 controls). We used unconditional logistic regression to estimate odds ratios (OR) and 95% confidence intervals (CI). We evaluated interactions with factors that may affect N-nitroso compound formation (i.e., red meat, vitamin C, smoking), and with water intake. RESULTS: Average drinking water nitrate concentration above the 95th percentile (>2.07 mg/L) compared with the lowest quartile (≤0.21 mg/L) was associated with bladder cancer (OR = 1.5, 95% CI = 0.97, 2.3; P trend = 0.01); the association was similar for average daily drinking water nitrate intake. We observed positive associations for dietary nitrate and nitrite intakes from processed meat (highest versus lowest quintile OR for nitrate = 1.4, 95% CI = 1.0, 2.0; P trend = 0.04; OR for nitrite = 1.5, 95% CI = 1.0, 2.1; P trend = 0.04, respectively), but not other dietary sources. We observed positive interactions between drinking water nitrate and red meat (P-interaction 0.05) and processed red meat (0.07). CONCLUSIONS: Our results suggest the importance of both drinking water and dietary nitrate sources as risk factors for bladder cancer.

Validity of retrospective occupational exposure estimates of lead and manganese in a case-control study.

OBJECTIVES: The validity of surrogate measures of retrospective occupational exposure in population-based epidemiological studies has rarely been evaluated. Using toenail samples as bioindicators of exposure, we assessed whether work tasks and expert assessments of occupational metal exposure obtained from personal interviews were associated with lead and manganese concentrations. METHODS: We selected 609 controls from a case-control study of bladder cancer in New England who had held a job for ≥1 year 8-24 months prior to toenail collection. We evaluated associations between toenail metal concentrations and five tasks extracted from occupational questionnaires (grinding, painting, soldering, welding, working near engines) using linear regression models. For 139 subjects, we also evaluated associations between the toenail concentrations and exposure estimates from three experts. RESULTS: We observed a 1.9-fold increase (95% CI 1.4 to 2.5) in toenail lead concentrations with painting and 1.4-fold increase (95% CI 1.1 to 1.7) in manganese concentrations with working around engines and handling fuel. We observed significant trends with increasing frequency of both activities. For lead, significant trends were observed with the ratings from all three experts. Their average ratings showed the strongest association, with subjects rated as possibly or probably exposed to lead having concentrations that were 2.0 and 2.5 times higher, respectively, than in unexposed subjects (ptrend <0.001). Expert estimates were only weakly associated with manganese toenail concentrations. CONCLUSIONS: Our findings support the ability of experts to identify broad contrasts in previous occupational exposure to lead. The stronger associations with task frequency and expert assessments support using refined exposure characterisation whenever possible.

Potential effect modifiers of the arsenic-bladder cancer risk relationship.

Populations exposed to arsenic in drinking water have an increased bladder cancer risk and evidence suggests that several factors may modify arsenic metabolism, influencing disease risk. We evaluated whether the association between cumulative lifetime arsenic exposure from drinking water and bladder cancer risk was modified by factors that may impact arsenic metabolism in a population-based case-control study of 1,213 cases and 1,418 controls. Unconditional logistic regression was used to calculate odds ratios (ORs) and 95% confidence intervals (CIs) for the association between cumulative arsenic intake and bladder cancer stratified by age, sex, smoking status, body mass index (BMI), alcohol consumption and folate intake. P-values for interaction were computed using a likelihood ratio test. We observed no statistically significant multiplicative interactions although some variations in associations were notable across risk factors, particularly for smoking and BMI. Among former smokers and current smokers, those with the highest cumulative arsenic intake had elevated risks of bladder cancer (OR = 1.4, 95% CI: 0.96-2.0 and OR = 1.6, 95% CI: 0.91-3.0, respectively; while the OR among never smokers was 1.1, 95% CI: 0.6-1.9, p-interaction = 0.49). Among those classified as normal or overweight based on usual adult BMI, the highest level of cumulative arsenic intake was associated with elevated risks of bladder cancer (OR = 1.3, 95% CI: 0.89-2.0 and OR = 1.6, 95% CI: 1.1-2.4, respectively), while risk was not elevated among those who were obese (OR = 0.9, 95% CI: 0.4-1.8) (p-interaction = 0.14). Our study provides some limited evidence of modifying roles of age, sex, smoking, BMI, folate and alcohol on arsenic-related bladder cancer risk that requires confirmation in other, larger studies.

Identification of a novel susceptibility locus at 13q34 and refinement of the 20p12.2 region as a multi-signal locus associated with bladder cancer risk in individuals of European ancestry.

Candidate gene and genome-wide association studies (GWAS) have identified 15 independent genomic regions associated with bladder cancer risk. In search for additional susceptibility variants, we followed up on four promising single-nucleotide polymorphisms (SNPs) that had not achieved genome-wide significance in 6911 cases and 11 814 controls (rs6104690, rs4510656, rs5003154 and rs4907479, P < 1 × 10(-6)), using additional data from existing GWAS datasets and targeted genotyping for studies that did not have GWAS data. In a combined analysis, which included data on up to 15 058 cases and 286 270 controls, two SNPs achieved genome-wide statistical significance: rs6104690 in a gene desert at 20p12.2 (P = 2.19 × 10(-11)) and rs4907479 within the MCF2L gene at 13q34 (P = 3.3 × 10(-10)). Imputation and fine-mapping analyses were performed in these two regions for a subset of 5551 bladder cancer cases and 10 242 controls. Analyses at the 13q34 region suggest a single signal marked by rs4907479. In contrast, we detected two signals in the 20p12.2 region-the first signal is marked by rs6104690, and the second signal is marked by two moderately correlated SNPs (r(2) = 0.53), rs6108803 and the previously reported rs62185668. The second 20p12.2 signal is more strongly associated with the risk of muscle-invasive (T2-T4 stage) compared with non-muscle-invasive (Ta, T1 stage) bladder cancer (case-case P ≤ 0.02 for both rs62185668 and rs6108803). Functional analyses are needed to explore the biological mechanisms underlying these novel genetic associations with risk for bladder cancer.

Pooled study of occupational exposure to aromatic hydrocarbon solvents and risk of multiple myeloma.

OBJECTIVES: To investigate the association between occupational exposure to aromatic hydrocarbon solvents and risk of multiple myeloma (MM) in a large, consortium-based study. METHODS: We pooled data on 2854 cases and 10 743 controls from nine studies participating in the InterLymph consortium. Occupational exposures to benzene, toluene and xylene were assigned by a job-exposure matrix, coupled with 'correction' of exposure probability by self-reported or expert-assessed exposure from the individual studies. Cumulative intensity was calculated as the job-specific exposure intensity multiplied by job duration, summed across jobs. Associations were estimated using logistic regression, with inclusion of covariates for study matching factors and other potential confounders. We repeated our main analysis using random-effects meta-analysis to evaluate heterogeneity of effect. RESULTS: Benzene, toluene and xylene were each associated with MM. For the three solvents, the highest quartile of high-probability cumulative intensity exposure (vs unexposed) was associated with 42% to 63% increased risks of MM. Associations with toluene and xylene exposures were fairly consistent and robust to sensitivity analyses. The estimated effect for benzene was moderately heterogeneous between the studies. Each solvent's association with MM was stronger for exposure occurring within 20 years of diagnosis than with exposure lagged by more than 20 years. CONCLUSIONS: Our study adds important evidence for a role of aromatic hydrocarbon solvents in causation of MM. The difficulty in disentangling individual compounds in this group and a lack of data on potential carcinogenicity of toluene and xylene, in widespread current use, underscore a need for further epidemiological evaluation.

Genome-wide association study identifies multiple loci associated with bladder cancer risk.

Candidate gene and genome-wide association studies (GWAS) have identified 11 independent susceptibility loci associated with bladder cancer risk. To discover additional risk variants, we conducted a new GWAS of 2422 bladder cancer cases and 5751 controls, followed by a meta-analysis with two independently published bladder cancer GWAS, resulting in a combined analysis of 6911 cases and 11 814 controls of European descent. TaqMan genotyping of 13 promising single nucleotide polymorphisms with P < 1 × 10(-5) was pursued in a follow-up set of 801 cases and 1307 controls. Two new loci achieved genome-wide statistical significance: rs10936599 on 3q26.2 (P = 4.53 × 10(-9)) and rs907611 on 11p15.5 (P = 4.11 × 10(-8)). Two notable loci were also identified that approached genome-wide statistical significance: rs6104690 on 20p12.2 (P = 7.13 × 10(-7)) and rs4510656 on 6p22.3 (P = 6.98 × 10(-7)); these require further studies for confirmation. In conclusion, our study has identified new susceptibility alleles for bladder cancer risk that require fine-mapping and laboratory investigation, which could further understanding into the biological underpinnings of bladder carcinogenesis.

Imputation and subset-based association analysis across different cancer types identifies multiple independent risk loci in the TERT-CLPTM1L region on chromosome 5p15.33.

Genome-wide association studies (GWAS) have mapped risk alleles for at least 10 distinct cancers to a small region of 63 000 bp on chromosome 5p15.33. This region harbors the TERT and CLPTM1L genes; the former encodes the catalytic subunit of telomerase reverse transcriptase and the latter may play a role in apoptosis. To investigate further the genetic architecture of common susceptibility alleles in this region, we conducted an agnostic subset-based meta-analysis (association analysis based on subsets) across six distinct cancers in 34 248 cases and 45 036 controls. Based on sequential conditional analysis, we identified as many as six independent risk loci marked by common single-nucleotide polymorphisms: five in the TERT gene (Region 1: rs7726159, P = 2.10 × 10(-39); Region 3: rs2853677, P = 3.30 × 10(-36) and PConditional = 2.36 × 10(-8); Region 4: rs2736098, P = 3.87 × 10(-12) and PConditional = 5.19 × 10(-6), Region 5: rs13172201, P = 0.041 and PConditional = 2.04 × 10(-6); and Region 6: rs10069690, P = 7.49 × 10(-15) and PConditional = 5.35 × 10(-7)) and one in the neighboring CLPTM1L gene (Region 2: rs451360; P = 1.90 × 10(-18) and PConditional = 7.06 × 10(-16)). Between three and five cancers mapped to each independent locus with both risk-enhancing and protective effects. Allele-specific effects on DNA methylation were seen for a subset of risk loci, indicating that methylation and subsequent effects on gene expression may contribute to the biology of risk variants on 5p15.33. Our results provide strong support for extensive pleiotropy across this region of 5p15.33, to an extent not previously observed in other cancer susceptibility loci.

Computer-based coding of free-text job descriptions to efficiently identify occupations in epidemiological studies.

BACKGROUND: Mapping job titles to standardised occupation classification (SOC) codes is an important step in identifying occupational risk factors in epidemiological studies. Because manual coding is time-consuming and has moderate reliability, we developed an algorithm called SOCcer (Standardized Occupation Coding for Computer-assisted Epidemiologic Research) to assign SOC-2010 codes based on free-text job description components. METHODS: Job title and task-based classifiers were developed by comparing job descriptions to multiple sources linking job and task descriptions to SOC codes. An industry-based classifier was developed based on the SOC prevalence within an industry. These classifiers were used in a logistic model trained using 14 983 jobs with expert-assigned SOC codes to obtain empirical weights for an algorithm that scored each SOC/job description. We assigned the highest scoring SOC code to each job. SOCcer was validated in 2 occupational data sources by comparing SOC codes obtained from SOCcer to expert assigned SOC codes and lead exposure estimates obtained by linking SOC codes to a job-exposure matrix. RESULTS: For 11 991 case-control study jobs, SOCcer-assigned codes agreed with 44.5% and 76.3% of manually assigned codes at the 6-digit and 2-digit level, respectively. Agreement increased with the score, providing a mechanism to identify assignments needing review. Good agreement was observed between lead estimates based on SOCcer and manual SOC assignments (κ 0.6-0.8). Poorer performance was observed for inspection job descriptions, which included abbreviations and worksite-specific terminology. CONCLUSIONS: Although some manual coding will remain necessary, using SOCcer may improve the efficiency of incorporating occupation into large-scale epidemiological studies.

Combining Decision Rules from Classification Tree Models and Expert Assessment to Estimate Occupational Exposure to Diesel Exhaust for a Case-Control Study.

OBJECTIVES: To efficiently and reproducibly assess occupational diesel exhaust exposure in a Spanish case-control study, we examined the utility of applying decision rules that had been extracted from expert estimates and questionnaire response patterns using classification tree (CT) models from a similar US study. METHODS: First, previously extracted CT decision rules were used to obtain initial ordinal (0-3) estimates of the probability, intensity, and frequency of occupational exposure to diesel exhaust for the 10 182 jobs reported in a Spanish case-control study of bladder cancer. Second, two experts reviewed the CT estimates for 350 jobs randomly selected from strata based on each CT rule's agreement with the expert ratings in the original study [agreement rate, from 0 (no agreement) to 1 (perfect agreement)]. Their agreement with each other and with the CT estimates was calculated using weighted kappa (κ w) and guided our choice of jobs for subsequent expert review. Third, an expert review comprised all jobs with lower confidence (low-to-moderate agreement rates or discordant assignments, n = 931) and a subset of jobs with a moderate to high CT probability rating and with moderately high agreement rates (n = 511). Logistic regression was used to examine the likelihood that an expert provided a different estimate than the CT estimate based on the CT rule agreement rates, the CT ordinal rating, and the availability of a module with diesel-related questions. RESULTS: Agreement between estimates made by two experts and between estimates made by each of the experts and the CT estimates was very high for jobs with estimates that were determined by rules with high CT agreement rates (κ w: 0.81-0.90). For jobs with estimates based on rules with lower agreement rates, moderate agreement was observed between the two experts (κ w: 0.42-0.67) and poor-to-moderate agreement was observed between the experts and the CT estimates (κ w: 0.09-0.57). In total, the expert review of 1442 jobs changed 156 probability estimates, 128 intensity estimates, and 614 frequency estimates. The expert was more likely to provide a different estimate when the CT rule agreement rate was <0.8, when the CT ordinal ratings were low to moderate, or when a module with diesel questions was available. CONCLUSIONS: Our reliability assessment provided important insight into where to prioritize additional expert review; as a result, only 14% of the jobs underwent expert review, substantially reducing the exposure assessment burden. Overall, we found that we could efficiently, reproducibly, and reliably apply CT decision rules from one study to assess exposure in another study.

Is the inverse association between selenium and bladder cancer due to confounding by smoking?

Selenium has been linked to a reduced risk of bladder cancer in some studies. Smoking, a well-established risk factor for bladder cancer, has been associated with lower selenium levels in the body. We investigated the selenium-bladder cancer association in subjects from Maine, New Hampshire, and Vermont in the New England Bladder Cancer Case-Control Study. At interview (2001-2005), participants provided information on a variety of factors, including a comprehensive smoking history, and submitted toenail samples, from which we measured selenium levels. We estimated odds ratios and 95% confidence intervals among 1,058 cases and 1,271 controls using logistic regression. After controlling for smoking, we saw no evidence of an association between selenium levels and bladder cancer (for fourth quartile vs. first quartile, odds ratio (OR) = 0.98, 95% confidence interval (CI): 0.77, 1.25). When results were restricted to regular smokers, there appeared to be an inverse association (OR = 0.76, 95% CI: 0.58, 0.99); however, when pack-years of smoking were considered, this association was attenuated (OR = 0.91, 95% CI: 0.68, 1.20), indicating potential confounding by smoking. Despite some reports of an inverse association between selenium and bladder cancer overall, our results, combined with an in-depth evaluation of other studies, suggested that confounding from smoking intensity or duration could explain this association. Our study highlights the need to carefully evaluate the confounding association of smoking in the selenium-bladder cancer association.

Exposure-response relationships for select cancer and non-cancer health outcomes in a cohort of U.S. firefighters from San Francisco, Chicago and Philadelphia (1950-2009).

OBJECTIVES: To examine exposure-response relationships between surrogates of firefighting exposure and select outcomes among previously studied US career firefighters. METHODS: Eight cancer and four non-cancer outcomes were examined using conditional logistic regression. Incidence density sampling was used to match each case to 200 controls on attained age. Days accrued in firefighting assignments (exposed-days), run totals (fire-runs) and run times (fire-hours) were used as exposure surrogates. HRs comparing 75th and 25th centiles of lagged cumulative exposures were calculated using loglinear, linear, log-quadratic, power and restricted cubic spline general relative risk models. Piecewise constant models were used to examine risk differences by time since exposure, age at exposure and calendar period. RESULTS: Among 19,309 male firefighters eligible for the study, there were 1333 cancer deaths and 2609 cancer incidence cases. Significant positive associations between fire-hours and lung cancer mortality and incidence were evident. A similar relation between leukaemia mortality and fire-runs was also found. The lung cancer associations were nearly linear in cumulative exposure, while the association with leukaemia mortality was attenuated at higher exposure levels and greater for recent exposures. Significant negative associations were evident for the exposure surrogates and colorectal and prostate cancers, suggesting a healthy worker survivor effect possibly enhanced by medical screening. CONCLUSIONS: Lung cancer and leukaemia mortality risks were modestly increasing with firefighter exposures. These findings add to evidence of a causal association between firefighting and cancer. Nevertheless, small effects merit cautious interpretation. We plan to continue to follow the occurrence of disease and injury in this cohort.

Comparison of ordinal and nominal classification trees to predict ordinal expert-based occupational exposure estimates in a case-control study.

OBJECTIVES: To evaluate occupational exposures in case-control studies, exposure assessors typically review each job individually to assign exposure estimates. This process lacks transparency and does not provide a mechanism for recreating the decision rules in other studies. In our previous work, nominal (unordered categorical) classification trees (CTs) generally successfully predicted expert-assessed ordinal exposure estimates (i.e. none, low, medium, high) derived from occupational questionnaire responses, but room for improvement remained. Our objective was to determine if using recently developed ordinal CTs would improve the performance of nominal trees in predicting ordinal occupational diesel exhaust exposure estimates in a case-control study. METHODS: We used one nominal and four ordinal CT methods to predict expert-assessed probability, intensity, and frequency estimates of occupational diesel exhaust exposure (each categorized as none, low, medium, or high) derived from questionnaire responses for the 14983 jobs in the New England Bladder Cancer Study. To replicate the common use of a single tree, we applied each method to a single sample of 70% of the jobs, using 15% to test and 15% to validate each method. To characterize variability in performance, we conducted a resampling analysis that repeated the sample draws 100 times. We evaluated agreement between the tree predictions and expert estimates using Somers' d, which measures differences in terms of ordinal association between predicted and observed scores and can be interpreted similarly to a correlation coefficient. RESULTS: From the resampling analysis, compared with the nominal tree, an ordinal CT method that used a quadratic misclassification function and controlled tree size based on total misclassification cost had a slightly better predictive performance that was statistically significant for the frequency metric (Somers' d: nominal tree = 0.61; ordinal tree = 0.63) and similar performance for the probability (nominal = 0.65; ordinal = 0.66) and intensity (nominal = 0.65; ordinal = 0.65) metrics. The best ordinal CT predicted fewer cases of large disagreement with the expert assessments (i.e. no exposure predicted for a job with high exposure and vice versa) compared with the nominal tree across all of the exposure metrics. For example, the percent of jobs with expert-assigned high intensity of exposure that the model predicted as no exposure was 29% for the nominal tree and 22% for the best ordinal tree. CONCLUSIONS: The overall agreements were similar across CT models; however, the use of ordinal models reduced the magnitude of the discrepancy when disagreements occurred. As the best performing model can vary by situation, researchers should consider evaluating multiple CT methods to maximize the predictive performance within their data.

Using hierarchical cluster models to systematically identify groups of jobs with similar occupational questionnaire response patterns to assist rule-based expert exposure assessment in population-based studies.

OBJECTIVES: Rule-based expert exposure assessment based on questionnaire response patterns in population-based studies improves the transparency of the decisions. The number of unique response patterns, however, can be nearly equal to the number of jobs. An expert may reduce the number of patterns that need assessment using expert opinion, but each expert may identify different patterns of responses that identify an exposure scenario. Here, hierarchical clustering methods are proposed as a systematic data reduction step to reproducibly identify similar questionnaire response patterns prior to obtaining expert estimates. As a proof-of-concept, we used hierarchical clustering methods to identify groups of jobs (clusters) with similar responses to diesel exhaust-related questions and then evaluated whether the jobs within a cluster had similar (previously assessed) estimates of occupational diesel exhaust exposure. METHODS: Using the New England Bladder Cancer Study as a case study, we applied hierarchical cluster models to the diesel-related variables extracted from the occupational history and job- and industry-specific questionnaires (modules). Cluster models were separately developed for two subsets: (i) 5395 jobs with ≥1 variable extracted from the occupational history indicating a potential diesel exposure scenario, but without a module with diesel-related questions; and (ii) 5929 jobs with both occupational history and module responses to diesel-relevant questions. For each subset, we varied the numbers of clusters extracted from the cluster tree developed for each model from 100 to 1000 groups of jobs. Using previously made estimates of the probability (ordinal), intensity (µg m(-3) respirable elemental carbon), and frequency (hours per week) of occupational exposure to diesel exhaust, we examined the similarity of the exposure estimates for jobs within the same cluster in two ways. First, the clusters' homogeneity (defined as >75% with the same estimate) was examined compared to a dichotomized probability estimate (<5 versus ≥5%; <50 versus ≥50%). Second, for the ordinal probability metric and continuous intensity and frequency metrics, we calculated the intraclass correlation coefficients (ICCs) between each job's estimate and the mean estimate for all jobs within the cluster. RESULTS: Within-cluster homogeneity increased when more clusters were used. For example, ≥80% of the clusters were homogeneous when 500 clusters were used. Similarly, ICCs were generally above 0.7 when ≥200 clusters were used, indicating minimal within-cluster variability. The most within-cluster variability was observed for the frequency metric (ICCs from 0.4 to 0.8). We estimated that using an expert to assign exposure at the cluster-level assignment and then to review each job in non-homogeneous clusters would require ~2000 decisions per expert, in contrast to evaluating 4255 unique questionnaire patterns or 14983 individual jobs. CONCLUSIONS: This proof-of-concept shows that using cluster models as a data reduction step to identify jobs with similar response patterns prior to obtaining expert ratings has the potential to aid rule-based assessment by systematically reducing the number of exposure decisions needed. While promising, additional research is needed to quantify the actual reduction in exposure decisions and the resulting homogeneity of exposure estimates within clusters for an exposure assessment effort that obtains cluster-level expert assessments as part of the assessment process.

Common genetic variants in the PSCA gene influence gene expression and bladder cancer risk.

Genome-wide association studies have identified a SNP, rs2294008, on 8q24.3 within the prostate stem cell antigen (PSCA) gene, as a risk factor for bladder cancer. To fine-map this region, we imputed 642 SNPs within 100 Kb of rs2294008 in addition to 33 markers genotyped in one of the reported genome-wide association study in 8,652 subjects. A multivariable logistic regression model adjusted for rs2294008 revealed a unique signal, rs2978974 (r(2) = 0.02, D' = 0.19 with rs2294008). In the combined analysis of 5,393 cases and 7,324 controls, we detected a per-allele odds ratio (OR) = 1.11 [95% confidence interval (CI) = 1.06-1.17, P = 5.8 × 10(-5)] for rs2294008 and OR = 1.07 (95% CI = 1.02-1.13, P = 9.7 × 10(-3)) for rs2978974. The effect was stronger in carriers of both risk variants (OR = 1.24, 95% CI = 1.08-1.41, P = 1.8 × 10(-3)) and there was a significant multiplicative interaction (P = 0.035) between these two SNPs, which requires replication in future studies. The T risk allele of rs2294008 was associated with increased PSCA mRNA expression in two sets of bladder tumor samples (n = 36, P = 0.0007 and n = 34, P = 0.0054) and in normal bladder samples (n = 35, P = 0.0155), but rs2978974 was not associated with PSCA expression. SNP rs2978974 is located 10 Kb upstream of rs2294008, within an alternative untranslated first exon of PSCA. The non-risk allele G of rs2978974 showed strong interaction with nuclear proteins from five cell lines tested, implying a regulatory function. In conclusion, a joint effect of two PSCA SNPs, rs2294008 and rs2978974, suggests that both variants may be important for bladder cancer susceptibility, possibly through different mechanisms that influence the control of mRNA expression and interaction with regulatory factors.

Genome-wide interaction study of smoking and bladder cancer risk.

Bladder cancer is a complex disease with known environmental and genetic risk factors. We performed a genome-wide interaction study (GWAS) of smoking and bladder cancer risk based on primary scan data from 3002 cases and 4411 controls from the National Cancer Institute Bladder Cancer GWAS. Alternative methods were used to evaluate both additive and multiplicative interactions between individual single nucleotide polymorphisms (SNPs) and smoking exposure. SNPs with interaction P values < 5 × 10(-) (5) were evaluated further in an independent dataset of 2422 bladder cancer cases and 5751 controls. We identified 10 SNPs that showed association in a consistent manner with the initial dataset and in the combined dataset, providing evidence of interaction with tobacco use. Further, two of these novel SNPs showed strong evidence of association with bladder cancer in tobacco use subgroups that approached genome-wide significance. Specifically, rs1711973 (FOXF2) on 6p25.3 was a susceptibility SNP for never smokers [combined odds ratio (OR) = 1.34, 95% confidence interval (CI) = 1.20-1.50, P value = 5.18 × 10(-) (7)]; and rs12216499 (RSPH3-TAGAP-EZR) on 6q25.3 was a susceptibility SNP for ever smokers (combined OR = 0.75, 95% CI = 0.67-0.84, P value = 6.35 × 10(-) (7)). In our analysis of smoking and bladder cancer, the tests for multiplicative interaction seemed to more commonly identify susceptibility loci with associations in never smokers, whereas the additive interaction analysis identified more loci with associations among smokers-including the known smoking and NAT2 acetylation interaction. Our findings provide additional evidence of gene-environment interactions for tobacco and bladder cancer.

Mapping of the UGT1A locus identifies an uncommon coding variant that affects mRNA expression and protects from bladder cancer.

A recent genome-wide association study of bladder cancer identified the UGT1A gene cluster on chromosome 2q37.1 as a novel susceptibility locus. The UGT1A cluster encodes a family of UDP-glucuronosyltransferases (UGTs), which facilitate cellular detoxification and removal of aromatic amines. Bioactivated forms of aromatic amines found in tobacco smoke and industrial chemicals are the main risk factors for bladder cancer. The association within the UGT1A locus was detected by a single nucleotide polymorphism (SNP) rs11892031. Now, we performed detailed resequencing, imputation and genotyping in this region. We clarified the original genetic association detected by rs11892031 and identified an uncommon SNP rs17863783 that explained and strengthened the association in this region (allele frequency 0.014 in 4035 cases and 0.025 in 5284 controls, OR = 0.55, 95%CI = 0.44-0.69, P = 3.3 × 10(-7)). Rs17863783 is a synonymous coding variant Val209Val within the functional UGT1A6.1 splicing form, strongly expressed in the liver, kidney and bladder. We found the protective T allele of rs17863783 to be associated with increased mRNA expression of UGT1A6.1 in in-vitro exontrap assays and in human liver tissue samples. We suggest that rs17863783 may protect from bladder cancer by increasing the removal of carcinogens from bladder epithelium by the UGT1A6.1 protein. Our study shows an example of genetic and functional role of an uncommon protective genetic variant in a complex human disease, such as bladder cancer.

A prospective study of circulating adipokine levels and risk of multiple myeloma.

It has been hypothesized that the observed excess risk of multiple myeloma (MM) among obese persons could be the result of altered circulating levels of adipokines, polypeptide hormones with pro- and anti-inflammatory properties secreted by adipose tissue. We investigated whether circulating levels of leptin, total adiponectin, and high molecular weight adiponectin are associated with subsequent MM risk among 174 MM patients and 348 controls within the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. Inverse associations with MM were observed for total adiponectin (highest quartile vs lowest: odds ratio = 0.49; 95% CI = 0.26-0.93, P(trend) = .03) and high molecular weight adiponectin (0.44; 0.23-0.85, P(trend) = .01). These associations remained after restricting to MM patients diagnosed ∼ 8 years or more after blood collection. Leptin levels were not associated with MM risk. The results of this study, to our knowledge the first prospective investigation of circulating adipokines and MM, suggest that adiponectin may play an important role in obesity-related myelomagenesis.