Comparative genomic hybridization to detect variation in the copy number of large DNA segments.

Array comparative genomic hybridization (CGH) is an excellent tool to scan the genome for copy number variations (CNVs) when used conscientiously. This article is intended to provide an understanding of the basic principles of array CGH and the different options available to the user to design their array CGH experiments. Specifically, the six subsections discuss the different array platforms available, test and reference DNA preparation, reference DNA choice, the basics of hybridization, data processing, and our current understanding of CNVs in the human genome.

Rigorous and thorough bioinformatic analyses of olfactory receptor promoters confirm enrichment of O/E and homeodomain binding sites but reveal no new common motifs.

BACKGROUND: Mammalian olfactory receptors (ORs) are subject to a remarkable but poorly understood regime of transcriptional regulation, whereby individual olfactory neurons each express only one allele of a single member of the large OR gene family. RESULTS: We performed a rigorous search for enriched sequence motifs in the largest dataset of OR promoter regions analyzed to date. We combined measures of cross-species conservation with databases of known transcription factor binding sites and ab initio motif-finding algorithms. We found strong enrichment of binding sites for the O/E family of transcription factors and for homeodomain factors, both already known to be involved in the transcriptional control of ORs, but did not identify any novel enriched sequences. We also found that TATA-boxes are present in at least a subset of OR promoters. CONCLUSIONS: Our rigorous approach provides a template for the analysis of the regulation of large gene families and demonstrates some of the difficulties and pitfalls of such analyses. Although currently available bioinformatics methods cannot detect all transcriptional regulatory elements, our thorough analysis of OR promoters shows that in the case of this gene family, experimental approaches have probably already identified all the binding factors common to large fractions of OR promoters.

Chimpanzee and human Y chromosomes are remarkably divergent in structure and gene content.

The human Y chromosome began to evolve from an autosome hundreds of millions of years ago, acquiring a sex-determining function and undergoing a series of inversions that suppressed crossing over with the X chromosome. Little is known about the recent evolution of the Y chromosome because only the human Y chromosome has been fully sequenced. Prevailing theories hold that Y chromosomes evolve by gene loss, the pace of which slows over time, eventually leading to a paucity of genes, and stasis. These theories have been buttressed by partial sequence data from newly emergent plant and animal Y chromosomes, but they have not been tested in older, highly evolved Y chromosomes such as that of humans. Here we finished sequencing of the male-specific region of the Y chromosome (MSY) in our closest living relative, the chimpanzee, achieving levels of accuracy and completion previously reached for the human MSY. By comparing the MSYs of the two species we show that they differ radically in sequence structure and gene content, indicating rapid evolution during the past 6 million years. The chimpanzee MSY contains twice as many massive palindromes as the human MSY, yet it has lost large fractions of the MSY protein-coding genes and gene families present in the last common ancestor. We suggest that the extraordinary divergence of the chimpanzee and human MSYs was driven by four synergistic factors: the prominent role of the MSY in sperm production, 'genetic hitchhiking' effects in the absence of meiotic crossing over, frequent ectopic recombination within the MSY, and species differences in mating behaviour. Although genetic decay may be the principal dynamic in the evolution of newly emergent Y chromosomes, wholesale renovation is the paramount theme in the continuing evolution of chimpanzee, human and perhaps other older MSYs.

Extreme variability among mammalian V1R gene families.

We report an evolutionary analysis of the V1R gene family across 37 mammalian genomes. V1Rs comprise one of three chemosensory receptor families expressed in the vomeronasal organ, and contribute to pheromone detection. We first demonstrate that Trace Archive data can be used effectively to determine V1R family sizes and to obtain sequences of most V1R family members. Analyses of V1R sequences from trace data and genome assemblies show that species-specific expansions previously observed in only eight species were prevalent throughout mammalian evolution, resulting in "semi-private" V1R repertoires for most mammals. The largest families are found in mouse and platypus, whose V1R repertoires have been published previously, followed by mouse lemur and rabbit (approximately 215 and approximately 160 intact V1Rs, respectively). In contrast, two bat species and dolphin possess no functional V1Rs, only pseudogenes, and suffered inactivating mutations in the vomeronasal signal transduction gene Trpc2. We show that primate V1R decline happened prior to acquisition of trichromatic vision, earlier during evolution than was previously thought. We also show that it is extremely unlikely that decline of the dog V1R repertoire occurred in response to selective pressures imposed by humans during domestication. Functional repertoire sizes in each species correlate roughly with anatomical observations of vomeronasal organ size and quality; however, no single ecological correlate explains the very diverse fates of this gene family in different mammalian genomes. V1Rs provide one of the most extreme examples observed to date of massive gene duplication in some genomes, with loss of all functional genes in other species.

