Identification of ten variants associated with risk of estrogen-receptor-negative breast cancer

Most common breast cancer susceptibility variants have been identified through genome-wide association studies (GWAS) of predominantly estrogen receptor (ER)-positive disease. We conducted a GWAS using 21,468 ER-negative cases and 100,594 controls combined with 18,908 BRCA1 mutation carriers (9,414 with breast cancer), all of European origin. We identified independent associations at P < 5 × 10-8 with ten variants at nine new loci. At P < 0.05, we replicated associations with 10 of 11 variants previously reported in ER-negative disease or BRCA1 mutation carrier GWAS and observed consistent associations with ER-negative disease for 105 susceptibility variants identified by other studies. These 125 variants explain approximately 16% of the familial risk of this breast cancer subtype. There was high genetic correlation (0.72) between risk of ER-negative breast cancer and breast cancer risk for BRCA1 mutation carriers. These findings may lead to improved risk prediction and inform further fine-mapping and functional work to better understand the biological basis of ER-negative breast cancer.

Association analysis identifies 65 new breast cancer risk loci

Breast cancer risk is influenced by rare coding variants in susceptibility genes, such as BRCA1, and many common, mostly non-coding variants. However, much of the genetic contribution to breast cancer risk remains unknown. Here we report the results of a genome-wide association study of breast cancer in 122,977 cases and 105,974 controls of European ancestry and 14,068 cases and 13,104 controls of East Asian ancestry. We identified 65 new loci that are associated with overall breast cancer risk at P < 5 × 10-8. The majority of credible risk single-nucleotide polymorphisms in these loci fall in distal regulatory elements, and by integrating in silico data to predict target genes in breast cells at each locus, we demonstrate a strong overlap between candidate target genes and somatic driver genes in breast tumours. We also find that heritability of breast cancer due to all single-nucleotide polymorphisms in regulatory features was 2-5-fold enriched relative to the genome-wide average, with strong enrichment for particular transcription factor binding sites. These results provide further insight into genetic susceptibility to breast cancer and will improve the use of genetic risk scores for individualized screening and prevention.

The OncoArray Consortium: A Network for Understanding the Genetic Architecture of Common Cancers

BACKGROUND

Common cancers develop through a multistep process often including inherited susceptibility. Collaboration among multiple institutions, and funding from multiple sources, has allowed the development of an inexpensive genotyping microarray, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits.

METHODS

The OncoArray can be genotyped using a novel technology developed by Illumina to facilitate efficient genotyping. The consortium developed standard approaches for selecting SNPs for study, for quality control of markers, and for ancestry analysis. The array was genotyped at selected sites and with prespecified replicate samples to permit evaluation of genotyping accuracy among centers and by ethnic background.

RESULTS

The OncoArray consortium genotyped 447,705 samples. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis.

CONCLUSIONS

Results from these analyses will enable researchers to identify new susceptibility loci, perform fine-mapping of new or known loci associated with either single or multiple cancers, assess the degree of overlap in cancer causation and pleiotropic effects of loci that have been identified for disease-specific risk, and jointly model genetic, environmental, and lifestyle-related exposures.

IMPACT

Ongoing analyses will shed light on etiology and risk assessment for many types of cancer. Cancer Epidemiol Biomarkers Prev; 26(1); 126-35. ©2016 AACR.

Variations in 5-methylcytosine and 5-hydroxymethylcytosine among human brain, blood, and saliva using oxBS and the Infinium MethylationEPIC array

Investigating 5-methylcytosine (5mC) has led to many hypotheses regarding molecular mechanism underlying human diseases and disorders. Many of these studies, however, utilize bisulfite conversion alone, which cannot distinguish 5mC from its recently discovered oxidative product, 5-hydroxymethylcytosine (5hmC). Furthermore, previous array-based technologies do not have the necessary probes to adequately investigate both modifications simultaneously. In this manuscript, we used technical replicates of DNA from human brain, human blood, and human saliva, in combination with oxidative bisulfite conversion and Illumina's Infinium MethylationEPIC array, to analyze 5mC and 5hmC at more than 650 000 and 450 000 relevant loci, respectively, in the human genome. We show the presence of loci with detectable 5mC and 5hmC to be equally distributed across chromosomes and genomic features, while also being present in genomic regions with transcriptional regulatory properties. We also describe 2528 5hmC sites common across tissue types that show a strong association with immune-related functions. Lastly, in human brain, we show that 5hmC accounts for one-third of the total signal from bisulfite-converted data. As such, not only do our results confirm the efficacy and sensitivity of pairing oxidative bisulfite conversion and the EPIC array to detect 5mC and 5hmC in all three tissue types, but they also highlight the importance of dissociating 5hmC from 5mC in future studies related to cytosine modifications.

