High-throughput PRIME-editing screens identify functional DNA variants in the human genome.

Despite tremendous progress in detecting DNA variants associated with human disease, interpreting their functional impact in a high-throughput and single-base resolution manner remains challenging. Here, we develop a pooled prime-editing screen method, PRIME, that can be applied to characterize thousands of coding and non-coding variants in a single experiment with high reproducibility. To showcase its applications, we first identified essential nucleotides for a 716 bp MYC enhancer via PRIME-mediated single-base resolution analysis. Next, we applied PRIME to functionally characterize 1,304 genome-wide association study (GWAS)-identified non-coding variants associated with breast cancer and 3,699 variants from ClinVar. We discovered that 103 non-coding variants and 156 variants of uncertain significance are functional via affecting cell fitness. Collectively, we demonstrate that PRIME is capable of characterizing genetic variants at single-base resolution and scale, advancing accurate genome annotation for disease risk prediction, diagnosis, and therapeutic target identification.

Germline genetic contribution to the immune landscape of cancer.

Understanding the contribution of the host's genetic background to cancer immunity may lead to improved stratification for immunotherapy and to the identification of novel therapeutic targets. We investigated the effect of common and rare germline variants on 139 well-defined immune traits in ∼9000 cancer patients enrolled in TCGA. High heritability was observed for estimates of NK cell and T cell subset infiltration and for interferon signaling. Common variants of IFIH1, TMEM173 (STING1), and TMEM108 were associated with differential interferon signaling and variants mapping to RBL1 correlated with T cell subset abundance. Pathogenic or likely pathogenic variants in BRCA1 and in genes involved in telomere stabilization and Wnt-ß-catenin also acted as immune modulators. Our findings provide evidence for the impact of germline genetics on the composition and functional orientation of the tumor immune microenvironment. The curated datasets, variants, and genes identified provide a resource toward further understanding of tumor-immune interactions.

A comprehensive analysis of clinical and polygenic germline influences on somatic mutational burden.

Tumor mutational burden (TMB), the total number of somatic mutations in the tumor, and copy number burden (CNB), the corresponding measure of aneuploidy, are established fundamental somatic features and emerging biomarkers for immunotherapy. However, the genetic and non-genetic influences on TMB/CNB and, critically, the manner by which they influence patient outcomes remain poorly understood. Here, we present a large germline-somatic study of TMB/CNB with >23,000 individuals across 17 cancer types, of which 12,000 also have extensive clinical, treatment, and overall survival (OS) measurements available. We report dozens of clinical associations with TMB/CNB, observing older age and male sex to have a strong effect on TMB and weaker impact on CNB. We additionally identified significant germline influences on TMB/CNB, including fine-scale European ancestry and germline polygenic risk scores (PRSs) for smoking, tanning, white blood cell counts, and educational attainment. We quantify the causal effect of exposures on somatic mutational processes using Mendelian randomization. Many of the identified features associated with TMB/CNB were additionally associated with OS for individuals treated at a single tertiary cancer center. For individuals receiving immunotherapy, we observed a complex relationship between PRSs for educational attainment, self-reported college attainment, TMB, and survival, suggesting that the influence of this biomarker may be substantially modified by socioeconomic status. While the accumulation of somatic alterations is a stochastic process, our work demonstrates that it can be shaped by host characteristics including germline genetics.

The Immune Landscape of Cancer.

We performed an extensive immunogenomic analysis of more than 10,000 tumors comprising 33 diverse cancer types by utilizing data compiled by TCGA. Across cancer types, we identified six immune subtypes-wound healing, IFN-γ dominant, inflammatory, lymphocyte depleted, immunologically quiet, and TGF-ß dominant-characterized by differences in macrophage or lymphocyte signatures, Th1:Th2 cell ratio, extent of intratumoral heterogeneity, aneuploidy, extent of neoantigen load, overall cell proliferation, expression of immunomodulatory genes, and prognosis. Specific driver mutations correlated with lower (CTNNB1, NRAS, or IDH1) or higher (BRAF, TP53, or CASP8) leukocyte levels across all cancers. Multiple control modalities of the intracellular and extracellular networks (transcription, microRNAs, copy number, and epigenetic processes) were involved in tumor-immune cell interactions, both across and within immune subtypes. Our immunogenomics pipeline to characterize these heterogeneous tumors and the resulting data are intended to serve as a resource for future targeted studies to further advance the field.

