Case-Case Genome-Wide Analyses Identify Subtype-Informative Variants that Confer Risk for Breast Cancer.

Breast cancer includes several subtypes with distinct characteristic biological, pathological, and clinical features. Elucidating subtype-specific genetic etiology could provide insights into the heterogeneity of breast cancer to facilitate development of improved prevention and treatment approaches. Here, we conducted pairwise case-case comparisons among five breast cancer subtypes by applying a case-case GWAS (CC-GWAS) approach to summary statistics data of the Breast Cancer Association Consortium. The approach identified 13 statistically significant loci and eight suggestive loci, the majority of which were identified from comparisons between triple-negative breast cancer (TNBC) and luminal A breast cancer. Associations of lead variants in 12 loci remained statistically significant after accounting for previously reported breast cancer susceptibility variants, among which two were genome-wide significant. Fine mapping implicated putative functional/causal variants and risk genes at several loci, e.g., 3q26.31/TNFSF10, 8q22.3/NACAP1/GRHL2, and 8q23.3/LINC00536/TRPS1, for TNBC as compared to luminal cancer. Functional investigation further identified rs16867605 at 8q22.3 as a SNP that modulates enhancer activity of GRHL2. Subtype-informative polygenic risk scores (PRS) were derived, and patients with a high subtype-informative PRS had an up to 2-fold increased risk of being diagnosed with TNBC instead of luminal cancers. The CC-GWAS PRS remained statistically significant after adjusting for TNBC PRS derived from traditional case-control GWAS in The Cancer Genome Atlas and the African Ancestry Breast Cancer Genetic Consortium. The CC-GWAS PRS was also associated with overall survival and disease-specific survival among breast cancer patients. Overall, these findings have advanced our understanding of the genetic etiology of breast cancer subtypes, particularly for TNBC.

A risk variant for Barrett's esophagus and esophageal adenocarcinoma at chr8p23.1 affects enhancer activity and implicates multiple gene targets.

Nineteen genetic susceptibility loci for esophageal adenocarcinoma (EAC) and its precursor Barrett's esophagus (BE) have been identified through genome-wide association studies (GWAS). Clinical translation of such discoveries, however, has been hindered by the slow pace of discovery of functional/causal variants and gene targets at these loci. We previously developed a systematic informatics pipeline to prioritize candidate functional variants using functional potential scores, applied the pipeline to select high-scoring BE/EAC risk loci and validated a functional variant at chr19p13.11 (rs10423674). Here, we selected two additional prioritized loci for experimental interrogation: chr3p13/rs1522552 and chr8p23.1/rs55896564. Candidate enhancer regions encompassing these variants were evaluated using luciferase reporter assays in two EAC cell lines. One of the two regions tested exhibited allele-specific enhancer activity - 8p23.1/rs55896564. CRISPR-mediated deletion of the putative enhancer in EAC cell lines correlated with reduced expression of three candidate gene targets: B lymphocyte kinase (BLK), nei like DNA glycosylase 2 (NEIL2) and cathepsin B (CTSB). Expression quantitative trait locus (eQTL) mapping in normal esophagus and stomach revealed strong associations between the BE/EAC risk allele at rs55896564 (G) and lower expression of CTSB, a protease gene implicated in epithelial wound repair. These results further support the utility of functional potential scores for GWAS variant prioritization, and provide the first experimental evidence of a functional variant and risk enhancer at the 8p23.1 GWAS locus. Identification of CTSB, BLK and NEIL2 as candidate gene targets suggests that altered expression of these genes may underlie the genetic risk association at 8p23.1 with BE/EAC.

Prioritization and functional analysis of GWAS risk loci for Barrett's esophagus and esophageal adenocarcinoma.

