ShallowHRD: detection of homologous recombination deficiency from shallow whole genome sequencing.

SUMMARY: We introduce shallowHRD, a software tool to evaluate tumor homologous recombination deficiency (HRD) based on whole genome sequencing (WGS) at low coverage (shallow WGS or sWGS; ∼1X coverage). The tool, based on mining copy number alterations profile, implements a fast and straightforward procedure that shows 87.5% sensitivity and 90.5% specificity for HRD detection. shallowHRD could be instrumental in predicting response to poly(ADP-ribose) polymerase inhibitors, to which HRD tumors are selectively sensitive. shallowHRD displays efficiency comparable to most state-of-art approaches, is cost-effective, generates low-storable outputs and is also suitable for fixed-formalin paraffin embedded tissues. AVAILABILITY AND IMPLEMENTATION: shallowHRD R script and documentation are available at https://github.com/aeeckhou/shallowHRD. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Lack of evidence for CDK12 as an ovarian cancer predisposing gene.

CDK12 variants were investigated as a genetic susceptibility to ovarian cancer in a series of 416 unrelated and consecutive patients with ovarian carcinoma and who carry neither germline BRCA1 nor BRCA2 pathogenic variant. The presence of CDK12 variants was searched in germline DNA by massive parallel sequencing on pooled DNAs. The lack of detection of deleterious variants and the observed proportion of missense variants in the series of ovarian carcinoma patients as compared with all human populations strongly suggests that CDK12 is not an ovarian cancer predisposing gene.

Medullary Breast Carcinoma, a Triple-Negative Breast Cancer Associated with BCLG Overexpression.

Medullary breast carcinoma (MBC) is a rare subtype of triple-negative breast cancer with specific genomic features within the spectrum of basal-like carcinoma (BLC). In this study of 19 MBCs and 36 non-MBC BLCs, we refined the transcriptomic and genomic knowledge about this entity. Unsupervised and supervised analysis of transcriptomic profiles confirmed that MBC clearly differs from non-MBC BLC, with 92 genes overexpressed and 154 genes underexpressed in MBC compared with non-MBC BLC. Immunity-related pathways are the most differentially represented pathways in MBC compared with non-MBC BLC. The proapoptotic gene BCLG (official name BCL2L14) is by far the most intensely overexpressed gene in MBC. A quantitative RT-PCR validation study conducted in 526 breast tumors corresponding to all molecular subtypes documented the specificity of BCLG overexpression in MBC, which was confirmed at the protein level by immunohistochemistry. We also found that most MBCs belong to the immunomodulatory triple-negative breast cancer subtype. Using pan-genomic analysis, it was found that MBC harbors more losses of heterozygosity than non-MBC BLC. These observations corroborate the notion that MBC remains a distinct entity that could benefit from specific treatment strategies (such as deescalation or targeted therapy) adapted to this rare tumor type.

Ovarian Cancers Harboring Inactivating Mutations in CDK12 Display a Distinct Genomic Instability Pattern Characterized by Large Tandem Duplications.

CDK12 is a recurrently mutated gene in serous ovarian carcinoma, whose downregulation is associated with impaired expression of DNA damage repair genes and subsequent hypersensitivity to DNA-damaging agents and PARP1/2 inhibitors. In this study, we investigated the genomic landscape associated with CDK12 inactivation in patients with serous ovarian carcinoma. We show that CDK12 loss was consistently associated with a particular genomic instability pattern characterized by hundreds of tandem duplications of up to 10 megabases (Mb) in size. Tandem duplications were characterized by a bimodal (∼0.3 and ∼3 Mb) size distribution and overlapping microhomology at the breakpoints. This genomic instability, denoted as the CDK12 TD-plus phenotype, is remarkably distinct from other alteration patterns described in breast and ovarian cancers. The CDK12 TD-plus phenotype was associated with a greater than 10% gain in genomic content and occurred at a 3% to 4% rate in The Cancer Genome Atlas-derived and in-house cohorts of patients with serous ovarian carcinoma. Moreover, CDK12-inactivating mutations together with the TD-plus phenotype were also observed in prostate cancers. Our finding provides new insight toward deciphering the function of CDK12 in genome maintenance and oncogenesis. Cancer Res; 76(7); 1882-91. ©2016 AACR.

