In environmental epidemiology, measurements of exposure biomarkers often fall below the assay's limit of detection. Existing methods for handling this problem, including deletion, substitution, parametric regression, and multiple imputation, can perform poorly if the proportion of "nondetects" is high or parametric models are mis-specified. We propose an approach that treats the measured analyte as the modeled outcome, implying a role reversal when the analyte is a putative cause of a health outcome. Following a scale reversal as well, our approach uses Cox regression to model the analyte, with confounder adjustment. The method makes full use of quantifiable analyte measures, while appropriately treating nondetects as censored. Under the proportional hazards assumption, the hazard ratio for a binary health outcome is interpretable as an adjusted odds ratio: the odds for the outcome at any particular analyte concentration divided by the odds given a lower concentration. Our approach is broadly applicable to cohort studies, case-control studies (frequency matched or not), and cross-sectional studies conducted to identify determinants of exposure. We illustrate the method with cross-sectional survey data to assess sex as a determinant of 2,3,7,8-tetrachlorodibenzo-p-dioxin concentration and with prospective cohort data to assess the association between 2,4,4'-trichlorobiphenyl exposure and psychomotor development.

Background

Phthalate exposures are ubiquitous during pregnancy and may contribute to racial and ethnic disparities in preterm birth.

Objectives

We investigated race and ethnicity in the relationship between biomarkers of phthalate exposure and preterm birth by examining: a) how hypothetical reductions in racial and ethnic disparities in phthalate metabolites might reduce the probability of preterm birth; and b) exposure-response models stratified by race and ethnicity.

Methods

We pooled individual-level data on 6,045 pregnancies from 16 U.S. cohorts. We investigated covariate-adjusted differences in nine urinary phthalate metabolite concentrations by race and ethnicity [non-Hispanic White (White, 43%), non-Hispanic Black (Black, 13%), Hispanic/Latina (38%), and Asian/Pacific Islander (3%)]. Using g-computation, we estimated changes in the probability of preterm birth under hypothetical interventions to eliminate disparities in levels of urinary phthalate metabolites by proportionally lowering average concentrations in Black and Hispanic/Latina participants to be approximately equal to the averages in White participants. We also used race and ethnicity-stratified logistic regression to characterize associations between phthalate metabolites and preterm birth.

Results

In comparison with concentrations among White participants, adjusted mean phthalate metabolite concentrations were consistently higher among Black and Hispanic/Latina participants by 23%-148% and 4%-94%, respectively. Asian/Pacific Islander participants had metabolite levels that were similar to those of White participants. Hypothetical interventions to reduce disparities in metabolite mixtures were associated with lower probabilities of preterm birth for Black [13% relative reduction; 95% confidence interval (CI): -34%, 8.6%] and Hispanic/Latina (9% relative reduction; 95% CI: -19%, 0.8%) participants. Odds ratios for preterm birth in association with phthalate metabolites demonstrated heterogeneity by race and ethnicity for two individual metabolites (mono-n-butyl and monoisobutyl phthalate), with positive associations that were larger in magnitude observed among Black or Hispanic/Latina participants.

Conclusions

Phthalate metabolite concentrations differed substantially by race and ethnicity. Our results show hypothetical interventions to reduce population-level racial and ethnic disparities in biomarkers of phthalate exposure could potentially reduce the probability of preterm birth. https://doi.org/10.1289/EHP12831.

Background

Recently, the International Agency for Research on Cancer (IARC) Programme for the Evaluation of Carcinogenic Risks to Humans has been criticized for several of its evaluations, and also for the approach used to perform these evaluations. Some critics have claimed that failures of IARC Working Groups to recognize study weaknesses and biases of Working Group members have led to inappropriate classification of a number of agents as carcinogenic to humans.

Objectives

The authors of this Commentary are scientists from various disciplines relevant to the identification and hazard evaluation of human carcinogens. We examined criticisms of the IARC classification process to determine the validity of these concerns. Here, we present the results of that examination, review the history of IARC evaluations, and describe how the IARC evaluations are performed.

Discussion

We concluded that these recent criticisms are unconvincing. The procedures employed by IARC to assemble Working Groups of scientists from the various disciplines and the techniques followed to review the literature and perform hazard assessment of various agents provide a balanced evaluation and an appropriate indication of the weight of the evidence. Some disagreement by individual scientists to some evaluations is not evidence of process failure. The review process has been modified over time and will undoubtedly be altered in the future to improve the process. Any process can in theory be improved, and we would support continued review and improvement of the IARC processes. This does not mean, however, that the current procedures are flawed.