Genome-wide association studies (GWAS) have identified ~20 genetic susceptibility loci for esophageal adenocarcinoma (EAC), and its precursor, Barrett's esophagus (BE). Despite such advances, functional/causal variants and gene targets at these loci remain undefined, hindering clinical translation. A key challenge is that most causal variants map to non-coding regulatory regions such as enhancers, and typically, numerous potential candidate variants at GWAS loci require testing. We developed a systematic informatics pipeline for prioritizing candidate functional variants via integrative functional potential scores (FPS) consolidated from multi-omics annotations, and used this pipeline to identify two high-scoring variants for experimental interrogation: chr9q22.32/rs11789015 and chr19p13.11/rs10423674. Minimal candidate enhancer regions spanning these variants were evaluated using luciferase reporter assays in two EAC cell lines. One of the two variants tested (rs10423674) exhibited allele-specific enhancer activity. CRISPR-mediated deletion of the putative enhancer region in EAC cell lines correlated with reduced expression of two genes-CREB-regulated transcription coactivator 1 (CRTC1) and Cartilage oligomeric matrix protein (COMP); expression of five other genes remained unchanged (CRLF1, KLHL26, TMEM59L, UBA52, RFXANK). Expression quantitative trait locus mapping indicated that rs10423674 genotype correlated with CRTC1 and COMP expression in normal esophagus. This study represents the first experimental effort to bridge GWAS associations to biology in BE/EAC and supports the utility of FPS to guide variant prioritization. Our findings reveal a functional variant and candidate risk enhancer at chr19p13.11 and implicate CRTC1 and COMP as putative gene targets, suggesting that altered expression of these genes may underlie the BE/EAC risk association.

Multitrait genetic association analysis identifies 50 new risk loci for gastro-oesophageal reflux, seven new loci for Barrett's oesophagus and provides insights into clinical heterogeneity in reflux diagnosis.

OBJECTIVE: Gastro-oesophageal reflux disease (GERD) has heterogeneous aetiology primarily attributable to its symptom-based definitions. GERD genome-wide association studies (GWASs) have shown strong genetic overlaps with established risk factors such as obesity and depression. We hypothesised that the shared genetic architecture between GERD and these risk factors can be leveraged to (1) identify new GERD and Barrett's oesophagus (BE) risk loci and (2) explore potentially heterogeneous pathways leading to GERD and oesophageal complications. DESIGN: We applied multitrait GWAS models combining GERD (78 707 cases; 288 734 controls) and genetically correlated traits including education attainment, depression and body mass index. We also used multitrait analysis to identify BE risk loci. Top hits were replicated in 23andMe (462 753 GERD cases, 24 099 BE cases, 1 484 025 controls). We additionally dissected the GERD loci into obesity-driven and depression-driven subgroups. These subgroups were investigated to determine how they relate to tissue-specific gene expression and to risk of serious oesophageal disease (BE and/or oesophageal adenocarcinoma, EA). RESULTS: We identified 88 loci associated with GERD, with 59 replicating in 23andMe after multiple testing corrections. Our BE analysis identified seven novel loci. Additionally we showed that only the obesity-driven GERD loci (but not the depression-driven loci) were associated with genes enriched in oesophageal tissues and successfully predicted BE/EA. CONCLUSION: Our multitrait model identified many novel risk loci for GERD and BE. We present strong evidence for a genetic underpinning of disease heterogeneity in GERD and show that GERD loci associated with depressive symptoms are not strong predictors of BE/EA relative to obesity-driven GERD loci.

Clinical and Lifestyle-Related Prognostic Indicators among Esophageal Adenocarcinoma Patients Receiving Treatment at a Comprehensive Cancer Center.

PURPOSE: The incidence of esophageal adenocarcinoma (EAC) has risen substantially in recent decades, while the average 5-year survival remains only ~20%. Disease stage and treatment are the strongest prognostic factors. The role of lifestyle factors in relation to survival remains uncertain, with a handful of studies to date investigating associations with obesity, smoking, physical activity, diet, or medications. METHODS: This study included patients diagnosed with primary adenocarcinoma of the esophagus, gastroesophageal junction, or cardia (N = 371) at Roswell Park Comprehensive Cancer Center between 2003 and 2019. Leveraging extensive data abstracted from electronic medical records, epidemiologic questionnaires, and a tumor registry, we analyzed clinical, behavioral, and environmental exposures and evaluated stage-specific associations with survival. Survival distributions were visualized using Kaplan-Meier curves. Cox proportional hazards regression models adjusted for age, sex, stage, treatment, and comorbidities were used to estimate the association between each exposure and all-cause or cancer-specific mortality. RESULTS: Among patients presenting with localized/regional tumors (stages I-III), current smoking was associated with increased overall mortality risk (HR = 2.5 [1.42-4.53], p = 0.002), while current physical activity was linked to reduced risk (HR = 0.58 [0.35-0.96], p = 0.035). Among patients with stage IV disease, individuals reporting pre-diagnostic use of statins (HR = 0.62 [0.42-0.92], p = 0.018) or NSAIDs (HR = 0.61 [0.42-0.91], p = 0.016) had improved overall survival. Exploratory analyses suggested that high pre-diagnostic dietary consumption of broccoli, carrots, and fiber correlated with prolonged overall survival in patients with localized/regional disease. CONCLUSION: Our data suggest that lifestyle exposures may be differentially associated with EAC survival based on disease stage. Future investigation of larger, diverse patient cohorts is essential to validate these findings. Our results may help inform the development of lifestyle-based interventions to improve EAC prognosis and quality of life.

