Identification of outcome-oriented cut-offs for copy number alterations in multiple myeloma: predictive biomarkers with improved prognostic accuracy

Aim of the present study was to develop a statistical approach to define the best cut-off Copy number alterations (CNAs) calling from genomic data provided by high throughput experiments, able to predict a specific clinical end-point (early relapse, 18 months) in the context of Multiple Myeloma (MM). 743 newly diagnosed MM patients with SNPs array-derived genomic and clinical data were included in the study. CNAs were called both by a conventional (classic, CL) and an outcome-oriented (OO) method, and Progression Free Survival (PFS) hazard ratios of CNAs called by the two approaches were compared. The OO approach successfully identified patients at higher risk of relapse and the univariate survival analysis showed stronger prognostic effects for OO-defined high-risk alterations, as compared to that defined by CL approach, statistically significant for 12 CNAs. Overall, 155/743 patients relapsed within 18 months from the therapy start. A small number of OO-defined CNAs were significantly recurrent in early-relapsed patients (ER-CNAs) - amp1q, amp2p, del2p, del12p, del17p, del19p -. Two groups of patients were identified either carrying or not ≥1 ER-CNAs (249 vs. 494, respectively), the first one with significantly shorter PFS and overall survivals (OS) (PFS HR 2.15, p<0001; OS HR 2.37, p<0.0001). The risk of relapse defined by the presence of ≥1 ER-CNAs was independent from those conferred both by R-IIS 3 (HR=1.51; p=0.01) and by low quality (< stable disease) clinical response (HR=2.59 p=0.004). Notably, the type of induction therapy was not descriptive, suggesting that ER is strongly related to patients’ baseline genomic architecture. In conclusion, the OO- approach employed allowed to define CNAs-specific dynamic clonality cut-offs, improving the CNAs calls’ accuracy to identify MM patients with the highest probability to ER. As being outcome-dependent, the OO-approach is dynamic and might be adjusted according to the selected outcome variable of interest.

High frequency of submicroscopic chromosomal imbalances in patients with syndromic craniosynostosis detected by a combined approach of microsatellite segregation analysis, multiplex ligation-dependent probe amplification and array-based comparative genome hybridisation

We present the first comprehensive study, to our knowledge, on genomic chromosomal analysis in syndromic craniosynostosis. In total, 45 patients with craniosynostotic disorders were screened with a variety of methods including conventional karyotype, microsatellite segregation analysis, subtelomeric multiplex ligation-dependent probe amplification) and whole-genome array-based comparative genome hybridisation. Causative abnormalities were present in 42.2% (19/45) of the samples, and 27.8% (10/36) of the patients with normal conventional karyotype carried submicroscopic imbalances. Our results include a wide variety of imbalances and point to novel chromosomal regions associated with craniosynostosis. The high incidence of pure duplications or trisomies suggests that these are important mechanisms in craniosynostosis, particularly in cases involving the metopic suture.

Chromosome imbalances in syndromic hearing loss

The cause of hearing impairment has not been elucidated in a large proportion of patients. We screened by 1-Mb array-based comparative genomic hybridization (aCGH) 29 individuals with syndromic hearing impairment whose clinical features were not typical of known disorders. Rare chromosomal copy number changes were detected in eight patients, four de novo imbalances and four inherited from a normal parent. The de novo alterations define candidate chromosome segments likely to harbor dosage-sensitive genes related to hearing impairment, namely 1q23.3-q25.2, 2q22q23, 6p25.3 and 11q13.2-q13.4. The rare imbalances also present in normal parents might be casually associated with hearing impairment, but its role as a predisposition gene remains a possibility. Our results show that syndromic deafness is frequently associated with chromosome microimbalances (14-27%), and the use of aCGH for defining disease etiology is recommended.

Nonrecurrent MECP2 duplications mediated by genomic architecture-driven DNA breaks and break-induced replication repair

