Repeatability of methylation measures using a QIAseq targeted methyl panel and comparison with the Illumina HumanMethylation450 assay

OBJECTIVE

In previous studies using Illumina Infinium methylation arrays, we have identified DNA methylation marks associated with cancer predisposition and progression. In the present study, we have sought to find appropriate technology to both technically validate our data and expand our understanding of DNA methylation in these genomic regions. Here, we aimed to assess the repeatability of methylation measures made using QIAseq targeted methyl panel and to compare them with those obtained from the Illumina HumanMethylation450 (HM450K) assay. We included in the analysis high molecular weight DNA extracted from whole blood (WB) and DNA extracted from formalin-fixed paraffin-embedded tissues (FFPE).

RESULTS

The repeatability of QIAseq-methylation measures was assessed at 40 CpGs, using the Intraclass Correlation Coefficient (ICC). The mean ICCs and 95% confidence intervals (CI) were 0.72 (0.62-0.81), 0.59 (0.47-0.71) and 0.80 (0.73-0.88) for WB, FFPE and both sample types combined, respectively. For technical replicates measured using QIAseq and HM450K, the mean ICCs (95% CI) were 0.53 (0.39-0.68), 0.43 (0.31-0.56) and 0.70 (0.59-0.80), respectively. Bland-Altman plots indicated good agreement between QIAseq and HM450K measurements. These results demonstrate that the QIAseq targeted methyl panel produces reliable and reproducible methylation measurements across the 40 CpGs that were examined.

Mismatch repair gene pathogenic germline variants in a population-based cohort of breast cancer

The advent of gene panel testing is challenging the previous practice of using clinically defined cancer family syndromes to inform single-gene genetic screening. Individual and family cancer histories that would have previously indicated testing of a single gene or a small number of related genes are now, increasingly, leading to screening across gene panels that contain larger numbers of genes. We have applied a gene panel test that included four DNA mismatch repair (MMR) genes (MLH1, MSH2, MSH6 and PMS2) to an Australian population-based case-control-family study of breast cancer. Altogether, eight pathogenic variants in MMR genes were identified: six in 1421 case-families (0.4%, 4 MSH6 and 2 PMS2) and two in 833 control-families (0.2%, one each of MLH1 and MSH2). This testing highlights the current and future challenges for clinical genetics in the context of anticipated gene panel-based population-based screening that includes the MMR genes. This testing is likely to provide additional opportunities for cancer prevention via cascade testing for Lynch syndrome and precision medicine for breast cancer treatment.

Population-based estimates of the age-specific cumulative risk of breast cancer for pathogenic variants in CHEK2: Findings from the Australian Breast Cancer Family Registry

551

Background: Case-control studies of breast cancer have consistently shown that pathogenic variants in CHEK2 are associated with about a 3-fold increased risk of breast cancer. Information about the recurrent protein truncating variant CHEK2c.1100delC dominates this estimate. There have been no formal estimates of age-specific cumulative risk of breast cancer for all CHEK2 pathogenic (including likely pathogenic) variants combined. Methods: We conducted a genetic screen of CHEK2 in an Australian population-based case-control-family study of breast cancer. This study is focused on disease at an early age and participants were unselected for family history. The age-specific cumulative risk (penetrance) of breast cancer was estimated using segregation analysis. Results: The estimated hazard ratio for carriers of pathogenic CHEK2 variants (combined) was 4.9 (95% CI 2.5-9.5; p < 0.0001) relative to non-carriers. The HR for carriers of the CHEK2 c.1100delC variant was estimated to be 3.5 (95% CI 1.02-11.6) and the HR for carriers of all other CHEK2 variants combined was estimated to be 5.7 (95% CI 2.5-12.9). The age-specific cumulative risk of breast cancer was estimated to be 18% (95% CI 11-30%) and 33% (95% CI 21-48%) to age 60 and 80 years, respectively. Conclusions: These findings provide important information for the clinical management of breast cancer risk for women carrying pathogenic variants in CHEK2.

Abstract PS7-04: Population-based estimates of breast cancer risk for germline pathogenic variants identified by gene-panel testing: An Australian perspective

BRA-STRAP is an Australia-wide study of breast cancer predisposition that brings together gene-panel data from 30,000 adult Australian women of all ages, across the breast cancer risk spectrum, with and without a diagnosis of breast cancer. The “BRA-STRAP panel” includes 24 genes* that are involved in, or putatively associated with, predisposition to breast and/or ovarian cancer. Despite insufficient evidence for clinical translation for some of these genes, all 24 are commonly included on panel tests for breast cancer predisposition.

We present findings from the population-based case-control sub-study of BRA-STRAP, which involved 1451 women diagnosed with breast cancer and 857 age-matched controls participating in the Australian Breast Cancer Family Registry (ABCFR), and 6101 healthy, elderly Australian women enrolled in the ASPREE study. These analyses focus on rare genetic variants predicted to lead to loss of function and/or classified as pathogenic/likely pathogenic (P/LP) in ClinVar. Odds ratios (ORs) for their associations with breast cancer were estimated by aggregating genetic variants for each gene.