Inherited variants in the inner centromere protein (INCENP) gene of the chromosomal passenger complex contribute to the susceptibility of ER-negative breast cancer

The chromosomal passenger complex (CPC) plays a pivotal role in the regulation of cell division. Therefore, inherited CPC variability could influence tumor development. The present candidate gene approach investigates the relationship between single nucleotide polymorphisms (SNPs) in genes encoding key CPC components and breast cancer risk. Fifteen SNPs in four CPC genes (INCENP, AURKB, BIRC5 and CDCA8) were genotyped in 88 911 European women from 39 case-control studies of the Breast Cancer Association Consortium. Possible associations were investigated in fixed-effects meta-analyses. The synonymous SNP rs1675126 in exon 7 of INCENP was associated with overall breast cancer risk [per A allele odds ratio (OR) 0.95, 95% confidence interval (CI) 0.92-0.98, P = 0.007] and particularly with estrogen receptor (ER)-negative breast tumors (per A allele OR 0.89, 95% CI 0.83-0.95, P = 0.0005). SNPs not directly genotyped were imputed based on 1000 Genomes. The SNPs rs1047739 in the 3' untranslated region and rs144045115 downstream of INCENP showed the strongest association signals for overall (per T allele OR 1.03, 95% CI 1.00-1.06, P = 0.0009) and ER-negative breast cancer risk (per A allele OR 1.06, 95% CI 1.02-1.10, P = 0.0002). Two genotyped SNPs in BIRC5 were associated with familial breast cancer risk (top SNP rs2071214: per G allele OR 1.12, 95% CI 1.04-1.21, P = 0.002). The data suggest that INCENP in the CPC pathway contributes to ER-negative breast cancer susceptibility in the European population. In spite of a modest contribution of CPC-inherited variants to the total burden of sporadic and familial breast cancer, their potential as novel targets for breast cancer treatment should be further investigated.

Fine-mapping identifies multiple prostate cancer risk loci at 5p15, one of which associates with TERT expression

Associations between single nucleotide polymorphisms (SNPs) at 5p15 and multiple cancer types have been reported. We have previously shown evidence for a strong association between prostate cancer (PrCa) risk and rs2242652 at 5p15, intronic in the telomerase reverse transcriptase (TERT) gene that encodes TERT. To comprehensively evaluate the association between genetic variation across this region and PrCa, we performed a fine-mapping analysis by genotyping 134 SNPs using a custom Illumina iSelect array or Sequenom MassArray iPlex, followed by imputation of 1094 SNPs in 22 301 PrCa cases and 22 320 controls in The PRACTICAL consortium. Multiple stepwise logistic regression analysis identified four signals in the promoter or intronic regions of TERT that independently associated with PrCa risk. Gene expression analysis of normal prostate tissue showed evidence that SNPs within one of these regions also associated with TERT expression, providing a potential mechanism for predisposition to disease.

GWAS meta-analysis and replication identifies three new susceptibility loci for ovarian cancer

Genome-wide association studies (GWAS) have identified four susceptibility loci for epithelial ovarian cancer (EOC), with another two suggestive loci reaching near genome-wide significance. We pooled data from a GWAS conducted in North America with another GWAS from the UK. We selected the top 24,551 SNPs for inclusion on the iCOGS custom genotyping array. We performed follow-up genotyping in 18,174 individuals with EOC (cases) and 26,134 controls from 43 studies from the Ovarian Cancer Association Consortium. We validated the two loci at 3q25 and 17q21 that were previously found to have associations close to genome-wide significance and identified three loci newly associated with risk: two loci associated with all EOC subtypes at 8q21 (rs11782652, P = 5.5 × 10(-9)) and 10p12 (rs1243180, P = 1.8 × 10(-8)) and another locus specific to the serous subtype at 17q12 (rs757210, P = 8.1 × 10(-10)). An integrated molecular analysis of genes and regulatory regions at these loci provided evidence for functional mechanisms underlying susceptibility and implicated CHMP4C in the pathogenesis of ovarian cancer.

Identification of 23 new prostate cancer susceptibility loci using the iCOGS custom genotyping array

Prostate cancer is the most frequently diagnosed cancer in males in developed countries. To identify common prostate cancer susceptibility alleles, we genotyped 211,155 SNPs on a custom Illumina array (iCOGS) in blood DNA from 25,074 prostate cancer cases and 24,272 controls from the international PRACTICAL Consortium. Twenty-three new prostate cancer susceptibility loci were identified at genome-wide significance (P < 5 × 10(-8)). More than 70 prostate cancer susceptibility loci, explaining ∼30% of the familial risk for this disease, have now been identified. On the basis of combined risks conferred by the new and previously known risk loci, the top 1% of the risk distribution has a 4.7-fold higher risk than the average of the population being profiled. These results will facilitate population risk stratification for clinical studies.