Genomic hallmarks of homologous recombination deficiency in invasive breast carcinomas.

Therapeutic strategies targeting Homologous Recombination Deficiency (HRD) in breast cancer requires patient stratification. The LST (Large-scale State Transitions) genomic signature previously validated for triple-negative breast carcinomas (TNBC) was evaluated as biomarker of HRD in luminal (hormone receptor positive) and HER2-overexpressing (HER2+) tumors. The LST genomic signature related to the number of large-scale chromosomal breakpoints in SNP-array tumor profile was applied to identify HRD in in-house and TCGA sets of breast tumors, in which the status of BRCA1/2 and other genes was also investigated. In the in-house dataset, HRD was predicted in 5% (20/385) of sporadic tumors luminal or HER2+ by the LST genomic signature and the inactivation of BRCA1, BRCA2 or RAD51C confirmed this prediction in 75% (12/16) of the tested cases. In 14% (6/43) of tumors occurring in BRCA1/2 mutant carriers, the corresponding wild-type allele was retained emphasizing the importance of determining the tumor status. In the TCGA luminal and HER2+ subtypes HRD incidence was estimated at 5% (18/329, 95%CI: 5-8%) and 2% (1/59, 95%CI: 2-9%), respectively. In TNBC cisplatin-based neo-adjuvant clinical trials, HRD is shown to be a necessary condition for cisplatin sensitivity. This analysis demonstrates the high performance of the LST genomic signature for HRD detection in breast cancers, which suggests its potential as a biomarker for genetic testing and patient stratification for clinical trials evaluating platinum salts and PARP inhibitors.

First description of a sporadic breast cancer in a woman with BRCA1 germline mutation.

We describe the case of a woman carrying a germline pathogenic BRCA1 mutation diagnosed with a breast cancer overexpressing HER2. Clinical presentation of the tumor, HER2-positivity, genomic profile and loss of the mutated BRCA1 allele in tumor evidence that BRCA1 is not inactivated in this breast cancer. It represents the first biological demonstration for the existence of a sporadic HER2-positive breast cancer independent from BRCA loss of function in a woman carrier of a deleterious BRCA1 mutation. In a context where targeted therapies based on BRCA loss of function in the tumor are developed, such case could have direct implications.

Loss of heterozygosity at 13q13 and 14q32 predicts BRCA2 inactivation in luminal breast carcinomas.

BRCA2 is the major high-penetrance predisposition gene for luminal (estrogen receptor [ER] positive) breast cancers. However, many BRCA2 mutant carriers lack family history of breast/ovarian cancers and do not benefit from genetic testing. Specific genomic features associated with BRCA2 inactivation in tumors could help identify patients for whom a genetic test for BRCA2 may be proposed. A series of ER-positive invasive ductal carcinomas (IDCs) including 30 carriers of BRCA2 mutations and 215 control cases was studied by single-nucleotide polymorphism (SNP) arrays. Cases and controls were stratified by grade and HER2 status. Independently, 7 BRCA2 and 51 control cases were used for validation. Absolute copy number and Loss of heterozygosity (LOH) profiles were obtained from SNP arrays by the genome alteration print (GAP) method. BRCA2 tumors were observed to display a discriminatively greater number of chromosomal breaks calculated after filtering out and smoothing <3 Mb variations. This argues for a BRCA2-associated genomic instability responsible for long-segment aberrations. Co-occurrence of two genomic features-LOH of 13q13 and 14q32-was found to predict BRCA2 status with 90% of sensitivity and 87% of specificity in discovery series of high-grade HER2-negative IDCs and 100% of sensitivity and 88% of specificity in an independent series of 58 IDCs. Estimated positive predictive value was 17.2% (confidence interval: 6.7-33.5) in the whole series. In conclusion, the simplified BRCA2 classifier based on the co-occurrence of LOH at 13q13 and 14q32 could provide an indication to test for BRCA2 mutation in patients with ER-positive IDC.

Germline BAP1 mutations predispose to renal cell carcinomas.