A pleiotropic variant in DNAJB4 is associated with multiple myeloma risk.

Pleiotropy, which consists of a single gene or allelic variant affecting multiple unrelated traits, is common across cancers, with evidence for genome-wide significant loci shared across cancer and noncancer traits. This feature is particularly relevant in multiple myeloma (MM) because several susceptibility loci that have been identified to date are pleiotropic. Therefore, the aim of this study was to identify novel pleiotropic variants involved in MM risk using 28 684 independent single nucleotide polymorphisms (SNPs) from GWAS Catalog that reached a significant association (P < 5 × 10-8 ) with their respective trait. The selected SNPs were analyzed in 2434 MM cases and 3446 controls from the International Lymphoma Epidemiology Consortium (InterLymph). The 10 SNPs showing the strongest associations with MM risk in InterLymph were selected for replication in an independent set of 1955 MM cases and 1549 controls from the International Multiple Myeloma rESEarch (IMMEnSE) consortium and 418 MM cases and 147 282 controls from the FinnGen project. The combined analysis of the three studies identified an association between DNAJB4-rs34517439-A and an increased risk of developing MM (OR = 1.22, 95%CI 1.13-1.32, P = 4.81 × 10-7 ). rs34517439-A is associated with a modified expression of the FUBP1 gene, which encodes a multifunctional DNA and RNA-binding protein that it was observed to influence the regulation of various genes involved in cell cycle regulation, among which various oncogenes and oncosuppressors. In conclusion, with a pleiotropic scan approach we identified DNAJB4-rs34517439 as a potentially novel MM risk locus.

Novel insights into the whole-blood DNA methylome of asthma in ethnically diverse children and youth.

BACKGROUND: The epigenetic mechanisms of asthma remain largely understudied in African Americans and Hispanics/Latinos, two populations disproportionately affected by asthma. We aimed to identify markers, regions and processes with differential patterns of DNA methylation (DNAm) in whole blood by asthma status in ethnically diverse children and youth, and to assess their functional consequences. METHODS: DNAm levels were profiled with the Infinium MethylationEPIC or HumanMethylation450 BeadChip arrays among 1226 African Americans or Hispanics/Latinos and assessed for differential methylation per asthma status at the CpG and region (differentially methylated region (DMR)) level. Novel associations were validated in blood and/or nasal epithelium from ethnically diverse children and youth. The functional and biological implications of the markers identified were investigated by combining epigenomics with transcriptomics from study participants. RESULTS: 128 CpGs and 196 DMRs were differentially methylated after multiple testing corrections, including 92.3% and 92.8% novel associations, respectively. 41 CpGs were replicated in other Hispanics/Latinos, prioritising cg17647904 (NCOR2) and cg16412914 (AXIN1) as asthma DNAm markers. Significant DNAm markers were enriched in previous associations for asthma, fractional exhaled nitric oxide, bacterial infections, immune regulation or eosinophilia. Functional annotation highlighted epigenetically regulated gene networks involved in corticosteroid response, host defence and immune regulation. Several implicated genes are targets for approved or experimental drugs, including TNNC1 and NDUFA12. Many differentially methylated loci previously associated with asthma were validated in our study. CONCLUSIONS: We report novel whole-blood DNAm markers for asthma underlying key processes of the disease pathophysiology and confirm the transferability of previous asthma DNAm associations to ethnically diverse populations.

Promotores' perspectives on the virtual adaptation of a hereditary breast cancer education program.