Comparative analyses of chromosome alterations in soft-tissue metastases within and across patients with castration-resistant prostate cancer.

Androgen deprivation is the mainstay of therapy for progressive prostate cancer. Despite initial and dramatic tumor inhibition, most men eventually fail therapy and die of metastatic castration-resistant (CR) disease. Here, we characterize the profound degree of genomic alteration found in CR tumors using array comparative genomic hybridization (array CGH), gene expression arrays, and fluorescence in situ hybridization (FISH). Bycluster analysis, we show that the similarity of the genomic profiles from primary and metastatic tumors is driven by the patient. Using data adjusted for this similarity, we identify numerous high-frequency alterations in the CR tumors, such as 8p loss and chromosome 7 and 8q gain. By integrating array CGH and expression array data, we reveal genes whose correlated values suggest they are relevant to prostate cancer biology. We find alterations that are significantly associated with the metastases of specific organ sites, and others with CR tumors versus the tumors of patients with localized prostate cancer not treated with androgen deprivation. Within the high-frequency sites of loss in CR metastases, we find an overrepresentation of genes involved in cellular lipid metabolism, including PTEN. Finally, using FISH, we verify the presence of a gene fusion between TMPRSS2 and ERG suggested by chromosome 21 deletions detected by array CGH. We find the fusion in 54% of our CR tumors, and 81% of the fusion-positive tumors contain cells with multiple copies of the fusion. Our investigation lays the foundation for a better understanding of and possible therapeutic targets for CR disease, the poorly responsive and final stage of prostate cancer.

Tissue-specific variation in DNA methylation levels along human chromosome 1.

BACKGROUND: DNA methylation is a major epigenetic modification important for regulating gene expression and suppressing spurious transcription. Most methods to scan the genome in different tissues for differentially methylated sites have focused on the methylation of CpGs in CpG islands, which are concentrations of CpGs often associated with gene promoters. RESULTS: Here, we use a methylation profiling strategy that is predominantly responsive to methylation differences outside of CpG islands. The method compares the yield from two samples of size-selected fragments generated by a methylation-sensitive restriction enzyme. We then profile nine different normal tissues from two human donors relative to spleen using a custom array of genomic clones covering the euchromatic portion of human chromosome 1 and representing 8% of the human genome. We observe gross regional differences in methylation states across chromosome 1 between tissues from the same individual, with the most striking differences detected in the comparison of cerebellum and spleen. Profiles of the same tissue from different donors are strikingly similar, as are the profiles of different lobes of the brain. Comparing our results with published gene expression levels, we find that clones exhibiting extreme ratios reflecting low relative methylation are statistically enriched for genes with high expression ratios, and vice versa, in most pairs of tissues examined. CONCLUSION: The varied patterns of methylation differences detected between tissues by our methylation profiling method reinforce the potential functional significance of regional differences in methylation levels outside of CpG islands.

Comparative sequence analysis of primate subtelomeres originating from a chromosome fission event.

Subtelomeres are concentrations of interchromosomal segmental duplications capped by telomeric repeats at the ends of chromosomes. The nature of the segments shared by different sets of human subtelomeres reflects their high rate of recent interchromosomal exchange. Here, we characterize the rearrangements incurred by the 15q subtelomere after it arose from a chromosome fission event in the common ancestor of great apes. We used FISH, sequencing of genomic clones, and PCR to map the breakpoint of this fission and track the fate of flanking sequence in human, chimpanzee, gorilla, orangutan, and macaque genomes. The ancestral locus, a cluster of olfactory receptor (OR) genes, lies internally on macaque chromosome 7. Sequence originating from this fission site is split between the terminus of 15q and the pericentromere of 14q in the great apes. Numerous structural rearrangements, including interstitial deletions and transfers of material to or from other subtelomeres, occurred subsequent to the fission, such that each species has a unique 15q structure and unique collection of ORs derived from the fission locus. The most striking rearrangement involved transfer of at least 200 kb from the fission-site region to the end of chromosome 4q, where much still resides in chimpanzee and gorilla, but not in human. This gross structural difference places the subtelomeric defect underlying facioscapulohumeral muscular dystrophy (FSHD) much closer to the telomere in human 4q than in the hybrid 4q-15q subtelomere of chimpanzee.

