Polygenic risk scores for prediction of breast cancer risk in women of African ancestry: a cross-ancestry approach.

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.

Associations of minor histocompatibility antigens with outcomes following allogeneic hematopoietic cell transplantation.

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.

Genetic ancestry and population differences in levels of inflammatory cytokines in women: Role for evolutionary selection and environmental factors.

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.

gwasurvivr: an R package for genome-wide survival analysis.

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.

Exome chip analyses identify genes affecting mortality after HLA-matched unrelated-donor blood and marrow transplantation.

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.

Replication and validation of genetic polymorphisms associated with survival after allogeneic blood or marrow transplant.

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.

Associations between self-reported diet during treatment and chemotherapy-induced peripheral neuropathy in a cooperative group trial (S0221).

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.

A survey of microRNA single nucleotide polymorphisms identifies novel breast cancer susceptibility loci in a case-control, population-based study of African-American women.

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.

Genome-wide association studies in women of African ancestry identified 3q26.21 as a novel susceptibility locus for oestrogen receptor negative breast cancer.

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.

Genome-wide meta-analyses identifies novel taxane-induced peripheral neuropathy-associated loci.

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

Integration of VDR genome wide binding and GWAS genetic variation data reveals co-occurrence of VDR and NF-κB binding that is linked to immune phenotypes.

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.

Single nucleotide variants in metastasis-related genes are associated with breast cancer risk, by lymph node involvement and estrogen receptor status, in women with European and African ancestry.

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.

The Association of Peripheral Blood Regulatory T-Cell Concentrations With Epithelial Ovarian Cancer: A Brief Report.

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.

An exome-wide analysis of low frequency and rare variants in relation to risk of breast cancer in African American Women: the AMBER Consortium.

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.

Gene-based analysis of the fibroblast growth factor receptor signaling pathway in relation to breast cancer in African American women: the AMBER consortium.

We conducted gene-based analysis in 26 genes in the FGFR signaling pathway to identify genes carrying genetic variation affecting risk of breast cancer and the specific estrogen receptor (ER) subtypes. Tagging single-nucleotide polymorphisms (SNPs) for each gene were selected and genotyped on a customized Illumina Exome Array. Imputation was carried out using 1000 Genomes haplotypes. The analysis included 3237 SNPs in 3663 breast cancer cases (including 1983 ER-positive, and 1098 ER-negative) and 4687 controls from the African American Breast Cancer Epidemiology and Risk consortium, a collaborative project of four large studies of breast cancer in African American women (Carolina Breast Cancer Study, Black Women's Health Study, Women's Circle of Health Study, and Multiethnic Cohort). We used a multi-locus adaptive joint (AdaJoint) test to determine the association of each gene in the FGFR signaling pathway with overall breast cancer and ER subtypes. The FGF1 gene was significantly associated with risk of ER-negative breast cancer (P = 0.001). The FGFR2 gene was associated with risk of overall breast cancer (P = 0.002) and ER-positive breast cancer (P = 0.002). The FGF1 gene affects risk of ER-negative breast cancer in African American women. We confirmed the association of the FGFR2 gene with risk of overall and ER-positive breast cancer. These results highlight the importance of the FGFR signaling pathway in the pathogenesis of breast cancer, and suggest that different genes in the same pathway may be associated with different ER breast cancer subtypes.

Vitamin D receptor and RXR in the post-genomic era.

Following the elucidation of the human genome and components of the epigenome, it is timely to revisit what is known of vitamin D receptor (VDR) function. Early transcriptomic studies using microarray approaches focused on the protein coding mRNA that were regulated by the VDR, usually following treatment with ligand. These studies quickly established the approximate size and surprising diversity of the VDR transcriptome, revealing it to be highly heterogenous and cell type and time dependent. Investigators also considered VDR regulation of non-protein coding RNA and again, cell and time dependency was observed. Attempts to integrate mRNA and miRNA regulation patterns are beginning to reveal patterns of co-regulation and interaction that allow for greater control of mRNA expression, and the capacity to govern more complex cellular events. Alternative splicing in the trasncriptome has emerged as a critical process in transcriptional control and there is evidence of the VDR interacting with components of the splicesome. ChIP-Seq approaches have proved to be pivotal to reveal the diversity of the VDR binding choices across cell types and following treatment, and have revealed that the majority of these are non-canonical in nature. The underlying causes driving the diversity of VDR binding choices remain enigmatic. Finally, genetic variation has emerged as important to impact the transcription factor affinity towards genomic binding sites, and recently the impact of this on VDR function has begun to be considered.

Establishment of Definitions and Review Process for Consistent Adjudication of Cause-specific Mortality after Allogeneic Unrelated-donor Hematopoietic Cell Transplantation.