Quantitative global lipidomics analysis of patients with ovarian cancer versus benign adnexal mass.

Altered lipid metabolism has emerged as an important feature of ovarian cancer (OC), yet the translational potential of lipid metabolites to aid in diagnosis and triage remains unproven. We conducted a multi-level interrogation of lipid metabolic phenotypes in patients with adnexal masses, integrating quantitative lipidomics profiling of plasma and ascites with publicly-available tumor transcriptome data. Using Sciex Lipidyzer, we assessed concentrations of > 500 plasma lipids in two patient cohorts-(i) a pilot set of 100 women with OC (50) or benign tumor (50), and (ii) an independent set of 118 women with malignant (60) or benign (58) adnexal mass. 249 lipid species and several lipid classes were significantly reduced in cases versus controls in both cohorts (FDR < 0.05). 23 metabolites-triacylglycerols, phosphatidylcholines, cholesterol esters-were validated at Bonferroni significance (P < 9.16 × 10-5). Certain lipids exhibited greater alterations in early- (diacylglycerols) or late-stage (lysophospholipids) cases, and multiple lipids in plasma and ascites were positively correlated. Lipoprotein receptor gene expression differed markedly in OC versus benign tumors. Importantly, several plasma lipid species, such as DAG(16:1/18:1), improved the accuracy of CA125 in differentiating early-stage OC cases from benign controls, and conferred a 15-20% increase in specificity at 90% sensitivity in multivariate models adjusted for age and BMI. This study provides novel insight into systemic and local lipid metabolic differences between OC and benign disease, further implicating altered lipid uptake in OC biology, and advancing plasma lipid metabolites as a complementary class of circulating biomarkers for OC diagnosis and triage.

Genes Relevant to Tissue Response to Cancer Therapy Display Diurnal Variation in mRNA Expression in Human Oral Mucosa.

BACKGROUND: To address a critical gap for application of cancer chronotherapy of when would be the best time(s) for treating an individual cancer patient, we conducted a pilot study to characterize diurnal variations of gene expression in oral mucosal tissue, which is vulnerable to damage from cancer therapies. METHODS: We conducted RNA-seq assay on individual oral mucosal samples collected from 11 healthy volunteers every 4 hours (6 time points). Using a cosine-based method, we estimated the individual and average values of peak-time and amplitude for each gene. Correlations between gene expression peak-times and age was examined, adjusting for individual's sleep timing. RESULTS: Among candidate gene pathways that are relevant to treatment response, 7 of 16 genes (PER3, CIART, TEF, PER1, PER2, CRY2, ARNTL) involved in circadian regulation and 1 of 118 genes (WEE1) involved in cell cycle regulation achieved p-value ≤ 0.1 and relative amplitude>0.1. The average peak times were approximately 10:15 for PER3, CIART and TEF, 10:45 for PER1, 13:00 for WEE1, PER2 and CRY2, and 19:30 for ARNTL. Ranges in peak times across individuals differed by gene (e.g., 8 hours for PER1; 16.7 hours for WEE1). Older people had later peak times for PER1 (r = 0.77, p = 0.03) and PER3 (r = 0.69, p-value = 0.06). CONCLUSION: In oral mucosa, expression of some genes relevant to treatment response displayed diurnal variation. These genes may be candidates for development of biomarkers for optimizing individual timing of cancer therapy using non-invasively collected oral mucosa.

