Clinical relevance of polymorphic DNA copy number variation (CNV) in African American women with stage I-II breast cancer

Objectives. The chief goal of this study was to analyze copy number variation (CNV) in breast cancer tumors from 25 African American women with early stage breast cancer (BC) using molecular inversion probes (MIP) in order to: (1) compare the degree of CNV in tumors compared to normal lymph nodes, and (2) determine whether gains and/or losses of genes in specific chromosomes differ between pathologic subtypes of breast cancer defined by known prognostic markers, (3) determine whether gains/losses in CN are associated with known oncogenes or tumor suppressor genes, and (4) determine whether increased gains/losses in CN for specific chromosomes were associated with differences in breast cancer recurrence. ^ Methods. Twenty to 37 nanograms of DNA extracted from 25 formalin-fixed paraffin embedded (FFPE) tumor samples and matched normal lymph nodes were added to individual tubes. Oligonucleotide probes with recognition sequences at each terminus were hybridized with a genomic target sequence to form a circular structure. Probes are released from genomic DNA obtained from FFPE samples, and those which have been correctly "circularized" in the proper allele/nucleotide reaction combination are amplified using polymerase chain reaction (PCR) primers. Amplicons were fluorescently labeled and the tag sequences released from the genome homology regions by treatment with uracil-N-glycosylase to cleave the probe at the site where uracils are present, and detected using a complementary tag array developed by Affymetrix. ^ Results. Analysis of CN gains and losses from tumors and normal tissues showed marked differences in tumors with numerous chromosomes affected. Similar changes were not observed in normal lymph nodes. When tumors were stratified into four groups based on expression or lack of expression of the estrogen receptor and HER2/neu, distinct patterns of CNV for different chromosomes were observed. Gains or losses in CN for specific chromosomes correlated with amplifications/deletions of particular oncogenes or tumor suppressor genes (i.e. such as found on chromosome 17) known to be associated with aggressive tumor phenotype and poor prognosis. There was a trend for increases in CN observed for chromosome 17 to correlate inversely with time to recurrence of BC (p=0.14 for trend). CNV was also observed for chromosomes 5, 8, 10, 11, and 16, which are known sites for several breast cancer susceptibility alleles. ^ Conclusions. This study is the first to validate the MIP technique, to correlate differences in gene expression with known prognostic tumor markers, and to correlate significant increases/decreases in CN with known tumor markers associated with prognosis. The results of this study may have far reaching public health implications towards identifying new high-risk groups based on genomic differences in CNP, both with respect to prognosis and response to therapy, and to eventually identify new therapeutic targets for prevention and treatment of this disease. ^

Quantifying single gene copy number by measuring fluorescent probe lengths on combed genomic DNA

An approach was developed for the quantification of subtle gains and losses of genomic DNA. The approach relies on a process called molecular combing. Molecular combing consists of the extension and alignment of purified molecules of genomic DNA on a glass coverslip. It has the advantage that a large number of genomes can be combed per coverslip, which allows for a statistically adequate number of measurements to be made on the combed DNA. Consequently, a high-resolution approach to mapping and quantifying genomic alterations is possible. The approach consists of applying fluorescence hybridization to the combed DNA by using probes to identify the amplified region. Measurements then are made on the linear hybridization signals to ascertain the region's exact size. The reliability of the approach first was tested for low copy number amplifications by determining the copy number of chromosome 21 in a normal and trisomy 21 cell line. It then was tested for high copy number amplifications by quantifying the copy number of an oncogene amplified in the tumor cell line GTL-16. These results demonstrate that a wide range of amplifications can be accurately and reliably quantified. The sensitivity and resolution of the approach likewise was assessed by determining the copy number of a single allele (160 kb) alteration.

rrndb: the Ribosomal RNA Operon Copy Number Database

The Ribosomal RNA Operon Copy Number Database (rrndb) is an Internet-accessible database containing annotated information on rRNA operon copy number among prokaryotes. Gene redundancy is uncommon in prokaryotic genomes, yet the rRNA genes can vary from one to as many as 15 copies. Despite the widespread use of 16S rRNA gene sequences for identification of prokaryotes, information on the number and sequence of individual rRNA genes in a genome is not readily accessible. In an attempt to understand the evolutionary implications of rRNA operon redundancy, we have created a phylogenetically arranged report on rRNA gene copy number for a diverse collection of prokaryotic microorganisms. Each entry (organism) in the rrndb contains detailed information linked directly to external websites including the Ribosomal Database Project, GenBank, PubMed and several culture collections. Data contained in the rrndb will be valuable to researchers investigating microbial ecology and evolution using 16S rRNA gene sequences. The rrndb web site is directly accessible on the WWW at http://rrndb.cme.msu.edu.