Recurrent submicroscopic genomic copy number changes are the result of nonallelic homologous recombination (NAHR). Nonrecurrent aberrations, however, can result from different nonexclusive recombination-repair mechanisms. We previously described small microduplications at Xq28 containing MECP2 in four male patients with a severe neurological phenotype. Here, we report on the fine-mapping and breakpoint analysis of 16 unique microduplications. The size of the overlapping copy number changes varies between 0.3 and 2.3 Mb, and FISH analysis on three patients demonstrated a tandem orientation. Although eight of the 32 breakpoint regions coincide with low-copy repeats, none of the duplications are the result of NAHR. Bioinformatics analysis of the breakpoint regions demonstrated a 2.5-fold higher frequency of Alu interspersed repeats as compared with control regions, as well as a very high GC content (53%). Unexpectedly, we obtained the junction in only one patient by long-range PCR, which revealed nonhomologous end joining as the mechanism. Breakpoint analysis in two other patients by inverse PCR and subsequent array comparative genomic hybridization analysis demonstrated the presence of a second duplicated region more telomeric at Xq28, of which one copy was inserted in between the duplicated MECP2 regions. These data suggest a two-step mechanism in which part of Xq28 is first inserted near the MECP2 locus, followed by breakage-induced replication with strand invasion of the normal sister chromatid. Our results indicate that the mechanism by which copy number changes occur in regions with a complex genomic architecture can yield complex rearrangements.

Genomic imbalances associated with mullerian aplasia

Background: Aplasia of the mullerian ducts leads to absence of the uterine corpus, uterine cervix, and upper (superior) vagina. Patients with mullerian aplasia (MA) often exhibit additional clinical features such as renal, vertebral and cardiac defects. A number of different syndromes have been associated with MA, and in most cases its aetiology remains poorly understood. Objective and methods: 14 syndromic patients with MA and 46, XX G-banded karyotype were screened for DNA copy number changes by similar to 1 Mb whole genome bacterial artificial chromosome (BAC) array based comparative genomic hybridisation (CGH). The detected alterations were validated by an independent method and further mapped by high resolution oligo-arrays. Results: Submicroscopic genomic imbalances affecting the 1q21.1, 17q12, 22q11.21, and Xq21.31 chromosome regions were detected in four probands. Presence of the alterations in the normal mother of one patient suggests incomplete penetrance and/or variable expressivity. Conclusion: 4 of the 14 patients (29%) were found to have cryptic genomic alterations. The imbalances on 22q11.21 support recent findings by us and others that alterations in this chromosome region may result in impairment of mullerian duct development. The remaining imbalances indicate involvement of previously unknown chromosome regions in MA, and point specifically to LHX1 and KLHL4 as candidate genes.

Transmissable Dog Cancer Genome Reveals the Origin and History of an Ancient Cell Lineage

Canine transmissible venereal tumor (CTVT) is the oldest known somatic cell lineage. It is a transmissible cancer that propagates naturally in dogs. We sequenced the genomes of two CTVT tumors and found that CTVT has acquired 1.9 million somatic substitution mutations and bears evidence of exposure to ultraviolet light. CTVT is remarkably stable and lacks subclonal heterogeneity despite thousands of rearrangements, copy-number changes, and retrotransposon insertions. More than 10,000 genes carry nonsynonymous variants, and 646 genes have been lost. CTVT first arose in a dog with low genomic heterozygosity that may have lived about 11,000 years ago. The cancer spawned by this individual dispersed across continents about 500 years ago. Our results provide a genetic identikit of an ancient dog and demonstrate the robustness of mammalian somatic cells to survive for millennia despite a massive mutation burden.

Frequent downregulation of the transcription factor Foxa2 in lung cancer through epigenetic silencing

