Hatef Darabi

Genome-wide association study of follicular lymphoma identifies a risk locus at 6p21.32

To identify susceptibility loci for non-Hodgkin lymphoma subtypes, we conducted a three-stage genome-wide association study. We identified two variants associated with follicular lymphoma at 6p21.32 (rs10484561, combined P = 1.12 × 10−29 and rs7755224, combined P = 2.00 × 10−19; r2 = 1.0), supporting the idea that major histocompatibility complex genetic variation influences follicular lymphoma susceptibility. We also found confirmatory evidence of a previously reported association between chronic lymphocytic leukemia/small lymphocytic lymphoma and rs735665 (combined P = 4.24 × 10−9).

GWAS of follicular lymphoma reveals allelic heterogeneity at 6p21.32 and suggests shared genetic susceptibility with diffuse large B-cell lymphoma.

Non-Hodgkin lymphoma (NHL) represents a diverse group of hematological malignancies, of which follicular lymphoma (FL) is a prevalent subtype. A previous genome-wide association study has established a marker, rs10484561 in the human leukocyte antigen (HLA) class II region on 6p21.32 associated with increased FL risk. Here, in a three-stage genome-wide association study, starting with a genome-wide scan of 379 FL cases and 791 controls followed by validation in 1,049 cases and 5,790 controls, we identified a second independent FL–associated locus on 6p21.32, rs2647012 (ORcombined = 0.64, Pcombined = 2×10−21) located 962 bp away from rs10484561 (r2<0.1 in controls). After mutual adjustment, the associations at the two SNPs remained genome-wide significant (rs2647012:ORadjusted = 0.70, Padjusted = 4×10−12; rs10484561:ORadjusted = 1.64, Padjusted = 5×10−15). Haplotype and coalescence analyses indicated that rs2647012 arose on an evolutionarily distinct haplotype from that of rs10484561 and tags a novel allele with an opposite (protective) effect on FL risk. Moreover, in a follow-up analysis of the top 6 FL–associated SNPs in 4,449 cases of other NHL subtypes, rs10484561 was associated with risk of diffuse large B-cell lymphoma (ORcombined = 1.36, Pcombined = 1.4×10−7). Our results reveal the presence of allelic heterogeneity within the HLA class II region influencing FL susceptibility and indicate a possible shared genetic etiology with diffuse large B-cell lymphoma. These findings suggest that the HLA class II region plays a complex yet important role in NHL.

Breast cancer risk prediction and individualised screening based on common genetic variation and breast density measurement

IntroductionOver the last decade several breast cancer risk alleles have been identified which has led to an increased interest in individualised risk prediction for clinical purposes.MethodsWe investigate the performance of an up-to-date 18 breast cancer risk single-nucleotide polymorphisms (SNPs), together with mammographic percentage density (PD), body mass index (BMI) and clinical risk factors in predicting absolute risk of breast cancer, empirically, in a well characterised Swedish case-control study of postmenopausal women. We examined the efficiency of various prediction models at a population level for individualised screening by extending a recently proposed analytical approach for estimating number of cases captured.ResultsThe performance of a risk prediction model based on an initial set of seven breast cancer risk SNPs is improved by additionally including eleven more recently established breast cancer risk SNPs (P = 4.69 × 10-4). Adding mammographic PD, BMI and all 18 SNPs to a Swedish Gail model improved the discriminatory accuracy (the AUC statistic) from 55% to 62%. The net reclassification improvement was used to assess improvement in classification of women into low, intermediate, and high categories of 5-year risk (P = 8.93 × 10-9). For scenarios we considered, we estimated that an individualised screening strategy based on risk models incorporating clinical risk factors, mammographic density and SNPs, captures 10% more cases than a screening strategy using the same resources, based on age alone. Estimates of numbers of cases captured by screening stratified by age provide insight into how individualised screening programs might appear in practice.ConclusionsTaken together, genetic risk factors and mammographic density offer moderate improvements to clinical risk factor models for predicting breast cancer.

Multi-variant pathway association analysis reveals the importance of genetic determinants of estrogen metabolism in breast and endometrial cancer susceptibility.