Clinical trials commonly use adjudication committees to refine endpoints, but observational research or genome-wide association studies rarely do. Our goals were to establish definitions of cause-specific death after unrelated-donor allogeneic hematopoietic cell transplantation (URD-HCT), to estimate discordance between reported and adjudicated cause-specific death, and to identify factors contributing to inconsistency in cause-specific death determination. A consensus panel adjudicated cause-specific death in 1484 patients who died within 1 year after HCT, derived from 3532 acute leukemia or myelodysplasia patients after URD-HCT from 2000 to 2011 reported by 151 US transplant centers to the Center for International Blood and Marrow Transplant Research. Deaths were classified as disease-related or transplant-related. The panel agreed with >99% of deaths reported by centers as disease-related and 80% reported as transplant-related. Year of transplant (cohort effect) and disease status significantly influenced agreement between the panel and centers. Sensitivity analysis of deaths < 100 days post-transplant yielded the lowest agreement between the panel and centers for myelodysplastic syndrome patients. Standard predefined criteria for adjudicating cause-specific death led to consistent application to similar clinical scenarios and clearer delineation of cause-specific death categories. Other studies of competing events such as cancer-specific versus treatment-related mortality would benefit from our results. Our detailed algorithm should result in more consistent reporting of cause-specific death by centers.

Hormone-related pathways and risk of breast cancer subtypes in African American women.

We sought to investigate genetic variation in hormone pathways in relation to risk of overall and subtype-specific breast cancer in women of African ancestry (AA). Genotyping and imputation yielded data on 143,934 SNPs in 308 hormone-related genes for 3663 breast cancer cases (1098 ER-, 1983 ER+, 582 ER unknown) and 4687 controls from the African American Breast Cancer Epidemiology and Risk (AMBER) Consortium. AMBER includes data from four large studies of AA women: the Carolina Breast Cancer Study, the Women's Circle of Health Study, the Black Women's Health Study, and the Multiethnic Cohort Study. Pathway- and gene-based analyses were conducted, and single-SNP tests were run for the top genes. There were no strong associations at the pathway level. The most significantly associated genes were GHRH, CALM2, CETP, and AKR1C1 for overall breast cancer (gene-based nominal p ≤ 0.01); NR0B1, IGF2R, CALM2, CYP1B1, and GRB2 for ER+ breast cancer (p ≤ 0.02); and PGR, MAPK3, MAP3K1, and LHCGR for ER- disease (p ≤ 0.02). Single-SNP tests for SNPs with pairwise linkage disequilibrium r (2) < 0.8 in the top genes identified 12 common SNPs (in CALM2, CETP, NR0B1, IGF2R, CYP1B1, PGR, MAPK3, and MAP3K1) associated with overall or subtype-specific breast cancer after gene-level correction for multiple testing. Rs11571215 in PGR (progesterone receptor) was the SNP most strongly associated with ER- disease. We identified eight genes in hormone pathways that contain common variants associated with breast cancer in AA women after gene-level correction for multiple testing.

Cooperative behavior of the nuclear receptor superfamily and its deregulation in prostate cancer.

The current study aimed to assess the topology of the nuclear receptor (NR) superfamily in normal prostate epithelial cells and its distortion in prostate cancer. Both in vitro and in silico approaches were utilized to profile NRs expressed in non-malignant RWPE-1 cells, which were subsequently investigated by treating cells with 132 binary NR ligand combinations. Nine significant cooperative interactions emerged including both superadditive [22(R)-hydroxycholesterol and eicosatetraenoic acid] and subadditive [1α,25(OH)2D3 and chenodeoxycholic acid] cellular responses, which could be explained in part by cooperative control of cell-cycle progression and candidate gene expression. In addition, publicly available data were employed to assess NR expression in human prostate tissue. Common and significant loss of NR superfamily expression was established in publicly available data from prostate tumors, in part predicting parallel distortion of targeting microRNA. These findings suggest that the NR superfamily in the prostate cooperatively integrates signals from dietary, hormonal and metabolic cues, and is significantly distorted in prostate cancer.

Epigenetic disruption of the RARγ complex impairs its function to bookmark AR enhancer interactions required for enzalutamide sensitivity in prostate cancer.

The current study in prostate cancer (PCa) focused on the genomic mechanisms at the cross-roads of pro-differentiation signals and the emergence of lineage plasticity. We explored an understudied cistromic mechanism involving RARγ's ability to govern AR cistrome-transcriptome relationships, including those associated with more aggressive PCa features. The RARγ complex in PCa cell models was enriched for canonical cofactors, as well as proteins involved in RNA processing and bookmarking. Identifying the repertoire of miR-96 bound and regulated gene targets, including those recognition elements marked by m6A, revealed their significant enrichment in the RARγ complex. RARγ significantly enhanced the AR cistrome, particularly in active enhancers and super-enhancers, and overlapped with the binding of bookmarking factors. Furthermore, RARγ expression led to nucleosome-free chromatin enriched with H3K27ac, and significantly enhanced the AR cistrome in G2/M cells. RARγ functions also antagonized the transcriptional actions of the lineage master regulator ONECUT2. Similarly, gene programs regulated by either miR-96 or antagonized by RARγ were enriched in alternative lineages and more aggressive PCa phenotypes. Together these findings reveal an under-investigated role for RARγ, modulated by miR-96, to bookmark enhancer sites during mitosis. These sites are required by the AR to promote transcriptional competence, and emphasize luminal differentiation, while antagonizing ONECUT2.