Characterizing Genetic Susceptibility to Breast Cancer in Women of African Ancestry

Novel breast cancer susceptibility loci under linkage peaks identified in African ancestry consortia

BACKGROUND

Expansion of genome-wide association studies across population groups is needed to improve our understanding of shared and unique genetic contributions to breast cancer. We performed association and replication studies guided by a priori linkage findings from African ancestry (AA) relative pairs.

METHODS

We performed fixed-effect inverse-variance weighted meta-analysis under three significant AA breast cancer linkage peaks (3q26-27, 12q22-23, and 16q21-22) in 9241 AA cases and 10 193 AA controls. We examined associations with overall breast cancer as well as estrogen receptor (ER)-positive and negative subtypes (193,132 SNPs). We replicated associations in the African-ancestry Breast Cancer Genetic Consortium (AABCG).

RESULTS

In AA women, we identified two associations on chr12q for overall breast cancer (rs1420647, OR = 1.15, p = 2.50×10-6; rs12322371, OR = 1.14, p = 3.15×10-6), and one for ER-negative breast cancer (rs77006600, OR = 1.67, p = 3.51×10-6). On chr3, we identified two associations with ER-negative disease (rs184090918, OR = 3.70, p = 1.23×10-5; rs76959804, OR = 3.57, p = 1.77×10-5) and on chr16q we identified an association with ER-negative disease (rs34147411, OR = 1.62, p = 8.82×10-6). In the replication study, the chr3 associations were significant and effect sizes were larger (rs184090918, OR: 6.66, 95% CI: 1.43, 31.01; rs76959804, OR: 5.24, 95% CI: 1.70, 16.16).

CONCLUSION

The two chr3 SNPs are upstream to open chromatin ENSR00000710716, a regulatory feature that is actively regulated in mammary tissues, providing evidence that variants in this chr3 region may have a regulatory role in our target organ. Our study provides support for breast cancer variant discovery using prioritization based on linkage evidence.

Single nucleotide polymorphisms associated with female breast cancer susceptibility in Chinese population

Breast cancer is a complex disease influenced by both external and internal factors, among which genetic factors play a critical role. Single-nucleotide polymorphisms (SNPs) are major contributors to the heritability of breast cancer, and their frequencies vary across ethnic groups. In this study, we aimed to investigate the association between 34 SNPs identified in previous genome-wide association studies (GWAS) and overall breast cancer risk, as well as breast cancer subtypes, in the Chinese female population. To accomplish this, we conducted an extensive association analysis using the high-throughput Sequenom MassARRAY® platform in a case-control study comprising 1848 breast cancer patients and 709 healthy controls. Our analysis, which utilized the SNPassoc package in R based on chi-squared (χ2) test and genetic model analysis, identified significant associations between breast cancer risk and SNP rs12493607 (TGFBR2, risk allele C, OR = 1.28 [1.11-1.47], P = 0.0005), as well as a less conservatively significant association with rs4784227 (CASC16, risk allele T, OR = 1.24 [1.08-1.42], P = 0.0017) and rs2046210 (ESR1, risk allele A, OR = 1.50 [1.16-1.95], P = 0.0016). Furthermore, our stratified analyses revealed that rs12493607 was significantly associated with invasive carcinoma, estrogen receptor (ER)-positive, progesterone receptor (PR)-positive, HER2-negative, and young (aged younger than 45) breast cancer. SNP rs4784227 and rs3803662 (CASC16) were associated with invasive carcinoma and ER-positive breast cancer, while rs2046210 was linked to ductal carcinoma in situ, ER-negative, PR-negative, HER2-positive, and elder (aged more than 45) breast cancers. SNPs rs10484919 (ESR1) and rs1038304 (CCDC170) showed links to HER2-positive breast cancer, and rs616488 (PEX14) with premenopausal breast cancer. In summary, our study shed light on the relationship between SNPs and breast cancer susceptibility within a vast Chinese cohort, supporting the development of polygenetic risk scores for the Chinese population. These findings provide valuable insights into the genetic basis of breast cancer and have important implications for risk prediction, early detection, and personalized treatment of this disease.

m6A readers, writers, erasers, and the m6A epitranscriptome in breast cancer

Epitranscriptomic modification of RNA regulates human development, health, and disease. The true diversity of the transcriptome in breast cancer including chemical modification of transcribed RNA (epitranscriptomics) is not well understood due to limitations of technology and bioinformatic analysis. N-6-methyladenosine (m6A) is the most abundant epitranscriptomic modification of mRNA and regulates splicing, stability, translation, and intracellular localization of transcripts depending on m6A association with reader RNA-binding proteins. m6A methylation is catalyzed by the METTL3 complex and removed by specific m6A demethylase ALKBH5, with the role of FTO as an 'eraser' uncertain. In this review, we provide an overview of epitranscriptomics related to mRNA and focus on m6A in mRNA and its detection. We summarize current knowledge on altered levels of writers, readers, and erasers of m6A and their roles in breast cancer and their association with prognosis. We summarize studies identifying m6A peaks and sites in genes in breast cancer cells.