Despite the central role of estrogen exposure in breast and endometrial cancer development and numerous studies of genes in the estrogen metabolic pathway, polymorphisms within the pathway have not been consistently associated with these cancers. We posit that this is due to the complexity of multiple weak genetic effects within the metabolic pathway that can only be effectively detected through multi-variant analysis. We conducted a comprehensive association analysis of the estrogen metabolic pathway by interrogating 239 tagSNPs within 35 genes of the pathway in three tumor samples. The discovery sample consisted of 1,596 breast cancer cases, 719 endometrial cancer cases, and 1,730 controls from Sweden; and the validation sample included 2,245 breast cancer cases and 1,287 controls from Finland. We performed admixture maximum likelihood (AML)–based global tests to evaluate the cumulative effect from multiple SNPs within the whole metabolic pathway and three sub-pathways for androgen synthesis, androgen-to-estrogen conversion, and estrogen removal. In the discovery sample, although no single polymorphism was significant after correction for multiple testing, the pathway-based AML global test suggested association with both breast (p global = 0.034) and endometrial (p global = 0.052) cancers. Further testing revealed the association to be focused on polymorphisms within the androgen-to-estrogen conversion sub-pathway, for both breast (p global = 0.008) and endometrial cancer (p global = 0.014). The sub-pathway association was validated in the Finnish sample of breast cancer (p global = 0.015). Further tumor subtype analysis demonstrated that the association of the androgen-to-estrogen conversion sub-pathway was confined to postmenopausal women with sporadic estrogen receptor positive tumors (p global = 0.0003). Gene-based AML analysis suggested CYP19A1 and UGT2B4 to be the major players within the sub-pathway. Our study indicates that the composite genetic determinants related to the androgen–estrogen conversion are important for the induction of two hormone-associated cancers, particularly for the hormone-driven breast tumour subtypes.

ESR1 and EGF genetic variation in relation to breast cancer risk and survival

IntroductionOestrogen exposure is a central factor in the development of breast cancer. Oestrogen receptor alpha (ESR1) is the main mediator of oestrogen effect in breast epithelia and has also been shown to be activated by epidermal growth factor (EGF). We sought to determine if common genetic variation in the ESR1 and EGF genes affects breast cancer risk, tumour characteristics or breast cancer survival.MethodsWe genotyped 157 single nucleotide polymorphisms (SNPs) in ESR1 and 54 SNPs in EGF in 92 Swedish controls and selected haplotype tagging SNPs (tagSNPs) that could predict both single SNP and haplotype variation in the genes with an R2 of at least 0.8. The tagSNPs were genotyped in 1,590 breast cancer cases and 1,518 controls, and their association with breast cancer risk, tumour characteristics and survival were assessed using unconditional logistic regression models, Cox proportional hazard models and haplotype analysis.ResultsThe single tagSNP analysis did not reveal association evidence for breast cancer risk, tumour characteristics, or survival. A multi-locus analysis of five adjacent tagSNPs suggested a region in ESR1 (between rs3003925 and rs2144025) for association with breast cancer risk (p = 0.001), but the result did not withstand adjustment for multiple comparisons (p = 0.086). A similar region was also implicated by haplotype analyses, but its significance needs to be verified by follow-up analysis.ConclusionOur results do not support a strong association between common variants in the ESR1 and EGF genes and breast cancer risk, tumour characteristics or survival.

A genome-wide association scan on estrogen receptor-negative breast cancer.

IntroductionBreast cancer is a heterogeneous disease and may be characterized on the basis of whether estrogen receptors (ER) are expressed in the tumour cells. ER status of breast cancer is important clinically, and is used both as a prognostic indicator and treatment predictor. In this study, we focused on identifying genetic markers associated with ER-negative breast cancer risk.MethodsWe conducted a genome-wide association analysis of 285,984 single nucleotide polymorphisms (SNPs) genotyped in 617 ER-negative breast cancer cases and 4,583 controls. We also conducted a genome-wide pathway analysis on the discovery dataset using permutation-based tests on pre-defined pathways. The extent of shared polygenic variation between ER-negative and ER-positive breast cancers was assessed by relating risk scores, derived using ER-positive breast cancer samples, to disease state in independent, ER-negative breast cancer cases.ResultsAssociation with ER-negative breast cancer was not validated for any of the five most strongly associated SNPs followed up in independent studies (1,011 ER-negative breast cancer cases, 7,604 controls). However, an excess of small P-values for SNPs with known regulatory functions in cancer-related pathways was found (global P = 0.052). We found no evidence to suggest that ER-negative breast cancer shares a polygenic basis to disease with ER-positive breast cancer.ConclusionsER-negative breast cancer is a distinct breast cancer subtype that merits independent analyses. Given the clinical importance of this phenotype and the likelihood that genetic effect sizes are small, greater sample sizes and further studies are required to understand the etiology of ER-negative breast cancers.

Coding Variants at Hexa-allelic Amino Acid 13 of HLA-DRB1 Explain Independent SNP Associations with Follicular Lymphoma Risk

