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.

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.

Modeling stochastic gene expression: Implications for haploinsufficiency

There is increasing recognition that stochastic processes regulate highly predictable patterns of gene expression in developing organisms, but the implications of stochastic gene expression for understanding haploinsufficiency remain largely unexplored. We have used simulations of stochastic gene expression to illustrate that gene copy number and expression deactivation rates are important variables in achieving predictable outcomes. In gene expression systems with non-zero expression deactivation rates, diploid systems had a higher probability of uninterrupted gene expression than haploid systems and were more successful at maintaining gene product above a very low threshold. Systems with relatively rapid expression deactivation rates (unstable gene expression) had more predictable responses to a gradient of inducer than systems with slow or zero expression deactivation rates (stable gene expression), and diploid systems were more predictable than haploid, with or without dosage compensation. We suggest that null mutations of a single allele in a diploid organism could decrease the probability of gene expression and present the hypothesis that some haploinsufficiency syndromes might result from an increased susceptibility to stochastic delays of gene initiation or interruptions of gene expression.

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.

Genetic analysis using genomic representations

Analysis of the genetic changes in human tumors is often problematical because of the presence of normal stroma and the limited availability of pure tumor DNA. However, large amounts of highly reproducible “representations” of tumor and normal genomes can be made by PCR from nanogram amounts of restriction endonuclease cleaved DNA that has been ligated to oligonucleotide adaptors. We show here that representations are useful for many types of genetic analyses, including measuring relative gene copy number, loss of heterozygosity, and comparative genomic hybridization. Representations may be prepared even from sorted nuclei from fixed and archived tumor biopsies.

Phenotypic conversion of drug-resistant bacteria to drug sensitivity

Plasmids that contain synthetic genes coding for small oligoribonucleotides called external guide sequences (EGSs) have been introduced into strains of Escherichia coli harboring antibiotic resistance genes. The EGSs direct RNase P to cleave the mRNAs transcribed from these genes thereby converting the phenotype of drug-resistant cells to drug sensitivity. Increasing the EGS-to-target mRNA ratio by changing gene copy number or the number of EGSs complementary to different target sites enhances the efficiency of the conversion process. We demonstrate a general method for the efficient phenotypic conversion of drug-resistant bacterial cultures.

The embAB genes of Mycobacterium avium encode an arabinosyl transferase involved in cell wall arabinan biosynthesis that is the target for the antimycobacterial drug ethambutol.