Conclusions

The IARC Monographs have made, and continue to make, major contributions to the scientific underpinning for societal actions to improve the public's health.

Background

Genome-wide studies of gene-environment interactions (G×E) may identify variants associated with disease risk in conjunction with lifestyle/environmental exposures. We conducted a genome-wide G×E analysis of ~ 7.6 million common variants and seven lifestyle/environmental risk factors for breast cancer risk overall and for estrogen receptor positive (ER +) breast cancer.

Methods

Analyses were conducted using 72,285 breast cancer cases and 80,354 controls of European ancestry from the Breast Cancer Association Consortium. Gene-environment interactions were evaluated using standard unconditional logistic regression models and likelihood ratio tests for breast cancer risk overall and for ER + breast cancer. Bayesian False Discovery Probability was employed to assess the noteworthiness of each SNP-risk factor pairs.

Results

Assuming a 1 × 10^-5 prior probability of a true association for each SNP-risk factor pairs and a Bayesian False Discovery Probability < 15%, we identified two independent SNP-risk factor pairs: rs80018847(9p13)-LINGO2 and adult height in association with overall breast cancer risk (OR_int = 0.94, 95% CI 0.92-0.96), and rs4770552(13q12)-SPATA13 and age at menarche for ER + breast cancer risk (OR_int = 0.91, 95% CI 0.88-0.94).

Conclusions

Overall, the contribution of G×E interactions to the heritability of breast cancer is very small. At the population level, multiplicative G×E interactions do not make an important contribution to risk prediction in breast cancer.

A large number of variants identified through clinical genetic testing in disease susceptibility genes, are of uncertain significance (VUS). Following the recommendations of the American College of Medical Genetics and Genomics (ACMG) and Association for Molecular Pathology (AMP), the frequency in case-control datasets (PS4 criterion), can inform their interpretation. We present a novel case-control likelihood ratio-based method that incorporates gene-specific age-related penetrance. We demonstrate the utility of this method in the analysis of simulated and real datasets. In the analyses of simulated data, the likelihood ratio method was more powerful compared to other methods. Likelihood ratios were calculated for a case-control dataset of BRCA1 and BRCA2 variants from the Breast Cancer Association Consortium (BCAC), and compared with logistic regression results. A larger number of variants reached evidence in favor of pathogenicity, and a substantial number of variants had evidence against pathogenicity - findings that would not have been reached using other case-control analysis methods. Our novel method provides greater power to classify rare variants compared to classical case-control methods. As an initiative from the ENIGMA Analytical Working Group, we provide user-friendly scripts and pre-formatted excel calculators for implementation of the method for rare variants in BRCA1, BRCA2 and other high-risk genes with known penetrance.

Background

Low-frequency variants play an important role in breast cancer (BC) susceptibility. Gene-based methods can increase power by combining multiple variants in the same gene and help identify target genes.

Methods

We evaluated the potential of gene-based aggregation in the Breast Cancer Association Consortium cohorts including 83,471 cases and 59,199 controls. Low-frequency variants were aggregated for individual genes' coding and regulatory regions. Association results in European ancestry samples were compared to single-marker association results in the same cohort. Gene-based associations were also combined in meta-analysis across individuals with European, Asian, African, and Latin American and Hispanic ancestry.

Results

In European ancestry samples, 14 genes were significantly associated (q < 0.05) with BC. Of those, two genes, FMNL3 (P = 6.11 × 10^-6) and AC058822.1 (P = 1.47 × 10^-4), represent new associations. High FMNL3 expression has previously been linked to poor prognosis in several other cancers. Meta-analysis of samples with diverse ancestry discovered further associations including established candidate genes ESR1 and CBLB. Furthermore, literature review and database query found further support for a biologically plausible link with cancer for genes CBLB, FMNL3, FGFR2, LSP1, MAP3K1, and SRGAP2C.

Conclusions

Using extended gene-based aggregation tests including coding and regulatory variation, we report identification of plausible target genes for previously identified single-marker associations with BC as well as the discovery of novel genes implicated in BC development. Including multi ancestral cohorts in this study enabled the identification of otherwise missed disease associations as ESR1 (P = 1.31 × 10^-5), demonstrating the importance of diversifying study cohorts.

Background

Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear.

Methods

Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes.

Results

Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions.

Conclusion

This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.

Objectives

Physical inactivity and sedentary behaviour are associated with higher breast cancer risk in observational studies, but ascribing causality is difficult. Mendelian randomisation (MR) assesses causality by simulating randomised trial groups using genotype. We assessed whether lifelong physical activity or sedentary time, assessed using genotype, may be causally associated with breast cancer risk overall, pre/post-menopause, and by case-groups defined by tumour characteristics.