Simultaneous A8344G heteroplasmy and mitochondrial DNA copy number quantification in Myoclonus Epilepsy and Ragged-Red Fibers (MERRF) syndrome by a multiplex Molecular Beacon based real-time fluorescence PCR

The association of a particular mitochondrial DNA (mtDNA) mutation with different clinical phenotypes is a well-known feature of mitochondrial diseases. A simple genotype–phenotype correlation has not been found between mutation load and disease expression. Tissue and intercellular mosaicism as well as mtDNA copy number are thought to be responsible for the different clinical phenotypes. As disease expression of mitochondrial tRNA mutations is mostly in postmitotic tissues, studies to elucidate disease mechanisms need to be performed on patient material. Heteroplasmy quantitation and copy number estimation using small patient biopsy samples has not been reported before, mainly due to technical restrictions. In order to resolve this problem, we have developed a robust assay that utilizes Molecular Beacons to accurately quantify heteroplasmy levels and determine mtDNA copy number in small samples carrying the A8344G tRNALys mutation. It provides the methodological basis to investigate the role of heteroplasmy and mtDNA copy number in determining the clinical phenotypes.

Comprehensive copy number and gene expression profiling of the 17q23 amplicon in human breast cancer

The biological significance of DNA amplification in cancer is thought to be due to the selection of increased expression of a single or few important genes. However, systematic surveys of the copy number and expression of all genes within an amplified region of the genome have not been performed. Here we have used a combination of molecular, genomic, and microarray technologies to identify target genes for 17q23, a common region of amplification in breast cancers with poor prognosis. Construction of a 4-Mb genomic contig made it possible to define two common regions of amplification in breast cancer cell lines. Analysis of 184 primary breast tumors by fluorescence in situ hybridization on tissue microarrays validated these results with the highest amplification frequency (12.5%) observed for the distal region. Based on GeneMap'99 information, 17 known genes and 26 expressed sequence tags were localized to the contig. Analysis of genomic sequence identified 77 additional transcripts. A comprehensive analysis of expression levels of these transcripts in six breast cancer cell lines was carried out by using complementary DNA microarrays. The expression patterns varied from one cell line to another, and several overexpressed genes were identified. Of these, RPS6KB1, MUL, APPBP2, and TRAP240 as well as one uncharacterized expressed sequence tag were located in the two common amplified regions. In summary, comprehensive analysis of the 17q23 amplicon revealed a limited number of highly expressed genes that may contribute to the more aggressive clinical course observed in breast cancer patients with 17q23-amplified tumors.

Quantitative analysis of chromosomal CGH in human breast tumors associates copy number abnormalities with p53 status and patient survival

We present a general method for rigorously identifying correlations between variations in large-scale molecular profiles and outcomes and apply it to chromosomal comparative genomic hybridization data from a set of 52 breast tumors. We identify two loci where copy number abnormalities are correlated with poor survival outcome (gain at 8q24 and loss at 9q13). We also identify a relationship between abnormalities at two loci and the mutational status of p53. Gain at 8q24 and loss at 5q15-5q21 are linked with mutant p53. The 9q and 5q losses suggest the possibility of gene products involved in breast cancer progression. The analytical techniques are general and also are applicable to the analysis of array-based expression data.

High copy number I-Ab transgenes induce production of IgE through an interluekin 4-dependent mechanism.

To better understand the role of class II major histocompatibility complex molecules in both normal and autoimmune responses, we have produced a series of I-Ab transgenic mice. One of these transgenic constructs, designated NOD.PD, has the sequence of the NOD beta chain (Abeta(g7)) except at positions 56 and 57, where Pro-Asp replaces His-Ser. Several NOD.PD transgenic lines have been produced. One line of these mice carried a very high number of copies (>50) of the NOD.PD transgene. As has been described in other mice carrying high copy numbers of I-Ab transgenes, B-cell development was abnormal. The steady state numbers of mature B cells (IgM+/IgD(hi)) in the periphery were greatly reduced in transgenic mice compared to nontransgenic littermates. Surprisingly, rather than being accompanied by a generalized hypogammaglobulinemia, this B-cell deficiency was accompanied by elevated concentrations of IgG1 and IgE in the serum. Conversely, the levels of IgG2a were reduced in transgenic mice compared to nontransgenic littermates. Because this isotype pattern was characteristic of interleukin (IL)-4-induced class-switching, we then investigated the role of IL-4 in causing the observed phenotype. We crossed the high copy number transgenic mice with an IL-4-deficient strain of mice. As expected, the elevated levels of IgE in high copy number transgenic mice were eliminated when the IL-4 gene was inactivated. However, the reduction in the number of B cells was not ameliorated. These data indicate that the primary defect caused by the transgene was to reduce the number of B cells in these mice. This reduction was accompanied by a secondary increase in IL-4 production, which drove the remaining B cells toward the production of IgGl and IgE.

