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.

Analysis of genomic alterations in benign, atypical, and anaplastic meningiomas: Toward a genetic model of meningioma progression

Nineteen benign [World Health Organization (WHO) grade I; MI], 21 atypical (WHO grade II; MII), and 19 anaplastic (WHO grade III; MIII) sporadic meningiomas were screened for chromosomal imbalances by comparative genomic hybridization (CGH). These data were supplemented by molecular genetic analyses of selected chromosomal regions and genes. With increasing malignancy grade, a marked accumulation of genomic aberrations was observed; i.e., the numbers (mean ± SEM) of total alterations detected per tumor were 2.9 ± 0.7 for MI, 9.2 ± 1.2 for MII, and 13.3 ± 1.9 for MIII. The most frequent alteration detected in MI was loss on 22q (58%). In MII, aberrations most commonly identified were losses on 1p (76%), 22q (71%), 14q (43%), 18q (43%), 10 (38%), and 6q (33%), as well as gains on 20q (48%), 12q (43%), 15q (43%), 1q (33%), 9q (33%), and 17q (33%). In MIII, most of these alterations were found at similar frequencies. However, an increase in losses on 6q (53%), 10 (68%), and 14q (63%) was observed. In addition, 32% of MIII demonstrated loss on 9p. Homozygous deletions in the CDKN2A gene at 9p21 were found in 4 of 16 MIII (25%). Highly amplified DNA sequences were mapped to 12q13–q15 by CGH in 1 MII. Southern blot analysis of this tumor revealed amplification of CDK4 and MDM2. By CGH, DNA sequences from 17q were found to be amplified in 1 MII and 8 MIII, involving 17q23 in all cases. Despite the high frequency of chromosomal aberrations in the MII and MIII investigated, none of these tumors showed mutations in exons 5–8 of the TP53 gene. On the basis of the most common aberrations identified in the various malignancy grades, a model for the genomic alterations associated with meningioma progression is proposed.

Chromosomal aberrations in PARP−/− mice: Genome stabilization in immortalized cells by reintroduction of poly(ADP-ribose) polymerase cDNA

Depletion of poly(ADP-ribose) polymerase (PARP) increases the frequency of recombination, gene amplification, sister chromatid exchanges, and micronuclei formation in cells exposed to genotoxic agents, implicating PARP in the maintenance of genomic stability. Flow cytometric analysis now has revealed an unstable tetraploid population in immortalized fibroblasts derived from PARP−/− mice. Comparative genomic hybridization detected partial chromosomal gains in 4C5-ter, 5F-ter, and 14A1-C1 in PARP−/−mice and immortalized PARP−/−fibroblasts. Neither the chromosomal gains nor the tetraploid population were apparent in PARP−/− cells stably transfected with PARP cDNA [PARP−/−(+PARP)], indicating negative selection of cells with these genetic aberrations after reintroduction of PARP cDNA. Although the tumor suppressor p53 was not detectable in PARP−/− cells, p53 expression was partially restored in PARP−/− (+PARP) cells. Loss of 14D3-ter that encompasses the tumor suppressor gene Rb-1 in PARP−/− mice was associated with a reduction in retinoblastoma(Rb) expression; increased expression of the oncogene Jun was correlated with a gain in 4C5-ter that harbors this oncogene. These results further implicate PARP in the maintenance of genomic stability and suggest that altered expression of p53, Rb, and Jun, as well as undoubtedly many other proteins may be a result of genomic instability associated with PARP deficiency.

Spectral karyotyping suggests additional subsets of colorectal cancers characterized by pattern of chromosome rearrangement

The abundant chromosome abnormalities in most carcinomas are probably a reflection of genomic instability present in the tumor, so the pattern and variability of chromosome abnormalities will reflect the mechanism of instability combined with the effects of selection. Chromosome rearrangement was investigated in 17 colorectal carcinoma-derived cell lines. Comparative genomic hybridization showed that the chromosome changes were representative of those found in primary tumors. Spectral karyotyping (SKY) showed that translocations were very varied and mostly unbalanced, with no translocation occurring in more than three lines. At least three karyotype patterns could be distinguished. Some lines had few chromosome abnormalities: they all showed microsatellite instability, the replication error (RER)+ phenotype. Most lines had many chromosome abnormalities: at least seven showed a surprisingly consistent pattern, characterized by multiple unbalanced translocations and intermetaphase variation, with chromosome numbers around triploid, 6–16 structural aberrations, and similarities in gains and losses. Almost all of these were RER−, but one, LS411, was RER+. The line HCA7 showed a novel pattern, suggesting a third kind of genomic instability: multiple reciprocal translocations, with little numerical change or variability. This line was also RER+. The coexistence in one tumor of two kinds of genomic instability is to be expected if the underlying defects are selected for in tumor evolution.

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.