Extensive copy-number variation of the human olfactory receptor gene family.

As much as a quarter of the human genome has been reported to vary in copy number between individuals, including regions containing about half of the members of the olfactory receptor (OR) gene family. We have undertaken a detailed study of copy-number variation of ORs to elucidate the selective and mechanistic forces acting on this gene family and the true impact of copy-number variation on human OR repertoires. We argue that the properties of copy-number variants (CNVs) and other sets of large genomic regions violate the assumptions of statistical methods that are commonly used in the assessment of gene enrichment. Using more appropriate methods, we provide evidence that OR enrichment in CNVs is not due to positive selection but is because of OR preponderance in segmentally duplicated regions, which are known to be frequently copy-number variable, and because purifying selection against CNVs is lower in OR-containing regions than in regions containing essential genes. We also combine multiplex ligation-dependent probe amplification (MLPA) and PCR to assay the copy numbers of 37 candidate CNV ORs in a panel of approximately 50 human individuals. We confirm copy-number variation of 18 ORs but find no variation in this human-diversity panel for 16 other ORs, highlighting the caveat that reported intervals often overrepresent true CNVs. The copy-number variation we describe is likely to underpin significant variation in olfactory abilities among human individuals. Finally, we show that both homology-based and homology-independent processes have played a recent role in remodeling the OR family.

Genomic alterations indicate tumor origin and varied metastatic potential of disseminated cells from prostate cancer patients.

Disseminated epithelial cells can be isolated from the bone marrow of a far greater fraction of prostate-cancer patients than the fraction of patients who progress to metastatic disease. To provide a better understanding of these cells, we have characterized their genomic alterations. We first present an array comparative genomic hybridization method capable of detecting genomic changes in the small number of disseminated cells (10-20) that can typically be obtained from bone marrow aspirates of prostate-cancer patients. We show multiple regions of copy-number change, including alterations common in prostate cancer, such as 8p loss, 8q gain, and gain encompassing the androgen-receptor gene on Xq, in the disseminated cell pools from 11 metastatic patients. We found fewer and less striking genomic alterations in the 48 pools of disseminated cells from patients with organ-confined disease. However, we identify changes shared by these samples with their corresponding primary tumors and prostate-cancer alterations reported in the literature, evidence that these cells, like those in advanced disease, are disseminated tumor cells (DTC). We also show that DTCs from patients with advanced and localized disease share several abnormalities, including losses containing cell-adhesion genes and alterations reported to associate with progressive disease. These shared alterations might confer the capability to disseminate or establish secondary disease. Overall, the spectrum of genomic deviations is evidence for metastatic capacity in advanced-disease DTCs and for variation in that capacity in DTCs from localized disease. Our analysis lays the foundation for elucidation of the relationship between DTC genomic alterations and progressive prostate cancer.

A sequence-based survey of the complex structural organization of tumor genomes.

BACKGROUND: The genomes of many epithelial tumors exhibit extensive chromosomal rearrangements. All classes of genome rearrangements can be identified using end sequencing profiling, which relies on paired-end sequencing of cloned tumor genomes. RESULTS: In the present study brain, breast, ovary, and prostate tumors, along with three breast cancer cell lines, were surveyed using end sequencing profiling, yielding the largest available collection of sequence-ready tumor genome breakpoints and providing evidence that some rearrangements may be recurrent. Sequencing and fluorescence in situ hybridization confirmed translocations and complex tumor genome structures that include co-amplification and packaging of disparate genomic loci with associated molecular heterogeneity. Comparison of the tumor genomes suggests recurrent rearrangements. Some are likely to be novel structural polymorphisms, whereas others may be bona fide somatic rearrangements. A recurrent fusion transcript in breast tumors and a constitutional fusion transcript resulting from a segmental duplication were identified. Analysis of end sequences for single nucleotide polymorphisms revealed candidate somatic mutations and an elevated rate of novel single nucleotide polymorphisms in an ovarian tumor. CONCLUSION: These results suggest that the genomes of many epithelial tumors may be far more dynamic and complex than was previously appreciated and that genomic fusions, including fusion transcripts and proteins, may be common, possibly yielding tumor-specific biomarkers and therapeutic targets.