Purpose: We sought to determine the mechanisms of downregulation of the airway transcription factor Foxa2 in lung cancer and the expression status of Foxa2 in non-small-cell lung cancer (NSCLC). Methods: A series of 25 lung cancer cell lines were evaluated for Foxa2 protein expression, FOXA2 mRNA levels, FOXA2 mutations, FOXA2 copy number changes and for evidence of FOXA2 promoter hypermethylation. In addition, 32 NSCLCs were sequenced for FOXA2 mutations and 173 primary NSCLC tumors evaluated for Foxa2 expression using an immunohistochemical assay. Results: Out of the 25 cell lines, 13 (52%) had undetectable FOXA2 mRNA. The expression of FOXA2 mRNA and Foxa2 protein were congruent in 19/22 cells (p = 0.001). FOXA2 mutations were not identified in primary NSCLCs and were infrequent in cell lines. Focal or broad chromosomal deletions involving FOXA2 were not present. The promoter region of FOXA2 had evidence of hypermethylation, with an inverse correlation between FOXA2 mRNA expression and presence of CpG dinucleotide methylation (p < 0.0001). In primary NSCLC tumor specimens, there was a high frequency of either absence (42/173, 24.2%) or no/low expression (96/173,55.4%) of Foxa2. In 130 patients with stage I NSCLC there was a trend towards decreased survival in tumors with no/low expression of Foxa2 (HR of 1.6, 95%CI 0.9-3.1; p = 0.122). Conclusions: Loss of expression of Foxa2 is frequent in lung cancer cell lines and NSCLCs. The main mechanism of downregulation of Foxa2 is epigenetic silencing through promoter hypermethylation. Further elucidation of the involvement of Foxa2 and other airway transcription factors in the pathogenesis of lung cancer may identify novel therapeutic targets. (C) 2012 Elsevier Ireland Ltd. All rights reserved.

Li-Fraumeni-like syndrome associated with a large BRCA1 intragenic deletion

Background: Li-Fraumeni (LFS) and Li-Fraumeni-like (LFL) syndromes are associated to germline TP53 mutations, and are characterized by the development of central nervous system tumors, sarcomas, adrenocortical carcinomas, and other early-onset tumors. Due to the high frequency of breast cancer in LFS/LFL families, these syndromes clinically overlap with hereditary breast cancer (HBC). Germline point mutations in BRCA1, BRCA2, and TP53 genes are associated with high risk of breast cancer. Large rearrangements involving these genes are also implicated in the HBC phenotype. Methods: We have screened DNA copy number changes by MLPA on BRCA1, BRCA2, and TP53 genes in 23 breast cancer patients with a clinical diagnosis consistent with LFS/LFL; most of these families also met the clinical criteria for other HBC syndromes. Results: We found no DNA copy number alterations in the BRCA2 and TP53 genes, but we detected in one patient a 36.4 Kb BRCA1 microdeletion, confirmed and further mapped by array-CGH, encompassing exons 9-19. Breakpoints sequencing analysis suggests that this rearrangement was mediated by flanking Alu sequences. Conclusion: This is the first description of a germline intragenic BRCA1 deletion in a breast cancer patient with a family history consistent with both LFL and HBC syndromes. Our results show that large rearrangements in these known cancer predisposition genes occur, but are not a frequent cause of cancer susceptibility.

Duplications of BHLHA9 are associated with ectrodactyly and tibia hemimelia inherited in non-Mendelian fashion

Background Split-hand/foot malformation (SHFM)-also known as ectrodactyly-is a congenital disorder characterised by severe malformations of the distal limbs affecting the central rays of hands and/or feet. A distinct entity termed SHFLD presents with SHFM and long bone deficiency. Mouse models suggest that a defect of the central apical ectodermal ridge leads to the phenotype. Although six different loci/mutations (SHFM1-6) have been associated with SHFM, the underlying cause in a large number of cases is still unresolved. Methods High resolution array comparative genomic hybridisation (CGH) was performed in patients with SHFLD to detect copy number changes. Candidate genes were further evaluated for expression and function during limb development by whole mount in situ hybridisation and morpholino knock-down experiments. Results Array CGH showed microduplications on chromosome 17p13.3, a locus previously associated with SHFLD. Detailed analysis of 17 families revealed that this copy number variation serves as a susceptibility factor for a highly variable phenotype with reduced penetrance, particularly in females. Compared to other known causes for SHFLD 17p duplications appear to be the most frequent cause of SHFLD. A similar to 11.8 kb minimal critical region was identified encompassing a single gene, BHLHA9, a putative basic loop helix transcription factor. Whole mount in situ hybridisation showed expression restricted to the limb bud mesenchyme underlying the apical ectodermal ridge in mouse and zebrafish embryos. Knock down of bhlha9 in zebrafish resulted in shortening of the pectoral fins. Conclusions Genomic duplications encompassing BHLHA9 are associated with SHFLD and non-Mendelian inheritance characterised by a high degree of non-penetrance with sex bias. Knock-down of bhlha9 in zebrafish causes severe reduction defects of the pectoral fin, indicating a role for this gene in limb development.

