Fluorescence in situ hybridization analysis of keratinocyte growth factor gene amplification and dispersion in evolution of great apes and humans

[Functional analysis of retroviral endogenous inserts in the human genome evolution]

Retroelements, mobile elements produced in DNA by reverse transcription, comprise about 40% of the human genome. A small part of these elements appeared in the genome quite recently after the divergence of humans and chimpanzees had occurred. Evolutionarily young retroelements are represented by the members of four groups, SVA, Alu, L1, and the endogenous HERV-K (HML-2) virus. These retroelements could play a functional role in the course of the molecular evolution of human DNA. We comprehensively studied the contribution of human-specific endogenous viruses (hsERV) to the structural modifications and regulation of the human genome. We found that hsERV presented in 134 copies occupied about 330 000 bp of human DNA. They added to genomic sequences the copies of 50 functional retroviral genes as well as 134 potential promoters and enhancers, 50% of which are located in the regions adjacent to known genes, and 22% in gene introns. At least 67% of these elements are human-specific promoters in vivo. hsERV viruses regulate the activity of known protein-encoding genes by means of RNA interference, function as enhancers, and provide new polyadenylation signals for mRNA.

Identification of human specific gene duplications relative to other primates by array CGH and quantitative PCR

In order to identify human lineage specific (HLS) copy number differences (CNDs) compared to other primates, we performed pair wise comparisons (human vs. chimpanzee, gorilla and orangutan) by using cDNA array comparative genomic hybridization (CGH). A set of 23 genes with HLS duplications were identified, as well as other lineage differences in gene copy number specific of chimpanzee, gorilla and orangutan. Each species has gained more copies of specific genes rather than losing gene copies. Eleven of the 23 genes have only been observed to have undergone HLS duplication in Fortna et al. (2004) and in the present study. Then, seven of these 11 genes were analyzed by quantitative PCR in chimpanzee, gorilla and orangutan, as well as in other six primate species (Hylobates lar, Cercopithecus aethiops, Papio hamadryas, Macaca mulatta, Lagothrix lagothricha, and Saimiri sciureus). Six genes confirmed array CGH data, and four of them appeared to have bona fide HLS duplications (ABCB10, E2F6, CDH12, and TDG genes). We propose that these gene duplications have a potential to contribute to specific human phenotypes.

High divergence in primate-specific duplicated regions: human and chimpanzee chorionic gonadotropin beta genes

Background

Low nucleotide divergence between human and chimpanzee does not sufficiently explain the species-specific morphological, physiological and behavioral traits. As gene duplication is a major prerequisite for the emergence of new genes and novel biological processes, comparative studies of human and chimpanzee duplicated genes may assist in understanding the mechanisms behind primate evolution. We addressed the divergence between human and chimpanzee duplicated genomic regions by using Luteinizing Hormone Beta (LHB)/Chorionic Gonadotropin Beta (CGB) gene cluster as a model. The placental CGB genes that are essential for implantation have evolved from an ancestral pituitary LHB gene by duplications in the primate lineage.

Results

We shotgun sequenced and compared the human (45,165 bp) and chimpanzee (39,876 bp) LHB/CGB regions and hereby present evidence for structural variation resulting in discordant number of CGB genes (6 in human, 5 in chimpanzee). The scenario of species-specific parallel duplications was supported (i) as the most parsimonious solution requiring the least rearrangement events to explain the interspecies structural differences; (ii) by the phylogenetic trees constructed with fragments of intergenic regions; (iii) by the sequence similarity calculations. Across the orthologous regions of LHB/CGB cluster, substitutions and indels contributed approximately equally to the interspecies divergence and the distribution of nucleotide identity was correlated with the regional repeat content. Intraspecies gene conversion may have shaped the LHB/CGB gene cluster. The substitution divergence (1.8–2.59%) exceeded two-three fold the estimates for single-copy loci and the fraction of transversional mutations was increased compared to the unique sequences (43% versus ~30%). Despite the high sequence identity among LHB/CGB genes, there are signs of functional differentiation among the gene copies. Estimates for d_n/d_s rate ratio suggested a purifying selection on LHB and CGB8, and a positive evolution of CGB1.

Conclusion

If generalized, our data suggests that in addition to species-specific deletions and duplications, parallel duplication events may have contributed to genetic differences separating humans from their closest relatives. Compared to unique genomic segments, duplicated regions are characterized by high divergence promoted by intraspecies gene conversion and species-specific chromosomal rearrangements, including the alterations in gene copy number.

Understanding the recent evolution of the human genome: insights from human-chimpanzee genome comparisons