For the women diagnosed with breast cancer, the median age at diagnosis (inter-quartile range, IQR) was 40.0 (14.0) years and the overall frequency of P/LP variant carriers across all genes was 156/1451 (10.8%). The median age (IQR) of the ABCFR and ASPREE controls were 39.4 (14.9) and 73.9 (5.8) years, respectively. The frequencies of P/LP variant carriers were 33/857 (3.9%) and 268/6101 (4.4%) in the ABCFR and ASPREE controls, respectively. We combined both control datasets and, after adjusting for age and other potential confounders, the ORs associated with P/LP variants in BRCA1 and BRCA2 were 4.1 [95% confidence interval (CI): 1.8-10.2] and 2.9 [95% CI: 1.5-6], respectively. We also found that the OR for P/LP variants in ATM was 4.0 [95% CI: 1.5-10.4] and the OR for P/LP variants in PALB2 was 2.2 [95% CI: 0.75-5.7] although this did not reach statistical significance.

These results contribute to international efforts to refine the breast cancer risk estimates for genetic variants identified from population-based screening of unselected women using genes that are included on panel tests and thought to be potentially breast cancer predisposition genes.The case-control-family design of the ABCFR will also allow us to estimate the age specific cumulative risk (penetrance) of these genetic variants, which is important for genetic counselling and the clinical management of carrier families.

*ATM, BARD1, BRCA1, BRCA2, BRIP1, CDH1, CHEK2, FANCM, MLH1, MRE11A, MSH2, MSH6, MUTYH, NBN, NF1, PALB2, PMS2, PTEN, RAD50, RAD51C, RAD51D, RECQL, STK11 and TP53