Large-scale genotyping identifies 41 new loci associated with breast cancer risk

Breast cancer is the most common cancer among women. Common variants at 27 loci have been identified as associated with susceptibility to breast cancer, and these account for ∼9% of the familial risk of the disease. We report here a meta-analysis of 9 genome-wide association studies, including 10,052 breast cancer cases and 12,575 controls of European ancestry, from which we selected 29,807 SNPs for further genotyping. These SNPs were genotyped in 45,290 cases and 41,880 controls of European ancestry from 41 studies in the Breast Cancer Association Consortium (BCAC). The SNPs were genotyped as part of a collaborative genotyping experiment involving four consortia (Collaborative Oncological Gene-environment Study, COGS) and used a custom Illumina iSelect genotyping array, iCOGS, comprising more than 200,000 SNPs. We identified SNPs at 41 new breast cancer susceptibility loci at genome-wide significance (P < 5 × 10(-8)). Further analyses suggest that more than 1,000 additional loci are involved in breast cancer susceptibility.

Genome-wide association study in BRCA1 mutation carriers identifies novel loci associated with breast and ovarian cancer risk

BRCA1-associated breast and ovarian cancer risks can be modified by common genetic variants. To identify further cancer risk-modifying loci, we performed a multi-stage GWAS of 11,705 BRCA1 carriers (of whom 5,920 were diagnosed with breast and 1,839 were diagnosed with ovarian cancer), with a further replication in an additional sample of 2,646 BRCA1 carriers. We identified a novel breast cancer risk modifier locus at 1q32 for BRCA1 carriers (rs2290854, P = 2.7 × 10(-8), HR = 1.14, 95% CI: 1.09-1.20). In addition, we identified two novel ovarian cancer risk modifier loci: 17q21.31 (rs17631303, P = 1.4 × 10(-8), HR = 1.27, 95% CI: 1.17-1.38) and 4q32.3 (rs4691139, P = 3.4 × 10(-8), HR = 1.20, 95% CI: 1.17-1.38). The 4q32.3 locus was not associated with ovarian cancer risk in the general population or BRCA2 carriers, suggesting a BRCA1-specific association. The 17q21.31 locus was also associated with ovarian cancer risk in 8,211 BRCA2 carriers (P = 2×10(-4)). These loci may lead to an improved understanding of the etiology of breast and ovarian tumors in BRCA1 carriers. Based on the joint distribution of the known BRCA1 breast cancer risk-modifying loci, we estimated that the breast cancer lifetime risks for the 5% of BRCA1 carriers at lowest risk are 28%-50% compared to 81%-100% for the 5% at highest risk. Similarly, based on the known ovarian cancer risk-modifying loci, the 5% of BRCA1 carriers at lowest risk have an estimated lifetime risk of developing ovarian cancer of 28% or lower, whereas the 5% at highest risk will have a risk of 63% or higher. Such differences in risk may have important implications for risk prediction and clinical management for BRCA1 carriers.

Identification and molecular characterization of a new ovarian cancer susceptibility locus at 17q21.31

Epithelial ovarian cancer (EOC) has a heritable component that remains to be fully characterized. Most identified common susceptibility variants lie in non-protein-coding sequences. We hypothesized that variants in the 3' untranslated region at putative microRNA (miRNA)-binding sites represent functional targets that influence EOC susceptibility. Here, we evaluate the association between 767 miRNA-related single-nucleotide polymorphisms (miRSNPs) and EOC risk in 18,174 EOC cases and 26,134 controls from 43 studies genotyped through the Collaborative Oncological Gene-environment Study. We identify several miRSNPs associated with invasive serous EOC risk (odds ratio=1.12, P=10(-8)) mapping to an inversion polymorphism at 17q21.31. Additional genotyping of non-miRSNPs at 17q21.31 reveals stronger signals outside the inversion (P=10(-10)). Variation at 17q21.31 is associated with neurological diseases, and our collaboration is the first to report an association with EOC susceptibility. An integrated molecular analysis in this region provides evidence for ARHGAP27 and PLEKHM1 as candidate EOC susceptibility genes.

Epigenetic analysis leads to identification of HNF1B as a subtype-specific susceptibility gene for ovarian cancer

