Gene expression interactions in breast cancer survival.

<p>(<b>A</b>) Kaplan–Meier survival curves based on categorization of <i>HMMR</i> (probe NM_012484) and <i>AURKA</i> (NM_003600) expression in tertiles (low, medium or high expression). For simplicity, only the tertiles for “high” <i>AURKA</i> are shown. The tumours with high expression levels for both genes were not those with the poorest prognosis. (<b>B</b>) Kaplan–Meier survival curves based on categorization of <i>HMMR</i> (NM_012484) and <i>TUBG1</i> (NM_016437) expression in tertiles (low, medium or high expression). For simplicity, only the tertiles for “high” <i>HMMR</i> are shown. The cases with high expression levels for both genes were those with the poorest prognosis.</p

The <i>HMMR</i> locus and breast cancer risk in <i>BRCA1</i> mutation carriers.

<p>(<b>A</b>) Forest plots showing rs299290 HRs and 95% CIs (retrospective likelihood trend estimation) for participating countries (relatively small sample sets are not shown) ordered by sample size. Left and right panels show results for <i>BRCA1</i> and <i>BRCA2</i> mutation carriers, respectively. The sizes of the rectangles are proportional to the corresponding country/study precision. (<b>B</b>) The rs299290-containing region, including the genes, variation and regulatory evidence mentioned in HMECs. Exons are marked by black-filled rectangles and the direction of transcription is marked by arrows in the genomic structure. The chromosome 5 positions (base pairs (bp)) and linkage disequilibrium structure from Caucasian HapMap individuals are also shown.</p

Potential GxG associated with breast cancer risk in <i>BRCA1/2</i> mutation carriers.

<p>*Each estimate is derived from the interaction term of a Cox regression model.</p><p>Potential GxG associated with breast cancer risk in <i>BRCA1/2</i> mutation carriers.</p

p-values of association (−log10 scale) with breast cancer risk in <i>BRCA2</i> carriers for genotyped and imputed SNPs in the <i>NEIL2</i> gene.

<p>SNP rs1466785 is indicated with a purple arrow and the best causal imputed SNPs, rs804276 and rs804271 are indicated with a red arrow. Colors represent the pariwise r². Plot generated with LocusZoom <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1004256#pgen.1004256-Pruim1" target="_blank">[42]</a> (<a href="http://csg.sph.umich.edu/locuszoom/" target="_blank">http://csg.sph.umich.edu/locuszoom/</a>).</p

Associations with breast and ovarian cancer risk for SNPs observed at p-trend<0.05 in stage II of the experiment.

a<p>Hazard Ratio per allele (1 df) estimated from the retrospective likelihood analysis.</p>b<p>Hazard Ratio under the genotype specific models (2df) estimated from the retrospective likelihood analysis.</p>c<p>p-values were based on the score test.</p>d<p>HR per allele of 1.69 and p-trend of 1×10⁻⁴ for <i>BRCA2</i> mutation carriers in stage I of the study.</p>e<p>HR per allele of 1.43 and p-trend of 0.01 for <i>BRCA1</i> mutation carriers in stage I of the study.</p>f<p>HR per allele of 1.30 and p-trend of 0.03 for <i>BRCA1</i> mutation carriers in stage I of the study.</p>g<p>HR per allele of 0.64 and p-trend of 0.057 for <i>BRCA2</i> mutation carriers in stage I of the study.</p>h<p>HR per allele of 1.25 and p-trend of 0.04 for <i>BRCA1</i> mutation carriers in stage I of the study.</p>i<p>HR per allele of 1.25 and p-trend of 0.058 for <i>BRCA2</i> mutation carriers in stage I of the study.</p>j<p>rs3093926 did not yield results under the genotype specific model due to the low minor allele frequency.</p><p>Complete description of results from stage I are included in Supplementary <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1004256#pgen.1004256.s002" target="_blank">Table S1</a>.</p><p>Highlighted in bold are those SNPs showing strongest associations with breast or ovarian cancer risk (p<0.01).</p

Mapping of the 17q21 locus.