The sequencing of the chimpanzee genome and the comparison with its human counterpart have begun to reveal the spectrum of genetic changes that has accompanied human evolution. In addition to gross karyotypic rearrangements such as the fusion that formed human chromosome 2 and the human-specific pericentric inversions of chromosomes 1 and 18, there is considerable submicroscopic structural variation involving deletions, duplications, and inversions. Lineage-specific segmental duplications, detected by array comparative genomic hybridization and direct sequence comparison, have made a very significant contribution to this structural divergence, which is at least three-fold greater than that due to nucleotide substitutions. Since structural genomic changes may have given rise to irreversible functional differences between the diverging species, their detailed analysis could help to identify the biological processes that have accompanied speciation. To this end, interspecies comparisons have revealed numerous human-specific gains and losses of genes as well as changes in gene expression. The very considerable structural diversity (polymorphism) evident within both lineages has, however, hampered the analysis of the structural divergence between the human and chimpanzee genomes. The concomitant evaluation of genetic divergence and diversity at the nucleotide level has nevertheless served to identify many genes that have evolved under positive selection and may thus have been involved in the development of human lineage-specific traits. Genes that display signs of weak negative selection have also been identified and could represent candidate loci for complex genomic disorders. Here, we review recent progress in comparing the human and chimpanzee genomes and discuss how the differences detected have improved our understanding of the evolution of the human genome.

Duplication and extensive remodeling shaped POTE family genes encoding proteins containing ankyrin repeat and coiled coil domains

The POTE family genes encode a highly homologous group of primate-specific proteins that contain ankyrin repeats and coiled coil domains. At least 13 paralogous POTE family genes are found on 8 human chromosomes (2, 8, 13, 14, 15, 18, 21 and 22), which can be sorted into 3 groups based on sequence similarity. We identified by a database search a group of additional human ankyrin repeat domain proteins, of which ANKRD26 and ANKRD30A are the best characterized; these are more distant homologs of POTE family proteins. A comprehensive comparison of the genomic organization indicates that ANKRD26 has the genomic structure of the possible ancestor of ANKRD30A and all POTE family genes. Extensive remodeling involving segmental loss and internal duplication appears to have reshaped the ANKRD30A and POTE family genes after the primal duplication of the ancestor gene. We also identified a mouse homolog of human ANKRD26, but failed to find a mouse homolog that bears the structural characteristics of any of the POTE family of proteins. The mouse Ankrd26 may serve as a useful model for the study of the function of human ANKRD26, ANKRD30A and POTE family proteins.

An evolutionary history of the FGF superfamily

Fibroblast growth factors (FGF) are associated with multiple developmental and metabolic processes in triploblasts, and perhaps also in diploblasts. The evolution of the FGF superfamily has accompanied the major morphological and functional innovations of metazoan species. The study of FGFs throughout species shows that the FGF superfamily can be subdivided in eight families in present-day organisms and has evolved through phases of gene duplications and gene losses. At least two major expansions of the superfamily can be recognized: a first expansion increased the number of FGFs from one or few archeo-FGFs to eight proto-FGFs, prototypic of the eight families. A second expansion, which took place during euchordate evolution, is associated with genome duplications. It increased the number of members in the families. Subsequent losses reduced that number to the present-day figures.

Punctuated duplication seeding events during the evolution of human chromosome 2p11

Primate genomic sequence comparisons are becoming increasingly useful for elucidating the evolutionary history and organization of our own genome. Such studies are particularly informative within human pericentromeric regions--areas of particularly rapid change in genomic structure. Here, we present a systematic analysis of the evolutionary history of one approximately 700-kb region of 2p11, including the first autosomal transition from pericentromeric sequence to higher-order alpha-satellite DNA. We show that this region is composed of segmental duplications corresponding to 14 ancestral segments ranging in size from 4 kb to approximately 115 kb. These duplicons show 94%-98.5% sequence identity to their ancestral loci. Comparative FISH and phylogenetic analysis indicate that these duplicons are differentially distributed in human, chimpanzee, and gorilla genomes, whereas baboon has a single putative ancestral locus for all but one of the duplications. Our analysis supports a model where duplicative transposition events occurred during a narrow window of evolution after the separation of the human/ape lineage from the Old World monkeys (10-20 million years ago). Although dramatic secondary dispersal events occurred during the radiation of the human, chimpanzee, and gorilla lineages, duplicative transposition seeding events of new material to this particular pericentromeric region abruptly ceased after this time period. The multiplicity of initial duplicative transpositions prior to the separation of humans and great-apes suggests a punctuated model for the formation of highly duplicated pericentromeric regions within the human genome. The data further indicate that factors other than sequence are important determinants for such bursts of duplicative transposition from the euchromatin to pericentromeric regions.

Molecular evolution of the human chromosome 15 pericentromeric region

We present a detailed molecular evolutionary analysis of 1.2 Mb from the pericentromeric region of human 15q11. Sequence analysis indicates the region has been subject to extensive interchromosomal and intrachromosomal duplications during primate evolution. Comparative FISH analyses among non-human primates show remarkable quantitative and qualitative differences in the organization and duplication history of this region - including lineage-specific deletions and duplication expansions. Phylogenetic and comparative analyses reveal that the region is composed of at least 24 distinct segmental duplications or duplicons that have populated the pericentromeric regions of the human genome over the last 40 million years of human evolution. The value of combining both cytogenetic and experimental data in understanding the complex forces which have shaped these regions is discussed.