Development and Validation of a Breast Cancer Polygenic Risk Score on the Basis of Genetic Ancestry Composition

PURPOSE

Polygenic risk scores (PRSs) for breast cancer (BC) risk stratification have been developed primarily in women of European ancestry. Their application to women of non-European ancestry has lagged because of the lack of a formal approach to incorporate genetic ancestry and ancestry-dependent variant frequencies and effect sizes. Here, we propose a multiple-ancestry PRS (MA-PRS) that addresses these issues and may be useful in the development of equitable PRSs across other cancers and common diseases.

MATERIALS AND METHODS

Women referred for hereditary cancer testing were divided into consecutive cohorts for development (n = 189,230) and for independent validation (n = 89,126). Individual genetic composition as fractions of three reference ancestries (African, East Asian, and European) was determined from ancestry-informative single-nucleotide polymorphisms. The MA-PRS is a combination of three ancestry-specific PRSs on the basis of genetic ancestral composition. Stratification of risk was evaluated by multivariable logistic regression models controlling for family cancer history. Goodness-of-fit analysis compared expected with observed relative risks by quantiles of the MA-PRS distribution.

RESULTS

In independent validation, the MA-PRS was significantly associated with BC risk in the full cohort (odds ratio, 1.43; 95% CI, 1.40 to 1.46; P = 8.6 × 10-308) and within each major ancestry. The top decile of the MA-PRS consistently identified patients with two-fold increased risk of developing BC. Goodness-of-fit tests showed that the MA-PRS was well calibrated and predicted BC risk accurately in the tails of the distribution for both European and non-European women.

CONCLUSION

The MA-PRS uses genetic ancestral composition to expand the utility of polygenic risk prediction to non-European women. Inclusion of genetic ancestry in polygenic risk prediction presents an opportunity for more personalized treatment decisions for women of varying and mixed ancestries.

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.

Exploring the effects of genetic variation on gene regulation in cancer in the context of 3D genome structure

Background

Numerous genome-wide association studies (GWAS) conducted to date revealed genetic variants associated with various diseases, including breast and prostate cancers. Despite the availability of these large-scale data, relatively few variants have been functionally characterized, mainly because the majority of single-nucleotide polymorphisms (SNPs) map to the non-coding regions of the human genome. The functional characterization of these non-coding variants and the identification of their target genes remain challenging.

Results

In this communication, we explore the potential functional mechanisms of non-coding SNPs by integrating GWAS with the high-resolution chromosome conformation capture (Hi-C) data for breast and prostate cancers. We show that more genetic variants map to regulatory elements through the 3D genome structure than the 1D linear genome lacking physical chromatin interactions. Importantly, the association of enhancers, transcription factors, and their target genes with breast and prostate cancers tends to be higher when these regulatory elements are mapped to high-risk SNPs through spatial interactions compared to simply using a linear proximity. Finally, we demonstrate that topologically associating domains (TADs) carrying high-risk SNPs also contain gene regulatory elements whose association with cancer is generally higher than those belonging to control TADs containing no high-risk variants.

Conclusions

Our results suggest that many SNPs may contribute to the cancer development by affecting the expression of certain tumor-related genes through long-range chromatin interactions with gene regulatory elements. Integrating large-scale genetic datasets with the 3D genome structure offers an attractive and unique approach to systematically investigate the functional mechanisms of genetic variants in disease risk and progression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-01021-x.

Can the epigenome contribute to risk stratification for cancer onset?

The increasing burden of cancer requires identifying and protecting individuals at highest risk. The epigenome provides an indispensable complement to genetic alterations for a risk stratification approach for the following reasons: gene transcription necessary for cancer onset is directed by epigenetic modifications and many risk factors studied so far have been associated with alterations related to the epigenome. The risk level depends on the plasticity of the epigenome during phases of life particularly sensitive to environmental and dietary impacts. Modifications in the activity of DNA regulatory regions and altered chromatin compaction may accumulate, hence leading to the increase of cancer risk. Moreover, tissue architecture directs the unique organization of the epigenome for each tissue and cell type, which allows the epigenome to control cancer risk in specific organs. Investigations of epigenetic signatures of risk should help identify a continuum of alterations leading to a threshold beyond which the epigenome cannot maintain homeostasis. We propose that this threshold may be similar in the population for a given tissue, but the pace to reach this threshold will depend on the combination of germline inheritance and the risk and protective factors encountered, particularly during windows of epigenetic susceptibility, by individuals.