The antimycobacterial compound ethambutol [Emb; dextro-2,2'-(ethylenediimino)-di-1-butanol] is used to treat tuberculosis as well as disseminated infections caused by Mycobacterium avium. The critical target for Emb lies in the pathway for the biosynthesis of cell wall arabinogalactan, but the molecular mechanisms for drug action and resistance are unknown. The cellular target for Emb was sought using drug resistance, via target overexpression by a plasmid vector, as a selection tool. This strategy led to the cloning of the M. avium emb region which rendered the otherwise susceptible Mycobacterium smegmatis host resistant to Emb. This region contains three complete open reading frames (ORFs), embR, embA, and embB. The translationally coupled embA and embB genes are necessary and sufficient for an Emb-resistant phenotype which depends on gene copy number, and their putative novel membrane proteins are homologous to each other. The predicted protein encoded by embR, which is related to known transcriptional activators from Streptomyces, is expendable for the phenotypic expression of Emb resistance, but an intact divergent promoter region between embR and embAB is required. An Emb-sensitive cell-free assay for arabinan biosynthesis shows that overexpression of embAB is associated with high-level Emb-resistant arabinosyl transferase activity, and that embR appears to modulate the in vitro level of this activity. These data suggest that embAB encode the drug target of Emb, the arabinosyl transferase responsible for the polymerization of arabinose into the arabinan of arabinogalactan, and that overproduction of this Emb-sensitive target leads to Emb resistance.

The human granulocyte-macrophage colony- stimulating factor gene is autonomously regulated in vivo by an inducible tissue-specific enhancer

The granulocyte-macrophage colony-stimulating factor (GM-CSF) gene is part of a cytokine gene cluster and is directly linked to a conserved upstream inducible enhancer. Here we examined the in vitro and in vivo functions of the human GM-CSF enhancer and found that it was required for the correctly regulated expression of the GM-CSF gene. An inducible DNase I-hypersensitive site appeared within the enhancer in cell types such as T cells, myeloid cells, and endothelial cells that express GM-CSF, but not in nonexpressing cells. In a panel of transfected cells the human GM-CSF enhancer was activated in a tissue-specific manner in parallel with the endogenous gene. The in vivo function of the enhancer was examined in a transgenic mouse model that also addressed the issue of whether the GM-CSF locus was correctly regulated in isolation from other segments of the cytokine gene cluster. After correction for copy number the mean level of human GM-CSF expression in splenocytes from 11 lines of transgenic mice containing a 10.5-kb human GM-CSF transgene was indistinguishable from mouse GM-CSF expression (99% ± 56% SD). In contrast, a 9.8-kb transgene lacking just the enhancer had a significantly reduced (P = 0.004) and more variable level of activity (29% ± 89% SD). From these studies we conclude that the GM-CSF enhancer is required for the correct copy number-dependent expression of the human GM-CSF gene and that the GM-CSF gene is regulated independently from DNA elements associated with the closely linked IL-3 gene or other members of the cytokine gene cluster.

Homology-dependent Gene Silencing in Paramecium

Microinjection at high copy number of plasmids containing only the coding region of a gene into the Paramecium somatic macronucleus led to a marked reduction in the expression of the corresponding endogenous gene(s). The silencing effect, which is stably maintained throughout vegetative growth, has been observed for all Paramecium genes examined so far: a single-copy gene (ND7), as well as members of multigene families (centrin genes and trichocyst matrix protein genes) in which all closely related paralogous genes appeared to be affected. This phenomenon may be related to posttranscriptional gene silencing in transgenic plants and quelling in Neurospora and allows the efficient creation of specific mutant phenotypes thus providing a potentially powerful tool to study gene function in Paramecium. For the two multigene families that encode proteins that coassemble to build up complex subcellular structures the analysis presented herein provides the first experimental evidence that the members of these gene families are not functionally redundant.

Long-term expression of γ-globin mRNA in mouse erythrocytes from retrovirus vectors containing the human γ-globin gene fused to the ankyrin-1 promoter

Gene therapy for patients with hemoglobin disorders has been hampered by the inability of retrovirus vectors to transfer globin genes and their cis-acting regulatory sequences into hematopoietic stem cells without rearrangement. In addition, the expression from intact globin gene vectors has been variable in red blood cells due to position effects and retrovirus silencing. We hypothesized that by substituting the globin gene promoter for the promoter of another gene expressed in red blood cells, we could generate stable retrovirus vectors that would express globin at sufficient levels to treat hemoglobinopathies. Recently, we have shown that the human ankyrin (Ank) gene promoter directs position-independent, copy number-dependent expression of a linked γ-globin gene in transgenic mice. We inserted the Ank/Aγ-globin gene into retrovirus vectors that could transfer one or two copies of the Ank/Aγ-globin gene to target cells. Both vectors were stable, transferring only intact proviral sequences into primary mouse hematopoietic stem cells. Expression of Ank/Aγ-globin mRNA in mature red blood cells was 3% (single copy) and 8% (double copy) of the level of mouse α-globin mRNA. We conclude that these novel retrovirus vectors may be valuable for treating a variety of red cell disorders by gene replacement therapy including severe β-thalassemia if the level of expression can be further increased.

Determination of gene organization in individual haplotypes by analyzing single DNA fragments from single spermatozoa

To determine human Ig heavy chain variable region (VH) gene segment organization on individual homologous chromosomes, an efficient approach has been developed. Single spermatozoa were used as subjects for the study. Upon sperm lysis, VH regions in each sperm were randomly sheared into fragments by the random Brownian force. The fragments were separated from each other by aliquoting the lysate into a certain number of tubes. The gene segments in the VH1 and VH4 families in each tube were identified by denaturing gradient gel electrophoresis after PCR amplification. The polymorphic VH sequences were used to determine the parental origins of the analyzed sperm. VH segment organization in the parental haplotypes was determined by aligning the overlapping fragments from the spermatozoa with the corresponding haplotypes. Based on this comparison between the resulting haplotype maps and the composite map reported previously, the VH region on chromosome 14 could be subdivided into four portions. The numbers and compositions of the VH gene segments differ considerably among the maps in two portions, but are highly conserved in the other two. The data also indicate that the VH region on chromosome 15 may contain a large duplicated block with copy number varying among haplotypes. The approach used in the present study may be used to construct high-resolution haplotype maps without molecular cloning.

Overexpression of 5-methylcytosine DNA glycosylase in human embryonic kidney cells EcR293 demethylates the promoter of a hormone-regulated reporter gene

We have shown that the DNA demethylation complex isolated from chicken embryos has a G⋅T mismatch DNA glycosylase that also possesses 5-methylcytosine DNA glycosylase (5-MCDG) activity. Herein we show that human embryonic kidney cells stably transfected with 5-MCDG cDNA linked to a cytomegalovirus promoter overexpress 5-MCDG. A 15- to 20-fold overexpression of 5-MCDG results in the specific demethylation of a stably integrated ecdysone-retinoic acid responsive enhancer-promoter linked to a β-galactosidase reporter gene. Demethylation occurs in the absence of the ligand ponasterone A (an analogue of ecdysone). The state of methylation of the transgene was investigated by Southern blot analysis and by the bisulfite genomic sequencing reaction. Demethylation occurs downstream of the hormone response elements. No genome-wide demethylation was observed. The expression of an inactive mutant of 5-MCDG or the empty vector does not elicit any demethylation of the promoter-enhancer of the reporter gene. An increase in 5-MCDG activity does not influence the activity of DNA methyltransferase(s) when tested in vitro with a hemimethylated substrate. There is no change in the transgene copy number during selection of the clones with antibiotics. Immunoprecipitation combined with Western blot analysis showed that an antibody directed against 5-MCDG precipitates a complex containing the retinoid X receptor α. The association between retinoid receptor and 5-MCDG is not ligand dependent. These results suggest that a complex of the hormone receptor with 5-MCDG may target demethylation of the transgene in this system.