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.

Maternal Copy Number Variants (CNV) transmission to their Autism Spectrum Disorder (ASD) sons correlates with phenotypic traits

Background: Common neurodevelopmental disorder, global prevalence ~1 %; Persistent deficits in social communication and social interaction; restricted and repetitive behavior, interests, or activities; Highly heterogeneous clinical presentation; Male to female ratio ~4:1.

Genome-wide copy number variation (CNV) detection in Nelore cattle reveals highly frequent variants in genome regions harboring QTLs affecting production traits.

Background: Copy number variations (CNVs) have been shown to account for substantial portions of observed genomic variation and have been associated with qualitative and quantitative traits and the onset of disease in a number of species. Information from high-resolution studies to detect, characterize and estimate population-specific variant frequencies will facilitate the incorporation of CNVs in genomic studies to identify genes affecting traits of importance. Results: Genome-wide CNVs were detected in high-density single nucleotide polymorphism (SNP) genotyping data from 1,717 Nelore (Bos indicus) cattle, and in NGS data from eight key ancestral bulls. A total of 68,007 and 12,786 distinct CNVs were observed, respectively. Cross-comparisons of results obtained for the eight resequenced animals revealed that 92 % of the CNVs were observed in both datasets, while 62 % of all detected CNVs were observed to overlap with previously validated cattle copy number variant regions (CNVRs). Observed CNVs were used for obtaining breed-specific CNV frequencies and identification of CNVRs, which were subsequently used for gene annotation. A total of 688 of the detected CNVRs were observed to overlap with 286 non-redundant QTLs associated with important production traits in cattle. All of 34 CNVs previously reported to be associated with milk production traits in Holsteins were also observed in Nelore cattle. Comparisons of estimated frequencies of these CNVs in the two breeds revealed 14, 13, 6 and 14 regions in high (>20 %), low (<20 %) and divergent (NEL > HOL, NEL < HOL) frequencies, respectively. Conclusions: Obtained results significantly enriched the bovine CNV map and enabled the identification of variants that are potentially associated with traits under selection in Nelore cattle, particularly in genome regions harboring QTLs affecting production traits.

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.

Down-regulation Of Slc8a1 As A Putative Apoptosis-evasion Mechanism By Modulation Of Calcium Levels In Penile Carcinoma.

The SLC8A1 gene, which encodes the Na(+)/Ca(2+) exchanger, plays a key role in calcium homeostasis. Our previous gene expression oligoarray data revealed SLC8A1 underexpression in penile carcinoma (PeCa). The aim of this study was to investigate whether the dysregulation of SLC8A1 expression is associated with apoptosis and cell proliferation in PeCa, via modulation of calcium concentration. The underlying mechanisms of SLC8A1 underexpression were also explored, focusing on copy number alteration and microRNA. Transcript levels of SLC8A1 gene and miR-223 were evaluated by quantitative PCR, comparing PeCa samples with normal glans tissues. SLC8A1 copy number was evaluated by microarray-based comparative genomic hybridization (array-CGH). Caspase-3 and Ki-67 immunostaining, as well as calcium distribution by Laser Ablation Imaging Inductively Coupled Plasma Mass Spectrometry [LA(i)-ICP-MS], were investigated in both normal and tumor samples. Confirming our previous data, SLC8A1 underexpression was detected in PeCa samples (P=0.001) and was not associated with gene copy number loss. In contrast, overexpression of miR-223 (P=0.002) was inversely correlated with SLC8A1 (P=0.015, r=-0.426), its putative repressor. In addition, SLC8A1 underexpression was associated with decreased calcium distribution, high Ki-67 and low caspase-3 immunoexpression in PeCa when compared with normal tissues. Down-regulation of the SLC8A1 gene, most likely mediated by its regulator miR-223, can lead to reduced calcium levels in PeCa and, consequently, to suppression of apoptosis and increased tumor cell proliferation. These data suggest that the miR-223-NCX1-calcium-signaling axis may represent a potential therapeutic approach in PeCa.

Chromosomal gains and losses in ocular melanoma detected by comparative genomic hybridization in an Australian population-based study

To define the location of potential oncogenes and tumor suppressor genes in ocular melanoma we carried out comparative genomic hybridization (CGH) analysis on a population-based series of 25 formalin-fixed, paraffin-embedded primary tumors comprising 17 choroidal, 2 ciliary body, 4 iris, and 2 conjunctival melanomas. Twelve (48%) of the 25 melanomas showed no chromosomal changes and 13 (52%) had at least one chromosomal gain or loss. The mean number of CGH changes in all tumors was 3.3, with similar mean numbers of chromosomal gains (1.5) and losses (1.8). The highest number of chromosomal changes (i.e., nine) occurred in a conjunctival melanoma and included four changes not observed in tumors at any other ocular site (gains in 22q and 11p and losses in 6p and 17p). The most frequent gains in all primary ocular melanomas were on chromosome arm 8q (69%), 6p (31%) and 8p (23%) and the most frequent losses were on 6q (38%), 10q (23%), and 16q (23%). The most common pairing was gain in 8p and gain in 8q, implying a whole chromosome copy number increase; gains in 8p occurred only in conjunction with gains in 8q. The smallest regions of copy number alteration were mapped to gain of 8q21 and loss of 6q21, 10q21, and 16q22. Sublocalization of these chromosomal changes to single-band resolution should accelerate the identification of genes involved in the genesis of ocular melanoma.