Rare copy number variants observed in hereditary breast cancer cases disrupt genes in estrogen signaling and TP53 tumor suppression network

Breast cancer is the most common cancer in women in developed countries, and the contribution of genetic susceptibility to breast cancer development has been well-recognized. However, a great proportion of these hereditary predisposing factors still remain unidentified. To examine the contribution of rare copy number variants (CNVs) in breast cancer predisposition, high-resolution genome-wide scans were performed on genomic DNA of 103 BRCA1, BRCA2, and PALB2 mutation negative familial breast cancer cases and 128 geographically matched healthy female controls; for replication an independent cohort of 75 similarly mutation negative young breast cancer patients was used. All observed rare variants were confirmed by independent methods. The studied breast cancer cases showed a consistent increase in the frequency of rare CNVs when compared to controls. Furthermore, the biological networks of the disrupted genes differed between the two groups. In familial cases the observed mutations disrupted genes, which were significantly overrepresented in cellular functions related to maintenance of genomic integrity, including DNA double-strand break repair (P = 0.0211). Biological network analysis in the two independent breast cancer cohorts showed that the disrupted genes were closely related to estrogen signaling and TP53 centered tumor suppressor network. These results suggest that rare CNVs represent an alternative source of genetic variation influencing hereditary risk for breast cancer.     

The nature of intraspecific and interspecific genome size variation in taxonomically complex eyebrights

Background and aimsGenome size varies considerably across the diversity of plant life. Although genome size is, by definition, affected by genetic presence/absence variants, which are ubiquitous in population sequencing studies, genome size is often treated as an intrinsic property of a species. Here, we studied intra- and interspecific genome size variation in taxonomically complex British eyebrights (Euphrasia, Orobanchaceae). Our aim is to document genome size diversity and investigate underlying evolutionary processes shaping variation between individuals, populations and species.MethodsWe generated genome size data for 192 individuals of diploid and tetraploid Euphrasia and analysed genome size variation in relation to ploidy, taxonomy, population affiliation and geography. We further compared the genomic repeat content of 30 samples.Key resultsWe found considerable intraspecific genome size variation, and observed isolation-by-distance for genome size in outcrossing diploids. Tetraploid Euphrasia showed contrasting patterns, with genome size increasing with latitude in outcrossing Euphrasia arctica, but with little genome size variation in the highly selfing Euphrasia micrantha. Interspecific differences in genome size and the genomic proportions of repeat sequences were small.ConclusionsWe show the utility of treating genome size as the outcome of polygenic variation. Like other types of genetic variation, such as single nucleotide polymorphisms, genome size variation may be affected by ongoing hybridization and the extent of population subdivision. In addition to selection on associated traits, genome size is predicted to be affected indirectly by selection due to pleiotropy of the underlying presence/absence variants.     

Targeted Sequencing of the Mitochondrial Genome of Women at High Risk of Breast Cancer without Detectable Mutations in BRCA1/2

Breast Cancer is a complex multifactorial disease for which high-penetrance mutations have been identified. Approaches used to date have identified genomic features explaining about 50% of breast cancer heritability. A number of low- to medium penetrance alleles (per-allele odds ratio < 1.5 and 4.0, respectively) have been identified, suggesting that the remaining heritability is likely to be explained by the cumulative effect of such alleles and/or by rare high-penetrance alleles. Relatively few studies have specifically explored the mitochondrial genome for variants potentially implicated in breast cancer risk. For these reasons, we propose an exploration of the variability of the mitochondrial genome in individuals diagnosed with breast cancer, having a positive breast cancer family history but testing negative for BRCA1/2 pathogenic mutations. We sequenced the mitochondrial genome of 436 index breast cancer cases from the GENESIS study. As expected, no pathogenic genomic pattern common to the 436 women included in our study was observed. The mitochondrial genes MT-ATP6 and MT-CYB were observed to carry the highest number of variants in the study. The proteins encoded by these genes are involved in the structure of the mitochondrial respiration chain, and variants in these genes may impact reactive oxygen species production contributing to carcinogenesis. More functional and epidemiological studies are needed to further investigate to what extent variants identified may influence familial breast cancer risk.     

The Validation and Clinical Implementation of BRCAplus: A Comprehensive High-Risk Breast Cancer Diagnostic Assay