Keratinocyte growth factor/fibroblast growth factor 7, a homeostatic factor with therapeutic potential for epithelial protection and repair

Keratinocyte growth factor (KGF) is a paracrine-acting, epithelial mitogen produced by cells of mesenchymal origin. It is a member of the fibroblast growth factor (FGF) family, and acts exclusively through a subset of FGF receptor isoforms (FGFR2b) expressed predominantly by epithelial cells. The upregulation of KGF after epithelial injury suggested it had an important role in tissue repair. This hypothesis was reinforced by evidence that intestinal damage was worse and healing impaired in KGF null mice. Preclinical data from several animal models demonstrated that recombinant human KGF could enhance the regenerative capacity of epithelial tissues and protect them from a variety of toxic exposures. These beneficial effects are attributed to multiple mechanisms that collectively act to strengthen the integrity of the epithelial barrier, and include the stimulation of cell proliferation, migration, differentiation, survival, DNA repair, and induction of enzymes involved in the detoxification of reactive oxygen species. KGF is currently being evaluated in clinical trials to test its ability to ameliorate severe oral mucositis (OM) that results from cancer chemoradiotherapy. In a phase 3 trial involving patients who were treated with myeloablative chemoradiotherapy before autologous peripheral blood progenitor cell transplantation for hematologic malignancies, KGF significantly reduced both the incidence and duration of severe OM. Similar investigations are underway in patients being treated for solid tumors. On the basis of its success in ameliorating chemoradiotherapy-induced OM in humans and tissue damage in a variety of animal models, additional clinical applications of KGF are worthy of investigation.

Characterization of the segmental duplication LCR7-20 in the human genome

Our previous study described the amplification of a genomic sequence containing exon 9 of CFTR in the human genome. Here we report that this CFTR sequence is part of a large duplicated sequence unit, provisionally named LCR7-20. Through successive screening of two human chromosome 7-specific cosmid libraries to construct a cosmid contig, we assembled two sequenced BAC clones into a single contig containing a prototypic LCR7-20 unit. Subsequent searches of existing human genome sequences identified additional six copies of LCR7-20-like sequences with more than 90% sequence homology. Additional genomic clones containing LCR7-20-like sequences were then isolated from total genomic BAC and PAC libraries. Restriction fragment analysis and limited sequencing data indicated that there could be around 30 copies of LCR7-20-like sequences in the human genome and that the average region of homology could extend over 120 kb. As indicated by fluorescence in situ hybridization analysis, LCR7-20-like sequences are dispersed on different chromosomes, mainly in the centromeric and pericentromeric regions, and some may exist in tandem copies. Our study also indicates that many genomic regions containing LCR7-20's either have been misassembled or are missing in current versions of the human genome sequence.

Lineage-specific gene duplication and loss in human and great ape evolution

Given that gene duplication is a major driving force of evolutionary change and the key mechanism underlying the emergence of new genes and biological processes, this study sought to use a novel genome-wide approach to identify genes that have undergone lineage-specific duplications or contractions among several hominoid lineages. Interspecies cDNA array-based comparative genomic hybridization was used to individually compare copy number variation for 39,711 cDNAs, representing 29,619 human genes, across five hominoid species, including human. We identified 1,005 genes, either as isolated genes or in clusters positionally biased toward rearrangement-prone genomic regions, that produced relative hybridization signals unique to one or more of the hominoid lineages. Measured as a function of the evolutionary age of each lineage, genes showing copy number expansions were most pronounced in human (134) and include a number of genes thought to be involved in the structure and function of the brain. This work represents, to our knowledge, the first genome-wide gene-based survey of gene duplication across hominoid species. The genes identified here likely represent a significant majority of the major gene copy number changes that have occurred over the past 15 million years of human and great ape evolution and are likely to underlie some of the key phenotypic characteristics that distinguish these species.

This genome-wide analysis reports the major lineage-specific gene copy number changes that have occurred over the past 15 million years of human and great ape evolution

Juxtacentromeric region of human chromosome 21: a boundary between centromeric heterochromatin and euchromatic chromosome arms