Human subtelomeric WASH genes encode a new subclass of the WASP family.

Subtelomeres are duplication-rich, structurally variable regions of the human genome situated just proximal of telomeres. We report here that the most terminally located human subtelomeric genes encode a previously unrecognized third subclass of the Wiskott-Aldrich Syndrome Protein family, whose known members reorganize the actin cytoskeleton in response to extracellular stimuli. This new subclass, which we call WASH, is evolutionarily conserved in species as diverged as Entamoeba. We demonstrate that WASH is essential in Drosophila. WASH is widely expressed in human tissues, and human WASH protein colocalizes with actin in filopodia and lamellipodia. The VCA domain of human WASH promotes actin polymerization by the Arp2/3 complex in vitro. WASH duplicated to multiple chromosomal ends during primate evolution, with highest copy number reached in humans, whose WASH repertoires vary. Thus, human subtelomeres are not genetic junkyards, and WASH's location in these dynamic regions could have advantageous as well as pathologic consequences.

V2R gene families degenerated in primates, dog and cow, but expanded in opossum.

The V2R genes are expressed in the mammalian vomeronasal organ, and their products are involved in detecting pheromones. Here, we describe the evolution of the V2R gene family. We have found that the human, chimpanzee, macaque, cow and dog V2R gene families have completely degenerated. Each now contains 9-20 pseudogenes but no intact V2R genes. By contrast, opossum has approximately 90 intact V2R genes that mostly arose by duplication after opossum and rodent lineages diverged. One V2R gene subfamily with unusual biology evolved atypically, showing limited expansion in rodents and persistence of a single, albeit sometimes dysfunctional, ortholog in all other species examined.

Elevated rates of sister chromatid exchange at chromosome ends.

Chromosome ends are known hotspots of meiotic recombination and double-strand breaks. We monitored mitotic sister chromatid exchange (SCE) in telomeres and subtelomeres and found that 17% of all SCE occurs in the terminal 0.1% of the chromosome. Telomeres and subtelomeres are significantly enriched for SCEs, exhibiting rates of SCE per basepair that are at least 1,600 and 160 times greater, respectively, than elsewhere in the genome.

Human subtelomeres are hot spots of interchromosomal recombination and segmental duplication.

Human subtelomeres are polymorphic patchworks of interchromosomal segmental duplications at the ends of chromosomes. Here we provide evidence that these patchworks arose recently through repeated translocations between chromosome ends. We assess the relative contribution of the principal mechanisms of ectopic DNA repair to the formation of subtelomeric duplications and find that non-homologous end-joining predominates. Once subtelomeric duplications arise, they are prone to homology-based sequence transfers as shown by the incongruent phylogenetic relationships of neighbouring sections. Interchromosomal recombination of subtelomeres is a potent force for recent change. Cytogenetic and sequence analyses reveal that pieces of the subtelomeric patchwork have changed location and copy number with unprecedented frequency during primate evolution. Half of the known subtelomeric sequence has formed recently, through human-specific sequence transfers and duplications. Subtelomeric dynamics result in a gene duplication rate significantly higher than the genome average and could have both advantageous and pathological consequences in human biology. More generally, our analyses suggest an evolutionary cycle between segmental polymorphisms and genome rearrangements.

Interchromosomal segmental duplications explain the unusual structure of PRSS3, the gene for an inhibitor-resistant trypsinogen.