<p><i>Top 3 panels:</i> P-values of association (−log₁₀ scale) with ovarian cancer risk for genotyped and imputed SNPs (1000 Genomes Project CEU), by chromosome position (b.37) at the 17q21 region, for <i>BRCA1</i>, <i>BRCA2</i> mutation carriers and combined. Results based on the kinship-adjusted score test statistic (1 d.f.). <i>Fourth panel:</i> Genes in the region spanning (43.4–44.9 Mb, b.37) and the location of the most significant genotyped SNPs (in red font) and imputed SNPs (in black font). <i>Bottom panel:</i> Pairwise r² values for genotyped SNPs on iCOG array in the 17q21 region covering positions (43.4–44.9 Mb, b.37).</p

Predicted breast and ovarian cancer absolute risks for <i>BRCA1</i> mutation carriers at the 5^th, 10^th, 90^th, and 95^th percentiles of the combined SNP profile distributions.

<p>The minimum, maximum and average risks are also shown. Predicted cancer risks are based on the associations of known breast or ovarian cancer susceptibility loci (identified through GWAS) with cancer risk for <i>BRCA1</i> mutation carriers and loci identified through the present study. Breast cancer risks based on the associations with: 1q32, 10q25.3, 19p13, 6q25.1, 12p11, <i>TOX3</i>, 2q35, <i>LSP1</i>, <i>RAD51L1</i> (based on HR and minor allele frequency estimates from <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003212#pgen-1003212-t001" target="_blank">Table 1</a>, <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003212#pgen-1003212-t002" target="_blank">Table 2</a>, and <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003212#pgen.1003212.s016" target="_blank">Table S4</a>) and <i>TERT </i><a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003212#pgen.1003212-Bojesen1" target="_blank">[31]</a>. Ovarian cancer risks based on the associations with: 9p22, 8q24, 3q25, 17q21, 19p13 (<a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003212#pgen-1003212-t001" target="_blank">Table 1</a>) and 17q21.31, 4q32.3 (<a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003212#pgen-1003212-t002" target="_blank">Table 2</a>). Only the top SNP from each region was chosen. Average breast and ovarian cancer risks were obtained from published data <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003212#pgen.1003212-Antoniou10" target="_blank">[25]</a>. The methods for calculating the predicted risks have been described previously <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003212#pgen.1003212-Antoniou11" target="_blank">[28]</a>.</p

Study design for selection of the SNPs and genotyping of <i>BRCA1</i> samples.

<p>GWAS data from 2,727 <i>BRCA1</i> mutation carriers were analysed for associations with breast and ovarian cancer risk and 32,557 SNPs were selected for inclusion on the iCOGS array. A total of 11,705 <i>BRCA1</i> samples (after quality control (QC) checks) were genotyped on the 31,812 <i>BRCA1</i>-GWAS SNPs from the iCOGS array that passed QC. Of these samples, 2,387 had been genotyped at the SNP selection stage and are referred to as “stage 1” samples, whereas 9,318 samples were unique to the iCOGS study (“Stage 2” samples). Next, 17 SNPs that exhibited the most significant associations with breast and ovarian cancer were selected for genotyping in a third stage involving an additional 2,646 <i>BRCA1</i> samples (after QC).</p

Associations with SNPs at the novel 17q21 region with ovarian cancer risk for <i>BRCA1</i> and <i>BRCA2</i> mutation carriers.

*<p>HRs estimated under the single disease risk model.</p

Associations with breast or ovarian cancer risk for loci previously reported to be associated with cancer risk for <i>BRCA1</i> mutation carriers.

<p>Freq = frequency of allele 2 in unaffected <i>BRCA1</i> carriers.</p><p>HR = Per allele Hazard Ratio associated with allele 2, under a single disease risk model, unless specified.</p><p>r²: correlation between the SNP in the present study and the published SNP.</p>*<p>SNP not in <i>BRCA1</i> GWAS SNP allocation on iCOGS chip.</p>a<p>: rs13387042 was previously found to be associated only under the 2-df model.</p>b<p>: analysis under a competing risks model.</p