Search

The coronavirus disease 2019 (COVID-19) pandemic has demonstrated a clear need for high-throughput, multiplexed and sensitive assays for detecting severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and other respiratory viruses and their emerging variants. Here, we present a cost-effective virus and variant detection platform, called microfluidic Combinatorial Arrayed Reactions for Multiplexed Evaluation of Nucleic acids (mCARMEN), which combines CRISPR-based diagnostics and microfluidics with a streamlined workflow for clinical use. We developed the mCARMEN respiratory virus panel to test for up to 21 viruses, including SARS-CoV-2, other coronaviruses and both influenza strains, and demonstrated its diagnostic-grade performance on 525 patient specimens in an academic setting and 166 specimens in a clinical setting. We further developed an mCARMEN panel to enable the identification of 6 SARS-CoV-2 variant lineages, including Delta and Omicron, and evaluated it on 2,088 patient specimens with near-perfect concordance to sequencing-based variant classification. Lastly, we implemented a combined Cas13 and Cas12 approach that enables quantitative measurement of SARS-CoV-2 and influenza A viral copies in samples. The mCARMEN platform enables high-throughput surveillance of multiple viruses and variants simultaneously, enabling rapid detection of SARS-CoV-2 variants.

Repeated emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants with increased fitness underscores the value of rapid detection and characterization of new lineages. We have developed PyR0, a hierarchical Bayesian multinomial logistic regression model that infers relative prevalence of all viral lineages across geographic regions, detects lineages increasing in prevalence, and identifies mutations relevant to fitness. Applying PyR0 to all publicly available SARS-CoV-2 genomes, we identify numerous substitutions that increase fitness, including previously identified spike mutations and many nonspike mutations within the nucleocapsid and nonstructural proteins. PyR0 forecasts growth of new lineages from their mutational profile, ranks the fitness of lineages as new sequences become available, and prioritizes mutations of biological and public health concern for functional characterization.

An outbreak of over 1,000 COVID-19 cases in Provincetown, Massachusetts (MA), in July 2021-the first large outbreak mostly in vaccinated individuals in the US-prompted a comprehensive public health response, motivating changes to national masking recommendations and raising questions about infection and transmission among vaccinated individuals. To address these questions, we combined viral genomic and epidemiological data from 467 individuals, including 40% of outbreak-associated cases. The Delta variant accounted for 99% of cases in this dataset; it was introduced from at least 40 sources, but 83% of cases derived from a single source, likely through transmission across multiple settings over a short time rather than a single event. Genomic and epidemiological data supported multiple transmissions of Delta from and between fully vaccinated individuals. However, despite its magnitude, the outbreak had limited onward impact in MA and the US overall, likely due to high vaccination rates and a robust public health response.

Population-based estimates of breast cancer risk for carriers of pathogenic variants identified by gene-panel testing are urgently required. Most prior research has been based on women selected for high-risk features and more data is needed to make inference about breast cancer risk for women unselected for family history, an important consideration of population screening. We tested 1464 women diagnosed with breast cancer and 862 age-matched controls participating in the Australian Breast Cancer Family Study (ABCFS), and 6549 healthy, older Australian women enroled in the ASPirin in Reducing Events in the Elderly (ASPREE) study for rare germline variants using a 24-gene-panel. Odds ratios (ORs) were estimated using unconditional logistic regression adjusted for age and other potential confounders. We identified pathogenic variants in 11.1% of the ABCFS cases, 3.7% of the ABCFS controls and 2.2% of the ASPREE (control) participants. The estimated breast cancer OR [95% confidence interval] was 5.3 [2.1-16.2] for BRCA1, 4.0 [1.9-9.1] for BRCA2, 3.4 [1.4-8.4] for ATM and 4.3 [1.0-17.0] for PALB2. Our findings provide a population-based perspective to gene-panel testing for breast cancer predisposition and opportunities to improve predictors for identifying women who carry pathogenic variants in breast cancer predisposition genes.

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.