Deletion of Foxa1 in the mouse mammary gland results in abnormal accumulation of luminal progenitor cells: a link between reproductive factors and ER-/TNBC breast cancer?

In humans, parity without breastfeeding increases risk of estrogen receptor-negative (ER-) breast cancer and is associated with hypermethylation of FOXA1, a pioneer factor regulating lineage commitment of mammary gland luminal progenitor cells. We postulate that pregnancy-associated repression of FOXA1 results in the accumulation of aberrant, differentiation-arrested luminal progenitor cells which, following additional genetic and epigenetic insults, may give rise to ER- tumors. Consistent with this hypothesis, we show that deletion of Foxa1 in the mouse mammary gland results in a two-fold increase in the proportion of luminal progenitor cells and a reduction in mammary gland epithelial cells that stain positive for ER. These results provide compelling support for the notion that reduced Foxa1 expression is sufficient to alter mammary gland luminal cell fate determination in vivo, which could be a mechanism linking parity with ER- breast cancer.

Circulating CD14+ HLA-DRlo/- monocytic cells as a biomarker for epithelial ovarian cancer progression.

PROBLEM: Previous studies identified circulating CD14+ HLA-DRlo/- monocytic cells as an immune suppressive subset in solid malignancies, such as prostate, renal cell carcinoma, and pancreatic cancer. Such monocytic cells have been implicated not only in tumour progression but also as a potential barrier for immunotherapy. This study examined the relationship between the frequency of circulating monocytic cells and epithelial ovarian cancer (EOC) progression pre- and post-frontline chemotherapy, defined by disease stage, which is a leading prognostic factor for this malignancy. METHOD OF STUDY: Incident cases of 236 women with EOC were recruited and comprehensive flow cytometry was utilized to assess the frequency of peripheral blood CD33+ CD11b+ HLA-DR-/low CD14+ CD15- monocytic cells, henceforth termed CD14+ HLA-DRlo/- monocytic cells, prior to and after completion of frontline chemotherapy. Multivariable odds ratios (OR) were used to estimate the association between CD14+ HLA-DRlo/- monocytic cell percentages and disease stage. Wilcoxon signed-rank tests evaluated changes in these monocytic cell levels pre- and post-chemotherapy in a patient subset (n = 70). RESULTS: Patients with elevated frequencies of circulating CD14+ HLA-DRlo/- monocytic cells at diagnosis were at 3.33-fold greater odds of having advanced stage (III/IV) EOC (CI: 1.04-10.64), with a significant trend in increasing CD14+ HLA-DRlo/- monocytic cell levels (P = .04). There was a 2.02% median decrease of these monocytic cells post-chemotherapy among a subset of patients with advanced stage disease (P < .0001). CONCLUSION: These findings support the potential clinical relevance of CD14+ HLA-DRlo/- monocytic cells in EOC for prognosis and may indicate a non-invasive biomarker to measure disease progression.

FOXA1 Protein Expression in ER+ and ER- Breast Cancer in Relation to Parity and Breastfeeding in Black and White Women.

BACKGROUND: Forkhead box protein A1 (FOXA1) promotes luminal differentiation, and hypermethylation of the gene can be a mechanism of developing estrogen receptor-negative (ER-) breast cancer. We examined FOXA1 in breast tumor and adjacent normal tissue in relation to reproductive factors, particularly higher parity and no breastfeeding, that are associated with ER- tumors. METHODS: We performed IHC for FOXA1 in breast tumors (n = 1,329) and adjacent normal tissues (n = 298) in the Women's Circle of Health Study (949 Blacks and 380 Whites). Protein expression levels were summarized by histology (H) scores. Generalized linear models were used to assess FOXA1 protein expression in relation to reproductive factors by ER status. RESULTS: ER-positive (ER+) versus ER- tumors had higher FOXA1 protein expression (P < 0.001). FOXA1 expression was higher in tumor versus paired adjacent normal tissue in women with ER+ or non-triple-negative cancer (both P < 0.001), but not in those with ER- or triple-negative cancer. Higher number of births (1, 2, and 3+) was associated with lower FOXA1 protein expression in ER+ tumors [differences in H score, or ß = -8.5; 95% confidence interval (CI), -15.1 to -2.0], particularly among parous women who never breastfed (ß = -10.4; 95% CI, -19.7 to -1.0), but not among those who breastfed (ß = -7.5; 95% CI, -16.9 to 1.8). The associations for ER- tumors were similar, although they were not statistically significant. CONCLUSIONS: In this tumor-based study, higher parity was associated with lower FOXA1 expression in ER+ tumors, and breastfeeding may ameliorate the influence. IMPACT: These findings contribute to our understanding of FOXA1 methylation and breast cancer etiology.