The genetic cause of some familial nonsyndromic renal cell carcinomas (RCC) defined by at least two affected first-degree relatives is unknown. By combining whole-exome sequencing and tumor profiling in a family prone to cases of RCC, we identified a germline BAP1 mutation c.277A>G (p.Thr93Ala) as the probable genetic basis of RCC predisposition. This mutation segregated with all four RCC-affected relatives. Furthermore, BAP1 was found to be inactivated in RCC-affected individuals from this family. No BAP1 mutations were identified in 32 familial cases presenting with only RCC. We then screened for germline BAP1 deleterious mutations in familial aggregations of cancers within the spectrum of the recently described BAP1-associated tumor predisposition syndrome, including uveal melanoma, malignant pleural mesothelioma, and cutaneous melanoma. Among the 11 families that included individuals identified as carrying germline deleterious BAP1 mutations, 6 families presented with 9 RCC-affected individuals, demonstrating a significantly increased risk for RCC. This strongly argues that RCC belongs to the BAP1 syndrome and that BAP1 is a RCC-predisposition gene.

Ploidy and large-scale genomic instability consistently identify basal-like breast carcinomas with BRCA1/2 inactivation.

BRCA1 inactivation is a frequent event in basal-like breast carcinomas (BLC). However, BRCA1 can be inactivated by multiple mechanisms and determining its status is not a trivial issue. As an alternate approach, we profiled 65 BLC cases using single-nucleotide polymorphism arrays to define a signature of BRCA1-associated genomic instability. Large-scale state transitions (LST), defined as chromosomal break between adjacent regions of at least 10 Mb, were found to be a robust indicator of BRCA1 status in this setting. Two major ploidy-specific cutoffs in LST distributions were sufficient to distinguish highly rearranged BLCs with 85% of proven BRCA1-inactivated cases from less rearranged BLCs devoid of proven BRCA1-inactivated cases. The genomic signature we defined was validated in a second independent series of 55 primary BLC cases and 17 BLC-derived tumor cell lines. High numbers of LSTs resembling BRCA1-inactivated BLC were observed in 4 primary BLC cases and 2 BLC cell lines that harbored BRCA2 mutations. Overall, the genomic signature we defined predicted BRCA1/2 inactivation in BLCs with 100% sensitivity and 90% specificity (97% accuracy). This assay may ease the challenge of selecting patients for genetic testing or recruitment to clinical trials of novel emerging therapies that target DNA repair deficiencies in cancer.

A whole-genome massively parallel sequencing analysis of BRCA1 mutant oestrogen receptor-negative and -positive breast cancers.

BRCA1 encodes a tumour suppressor protein that plays pivotal roles in homologous recombination (HR) DNA repair, cell-cycle checkpoints, and transcriptional regulation. BRCA1 germline mutations confer a high risk of early-onset breast and ovarian cancer. In more than 80% of cases, tumours arising in BRCA1 germline mutation carriers are oestrogen receptor (ER)-negative; however, up to 15% are ER-positive. It has been suggested that BRCA1 ER-positive breast cancers constitute sporadic cancers arising in the context of a BRCA1 germline mutation rather than being causally related to BRCA1 loss-of-function. Whole-genome massively parallel sequencing of ER-positive and ER-negative BRCA1 breast cancers, and their respective germline DNAs, was used to characterize the genetic landscape of BRCA1 cancers at base-pair resolution. Only BRCA1 germline mutations, somatic loss of the wild-type allele, and TP53 somatic mutations were recurrently found in the index cases. BRCA1 breast cancers displayed a mutational signature consistent with that caused by lack of HR DNA repair in both ER-positive and ER-negative cases. Sequencing analysis of independent cohorts of hereditary BRCA1 and sporadic non-BRCA1 breast cancers for the presence of recurrent pathogenic mutations and/or homozygous deletions found in the index cases revealed that DAPK3, TMEM135, KIAA1797, PDE4D, and GATA4 are potential additional drivers of breast cancers. This study demonstrates that BRCA1 pathogenic germline mutations coupled with somatic loss of the wild-type allele are not sufficient for hereditary breast cancers to display an ER-negative phenotype, and has led to the identification of three potential novel breast cancer genes (ie DAPK3, TMEM135, and GATA4).

Leukemic phase of follicular lymphomas: an atypical presentation.