PCR-Based Analysis of Mitochondrial DNA Copy Number, Mitochondrial DNA Damage, and Nuclear DNA Damage.

Because of the role that DNA damage and depletion play in human disease, it is important to develop and improve tools to assess these endpoints. This unit describes PCR-based methods to measure nuclear and mitochondrial DNA damage and copy number. Long amplicon quantitative polymerase chain reaction (LA-QPCR) is used to detect DNA damage by measuring the number of polymerase-inhibiting lesions present based on the amount of PCR amplification; real-time PCR (RT-PCR) is used to calculate genome content. In this unit, we provide step-by-step instructions to perform these assays in Homo sapiens, Mus musculus, Rattus norvegicus, Caenorhabditis elegans, Drosophila melanogaster, Danio rerio, Oryzias latipes, Fundulus grandis, and Fundulus heteroclitus, and discuss the advantages and disadvantages of these assays.

A compendium of myeloma-associated chromosomal copy number abnormalities and their prognostic value.

To obtain a comprehensive genomic profile of presenting multiple myeloma cases we performed high-resolution single nucleotide polymorphism mapping array analysis in 114 samples alongside 258 samples analyzed by U133 Plus 2.0 expression array (Affymetrix). We examined DNA copy number alterations and loss of heterozygosity (LOH) to define the spectrum of minimally deleted regions in which relevant genes of interest can be found. The most frequent deletions are located at 1p (30%), 6q (33%), 8p (25%), 12p (15%), 13q (59%), 14q (39%), 16q (35%), 17p (7%), 20 (12%), and 22 (18%). In addition, copy number-neutral LOH, or uniparental disomy, was also prevalent on 1q (8%), 16q (9%), and X (20%), and was associated with regions of gain and loss. Based on fluorescence in situ hybridization and expression quartile analysis, genes of prognostic importance were found to be located at 1p (FAF1, CDKN2C), 1q (ANP32E), and 17p (TP53). In addition, we identified common homozygously deleted genes that have functions relevant to myeloma biology. Taken together, these analyses indicate that the crucial pathways in myeloma pathogenesis include the nuclear factor-κB pathway, apoptosis, cell-cycle regulation, Wnt signaling, and histone modifications. This study was registered at http://isrctn.org as ISRCTN68454111.

Small-plasmid-mediated antibiotic resistance is enhanced by increases in plasmid copy number and bacterial fitness

Plasmids play a key role in the horizontal spread of antibiotic resistance determinants among bacterial pathogens. When an antibiotic resistance plasmid arrives in a new bacterial host, it produces a fitness cost, causing a competitive disadvantage for the plasmid-bearing bacterium in the absence of antibiotics. On the other hand, in the presence of antibiotics, the plasmid promotes the survival of the clone. The adaptations experienced by plasmid and bacterium in the presence of antibiotics during the first generations of coexistence will be crucial for the progress of the infection and the maintenance of plasmid-mediated resistance once the treatment is over. Here we developed a model system using the human pathogen Haemophilus influenzae carrying the small plasmid pB1000 conferring resistance to β-lactam antibiotics to investigate host and plasmid adaptations in the course of a simulated ampicillin therapy. Our results proved that plasmid-bearing clones compensated for the fitness disadvantage during the first 100 generations of plasmid-host adaptation. In addition, ampicillin treatment was associated with an increase in pB1000 copy number. The augmentation in both bacterial fitness and plasmid copy number gave rise to H. influenzae populations with higher ampicillin resistance levels. In conclusion, we show here that the modulations in bacterial fitness and plasmid copy number help a plasmid-bearing bacterium to adapt during antibiotic therapy, promoting both the survival of the host and the spread of the plasmid.

TiMMing: developing an innovative suite of bioinformatic tools to harmonize and track the origin of copy number alterations in the evolutive history of multiple myeloma