We have analysed the genomic structure and transcriptional activity of a 2.3-Mb genomic sequence in the juxtacentromeric region of human chromosome 21. Our work shows that this region comprises two different chromosome domains. The 1.5-Mb proximal domain: (i) is a patchwork of chromosome duplications; (ii) shares sequence similarity with several chromosomes; (iii) contains several gene fragments (truncated genes having an intron/exon structure) intermingled with retrotransposed pseudogenes; and (iv) harbours two genes (TPTE and BAGE2) that belong to gene families and have a cancer and/or testis expression profile. The TPTE gene family was generated before the branching of Old World monkeys from the great ape lineage, by intra- and interchromosome duplications of the ancestral TPTE gene mapping to phylogenetic chromosome XIII. By contrast, the 0.8-Mb distal domain: (i) is devoid of chromosome duplications; (ii) has a chromosome 21-specific sequence; (iii) contains no gene fragments and only one retrotransposed pseudogene; and (iv) harbours six genes including housekeeping genes. G-rich sequences commonly associated with duplication termini cluster at the boundary between the two chromosome domains. These structural and transcriptional features lead us to suggest that the proximal domain has heterochromatic properties, whereas the distal domain has euchromatic properties.

POTE, a highly homologous gene family located on numerous chromosomes and expressed in prostate, ovary, testis, placenta, and prostate cancer

We have identified a gene located on chromosomes 21 that is expressed in normal and neoplastic prostate, and in normal testis, ovary, and placenta. We name this gene POTE (expressed in prostate, ovary, testis, and placenta). The POTE gene has 11 exons and 10 introns and spans approximately equal 32 kb of chromosome 21q11.2 region. The 1.83-kb mRNA of POTE encodes a protein of 66 kDa. Ten paralogs of the gene have been found dispersed among eight different chromosomes (2, 8, 13, 14, 15, 18, 21, and 22) with preservation of ORFs and splice junctions. The synonymous:nonsynonymous ratio indicates that the genes were duplicated rather recently but are diverging at a rate faster than the average for other paralogous genes. In prostate and in testis, at least five different paralogs are expressed. In situ hybridization shows that POTE is expressed in basal and terminal cells of normal prostate epithelium. It is also expressed in some prostate cancers and in the LnCAP prostate cancer cell line. The POTE protein contains seven ankyrin repeats between amino acids 140 and 380. Expression of POTE in prostate cancer and its undetectable expression in normal essential tissues make POTE a candidate for the immunotherapy of prostate cancer. The existence of a large number of closely related but rapidly diverging members, their location on multiple chromosomes and their limited expression pattern suggest an important role for the POTE gene family in reproductive processes.

Organisation of the pericentromeric region of chromosome 15: at least four partial gene copies are amplified in patients with a proximal duplication of 15q

Clinical cytogenetic laboratories frequently identify an apparent duplication of proximal 15q that does not involve probes within the PWS/AS critical region and is not associated with any consistent phenotype. Previous mapping data placed several pseudogenes, NF1, IgH D/V, and GABRA5 in the pericentromeric region of proximal 15q. Recent studies have shown that these pseudogene sequences have increased copy numbers in subjects with apparent duplications of proximal 15q. To determine the extent of variation in a control population, we analysed NF1 and IgH D pseudogene copy number in interphase nuclei from 20 cytogenetically normal subjects by FISH. Both loci are polymorphic in controls, ranging from 1-4 signals for NF1 and 1-3 signals for IgH D. Eight subjects with apparent duplications, examined by the same method, showed significantly increased NF1 copy number (5-10 signals). IgH D copy number was also increased in 6/8 of these patients (4-9 signals). We identified a fourth pseudogene, BCL8A, which maps to the pericentromeric region and is coamplified along with the NF1 sequences. Interphase FISH ordering experiments show that IgH D lies closest to the centromere, while BCL8A is the most distal locus in this pseudogene array; the total size of the amplicon is estimated at approximately 1 Mb. The duplicated chromosome was inherited from either sex parent, indicating no parent of origin effect, and no consistent phenotype was present. FISH analysis with one or more of these probes is therefore useful in discriminating polymorphic amplification of proximal pseudogene sequences from clinically significant duplications of 15q.

The nucleotide sequences of the parathyroid gene in primates (suborder Anthropoidea)

Nucleotide sequences of the parathyroid (PTH) gene of 12 species of primates belonging to suborder Anthropoidea were examined. The PTH gene contains one intron that separates two exons that code the sequence of prepro and PTH, respectively. The intron of the PTH gene in Cebus apella, Callithrix jacchus, and Saguinus oedipus was 102 bp long, whereas a 103-bp intron was observed in the remaining species. Phylogenetic analysis using the nucleotide sequences of PTH revealed that these 12 species of primates of suborder Anthropoidea could be divided into two groups of the infraorder Platyrrhini (C. apella, C. jacchus, and S. oedipus) and the infraorder Catarrhini (Macaca fascicularis, Macaca fuscata, Cercopithecus aethiops, Papio hamadryas, Presbytes obscura, Hylobates lar, Pongo pygmaeus, Pan troglodytes, and Pan paniscus). The latter infraorder could be further subdivided into two subgroups belonging to the superfamily Cercopithecoidea (M. fascicularis, M. fuscata, C. aethiops, P. hamadryas, and P. obscura) and the superfamily Hominoidea (H. lar, P. pygmaeus, P. troglodytes, and P. paniscus). The deduced amino acid sequences of PTH gene between 12 species of nonhuman primates and human revealed no amino acid substitution in mature PTH among orangutans, chimpanzees, and humans. The results indicated that the PTH gene is very conserved among primates, especially between great apes and humans. The apes are the most suitable animals to be used for studying the bone metabolism and applying the knowledge to clinical use in humans.