Follicular lymphomas (FLs) are frequent B-cell derived malignancies, generally demonstrating an indolent evolution. Although circulating FL cells may be detected by high-resolution analysis, bloodstream involvement by FL cells is unusual. We observed in 10 patients a leukemic phase of FL at the onset of the disease. Six of them had concomitant lymph node involvement and all of them required treatment at diagnosis due to a high tumor burden, whereas four patients had pure FL-cell leukemia, which was associated with a more indolent clinical outcome. The detection of a leukemic phase should therefore be studied as a potential prognosis marker in further studies.

Syntenic relationships between genomic profiles of fiber-induced murine and human malignant mesothelioma.

Malignant mesothelioma (MM) is an aggressive tumor with a poor prognosis mainly linked to past asbestos exposure. Murine models of MM based on fiber exposure have been developed to elucidate the mechanism of mesothelioma formation. Genomic alterations in murine MM have now been partially characterized. To gain insight into the pathophysiology of mesothelioma, 16 murine and 35 human mesotheliomas were characterized by array-comparative genomic hybridization and were screened for common genomic alterations. Alteration of the 9p21 human region, often by biallelic deletion, was the most frequent alteration in both species, in agreement with the CDKN2A/CDKN2B locus deletion in human disease and murine models. Other shared aberrations were losses of 1p36.3-p35 and 13q14-q33 and gains of 5p15.3-p13 regions. However, some differences were noted, such as absence of recurrent alterations in mouse regions corresponding to human chromosome 22. Comparison between altered recurrent regions in asbestos-exposed and non-asbestos-exposed patients showed a significant difference in the 14q11.2-q21 region, which was also lost in fiber-induced murine mesothelioma. A correlation was also demonstrated between genomic instability and tumorigenicity of human mesothelioma xenografts in nude mice. Overall, these data show similarities between murine and human disease, and contribute to the understanding of the influence of fibers in the pathogenesis of mesothelioma and validation of the murine model for preclinical testing.

Genome Alteration Print (GAP): a tool to visualize and mine complex cancer genomic profiles obtained by SNP arrays.

We describe a method for automatic detection of absolute segmental copy numbers and genotype status in complex cancer genome profiles measured with single-nucleotide polymorphism (SNP) arrays. The method is based on pattern recognition of segmented and smoothed copy number and allelic imbalance profiles. Assignments were verified by DNA indexes of primary tumors and karyotypes of cell lines. The method performs well even for poor-quality data, low tumor content, and highly rearranged tumor genomes.

Beta-catenin status in paediatric medulloblastomas: correlation of immunohistochemical expression with mutational status, genetic profiles, and clinical characteristics.

Medulloblastoma is the most frequent malignant paediatric brain tumour. The activation of the Wnt/beta-catenin pathway occurs in 10-15% of medulloblastomas and has been recently described as a marker for favourable patient outcome. We report a series of 72 paediatric medulloblastomas evaluated for beta-catenin protein expression, CTNNB1 mutations, and comparative genomic hybridization. Gene expression profiles were also available in a subset of 40 cases. Immunostaining of beta-catenin showed extensive nuclear staining (>50% of the tumour cells) in six cases and focal nuclear staining (<10% of cells) in three cases. The other cases either exhibited a signal strictly limited to the cytoplasm (58 cases) or were negative (five cases). CTNNB1 mutations were detected in all beta-catenin extensively nucleopositive cases. The expression profiles of these cases documented strong activation of the Wnt/beta-catenin pathway. Remarkably, five out of these six tumours showed a complete loss of chromosome 6. In contrast, cases with focal nuclear beta-catenin staining, as well as tumours with negative or cytoplasmic staining, never demonstrated CTNNB1 mutation, Wnt/beta-catenin pathway activation or chromosome 6 loss. Patients with extensive nuclear staining were significantly older at diagnosis and were in continuous complete remission after a mean follow-up of 75.7 months (range 27.5-121.2 months) from diagnosis. All three patients with focal nuclear staining of beta-catenin died within 36 months from diagnosis. Altogether, these data confirm and extend previous observations that CTNNB1-mutated tumours represent a distinct molecular subgroup of medulloblastomas with favourable outcome, indicating that therapy de-escalation should be considered. International consensus on the definition criteria of this distinct medulloblastoma subgroup should be achieved.

High frequency of TP53 mutation in BRCA1 and sporadic basal-like carcinomas but not in BRCA1 luminal breast tumors.