Genome-wide DNA profiling of marginal zone lymphomas identifies subtype-specific lesions with an impact on the clinical outcome

Marginal zone B-cell lymphomas (MZLs) have been divided into 3 distinct subtypes (extranodal MZLs of mucosa-associated lymphoid tissue [MALT] type, nodal MZLs, and splenic MZLs). Nevertheless, the relationship between the subtypes is still unclear. We performed a comprehensive analysis of genomic DNA copy number changes in a very large series of MZL cases with the aim of addressing this question. Samples from 218 MZL patients (25 nodal, 57 MALT, 134 splenic, and 2 not better specified MZLs) were analyzed with the Affymetrix Human Mapping 250K SNP arrays, and the data combined with matched gene expression in 33 of 218 cases. MALT lymphoma presented significantly more frequently gains at 3p, 6p, 18p, and del(6q23) (TNFAIP3/A20), whereas splenic MZLs was associated with del(7q31), del(8p). Nodal MZLs did not show statistically significant differences compared with MALT lymphoma while lacking the splenic MZLs-related 7q losses. Gains of 3q and 18q were common to all 3 subtypes. del(8p) was often present together with del(17p) (TP53). Although del(17p) did not determine a worse outcome and del(8p) was only of borderline significance, the presence of both deletions had a highly significant negative impact on the outcome of splenic MZLs.

In vivo function of mutated spliced leader RNAs in Caenorhabditis elegans

The role of spliced leader RNA (SL RNA) in trans-splicing in Caenorhabditis elegans has been studied through a combination of in vitro mutagenesis and in vivo complementation of rrs-1 mutant nematodes, which lack endogenous SL1 RNA. Three classes of mutant SL1 RNAs have been found—those that rescue the lethal phenotype at low concentration of transforming DNA, those that rescue at high but not low concentration, and those that do not rescue at all. These studies showed that some mutations in the otherwise highly conserved 22-nt spliced leader are tolerated for splicing and post-splicing events. A longer spliced leader also can be tolerated but only when present in high copy number. Changes in the first 16 nucleotides result in the appearance of no SL RNA, consistent with the in vitro studies by others showing that the SL1 RNA promoter partly resides within the spliced leader sequence.

The molecular genetics of breast cancer: The contribution of comparative genomic hybridization

Comparative genomic hybridization (CGH) has been the technique of choice over the last 10 years for mapping DNA copy number changes in human tumors. Here we review the literature to demonstrate how CGH has contributed to the comprehension of molecular aspects of breast tumorigenesis. At least two distinct molecular pathways of breast cancer have been characterized that show a strong correlation with histological grade. It seems that grade I invasive ductal carcinomas (IDCs) arise from well-differentiated ductal carcinoma in situ (DCIS), whereas grade III IDCs come from poorly differentiated DCIS. In addition, dedifferentiation from a low- to a high-grade breast cancer has proven an unlikely phenomenon. CGH has been instrumental in dissecting distinct molecular pathways toward breast malignancy and in establishing a direct relationship between genotype and clinical pathological features. (C) 2005 Elsevier GrnbH. All rights reserved.

Unlocking pathology archives for molecular genetic studies: a reliable method to generate probes for chromogenic and fluorescent in situ hybridization