Human paralogs of KIAA0187 were created through independent pericentromeric-directed and chromosome-specific duplication mechanisms

KIAA0187 is a gene of unknown function that maps to 10q11 and has been subject to recent duplication events. Here we analyze 18 human paralogs of this gene and show that paralogs of exons 14-23 were formed through satellite-associated pericentromeric-directed duplication, whereas paralogs of exons 1-9 were created via chromosome-specific satellite-independent duplications. In silico, Northern, and RT-PCR analyses indicate that nine paralogs are transcribed, including four in which KIAA0187 exons are spliced onto novel sequences. Despite this, no new genes appear to have been created by these events. The chromosome 10 paralogs map to 10q11, 10q22, 10q23.1, and 10q23.3, forming part of a complex family of chromosome-specific repeats that includes GLUD1, Cathepsin L, and KIAA1099 pseudogenes. Phylogenetic analyses and comparative FISH indicates that the 10q23.1 and 10q23.3 repeats were created in 10q11 and relocated by a paracentric inversion 13 to 27 Myr ago. Furthermore, the most recent duplications, involving the KIAA1099 pseudogenes, have largely been confined to 10q11. These results indicate a simple model for the evolution of this repeat family, involving multiple rounds of centromere-proximal duplication and dispersal through intrachromosomal rearrangement. However, more complex events must be invoked to account for high sequence identity between some paralogs.

Human-specific duplication and mosaic transcripts: the recent paralogous structure of chromosome 22

In recent decades, comparative chromosomal banding, chromosome painting, and gene-order studies have shown strong conservation of gross chromosome structure and gene order in mammals. However, findings from the human genome sequence suggest an unprecedented degree of recent (<35 million years ago) segmental duplication. This dynamism of segmental duplications has important implications in disease and evolution. Here we present a chromosome-wide view of the structure and evolution of the most highly homologous duplications (> or = 1 kb and > or = 90%) on chromosome 22. Overall, 10.8% (3.7/33.8 Mb) of chromosome 22 is duplicated, with an average sequence identity of 95.4%. To organize the duplications into tractable units, intron-exon structure and well-defined duplication boundaries were used to define 78 duplicated modules (minimally shared evolutionary segments) with 157 copies on chromosome 22. Analysis of these modules provides evidence for the creation or modification of 11 novel transcripts. Comparative FISH analyses of human, chimpanzee, gorilla, orangutan, and macaque reveal qualitative and quantitative differences in the distribution of these duplications--consistent with their recent origin. Several duplications appear to be human specific, including a approximately 400-kb duplication (99.4%-99.8% sequence identity) that transposed from chromosome 14 to the most proximal pericentromeric region of chromosome 22. Experimental and in silico data further support a pericentromeric gradient of duplications where the most recent duplications transpose adjacent to the centromere. Taken together, these data suggest that segmental duplications have been an ongoing process of primate genome evolution, contributing to recent gene innovation and the dynamic transformation of genome architecture within and among closely related species.

Recent duplication, domain accretion and the dynamic mutation of the human genome

An estimated 5% of the human genome consists of interspersed duplications that have arisen over the past 35 million years of evolution. Two categories of such recently duplicated segments can be distinguished: segmental duplications between nonhomologous chromosomes (transchromosomal duplications) and duplications mainly restricted to a particular chromosome (chromosome-specific duplications). Many of these duplications exhibit an extraordinarily high degree of sequence identity at the nucleotide level (>95%) and span large genomic distances (1-100 kb). Preliminary analyses indicate that these same regions are targets for rapid evolutionary turnover among the genomes of closely related primates. The dynamic nature of these regions because of recurrent chromosomal rearrangement, and their ability to create fusion genes from juxtaposed cassettes suggest that duplicative transposition was an important force in the evolution of our genome.

Segmental duplications: an 'expanding' role in genomic instability and disease

The knowledge that specific genetic diseases are caused by recurrent chromosomal aberrations has indicated that genomic instability might be directly related to the structure of the regions involved. The sequencing of the human genome has directed significant attention towards understanding the molecular basis of such recombination 'hot spots'. Segmental duplications have emerged as a significant factor in the aetiology of disorders that are caused by abnormal gene dosage. These observations bring us closer to understanding the mechanisms and consequences of genomic rearrangement.

Segmental duplications: organization and impact within the current human genome project assembly