Breast tumors with a germ-line mutation of BRCA1 (BRCA1 tumors) and basal-like carcinoma (BLC) are associated with a high rate of TP53 mutation. Because BRCA1 tumors frequently display a basal-like phenotype, this study was designed to determine whether TP53 mutations are correlated with the hereditary BRCA1 mutated status or the particular phenotype of these tumors. The TP53 gene status was first investigated in a series of 35 BRCA1 BLCs using immunohistochemistry, direct sequencing of the coding sequence, and functional analysis of separated alleles in yeast, and compared with the TP53 status in a series of 38 sporadic (nonhereditary) BLCs. Using this sensitive approach, TP53 was found to be frequently mutated in both BRCA1 (34 of 35, 97%) and sporadic (35 of 38, 92%) BLCs. However, the spectrum of mutation was different, particularly with a higher rate of complex mutations, such as insertion/deletion, in BRCA1 BLCs than in the sporadic group [14 of 33 (42%) and 3 of 34 (9%), [corrected] respectively; P = 0.002]. Secondly, the incidence of TP53 mutations was analyzed in 19 BRCA1 luminal tumors using the same strategy. Interestingly, only 10 of these 19 tumors were mutated (53%), a frequency similar to that found in grade-matched sporadic luminal tumors. In conclusion, TP53 mutation is highly recurrent in BLCs independently of BRCA1 status, but not a common feature of BRCA1 luminal tumors.

Genetic heterogeneity versus molecular analysis of prion susceptibility in neuroblasma N2a sublines.

The neuroblastoma-derived cell line N2a is permissive to certain prion strains but resistant sublines unable to accumulate the pathological proteinase-K resistant form of the prion protein can be isolated. We compared for gene expression and phenotypes different N2a sublines that were susceptible or resistant to the 22L prion strain. Karyotypes and comparative genomic hybridization arrays revealed chromosomal imbalances but did not demonstrate a characteristic profile of genomic alterations linked to prion susceptibility. Likewise, we showed that this phenotype was not dependent on the binding of PrPres, the expression of the prion protein gene, or on its primary sequence. We completed this analysis by looking using real-time quantitative PCR at the expression of a set of genes encoding proteins linked to prion biology. None of the candidates could account by itself for the infection phenotype, nevertheless sublines had distinct transcriptional profiles. Taken together, our results do not support a role for specific genomic abnormalities and possible candidate proteins in N2a prion susceptibility. They also reveal genetic heterogeneity among the sublines and serve as a guidance for further investigation into the molecular mechanisms of prion infection.

BAC array CGH distinguishes mutually exclusive alterations that define clinicogenetic subtypes of gliomas.

The pathological classification of gliomas constitutes a critical step of the clinical management of patients, yet it is frequently challenging. To assess the relationship between genetic abnormalities and clinicopathological characteristics, we have performed a genetic and clinical analysis of a series of gliomas. A total of 112 gliomas were analyzed by comparative genomic hybridization on a BAC array with a 1 megabase resolution. Altered regions were identified and correlation analysis enabled to retrieve significant associations and exclusions. Whole chromosomes (chrs) 1p and 19q losses with centromeric breakpoints and EGFR high level amplification were found to be mutually exclusive, permitting identification of 3 distinct, nonoverlapping groups of tumors with striking clinicopathological differences. Type A tumors with chrs 1p and 19q co-deletion exhibited an oligodendroglial phenotype and a longer patient survival. Type B tumors were characterized by EGFR amplification. They harbored a WHO high grade of malignancy and a short patient survival. Finally, type C tumors displayed none of the previous patterns but the presence of chr 7 gain, chr 9p deletion and/or chr 10 loss. It included astrocytic tumors in patients younger than in type B and whose prognosis was highly dependent upon the number of alterations. A multivariate analysis based on a Cox model shows that age, WHO grade and genomic type provide complementary prognostic informations. Finally, our results highlight the potential of a whole-genome analysis as an additional diagnostic in cases of unclear conventional genetic findings.

X inactive-specific transcript RNA coating and genetic instability of the X chromosome in BRCA1 breast tumors.