Chromogenic (CISH) and fluorescent ( FISH) in situ hybridization have emerged as reliable techniques to identify amplifications and chromosomal translocations. CISH provides a spatial distribution of gene copy number changes in tumour tissue and allows a direct correlation between copy number changes and the morphological features of neoplastic cells. However, the limited number of commercially available gene probes has hindered the use of this technique. We have devised a protocol to generate probes for CISH that can be applied to formalin-fixed, paraffin-embedded tissue sections (FFPETS). Bacterial artificial chromosomes ( BACs) containing fragments of human DNA which map to specific genomic regions of interest are amplified with phi 29 polymerase and random primer labelled with biotin. The genomic location of these can be readily confirmed by BAC end pair sequencing and FISH mapping on normal lymphocyte metaphase spreads. To demonstrate the reliability of the probes generated with this protocol, four strategies were employed: (i) probes mapping to cyclin D1 (CCND1) were generated and their performance was compared with that of a commercially available probe for the same gene in a series of 10 FFPETS of breast cancer samples of which five harboured CCND1 amplification; (ii) probes targeting cyclin-dependent kinase 4 were used to validate an amplification identified by microarray-based comparative genomic hybridization (aCGH) in a pleomorphic adenoma; (iii) probes targeting fibroblast growth factor receptor 1 and CCND1 were used to validate amplifications mapping to these regions, as defined by aCGH, in an invasive lobular breast carcinoma with FISH and CISH; and (iv) gene-specific probes for ETV6 and NTRK3 were used to demonstrate the presence of t(12; 15)(p12; q25) translocation in a case of breast secretory carcinoma with dual colour FISH. In summary, this protocol enables the generation of probes mapping to any gene of interest that can be applied to FFPETS, allowing correlation of morphological features with gene copy number.

Microarray-based comparative genomic hybridisation of breast cancer patients receiving neoadjuvant chemotherapy

We analysed the molecular genetic profiles of breast cancer samples before and after neoadjuvant chemotherapy with combination doxorubicin and cyclophosphamide (AC). DNA was obtained from microdissected frozen breast core biopsies from 44 patients before chemotherapy. Additional samples were obtained before the second course of chemotherapy (D21) and after the completion of the treatment (surgical specimens) in 17 and 21 patients, respectively. Microarray-based comparative genome hybridisation was performed using a platform containing approx5800 bacterial artificial chromosome clones (genome-wide resolution: 0.9 Mb). Analysis of the 44 pretreatment biopsies revealed that losses of 4p, 4q, 5q, 12q13.11–12q13.12, 17p11.2 and 17q11.2; and gains of 1p, 2p, 7q, 9p, 11q, 19p and 19q were significantly associated with oestrogen receptor negativity. 16q21–q22.1 losses were associated with lobular and 8q24 gains with ductal types. Losses of 5q33.3–q4 and 18p11.31 and gains of 6p25.1–p25.2 and Xp11.4 were associated with HER2 amplification. No correlations between DNA copy number changes and clinical response to AC were found. Microarray-based comparative genome hybridisation analysis of matched pretreatment and D21 biopsies failed to identify statistically significant differences, whereas a comparison between matched pretreatment and surgical samples revealed a statistically significant acquired copy number gain on 11p15.2–11p15.5. The modest chemotherapy-driven genomic changes, despite profound loss of cell numbers, suggest that there is little therapeutic selection of resistant non-modal cell lineages.

Characterising the TP53-deleted subgroup of chronic lymphocytic leukemia: an analysis of additional cytogenetic abnormalities detected by interphase fluorescence in situ hybridisation and array-based comparative genomic hybridisation.

Deletion of the TP53 gene on chromosome 17p13.1 is the prognostic factor associated with the shortest survival in CLL. We used array-based comparative genomic hybridisation (arrayCGH) to identify additional DNA copy number changes in peripheral blood samples from 74 LRF CLL4 trial patients, 37 with >or=5% and 37 without TP53-deleted cells. ArrayCGH reliably detected deletions on 17p, including the TP53 locus, in cases with >or=50%TP53-deleted cells detected by fluorescence in situ hybridisation, plus seven additional cases with deleted regions on 17p excluding TP53. Losses on chromosomal regions 18p and/or 20p were found exclusively in cases with >or=5%TP53-deleted cells (por=5%TP53-deleted cases (p=0.02). In particular, amplification of 2p and deletion of 6q were both more frequent. Cases with >20%TP53-deleted cells had the worst prognosis in the LRF CLL4 trial.