Gain of chromosome 3q defines the transition from severe dysplasia to invasive carcinoma of the uterine cervix.

We have chosen tumors of the uterine cervix as a model system to identify chromosomal aberrations that occur during carcinogenesis. A phenotype/genotype correlation was established in defined regions of archived, formalin-fixed, and hematoxylin/eosin-stained tissue sections that were dissected from normal cervical epithelium (n = 3), from mild (n = 4), moderate (n = 6), and severe dysplasias/carcinomas in situ (CIS) (n = 13), and from invasive carcinomas (n = 10) and investigated by comparative genomic hybridization. The same tissues were analyzed for DNA ploidy, proliferative activity, and the presence of human papillomavirus (HPV) sequences. The results show that an increase in proliferative activity and tetraploidization had occurred already in mildly dysplastic lesions. No recurrent chromosomal aberrations were observed in DNA extracted from normal epithelium or from mild and moderate dysplasias, indicating that the tetraploidization precedes the loss or gain of specific chromosomes. A gain of chromosome 3q became visible in one of the severe dysplasias/CIS. Notably, chromosome 3q was overrepresented in 90% of the carcinomas and was also found to have undergone a high-level copy-number increase (amplification). We therefore conclude that the gain of chromosome 3q that occurs in HPV16-infected, aneuploid cells represents a pivotal genetic aberration at the transition from severe dysplasia/CIS to invasive cervical carcinoma.

Comparative analysis of the gene-dense ACHE/TFR2 region on human chromosome 7q22 with the orthologous region on mouse chromosome 5

Chromosome 7q22 has been the focus of many cytogenetic and molecular studies aimed at delineating regions commonly deleted in myeloid leukemias and myelodysplastic syndromes. We have compared a gene-dense, GC-rich sub-region of 7q22 with the orthologous region on mouse chromosome 5. A physical map of 640 kb of genomic DNA from mouse chromosome 5 was derived from a series of overlapping bacterial artificial chromosomes. A 296 kb segment from the physical map, spanning Ache to Tfr2, was compared with 267 kb of human sequence. We identified a conserved linkage of 12 genes including an open reading frame flanked by Ache and Asr2, a novel cation-chloride cotransporter interacting protein Cip1, Ephb4, Zan and Perq1. While some of these genes have been previously described, in each case we present new data derived from our comparative sequence analysis. Adjacent unfinished sequence data from the mouse contains an orthologous block of 10 additional genes including three novel cDNA sequences that we subsequently mapped to human 7q22. Methods for displaying comparative genomic information, including unfinished sequence data, are becoming increasingly important. We supplement our printed comparative analysis with a new, Web-based program called Laj (local alignments with java). Laj provides interactive access to archived pairwise sequence alignments via the WWW. It displays synchronized views of a dot-plot, a percent identity plot, a nucleotide-level local alignment and a variety of relevant annotations. Our mouse–human comparison can be viewed at http://web.uvic.ca/~bioweb/laj.html. Laj is available at http://bio.cse.psu.edu/, along with online documentation and additional examples of annotated genomic regions.

Genomic evidence for a complete sexual cycle in Candida albicans

Candida albicans is a diploid fungus that has become a medically important opportunistic pathogen in immunocompromised individuals. We have sequenced the C. albicans genome to 10.4-fold coverage and performed a comparative genomic analysis between C. albicans and Saccharomyces cerevisiae with the objective of assessing whether Candida possesses a genetic repertoire that could support a complete sexual cycle. Analyzing over 500 genes important for sexual differentiation in S. cerevisiae, we find many homologues of genes that are implicated in the initiation of meiosis, chromosome recombination, and the formation of synaptonemal complexes. However, others are striking in their absence. C. albicans seems to have homologues of all of the elements of a functional pheromone response pathway involved in mating in S. cerevisiae but lacks many homologues of S. cerevisiae genes for meiosis. Other meiotic gene homologues in organisms ranging from filamentous fungi to Drosophila melanogaster and Caenorhabditis elegans were also found in the C. albicans genome, suggesting potential alternative mechanisms of genetic exchange.