Identifying the link between chemical exposures and breast cancer in African American women via integrated in vitro and exposure biomarker data

Among women, breast cancer is the most prevalent form of cancer worldwide and has the second highest mortality rate of any cancer in the United States. The breast cancer related death rate is 40 % higher in non-Hispanic Black women compared to non-Hispanic White women. The incidence of triple negative breast cancer (TNBC), an aggressive subtype of breast cancer for which there is no targeted therapy, is also approximately three times higher for Black, relative to, White women. The drivers of these differences are poorly understood. Here, we aimed to identify chemical exposures which play a role in breast cancer disparities. Using chemical biomonitoring data from the National Health and Nutrition Examination Survey (NHANES) and biological activity data from the EPA's ToxCast program, we assessed the toxicological profiles of chemicals to which US Black women are disproportionately exposed. We conducted a literature search to identify breast cancer targets in ToxCast to analyze the response of chemicals with exposure disparities in these assays. Forty-three chemical biomarkers are significantly higher in Black women. Investigation of these chemicals in ToxCast resulted in 32,683 assays for analysis, 5172 of which contained nonzero values for the concentration at which the dose-response fitted model reaches the cutoff considered "active". Of these chemicals BPA, PFOS, and thiram are most comprehensively assayed. 2,5-dichlorophenol, 1,4-dichlorobenzene, and methyl and propyl parabens had higher biomarker concentrations in Black women and moderate testing and activity in ToxCast. The distribution of active concentrations for these chemicals in ToxCast assays are comparable to biomarker concentrations in Black women NHANES participants. Through this integrated analysis, we identify that multiple chemicals, including thiram, propylparaben, and p,p' DDE, have disproportionate exposures in Black women and have breast cancer associated biological activity at human exposure relevant doses.

Evaluating Polygenic Risk Scores for Breast Cancer in Women of African Ancestry

BACKGROUND

Polygenic risk scores (PRSs) have been demonstrated to identify women of European, Asian, and Latino ancestry at elevated risk of developing breast cancer (BC). We evaluated the performance of existing PRSs trained in European ancestry populations among women of African ancestry.

METHODS

We assembled genotype data for women of African ancestry, including 9241 case subjects and 10 193 control subjects. We evaluated associations of 179- and 313-variant PRSs with overall and subtype-specific BC risk. PRS discriminatory accuracy was assessed using area under the receiver operating characteristic curve. We also evaluated a recalibrated PRS, replacing the index variant with variants in each region that better captured risk in women of African ancestry and estimated lifetime absolute risk of BC in African Americans by PRS category.

RESULTS

For overall BC, the odds ratio per SD of the 313-variant PRS (PRS313) was 1.27 (95% confidence interval [CI] = 1.23 to 1.31), with an area under the receiver operating characteristic curve of 0.571 (95% CI = 0.562 to 0.579). Compared with women with average risk (40th-60th PRS percentile), women in the top decile of PRS313 had a 1.54-fold increased risk (95% CI = 1.38-fold to 1.72-fold). By age 85 years, the absolute risk of overall BC was 19.6% for African American women in the top 1% of PRS313 and 6.7% for those in the lowest 1%. The recalibrated PRS did not improve BC risk prediction.

CONCLUSION

The PRSs stratify BC risk in women of African ancestry, with attenuated performance compared with that reported in European, Asian, and Latina populations. Future work is needed to improve BC risk stratification for women of African ancestry.

Understanding Racial Disparities in Gastrointestinal Cancer Outcomes: Lack of Surgery Contributes to Lower Survival in African American Patients

BACKGROUND

Race/ethnicity-related differences in rates of cancer surgery and cancer mortality have been observed for gastrointestinal (GI) cancers. This study aims to estimate the extent to which differences in receipt of surgery explain racial/ethnic disparities in cancer survival.

METHODS

The National Cancer Database was used to obtain data for patients diagnosed with stage I-III mid-esophageal, distal esophagus/gastric cardia (DEGC), noncardia gastric, pancreatic, and colorectal cancer in years 2004-2015. Mediation analysis was used to identify variables influencing the relationship between race/ethnicity and mortality, including surgery.

RESULTS