Segmental duplications play fundamental roles in both genomic disease and gene evolution. To understand their organization within the human genome, we have developed the computational tools and methods necessary to detect identity between long stretches of genomic sequence despite the presence of high copy repeats and large insertion-deletions. Here we present our analysis of the most recent genome assembly (January 2001) in which we focus on the global organization of these segments and the role they play in the whole-genome assembly process. Initially, we considered only large recent duplication events that fell well-below levels of draft sequencing error (alignments 90%-98% similar and > or =1 kb in length). Duplications (90%-98%; > or =1 kb) comprise 3.6% of all human sequence. These duplications show clustering and up to 10-fold enrichment within pericentromeric and subtelomeric regions. In terms of assembly, duplicated sequences were found to be over-represented in unordered and unassigned contigs indicating that duplicated sequences are difficult to assign to their proper position. To assess coverage of these regions within the genome, we selected BACs containing interchromosomal duplications and characterized their duplication pattern by FISH. Only 47% (106/224) of chromosomes positive by FISH had a corresponding chromosomal position by comparison. We present data that indicate that this is attributable to misassembly, misassignment, and/or decreased sequencing coverage within duplicated regions. Surprisingly, if we consider putative duplications >98% identity, we identify 10.6% (286 Mb) of the current assembly as paralogous. The majority of these alignments, we believe, represent unmerged overlaps within unique regions. Taken together the above data indicate that segmental duplications represent a significant impediment to accurate human genome assembly, requiring the development of specialized techniques to finish these exceptional regions of the genome. The identification and characterization of these highly duplicated regions represents an important step in the complete sequencing of a human reference genome.

Segmental Duplications: Organization and Impact Within the Current Human Genome Project Assembly

Segmental duplications play fundamental roles in both genomic disease and gene evolution. To understand their organization within the human genome, we have developed the computational tools and methods necessary to detect identity between long stretches of genomic sequence despite the presence of high copy repeats and large insertion-deletions. Here we present our analysis of the most recent genome assembly (January 2001) in which we focus on the global organization of these segments and the role they play in the whole-genome assembly process. Initially, we considered only large recent duplication events that fell well-below levels of draft sequencing error (alignments 90%–98% similar and ≥1 kb in length). Duplications (90%–98%; ≥1 kb) comprise 3.6% of all human sequence. These duplications show clustering and up to 10-fold enrichment within pericentromeric and subtelomeric regions. In terms of assembly, duplicated sequences were found to be over-represented in unordered and unassigned contigs indicating that duplicated sequences are difficult to assign to their proper position. To assess coverage of these regions within the genome, we selected BACs containing interchromosomal duplications and characterized their duplication pattern by FISH. Only 47% (106/224) of chromosomes positive by FISH had a corresponding chromosomal position by BLAST comparison. We present data that indicate that this is attributable to misassembly, misassignment, and/or decreased sequencing coverage within duplicated regions. Surprisingly, if we consider putative duplications >98% identity, we identify 10.6% (286 Mb) of the current assembly as paralogous. The majority of these alignments, we believe, represent unmerged overlaps within unique regions. Taken together the above data indicate that segmental duplications represent a significant impediment to accurate human genome assembly, requiring the development of specialized techniques to finish these exceptional regions of the genome. The identification and characterization of these highly duplicated regions represents an important step in the complete sequencing of a human reference genome.

Initial sequencing and analysis of the human genome

The human genome holds an extraordinary trove of information about human development, physiology, medicine and evolution. Here we report the results of an international collaboration to produce and make freely available a draft sequence of the human genome. We also present an initial analysis of the data, describing some of the insights that can be gleaned from the sequence.

Genetic differences between humans and great apes

The remarkable similarity among the genomes of humans and the African great apes could warrant their classification together as a single genus. However, whereas there are many similarities in the biology, life history, and behavior of humans and great apes, there are also many striking differences that need to be explained. The complete sequencing of the human genome creates an opportunity to ask which genes are involved in those differences. A logical approach would be to use the chimpanzee genome for comparison and the other great ape genomes for confirmation. Until such a great ape genome project can become reality, the next best approach must be educated guesses of where the genetic differences may lie and a careful analysis of differences that we do know about. Our group recently discovered a human-specific inactivating mutation in the CMP-sialic acid hydroxylase gene, which results in the loss of expression of a common mammalian cell-surface sugar throughout all cells in the human body. We are currently investigating the implications of this difference for a variety of issues relevant to humans, ranging from pathogen susceptibility to brain development. Evaluating the uniqueness of this finding has also led us to explore the existing literature on the broader issue of genetic differences between humans and great apes. The aim of this brief review is to consider a listing of currently known genetic differences between humans and great apes and to suggest avenues for future research. The differences reported between human and great ape genomes include cytogenetic differences, differences in the type and number of repetitive genomic DNA and transposable elements, abundance and distribution of endogenous retroviruses, the presence and extent of allelic polymorphisms, specific gene inactivation events, gene sequence differences, gene duplications, single nucleotide polymorphisms, gene expression differences, and messenger RNA splicing variations. Evaluation of the reported findings in all these categories indicates that the CMP-sialic hydroxylase mutation is the only one that has so far been shown to result in a global biochemical and structural difference between humans and great apes. Several of the other known genetic dissimilarities deserve more exploration at the functional level. Among the areas of focus for the future should be genes affecting development, mental maturation, reproductive biology, and other aspects of life history. The approaches taken should include both going from the genome up to the adaptive potential of the organisms and going from novel adaptive regimes down to the relevant repercussions in the genome. Also, as much as we desire a simple genetic explanation for the human phenomenon, it is much more probable that our evolution occurred in multiple genetic steps, many of which must have left detectable footprints in our genomes. Ultimately, we need to know the exact number of genetic steps, the order in which they occurred, and the temporal, spatial, environmental, and cultural contexts that determined their impact on human evolution.