Breast cancer is the most commonly diagnosed cancer in women, with 10% of disease attributed to hereditary factors. Although BRCA1 and BRCA2 account for a high percentage of hereditary cases, there are more than 25 susceptibility genes that differentially impact the risk for breast cancer. Traditionally, germline testing for breast cancer was performed by Sanger dideoxy terminator sequencing in a reflexive manner, beginning with BRCA1 and BRCA2. The introduction of next-generation sequencing (NGS) has enabled the simultaneous testing of all genes implicated in breast cancer resulting in diagnostic labs offering large, comprehensive gene panels. However, some physicians prefer to only test for those genes in which established surveillance and treatment protocol exists. The NGS based BRCAplus test utilizes a custom tiled PCR based target enrichment design and bioinformatics pipeline coupled with array comparative genomic hybridization (aCGH) to identify mutations in the six high-risk genes: BRCA1, BRCA2, PTEN, TP53, CDH1, and STK11. Validation of the assay with 250 previously characterized samples resulted in 100% detection of 3,025 known variants and analytical specificity of 99.99%. Analysis of the clinical performance of the first 3,000 BRCAplus samples referred for testing revealed an average coverage greater than 9,000X per target base pair resulting in excellent specificity and the sensitivity to detect low level mosaicism and allele-drop out. The unique design of the assay enabled the detection of pathogenic mutations missed by previous testing. With the abundance of NGS diagnostic tests being released, it is essential that clinicians understand the advantages and limitations of different test designs.     

Exome Sequencing of Germline DNA from Non-BRCA1/2 Familial Breast Cancer Cases Selected on the Basis of aCGH Tumor Profiling

The bulk of familial breast cancer risk (∼70%) cannot be explained by mutations in the known predisposition genes, primarily BRCA1 and BRCA2. Underlying genetic heterogeneity in these cases is the probable explanation for the failure of all attempts to identify further high-risk alleles. While exome sequencing of non-BRCA1/2 breast cancer cases is a promising strategy to detect new high-risk genes, rational approaches to the rigorous pre-selection of cases are needed to reduce heterogeneity. We selected six families in which the tumours of multiple cases showed a specific genomic profile on array comparative genomic hybridization (aCGH). Linkage analysis in these families revealed a region on chromosome 4 with a LOD score of 2.49 under homogeneity. We then analysed the germline DNA of two patients from each family using exome sequencing. Initially focusing on the linkage region, no potentially pathogenic variants could be identified in more than one family. Variants outside the linkage region were then analysed, and we detected multiple possibly pathogenic variants in genes that encode DNA integrity maintenance proteins. However, further analysis led to the rejection of all variants due to poor co-segregation or a relatively high allele frequency in a control population. We concluded that using CGH results to focus on a sub-set of families for sequencing analysis did not enable us to identify a common genetic change responsible for the aggregation of breast cancer in these families. Our data also support the emerging view that non-BRCA1/2 hereditary breast cancer families have a very heterogeneous genetic basis.     

Whole Exome Sequencing Suggests Much of Non-BRCA1/BRCA2 Familial Breast Cancer Is Due to Moderate and Low Penetrance Susceptibility Alleles

The identification of the two most prevalent susceptibility genes in breast cancer, BRCA1 and BRCA2, was the beginning of a sustained effort to uncover new genes explaining the missing heritability in this disease. Today, additional high, moderate and low penetrance genes have been identified in breast cancer, such as P53, PTEN, STK11, PALB2 or ATM, globally accounting for around 35 percent of the familial cases. In the present study we used massively parallel sequencing to analyze 7 BRCA1/BRCA2 negative families, each having at least 6 affected women with breast cancer (between 6 and 10) diagnosed under the age of 60 across generations. After extensive filtering, Sanger sequencing validation and co-segregation studies, variants were prioritized through either control-population studies, including up to 750 healthy individuals, or case-control assays comprising approximately 5300 samples. As a result, a known moderate susceptibility indel variant (CHEK2 1100delC) and a catalogue of 11 rare variants presenting signs of association with breast cancer were identified. All the affected genes are involved in important cellular mechanisms like DNA repair, cell proliferation and survival or cell cycle regulation. This study highlights the need to investigate the role of rare variants in familial cancer development by means of novel high throughput analysis strategies optimized for genetically heterogeneous scenarios. Even considering the intrinsic limitations of exome resequencing studies, our findings support the hypothesis that the majority of non-BRCA1/BRCA2 breast cancer families might be explained by the action of moderate and/or low penetrance susceptibility alleles.     

Direct Visualization of the Highly Polymorphic RNU2 Locus in Proximity to the BRCA1 Gene