PURPOSE: To characterize the spectrum of BRCA1 and BRCA2 pathogenic germline variants in women from south-west Poland and west Ukraine affected with breast or ovarian cancer. Testing in women at high risk of breast and ovarian cancer in these regions is currently mainly limited to founder mutations. METHODS: Unrelated women affected with breast and/or ovarian cancer from Poland (n = 337) and Ukraine (n = 123) were screened by targeted sequencing. Excluded from targeted sequencing were 34 Polish women who had previously been identified as carrying a founder mutation in BRCA1. No prior testing had been conducted among the Ukrainian women. Thus, this study screened BRCA1 and BRCA2 in the germline DNA of 426 women in total. RESULTS: We identified 31 and 18 women as carriers of pathogenic/likely pathogenic (P/LP) genetic variants in BRCA1 and BRCA2, respectively. We observed five BRCA1 and eight BRCA2 P/LP variants (13/337, 3.9%) in the Polish women. Combined with the 34/337 (10.1%) founder variants identified prior to this study, the overall P/LP variant frequency in the Polish women was thus 14% (47/337). Among the Ukrainian women, 16/123 (13%) women were identified as carrying a founder mutation and 20/123 (16.3%) were found to carry non-founder P/LP variants (10 in BRCA1 and 10 in BRCA2). CONCLUSIONS: These results indicate that genetic testing in women at high risk of breast and ovarian cancer in Poland and Ukraine should not be limited to founder mutations. Extended testing will enhance risk stratification and management for these women and their families.

DNA methylation influences predisposition, development and prognosis for many diseases, including cancer. However, it is not uncommon to encounter samples with incorrect sex labelling or atypical sex chromosome arrangement. Sex is one of the strongest influencers of the genomic distribution of DNA methylation and, therefore, correct assignment of sex and filtering of abnormal samples are essential for the quality control of study data. Differences in sex chromosome copy numbers between sexes and X-chromosome inactivation in females result in distinctive sex-specific patterns in the distribution of DNA methylation levels. In this study, we present a software tool, sEst, which incorporates clustering analysis to infer sex and to detect sex-chromosome abnormalities from DNA methylation microarray data. Testing with two publicly available datasets demonstrated that sEst not only correctly inferred the sex of the test samples, but also identified mislabelled samples and samples with potential sex-chromosome abnormalities, such as Klinefelter syndrome and Turner syndrome, the latter being a feature not offered by existing methods. Considering that sex and the sex-chromosome abnormalities can have large effects on many phenotypes, including diseases, our method can make a significant contribution to DNA methylation studies that are based on microarray platforms.

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.

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.

Loss-of-function germline mutations in the PALB2 gene are associated with an increase of breast cancer risk. The purpose of this study was to characterise the spectrum of PALB2 mutations in women affected with breast or ovarian cancer from South-West Poland and West Ukraine. 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 PALB2 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. Targeted-sequencing identified two frameshift deletions: PALB2:c.509_510del; p.R170Ifs in three women affected with breast cancer and PALB2:c.172_175del;p.Q60Rfs in one woman affected with ovarian cancer. A number of other previously described missense (some predicted to be damaging by PolyPhen-2 and CADD) and synonymous mutations were also identified in this population. This study is consistent with previous reports that PALB2:c.509_510del and PALB2:c.172_175del are recurrent mutations associated with breast cancer predisposition in Polish women with a family history of the disease. Our study contributes to the accumulating evidence indicating that PALB2 should be included in genetic testing for breast cancer susceptibility in these populations to enhance risk assessment and management of women at high-risk of developing breast cancer. This data could also contribute to ongoing work that is assessing the possible association between ovarian cancer risk and PALB2 mutations for which there is currently no evidence.

BACKGROUND: Genetic variant effect prediction algorithms are used extensively in clinical genomics and research to determine the likely consequences of amino acid substitutions on protein function. It is vital that we better understand their accuracies and limitations because published performance metrics are confounded by serious problems of circularity and error propagation. Here, we derive three independent, functionally determined human mutation datasets, UniFun, BRCA1-DMS and TP53-TA, and employ them, alongside previously described datasets, to assess the pre-eminent variant effect prediction tools. RESULTS: Apparent accuracies of variant effect prediction tools were influenced significantly by the benchmarking dataset. Benchmarking with the assay-determined datasets UniFun and BRCA1-DMS yielded areas under the receiver operating characteristic curves in the modest ranges of 0.52 to 0.63 and 0.54 to 0.75, respectively, considerably lower than observed for other, potentially more conflicted datasets. CONCLUSIONS: These results raise concerns about how such algorithms should be employed, particularly in a clinical setting. Contemporary variant effect prediction tools are unlikely to be as accurate at the general prediction of functional impacts on proteins as reported prior. Use of functional assay-based datasets that avoid prior dependencies promises to be valuable for the ongoing development and accurate benchmarking of such tools.