Fibroblast growth factors

SUMMARY

Fibroblast growth factors (FGFs) make up a large family of polypeptide growth factors that are found in organisms ranging from nematodes to humans. In vertebrates, the 22 members of the FGF family range in molecular mass from 17 to 34 kDa and share 13-71% amino acid identity. Between vertebrate species, FGFs are highly conserved in both gene structure and amino-acid sequence. FGFs have a high affinity for heparan sulfate proteoglycans and require heparan sulfate to activate one of four cell-surface FGF receptors. During embryonic development, FGFs have diverse roles in regulating cell proliferation, migration and differentiation. In the adult organism, FGFs are homeostatic factors and function in tissue repair and response to injury. When inappropriately expressed, some FGFs can contribute to the pathogenesis of cancer. A subset of the FGF family, expressed in adult tissue, is important for neuronal signal transduction in the central and peripheral nervous systems.

The mosaic structure of human pericentromeric DNA: a strategy for characterizing complex regions of the human genome

The pericentromeric regions of human chromosomes pose particular problems for both mapping and sequencing. These difficulties are due, in large part, to the presence of duplicated genomic segments that are distributed among multiple human chromosomes. To ensure contiguity of genomic sequence in these regions, we designed a sequence-based strategy to characterize different pericentromeric regions using a single (162 kb) 2p11 seed sequence as a point of reference. Molecular and cytogenetic techniques were first used to construct a paralogy map that delineated the interchromosomal distribution of duplicated segments throughout the human genome. Monochromosomal hybrid DNAs were PCR amplified by primer pairs designed to the 2p11 reference sequence. The PCR products were directly sequenced and used to develop a catalog of sequence tags for each duplicon for each chromosome. A total of 685 paralogous sequence variants were generated by sequencing 34.7 kb of paralogous pericentromeric sequence. Using PCR products as hybridization probes, we were able to identify 702 human BAC clones, of which a subset, 107 clones, were analyzed at the sequence level. We used diagnostic paralogous sequence variants to assign 65 of these BACs to at least 9 chromosomal pericentromeric regions: 1q12, 2p11, 9p11/q12, 10p11, 14q11, 15q11, 16p11, 17p11, and 22q11. Comparisons with existing sequence and physical maps for the human genome suggest that many of these BACs map to regions of the genome with sequence gaps. Our analysis indicates that large portions of pericentromeric DNA are virtually devoid of unique sequences. Instead, they consist of a mosaic of different genomic segments that have had different propensities for duplication. These biologic properties may be exploited for the rapid characterization of, not only pericentromeric DNA, but also other complex paralogous regions of the human genome.

Structure of chromosomal duplicons and their role in mediating human genomic disorders

Chromosome-specific low-copy repeats, or duplicons, occur in multiple regions of the human genome. Homologous recombination between different duplicon copies leads to chromosomal rearrangements, such as deletions, duplications, inversions, and inverted duplications, depending on the orientation of the recombining duplicons. When such rearrangements cause dosage imbalance of a developmentally important gene(s), genetic diseases now termed genomic disorders result, at a frequency of 0.7-1/1000 births. Duplicons can have simple or very complex structures, with variation in copy number from 2 to >10 repeats, and each varying in size from a few kilobases in length to hundreds of kilobases. Analysis of the different duplicons involved in human genomic disorders identifies features that may predispose to recombination, including large size and high sequence identity between the recombining copies, putative recombination promoting features, and the presence of multiple genes/pseudogenes that may include genes expressed in germ cells. Most of the chromosome rearrangements involve duplicons near pericentromeric regions, which may relate to the propensity of such regions to accumulate duplicons. Detailed analyses of the structure, polymorphic variation, and mechanisms of recombination in genomic disorders, as well as the evolutionary origin of various duplicons will further our understanding of the structure, function, and fluidity of the human genome.

CAGGG repeats and the pericentromeric duplication of the hominoid genome