Breast cancer is the most common cancer in women in the United States (U.S.) and the leading cause of cancer related death among U.S. Hispanics/Latinas (H/Ls). H/Ls have lower rates of screening and longer time to follow up after an abnormal mammogram. We developed a comprehensive community health educator (promotores)-led education and risk identification program for Spanish-speaking H/Ls in California to increase mammography screening, genetic testing, and the understanding of the impact of family history on cancer risk. Due to COVID-19, we adapted the program to a virtual platform. The experience of transforming the program to a virtual platform provided unique opportunities for collaboration between researchers, community partners, and participants. Promotores are major partners in community based participatory research and in the provision of health care services, but their voices are often excluded from scientific reports. This commentary is an effort to provide a platform for promotores to share their experiences and for the readers to understand their approach in bridging the gap between health care services and communities.

Differences in somatic TP53 mutation type in breast tumors by race and receptor status.

PURPOSE: Somatic driver mutations in TP53 are associated with triple-negative breast cancer (TNBC) and poorer outcomes. Breast cancers in women of African ancestry (AA) are more likely to be TNBC and have somatic TP53 mutations than cancers in non-Hispanic White (NHW) women. Missense driver mutations in TP53 have varied functional impact including loss-of-function (LOF) or gain-of-function (GOF) activity, and dominant negative (DNE) effects. We aimed to determine if there were differences in somatic TP53 mutation types by patient ancestry or TNBC status. METHODS: We identified breast cancer datasets with somatic TP53 mutation data, ancestry, age, and hormone receptor status. Mutations were classified for functional impact using published data and type of mutation. We assessed differences using Fisher's exact test. RESULTS: From 96 breast cancer studies, we identified 2964 women with somatic TP53 mutations: 715 (24.1%) Asian, 258 (8.7%) AA, 1931 (65.2%) NHW, and 60 (2%) Latina. The distribution of TP53 mutation type was similar by ancestry. However, 35.8% of tumors from NHW individuals had GOF mutations compared to 29% from AA individuals (p = 0.04). Mutations with DNE activity were positively associated with TNBC (OR 1.37, p = 0.03) and estrogen receptor (ER) negative status (OR 1.38; p = 0.005). CONCLUSIONS: Somatic TP53 mutation types did not differ by ancestry overall, but GOF mutations were more common in NHW women than AA women. ER-negative and TNBC tumors are less likely to have DNE+ TP53 mutations which could reflect biological processes. Larger cohorts and functional studies are needed to further elucidate these findings.

Assessment of genetic susceptibility to multiple primary cancers through whole-exome sequencing in two large multi-ancestry studies.

BACKGROUND: Up to one of every six individuals diagnosed with one cancer will be diagnosed with a second primary cancer in their lifetime. Genetic factors contributing to the development of multiple primary cancers, beyond known cancer syndromes, have been underexplored. METHODS: To characterize genetic susceptibility to multiple cancers, we conducted a pan-cancer, whole-exome sequencing study of individuals drawn from two large multi-ancestry populations (6429 cases, 165,853 controls). We created two groupings of individuals diagnosed with multiple primary cancers: (1) an overall combined set with at least two cancers across any of 36 organ sites and (2) cancer-specific sets defined by an index cancer at one of 16 organ sites with at least 50 cases from each study population. We then investigated whether variants identified from exome sequencing were associated with these sets of multiple cancer cases in comparison to individuals with one and, separately, no cancers. RESULTS: We identified 22 variant-phenotype associations, 10 of which have not been previously discovered and were significantly overrepresented among individuals with multiple cancers, compared to those with a single cancer. CONCLUSIONS: Overall, we describe variants and genes that may play a fundamental role in the development of multiple primary cancers and improve our understanding of shared mechanisms underlying carcinogenesis.

Expression quantitative trait loci of genes predicting outcome are associated with survival of multiple myeloma patients.