Survival disparities among racial/ethnic groups of women with ovarian cancer: An update on data from the Surveillance, Epidemiology and End Results (SEER) registry.

OBJECTIVE: Update information on racial disparities in ovarian cancer survival from the Surveillance, Epidemiology, and End Results (SEER) Program. METHODS: Data on women with epithelial ovarian cancer from the SEER Program between 1995-2015 were collected including; patient ID, age at diagnosis, year of diagnosis, surgery, chemotherapy, radiation, insurance status, region of registry, tumor grade, tumor histology, tumor summary stage, survival months, race/ethnicity, and vital status. Multivariable analyses were performed to examine overall survival, differences in survival by age at diagnosis, by year of diagnosis, risk of not receiving surgery, and risk of 12-month death across racial/ethnic groups. RESULTS: Non-Hispanic black women (n = 4261) had an increased risk of overall mortality (HR = 1.28, CI: 1.23-1.33) when compared to non-Hispanic white women (n = 47,475), which appears more pronounced among women diagnosed under age 50. Hispanic women (n = 7052) had no difference in survival when compared to non-Hispanic white women (HR = 1.03, CI: 0.99-1.07). Non-Hispanic Asian/PI women (n = 5008) exhibited slightly reduced risk (HR = 0.95, CI: 0.91-0.99) when compared to non-Hispanic white women. Risk of not receiving surgical intervention remains high among non-Hispanic black women and Hispanic women, when compared to non-Hispanic white women. Non-Hispanic black women, non-Hispanic Asian/PI women, and Hispanic women were all at significantly greater risk of dying within the first 12 months of cancer diagnosis when compared to non-Hispanic white women. CONCLUSION: Disparities in survival remain across various racial/ethnic groups, when compared to non-Hispanic white women with ovarian cancer. These disparities should continue to be examined in an effort to decrease such gaps.

Differences in microRNA expression in breast cancer between women of African and European ancestry.

Breast cancer is a heterogeneous disease, characterized by molecularly and phenotypically distinct tumor subtypes, linked to disparate clinical outcomes. American women of African ancestry (AA) are more likely than those of European ancestry (EA) to be diagnosed with aggressive, estrogen receptor negative (ER-) or triple negative breast cancer, and to die of this disease. However, the underlying causes of AA predisposition to ER-/triple negative breast cancer are still largely unknown. In this study, we performed high-throughput whole-genome miRNA expression profiling in breast tissue samples from both AA and EA women. A number of differentially expressed miRNAs, i.e., DEmiRs defined as >2-fold change in expression and false discovery rate <0.05, were identified as up- or downregulated by tumor ER status or by ancestry. We found that among 102 ER-subtype-related DEmiRs identified in breast tumors, the majority of these DEmiRs were race specific, with only 23 DEmiRs shared in tumors from both AAs and EAs; this finding indicates that there are unique subsets of miRNAs differentially expressed between ER- and ER positive tumors within AAs versus EAs. Our overall results support the notion that miRNA expression patterns may differ not only by tumor subtype but by ancestry, indicating differences in tumor biology and heterogeneity of breast cancer between AAs and EAs. These results will provide the basis for further functional analysis to elucidate biological differences between AAs and EAs and to help develop targeted treatment strategies to reduce disparities in breast cancer.

Whole-genome sequencing of esophageal adenocarcinoma in Chinese patients reveals distinct mutational signatures and genomic alterations.

While the incidence of esophageal adenocarcinoma (EAC) has risen drastically in Western countries over the last 40 years, a similar trend has not been observed for EAC in China. Here, we analyzed mutational spectrum, copy number alterations, and structural variants from whole-genome sequencing of 10 Chinese EAC tumor samples and their matched normal samples, and compared them to previously reported EAC tumor specimens from Western countries. The mutational burden in Chinese EAC was significantly lower than that found in EAC from Western countries. The hallmark A>C mutational signature observed at high frequency in EAC from Western countries, which has been linked to acid reflux, is completely absent in Chinese samples. Furthermore, none of the Chinese samples showed evidence of chromothripsis and genome doubling that are often found in EAC from Western countries. In summary, Chinese EAC tumor samples had distinct genomic profiles and signatures, suggesting that EAC in Chinese individuals may arise from a different etiological pathway.

