ABSTRACT
Polygenic risk scores (PRSs) are useful for predicting breast cancer risk, but the prediction accuracy of existing PRSs in women of African ancestry (AA) remains relatively low. We aim to develop optimal PRSs for the prediction of overall and estrogen receptor (ER) subtype-specific breast cancer risk in AA women. The AA dataset comprised 9235 cases and 10 184 controls from four genome-wide association study (GWAS) consortia and a GWAS study in Ghana. We randomly divided samples into training and validation sets. We built PRSs using individual-level AA data by a forward stepwise logistic regression and then developed joint PRSs that combined (1) the PRSs built in the AA training dataset and (2) a 313-variant PRS previously developed in women of European ancestry. PRSs were evaluated in the AA validation set. For overall breast cancer, the odds ratio per standard deviation of the joint PRS in the validation set was 1.34 [95% confidence interval (CI): 1.27-1.42] with the area under receiver operating characteristic curve (AUC) of 0.581. Compared with women with average risk (40th-60th PRS percentile), women in the top decile of the PRS had a 1.98-fold increased risk (95% CI: 1.63-2.39). For PRSs of ER-positive and ER-negative breast cancer, the AUCs were 0.608 and 0.576, respectively. Compared with existing methods, the proposed joint PRSs can improve prediction of breast cancer risk in AA women.
Subject(s)
Breast Neoplasms , Genome-Wide Association Study , Breast Neoplasms/genetics , Female , Genetic Predisposition to Disease , Humans , Multifactorial Inheritance/genetics , Receptors, Estrogen/genetics , Risk FactorsABSTRACT
The role of minor histocompatibility antigens (mHAs) in mediating graft versus leukemia and graft versus host disease (GvHD) following allogeneic hematopoietic cell transplantation (alloHCT) is recognized but not well-characterized. By implementing improved methods for mHA prediction in two large patient cohorts, this study aimed to comprehensively explore the role of mHAs in alloHCT by analyzing whether (1) the number of predicted mHAs, or (2) individual mHAs are associated with clinical outcomes. The study population consisted of 2249 donor-recipient pairs treated for acute myeloid leukemia and myelodysplastic syndrome with alloHCT. A Cox proportional hazard model showed that patients with a class I mHA count greater than the population median had an increased hazard of GvHD mortality (hazard ratio [HR] = 1.39, 95% confidence interval [CI] = 1.01, 1.77, p = .046). Competing risk analyses identified the class I mHAs DLRCKYISL (GSTP), WEHGPTSLL (CRISPLD2), and STSPTTNVL (SERPINF2) were associated with increased GVHD mortality (HR = 2.84, 95% CI = 1.52, 5.31, p = .01), decreased leukemia-free survival (LFS) (HR = 1.94, 95% CI = 1.27, 2.95, p = .044), and increased disease-related mortality (DRM) (HR = 2.32, 95% CI = 1.5, 3.6, p = .008), respectively. One class II mHA YQEIAAIPSAGRERQ (TACC2) was associated with increased risk of treatment-related mortality (TRM) (HR = 3.05, 95% CI = 1.75, 5.31, p = .02). WEHGPTSLL and STSPTTNVL were both present within HLA haplotype B*40:01-C*03:04 and showed a positive dose-response relationship with increased all-cause mortality and DRM and decreased LFS, indicating these two mHAs contribute to the risk of mortality in an additive manner. Our study reports the first large-scale investigation of the associations of predicted mHA peptides with clinical outcomes following alloHCT.
Subject(s)
Graft vs Host Disease , Hematopoietic Stem Cell Transplantation , Leukemia, Myeloid, Acute , Humans , Minor Histocompatibility Antigens/genetics , Transplantation, Homologous/adverse effects , Hematopoietic Stem Cell Transplantation/adverse effects , Hematopoietic Stem Cell Transplantation/methods , Leukemia, Myeloid, Acute/therapy , Retrospective StudiesABSTRACT
BACKGROUND: Gauging fitness remains a challenge among older adults with hematologic malignancies, and interventions to restore function are lacking. We pilot a structured exercise intervention and novel biologic correlates of aging using epigenetic clocks and markers of immunosenescence to evaluate changes in function and clinical outcomes. METHODS: Older adults (n=30) with hematologic malignancy actively receiving treatment were screened and enrolled in a 6-month exercise intervention, the Otago Exercise Programme (OEP). The impact of the OEP on geriatric assessment metrics and health-related quality of life were captured. Clinical outcomes of overall survival and hospital utilization (inpatient length of stay and emergency department use) in relationship to geriatric deficits were analyzed. RESULTS: Older adults (median age, 75.5 years [range, 62-83 years]) actively receiving treatment were enrolled in the OEP. Instrumental activities of daily living and physical health scores (PHS) increased significantly with the OEP intervention (median PHS: visit 1, 55 [range, 0-100]; visit 2, 70 [range, 30-100]; P<.01). Patient-reported Karnofsky performance status increased significantly, and the improvement was sustained (median [range]: visit 1, 80 [40-100]; visit 3, 90 [50-100]; P=.05). Quality of life (Patient-Reported Outcome Measurement Information System [PROMIS]) improved significantly by the end of the 6-month period (median [range]: visit 1, 32.4 [19.9-47.7]; visit 3, 36.2 [19.9-47.7]; P=.01]. Enhanced measures of gait speed and balance, using the Short Physical Performance Battery scores, were associated with a 20% decrease in risk of death (hazard ratio, 0.80; 95% CI, 0.65-0.97; P=.03) and a shorter hospital length of stay (decrease of 1.29 days; 95% CI, -2.46 to -0.13; P=.03). Peripheral blood immunosenescent markers were analyzed in relationship to clinical frailty and reports of mPhenoAge epigenetic analysis are preliminarily reported. Chronologic age had no relationship to overall survival, length of stay, or emergency department utilization. CONCLUSIONS: The OEP was effective in improving quality of life, and geriatric tools predicted survival and hospital utilization among older adults with hematologic malignancies.
Subject(s)
Activities of Daily Living , Hematologic Neoplasms , Aged , Aging , Geriatric Assessment , Hematologic Neoplasms/therapy , Humans , Phenotype , Quality of LifeABSTRACT
Selection pressure due to exposure to infectious pathogens endemic to Africa may explain distinct genetic variations in immune response genes. However, the impact of those genetic variations on human immunity remains understudied, especially within the context of modern lifestyles and living environments, which are drastically different from early humans in sub Saharan Africa. There are few data on population differences in constitutional immune environment, where genetic ancestry and environment are likely two primary sources of variation. In a study integrating genetic, molecular and epidemiologic data, we examined population differences in plasma levels of 14 cytokines involved in innate and adaptive immunity, including those implicated in chronic inflammation, and possible contributing factors to such differences, in 914 AA and 855 EA women. We observed significant differences in 7 cytokines, including higher plasma levels of CCL2, CCL11, IL4 and IL10 in EAs and higher levels of IL1RA and IFNα2 in AAs. Analyses of a wide range of demographic and lifestyle factors showed significant impact, with age, education level, obesity, smoking, and alcohol intake, accounting for some, but not all, observed population differences for the cytokines examined. Levels of two pro-inflammatory chemokines, CCL2 and CCL11, were strongly associated with percent of African ancestry among AAs. Through admixture mapping, the signal was pinpointed to local ancestry at 1q23, with fine-mapping analysis refined to the Duffy-null allele of rs2814778. In AA women, this variant was a major determinant of systemic levels of CCL2 (p = 1.1e-58) and CCL11 (p = 2.2e-110), accounting for 19% and 40% of the phenotypic variance, respectively. Our data reveal strong ancestral footprints in inflammatory chemokine regulation. The Duffy-null allele may indicate a loss of the buffering function for chemokine levels. The substantial immune differences by ancestry may have broad implications to health disparities between AA and EA populations.
Subject(s)
Adaptation, Biological/genetics , Cytokines/genetics , Gene Expression Regulation , Genetic Variation , Selection, Genetic , Adaptive Immunity/genetics , Adult , Alleles , Biological Evolution , Black People/genetics , Cytokines/blood , Duffy Blood-Group System/genetics , Environment , Female , Gene Frequency , Health Status Disparities , Healthy Volunteers , Humans , Immunity, Innate/genetics , Middle Aged , White People/geneticsABSTRACT
SUMMARY: To address the limited software options for performing survival analyses with millions of SNPs, we developed gwasurvivr, an R/Bioconductor package with a simple interface for conducting genome-wide survival analyses using VCF (outputted from Michigan or Sanger imputation servers), IMPUTE2 or PLINK files. To decrease the number of iterations needed for convergence when optimizing the parameter estimates in the Cox model, we modified the R package survival; covariates in the model are first fit without the SNP, and those parameter estimates are used as initial points. We benchmarked gwasurvivr with other software capable of conducting genome-wide survival analysis (genipe, SurvivalGWAS_SV and GWASTools). gwasurvivr is significantly faster and shows better scalability as sample size, number of SNPs and number of covariates increases. AVAILABILITY AND IMPLEMENTATION: gwasurvivr, including source code, documentation and vignette are available at: http://bioconductor.org/packages/gwasurvivr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Genome , Software , Polymorphism, Single Nucleotide , Survival AnalysisABSTRACT
Although survival outcomes have significantly improved, up to 40% of patients die within 1 year of HLA-matched unrelated-donor blood and marrow transplantation (BMT). To identify non-HLA genetic contributors to mortality after BMT, we performed the first exome-wide association study in the DISCOVeRY-BMT cohorts using the Illumina HumanExome BeadChip. This study includes 2473 patients with acute myeloid leukemia, acute lymphoblastic leukemia, or myelodysplastic syndrome and 2221 10/10 HLA-matched donors treated from 2000 to 2011. Single-variant and gene-level analyses were performed on overall survival (OS), transplantation-related mortality (TRM), and disease-related mortality (DRM). Genotype mismatches between recipients and donors in a rare nonsynonymous variant of testis-expressed gene TEX38 significantly increased risk of TRM, which was more dramatic when either the recipient or donor was female. Using the SKAT-O test to evaluate gene-level effects, variant genotypes of OR51D1 in recipients were significantly associated with OS and TRM. In donors, 4 (ALPP, EMID1, SLC44A5, LRP1), 1 (HHAT), and 2 genes (LYZL4, NT5E) were significantly associated with OS, TRM, and DRM, respectively. Inspection of NT5E crystal structures showed 4 of the associated variants affected the enzyme structure and likely decreased the catalytic efficiency of the enzyme. Further confirmation of these findings and additional functional studies may provide individualized risk prediction and prognosis, as well as alternative donor selection strategies.
Subject(s)
Bone Marrow Transplantation/mortality , Exome , Leukemia, Myeloid, Acute/therapy , Myelodysplastic Syndromes/therapy , Precursor Cell Lymphoblastic Leukemia-Lymphoma/therapy , 5'-Nucleotidase/genetics , Adolescent , Adult , Bone Marrow Transplantation/methods , Female , GPI-Linked Proteins/genetics , Genotype , Histocompatibility Testing , Humans , Leukemia, Myeloid, Acute/genetics , Male , Middle Aged , Models, Molecular , Myelodysplastic Syndromes/genetics , Precursor Cell Lymphoblastic Leukemia-Lymphoma/genetics , Transplantation, Homologous/methods , Transplantation, Homologous/mortality , Treatment Outcome , Unrelated Donors , Young AdultABSTRACT
Multiple candidate gene-association studies of non-HLA single-nucleotide polymorphisms (SNPs) and outcomes after blood or marrow transplant (BMT) have been conducted. We identified 70 publications reporting 45 SNPs in 36 genes significantly associated with disease-related mortality, progression-free survival, transplant-related mortality, and/or overall survival after BMT. Replication and validation of these SNP associations were performed using DISCOVeRY-BMT (Determining the Influence of Susceptibility COnveying Variants Related to one-Year mortality after BMT), a well-powered genome-wide association study consisting of 2 cohorts, totaling 2888 BMT recipients with acute myeloid leukemia, acute lymphoblastic leukemia, or myelodysplastic syndrome, and their HLA-matched unrelated donors, reported to the Center for International Blood and Marrow Transplant Research. Gene-based tests were used to assess the aggregate effect of SNPs on outcome. None of the previously reported significant SNPs replicated at P < .05 in DISCOVeRY-BMT. Validation analyses showed association with one previously reported donor SNP at P < .05 and survival; more associations would be anticipated by chance alone. No gene-based tests were significant at P < .05. Functional annotation with publicly available data shows these candidate SNPs most likely do not have biochemical function; only 13% of candidate SNPs correlate with gene expression or are predicted to impact transcription factor binding. Of these, half do not impact the candidate gene of interest; the other half correlate with expression of multiple genes. These findings emphasize the peril of pursing candidate approaches and the importance of adequately powered tests of unbiased genome-wide associations with BMT clinical outcomes given the ultimate goal of improving patient outcomes.
Subject(s)
Bone Marrow Transplantation/mortality , Disease-Free Survival , Genome-Wide Association Study , Polymorphism, Single Nucleotide , Stem Cell Transplantation/mortality , Validation Studies as Topic , Allografts , Humans , Leukemia, Myeloid, Acute/genetics , Leukemia, Myeloid, Acute/mortality , Leukemia, Myeloid, Acute/therapy , Myelodysplastic Syndromes/genetics , Myelodysplastic Syndromes/mortality , Myelodysplastic Syndromes/therapy , Precursor Cell Lymphoblastic Leukemia-Lymphoma/genetics , Precursor Cell Lymphoblastic Leukemia-Lymphoma/mortality , Precursor Cell Lymphoblastic Leukemia-Lymphoma/therapyABSTRACT
BACKGROUND: The pathophysiology of chemotherapy-induced peripheral neuropathy (CIPN) is not well understood. Currently, dose reduction is the only recommendation for alleviating symptoms, often leading to premature treatment cessation. The primary aim of this analysis was to determine the association between components of diet during taxane treatment for breast cancer and change in CIPN symptoms over treatment. METHODS: Women with stage II or III invasive breast cancer were enrolled into an ancillary study to the North American Breast Cancer Intergroup phase III trial (S0221) led by the Southwest Oncology Group (SWOG). Questionnaires including a food frequency questionnaire and the Functional Assessment of Cancer Treatment Gynecologic Oncology Group-Neurotoxicity were administered to assess diet and neuropathic conditions at baseline and during chemotherapy. Ordinal regression was used to estimate odds ratios (ORs) for associations between various food groups and change in neuropathy score (< 10%, 10-30%, > 30%) (n = 900). RESULTS: The odds of worse neuropathy decreased by 21% for each increase in tertile of grain consumption (OR = 0.79, 95% CI 0.66-0.94, p = 0.009). We also observed a nominal 19% increase with higher consumption of citrus fruits (OR = 1.19, 95% CI 1.01-1.40, p = 0.05). CONCLUSIONS: Distinguishing between those who experienced a moderate and a severe change in neuropathy, we found that citrus fruit and grain consumption may play a role in the severity of symptoms. Since there are no existing dietary recommendations for the management of CIPN, further research is needed to investigate whether there may be certain foods that could worsen or alleviate neuropathy symptoms associated with treatment for breast cancer. TRIAL REGISTRATION: ClinicalTrials.gov, NCT03413761 . Registered retrospectively on 29 January 2018.
Subject(s)
Antineoplastic Agents/adverse effects , Breast Neoplasms/drug therapy , Bridged-Ring Compounds/adverse effects , Diet/statistics & numerical data , Peripheral Nervous System Diseases/epidemiology , Taxoids/adverse effects , Adult , Female , Humans , Middle Aged , Peripheral Nervous System Diseases/chemically induced , Peripheral Nervous System Diseases/diet therapy , Peripheral Nervous System Diseases/prevention & control , Quality of Life , Self Report/statistics & numerical dataABSTRACT
BACKGROUND: MicroRNAs (miRNAs) regulate gene expression and influence cancer. Primary transcripts of miRNAs (pri-miRNAs) are poorly annotated and little is known about the role of germline variation in miRNA genes and breast cancer (BC). We sought to identify germline miRNA variants associated with BC risk and tumor subtype among African-American (AA) women. METHODS: Under the African American Breast Cancer Epidemiology and Risk (AMBER) Consortium, genotyping and imputed data from four studies on BC in AA women were combined into a final dataset containing 224,188 miRNA gene single nucleotide polymorphisms (SNPs) for 8350 women: 3663 cases and 4687 controls. The primary miRNA sequence was identified for 566 miRNA genes expressed in Encyclopedia of DNA Elements (ENCODE) Tier 1 cell types and human pancreatic islets. Association analysis was conducted using logistic regression for BC status overall and by tumor subtype. RESULTS: A novel BC signal was localized to an 8.6-kb region of 17q25.3 by four SNPs (rs9913477, rs1428882938, rs28585511, and rs7502931) and remained statistically significant after multiple test correction (odds ratio (OR) = 1.44, 95% confidence interval (CI) = 1.26-1.65; p = 3.15 × 10-7; false discovery rate (FDR) = 0.03). These SNPs reside in a genomic location that includes both the predicted primary transcript of the noncoding miRNA gene MIR3065 and the first intron of the gene for brain-specific angiogenesis inhibitor 1-associated protein 2 (BAIAP2). Furthermore, miRNA-associated SNPs on chromosomes 1p32.3, 5q32, and 3p25.1 were the strongest signals for hormone receptor, luminal versus basal-like, and HER2 enrichment status, respectively. A second phase of genotyping (1397 BC cases, 2418 controls) that included two SNPs in the 8.6-kb region was used for validation and meta-analysis. While neither rs4969239 nor rs9913477 was validated, when meta-analyzed with the original dataset their association with BC remained directionally consistent (OR = 1.29, 95% CI = 1.16-1.44 (p = 4.18 × 10-6) and OR = 1.33, 95% CI = 1.17-1.51 (p = 1.6 × 10-5), respectively). CONCLUSION: Germline genetic variation indicates that MIR3065 may play an important role in BC development and heterogeneity among AA women. Further investigation to determine the potential functional effects of these SNPs is warranted. This study contributes to our understanding of BC risk in AA women and highlights the complexity in evaluating variation in gene-dense regions of the human genome.
Subject(s)
Breast Neoplasms/genetics , Genetic Association Studies , Genetic Predisposition to Disease , MicroRNAs/genetics , Black or African American/genetics , Aged , Female , Gene Expression Regulation, Neoplastic , Genome, Human/genetics , Germ-Line Mutation , Humans , Middle Aged , Nerve Tissue Proteins/genetics , Polymorphism, Single Nucleotide/genetics , Receptor, ErbB-2/geneticsABSTRACT
Multiple breast cancer loci have been identified in previous genome-wide association studies, but they were mainly conducted in populations of European ancestry. Women of African ancestry are more likely to have young-onset and oestrogen receptor (ER) negative breast cancer for reasons that are unknown and understudied. To identify genetic risk factors for breast cancer in women of African descent, we conducted a meta-analysis of two genome-wide association studies of breast cancer; one study consists of 1,657 cases and 2,029 controls genotyped with Illumina's HumanOmni2.5 BeadChip and the other study included 3,016 cases and 2,745 controls genotyped using Illumina Human1M-Duo BeadChip. The top 18,376 single nucleotide polymorphisms (SNP) from the meta-analysis were replicated in the third study that consists of 1,984 African Americans cases and 2,939 controls. We found that SNP rs13074711, 26.5 Kb upstream of TNFSF10 at 3q26.21, was significantly associated with risk of oestrogen receptor (ER)-negative breast cancer (odds ratio [OR]=1.29, 95% CI: 1.18-1.40; P = 1.8 × 10 − 8). Functional annotations suggest that the TNFSF10 gene may be involved in breast cancer aetiology, but further functional experiments are needed. In addition, we confirmed SNP rs10069690 was the best indicator for ER-negative breast cancer at 5p15.33 (OR = 1.30; P = 2.4 × 10 − 10) and identified rs12998806 as the best indicator for ER-positive breast cancer at 2q35 (OR = 1.34; P = 2.2 × 10 − 8) for women of African ancestry. These findings demonstrated additional susceptibility alleles for breast cancer can be revealed in diverse populations and have important public health implications in building race/ethnicity-specific risk prediction model for breast cancer.
Subject(s)
Breast Neoplasms/genetics , Chromosomes, Human, Pair 3/genetics , Black or African American/genetics , Alleles , Black People/genetics , Case-Control Studies , Female , Gene Frequency/genetics , Genetic Loci , Genetic Predisposition to Disease/genetics , Genome-Wide Association Study/methods , Humans , Polymorphism, Single Nucleotide/genetics , Receptors, Estrogen/genetics , Risk Factors , TNF-Related Apoptosis-Inducing Ligand/genetics , TNF-Related Apoptosis-Inducing Ligand/metabolismABSTRACT
OBJECTIVE: Taxane containing chemotherapy extends survival for breast cancer patients. However, taxane-induced peripheral neuropathy (TIPN) cannot be predicted, prevented or effectively treated. Using genome-wide analyses, we sought to identify common risk variants for TIPN. PATIENTS AND METHODS: Women with high-risk breast cancer enrolled in SWOG 0221 were genotyped using the Illumina 1M chip. Genome-wide analyses were performed in relation to ≥grade 3 Common Terminology Criteria for Adverse Events (CTCAE) neuropathy in European and African Americans. Data were meta-analyzed with GW associations of CTCAE ≥grade 3 versus Subject(s)
Breast Neoplasms/drug therapy
, Bridged-Ring Compounds/adverse effects
, Genetic Predisposition to Disease
, Peripheral Nervous System Diseases/genetics
, Taxoids/adverse effects
, Black or African American/genetics
, Breast Neoplasms/complications
, Breast Neoplasms/genetics
, Breast Neoplasms/pathology
, Bridged-Ring Compounds/therapeutic use
, Female
, Genome-Wide Association Study
, Genomics
, Genotype
, Humans
, Peripheral Nervous System Diseases/pathology
, Polymorphism, Single Nucleotide
, Taxoids/therapeutic use
, White People/genetics
ABSTRACT
Previous studies have identified single-nucleotide polymorphisms (SNPs) associated with the risk of chronic graft-versus-host disease (GVHD) after allogeneic hematopoietic cell transplantation. The current study determined whether these associations could be replicated in large cohorts of donors and recipients. Each SNP was tested with cohorts of patients having the same donor type (HLA-matched related, unrelated, or both) reported in the original publication, and testing was limited to the same genome (recipient or donor) and genetic model (dominant, recessive, or allelic) reported in the original study. The 21 SNPs reported in this study represent 19 genes, and the analysis encompassed 22 SNP association tests. The hazard ratio (HR) point estimates and risk ratio point estimates corresponding to odds ratios in previous studies consistently fall outside the 95% confidence intervals of HR estimates in the current study. Despite the large size of the cohorts available for the current study, the 95% confidence intervals for most HRs did not exclude 1.0. Three SNPs representing CTLA4, HPSE, and IL1R1 showed evidence of association with the risk of chronic GVHD in unrelated donor-recipient pairs from 1 cohort, but none of these associations was replicated when tested in unrelated donor-recipient pairs from an independent cohort. Two SNPs representing CCR6 and FGFR1OP showed possible associations with the risk of chronic GVHD in related donor-recipient pairs but not in unrelated donor-recipient pairs. These results remain to be tested for replication in other cohorts of related donor-recipient pairs.
Subject(s)
Graft vs Host Disease/genetics , Polymorphism, Single Nucleotide , Adolescent , Adult , Aged , Child , Child, Preschool , Chronic Disease , Cohort Studies , Female , Genetic Association Studies , Hematopoietic Stem Cell Transplantation , Histocompatibility Testing , Humans , Infant , Infant, Newborn , Male , Middle Aged , Transplantation, Homologous , Unrelated Donors , Young AdultABSTRACT
BACKGROUND: The nuclear hormone receptor superfamily acts as a genomic sensor of diverse signals. Their actions are often intertwined with other transcription factors. Nuclear hormone receptors are targets for many therapeutic drugs, and include the vitamin D receptor (VDR). VDR signaling is pleotropic, being implicated in calcaemic function, antibacterial actions, growth control, immunomodulation and anti-cancer actions. Specifically, we hypothesized that the biologically significant relationships between the VDR transcriptome and phenotype-associated biology could be discovered by integrating the known VDR transcription factor binding sites and all published trait- and disease-associated SNPs. By integrating VDR genome-wide binding data (ChIP-seq) with the National Human Genome Research Institute (NHGRI) GWAS catalog of SNPs we would see where and which target gene interactions and pathways are impacted by inherited genetic variation in VDR binding sites, indicating which of VDR's multiple functions are most biologically significant. RESULTS: To examine how genetic variation impacts VDR function we overlapped 23,409 VDR genomic binding peaks from six VDR ChIP-seq datasets with 191,482 SNPs, derived from GWAS-significant SNPs (Lead SNPs) and their correlated variants (r 2 > 0.8) from HapMap3 and the 1000 genomes project. In total, 574 SNPs (71 Lead and 503 SNPs in linkage disequilibrium with Lead SNPs) were present at VDR binding loci and associated with 211 phenotypes. For each phenotype a hypergeometric test was used to determine if SNPs were enriched at VDR binding sites. Bonferroni correction for multiple testing across the 211 phenotypes yielded 42 SNPs that were either disease- or phenotype-associated with seven predominately immune related including self-reported allergy; esophageal cancer was the only cancer phenotype. Motif analyses revealed that only two of these 42 SNPs reside within a canonical VDR binding site (DR3 motif), and that 1/3 of the 42 SNPs significantly impacted binding and gene regulation by other transcription factors, including NF-κB. This suggests a plausible link for the potential cross-talk between VDR and NF-κB. CONCLUSIONS: These analyses showed that VDR peaks are enriched for SNPs associated with immune phenotypes suggesting that VDR immunomodulatory functions are amongst its most important actions. The enrichment of genetic variation in non-DR3 motifs suggests a significant role for the VDR to bind in multimeric complexes containing other transcription factors that are the primary DNA binding component. Our work provides a framework for the combination of ChIP-seq and GWAS findings to provide insight into the underlying phenotype-associated biology of a given transcription factor.
Subject(s)
Genome-Wide Association Study , Immunity/genetics , NF-kappa B/metabolism , Phenotype , Polymorphism, Single Nucleotide , Receptors, Calcitriol/genetics , Receptors, Calcitriol/metabolism , Cell Line , Genomics , Humans , Linkage Disequilibrium , Protein Binding , Transcription Factors/metabolismABSTRACT
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.
Subject(s)
Black or African American/genetics , Breast Neoplasms/ethnology , DNA Methylation , Hepatocyte Nuclear Factor 3-alpha/genetics , Hepatocyte Nuclear Factor 3-alpha/metabolism , Parity/genetics , Breast Feeding , Breast Neoplasms/genetics , Breast Neoplasms/metabolism , Down-Regulation , Female , Genetic Association Studies , High-Throughput Nucleotide Sequencing , Humans , Linear Models , Receptors, Estrogen/metabolism , Sequence Analysis, DNA , Sequence Analysis, RNA , White People/geneticsABSTRACT
Single nucleotide polymorphisms (SNPs) in pathways influencing lymph node (LN) metastasis and estrogen receptor (ER) status in breast cancer may partially explain inter-patient variability in prognosis. We examined 154 SNPs in 12 metastasis-related genes for associations with breast cancer risk, stratified by LN and ER status, in European-American (EA) and African-American (AA) women. Two-thousand six hundred and seventy-one women enrolled in the Women's Circle of Health Study were genotyped. Pathway analyses were conducted using the adaptive rank truncated product (ARTP) method, with pARTP ≤ 0.10 as significant. Multi-allelic risk scores were created for the ARTP-significant gene(s). Single-SNP and risk score associations were modeled using logistic regression, with false discovery rate (FDR) P-value adjustment. Although single-SNP associations were not significant at pFDR < 0.05, several genes were significant in the ARTP analyses. In AA women, significant ARTP gene-level associations included CDH1 with LN+ (pARTP = 0.10; multi-allelic OR = 1.13, 95%CI 1.07-1.19, pFDR = 0.0003) and SIPA1 with ER- breast cancer (pARTP = 0.10; multi-allelic OR = 1.16, 95%CI 1.02-1.31, pFDR = 0.03). In EA women, MTA2 was associated with overall breast cancer risk (pARTP = 0.004), regardless of ER status, and with LN- disease (pARTP = 0.01). Also significant were SATB1 in ER- (pARTP = 0.03; multi-allelic OR = 1.12, 95%CI 1.05-1.20, pFDR = 0.003) and KISS1 in LN- (pARTP = 0.10; multi-allelic OR = 1.18, 95%CI 1.08-1.29, pFDR = 0.002) analyses. Among LN+ cases, significant ARTP associations were observed for SNAI1, CD82, NME1, and CTNNB1 (multi-allelic OR = 1.09, 95%CI 1.04-1.14, pFDR = 0.001). Our findings suggest that variants in several metastasis genes may affect breast cancer risk by LN or ER status, although verification in larger studies is required. © 2016 Wiley Periodicals, Inc.
Subject(s)
Black People/genetics , Breast Neoplasms/genetics , Lymph Nodes/pathology , Polymorphism, Single Nucleotide , Receptors, Estrogen/metabolism , White People/genetics , Adult , Aged , Breast Neoplasms/metabolism , Female , Genetic Predisposition to Disease , Humans , Middle Aged , Neoplasm Metastasis , Young AdultABSTRACT
OBJECTIVE: There is a mounting body of evidence demonstrating higher percentages of regulatory T (Treg) cells in the peripheral blood of patients with cancer in comparison to healthy controls, but there is a paucity of epidemiological literature characterizing circulating Treg cells among patients with epithelial ovarian cancer (EOC). To investigate the role of peripheral Treg cells in ovarian neoplasms, we conducted a case-control study to characterize circulating concentrations of Treg cells among patients with EOC, women with benign ovarian conditions, and healthy controls without a history of cancer. MATERIALS AND METHODS: Participants were identified for inclusion due to their participation in the Data Bank and BioRepository program at Roswell Park Cancer Institute in Buffalo, NY. Patients included 71 women with a primary diagnosis of EOC and 195 women with a diagnosis of benign ovarian conditions. Controls included 101 age- and race-matched women without a history of cancer. Nonfasting, pretreatment peripheral blood levels of CD3+CD4+CD25+FOXP3+ Treg cells were measured using flow cytometric analyses and expressed as a percentage of total CD3+ cells and as a percentage of total CD3+CD4+ cells. RESULTS: Compared to healthy controls and women with benign ovarian conditions, patients with EOC had significantly higher frequency of Treg cells (P < 0.04). In multivariable logistic regression analyses using Treg frequency expressed as a percentage of CD+3 cells, we observed a significant positive association between Treg cell percentage and EOC risk, with each 1% increase associated with a 37% increased risk of EOC (odds ratio, 1.37; 95% confidence interval, 1.04-1.80). We observed a similar trend when Treg frequency was expressed as a percentage of CD3+CD+4 cells (odds ratio, 1.22; 95% confidence interval, 0.99-1.49). CONCLUSIONS: The current study provides support that peripheral Treg cell frequency is elevated in patients with EOC in comparison to women with benign ovarian conditions and healthy controls.
Subject(s)
Neoplasms, Glandular and Epithelial/blood , Ovarian Neoplasms/blood , T-Lymphocytes, Regulatory/pathology , Age Factors , Carcinoma, Ovarian Epithelial , Case-Control Studies , Female , Humans , Middle Aged , Neoplasms, Glandular and Epithelial/immunology , Ovarian Neoplasms/immunology , T-Lymphocytes, Regulatory/immunologyABSTRACT
A large percentage of breast cancer heritability remains unaccounted for, and most of the known susceptibility loci have been established in European and Asian populations. Rare variants may contribute to the unexplained heritability of this disease, including in women of African ancestry (AA). We conducted an exome-wide analysis of rare variants in relation to risk of overall and subtype-specific breast cancer in the African American Breast Cancer Epidemiology and Risk (AMBER) Consortium, which includes data from four large studies of AA women. Genotyping on the Illumina Human Exome Beadchip yielded data for 170 812 SNPs and 8287 subjects: 3629 cases (1093 estrogen receptor negative (ER-), 1968 ER+, 568 ER unknown) and 4658 controls, the largest exome chip study to date for AA breast cancer. Pooled gene-based association analyses were performed using the unified optimal sequence kernel association test (SKAT-O) for variants with minor allele frequency (MAF) ≤ 5%. In addition, each variant with MAF >0.5% was tested for association using logistic regression. There were no significant associations with overall breast cancer. However, a novel gene, FBXL22 (P = 8.2×10(-6)), and a gene previously identified in GWAS of European ancestry populations, PDE4D (P = 1.2×10(-6)), were significantly associated with ER- breast cancer after correction for multiple testing. Cases with the associated rare variants were also negative for progesterone and human epidermal growth factor receptors-thus, triple-negative cancer. Replication is required to confirm these gene-level associations, which are based on very small counts at extremely rare SNPs.
Subject(s)
Black or African American/genetics , Breast Neoplasms/genetics , Adult , Aged , Breast Neoplasms/chemistry , Breast Neoplasms/etiology , Cyclic Nucleotide Phosphodiesterases, Type 4/genetics , Exome , F-Box Proteins/genetics , Female , Genome-Wide Association Study , Humans , Logistic Models , Middle Aged , Polymorphism, Single Nucleotide , Receptors, Estrogen/analysis , RiskABSTRACT
The phosphatidylinositol 3-kinase-AKT-mammalian target of rapamycin (mTOR) pathway has been implicated in breast carcinogenesis. However, there has been no large-scale investigation of genetic variants in the mTOR pathway and breast cancer risk. We examined 28847 single-nucleotide polymorphisms (SNPs) in 61 mTOR pathway genes in the African American Breast Cancer Epidemiology and Risk consortium of 3663 cases [1983 estrogen receptor-positive (ER+) and 1098 ER-negative (ER-)] and 4687 controls. Gene-level analyses were conducted using the adaptive rank truncated product (ARTP) test for 10773 SNPs that were not highly correlated (r (2) < 0.8), and SNP-level analyses were conducted with logistic regression. Among genes that were prioritized (nominal P < 0.05, ARTP tests), associations were observed for intronic SNPs TSC2 rs181088346 [odds ratio (OR) of each copy of variant allele = 0.77, 95% confidence interval (CI) = 0.65-0.88 for all breast cancer] and BRAF rs114729114 (OR = 1.53, 95% CI = 1.24-1.91 for all breast cancer and OR = 2.03, 95% CI = 1.50-2.76 for ER- tumors). For ER- tumors, intronic SNPs PGF rs11542848 (OR = 1.38, 95% CI = 1.15-1.66) and rs61759375 (OR = 1.34, 95% CI = 1.14-1.57) and MAPK3 rs78564187 (OR = 1.26, 95% CI = 1.11-1.43) were associated with increased risk. These SNPs were significant at a gene-wide level (Bonferroni-corrected P < 0.05). The variant allele of RPS6KB2 rs35363135, a synonymous coding SNP, was more likely to be observed in ER- than ER+ tumors (OR = 1.18, 95% CI = 1.05-1.31, gene-wide Bonferroni-corrected P = 0.06). In conclusion, specific mTOR pathway genes are potentially important to breast cancer risk and to the ER negativity in African American women.
Subject(s)
Black or African American/genetics , Breast Neoplasms/ethnology , Breast Neoplasms/genetics , TOR Serine-Threonine Kinases/genetics , Adult , Aged , Breast Neoplasms/metabolism , Carcinoma in Situ/ethnology , Carcinoma in Situ/genetics , Carcinoma in Situ/metabolism , Carcinoma, Ductal, Breast/ethnology , Carcinoma, Ductal, Breast/genetics , Carcinoma, Ductal, Breast/metabolism , Case-Control Studies , Female , Genetic Predisposition to Disease , Genetic Variation , Humans , Middle Aged , Neoplasm Invasiveness , Polymorphism, Single Nucleotide , Receptors, Estrogen/metabolism , Signal Transduction/genetics , TOR Serine-Threonine Kinases/metabolism , Young AdultABSTRACT
The Hippo signaling pathway regulates cellular proliferation and survival, thus exerting profound effects on normal cell fate and tumorigenesis. Dysfunction of the Hippo pathway components has been linked with breast cancer stem cell regulation, as well as breast tumor progression and metastasis. TAZ, a key component of the Hippo pathway, is highly expressed in triple negative breast cancer; however, the associations of genetic variations in this important pathway with breast cancer risk remain largely unexplored. Here, we analyzed 8309 germline variants in 15 genes from the Hippo pathway with a total of 3663 cases and 4687 controls from the African American Breast Cancer Epidemiology and Risk Consortium. Odds ratios (ORs) were estimated using logistic regression for overall breast cancer, by estrogen receptor (ER) status (1983 ER positive and 1098 ER negative), and for case-only analyses by ER status. The Hippo signaling pathway was significantly associated with ER-negative breast cancer (pathway level P = 0.02). Gene-based analyses revealed that CDH1 was responsible for the pathway association (P < 0.01), with rs4783673 in CDH1 statistically significant after gene-level adjustment for multiple comparisons (P = 9.2×10(-5), corrected P = 0.02). rs142697907 in PTPN14 was associated with ER-positive breast cancer and rs2456773 in CDK1 with ER-negativity in case-only analysis after gene-level correction for multiple comparisons (corrected P < 0.05). In conclusion, common genetic variations in the Hippo signaling pathway may contribute to both ER-negative and ER+ breast cancer risk in AA women.
Subject(s)
CDC2 Protein Kinase/genetics , Cadherins/genetics , Estrogen Receptor alpha/genetics , Protein Tyrosine Phosphatases, Non-Receptor/genetics , Triple Negative Breast Neoplasms/genetics , Antigens, CD , Female , Gene Expression Regulation, Neoplastic , Genetic Association Studies , Genetic Predisposition to Disease , Hippo Signaling Pathway , Humans , Neoplasm Metastasis , Neoplasm Proteins/genetics , Neoplastic Stem Cells/metabolism , Neoplastic Stem Cells/pathology , Polymorphism, Single Nucleotide , Protein Serine-Threonine Kinases/genetics , Risk Factors , Signal Transduction , Triple Negative Breast Neoplasms/pathologyABSTRACT
Studies of genetic variations in vitamin D-related pathways and breast cancer risk have been conducted mostly in populations of European ancestry, and only sparsely in African Americans (AA), who are known for a high prevalence of vitamin D deficiency. We analyzed 24,445 germline variants in 63 genes from vitamin D-related pathways in the African American Breast Cancer Epidemiology and Risk (AMBER) consortium, including 3,663 breast cancer cases and 4,687 controls. Odds ratios (OR) were derived from logistic regression models for overall breast cancer, by estrogen receptor (ER) status (1,983 ER positive and 1,098 ER negative), and for case-only analyses of ER status. None of the three vitamin D-related pathways were associated with breast cancer risk overall or by ER status. Gene-level analyses identified associations with risk for several genes at a nominal p ≤ 0.05, particularly for ER- breast cancer, including rs4647707 in DDB2. In case-only analyses, vitamin D metabolism and signaling pathways were associated with ER- cancer (pathway-level p = 0.02), driven by a single gene CASR (gene-level p = 0.001). The top SNP in CASR was rs112594756 (p = 7 × 10(-5), gene-wide corrected p = 0.01), followed by a second signal from a nearby SNP rs6799828 (p = 1 × 10(-4), corrected p = 0.03). In summary, several variants in vitamin D pathways were associated with breast cancer risk in AA women. In addition, CASR may be related to tumor ER status, supporting a role of vitamin D or calcium in modifying breast cancer phenotypes.