Gene expression profiling can be used for predicting survival in multiple myeloma (MM) and identifying patients who will benefit from particular types of therapy. Some germline single nucleotide polymorphisms (SNPs) act as expression quantitative trait loci (eQTLs) showing strong associations with gene expression levels. We performed an association study to test whether eQTLs of genes reported to be associated with prognosis of MM patients are directly associated with measures of adverse outcome. Using the genotype-tissue expression portal, we identified a total of 16 candidate genes with at least one eQTL SNP associated with their expression with P < 10-7 either in EBV-transformed B-lymphocytes or whole blood. We genotyped the resulting 22 SNPs in 1327 MM cases from the International Multiple Myeloma rESEarch (IMMEnSE) consortium and examined their association with overall survival (OS) and progression-free survival (PFS), adjusting for age, sex, country of origin and disease stage. Three polymorphisms in two genes (TBRG4-rs1992292, TBRG4-rs2287535 and ENTPD1-rs2153913) showed associations with OS at P < .05, with the former two also associated with PFS. The associations of two polymorphisms in TBRG4 with OS were replicated in 1277 MM cases from the International Lymphoma Epidemiology (InterLymph) Consortium. A meta-analysis of the data from IMMEnSE and InterLymph (2579 cases) showed that TBRG4-rs1992292 is associated with OS (hazard ratio = 1.14, 95% confidence interval 1.04-1.26, P = .007). In conclusion, we found biologically a plausible association between a SNP in TBRG4 and OS of MM patients.

Cancer health disparities in racial/ethnic minorities in the United States.

There are well-established disparities in cancer incidence and outcomes by race/ethnicity that result from the interplay between structural, socioeconomic, socio-environmental, behavioural and biological factors. However, large research studies designed to investigate factors contributing to cancer aetiology and progression have mainly focused on populations of European origin. The limitations in clinicopathological and genetic data, as well as the reduced availability of biospecimens from diverse populations, contribute to the knowledge gap and have the potential to widen cancer health disparities. In this review, we summarise reported disparities and associated factors in the United States of America (USA) for the most common cancers (breast, prostate, lung and colon), and for a subset of other cancers that highlight the complexity of disparities (gastric, liver, pancreas and leukaemia). We focus on populations commonly identified and referred to as racial/ethnic minorities in the USA-African Americans/Blacks, American Indians and Alaska Natives, Asians, Native Hawaiians/other Pacific Islanders and Hispanics/Latinos. We conclude that even though substantial progress has been made in understanding the factors underlying cancer health disparities, marked inequities persist. Additional efforts are needed to include participants from diverse populations in the research of cancer aetiology, biology and treatment. Furthermore, to eliminate cancer health disparities, it will be necessary to facilitate access to, and utilisation of, health services to all individuals, and to address structural inequities, including racism, that disproportionally affect racial/ethnic minorities in the USA.

Whole-Genome Sequencing Identifies Novel Functional Loci Associated with Lung Function in Puerto Rican Youth.

Rationale: Puerto Ricans have the highest childhood asthma prevalence in the United States (23.6%); however, the etiology is uncertain.Objectives: In this study, we sought to uncover the genetic architecture of lung function in Puerto Rican youth with and without asthma who were recruited from the island (n = 836).Methods: We used admixture-mapping and whole-genome sequencing data to discover genomic regions associated with lung function. Functional roles of the prioritized candidate SNPs were examined with chromatin immunoprecipitation sequencing, RNA sequencing, and expression quantitative trait loci data.Measurements and Main Results: We discovered a genomic region at 1q32 that was significantly associated with a 0.12-L decrease in the lung volume of exhaled air (95% confidence interval, -0.17 to -0.07; P = 6.62 × 10-8) with each allele of African ancestry. Within this region, two SNPs were expression quantitative trait loci of TMEM9 in nasal airway epithelial cells and MROH3P in esophagus mucosa. The minor alleles of these SNPs were associated with significantly decreased lung function and decreased TMEM9 gene expression. Another admixture-mapping peak was observed on chromosome 5q35.1, indicating that each Native American ancestry allele was associated with a 0.15-L increase in lung function (95% confidence interval, 0.08-0.21; P = 5.03 × 10-6). The region-based association tests identified four suggestive windows that harbored candidate rare variants associated with lung function.Conclusions: We identified common and rare genetic variants that may play a critical role in lung function among Puerto Rican youth. We independently validated an inflammatory pathway that could potentially be used to develop more targeted treatments and interventions for patients with asthma.