Methods

We performed two-sample inverse-variance-weighted MR using individual-level Breast Cancer Association Consortium case-control data from 130 957 European-ancestry women (69 838 invasive cases), and published UK Biobank data (n=91 105-377 234). Genetic instruments were single nucleotide polymorphisms (SNPs) associated in UK Biobank with wrist-worn accelerometer-measured overall physical activity (n_snps=5) or sedentary time (n_snps=6), or accelerometer-measured (n_snps=1) or self-reported (n_snps=5) vigorous physical activity.

Results

Greater genetically-predicted overall activity was associated with lower breast cancer overall risk (OR=0.59; 95% confidence interval (CI) 0.42 to 0.83 per-standard deviation (SD;~8 milligravities acceleration)) and for most case-groups. Genetically-predicted vigorous activity was associated with lower risk of pre/perimenopausal breast cancer (OR=0.62; 95% CI 0.45 to 0.87,≥3 vs. 0 self-reported days/week), with consistent estimates for most case-groups. Greater genetically-predicted sedentary time was associated with higher hormone-receptor-negative tumour risk (OR=1.77; 95% CI 1.07 to 2.92 per-SD (~7% time spent sedentary)), with elevated estimates for most case-groups. Results were robust to sensitivity analyses examining pleiotropy (including weighted-median-MR, MR-Egger).

Conclusion

Our study provides strong evidence that greater overall physical activity, greater vigorous activity, and lower sedentary time are likely to reduce breast cancer risk. More widespread adoption of active lifestyles may reduce the burden from the most common cancer in women.

A combination of genetic and functional approaches has identified three independent breast cancer risk loci at 2q35. A recent fine-scale mapping analysis to refine these associations resulted in 1 (signal 1), 5 (signal 2), and 42 (signal 3) credible causal variants at these loci. We used publicly available in silico DNase I and ChIP-seq data with in vitro reporter gene and CRISPR assays to annotate signals 2 and 3. We identified putative regulatory elements that enhanced cell-type-specific transcription from the IGFBP5 promoter at both signals (30- to 40-fold increased expression by the putative regulatory element at signal 2, 2- to 3-fold by the putative regulatory element at signal 3). We further identified one of the five credible causal variants at signal 2, a 1.4 kb deletion (esv3594306), as the likely causal variant; the deletion allele of this variant was associated with an average additional increase in IGFBP5 expression of 1.3-fold (MCF-7) and 2.2-fold (T-47D). We propose a model in which the deletion allele of esv3594306 juxtaposes two transcription factor binding regions (annotated by estrogen receptor alpha ChIP-seq peaks) to generate a single extended regulatory element. This regulatory element increases cell-type-specific expression of the tumor suppressor gene IGFBP5 and, thereby, reduces risk of estrogen receptor-positive breast cancer (odds ratio = 0.77, 95% CI 0.74-0.81, p = 3.1 × 10^-31).

Background

Epidemiological studies provide strong evidence for a role of endogenous sex hormones in the aetiology of breast cancer. The aim of this analysis was to identify genetic variants that are associated with urinary sex-hormone levels and breast cancer risk.

Methods

We carried out a genome-wide association study of urinary oestrone-3-glucuronide and pregnanediol-3-glucuronide levels in 560 premenopausal women, with additional analysis of progesterone levels in 298 premenopausal women. To test for the association with breast cancer risk, we carried out follow-up genotyping in 90,916 cases and 89,893 controls from the Breast Cancer Association Consortium. All women were of European ancestry.

Results

For pregnanediol-3-glucuronide, there were no genome-wide significant associations; for oestrone-3-glucuronide, we identified a single peak mapping to the CYP3A locus, annotated by rs45446698. The minor rs45446698-C allele was associated with lower oestrone-3-glucuronide (-49.2%, 95% CI -56.1% to -41.1%, P = 3.1 × 10^-18); in follow-up analyses, rs45446698-C was also associated with lower progesterone (-26.7%, 95% CI -39.4% to -11.6%, P = 0.001) and reduced risk of oestrogen and progesterone receptor-positive breast cancer (OR = 0.86, 95% CI 0.82-0.91, P = 6.9 × 10^-8).

Conclusions

The CYP3A7*1C allele is associated with reduced risk of hormone receptor-positive breast cancer possibly mediated via an effect on the metabolism of endogenous sex hormones in premenopausal women.