A total of 600,063 patients were included in the study: 3.5% mid-esophageal, 12.4% DEGC, 4.9% noncardia gastric, 17.0% pancreatic, 40.1% colon, and 22.0% rectal cancers. The operative rates for Black patients were low relative to White patients, with absolute differences of 21.0%, 19.9%, 2.3%, 8.3%, 1.6%, and 7.7%. Adjustment for age, stage, and comorbidities revealed even lower odds of receiving surgery for Black patients compared with White patients. The observed HRs for Black patients compared with White patients ranged from 1.01 to 1.42. Mediation analysis showed that receipt of surgery and socioeconomic factors had greatest influence on the survival disparity.

CONCLUSIONS

The results of this study indicate that Black patients appear to be undertreated compared with White patients for GI cancers. The disproportionately low operative rates contribute to the known survival disparity between Black and White patients.

IMPACT

Interventions to reduce barriers to surgery for Black patients should be promoted to reduce disparities in GI cancer outcomes.See related commentary by Hébert, p. 438.

Exploring the effects of genetic variation on gene regulation in cancer in the context of 3D genome structure.. [DOI: 10.1101/2020.10.06.328567] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Numerous genome-wide association studies (GWAS) conducted to date revealed genetic variants associated with various diseases, including breast and prostate cancers. Despite the availability of these large-scale data, relatively few variants have been functionally characterized, mainly because the majority of single-nucleotide polymorphisms (SNPs) map to the non-coding regions of the human genome. The functional characterization of these non-coding variants and the identification of their target genes remain challenging. In this communication, we explore the potential functional mechanisms of non-coding SNPs by integrating GWAS with the high-resolution chromosome conformation capture (Hi-C) data for breast and prostate cancers. We show that more genetic variants map to regulatory elements through the 3D genome structure than the 1D linear genome lacking physical chromatin interactions. Importantly, the association of enhancers, transcription factors, and their target genes with breast and prostate cancers tends to be higher when these regulatory elements are mapped to high-risk SNPs through spatial interactions compared to simply using a linear proximity. Finally, we demonstrate that topologically associating domains (TADs) carrying high-risk SNPs also contain gene regulatory elements whose association with cancer is generally higher than those belonging to control TADs containing no high-risk variants. Our results suggest that many SNPs may contribute to the cancer development by affecting the expression of certain tumor-related genes through long-range chromatin interactions with gene regulatory elements. Integrating large-scale genetic datasets with the 3D genome structure offers an attractive and unique approach to systematically investigate the functional mechanisms of genetic variants in disease risk and progression.

Identification of novel epithelial ovarian cancer loci in women of African ancestry

Women of African ancestry have lower incidence of epithelial ovarian cancer (EOC) yet worse survival compared to women of European ancestry. We conducted a genome-wide association study in African ancestry women with 755 EOC cases, including 537 high-grade serous ovarian carcinomas (HGSOC) and 1,235 controls. We identified four novel loci with suggestive evidence of association with EOC (p < 1 × 10-6 ), including rs4525119 (intronic to AKR1C3), rs7643459 (intronic to LOC101927394), rs4286604 (12 kb 3' of UGT2A2) and rs142091544 (5 kb 5' of WWC1). For HGSOC, we identified six loci with suggestive evidence of association including rs37792 (132 kb 5' of follistatin [FST]), rs57403204 (81 kb 3' of MAGEC1), rs79079890 (LOC105376360 intronic), rs66459581 (5 kb 5' of PRPSAP1), rs116046250 (GABRG3 intronic) and rs192876988 (32 kb 3' of GK2). Among the identified variants, two are near genes known to regulate hormones and diseases of the ovary (AKR1C3 and FST), and two are linked to cancer (AKR1C3 and MAGEC1). In follow-up studies of the 10 identified variants, the GK2 region SNP, rs192876988, showed an inverse association with EOC in European ancestry women (p = 0.002), increased risk of ER positive breast cancer in African ancestry women (p = 0.027) and decreased expression of GK2 in HGSOC tissue from African ancestry women (p = 0.004). A European ancestry-derived polygenic risk score showed positive associations with EOC and HGSOC in women of African ancestry suggesting shared genetic architecture. Our investigation presents evidence of variants for EOC shared among European and African ancestry women and identifies novel EOC risk loci in women of African ancestry.

Using 3D epigenomic maps of primary olfactory neuronal cells from living individuals to understand gene regulation