Citation Format: Tu Nguyen-Dumont, James Dowty, Katherine Tucker, Judy Kirk, Paul James, Alison Trainer, Ingrid Winship, Nicholas Pachter, Nicola Poplawski, Scott Grist, Daniel J Park, Anne-Laure Renault, Fleur Hammet, Maryam Mahmoodi, Helen Tsimiklis, Jason A Steen, Derrick Theys, Amanda Rewse, Amanda Willis, April Morrow, Catherine Speechly, Rebecca Harris, Moeen Riaz, Robert Sebra, Eric Schadt, Paul Lacaze, John McNeil, John L Hopper, Melissa C Southey. Population-based estimates of breast cancer risk for germline pathogenic variants identified by gene-panel testing: An Australian perspective [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PS7-04.

Rare germline genetic variants and risk of aggressive prostate cancer

Few genetic risk factors have been demonstrated to be specifically associated with aggressive prostate cancer (PrCa). Here, we report a case-case study of PrCa comparing the prevalence of germline pathogenic/likely pathogenic (P/LP) genetic variants in 787 men with aggressive disease and 769 with nonaggressive disease. Overall, we observed P/LP variants in 11.4% of men with aggressive PrCa and 9.8% of men with nonaggressive PrCa (two-tailed Fisher's exact tests, P = .28). The proportion of BRCA2 and ATM P/LP variant carriers in men with aggressive PrCa exceeded that observed in men with nonaggressive PrCa; 18/787 carriers (2.3%) and 4/769 carriers (0.5%), P = .004, and 14/787 carriers (0.02%) and 5/769 carriers (0.01%), P = .06, respectively. Our findings contribute to the extensive international effort to interpret the genetic variation identified in genes included on gene-panel tests, for which there is currently an insufficient evidence-base for clinical translation in the context of PrCa risk.

Hi-Plex for Simple, Accurate, and Cost-Effective Amplicon-based Targeted DNA Sequencing

Hi-Plex is a suite of methods to enable simple, accurate, and cost-effective highly multiplex PCR-based targeted sequencing (Nguyen-Dumont et al., Biotechniques 58:33-36, 2015). At its core is the principle of using gene-specific primers (GSPs) to "seed" (or target) the reaction and universal primers to "drive" the majority of the reaction. In this manner, effects on amplification efficiencies across the target amplicons can, to a large extent, be restricted to early seeding cycles. Product sizes are defined within a relatively narrow range to enable high-specificity size selection, replication uniformity across target sites (including in the context of fragmented input DNA such as that derived from fixed tumor specimens (Nguyen-Dumont et al., Biotechniques 55:69-74, 2013; Nguyen-Dumont et al., Anal Biochem 470:48-51, 2015), and application of high-specificity genetic variant calling algorithms (Pope et al., Source Code Biol Med 9:3, 2014; Park et al., BMC Bioinformatics 17:165, 2016). Hi-Plex offers a streamlined workflow that is suitable for testing large numbers of specimens without the need for automation.

Is RNASEL:p.Glu265* a modifier of early-onset breast cancer risk for carriers of high-risk mutations?

Background

Breast cancer risk for BRCA1 and BRCA2 pathogenic mutation carriers is modified by risk factors that cluster in families, including genetic modifiers of risk. We considered genetic modifiers of risk for carriers of high-risk mutations in other breast cancer susceptibility genes.

Methods

In a family known to carry the high-risk mutation PALB2:c.3113G>A (p.Trp1038*), whole-exome sequencing was performed on germline DNA from four affected women, three of whom were mutation carriers.

Results

RNASEL:p.Glu265* was identified in one of the PALB2 carriers who had two primary invasive breast cancer diagnoses before 50 years. Gene-panel testing of BRCA1, BRCA2, PALB2 and RNASEL in the Australian Breast Cancer Family Registry identified five carriers of RNASEL:p.Glu265* in 591 early onset breast cancer cases. Three of the five women (60%) carrying RNASEL:p.Glu265* also carried a pathogenic mutation in a breast cancer susceptibility gene compared with 30 carriers of pathogenic mutations in the 586 non-carriers of RNASEL:p.Glu265* (5%) (p < 0.002). Taqman genotyping demonstrated that the allele frequency of RNASEL:p.Glu265* was similar in affected and unaffected Australian women, consistent with other populations.

Conclusion

Our study suggests that RNASEL:p.Glu265* may be a genetic modifier of risk for early-onset breast cancer predisposition in carriers of high-risk mutations. Much larger case-case and case-control studies are warranted to test the association observed in this report.

FANCM and RECQL genetic variants and breast cancer susceptibility: relevance to South Poland and West Ukraine

Background

FANCM and RECQL have recently been reported as breast cancer susceptibility genes and it has been suggested that they should be included on gene panel tests for breast cancer predisposition. However, the clinical value of testing for mutations in RECQL and FANCM remains to be determined. In this study, we have characterised the spectrum of FANCM and RECQL mutations in women affected with breast or ovarian cancer from South-West Poland and West Ukraine.

Methods

We applied Hi-Plex, an amplicon-based enrichment method for targeted massively parallel sequencing, to screen the coding exons and proximal intron-exon junctions of FANCM and RECQL in germline DNA from unrelated women affected with breast cancer (n = 338) and ovarian cancer (n = 89) from Poland (n = 304) and Ukraine (n = 123). These women were at high-risk of carrying a genetic predisposition to breast and/or ovarian cancer due to a family history and/or early-onset disease.

Results

Among 427 women screened, we identified one carrier of the FANCM:c.1972C > T nonsense mutation (0.23%), and two carriers of the frameshift insertion FANCM:c.1491dup (0.47%). None of the variants we observed in RECQL were predicted to be loss-of-function mutations by standard variant effect prediction tools.

Conclusions

Our study of the Polish and Ukrainian populations has identified a carrier frequency of truncating mutations in FANCM consistent with previous reports. Although initial reports suggesting that mutations in RECQL could be associated with increased breast cancer risk included women from Poland and identified the RECQL:c.1667_1667 + 3delAGTA mutation in 0.23–0.35% of breast cancer cases, we did not observe any carriers in our study cohort. Continued screening, both in research and diagnostic settings, will enable the accumulation of data that is needed to establish the clinical utility of including RECQL and FANCM on gene panel tests.

Electronic supplementary material

The online version of this article (10.1186/s12881-018-0524-x) contains supplementary material, which is available to authorized users.

A high-plex PCR approach for massively parallel sequencing

Current methods for targeted massively parallel sequencing (MPS) have several drawbacks, including limited design flexibility, expense, and protocol complexity, which restrict their application to settings involving modest target size and requiring low cost and high throughput. To address this, we have developed Hi-Plex, a PCR-MPS strategy intended for high-throughput screening of multiple genomic target regions that integrates simple, automated primer design software to control product size. Featuring permissive thermocycling conditions and clamp bias reduction, our protocol is simple, cost- and time-effective, uses readily available reagents, does not require expensive instrumentation, and requires minimal optimization. In a 60-plex assay targeting the breast cancer predisposition genes PALB2 and XRCC2, we applied Hi-Plex to 100 ng LCL-derived DNA, and 100 ng and 25 ng FFPE tumor-derived DNA. Altogether, at least 86.94% of the human genome-mapped reads were on target, and 100% of targeted amplicons were represented within 25-fold of the mean. Using 25 ng FFPE-derived DNA, 95.14% of mapped reads were on-target and relative representation ranged from 10.1-fold lower to 5.8-fold higher than the mean. These results were obtained using only the initial automatically-designed primers present in equal concentration. Hi-Plex represents a powerful new approach for screening panels of genomic target regions.

Mutation screening of PALB2 in clinically ascertained families from the Breast Cancer Family Registry

Loss-of-function mutations in PALB2 are associated with an increased risk of breast cancer, with recent data showing that female breast cancer risks for PALB2 mutation carriers are comparable in magnitude to those for BRCA2 mutation carriers. This study applied targeted massively parallel sequencing to characterize the mutation spectrum of PALB2 in probands attending breast cancer genetics clinics in the USA. The coding regions and proximal intron-exon junctions of PALB2 were screened in probands not known to carry a mutation in BRCA1 or BCRA2 from 1,250 families enrolled through familial cancer clinics by the Breast Cancer Family Registry. Mutation screening was performed using Hi-Plex, an amplicon-based targeted massively parallel sequencing platform. Screening of PALB2 was successful in 1,240/1,250 probands and identified nine women with protein-truncating mutations (three nonsense mutations and five frameshift mutations). Four of the 33 missense variants were predicted to be deleterious to protein function by in silico analysis using two different programs. Analysis of tumors from carriers of truncating mutations revealed that the majority were high histological grade, invasive ductal carcinomas. Young onset was apparent in most families, with 19 breast cancers under 50 years of age, including eight under the age of 40 years. Our data demonstrate the utility of Hi-Plex in the context of high-throughput testing for rare genetic mutations and provide additional timely information about the nature and prevalence of PALB2 mutations, to enhance risk assessment and risk management of women at high risk of cancer attending clinical genetic services.

Hi-Plex targeted sequencing is effective using DNA derived from archival dried blood spots

Many genetic epidemiology resources have collected dried blood spots (predominantly as Guthrie Cards) as an economical and efficient means of archiving sources of DNA, conferring great value to genetic screening methods that are compatible with this medium. We applied Hi-Plex to screen the breast cancer predisposition gene PALB2 in 93 Guthrie Card-derived DNA specimens previously characterized for PALB2 genetic variants via DNA derived from lymphoblastoid cell lines, whole blood, and buffy coat. Of the 93 archival Guthrie Card-derived DNAs, 92 (99%) were processed successfully and sequenced using approximately half of a MiSeq run. From these 92 DNAs, all 59 known variants were detected and no false-positive variant calls were yielded. Fully 98.13% of amplicons (5417/5520) were represented within 15-fold of the median coverage (2786 reads), and 99.98% of amplicons (5519/5520) were represented at a depth of 10 read-pairs or greater. With Hi-Plex, we show for the first time that a High-Plex amplicon-based massively parallel sequencing (MPS) system can be applied effectively to DNA prepared from dried blood spot archival specimens and, as such, can dramatically increase the scopes of both method and resource.

Rare mutations in RINT1 predispose carriers to breast and Lynch syndrome-spectrum cancers

Approximately half of the familial aggregation of breast cancer remains unexplained. A multiple-case breast cancer family exome-sequencing study identified three likely pathogenic mutations in RINT1 (NM_021930.4) not present in public sequencing databases: RINT1 c.343C>T (p.Q115X), c.1132_1134del (p.M378del), and c.1207G>T (p.D403Y). On the basis of this finding, a population-based case-control mutation-screening study was conducted that identified 29 carriers of rare (minor allele frequency < 0.5%), likely pathogenic variants: 23 in 1,313 early-onset breast cancer cases and six in 1,123 frequency-matched controls [OR, 3.24; 95% confidence interval (CI), 1.29-8.17; P = 0.013]. RINT1 mutation screening of probands from 798 multiple-case breast cancer families identified four additional carriers of rare genetic variants. Analysis of the incidence of first primary cancers in families of women carrying RINT1 mutations estimated that carriers were at increased risk of Lynch syndrome-spectrum cancers [standardized incidence ratio (SIR), 3.35; 95% CI, 1.7-6.0; P = 0.005], particularly for relatives diagnosed with cancer under the age of 60 years (SIR, 10.9; 95% CI, 4.7-21; P = 0.0003).

SIGNIFICANCE

The work described in this study adds RINT1 to the growing list of genes in which rare sequence variants are associated with intermediate levels of breast cancer risk. Given that RINT1 is also associated with a spectrum of cancers with mismatch repair defects, these findings have clinical applications and raise interesting biological questions.

ROVER variant caller: read-pair overlap considerate variant-calling software applied to PCR-based massively parallel sequencing datasets

Background

We recently described Hi-Plex, a highly multiplexed PCR-based target-enrichment system for massively parallel sequencing (MPS), which allows the uniform definition of library size so that subsequent paired-end sequencing can achieve complete overlap of read pairs. Variant calling from Hi-Plex-derived datasets can thus rely on the identification of variants appearing in both reads of read-pairs, permitting stringent filtering of sequencing chemistry-induced errors. These principles underly ROVER software (derived from Read Overlap PCR-MPS variant caller), which we have recently used to report the screening for genetic mutations in the breast cancer predisposition gene PALB2. Here, we describe the algorithms underlying ROVER and its usage.

Results

ROVER enables users to quickly and accurately identify genetic variants from PCR-targeted, overlapping paired-end MPS datasets. The open-source availability of the software and threshold tailorability enables broad access for a range of PCR-MPS users.

Methods

ROVER is implemented in Python and runs on all popular POSIX-like operating systems (Linux, OS X). The software accepts a tab-delimited text file listing the coordinates of the target-specific primers used for targeted enrichment based on a specified genome-build. It also accepts aligned sequence files resulting from mapping to the same genome-build. ROVER identifies the amplicon a given read-pair represents and removes the primer sequences by using the mapping co-ordinates and primer co-ordinates. It considers overlapping read-pairs with respect to primer-intervening sequence. Only when a variant is observed in both reads of a read-pair does the signal contribute to a tally of read-pairs containing or not containing the variant. A user-defined threshold informs the minimum number of, and proportion of, read-pairs a variant must be observed in for a ‘call’ to be made. ROVER also reports the depth of coverage across amplicons to facilitate the identification of any regions that may require further screening.

Conclusions

ROVER can facilitate rapid and accurate genetic variant calling for a broad range of PCR-MPS users.

Hi-Plex for high-throughput mutation screening: application to the breast cancer susceptibility gene PALB2

Background

Massively parallel sequencing (MPS) has revolutionised biomedical research and offers enormous capacity for clinical application. We previously reported Hi-Plex, a streamlined highly-multiplexed PCR-MPS approach, allowing a given library to be sequenced with both the Ion Torrent and TruSeq chemistries. Comparable sequencing efficiency was achieved using material derived from lymphoblastoid cell lines and formalin-fixed paraffin-embedded tumour.

Methods

Here, we report high-throughput application of Hi-Plex by performing blinded mutation screening of the coding regions of the breast cancer susceptibility gene PALB2 on a set of 95 blood-derived DNA samples that had previously been screened using Sanger sequencing and high-resolution melting curve analysis (n = 90), or genotyped by Taqman probe-based assays (n = 5). Hi-Plex libraries were prepared simultaneously using relatively inexpensive, readily available reagents in a simple half-day protocol followed by MPS on a single MiSeq run.

Results

We observed that 99.93% of amplicons were represented at ≥10X coverage. All 56 previously identified variant calls were detected and no false positive calls were assigned. Four additional variant calls were made and confirmed upon re-analysis of previous data or subsequent Sanger sequencing.

Conclusions

These results support Hi-Plex as a powerful approach for rapid, cost-effective and accurate high-throughput mutation screening. They further demonstrate that Hi-Plex methods are suitable for and can meet the demands of high-throughput genetic testing in research and clinical settings.

Cross-platform compatibility of Hi-Plex, a streamlined approach for targeted massively parallel sequencing

Although per-base sequencing costs have decreased during recent years, library preparation for targeted massively parallel sequencing remains constrained by high reagent cost, limited design flexibility, and protocol complexity. To address these limitations, we previously developed Hi-Plex, a polymerase chain reaction (PCR) massively parallel sequencing strategy for screening panels of genomic target regions. Here, we demonstrate that Hi-Plex applied with hybrid adapters can generate a library suitable for sequencing with both the Ion Torrent and the TruSeq chemistries and that adjusting primer concentrations improves coverage uniformity. These results expand Hi-Plex capabilities as an accurate, affordable, flexible, and rapid approach for various genetic screening applications.

Using SNP genotypes to improve the discrimination of a simple breast cancer risk prediction model

It has been shown that, for women aged 50 years or older, the discriminatory accuracy of the Breast Cancer Risk Prediction Tool (BCRAT) can be modestly improved by the inclusion of information on common single nucleotide polymorphisms (SNPs) that are associated with increased breast cancer risk. We aimed to determine whether a similar improvement is seen for earlier onset disease. We used the Australian Breast Cancer Family Registry to study a population-based sample of 962 cases aged 35-59 years, and 463 controls frequency matched for age and for whom genotyping data was available. Overall, the inclusion of data on seven SNPs improved the area under the receiver operating characteristic curve (AUC) from 0.58 (95 % confidence interval [CI] 0.55-0.61) for BCRAT alone to 0.61 (95 % CI 0.58-0.64) for BCRAT and SNP data combined (p < 0.001). For women aged 35-39 years at interview, the corresponding improvement in AUC was from 0.61 (95 % CI 0.56-0.66) to 0.65 (95 % CI 0.60-0.70; p = 0.03), while for women aged 40-49 years at diagnosis, the AUC improved from 0.61 (95 % CI 0.55-0.66) to 0.63 (95 % CI 0.57-0.69; p = 0.04). Using previously used classifications of low, intermediate and high risk, 2.1 % of cases and none of the controls aged 35-39 years, and 10.9 % of cases and 4.0 % of controls aged 40-49 years were classified into a higher risk group. Including information on seven SNPs associated with breast cancer risk, improves the discriminatory accuracy of BCRAT for women aged 35-39 years and 40-49 years. Given, the low absolute risk for women in these age groups, only a small proportion are reclassified into a higher category for predicted 5-year risk of breast cancer.

Genome-wide association studies identify four ER negative-specific breast cancer risk loci

Estrogen receptor (ER)-negative tumors represent 20-30% of all breast cancers, with a higher proportion occurring in younger women and women of African ancestry. The etiology and clinical behavior of ER-negative tumors are different from those of tumors expressing ER (ER positive), including differences in genetic predisposition. To identify susceptibility loci specific to ER-negative disease, we combined in a meta-analysis 3 genome-wide association studies of 4,193 ER-negative breast cancer cases and 35,194 controls with a series of 40 follow-up studies (6,514 cases and 41,455 controls), genotyped using a custom Illumina array, iCOGS, developed by the Collaborative Oncological Gene-environment Study (COGS). SNPs at four loci, 1q32.1 (MDM4, P = 2.1 × 10(-12) and LGR6, P = 1.4 × 10(-8)), 2p24.1 (P = 4.6 × 10(-8)) and 16q12.2 (FTO, P = 4.0 × 10(-8)), were associated with ER-negative but not ER-positive breast cancer (P > 0.05). These findings provide further evidence for distinct etiological pathways associated with invasive ER-positive and ER-negative breast cancers.

FAVR (Filtering and Annotation of Variants that are Rare): methods to facilitate the analysis of rare germline genetic variants from massively parallel sequencing datasets

Background

Characterising genetic diversity through the analysis of massively parallel sequencing (MPS) data offers enormous potential to significantly improve our understanding of the genetic basis for observed phenotypes, including predisposition to and progression of complex human disease. Great challenges remain in resolving genetic variants that are genuine from the millions of artefactual signals.

Results

FAVR is a suite of new methods designed to work with commonly used MPS analysis pipelines to assist in the resolution of some of the issues related to the analysis of the vast amount of resulting data, with a focus on relatively rare genetic variants. To the best of our knowledge, no equivalent method has previously been described. The most important and novel aspect of FAVR is the use of signatures in comparator sequence alignment files during variant filtering, and annotation of variants potentially shared between individuals. The FAVR methods use these signatures to facilitate filtering of (i) platform and/or mapping-specific artefacts, (ii) common genetic variants, and, where relevant, (iii) artefacts derived from imbalanced paired-end sequencing, as well as annotation of genetic variants based on evidence of co-occurrence in individuals. We applied conventional variant calling applied to whole-exome sequencing datasets, produced using both SOLiD and TruSeq chemistries, with or without downstream processing by FAVR methods. We demonstrate a 3-fold smaller rare single nucleotide variant shortlist with no detected reduction in sensitivity. This analysis included Sanger sequencing of rare variant signals not evident in dbSNP131, assessment of known variant signal preservation, and comparison of observed and expected rare variant numbers across a range of first cousin pairs. The principles described herein were applied in our recent publication identifying XRCC2 as a new breast cancer risk gene and have been made publically available as a suite of software tools.

Conclusions

FAVR is a platform-agnostic suite of methods that significantly enhances the analysis of large volumes of sequencing data for the study of rare genetic variants and their influence on phenotypes.

Breast cancer risk and 6q22.33: combined results from Breast Cancer Association Consortium and Consortium of Investigators on Modifiers of BRCA1/2

Recently, a locus on chromosome 6q22.33 (rs2180341) was reported to be associated with increased breast cancer risk in the Ashkenazi Jewish (AJ) population, and this association was also observed in populations of non-AJ European ancestry. In the present study, we performed a large replication analysis of rs2180341 using data from 31,428 invasive breast cancer cases and 34,700 controls collected from 25 studies in the Breast Cancer Association Consortium (BCAC). In addition, we evaluated whether rs2180341 modifies breast cancer risk in 3,361 BRCA1 and 2,020 BRCA2 carriers from 11 centers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Based on the BCAC data from women of European ancestry, we found evidence for a weak association with breast cancer risk for rs2180341 (per-allele odds ratio (OR) = 1.03, 95% CI 1.00-1.06, p = 0.023). There was evidence for heterogeneity in the ORs among studies (I(2) = 49.3%; p = <0.004). In CIMBA, we observed an inverse association with the minor allele of rs2180341 and breast cancer risk in BRCA1 mutation carriers (per-allele OR = 0.89, 95%CI 0.80-1.00, p = 0.048), indicating a potential protective effect of this allele. These data suggest that that 6q22.33 confers a weak effect on breast cancer risk.

Associations of common variants at 1p11.2 and 14q24.1 (RAD51L1) with breast cancer risk and heterogeneity by tumor subtype: findings from the Breast Cancer Association Consortium

A genome-wide association study (GWAS) identified single-nucleotide polymorphisms (SNPs) at 1p11.2 and 14q24.1 (RAD51L1) as breast cancer susceptibility loci. The initial GWAS suggested stronger effects for both loci for estrogen receptor (ER)-positive tumors. Using data from the Breast Cancer Association Consortium (BCAC), we sought to determine whether risks differ by ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2), grade, node status, tumor size, and ductal or lobular morphology. We genotyped rs11249433 at 1p.11.2, and two highly correlated SNPs rs999737 and rs10483813 (r(2)= 0.98) at 14q24.1 (RAD51L1), for up to 46 036 invasive breast cancer cases and 46 930 controls from 39 studies. Analyses by tumor characteristics focused on subjects reporting to be white women of European ancestry and were based on 25 458 cases, of which 87% had ER data. The SNP at 1p11.2 showed significantly stronger associations with ER-positive tumors [per-allele odds ratio (OR) for ER-positive tumors was 1.13, 95% CI = 1.10-1.16 and, for ER-negative tumors, OR was 1.03, 95% CI = 0.98-1.07, case-only P-heterogeneity = 7.6 × 10(-5)]. The association with ER-positive tumors was stronger for tumors of lower grade (case-only P= 6.7 × 10(-3)) and lobular histology (case-only P= 0.01). SNPs at 14q24.1 were associated with risk for most tumor subtypes evaluated, including triple-negative breast cancers, which has not been described previously. Our results underscore the need for large pooling efforts with tumor pathology data to help refine risk estimates for SNP associations with susceptibility to different subtypes of breast cancer.

Low penetrance breast cancer susceptibility loci are associated with specific breast tumor subtypes: findings from the Breast Cancer Association Consortium

Breast cancers demonstrate substantial biological, clinical and etiological heterogeneity. We investigated breast cancer risk associations of eight susceptibility loci identified in GWAS and two putative susceptibility loci in candidate genes in relation to specific breast tumor subtypes. Subtypes were defined by five markers (ER, PR, HER2, CK5/6, EGFR) and other pathological and clinical features. Analyses included up to 30 040 invasive breast cancer cases and 53 692 controls from 31 studies within the Breast Cancer Association Consortium. We confirmed previous reports of stronger associations with ER+ than ER- tumors for six of the eight loci identified in GWAS: rs2981582 (10q26) (P-heterogeneity = 6.1 × 10(-18)), rs3803662 (16q12) (P = 3.7 × 10(-5)), rs13281615 (8q24) (P = 0.002), rs13387042 (2q35) (P = 0.006), rs4973768 (3p24) (P = 0.003) and rs6504950 (17q23) (P = 0.002). The two candidate loci, CASP8 (rs1045485, rs17468277) and TGFB1 (rs1982073), were most strongly related with the risk of PR negative tumors (P = 5.1 × 10(-6) and P = 4.1 × 10(-4), respectively), as previously suggested. Four of the eight loci identified in GWAS were associated with triple negative tumors (P ≤ 0.016): rs3803662 (16q12), rs889312 (5q11), rs3817198 (11p15) and rs13387042 (2q35); however, only two of them (16q12 and 2q35) were associated with tumors with the core basal phenotype (P ≤ 0.002). These analyses are consistent with different biological origins of breast cancers, and indicate that tumor stratification might help in the identification and characterization of novel risk factors for breast cancer subtypes. This may eventually result in further improvements in prevention, early detection and treatment.

Associations of breast cancer risk factors with tumor subtypes: a pooled analysis from the Breast Cancer Association Consortium studies

BACKGROUND

Previous studies have suggested that breast cancer risk factors are associated with estrogen receptor (ER) and progesterone receptor (PR) expression status of the tumors.

METHODS

We pooled tumor marker and epidemiological risk factor data from 35,568 invasive breast cancer case patients from 34 studies participating in the Breast Cancer Association Consortium. Logistic regression models were used in case-case analyses to estimate associations between epidemiological risk factors and tumor subtypes, and case-control analyses to estimate associations between epidemiological risk factors and the risk of developing specific tumor subtypes in 12 population-based studies. All statistical tests were two-sided.

RESULTS

In case-case analyses, of the epidemiological risk factors examined, early age at menarche (≤12 years) was less frequent in case patients with PR(-) than PR(+) tumors (P = .001). Nulliparity (P = 3 × 10(-6)) and increasing age at first birth (P = 2 × 10(-9)) were less frequent in ER(-) than in ER(+) tumors. Obesity (body mass index [BMI] ≥ 30 kg/m(2)) in younger women (≤50 years) was more frequent in ER(-)/PR(-) than in ER(+)/PR(+) tumors (P = 1 × 10(-7)), whereas obesity in older women (>50 years) was less frequent in PR(-) than in PR(+) tumors (P = 6 × 10(-4)). The triple-negative (ER(-)/PR(-)/HER2(-)) or core basal phenotype (CBP; triple-negative and cytokeratins [CK]5/6(+) and/or epidermal growth factor receptor [EGFR](+)) accounted for much of the heterogeneity in parity-related variables and BMI in younger women. Case-control analyses showed that nulliparity, increasing age at first birth, and obesity in younger women showed the expected associations with the risk of ER(+) or PR(+) tumors but not triple-negative (nulliparity vs parity, odds ratio [OR] = 0.94, 95% confidence interval [CI] = 0.75 to 1.19, P = .61; 5-year increase in age at first full-term birth, OR = 0.95, 95% CI = 0.86 to 1.05, P = .34; obesity in younger women, OR = 1.36, 95% CI = 0.95 to 1.94, P = .09) or CBP tumors.

CONCLUSIONS

This study shows that reproductive factors and BMI are most clearly associated with hormone receptor-positive tumors and suggest that triple-negative or CBP tumors may have distinct etiology.

Candidate tumor-suppressor genes on chromosome arm 8p in early-onset and high-grade breast cancers

Loss of genetic material from chromosome arm 8p occurs commonly in breast carcinomas, suggesting that this region is the site of one or more tumor-suppressor genes (TSGs). Comparative genomic hybridization analysis showed that 8p loss is more common in breast cancers from pre-menopausal compared with post-menopausal patients, as well as in high-grade breast cancers, regardless of the menopausal status. Subsequent high-resolution gene expression profiling of genes mapped to chromosome arm 8p, on an extended cohort of clinical tumor samples, indicated a similar dichotomy of breast cancer clinicopathologic types. Some of these genes showed differential downregulation in early-onset and later-onset, high-grade cancers compared with lower-grade, later-onset cancers. Three such genes were analysed further by in situ technologies, performed on tissue microarrays representing breast tumor and normal tissue samples. PCM1, which encodes a centrosomal protein, and DUSP4/MKP-2, which encodes a MAP kinase phosphatase, both showed frequent gene and protein loss in carcinomas. In contrast, there was an excess of cases showing loss of expression in the absence of reduced gene copy number of SFRP1, which encodes a dominant-negative receptor for Wnt-family ligands. These candidate TSGs may constitute some of the molecular drivers of chromosome arm 8p loss in breast carcinogenesis.

Detailed gene copy number and RNA expression analysis of the 17q12-23 region in primary breast cancers

Chromosome region 17q12-23 commonly shows an increase in DNA copy number in breast cancers, suggesting that several oncogenes are located at this site. We performed a high-resolution expression array and comparative genomic hybridization analysis of genes mapped to the entire 17q12-23 region, to identify novel candidate oncogenes. We identified 24 genes that showed significant overexpression in breast cancers with gain of 17q12-23, compared to cancers without gain. These genes included previously identified oncogenes, together with several novel candidate oncogenes. FISH analysis using specific gene probes hybridized to tissue arrays confirmed the underlying amplification of overexpressed genes. This high-resolution analysis of the 17q12-23 region indicates that several established and novel candidate oncogenes, including a Wnt-signaling pathway member, are amplified and overexpressed within individual primary breast cancer samples. We were also able to confirm the presence of two apparently separate and reciprocally amplified groups of genes within this region. Investigation of these genes and their functional interactions will facilitate our understanding of breast oncogenesis and optimal management of this disease.

Genetic aberrations detected by comparative genomic hybridisation in vulvar cancers

Squamous cell carcinoma of the vulva is a disease of significant clinical importance, which arises in the presence or absence of human papillomavirus. We used comparative genomic hybridisation to document non-random chromosomal gains and losses within human papillomavirus positive and negative vulvar cancers. Gain of 3q was significantly more common in human papillomavirus-positive cancers compared to human papillomavirus-negative cancers. The smallest area of gain was 3q22-25, a chromosome region which is frequently gained in other human papillomavirus-related cancers. Chromosome 8q was more commonly gained in human papillomavirus-negative compared to human papillomavirus-positive cancers. 8q21 was the smallest region of gain, which has been identified in other, non-human papillomavirus-related cancers. Chromosome arms 3p and 11q were lost in both categories of vulvar cancer. This study has demonstrated chromosome locations important in the development of vulvar squamous cell carcinoma. Additionally, taken together with previous studies of human papillomavirus-positive cancers of other anogenital sites, the data indicate that one or more oncogenes important in the development and progression of human papillomavirus-induced carcinomas are located on 3q. The different genetic changes seen in human papillomavirus-positive and negative vulvar squamous cell carcinomas support the clinicopathological data indicating that these are different cancer types.

Distinct molecular pathogeneses of early-onset breast cancers in BRCA1 and BRCA2 mutation carriers: a population-based study

Breast cancers arising in women with and without a germline mutation in the BRCA1 or BRCA2 gene display different histological features, which suggests unique mechanisms of molecular pathogenesis: We used a molecular pathological analysis to define the genetic abnormalities relevant to these specific pathogeneses. Tumor material was studied from 40 women with breast cancer diagnosed before 40 years of age, sampled from a population-based study and stratified by BRCA1 and BRCA2 germline mutation status. Cases were not selected for family history or ethnic origin, and none were known to be genetically related. Thus, germline mutation itself is likely to impact on the molecular pathogenesis of these tumors, with no substantial influence due to modifying genetic or environmental factors. Breast cancers occurring in BRCA1 mutation carriers had significantly higher levels of p53 expression, including the preinvasive (carcinoma in situ) stage of disease, compared with cancers occurring in BRCA2 mutation carriers or women with no detectable germline mutation. These cancers also had a higher proliferation rate as measured by Ki-67 antibody. Expression of the prognostic factors c-erbB-2, cyclin D1, and estrogen receptor was significantly less common in BRCA1 mutation carriers. Lower levels of cyclin D1 were also found in cancers from BRCA2 mutation carriers compared with non-mutation carriers. Direct p53 mutation analysis revealed mutations in 18% of all of the early-onset breast cancers within the study and included rare insertion and deletional mutations in cancers from BRCA1 mutation carriers. Our data indicate that a BRCA1 breast cancer phenotype may be recognized by an exceptionally high proliferation rate and early and frequent p53 overexpression but infrequent selection for overexpression of several other prognostic factor proteins known to be involved in breast oncogenesis. In contrast, breast cancers arising in BRCA2 mutation carriers have a more heterogeneous phenotypic profile.

Chromosomal localization of the human P2y6 purinoceptor gene and phylogenetic analysis of the P2y purinoceptor family

The G-protein-coupled P2Y purinoceptors mediate a variety of physiological effects in response to extracellular nucleotides. With the recent discovery of several new members from a variety of species, the P2Y purinoceptor family now encompasses types P2Y1 to P2Y6. By fluorescence in situ hybridization and utilization of the National Center for Biotechnology Information (NCBI) database, the human P2Y6 gene was localized to chromosome 11q13.5, between polymorphic markers D11S1314 and D11S916. NCBI database analysis of the remaining human P2Y purinoceptor genes revealed that P2Y2 and P2Y6 mapped to within less than 4 cM, and thus constitute the first described chromosomal clustering of this gene family. Phylogenetic analysis of the P2Y purinoceptor family demonstrated the presence of five evolutionary branches and suggests the occurrence of an ancient gene duplication event.

Molecular cloning and sequencing of a novel human P2 nucleotide receptor

A novel human P2 nucleotide receptor has been cloned from a T-cell cDNA library. The predicted amino acid sequence shows characteristics of a G-protein-coupled receptor, and shares 88% homology with a recently characterised rat P2 nucleotide receptor sequence. Distinctive features include an extremely short cytoplasmic tail with only one putative protein kinase C phosphorylation site. Northern blot analysis revealed a 1.9 kb transcript expressed in the placenta.