Salivary metabolite profiling distinguishes patients with oral cavity squamous cell carcinoma from normal controls.

Oral cavity squamous cell carcinoma (OCC) and oropharyngeal squamous cell carcinoma (OPC) are among the most common cancers worldwide and are associated with high mortality and morbidity. The purpose of this study is to identify potential biomarkers to distinguish OCC/OPC from normal controls and to distinguish OCC patients with and without nodal metastasis. We tested saliva samples from 101 OCC, 58 OPC, and 35 normal controls using four analytical platforms (NMR, targeted aqueous by LC-MS/MS, global aqueous and global lipidomics by LC-Q-TOF). Samples from OCC and normal controls were divided into discovery and validation sets. Using linear regression adjusting for age, sex, race and experimental batches, we found the levels of two metabolites (glycine and proline) to be significantly different between OCC and controls (FDR < 0.1 for both discovery and validation sets) but did not find any appreciable differences in metabolite levels between OPC and controls or between OCC with and without nodal metastasis. Four metabolites, including glycine, proline, citrulline, and ornithine were associated with early stage OCC in both discovery and validation sets. Further study is warranted to confirm these results in the development of salivary metabolites as diagnostic markers.

Recommendation to use exact P-values in biomarker discovery research in place of approximate P-values.

BACKGROUND: Biomarker candidates are often ranked using P-values. Standard P-value calculations use normal or logit-normal approximations, which may not be correct for small P-values and small sample sizes common in discovery research. METHODS: We compared exact P-values, correct by definition, with logit-normal approximations in a simulated study of 40 cases and 160 controls. The key measure of biomarker performance was sensitivity at 90% specificity. Data for 3000 uninformative false markers and 30 informative true markers were generated randomly. We also analyzed real data for 2371 plasma protein markers measured in 121 breast cancer cases and 121 controls. RESULTS: In our simulation, using the same discovery criterion, exact P-values led to discovery of 24 true and 82 false biomarkers, while logit-normal approximate P-values yielded 20 true and 106 false biomarkers. The estimated true discovery rate was substantially off for approximate P-values: logit-normal estimated 42 but found 20. The exact method estimated 22, very close to 24, which was the actual number of true discoveries. Although these results are based on one specific simulation, qualitatively similar results were obtained from 10 random repetitions. With real data, ranking candidate biomarkers by exact P-values, versus approximate P-values, resulted in a very different ordering of these markers. CONCLUSIONS: Exact P-values, which correspond to permutation tests with non-parametric rank statistics such as empirical ROC statistics, are preferred over approximate P-values. Approximate P-values can lead to inappropriate biomarker selection rules and incorrect conclusions. IMPACT: Exact P-values in place of approximate P-values in discovery research may improve the yield of biomarkers that validate clinically.

Interactions Between Genetic Variants and Environmental Factors Affect Risk of Esophageal Adenocarcinoma and Barrett's Esophagus.