As part of PsychENCODE, we developed a three-dimensional (3D) epigenomic map of primary cultured neuronal cells derived from olfactory neuroepithelium (CNON). We mapped topologically associating domains and high-resolution chromatin interactions using Hi-C and identified regulatory elements using chromatin immunoprecipitation and nucleosome positioning assays. Using epigenomic datasets from biopsies of 63 living individuals, we found that epigenetic marks at distal regulatory elements are more variable than marks at proximal regulatory elements. By integrating genotype and metadata, we identified enhancers that have different levels corresponding to differences in genetic variation, gender, smoking, and schizophrenia. Motif searches revealed that many CNON enhancers are bound by neuronal-related transcription factors. Last, we combined 3D epigenomic maps and gene expression profiles to predict enhancer-target gene interactions on a genome-wide scale. This study not only provides a framework for understanding individual epigenetic variation using a primary cell model system but also contributes valuable data resources for epigenomic studies of neuronal epithelium.

Validation of a genetic risk score for Arkansas women of color

African American women in the state of Arkansas have high breast cancer mortality rates. Breast cancer risk assessment tools developed for African American underestimate breast cancer risk. Combining African American breast cancer associated single-nucleotide polymorphisms (SNPs) into breast cancer risk algorithms may improve individualized estimates of a woman's risk of developing breast cancer and enable improved recommendation of screening and chemoprevention for women at high risk. The goal of this study was to confirm with an independent dataset consisting of Arkansas women of color, whether a genetic risk score derived from common breast cancer susceptibility SNPs can be combined with a clinical risk estimate provided by the Breast Cancer Risk Assessment Tool (BCRAT) to produce a more accurate individualized breast cancer risk estimate. A population-based cohort of African American women representative of Arkansas consisted of 319 cases and 559 controls for this study. Five-year and lifetime risks from the BCRAT were measured and combined with a risk score based on 75 independent susceptibility SNPs in African American women. We used the odds ratio (OR) per adjusted standard deviation to evaluate the improvement in risk estimates produced by combining the polygenic risk score (PRS) with 5-year and lifetime risk scores estimated using BCRAT. For 5-year risk OR per standard deviation increased from 1.84 to 2.08 with the addition of the polygenic risk score and from 1.79 to 2.07 for the lifetime risk score. Reclassification analysis indicated that 13% of cases had their 5-year risk increased above the 1.66% guideline threshold (NRI = 0.020 (95% CI -0.040, 0.080)) and 6.3% of cases had their lifetime risk increased above the 20% guideline threshold by the addition of the polygenic risk score (NRI = 0.034 (95% CI 0.000, 0.070)). Our data confirmed that discriminatory accuracy of BCRAT is improved for African American women in Arkansas with the inclusion of specific SNP breast cancer risk alleles.

Bridging the Data Gap in Breast Cancer Risk Assessment to Enable Widespread Clinical Implementation across the Multiethnic Landscape of the US

Breast cancer remains the second leading cause of cancer death among women and is the most commonly diagnosed cancer in women. Breast cancer risk assessment has been clinically available for nearly 30 years yet is under-utilized in practice for multiple reasons. Incorporation of polygenic risk as well as breast density measurements, promise to increase the accuracy of risk assessment. With that comes the hope that both prevention and screening become more personalized and thus more effective. Incidence rates have been static over the past 15 years and have even increased slightly in African American and Asian/Pacific Islander populations despite the robust data on breast cancer risk reduction measures that exist. Current challenges in reducing breast cancer incidence begin with robust data curation that allows for appropriate risk stratification across our multiethnic population and conclude with the implementation of prevention strategies within our fractured healthcare system.

Genetic variation in the Hippo pathway and breast cancer risk in women of African ancestry

Gene expression changes within the Hippo pathway were found to be associated with large tumor size and metastasis in breast cancer. The combined effect of genetic variants in genes of this pathway may have a causal role in breast cancer development. We examined 7086 SNPs that were not highly correlated (r²  < 0.8) in 35 Hippo pathway genes using data from the genome-wide association study of breast cancer from the Root Consortium, which includes 3686 participants of African ancestry from Nigeria, United States of America, and Barbados: 1657 cases (403 estrogen receptor-positive [ER+], 374 ER-) and 2029 controls. Gene-level analyses were conducted using improved AdaJoint test for large-scale genetic association studies adjusting for age, study site and the first four eigenvectors from the principal component analysis. SNP-level analyses were conducted with logistic regression. The Hippo pathway was significantly associated with risk of ER+ breast cancer (pathway-level P = 0.019), with WWC1 (P_adj = 0.04) being the leading gene. The pathway-level significance was lost without WWC1 (P = 0.12). rs147106204 in the WWC1 gene was the most statistically significant SNP after gene-level adjustment for multiple comparisons (OR = 0.53, 95%CI = 0.41-0.70, P_adj  = 0.025). We found evidence of an association between genetic variations in the Hippo pathway and ER+ breast cancer. Moreover, WWC1 was identified as the most important genetic susceptibility locus highlighting the importance of genetic epidemiology studies of breast cancer in understudied populations.