Constitutional biological processes involve the generation of DNA double-strand breaks (DSBs). The production of such breaks and their subsequent resolution are also highly relevant to neurodegenerative diseases and cancer, in which extensive DNA fragmentation has been described Stephens et al. (2011), Blondet et al. (2001). Tchurikov et al. Tchurikov et al. (2011, 2013) have reported previously that frequent sites of DSBs occur in chromosomal domains involved in the co-ordinated expression of genes. This group report that hot spots of DSBs in human HEK293T cells often coincide with H3K4me3 marks, associated with active transcription Kravatsky et al. (2015) and that frequent sites of DNA double-strand breakage are likely to be relevant to cancer genomics Tchurikov et al. (2013, 2016) . Recently, they applied a RAFT (rapid amplification of forum termini) protocol that selects for blunt-ended DSB sites and mapped these to the human genome within defined co-ordinate 'windows'. In this paper, we re-analyse public RAFT data to derive sites of DSBs at the single-nucleotide level across the built genome for human HEK293T cells (https://figshare.com/s/35220b2b79eaaaf64ed8). This refined mapping, combined with accessory ENCODE data tracks and ribosomal DNA-related sequence annotations, will likely be of value for the design of clinically relevant targeted assays such as those for cancer susceptibility, diagnosis, treatment-matching and prognostication.

Neural injury triggers swift responses from glia, including glial migration and phagocytic clearance of damaged neurons. The transcriptional programs governing these complex innate glial immune responses are still unclear. Here, we describe a novel injury assay in adult Drosophila that elicits widespread glial responses in the ventral nerve cord (VNC). We profiled injury-induced changes in VNC gene expression by RNA sequencing (RNA-seq) and found that responsive genes fall into diverse signaling classes. One factor, matrix metalloproteinase-1 (MMP-1), is induced in Drosophila ensheathing glia responding to severed axons. Interestingly, glial induction of MMP-1 requires the highly conserved engulfment receptor Draper, as well as AP-1 and STAT92E. In MMP-1 depleted flies, glia do not properly infiltrate neuropil regions after axotomy and, as a consequence, fail to clear degenerating axonal debris. This work identifies Draper-dependent activation of MMP-1 as a novel cascade required for proper glial clearance of severed axons.

An apparently balanced t(2;3)(q37.3;q13.2) translocation that appears to segregate with renal cell carcinoma (RCC) has indicated potential areas to search for the elusive genetic basis of clear cell RCC. We applied Hi-Plex targeted sequencing to analyse germline DNA from 479 individuals affected with clear cell RCC for this breakpoint translocation and genetic variants in neighbouring genes on chromosome 2, ACKR3 and COPS8. While only synonymous variants were found in COPS8, one of the missense variants in ACKR3:c.892C>T, observed in 4/479 individuals screened (0.8%), was predicted likely to damage ACKR3 function. Identification of causal genes for RCC has potential clinical utility, where risk assessment and risk management can offer better outcomes, with surveillance for at-risk relatives and nephron sparing surgery through earlier intervention.

BACKGROUND: Global DNA methylation has been reported to be associated with urothelial cell carcinoma (UCC) by studies using blood samples collected at diagnosis. Using the Illumina HumanMethylation450 assay, we derived genome-wide measures of blood DNA methylation and assessed them for their prospective association with UCC risk. METHODS: We used 439 case-control pairs from the Melbourne Collaborative Cohort Study matched on age, sex, country of birth, DNA sample type, and collection period. Conditional logistic regression was used to compute odds ratios (OR) of UCC risk per s.d. of each genome-wide measure of DNA methylation and 95% confidence intervals (CIs), adjusted for potential confounders. We also investigated associations by disease subtype, sex, smoking, and time since blood collection. RESULTS: The risk of superficial UCC was decreased for individuals with higher levels of our genome-wide DNA methylation measure (OR=0.71, 95% CI: 0.54-0.94; P=0.02). This association was particularly strong for current smokers at sample collection (OR=0.47, 95% CI: 0.27-0.83). Intermediate levels of our genome-wide measure were associated with decreased risk of invasive UCC. Some variation was observed between UCC subtypes and the location and regulatory function of the CpGs included in the genome-wide measures of methylation. CONCLUSIONS: Higher levels of our genome-wide DNA methylation measure were associated with decreased risk of superficial UCC and intermediate levels were associated with reduced risk of invasive disease. These findings require replication by other prospective studies.

BACKGROUND: Previously, we described ROVER, a DNA variant caller which identifies genetic variants from PCR-targeted massively parallel sequencing (MPS) datasets generated by the Hi-Plex protocol. ROVER permits stringent filtering of sequencing chemistry-induced errors by requiring reported variants to appear in both reads of overlapping pairs above certain thresholds of occurrence. ROVER was developed in tandem with Hi-Plex and has been used successfully to screen for genetic mutations in the breast cancer predisposition gene PALB2. ROVER is applied to MPS data in BAM format and, therefore, relies on sequence reads being mapped to a reference genome. In this paper, we describe an improvement to ROVER, called UNDR ROVER (Unmapped primer-Directed ROVER), which accepts MPS data in FASTQ format, avoiding the need for a computationally expensive mapping stage. It does so by taking advantage of the location-specific nature of PCR-targeted MPS data. RESULTS: The UNDR ROVER algorithm achieves the same stringent variant calling as its predecessor with a significant runtime performance improvement. In one indicative sequencing experiment, UNDR ROVER (in its fastest mode) required 8-fold less sequential computation time than the ROVER pipeline and 13-fold less sequential computation time than a variant calling pipeline based on the popular GATK tool. UNDR ROVER is implemented in Python and runs on all popular POSIX-like operating systems (Linux, OS X). It requires as input a tab-delimited format file containing primer sequence information, a FASTA format file containing the reference genome sequence, and paired FASTQ files containing sequence reads. Primer sequences at the 5' end of reads associate read-pairs with their targeted amplicon and, thus, their expected corresponding coordinates in the reference genome. The primer-intervening sequence of each read is compared against the reference sequence from the same location and variants are identified using the same algorithm a

DNA breakage arises during a variety of biological processes, including transcription, replication and genome rearrangements. In the context of disease, extensive fragmentation of DNA has been described in cancer cells and during early stages of neurodegeneration (Stephens et al., 2011 Stephens et al. (2011) [5]; Blondet et al., 2001 Blondet et al. (2001) [1]). Stults et al. (2009) Stults et al. (2009) [6] reported that human rDNA gene clusters are hotspots for recombination and that rDNA restructuring is among the most common chromosomal alterations in adult solid tumours. As such, analysis of rDNA regions is likely to have significant prognostic and predictive value, clinically. Tchurikov et al. (2015a, 2016) Tchurikov et al. (2015a, 2016) [7], [9] have made major advances in this direction, reporting that sites of human genome double-strand breaks (DSBs) occur frequently at sites in rDNA that are tightly linked with active transcription - the authors used a RAFT (rapid amplification of forum termini) protocol that selects for blunt-ended sites. They reported the relative frequency of these rDNA DSBs within defined co-ordinate 'windows' of varying size and made these data (as well as the relevant 'raw' sequencing information) available to the public (Tchurikov et al., 2015b). Assay designs targeting rDNA DSB hotspots will benefit greatly from the publication of break sites at greater resolution. Here, we re-analyse public RAFT data and make available rDNA DSB co-ordinates to the single-nucleotide level.

Early-onset breast cancer (EOBC) causes substantial loss of life and productivity, creating a major burden among women worldwide. We analyzed 1,265,548 Hapmap3 single-nucleotide polymorphisms (SNP) among a discovery set of 3,523 EOBC incident cases and 2,702 population control women ages ≤ 51 years. The SNPs with smallest P values were examined in a replication set of 3,470 EOBC cases and 5,475 control women. We also tested EOBC association with 19,684 genes by annotating each gene with putative functional SNPs, and then combining their P values to obtain a gene-based P value. We examined the gene with smallest P value for replication in 1,145 breast cancer cases and 1,142 control women. The combined discovery and replication sets identified 72 new SNPs associated with EOBC (P < 4 × 10(-8)) located in six genomic regions previously reported to contain SNPs associated largely with later-onset breast cancer (LOBC). SNP rs2229882 and 10 other SNPs on chromosome 5q11.2 remained associated (P < 6 × 10(-4)) after adjustment for the strongest published SNPs in the region. Thirty-two of the 82 currently known LOBC SNPs were associated with EOBC (P < 0.05). Low power is likely responsible for the remaining 50 unassociated known LOBC SNPs. The gene-based analysis identified an association between breast cancer and the phosphofructokinase-muscle (PFKM) gene on chromosome 12q13.11 that met the genome-wide gene-based threshold of 2.5 × 10(-6). In conclusion, EOBC and LOBC seem to have similar genetic etiologies; the 5q11.2 region may contain multiple distinct breast cancer loci; and the PFKM gene region is worthy of further investigation. These findings should enhance our understanding of the etiology of breast cancer.

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 identificatio

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

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.

Linkage analysis, positional cloning, candidate gene mutation scanning and genome-wide association study approaches have all contributed significantly to our understanding of the underlying genetic architecture of breast cancer. Taken together, these approaches have identified genetic variation that explains approximately 30% of the overall familial risk of breast cancer, implying that more, and likely rarer, genetic susceptibility alleles remain to be discovered.

BACKGROUND: Population-based studies of breast cancer have estimated that at least some PALB2 mutations are associated with high breast cancer risk. For women carrying PALB2 mutations, knowing their carrier status could be useful in directing them towards effective cancer risk management and therapeutic strategies. We sought to determine whether morphological features of breast tumours can predict PALB2 germline mutation status. METHODS: Systematic pathology review was conducted on breast tumours from 28 female carriers of PALB2 mutations (non-carriers of other known high-risk mutations, recruited through various resources with varying ascertainment) and on breast tumours from a population-based sample of 828 Australian women diagnosed before the age of 60 years (which included 40 BRCA1 and 18 BRCA2 mutation carriers). Tumour morphological features of the 28 PALB2 mutation carriers were compared with those of 770 women without high-risk mutations. RESULTS: Tumours arising in PALB2 mutation carriers were associated with minimal sclerosis (odds ratio (OR)=19.7; 95% confidence interval (CI)=6.0-64.6; P=5 × 10(-7)). Minimal sclerosis was also a feature that distinguished PALB2 mutation carriers from BRCA1 (P=0.05) and BRCA2 (P=0.04) mutation carriers. CONCLUSION: This study identified minimal sclerosis to be a predictor of germline PALB2 mutation status. Morphological review can therefore facilitate the identification of women most likely to carry mutations in PALB2.

BACKGROUND: Dried blood (Guthrie card) spots provide an efficient way to collect and store blood specimens. DNA from this source has been utilised for a number of molecular analyses including genome-wide association studies, but only few studies have tested the feasibility of using it for epigenetic applications, particularly at a genome-wide level. RESULTS: In this study, we demonstrate the successful use of DNA isolated from archived dried blood spots for the Infinium HumanMethylation450 Beadchip, along with DNA from matched frozen buffy coats. We obtained high quality and reproducible genome-wide DNA methylation profiles using both sample types. We also report high correlations (r > 0.9907) between DNA obtained from matched dried blood spots and frozen buffy coats, sufficient to distinguish between unrelated individuals. CONCLUSIONS: We, thus, demonstrate that DNA from archived dried blood spots is suitable for genome-wide DNA methylation profiling.

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.

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.

The 6q25.1 locus was first identified via a genome-wide association study (GWAS) in Chinese women and marked by single nucleotide polymorphism (SNP) rs2046210, approximately 180 Kb upstream of ESR1. There have been conflicting reports about the association of this locus with breast cancer in Europeans, and a GWAS in Europeans identified a different SNP, tagged here by rs12662670. We examined the associations of both SNPs in up to 61,689 cases and 58,822 controls from forty-four studies collaborating in the Breast Cancer Association Consortium, of which four studies were of Asian and 39 of European descent. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI). Case-only analyses were used to compare SNP effects in Estrogen Receptor positive (ER+) versus negative (ER-) tumours. Models including both SNPs were fitted to investigate whether the SNP effects were independent. Both SNPs are significantly associated with breast cancer risk in both ethnic groups. Per-allele ORs are higher in Asian than in European studies [rs2046210: OR (A/G) = 1.36 (95% CI 1.26-1.48), p = 7.6 × 10(-14) in Asians and 1.09 (95% CI 1.07-1.11), p = 6.8 × 10(-18) in Europeans. rs12662670: OR (G/T) = 1.29 (95% CI 1.19-1.41), p = 1.2 × 10(-9) in Asians and 1.12 (95% CI 1.08-1.17), p = 3.8 × 10(-9) in Europeans]. SNP rs2046210 is associated with a significantly greater risk of ER- than ER+ tumours in Europeans [OR (ER-) = 1.20 (95% CI 1.15-1.25), p = 1.8 × 10(-17) versus OR (ER+) = 1.07 (95% CI 1.04-1.1), p = 1.3 × 10(-7), p(heterogeneity) = 5.1 × 10(-6)]. In these Asian studies, by contrast, there is no clear evidence of a differential association by tumour receptor status. Each SNP is associated with risk after adjustment for the other SNP. These results suggest the presence of two variants at 6q25.1 each independently associated with breast cancer risk in Asians and in Europeans. Of these two, the one tagged by rs2046210 is associated with a greater risk of