Homo sapiens possess several trypsinogen or trypsinogen-like genes of which three (PRSS1, PRSS2, and PRSS3) produce functional trypsins in the digestive tract. PRSS1 and PRSS2 are located on chromosome 7q35, while PRSS3 is found on chromosome 9p13. Here, we report a variation of the theme of new gene creation by duplication: the PRSS3 gene was formed by segmental duplications originating from chromosomes 7q35 and 11q24. As a result, PRSS3 transcripts display two variants of exon 1. The PRSS3 transcript whose gene organization most resembles PRSS1 and PRSS2 encodes a functional protein originally named mesotrypsinogen. The other variant is a fusion transcript, called trypsinogen IV. We show that the first exon of trypsinogen IV is derived from the noncoding first exon of LOC120224, a chromosome 11 gene. LOC120224 codes for a widely conserved transmembrane protein of unknown function. Comparative analyses suggest that these interchromosomal duplications occurred after the divergence of Old World monkeys and hominids. PRSS3 transcripts consist of a mixed population of mRNAs, some expressed in the pancreas and encoding an apparently functional trypsinogen and others of unknown function expressed in brain and a variety of other tissues. Analysis of the selection pressures acting on the trypsinogen gene family shows that, while the apparently functional genes are under mild to strong purifying selection overall, a few residues appear under positive selection. These residues could be involved in interactions with inhibitors.

Divergent V1R repertoires in five species: Amplification in rodents, decimation in primates, and a surprisingly small repertoire in dogs.

The V1R gene family comprises one of two types of putative pheromone receptors expressed in the mammalian vomeronasal organ (VNO). We searched the most recent mouse, rat, dog, chimpanzee, and human genome sequence assemblies to compile a near-complete repertoire of V1R genes for each species. Dog, human, and chimpanzee have very few intact V1Rs (8, 2, and 0, respectively) compared to more than a hundred intact V1Rs in each of the rat (106) and mouse (165) genomes. We also provide the first description of the diversity of V1R pseudogenes in these species. We identify at least 165 pseudogenes in mouse, 110 in rat, 102 in chimpanzee, 115 in human, and 54 in dog. Primate and dog pseudogenes are distributed among almost all V1R subfamilies seen in rodents, indicating that the common ancestor of these species had a diverse V1R repertoire. We find that V1R genes were subject to strikingly different fates in different species and in different subfamilies. In rodents, some subfamilies remained relatively stable or underwent roughly equivalent expansion in mouse and rat; other subfamilies expanded in one species but not the other. The small number of intact V1Rs in the dog genome is unexpected given the presumption that dogs, like rodents, have a functional VNO, and a complex system of pheromone-based behaviors. We identify an intact transient receptor potential channel 2beta in the dog genome, consistent with a functional VNO in dogs. The diminished V1R repertoire in dogs raises questions about the relative contributions of V1Rs versus other candidate pheromone receptor genes in the establishment of complex pheromone systems in mammals.

Array comparative genomic hybridization analysis of genomic alterations in breast cancer subtypes.

In this study, we performed high-resolution array comparative genomic hybridization with an array of 4153 bacterial artificial chromosome clones to assess copy number changes in 44 archival breast cancers. The tumors were flow sorted to exclude non-tumor DNA and increase our ability to detect gene copy number changes. In these tumors, losses were more frequent than gains, and gains in 1q and loss in 16q were the most frequent alterations. We compared gene copy number changes in the tumors based on histologic subtype and estrogen receptor (ER) status, i.e., ER-negative infiltrating ductal carcinoma, ER-positive infiltrating ductal carcinoma, and ER-positive infiltrating lobular carcinoma. We observed a consistent association between loss in regions of 5q and ER-negative infiltrating ductal carcinoma, as well as more frequent loss in 4p16, 8p23, 8p21, 10q25, and 17p11.2 in ER-negative infiltrating ductal carcinoma compared with ER-positive infiltrating ductal carcinoma (adjusted P values < or = 0.05). We also observed high-level amplifications in ER-negative infiltrating ductal carcinoma in regions of 8q24 and 17q12 encompassing the c-myc and c-erbB-2 genes and apparent homozygous deletions in 3p21, 5q33, 8p23, 8p21, 9q34, 16q24, and 19q13. ER-positive infiltrating ductal carcinoma showed a higher frequency of gain in 16p13 and loss in 16q21 than ER-negative infiltrating ductal carcinoma. Correlation analysis highlighted regions of change commonly seen together in ER-negative infiltrating ductal carcinoma. ER-positive infiltrating lobular carcinoma differed from ER-positive infiltrating ductal carcinoma in the frequency of gain in 1q and loss in 11q and showed high-level amplifications in 1q32, 8p23, 11q13, and 11q14. These results indicate that array comparative genomic hybridization can identify significant differences in the genomic alterations between subtypes of breast cancer.