Cis-eQTL-based trans-ethnic meta-analysis reveals novel genes associated with breast cancer risk.

Breast cancer is the most common solid organ malignancy and the most frequent cause of cancer death among women worldwide. Previous research has yielded insights into its genetic etiology, but there remains a gap in the understanding of genetic factors that contribute to risk, and particularly in the biological mechanisms by which genetic variation modulates risk. The National Cancer Institute's "Up for a Challenge" (U4C) competition provided an opportunity to further elucidate the genetic basis of the disease. Our group leveraged the seven datasets made available by the U4C organizers and data from the publicly available UK Biobank cohort to examine associations between imputed gene expression and breast cancer risk. In particular, we used reference datasets describing the breast tissue and whole blood transcriptomes to impute expression levels in breast cancer cases and controls. In trans-ethnic meta-analyses of U4C and UK Biobank data, we found significant associations between breast cancer risk and the expression of RCCD1 (joint p-value: 3.6x10-06) and DHODH (p-value: 7.1x10-06) in breast tissue, as well as a suggestive association for ANKLE1 (p-value: 9.3x10-05). Expression of RCCD1 in whole blood was also suggestively associated with disease risk (p-value: 1.2x10-05), as were expression of ACAP1 (p-value: 1.9x10-05) and LRRC25 (p-value: 5.2x10-05). While genome-wide association studies (GWAS) have implicated RCCD1 and ANKLE1 in breast cancer risk, they have not identified the remaining three genes. Among the genetic variants that contributed to the predicted expression of the five genes, we found 23 nominally (p-value < 0.05) associated with breast cancer risk, among which 15 are not in high linkage disequilibrium with risk variants previously identified by GWAS. In summary, we used a transcriptome-based approach to investigate the genetic underpinnings of breast carcinogenesis. This approach provided an avenue for deciphering the functional relevance of genes and genetic variants involved in breast cancer.

Assessing the performance of in silico methods for predicting the pathogenicity of variants in the gene CHEK2, among Hispanic females with breast cancer.

The availability of disease-specific genomic data is critical for developing new computational methods that predict the pathogenicity of human variants and advance the field of precision medicine. However, the lack of gold standards to properly train and benchmark such methods is one of the greatest challenges in the field. In response to this challenge, the scientific community is invited to participate in the Critical Assessment for Genome Interpretation (CAGI), where unpublished disease variants are available for classification by in silico methods. As part of the CAGI-5 challenge, we evaluated the performance of 18 submissions and three additional methods in predicting the pathogenicity of single nucleotide variants (SNVs) in checkpoint kinase 2 (CHEK2) for cases of breast cancer in Hispanic females. As part of the assessment, the efficacy of the analysis method and the setup of the challenge were also considered. The results indicated that though the challenge could benefit from additional participant data, the combined generalized linear model analysis and odds of pathogenicity analysis provided a framework to evaluate the methods submitted for SNV pathogenicity identification and for comparison to other available methods. The outcome of this challenge and the approaches used can help guide further advancements in identifying SNV-disease relationships.

Identification of novel common breast cancer risk variants at the 6q25 locus among Latinas.