BACKGROUND & AIMS: Genome-wide association studies (GWAS) have identified more than 20 susceptibility loci for esophageal adenocarcinoma (EA) and Barrett's esophagus (BE). However, variants in these loci account for a small fraction of cases of EA and BE. Genetic factors might interact with environmental factors to affect risk of EA and BE. We aimed to identify single nucleotide polymorphisms (SNPs) that may modify the associations of body mass index (BMI), smoking, and gastroesophageal reflux disease (GERD), with risks of EA and BE. METHODS: We collected data on single BMI measurements, smoking status, and symptoms of GERD from 2284 patients with EA, 3104 patients with BE, and 2182 healthy individuals (controls) participating in the Barrett's and Esophageal Adenocarcinoma Consortium GWAS, the UK Barrett's Esophagus Gene Study, and the UK Stomach and Oesophageal Cancer Study. We analyzed 993,501 SNPs in DNA samples of all study subjects. We used standard case-control logistic regression to test for gene-environment interactions. RESULTS: For EA, rs13429103 at chromosome 2p25.1, near the RNF144A-LOC339788 gene, showed a borderline significant interaction with smoking status (P = 2.18×10-7). Ever smoking was associated with an almost 12-fold increase in risk of EA among individuals with rs13429103-AA genotype (odds ratio=11.82; 95% CI, 4.03-34.67). Three SNPs (rs12465911, rs2341926, rs13396805) at chromosome 2q23.3, near the RND3-RBM43 gene, interacted with GERD symptoms (P = 1.70×10-7, P = 1.83×10-7, and P = 3.58×10-7, respectively) to affect risk of EA. For BE, rs491603 at chromosome 1p34.3, near the EIF2C3 gene, and rs11631094 at chromosome 15q14, at the SLC12A6 gene, interacted with BMI (P = 4.44×10-7) and pack-years of smoking history (P = 2.82×10-7), respectively. CONCLUSION: The associations of BMI, smoking, and GERD symptoms with risks of EA and BE appear to vary with SNPs at chromosomes 1, 2, and 15. Validation of these suggestive interactions is warranted.

Determining Risk of Barrett's Esophagus and Esophageal Adenocarcinoma Based on Epidemiologic Factors and Genetic Variants.

BACKGROUND & AIMS: We developed comprehensive models to determine risk of Barrett's esophagus (BE) or esophageal adenocarcinoma (EAC) based on genetic and non-genetic factors. METHODS: We used pooled data from 3288 patients with BE, 2511 patients with EAC, and 2177 individuals without either (controls) from participants in the international Barrett's and EAC consortium as well as the United Kingdom's BE gene study and stomach and esophageal cancer study. We collected data on 23 genetic variants associated with risk for BE or EAC, and constructed a polygenic risk score (PRS) for cases and controls by summing the risk allele counts for the variants weighted by their natural log-transformed effect estimates (odds ratios) extracted from genome-wide association studies. We also collected data on demographic and lifestyle factors (age, sex, smoking, body mass index, use of nonsteroidal anti-inflammatory drugs) and symptoms of gastroesophageal reflux disease (GERD). Risk models with various combinations of non-genetic factors and the PRS were compared for their accuracy in identifying patients with BE or EAC using the area under the receiver operating characteristic curve (AUC) analysis. RESULTS: Individuals in the highest quartile of risk, based on genetic factors (PRS), had a 2-fold higher risk of BE (odds ratio, 2.22; 95% confidence interval, 1.89-2.60) or EAC (odds ratio, 2.46; 95% confidence interval, 2.07-2.92) than individual in the lowest quartile of risk based on PRS. Risk models developed based on only demographic or lifestyle factors or GERD symptoms identified patients with BE or EAC with AUC values ranging from 0.637 to 0.667. Combining data on demographic or lifestyle factors with data on GERD symptoms identified patients with BE with an AUC of 0.793 and patients with EAC with an AUC of 0.745. Including PRSs with these data only minimally increased the AUC values for BE (to 0.799) and EAC (to 0.754). Including the PRSs in the model developed based on non-genetic factors resulted in a net reclassification improvement for BE of 3.0% and for EAC of 5.6%. CONCLUSIONS: We used data from 3 large databases of patients from studies of BE or EAC to develop a risk prediction model based on genetic, clinical, and demographic/lifestyle factors. We identified a PRS that increases discrimination and net reclassification of individuals with vs without BE and EAC. However, the absolute magnitude of improvement is not sufficient to justify its clinical use.

FOXA1 hypermethylation: link between parity and ER-negative breast cancer in African American women?