HNF1B is overexpressed in clear cell epithelial ovarian cancer, and we observed epigenetic silencing in serous epithelial ovarian cancer, leading us to hypothesize that variation in this gene differentially associates with epithelial ovarian cancer risk according to histological subtype. Here we comprehensively map variation in HNF1B with respect to epithelial ovarian cancer risk and analyse DNA methylation and expression profiles across histological subtypes. Different single-nucleotide polymorphisms associate with invasive serous (rs7405776 odds ratio (OR)=1.13, P=3.1 × 10(-10)) and clear cell (rs11651755 OR=0.77, P=1.6 × 10(-8)) epithelial ovarian cancer. Risk alleles for the serous subtype associate with higher HNF1B-promoter methylation in these tumours. Unmethylated, expressed HNF1B, primarily present in clear cell tumours, coincides with a CpG island methylator phenotype affecting numerous other promoters throughout the genome. Different variants in HNF1B associate with risk of serous and clear cell epithelial ovarian cancer; DNA methylation and expression patterns are also notably distinct between these subtypes. These findings underscore distinct mechanisms driving different epithelial ovarian cancer histological subtypes.

Genetic variant near IRS1 is associated with type 2 diabetes, insulin resistance and hyperinsulinemia

Genome-wide association studies have identified common variants that only partially explain the genetic risk for type 2 diabetes (T2D). Using genome-wide association data from 1,376 French individuals, we identified 16,360 SNPs nominally associated with T2D and studied these SNPs in an independent sample of 4,977 French individuals. We then selected the 28 best hits for replication in 7,698 Danish subjects and identified 4 SNPs showing strong association with T2D, one of which (rs2943641, P = 9.3 x 10(-12), OR = 1.19) was located adjacent to the insulin receptor substrate 1 gene (IRS1). Unlike previously reported T2D risk loci, which predominantly associate with impaired beta cell function, the C allele of rs2943641 was associated with insulin resistance and hyperinsulinemia in 14,358 French, Danish and Finnish participants from population-based cohorts; this allele was also associated with reduced basal levels of IRS1 protein and decreased insulin induction of IRS1-associated phosphatidylinositol-3-OH kinase activity in human skeletal muscle biopsies.

Identification of susceptibility genes for complex diseases using pooling-based genome-wide association scans

The success of genome-wide association studies (GWAS) to identify risk loci of complex diseases is now well-established. One persistent major hurdle is the cost of those studies, which make them beyond the reach of most research groups. Performing GWAS on pools of DNA samples may be an effective strategy to reduce the costs of these studies. In this study, we performed pooling-based GWAS with more than 550,000 SNPs in two case-control cohorts consisting of patients with Type II diabetes (T2DM) and with chronic rhinosinusitis (CRS). In the T2DM study, the results of the pooling experiment were compared to individual genotypes obtained from a previously published GWAS. TCF7L2 and HHEX SNPs associated with T2DM by the traditional GWAS were among the top ranked SNPs in the pooling experiment. This dataset was also used to refine the best strategy to correctly identify SNPs that will remain significant based on individual genotyping. In the CRS study, the top hits from the pooling-based GWAS located within ten kilobases of known genes were validated by individual genotyping of 1,536 SNPs. Forty-one percent (598 out of the 1,457 SNPs that passed quality control) were associated with CRS at a nominal P value of 0.05, confirming the potential of pooling-based GWAS to identify SNPs that differ in allele frequencies between two groups of subjects. Overall, our results demonstrate that a pooling experiment on high-density genotyping arrays can accurately determine the minor allelic frequency as compared to individual genotyping and produce a list of top ranked SNPs that captures genuine allelic differences between a group of cases and controls. The low cost associated with a pooling-based GWAS clearly justifies its use in screening for genetic determinants of complex diseases.

Genetic control of alternative splicing in the TAP2 gene: possible implication in the genetics of type 1 diabetes