BACKGROUND: Breast cancer is a partially heritable trait and genome-wide association studies (GWAS) have identified over 180 common genetic variants associated with breast cancer. We have previously performed breast cancer GWAS in Latinas and identified a strongly protective single nucleotide polymorphism (SNP) at 6q25, with the protective minor allele originating from indigenous American ancestry. Here we report on fine mapping of the 6q25 locus in an expanded sample of Latinas. METHODS: We performed GWAS in 2385 cases and 6416 controls who were either US Latinas or Mexican women. We replicated the top SNPs in 2412 cases and 1620 controls of US Latina, Mexican, and Colombian women. In addition, we validated the top novel variants in studies of African, Asian and European ancestry. In each dataset we used logistic regression models to test the association between SNPs and breast cancer risk and corrected for genetic ancestry using either principal components or genetic ancestry inferred from ancestry informative markers using a model-based approach. RESULTS: We identified a novel set of SNPs at the 6q25 locus associated with genome-wide levels of significance (p = 3.3 × 10- 8 - 6.0 × 10- 9) not in linkage disequilibrium (LD) with variants previously reported at this locus. These SNPs were in high LD (r2 > 0.9) with each other, with the top SNP, rs3778609, associated with breast cancer with an odds ratio (OR) and 95% confidence interval (95% CI) of 0.76 (0.70-0.84). In a replication in women of Latin American origin, we also observed a consistent effect (OR 0.88; 95% CI 0.78-0.99; p = 0.037). We also performed a meta-analysis of these SNPs in East Asians, African ancestry and European ancestry populations and also observed a consistent effect (rs3778609, OR 0.95; 95% CI 0.91-0.97; p = 0.0017). CONCLUSION: Our study adds to evidence about the importance of the 6q25 locus for breast cancer susceptibility. Our finding also highlights the utility of performing additional searches for genetic variants for breast cancer in non-European populations.

Pathogenic and likely pathogenic variants in PALB2, CHEK2, and other known breast cancer susceptibility genes among 1054 BRCA-negative Hispanics with breast cancer.

BACKGROUND: Breast cancer (BC) is the most common cancer and related cause of mortality among Hispanics, yet susceptibility has been understudied. BRCA1 and BRCA2 (BRCA) mutations explain less than one-half of hereditary BC, and the proportion associated with other BC susceptibility genes is unknown. METHODS: Germline DNA from 1054 BRCA-mutation-negative Hispanic women with hereditary BC (BC diagnosed at age <51 years, bilateral BC, breast and ovarian cancer, or BC diagnosed at ages 51-70 years with ≥2 first-degree or second-degree relatives who had BC diagnosed at age <70 years), 312 local controls, and 887 multiethnic cohort controls was sequenced and analyzed for 12 known and suspected, high-penetrance and moderate-penetrance cancer susceptibility genes (ataxia telangiectasia mutated [ATM], breast cancer 1 interacting protein C-terminal helicase 1 [BRIP1], cadherin 1 [CDH1], checkpoint kinase 2 [CHEK2], nibrin [NBN], neurofibromatosis type 1 [NF1], partner and localizer of BRCA2 [PALB2], phosphatase and tensin homolog [PTEN], RAD51 paralog 3 [RAD51C], RAD51D, serine/threonine kinase 11 [STK11], and TP53). RESULTS: Forty-nine (4.6%) pathogenic or likely pathogenic variants (PVs) in 47 of 1054 participants (4.5%), including 21 truncating frameshift, 20 missense, 5 nonsense, and 4 splice variants, were identified in CHEK2 (n = 20), PALB2 (n = 18), ATM (n = 5), TP53 (n = 3), BRIP1 (n = 2), and CDH1 and NF1 (both n = 1) and none were identified in NBN, PTEN, STK11, RAD51C, or RAD51D. Nine participants carried the PALB2 c.2167_2168del PV (0.85%), and 14 carried the CHEK2 c.707T>C PV (1.32%). CONCLUSIONS: Of 1054 BRCA-negative, high-risk Hispanic women, 4.5% carried a PV in a cancer susceptibility gene, increasing understanding of hereditary BC in this population. Recurrent PVs in PALB2 and CHEK2 represented 47% (23 of 49) of the total, suggesting a founder effect. Accurate classification of variants was enabled by carefully controlling for ancestry and the increased identification of at-risk Hispanics for screening and prevention.