Multiple Myeloma (MM) is a hematologic cancer with heterogeneous and complex genomic landscape, where Copy Number Alterations (CNAs) play a key role in the disease's pathogenesis and prognosis. It is of biological and clinical interest to study the temporal occurrence of early alterations, as they play a disease "driver" function by deregulating key tumor pathways. This study presents an innovative bioinformatic tools suite created for harmonizing and tracing the origin of CNAs throughout the evolutionary history of MM. To this aim, large cohorts of newly-diagnosed MM (NDMM, N=1582) and Smoldering-MM (SMM, N=282) were aggregated. The tools developed in this study enable the harmonization of CNAs as obtained from different genomic platforms in such a way that a high statistical power can be obtained. By doing so, the high numerosity of those cohorts was harnessed for the identification of novel genes characterized as "driver" (NFKB2, NOTCH2, MAX, EVI5 and MYC-ME2-enhancer), and the generation of an innovative timing model, implemented with a statistical method to introduce confidence intervals in the CNAs-calls. By applying this model on both NDMM and SMM cohorts, it was possible to identify specific CNAs (1q(CKS1B)amp, 13q(RB1)del, 11q(CCND1)amp and 14q(MAX)del) and categorize them as "early"/ "driver" events. A high level of precision was guaranteed by the narrow confidence intervals in the timing estimates. These CNAs were proposed as critical MM alterations, which play a foundational role in the evolutionary history of both SMM and NDMM. Finally, a multivariate survival model was able to identify the independent genomic alterations with the greatest effect on patients’ survival, including RB1-del, CKS1B-amp, MYC-amp, NOTCH2-amp and TRAF3-del/mut. In conclusion, the alterations that were identified as both "early-drivers” and correlated with patients’ survival were proposed as biomarkers that, if included in wider survival models, could provide a better disease stratification and an improved prognosis definition.

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.

The 16p11.2 locus modulates brain structures common to autism, schizophrenia and obesity

Anatomical structures and mechanisms linking genes to neuropsychiatric disorders are not deciphered. Reciprocal copy number variants at the 16p11.2 BP4-BP5 locus offer a unique opportunity to study the intermediate phenotypes in carriers at high risk for autism spectrum disorder (ASD) or schizophrenia (SZ). We investigated the variation in brain anatomy in 16p11.2 deletion and duplication carriers. Beyond gene dosage effects on global brain metrics, we show that the number of genomic copies negatively correlated to the gray matter volume and white matter tissue properties in cortico-subcortical regions implicated in reward, language and social cognition. Despite the near absence of ASD or SZ diagnoses in our 16p11.2 cohort, the pattern of brain anatomy changes in carriers spatially overlaps with the well-established structural abnormalities in ASD and SZ. Using measures of peripheral mRNA levels, we confirm our genomic copy number findings. This combined molecular, neuroimaging and clinical approach, applied to larger datasets, will help interpret the relative contributions of genes to neuropsychiatric conditions by measuring their effect on local brain anatomy.Molecular Psychiatry advance online publication, 25 November 2014; doi:10.1038/mp.2014.145.

A higher mutational burden in females supports a "female protective model" in neurodevelopmental disorders.

Increased male prevalence has been repeatedly reported in several neurodevelopmental disorders (NDs), leading to the concept of a "female protective model." We investigated the molecular basis of this sex-based difference in liability and demonstrated an excess of deleterious autosomal copy-number variants (CNVs) in females compared to males (odds ratio [OR] = 1.46, p = 8 × 10(-10)) in a cohort of 15,585 probands ascertained for NDs. In an independent autism spectrum disorder (ASD) cohort of 762 families, we found a 3-fold increase in deleterious autosomal CNVs (p = 7 × 10(-4)) and an excess of private deleterious single-nucleotide variants (SNVs) in female compared to male probands (OR = 1.34, p = 0.03). We also showed that the deleteriousness of autosomal SNVs was significantly higher in female probands (p = 0.0006). A similar bias was observed in parents of probands ascertained for NDs. Deleterious CNVs (>400 kb) were maternally inherited more often (up to 64%, p = 10(-15)) than small CNVs < 400 kb (OR = 1.45, p = 0.0003). In the ASD cohort, increased maternal transmission was also observed for deleterious CNVs and SNVs. Although ASD females showed higher mutational burden and lower cognition, the excess mutational burden remained, even after adjustment for those cognitive differences. These results strongly suggest that females have an increased etiological burden unlinked to rare deleterious variants on the X chromosome. Carefully phenotyped and genotyped cohorts will be required for identifying the symptoms, which show gender-specific liability to mutational burden.