Identification among breast tumors of those arising in a hereditary BRCA1 context remains a medical challenge. Abnormalities in X chromosome copy number and in the epigenetic stability of the inactive X chromosome (Xi) have been proposed to characterize BRCA1 breast tumors. In particular, it has been proposed that loss of BRCA1 function can lead to loss of X inactive-specific transcript (XIST) RNA association with the Xi. However, few studies have addressed this issue in a sufficiently large series of BRCA1 primary tumors. Here we assess X-chromosome status using single-cell (RNA and DNA fluorescence in situ hybridization) and global genomic (array-comparative genomic hybridization and allelotyping) approaches on a series of 11 well-defined BRCA1 tumors. We show that many or most cells of the tumors contain one or more XIST RNA domains. Furthermore, the number of XIST RNA domains per cell varied considerably even within a single tumor. Frequent X-chromosome allelic and copy number aberrations were found, in agreement with aberrant XIST RNA domain numbers. In summary, by combining multiple approaches to assess the genetics and epigenetics of a large series of BRCA1 primary tumors, we can conclude definitively that BRCA1 is not required for XIST RNA coating of the X chromosome. The intratumoral and intertumoral variability in XIST RNA domain number in BRCA1 tumors correlates with chromosomal genetic abnormalities, including gains, losses, reduplications, and rearrangements of the X-chromosome. Finally, we also show the necessity for combined global and single-cell approaches in the assessment of tumors with such a high degree of heterogeneity.

VAMP: visualization and analysis of array-CGH, transcriptome and other molecular profiles.

MOTIVATION: Microarray-based CGH (Comparative Genomic Hybridization), transcriptome arrays and other large-scale genomic technologies are now routinely used to generate a vast amount of genomic profiles. Exploratory analysis of this data is crucial in helping to understand the data and to help form biological hypotheses. This step requires visualization of the data in a meaningful way to visualize the results and to perform first level analyses. RESULTS: We have developed a graphical user interface for visualization and first level analysis of molecular profiles. It is currently in use at the Institut Curie for cancer research projects involving CGH arrays, transcriptome arrays, SNP (single nucleotide polymorphism) arrays, loss of heterozygosity results (LOH), and Chromatin ImmunoPrecipitation arrays (ChIP chips). The interface offers the possibility of studying these different types of information in a consistent way. Several views are proposed, such as the classical CGH karyotype view or genome-wide multi-tumor comparison. Many functionalities for analyzing CGH data are provided by the interface, including looking for recurrent regions of alterations, confrontation to transcriptome data or clinical information, and clustering. Our tool consists of PHP scripts and of an applet written in Java. It can be run on public datasets at http://bioinfo.curie.fr/vamp AVAILABILITY: The VAMP software (Visualization and Analysis of array-CGH,transcriptome and other Molecular Profiles) is available upon request. It can be tested on public datasets at http://bioinfo.curie.fr/vamp. The documentation is available at http://bioinfo.curie.fr/vamp/doc.

Spatial normalization of array-CGH data.

BACKGROUND: Array-based comparative genomic hybridization (array-CGH) is a recently developed technique for analyzing changes in DNA copy number. As in all microarray analyses, normalization is required to correct for experimental artifacts while preserving the true biological signal. We investigated various sources of systematic variation in array-CGH data and identified two distinct types of spatial effect of no biological relevance as the predominant experimental artifacts: continuous spatial gradients and local spatial bias. Local spatial bias affects a large proportion of arrays, and has not previously been considered in array-CGH experiments. RESULTS: We show that existing normalization techniques do not correct these spatial effects properly. We therefore developed an automatic method for the spatial normalization of array-CGH data. This method makes it possible to delineate and to eliminate and/or correct areas affected by spatial bias. It is based on the combination of a spatial segmentation algorithm called NEM (Neighborhood Expectation Maximization) and spatial trend estimation. We defined quality criteria for array-CGH data, demonstrating significant improvements in data quality with our method for three data sets coming from two different platforms (198, 175 and 26 BAC-arrays). CONCLUSION: We have designed an automatic algorithm for the spatial normalization of BAC CGH-array data, preventing the misinterpretation of experimental artifacts as biologically relevant outliers in the genomic profile. This algorithm is implemented in the R package MANOR (Micro-Array NORmalization), which is described at http://bioinfo.curie.fr/projects/manor and available from the Bioconductor site http://www.bioconductor.org. It can also be tested on the CAPweb bioinformatics platform at http://bioinfo.curie.fr/CAPweb.