Non-Hodgkin lymphoma represents a diverse group of blood malignancies, of which follicular lymphoma (FL) is a common subtype. Previous genome-wide association studies (GWASs) have identified in the human leukocyte antigen (HLA) class II region multiple independent SNPs that are significantly associated with FL risk. To dissect these signals and determine whether coding variants in HLA genes are responsible for the associations, we conducted imputation, HLA typing, and sequencing in three independent populations for a total of 689 cases and 2,446 controls. We identified a hexa-allelic amino acid polymorphism at position 13 of the HLA-DR beta chain that showed the strongest association with FL within the major histocompatibility complex (MHC) region (multiallelic p = 2.3 × 10⁻¹⁵). Out of six possible amino acids that occurred at that position within the population, we classified two as high risk (Tyr and Phe), two as low risk (Ser and Arg), and two as moderate risk (His and Gly). There was a 4.2-fold difference in risk (95% confidence interval = 2.9-6.1) between subjects carrying two alleles encoding high-risk amino acids and those carrying two alleles encoding low-risk amino acids (p = 1.01 × 10⁻¹⁴). This coding variant might explain the complex SNP associations identified by GWASs and suggests a common HLA-DR antigen-driven mechanism for the pathogenesis of FL and rheumatoid arthritis.

Five endometrial cancer risk loci identified through genome-wide association analysis

We conducted a meta-analysis of three endometrial cancer genome-wide association studies (GWAS) and two follow-up phases totaling 7,737 endometrial cancer cases and 37,144 controls of European ancestry. Genome-wide imputation and meta-analysis identified five new risk loci of genome-wide significance at likely regulatory regions on chromosomes 13q22.1 (rs11841589, near KLF5), 6q22.31 (rs13328298, in LOC643623 and near HEY2 and NCOA7), 8q24.21 (rs4733613, telomeric to MYC), 15q15.1 (rs937213, in EIF2AK4, near BMF) and 14q32.33 (rs2498796, in AKT1, near SIVA1). We also found a second independent 8q24.21 signal (rs17232730). Functional studies of the 13q22.1 locus showed that rs9600103 (pairwise r2 = 0.98 with rs11841589) is located in a region of active chromatin that interacts with the KLF5 promoter region. The rs9600103[T] allele that is protective in endometrial cancer suppressed gene expression in vitro, suggesting that regulation of the expression of KLF5, a gene linked to uterine development, is implicated in tumorigenesis. These findings provide enhanced insight into the genetic and biological basis of endometrial cancer.

CYP19A1 fine-mapping and Mendelian randomization: estradiol is causal for endometrial cancer

Candidate gene studies have reported CYP19A1 variants to be associated with endometrial cancer and with estradiol (E2) concentrations. We analyzed 2937 single nucleotide polymorphisms (SNPs) in 6608 endometrial cancer cases and 37 925 controls and report the first genome wide-significant association between endometrial cancer and a CYP19A1 SNP (rs727479 in intron 2, P=4.8×10−11). SNP rs727479 was also among those most strongly associated with circulating E2 concentrations in 2767 post-menopausal controls (P=7.4×10−8). The observed endometrial cancer odds ratio per rs727479 A-allele (1.15, CI=1.11–1.21) is compatible with that predicted by the observed effect on E2 concentrations (1.09, CI=1.03–1.21), consistent with the hypothesis that endometrial cancer risk is driven by E2. From 28 candidate-causal SNPs, 12 co-located with three putative gene-regulatory elements and their risk alleles associated with higher CYP19A1 expression in bioinformatical analyses. For both phenotypes, the associations with rs727479 were stronger among women with a higher BMI (Pinteraction=0.034 and 0.066 respectively), suggesting a biologically plausible gene-environment interaction.

Identification of new genetic susceptibility Loci for breast cancer through consideration of gene-environment interactions

Genes that alter disease risk only in combination with certain environmental exposures may not be detected in genetic association analysis. By using methods accounting for gene‐environment (G × E) interaction, we aimed to identify novel genetic loci associated with breast cancer risk. Up to 34,475 cases and 34,786 controls of European ancestry from up to 23 studies in the Breast Cancer Association Consortium were included. Overall, 71,527 single nucleotide polymorphisms (SNPs), enriched for association with breast cancer, were tested for interaction with 10 environmental risk factors using three recently proposed hybrid methods and a joint test of association and interaction. Analyses were adjusted for age, study, population stratification, and confounding factors as applicable. Three SNPs in two independent loci showed statistically significant association: SNPs rs10483028 and rs2242714 in perfect linkage disequilibrium on chromosome 21 and rs12197388 in ARID1B on chromosome 6. While rs12197388 was identified using the joint test with parity and with age at menarche (P‐values = 3 × 10−07), the variants on chromosome 21 q22.12, which showed interaction with adult body mass index (BMI) in 8,891 postmenopausal women, were identified by all methods applied. SNP rs10483028 was associated with breast cancer in women with a BMI below 25 kg/m2 (OR = 1.26, 95% CI 1.15–1.38) but not in women with a BMI of 30 kg/m2 or higher (OR = 0.89, 95% CI 0.72–1.11, P for interaction = 3.2 × 10−05). Our findings confirm comparable power of the recent methods for detecting G × E interaction and the utility of using G × E interaction analyses to identify new susceptibility loci.