The 6q25.1 locus was first identified via a genome-wide association study (GWAS) in Chinese women and marked by single nucleotide polymorphism (SNP) rs2046210, approximately 180 Kb upstream of ESR1. There have been conflicting reports about the association of this locus with breast cancer in Europeans, and a GWAS in Europeans identified a different SNP, tagged here by rs12662670. We examined the associations of both SNPs in up to 61,689 cases and 58,822 controls from forty-four studies collaborating in the Breast Cancer Association Consortium, of which four studies were of Asian and 39 of European descent. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI). Case-only analyses were used to compare SNP effects in Estrogen Receptor positive (ER+) versus negative (ER-) tumours. Models including both SNPs were fitted to investigate whether the SNP effects were independent. Both SNPs are significantly associated with breast cancer risk in both ethnic groups. Per-allele ORs are higher in Asian than in European studies [rs2046210: OR (A/G) = 1.36 (95% CI 1.26-1.48), p = 7.6x10(-14) in Asians and 1.09 (95% CI 1.07-1.11), p = 6.8x10(-18) in Europeans. rs12662670: OR (G/T) = 1.29 (95% CI 1.19-1.41), p = 1.2x10(-9) in Asians and 1.12 (95% CI 1.08-1.17), p = 3.8x10(-9) in Europeans]. SNP rs2046210 is associated with a significantly greater risk of ER- than ER+ tumours in Europeans [OR (ER-) = 1.20 (95% CI 1.15-1.25), p = 1.8x10(-17) versus OR (ER+) = 1.07 (95% CI 1.04-1.1), p = 1.3x10(-7), p(heterogeneity) = 5.1x10(-6)]. In these Asian studies, by contrast, there is no clear evidence of a differential association by tumour receptor status. Each SNP is associated with risk after adjustment for the other SNP. These results suggest the presence of two variants at 6q25.1 each independently associated with breast cancer risk in Asians and in Europeans. Of these two, the one tagged by rs2046210 is associated with a greater risk of ER- tumours.