Structure of the gene for porcine peptide antibiotic PR-39, a cathelin gene family member: comparative mapping of the locus for the human peptide antibiotic FALL-39.

PR-39 is a porcine 39-aa peptide antibiotic composed of 49% proline and 24% arginine, with an activity against Gram-negative bacteria comparable to that of tetracycline. In Escherichia coli, it inhibits DNA and protein synthesis. PR-39 was originally isolated from pig small intestine, but subsequent cDNA cloning showed that the gene is expressed in the bone marrow. The open reading frame of the clone showed that PR-39 is made as 173-aa precursor whose proregion belongs to the cathelin family. The PR39 gene, which is rather compact and spans only 1784 bp has now been sequenced. The coding information is split into four exons. The first exon contains the signal sequence of 29 residues and the first 37 residues of the cathelin propart. Exons 2 and 3 contain only cathelin information, while exon 4 codes for the four C-terminal cathelin residues and the mature PR-39 peptide extended by three residues. The sequenced upstream region (1183 bp) contains four potential recognition sites for NF-IL6 and three for APRF, transcription factors known to regulate genes for both cytokines and acute phase response factors. Genomic hybridizations revealed a fairly high level of restriction fragment length polymorphism and indicated that there are at least two copies of the PR39 gene in the pig genome. PR39 was mapped to pig chromosome 13 by linkage and in situ hybridization mapping. The gene for the human peptide antibiotic FALL-39 (also a member of the cathelin family) was mapped to human chromosome 3, which is homologous to pig chromosome 13.

Exploring drug-induced alterations in gene expression in Mycobacterium tuberculosis by microarray hybridization

Tuberculosis is a chronic infectious disease that is transmitted by cough-propelled droplets that carry the etiologic bacterium, Mycobacterium tuberculosis. Although currently available drugs kill most isolates of M. tuberculosis, strains resistant to each of these have emerged, and multiply resistant strains are increasingly widespread. The growing problem of drug resistance combined with a global incidence of seven million new cases per year underscore the urgent need for new antituberculosis therapies. The recent publication of the complete sequence of the M. tuberculosis genome has made possible, for the first time, a comprehensive genomic approach to the biology of this organism and to the drug discovery process. We used a DNA microarray containing 97% of the ORFs predicted from this sequence to monitor changes in M. tuberculosis gene expression in response to the antituberculous drug isoniazid. Here we show that isoniazid induced several genes that encode proteins physiologically relevant to the drug’s mode of action, including an operonic cluster of five genes encoding type II fatty acid synthase enzymes and fbpC, which encodes trehalose dimycolyl transferase. Other genes, not apparently within directly affected biosynthetic pathways, also were induced. These genes, efpA, fadE23, fadE24, and ahpC, likely mediate processes that are linked to the toxic consequences of the drug. Insights gained from this approach may define new drug targets and suggest new methods for identifying compounds that inhibit those targets.

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.

t(11;22)(q23;q11.2) in acute myeloid leukemia of infant twins fuses MLL with hCDCrel, a cell division cycle gene in the genomic region of deletion in DiGeorge and velocardiofacial syndromes

We examined the MLL genomic translocation breakpoint in acute myeloid leukemia of infant twins. Southern blot analysis in both cases showed two identical MLL gene rearrangements indicating chromosomal translocation. The rearrangements were detectable in the second twin before signs of clinical disease and the intensity relative to the normal fragment indicated that the translocation was not constitutional. Fluorescence in situ hybridization with an MLL-specific probe and karyotype analyses suggested t(11;22)(q23;q11.2) disrupting MLL. Known 5′ sequence from MLL but unknown 3′ sequence from chromosome band 22q11.2 formed the breakpoint junction on the der(11) chromosome. We used panhandle variant PCR to clone the translocation breakpoint. By ligating a single-stranded oligonucleotide that was homologous to known 5′ MLL genomic sequence to the 5′ ends of BamHI-digested DNA through a bridging oligonucleotide, we formed the stem–loop template for panhandle variant PCR which yielded products of 3.9 kb. The MLL genomic breakpoint was in intron 7. The sequence of the partner DNA from band 22q11.2 was identical to the hCDCrel (human cell division cycle related) gene that maps to the region commonly deleted in DiGeorge and velocardiofacial syndromes. Both MLL and hCDCrel contained homologous CT, TTTGTG, and GAA sequences within a few base pairs of their respective breakpoints, which may have been important in uniting these two genes by translocation. Reverse transcriptase-PCR amplified an in-frame fusion of MLL exon 7 to hCDCrel exon 3, indicating that an MLL-hCDCrel chimeric mRNA had been transcribed. Panhandle variant PCR is a powerful strategy for cloning translocation breakpoints where the partner gene is undetermined. This application of the method identified a region of chromosome band 22q11.2 involved in both leukemia and a constitutional disorder.