Although the breast cancer susceptibility gene BRCA1 is one of the most extensively characterized genetic loci, much less is known about its upstream variable number tandem repeat element, the RNU2 locus. RNU2 encodes the U2 small nuclear RNA, an essential splicing element, but this locus is missing from the human genome assembly due to the inherent difficulty in the assembly of repetitive sequences. To fill the gap between RNU2 and BRCA1, we have reconstructed the physical map of this region by re-examining genomic clone sequences of public databases, which allowed us to precisely localize the RNU2 array 124 kb telomeric to BRCA1. We measured by performing FISH analyses on combed DNA for the first time the exact number of repeats carried by each of the two alleles in 41 individuals and found a range of 6-82 copies and a level of heterozygosity of 98%. The precise localisation of the RNU2 locus in the genome reference assembly and the implementation of a new technical tool to study it will make the detailed exploration of this locus possible. This recently neglected macrosatellite could be valuable for evaluating the potential role of structural variations in disease due to its location next to a major cancer susceptibility gene.     

Secondary variants in individuals undergoing exome sequencing: screening of 572 individuals identifies high-penetrance mutations in cancer-susceptibility genes

Genome- and exome-sequencing costs are continuing to fall, and many individuals are undergoing these assessments as research participants and patients. The issue of secondary (so-called incidental) findings in exome analysis is controversial, and data are needed on methods of detection and their frequency. We piloted secondary variant detection by analyzing exomes for mutations in cancer-susceptibility syndromes in subjects ascertained for atherosclerosis phenotypes. We performed exome sequencing on 572 ClinSeq participants, and in 37 genes, we interpreted variants that cause high-penetrance cancer syndromes by using an algorithm that filtered results on the basis of mutation type, quality, and frequency and that filtered mutation-database entries on the basis of defined categories of causation. We identified 454 sequence variants that differed from the human reference. Exclusions were made on the basis of sequence quality (26 variants) and high frequency in the cohort (77 variants) or dbSNP (17 variants), leaving 334 variants of potential clinical importance. These were further filtered on the basis of curation of literature reports. Seven participants, four of whom were of Ashkenazi Jewish descent and three of whom did not meet family-history-based referral criteria, had deleterious BRCA1 or BRCA2 mutations. One participant had a deleterious SDHC mutation, which causes paragangliomas. Exome sequencing, coupled with multidisciplinary interpretation, detected clinically important mutations in cancer-susceptibility genes; four of such mutations were in individuals without a significant family history of disease. We conclude that secondary variants of high clinical importance will be detected at an appreciable frequency in exomes, and we suggest that priority be given to the development of more efficient modes of interpretation with trials in larger patient groups.     

Prevalence of 185delAG and 5382insC mutations in BRCA1, and 6174delT in BRCA2 in women of Ashkenazi Jewish origin in southern Brazil

Certain mutations in BRCA1 and BRCA2 genes are frequent in the Ashkenazi Jewish population. Several factors contribute to this increased frequency, including consanguineous marriages and an event known as a “bottleneck”, which occurred in the past and caused a drastic reduction in the genetic variability of this population. Several studies were performed over the years in an attempt to elucidate the role of BRCA1 and BRCA2 genes in susceptibility to breast cancer. The aim of this study was to estimate the carrier frequency of certain common mutations in the BRCA1 (185delAG and 5382insC) and BRCA2 (6174delT) genes in an Ashkenazi Jewish population from Porto Alegre, Brazil. Molecular analyses were done by PCR followed by RFLP (ACRS). The carrier frequencies for BRCA1 185delAG and 5382insC were 0.78 and 0 respectively, and 0.4 for the BRCA2 6174deT mutation. These findings are similar to those of some prior studies but differ from others, possibly due to excluding individuals with a personal or family history of cancer. Our sample was drawn from the community group and included individuals with or without a family or personal history of cancer. Furthermore, increased dispersion among Ashkenazi subpopulations may be the result of strong genetic drift and/or admixture. It is therefore necessary to consider the effects of local admixture on the mismatch distributions of various Jewish populations.     

RAD51 and Breast Cancer Susceptibility: No Evidence for Rare Variant Association in the Breast Cancer Family Registry Study

Background

Although inherited breast cancer has been associated with germline mutations in genes that are functionally involved in the DNA homologous recombination repair (HRR) pathway, including BRCA1, BRCA2, TP53, ATM, BRIP1, CHEK2 and PALB2, about 70% of breast cancer heritability remains unexplained. Because of their critical functions in maintaining genome integrity and already well-established associations with breast cancer susceptibility, it is likely that additional genes involved in the HRR pathway harbor sequence variants associated with increased risk of breast cancer. RAD51 plays a central biological function in DNA repair and despite the fact that rare, likely dysfunctional variants in three of its five paralogs, RAD51C, RAD51D, and XRCC2, have been associated with breast and/or ovarian cancer risk, no population-based case-control mutation screening data are available for the RAD51 gene. We thus postulated that RAD51 could harbor rare germline mutations that confer increased risk of breast cancer.