NTRODUCTION: As a group, women who carry germline mutations in partner and localizer of breast cancer 2 susceptibility protein (PALB2) are at increased risk of breast cancer. Little is known about by how much or whether risk differs by mutation or family history, owing to the paucity of studies of cases unselected for family history. METHODS: We screened 1,403 case probands for PALB2 mutations in a population-based study of Australian women with invasive breast cancer stratified by age at onset. The age-specific risk of breast cancer was estimated from the cancer histories of first- and second-degree relatives of mutation-carrying probands using a modified segregation analysis that included a polygenic modifier and was conditioned on the carrier case proband. Further screening for PALB2 c.3113G > A (W1038X) was conducted for 779 families with multiple cases of breast cancer ascertained through family cancer clinics in Australia and New Zealand and 764 population-based controls. RESULTS: We found five independent case probands in the population-based sample with the protein-truncating mutation PALB2 c.3113G > A (W1038X); 2 of 695 were diagnosed before age 40 years and 3 of 708 were diagnosed when between ages 40 and 59 years. Both of the two early-onset carrier case probands had very strong family histories of breast cancer. Further testing found that the mutation segregated with breast cancer in these families. No c.3113G > A (W1038X) carriers were found in 764 population-based unaffected controls. The hazard ratio was estimated to be 30.1 (95% confidence interval (CI), 7.5 to 120; P < 0.0001), and the corresponding cumulative risk estimates were 49% (95% CI, 15 to 93) to age 50 and 91% (95% CI, 44 to 100) to age 70. We found another eight families carrying this mutation in 779 families with multiple cases of breast cancer ascertained through family cancer clinics. CONCLUSIONS: The PALB2 c.3113G > A mutation appears to be associated with substantial risks of breast