Whole-Genome Sequencing of Pharmacogenetic Drug Response in Racially Diverse Children with Asthma.

RATIONALE: Albuterol, a bronchodilator medication, is the first-line therapy for asthma worldwide. There are significant racial/ethnic differences in albuterol drug response. OBJECTIVES: To identify genetic variants important for bronchodilator drug response (BDR) in racially diverse children. METHODS: We performed the first whole-genome sequencing pharmacogenetics study from 1,441 children with asthma from the tails of the BDR distribution to identify genetic association with BDR. MEASUREMENTS AND MAIN RESULTS: We identified population-specific and shared genetic variants associated with BDR, including genome-wide significant (P < 3.53 × 10-7) and suggestive (P < 7.06 × 10-6) loci near genes previously associated with lung capacity (DNAH5), immunity (NFKB1 and PLCB1), and ß-adrenergic signaling (ADAMTS3 and COX18). Functional analyses of the BDR-associated SNP in NFKB1 revealed potential regulatory function in bronchial smooth muscle cells. The SNP is also an expression quantitative trait locus for a neighboring gene, SLC39A8. The lack of other asthma study populations with BDR and whole-genome sequencing data on minority children makes it impossible to perform replication of our rare variant associations. Minority underrepresentation also poses significant challenges to identify age-matched and population-matched cohorts of sufficient sample size for replication of our common variant findings. CONCLUSIONS: The lack of minority data, despite a collaboration of eight universities and 13 individual laboratories, highlights the urgent need for a dedicated national effort to prioritize diversity in research. Our study expands the understanding of pharmacogenetic analyses in racially/ethnically diverse populations and advances the foundation for precision medicine in at-risk and understudied minority populations.

Thyroid hormones and breast cancer association according to menopausal status and body mass index.

BACKGROUND: Thyroxine (T4) has been positively associated with tumor cell proliferation, while the effect of triiodothyronine (T3) on cell proliferation has not been well-established because it differs according to the type of cell line used. In Mexico, it has been reported that 14.5% of adult women have some type of thyroid dysfunction and abnormalities in thyroid function tests have been observed in a variety of non-thyroidal illnesses, including breast cancer (BC). These abnormalities might change with body mass index (BMI) because thyroid hormones are involved in the regulation of various metabolic pathways and probably by menopausal status because obesity has been negatively associated with BC in premenopausal women and has been positively associated with BC in postmenopausal women. METHODS: To assess the association between serum thyroid hormone concentration (T4 and T3) and BC and the influence of obesity as an effect modifier of this relationship in premenopausal and postmenopausal women, we measured serum thyroid hormone and thyroid antibody levels in 682 patients with incident breast cancer (cases) and 731 controls, who participated in a population-based case-control study performed from 2004 to 2007 in three states of Mexico. We tested the association of total T4 (TT4) and total T3 (TT3) stratifying by menopausal status and body mass index (BMI), and adjusted for other health and demographic risk factors using logistic regressions models. RESULTS: Higher serum total T4 (TT4) concentrations were associated with BC in both premenopausal (odds ratio (OR) per standard deviation = 5.98, 95% CI 3.01-11.90) and postmenopausal women (OR per standard deviation = 2.81, 95% CI 2.17-3.65). In premenopausal women, the effect of TT4 decreased as BMI increased while the opposite was observed in postmenopausal women. The significance of the effect modification was marginal (p = 0.059) in postmenopausal women and was not significant in premenopausal women (p = 0.22). Lower TT3 concentrations were associated with BC in both premenopausal and postmenopausal women and no effect modification was observed. CONCLUSIONS: There is a strong association between BC and serum concentrations of TT3 and TT4; this needs to be further investigated to understand why it happens and how important it is to consider these alterations in treatment.