Gene duplication is one of the primary forces of evolutionary change. We present data from three different pericentromeric regions of human chromosomes, which indicate that such regions of the genome have been sites of recent genomic duplication. This form of duplication has involved the evolutionary movement of segments of genomic material, including both intronic and exonic sequence, from diverse regions of the genome toward the pericentromeric regions. Sequence analyses of the target sites of duplication have identified a novel class of interspersed GC-rich repeats located precisely at the boundaries of duplication. Estimates of the evolutionary age of these duplications indicate that they have occurred between 10 and 25 mya. In contrast, comparative analyses confirm that the GC-rich pericentromeric repeats have existed within the pericentromeric regions of primate chromosomes before the divergence of the cercopithecoid and hominoid lineages ( approximately 30 mya). These data provide molecular evidence for considerable interchromosomal duplication of genic segments during the evolution of the hominoid genome and strongly implicate GC-rich repeat elements as playing a direct role in the pericentromeric localization of these events

Juxta-centromeric region of human chromosome 21 is enriched for pseudogenes and gene fragments

A physical map including four pseudogenes and 10 gene fragments and spanning 500 kb in the juxta-centromeric region of the long arm of human chromosome 21 is presented. cDNA fragments isolated from a selected cDNA library were characterized and mapped to the 831B6 YAC and to two BAC contigs that cover 250 kb of the region. An 85 kb genomic sequence located in the proximal region of the map was analyzed for putative exons. Four pseudogenes were found, including psiIGSF3, psiEIF3, psiGCT-rel whose functional copies map to chromosome 1p13, chromosome 2 and chromosome 22q11, respectively. The TTLL1 pseudogene corresponds to a new gene whose functional copy maps to chromosome 22q13. Ten gene fragments represent novel sequences that have related sequences on different human chromosomes and show 97-100% nucleotide identity to chromosome 21. These may correspond to pseudogenes on chromosome 21 and to functional genes in other chromosomes. The 85 kb genomic sequence was analyzed also for GC content, CpG islands, and repetitive sequence distribution. A GC-poor L isochore spanning 40 kb from satellite 1 was observed in the most centromeric region, next to a GC-rich H isochore that is a candidate region for the presence of functional genes. The pericentric duplication of a 7.8 kb region that is derived from the 22q13 chromosome band is described. We showed that the juxta-centromeric region of human chromosome 21 is enriched for retrotransposed pseudogenes and gene fragments transferred by interchromosome duplications, but we do not rule out the possibility that the region harbors functional genes also.

Genome duplications and other features in 12 Mb of DNA sequence from human chromosome 16p and 16q

Several publicly funded large-scale sequencing efforts have been initiated with the goal of completing the first reference human genome sequence by the year 2005. Here we present the results of analysis of 11.8 Mb of genomic sequence from chromosome 16. The apparent gene density varies throughout the region, but the number of genes predicted (84) suggests that this is a gene-poor region. This result may also suggest that the total number of human genes is likely to be at the lower end of published estimates. One of the most interesting aspects of this region of the genome is the presence of highly homologous, recently duplicated tracts of sequence distributed throughout the p-arm. Such duplications have implications for mapping and gene analysis as well as the predisposition to recurrent chromosomal structural rearrangements associated with genetic disease.

Complex beta-satellite repeat structures and the expansion of the zinc finger gene cluster in 19p12

We investigated the organization, architecture, and evolution of the largest cluster ( approximately 4 Mb) of Krüppel-associated box zinc finger (KRAB-ZNF) genes located in cytogenetic band interval 19p12. A highly integrated physical map ( approximately 700 kb) of overlapping cosmid and BAC clones was developed between genetic STS markers D19S454 and D19S269. Using ZNF91 exon-specific probes to interrogate a detailed EcoRI restriction map of the region, ZNF genes were found to be distributed in a head-to-tail fashion throughout the region with an average density of one ZNF duplicon every 150-180 kb of genomic distance. Sequence analysis of 208,967 bp of this region indicated the presence of two putative ZNF genes: one consisting of a novel member of this gene family (ZNF208) expressed ubiquitously in all tissues examined and the other representing a nonprocessed pseudogene (ZNF209), located 450 kb proximal to ZNF208. Large blocks of ( approximately 25-kb) inverted beta-satellite repeats with a remarkably symmetrical higher order repeat structure were found to bracket the functional ZNF gene. Hybridization analysis using the beta-satellite repeat as a probe indicates that beta-satellite interspersion between ZNF gene cassettes is a general property for 1.5 Mb of the ZNF gene cluster in 19p12. Both molecular clock data as well as a retroposon-mapping molecular fossil approach indicate that this ZNF cluster arose early during primate evolution (approximately 50 million years ago). We propose an evolutionary model in which heteromorphic pericentromeric repeat structures such as the beta satellites have been coopted to accommodate rapid expansion of a large gene family over a short period of evolutionary time. [The sequence data described in this paper have been submitted to GenBank under accession nos. AC003973 and AC004004.]