Methodology/Principal Findings

We screened the coding exons and proximal splice junction regions of the gene for germline sequence variation in 1,330 early-onset breast cancer cases and 1,123 controls from the Breast Cancer Family Registry, using the same population-based sampling and analytical strategy that we developed for assessment of rare sequence variants in ATM and CHEK2. In total, 12 distinct very rare or private variants were characterized in RAD51, with 10 cases (0.75%) and 9 controls (0.80%) carrying such a variant. Variants were either likely neutral missense substitutions (3), silent substitutions (4) or non-coding substitutions (5) that were predicted to have little effect on efficiency of the splicing machinery.

Conclusion

Altogether, our data suggest that RAD51 tolerates so little dysfunctional sequence variation that rare variants in the gene contribute little, if anything, to breast cancer susceptibility.     

BRCA1 and BRCA2 missense variants of high and low clinical significance influence lymphoblastoid cell line post-irradiation gene expression

The functional consequences of missense variants in disease genes are difficult to predict. We assessed if gene expression profiles could distinguish between BRCA1 or BRCA2 pathogenic truncating and missense mutation carriers and familial breast cancer cases whose disease was not attributable to BRCA1 or BRCA2 mutations (BRCAX cases). 72 cell lines from affected women in high-risk breast ovarian families were assayed after exposure to ionising irradiation, including 23 BRCA1 carriers, 22 BRCA2 carriers, and 27 BRCAX individuals. A subset of 10 BRCAX individuals carried rare BRCA1/2 sequence variants considered to be of low clinical significance (LCS). BRCA1 and BRCA2 mutation carriers had similar expression profiles, with some subclustering of missense mutation carriers. The majority of BRCAX individuals formed a distinct cluster, but BRCAX individuals with LCS variants had expression profiles similar to BRCA1/2 mutation carriers. Gaussian Process Classifier predicted BRCA1, BRCA2 and BRCAX status, with a maximum of 62% accuracy, and prediction accuracy decreased with inclusion of BRCAX samples carrying an LCS variant, and inclusion of pathogenic missense carriers. Similarly, prediction of mutation status with gene lists derived using Support Vector Machines was good for BRCAX samples without an LCS variant (82–94%), poor for BRCAX with an LCS (40–50%), and improved for pathogenic BRCA1/2 mutation carriers when the gene list used for prediction was appropriate to mutation effect being tested (71–100%). This study indicates that mutation effect, and presence of rare variants possibly associated with a low risk of cancer, must be considered in the development of array-based assays of variant pathogenicity.     

Whole genome sequencing of Turkish genomes reveals functional private alleles and impact of genetic interactions with Europe,Asia and Africa

Background

Turkey is a crossroads of major population movements throughout history and has been a hotspot of cultural interactions. Several studies have investigated the complex population history of Turkey through a limited set of genetic markers. However, to date, there have been no studies to assess the genetic variation at the whole genome level using whole genome sequencing. Here, we present whole genome sequences of 16 Turkish individuals resequenced at high coverage (32 × -48×).

Results

We show that the genetic variation of the contemporary Turkish population clusters with South European populations, as expected, but also shows signatures of relatively recent contribution from ancestral East Asian populations. In addition, we document a significant enrichment of non-synonymous private alleles, consistent with recent observations in European populations. A number of variants associated with skin color and total cholesterol levels show frequency differentiation between the Turkish populations and European populations. Furthermore, we have analyzed the 17q21.31 inversion polymorphism region (MAPT locus) and found increased allele frequency of 31.25% for H1/H2 inversion polymorphism when compared to European populations that show about 25% of allele frequency.

Conclusion

This study provides the first map of common genetic variation from 16 western Asian individuals and thus helps fill an important geographical gap in analyzing natural human variation and human migration. Our data will help develop population-specific experimental designs for studies investigating disease associations and demographic history in Turkey.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-963) contains supplementary material, which is available to authorized users.     

Variants in CHEK2 other than 1100delC do not make a major contribution to breast cancer susceptibility