Species specificity in rodent pheromone receptor repertoires.

The mouse V1R putative pheromone receptor gene family consists of at least 137 intact genes clustered at multiple chromosomal locations in the genome. Species-specific pheromone receptor repertoires may partly explain species-specific social behavior. We conducted a genomic analysis of an orthologous pair of mouse and rat V1R gene clusters to test for species specificity in rodent pheromone systems. Mouse and rat have lineage-specific V1R repertoires in each of three major subfamilies at these loci as a result of postspeciation duplications, gene loss, and gene conversions. The onset of this diversification roughly coincides with a wave of Line1 (L1) retrotranspositions into the two loci. We propose that L1 activity has facilitated postspeciation V1R duplications and gene conversions. In addition, we find extensive homology among putative V1R promoter regions in both species. We propose a regulatory model in which promoter homogenization could ensure that V1R genes are equally competitive for a limiting transcriptional structure to account for mutually exclusive V1R expression in vomeronasal neurons.

Odorant receptor expressed sequence tags demonstrate olfactory expression of over 400 genes, extensive alternate splicing and unequal expression levels.

BACKGROUND: The olfactory receptor gene family is one of the largest in the mammalian genome. Previous computational analyses have identified approximately 1,500 mouse olfactory receptors, but experimental evidence confirming olfactory function is available for very few olfactory receptors. We therefore screened a mouse olfactory epithelium cDNA library to obtain olfactory receptor expressed sequence tags, providing evidence of olfactory function for many additional olfactory receptors, as well as identifying gene structure and putative promoter regions. RESULTS: We identified more than 1,200 odorant receptor cDNAs representing more than 400 genes. Using real-time PCR to confirm expression level differences suggested by our screen, we find that transcript levels in the olfactory epithelium can differ between olfactory receptors by up to 300-fold. Differences for one gene pair are apparently due to both unequal numbers of expressing cells and unequal transcript levels per expressing cell. At least two-thirds of olfactory receptors exhibit multiple transcriptional variants, with alternative isoforms of both 5' and 3' untranslated regions. Some transcripts (5%) utilize splice sites within the coding region, contrary to the stereotyped olfactory receptor gene structure. Most atypical transcripts encode nonfunctional olfactory receptors, but can occasionally increase receptor diversity. CONCLUSIONS: Our cDNA collection confirms olfactory function of over one-third of the intact mouse olfactory receptors. Most of these genes were previously annotated as olfactory receptors based solely on sequence similarity. Our finding that different olfactory receptors have different expression levels is intriguing given the one-neuron, one-gene expression regime of olfactory receptors. We provide 5' untranslated region sequences and candidate promoter regions for more than 300 olfactory receptors, valuable resources for computational regulatory motif searches and for designing olfactory receptor microarrays and other experimental probes.

Complex evolution of 7E olfactory receptor genes in segmental duplications.

Large segmental duplications (SDs) constitute at least 3.6% of the human genome and have increased its size, complexity, and diversity. SDs can mediate ectopic sequence exchange resulting in gross chromosomal rearrangements that could contribute to speciation and disease. We have identified and evaluated a subset of human SDs that harbor an 88-member subfamily of olfactory receptor (OR)-like genes called the 7Es. At least 92% of these genes appear to be pseudogenes when compared to other OR genes. The 7E-containing SDs (7E SDs) have duplicated to at least 35 regions of the genome via intra- and interchromosomal duplication events. In contrast to many human SDs, the 7E SDs are not biased towards pericentromeric or subtelomeric regions. We find evidence for gene conversion among 7E genes and larger sequence exchange between 7E SDs, supporting the hypothesis that long, highly similar stretches of DNA facilitate ectopic interactions. The complex structure and history of the 7E SDs necessitates extension of the current model of large-scale DNA duplication. Despite their appearance as pseudogenes, some 7E genes exhibit a signature of purifying selection, and at least one 7E gene is expressed.