The transporter 2, ATP-binding cassette, subfamily B (TAP2) is involved in the transport of antigenic peptides to HLA molecules. Coding TAP2 polymorphisms shows a strong association with type 1 diabetes, but it is not clear whether this association may be entirely due to linkage disequilibrium with HLA DR and DQ. Functionally, rat Tap2 nonsynonymous single-nucleotide polymorphisms (nsSNPs) confer differential selectivity for antigenic peptides, but this was not shown to be the case for human TAP2 nsSNPs. In the human, differential peptide selectivity is rather conferred by two splicing isoforms with alternative carboxy terminals. Here, we tested the hypothesis that alleles at the coding SNPs favor different splicing isoforms, thus determining peptide selectivity indirectly. This may be the basis for independent contribution to the type 1 diabetes association. In RNA from heterozygous lymphoblastoid lines, we measured the relative abundance of each SNP haplotype in each isoform. In isoform NM_000544, the G (Ala) allele at 665 Thr>Ala (rs241447) is more than twice as abundant as A (Thr) (GA = 2.2 +/- 0.4, P = 1.5 x 10(-4)), while isoform NM_018833 is derived almost exclusively from chromosomes carrying A (AG = 18.1 +/- 5.6, P = 2.04 x 10(-7)). In 889 Canadian children with type 1 diabetes, differential transmission of parental TAP2 alleles persisted (P = 0.011) when analysis was confined to chromosomes carrying only DQ*02 alleles, which mark a conserved DR-DQ haplotype, thus eliminating most of the variation at DR-DQ. Thus, we present evidence of TAP2 association with type 1 diabetes that is independent of HLA DR-DQ and describe a plausible functional mechanism based on allele dependence of splicing into isoforms known to have differential peptide selectivities.

No association of type 1 diabetes with a functional polymorphism of the LRAP gene

AIMS/HYPOTHESIS

In a recent linkage analysis of genome-wide gene-expression patterns in lymphoblastoid lines, the gene encoding the leukocyte-derived arginine aminopeptidase (LRAP) was identified as having its expression levels modulated by one of the most pronounced effects in cis, mapping to a single-nucleotide polymorphism (rs2762). As this enzyme has an important role in processing antigenic peptides for the HLA I molecules, a variant that drastically modulates its expression levels might affect risk of autoimmunity. This study was designed to confirm that LRAP expression in B-cell derived lines is controlled by a haplotype marked by rs2762 and to see whether this would be the basis of an association with type 1 diabetes (T1D).

METHODS

Single-nucleotide primer extension (SNuPE) was adapted to determine the haplotype-specific expression. Genetic association was tested in 892 nuclear families with one T1D-affected offspring and two parents (2676 individuals).

RESULTS

All nine heterozygous RNA samples showed an eight-fold higher level of one haplotype over the other (7.97+/-0.99, p=1.33x10(-9)). However, no association of rs2762 with T1D was found by the transmission disequilibrium test (transmission ratio A/G=377/388, p=0.69).

CONCLUSIONS/INTERPRETATION

The genetically determined LRAP expression does not play significant roles in T1D. However, this dramatic genetic effect on LRAP expression justifies further investigation of association with other phenotypes, especially autoimmune and related to host defense to specific pathogens.

Lack of association of type 1 diabetes with the IL4R gene

AIMS/HYPOTHESIS

The association between IL4R and type 1 diabetes has been tested in many studies in recent years, with contradictory results. The aim of this study was to re-evaluate the genetic association in type 1 diabetic nuclear families of mixed European background.

SUBJECTS, MATERIALS AND METHODS

We genotyped six non-synonymous single-nucleotide polymorphisms (SNPs) of the IL4R gene in 830 nuclear families as specified above, including a French Canadian subset.

RESULTS

No association between type 1 diabetes and any SNP or haplotype was found by the transmission disequilibrium test.

CONCLUSIONS/INTERPRETATION

Previous positive reports may be due to population stratification as, according to HapMap data, allele frequencies in the IL4R region vary considerably by ethnicity.

Type 1 diabetes and the OAS gene cluster: association with splicing polymorphism or haplotype?

BACKGROUND

The 2',5'-oligoadenylate synthetase genes (OAS1, OAS2, and OAS3) map to human chromosome 12q24 and encode a family of enzymes pivotal to innate antiviral defence. Recently, the minor allele of an OAS1 single nucleotide polymorphism (SNP) that alters splicing (rs10774671) was found to be associated with increased enzymatic activity and, in a case-sibling control study, with type 1 diabetes (T1D).

METHODS

We have confirmed this T1D association in 784 nuclear families (two parents and at least one affected offspring) by the transmission disequilibrium test (TDT; G:A = 386:329, p = 0.033). However, because of linkage disequilibrium within OAS1 and with the other two OAS genes, functional attribution of the association to this SNP cannot be assumed. To help answer this question, we also genotyped two non-synonymous SNPs in OAS1 exons 3 and 7.

RESULTS

All three SNPs showed significant transmission distortion. Three of the eight possible haplotypes accounted for 98.4% of parental chromosomes and two of them carried the non-predisposing A allele at rs10774671. Parents heterozygous for these two haplotypes showed significant transmission distortion (p = 0.009) despite being homozygous at rs10774671.

CONCLUSIONS

We confirm the T1D association with rs10774671, but we conclude that it cannot be attributed (solely) to the splicing variant rs10774671. A serine/glycine substitution in OAS1 exon 3 is more likely a functional variant.