We recently reported that a sequence variant in the cell-cycle-checkpoint kinase CHEK2 (CHEK2 1100delC) is a low-penetrance breast cancer-susceptibility allele in noncarriers of BRCA1 or BRCA2 mutations. To investigate whether other CHEK2 variants confer susceptibility to breast cancer, we screened the full CHEK2 coding sequence in BRCA1/2-negative breast cancer cases from 89 pedigrees with three or more cases of breast cancer. We identified one novel germline variant, R117G, in two separate families. To evaluate the possible association of R117G and two germline variants reported elsewhere, R145W and I157T with breast cancer, we screened 737 BRCA1/2-negative familial breast cancer cases from 605 families, 459 BRCA1/2-positive cases from 335 families, and 723 controls from the United Kingdom, the Netherlands, and North America. All three variants were rare in all groups, and none occurred at significantly elevated frequency in familial breast cancer cases compared with controls. These results indicate that 1100delC may be the only CHEK2 allele that makes an appreciable contribution to breast cancer susceptibility.     

Identification of Novel Candidate Genes for Early-Onset Colorectal Cancer Susceptibility

Approximately 25–30% of colorectal cancer (CRC) cases are expected to result from a genetic predisposition, but in only 5–10% of these cases highly penetrant germline mutations are found. The remaining CRC heritability is still unexplained, and may be caused by a hitherto-undefined set of rare variants with a moderately penetrant risk. Here we aimed to identify novel risk factors for early-onset CRC using whole-exome sequencing, which was performed on a cohort of CRC individuals (n = 55) with a disease onset before 45 years of age. We searched for genes that were recurrently affected by rare variants (minor allele frequency ≤0.001) with potentially damaging effects and, subsequently, re-sequenced the candidate genes in a replication cohort of 174 early-onset or familial CRC individuals. Two functionally relevant genes with low frequency variants with potentially damaging effects, PTPN12 and LRP6, were found in at least three individuals. The protein tyrosine phosphatase PTP-PEST, encoded by PTPN12, is a regulator of cell motility and LRP6 is a component of the WNT-FZD-LRP5-LRP6 complex that triggers WNT signaling. All variants in LRP6 were identified in individuals with an extremely early-onset of the disease (≤30 years of age), and two of the three variants showed increased WNT signaling activity in vitro. In conclusion, we present PTPN12 and LRP6 as novel candidates contributing to the heterogeneous susceptibility to CRC.     

The functional impact of BRCA1 BRCT domain variants using multiplexed DNA double-strand break repair assays

Pathogenic variants in BRCA1 are associated with a greatly increased risk of hereditary breast and ovarian cancer (HBOC). With the increased availability and affordability of genetic testing, many individuals have been identified with BRCA1 variants of uncertain significance (VUSs), which are individually detected in the population too infrequently to ascertain a clinical risk. Functional assays can be used to experimentally assess the effects of these variants. In this study, we used multiplexed DNA repair assays of variants in the BRCA1 carboxyl terminus to functionally characterize 2,271 variants for homology-directed repair function (HDR) and 1,427 variants for cisplatin resistance (CR). We found a high level of consistent results (Pearson’s r = 0.74) in the two multiplexed functional assays with non-functional variants located within regions of the BRCA1 protein necessary for its tumor suppression activity. In addition, functional categorizations of variants tested in the multiplex HDR and CR assays correlated with known clinical significance and with other functional assays for BRCA1 (Pearson’s r = 0.53 to 0.71). The results of the multiplex HDR and CR assays are useful resources for characterizing large numbers of BRCA1 VUSs.     

Personalized Oncogenomics: Clinical Experience with Malignant Peritoneal Mesothelioma Using Whole Genome Sequencing

Peritoneal mesothelioma is a rare and sometimes lethal malignancy that presents a clinical challenge for both diagnosis and management. Recent studies have led to a better understanding of the molecular biology of peritoneal mesothelioma. Translation of the emerging data into better treatments and outcome is needed. From two patients with peritoneal mesothelioma, we derived whole genome sequences, RNA expression profiles, and targeted deep sequencing data. Molecular data were made available for translation into a clinical treatment plan. Treatment responses and outcomes were later examined in the context of molecular findings. Molecular studies presented here provide the first reported whole genome sequences of peritoneal mesothelioma. Mutations in known mesothelioma-related genes NF2, CDKN2A, LATS2, amongst others, were identified. Activation of MET-related signaling pathways was demonstrated in both cases. A hypermutated phenotype was observed in one case (434 vs. 18 single nucleotide variants) and was associated with a favourable outcome despite sarcomatoid histology and multifocal disease. This study represents the first report of whole genome analyses of peritoneal mesothelioma, a key step in the understanding and treatment of this disease.