BACKGROUND: Reproductive factors, particularly parity, have differential effects on breast cancer risk according to estrogen receptor (ER) status, especially among African American (AA) women. One mechanism could be through DNA methylation, leading to altered expression levels of genes important in cell fate decisions. METHODS: Using the Illumina 450K BeadChip, we compared DNA methylation levels in paraffin-archived tumor samples from 383 AA and 350 European American (EA) women in the Women's Circle of Health Study (WCHS). We combined 450K profiles with RNA-seq data and prioritized genes based on differential methylation by race, correlation between methylation and gene expression, and biological function. We measured tumor protein expression and assessed its relationship to DNA methylation. We evaluated associations between reproductive characteristics and DNA methylation using linear regression. RESULTS: 410 loci were differentially methylated by race, with the majority unique to ER- tumors. FOXA1 was hypermethylated in tumors from AA versus EA women with ER- cancer, and increased DNA methylation correlated with reduced RNA and protein expression. Importantly, parity was positively associated with FOXA1 methylation among AA women with ER- tumors (P = 0.022), as was number of births (P = 0.026), particularly among those who did not breastfeed (P = 0.008). These same relationships were not observed among EA women, although statistical power was more limited. CONCLUSIONS: Methylation and expression of FOXA1 is likely impacted by parity and breastfeeding. Because FOXA1 regulates a luminal gene expression signature in progenitor cells and represses the basal phenotype, this could be a mechanism that links these reproductive exposures with ER- breast cancer.

Candidate serum metabolite biomarkers for differentiating gastroesophageal reflux disease, Barrett's esophagus, and high-grade dysplasia/esophageal adenocarcinoma.

INTRODUCTION/OBJECTIVES: Incidence of esophageal adenocarcinoma (EA), an often fatal cancer, has increased sharply over recent decades. Several important risk factors (reflux, obesity, smoking) have been identified for EA and its precursor, Barrett's esophagus (BE), but a key challenge remains identifying individuals at highest risk, since most with reflux do not develop BE, and most with BE do not progress to cancer. Metabolomics represents an emerging approach for identifying novel biomarkers associated with cancer development. METHODS: We used targeted liquid chromatography-mass spectrometry (LC-MS) to profile 57 metabolites in 322 serum specimens derived from individuals with gastroesophageal reflux disease (GERD), BE, high-grade dysplasia (HGD), or EA, drawn from two well-annotated epidemiologic parent studies. RESULTS: Multiple metabolites differed significantly (P<0.05) between BE versus GERD (n=9), and between HGD/EA versus BE (n=4). Several top candidates (FDR q≤0.15), including urate, homocysteine, and 3-nitrotyrosine, are linked to inflammatory processes, which may contribute to BE/EA pathogenesis. Multivariate modeling achieved moderate discrimination between HGD/EA and BE (AUC=0.75), with less pronounced separation for BE versus GERD (AUC=0.64). CONCLUSION: Serum metabolite differences can be detected between individuals with GERD versus BE, and between those with BE versus HGD/EA, and may help differentiate patients at different stages of progression to EA.

Germline variation in inflammation-related pathways and risk of Barrett's oesophagus and oesophageal adenocarcinoma.

OBJECTIVE: Oesophageal adenocarcinoma (OA) incidence has risen sharply in Western countries over recent decades. Local and systemic inflammation is considered an important contributor to OA pathogenesis. Established risk factors for OA and its precursor, Barrett's oesophagus (BE), include symptomatic reflux, obesity and smoking. The role of inherited genetic susceptibility remains an area of active investigation. Here, we explore whether germline variation related to inflammatory processes influences susceptibility to BE/OA. DESIGN: We used data from a genomewide association study of 2515 OA cases, 3295 BE cases and 3207 controls. Our analysis included 7863 single-nucleotide polymorphisms (SNPs) in 449 genes assigned to five pathways: cyclooxygenase (COX), cytokine signalling, oxidative stress, human leucocyte antigen and nuclear factor-κB. A principal components-based analytic framework was employed to evaluate pathway-level and gene-level associations with disease risk. RESULTS: We identified a significant signal for the COX pathway in relation to BE risk (p=0.0059, false discovery rate q=0.03), and in gene-level analyses found an association with microsomal glutathione-S-transferase 1 (MGST1); (p=0.0005, q=0.005). Assessment of 36 MGST1 SNPs identified 14 variants associated with elevated BE risk (q<0.05). Four of these were subsequently confirmed (p<5.5×10-5) in a meta-analysis encompassing an independent set of 1851 BE cases and 3496 controls, and are known strong expression quantitative trait loci for MGST1. Three such variants were associated with similar elevations in OA risk. CONCLUSIONS: This study provides the most comprehensive evaluation of inflammation-related germline variation in relation to risk of BE/OA and suggests that variants in MGST1 influence disease susceptibility.