Comparative genetics in the grasses

Genetic mapping of wheat, maize, and rice and other grass species with common DNA probes has revealed remarkable conservation of gene content and gene order over the 60 million years of radiation of Poaceae. The linear organization of genes in some nine different genomes differing in basic chromosome number from 5 to 12 and nuclear DNA amount from 400 to 6,000 Mb, can be described in terms of only 25 “rice linkage blocks.” The extent to which this intergenomic colinearity is confounded at the micro level by gene duplication and micro-rearrangements is still an open question. Nevertheless, it is clear that the elucidation of the organization of the economically important grasses with larger genomes, such as maize (2n = 10, 4,500 Mb DNA), will, to a greater or lesser extent, be predicted from sequence analysis of smaller genomes such as rice, with only 400 Mb, which in turn may be greatly aided by knowledge of the entire sequence of Arabidopsis, which may be available as soon as the turn of the century. Comparative genetics will provide the key to unlock the genomic secrets of crop plants with bigger genomes than Homo sapiens.

The role of genetic and genomic attributes in the success of polyploids

In 1950, G. Ledyard Stebbins devoted two chapters of his book Variation and Evolution in Plants (Columbia Univ. Press, New York) to polyploidy, one on occurrence and nature and one on distribution and significance. Fifty years later, many of the questions Stebbins posed have not been answered, and many new questions have arisen. In this paper, we review some of the genetic attributes of polyploids that have been suggested to account for the tremendous success of polyploid plants. Based on a limited number of studies, we conclude: (i) Polyploids, both individuals and populations, generally maintain higher levels of heterozygosity than do their diploid progenitors. (ii) Polyploids exhibit less inbreeding depression than do their diploid parents and can therefore tolerate higher levels of selfing; polyploid ferns indeed have higher levels of selfing than do their diploid parents, but polyploid angiosperms do not differ in outcrossing rates from their diploid parents. (iii) Most polyploid species are polyphyletic, having formed recurrently from genetically different diploid parents. This mode of formation incorporates genetic diversity from multiple progenitor populations into the polyploid “species”; thus, genetic diversity in polyploid species is much higher than expected by models of polyploid formation involving a single origin. (iv) Genome rearrangement may be a common attribute of polyploids, based on evidence from genome in situ hybridization (GISH), restriction fragment length polymorphism (RFLP) analysis, and chromosome mapping. (v) Several groups of plants may be ancient polyploids, with large regions of homologous DNA. These duplicated genes and genomes can undergo divergent evolution and evolve new functions. These genetic and genomic attributes of polyploids may have both biochemical and ecological benefits that contribute to the success of polyploids in nature.

Genomic analysis of orthologous mouse and human olfactory receptor loci

Olfactory receptor (OR) genes represent ≈1% of genomic coding sequence in mammals, and these genes are clustered on multiple chromosomes in both the mouse and human genomes. We have taken a comparative genomics approach to identify features that may be involved in the dynamic evolution of this gene family and in the transcriptional control that results in a single OR gene expressed per olfactory neuron. We sequenced ≈350 kb of the murine P2 OR cluster and used synteny, gene linkage, and phylogenetic analysis to identify and sequence ≈111 kb of an orthologous cluster in the human genome. In total, 18 mouse and 8 human OR genes were identified, including 7 orthologs that appear to be functional in both species. Noncoding homology is evident between orthologs and generally is confined within the transcriptional unit. We find no evidence for common regulatory features shared among paralogs, and promoter regions generally do not contain strong promoter motifs. We discuss these observations, as well as OR clustering, in the context of evolutionary expansion and transcriptional regulation of OR repertoires.