- Open Access
Genome-wide detection of segmental duplications and potential assembly errors in the human genome sequence
Genome Biology volume 4, Article number: R25 (2003)
Previous studies have suggested that recent segmental duplications, which are often involved in chromosome rearrangements underlying genomic disease, account for some 5% of the human genome. We have developed rapid computational heuristics based on BLAST analysis to detect segmental duplications, as well as regions containing potential sequence misassignments in the human genome assemblies.
Our analysis of the June 2002 public human genome assembly revealed that 107.4 of 3,043.1 megabases (Mb) (3.53%) of sequence contained segmental duplications, each with size equal or more than 5 kb and 90% identity. We have also detected that 38.9 Mb (1.28%) of sequence within this assembly is likely to be involved in sequence misassignment errors. Furthermore, we have identified a significant subset (199,965 of 2,327,473 or 8.6%) of single-nucleotide polymorphisms (SNPs) in the public databases that are not true SNPs but are potential paralogous sequence variants.
Using two distinct computational approaches, we have identified most of the sequences in the human genome that have undergone recent segmental duplications. Near-identical segmental duplications present a major challenge to the completion of the human genome sequence. Potential sequence misassignments detected in this study would require additional efforts to resolve.
Segments of DNA with near-identical sequence (segmental duplications or duplicons) in the human genome can be hot spots or predisposition sites for the occurrence of non-allelic homologous recombination or unequal crossing-over leading to genomic mutations such as deletion , duplication , inversion  or translocation [3, 4]. These structural alterations, in turn, can cause dosage imbalance of genetic material or lead to the generation of new gene products resulting in diseases defined as genomic disorders .
Previous studies to identify segmental duplications in the human genome have analyzed older versions of the genome assembly, which contained higher amounts of unfinished sequence and incorrectly mapped regions, and have used different computational approaches all performed by the same group [6–8]. With the human genome sequence now nearing completion, we have examined its content for segmental duplications using two distinct computational methods. In the first, we utilized the rapid BLAST2  algorithms that allow direct chromosomal-wide sequence comparisons to be made. All BLAST results reported in table formats can be subsequently grouped, parsed and analyzed for the detection of duplicated sequences. In addition, we have shown previously that there is a strong correlation between ambiguously mapped SNPs (ambSNPs), as well as the density of SNPs, and segmental duplications . AmbSNPs are SNPs that were annotated to map to two locations on a particular chromosome in the NCBI dbSNP. A subset of these ambSNPs are not true SNPs but are likely to be computer-generated nucleotide mismatches from paralogous copies of duplicated sequences and should be more appropriately labeled as paralogous sequence variants (PSVs) . Another subset is likely to be false ambSNPs of genomic sequences that have been misassigned in genome assemblies. Here, we report our analysis of all potential PSVs in the human genome and their correlation with segmental duplications as detected by our BLAST analysis. Furthermore, we provide a critical assessment on the three latest human genome assemblies from our analysis of sequence misassignments as identified from this study.
Results and discussion
Human genome segmental duplication content
On the basis of the June 2002 (NCBI Build 30) human genome assembly, a total of 107.4 Mb (3.53%) of the human genome content (3,043.1 Mb) were found to be involved in recent segmental duplications by our BLAST analysis criteria (Table 1). This content is composed of more than 1,530 distinct intrachromosomal segmental duplications (80.3 Mb or 2.64% of the total genome, Figure 1) and 1,637 distinct interchromosomal duplications (43.8 Mb or 1.44% of the total genome). In addition, 29% of all duplications are located in unfinished regions of the current genome assembly. Our results are shown using the Generic Genome Browser [11, 12]. We have also found that 38% of the duplications (52.3 Mb) can be considered as tandem duplications - defined here as two related duplicons separated by less than 200 kb.
In this study, we only analyzed large (size ≥ 5 kb) and recent (sequence identity ≥ 90%) duplications because we can achieve higher confidence and to prioritize those regions for their potential involvement in diseases. Previously, Bailey and colleagues  reported a total of 5.2% of the human genome involved in recent segmental duplications. The 1.6% discrepancy between our findings could be due to the difference in our detection criteria (size cutoff of 5 kb used in this study versus 1 kb used in Bailey et al. ). Moreover, we have identified 38.9 Mb of sequences (1.28% of the June 2002 genome assembly) likely to be artifactual duplications resulting from sequence misassignment errors present in the assembly. By comparing our results with those published previously , we found that 482/2579 clones that we identified to be involved in duplication were novel.
The molecular mechanism by which segmental duplications are created is still unclear at the moment. A recent report has suggested that Alu repeat clusters had a role as mediators of recurrent chromosomal rearrangements . We have examined whether elevated amounts of repetitive elements could be found in duplicon junctions. We inspected all duplication borders from our results and calculated the occurrence of different repeat types within the 500 bp window outside each duplicon junction. The whole-genome average frequencies were determined by sampling random 500 bp windows across the genome (excluding gap regions). Overall, we found that there are significant enrichment (or relative fold increase) for the presence of small ribonucleoprotein RNA (srpRNA), satellite, long terminal (LTR) and SINE/Alu repeats (see Additional data file). In addition, our data also showed that for some chromosomes the amount of duplicated sequence is higher in the pericentromeric and subtelomeric regions of chromosomes (Figure 1), supporting the hypothesis that these repeat-dense regions have made an important contribution to the evolution of the human genome .
Regions containing recently occurring segmental duplications can harbor rapidly evolving hominoid-specific genes, as well as novel gene families that are unique to primates [15, 16]. Using the National Center for Biotechnology Information RefSeq annotation, we identified 1,152 human genes that were mapped to duplicated regions. Of these, 475 genes were fully contained within duplicated regions and were best candidates for recent whole-gene duplication. We have carried out functional analysis of these 475 genes using the Gene Ontology Consortium database  and found that there is a significant increase in gene duplications for genes involved in immune defense (antibodies, blood-group antigens) and reproduction (pregnancy, sex differentiation) (see Additional data file).
Sequence misassignment errors in the human genome sequence assembly
We were aware that in silico detection methods, such as the ones used in this study, would not allow us to distinguish completely true duplications from artifactual duplications arising from misassigned sequences, especially in cases where sequence identity between two detected duplications exceeded 99.5% over a substantial length (> 5 kb) in regions composed of draft sequences. Although a small proportion of such results (duplications with > 99.5% identity) might represent unfinished regions of the genome that contain true duplications that have arisen very recently in the evolution of the human genome (such as the large and nearly perfect palindromic repeats located in the AZFc region on chromosome Yq11.223 involved in male infertility ), we suspect that most of the duplications (> 99.5% identity and contain draft sequences) are in fact sequence misassignment errors in the genome assembly. An explanation for such errors would be when two identical sequences belonging to the same genomic location were misassigned to distinct regions in the genome assembly.
We have used the NCBI e-PCR  to evaluate our results (potential sequence misassignment errors) from the June 2002 human genome assembly. Using some of the largest interchromosomal misassignment errors detected in our study, we found that none of the STS markers located within these misassigned sequences maps to their incorrectly assigned chromosomes. For example, AC121339 is incorrectly mapped to 3q13.13 in the June 2002 genome assembly, as supported by a consensus number of chromosome X sequence-tagged site (STS) markers (Figure 2, Table 2).
From this genome assembly, we identified that a total of 38.9 Mb of sequences, representing 1.28% of the total sequence content, are involved in such potential errors (a full list of potentially misassigned sequences can be obtained from ) that would require additional effort and further sequencing to achieve resolutions. We also analyzed an additional two previous versions of the human genome assemblies, December 2001 and April 2002, and our results showed that there has been a dramatic reduction in potential errors in the latest human genome assembly compared to the two previous genome assemblies (Table 3). Furthermore, we examined the distribution of the amount of duplications in five different categories on the basis of their level of sequence identity to each other (Table 3). We observed a large reduction in duplications that fall within the 98-100% category, supporting the fact that the genome assemblies continue to improve and correct errors made. In addition, our data showed that there have been major improvements for chromosomes 5, 6, 7, 13, 14 and 19 over the last three genome assemblies. And for chromosomes that had reached finished status, such as chromosomes 20, 21 and 22, the number of errors was negligible.
Paralogous sequence variants in the human genome
We have previously shown a strong correlation between ambSNPs with segmental duplications . AmbSNPs are SNPs that were annotated to map to two locations on a particular chromosome in the NCBI dbSNP. Here we show on a genome scale that ambSNPs most specifically correlate with intrachromosomal segmental duplications, suggesting they are paralogous sequence variants (PSVs) (Figure 1). These PSVs were perhaps mistakenly introduced into dbSNP by automated in silico-generated analysis, arising from nucleotide mismatches in paralogous copies of duplicated sequences. Overall, a surprisingly high proportion, 8.6% (199,965 of 2,327,473), of the refSNPs were annotated as ambSNPs from dbSNP (Build 108). A significant number of the ambSNPs (139,974 of 199,965 or 70.0%) are located within duplicated regions as identified by BLAST and should be regarded as PSVs (see Additional data file).
The non-identification by BLAST analysis of regions that contain ambSNPs could be due to one of three possibilities. First, the duplicated copy(s) could have been removed from the sequence assembly or the two have been conflated, that is, mistakenly thought to be the same sequence owing to their high sequence similarity. A second possibility is that the duplication is smaller than 5 kb and was excluded in our BLAST analysis. A third possibility is that a collection of ambSNPs could have been generated from misassigned sequences (identical sequences but misassigned to two difference locations in the genome) in older assembly builds due to sequencing errors or true SNPs (with high polymorphism rate) in the sequence. We also observed that the density of ambSNPs generally correlates with the size of the putative duplication, although this might be affected by the level of sequence identity between duplications. For example, two duplicated sequences sharing 98% sequence identity compared to 95% over the same length might contain fewer PSVs as the number of base-pair mismatches would be fewer in the former. In addition, we observed that regions identified by our BLAST method but that do not contain ambSNPs often correspond to artifactual duplications generated from assembly errors.
Duplicons related to genomic disorders
The size, orientation, and contents of segmental duplications are highly variable and most of them show great organizational complexity. This is perhaps due to successive transposition and rearrangement events leading to the creation of segmental duplications . In many cases, a contiguous duplicon is organized into multiple modules with different orientations and sizes. For example, one of the largest segmental duplicons detected in this study was 359 kb in size at the Williams-Beuren locus on 7q11.23 [20, 21]. In this case, the two duplicons are separated by 1.6 Mb of intervening sequence with the telomeric duplicon comprising several separate smaller modules as compared to the primary duplicon. The results presented in our study (provided in tables available at ) would also allow rapid identification of new duplicons that are potentially responsible for chromosome rearrangements and genomic disorders. For example, the location of the duplicons on chromosomes 9q34/22 q11 that have been suggested to mediate recombination leading to the Philadelphia chromosome  was identified in our analysis, as were other medically relevant chromosomal regions (Table 4) [22, 23].
The characterization of most large segmental duplications is complicated by the fact that many of them (29% of all duplications) are only represented as draft sequences from the current genome assembly. Despite the fact that both BLAST and PSVs analyses allowed us to identify most segmental duplications involved in known genomic disorder mutations (Table 4), estimations of the size of rearranged regions were different from those previously reported . In fact, with the exception of several small duplications and the segmental duplications on chromosome 22 , other regions containing duplications involved in genomic disorders were often erroneously assembled and misplaced. Furthermore, we have searched the Celera human genome C3 (publicly released version ) and C4 (subscription-based version) sequence assemblies for large duplications found on chromosome 7. We observed that most of them were not represented in large scaffolds, but instead were located in their sequence gaps, or only partially found at ends of scaffolds leading into gaps (see Table 4) . This suggests that the whole-genome assembly approach  alone might not be able to finish such duplicated regions in mammalian genomes.
We have used two different computational approaches to identify the locations of all recent segmental duplications in the current human genome draft sequence. The fidelity of the results reflects the quality of the assembly examined and the parameters used. In addition, our approach has detected numerous potential sequence misassignment errors in the current genome annotation, allowing rapid error detection in future sequence assemblies. The segmental duplication map of the human genome should serve as a guide for investigation of the role of duplications in genomic disorders, as well as their contributions to normal human genomic variability [2, 3, 27]. It is clear that genomic regions containing segmental duplications present a major challenge to the completion of the human genome sequence by April 2003. Focused efforts including targeted sequencing of allele-specific clones, high-resolution fluorescence in situ hybridization, and expert curation would be required to validate the actual (or proposed) organization of these complex regions as well as to complete the human genome reference sequence.
Materials and methods
Genome sequence and chromosome-wide BLAST
We obtained the December 2001, April 2002, and June 2002 (NCBI Build 28, 29 and 30 respectively) human genome assemblies through the University of California, Santa Cruz Human Genome Browser . All chromosome sequences were lower-case masked for highly repetitive elements by RepeatMasker (A.F.A. Smit and P. Green, unpublished). For each assembly build, each of the 25 masked chromosome sequences (including one unmapped chromosome sequence 'ChrUn') was compared against itself by chromosome-wide BLAST2  to detect intrachromosomal segmental duplications (25 comparisons made), as well as pairwise comparisons to each of the other 24 chromosomes to detect interchromosomal segmental duplications (600 comparisons made). All BLAST results were subsequently parsed to eliminate low-quality and fragmented alignments under the following criteria: BLAST results having ≥ 90% sequence identity, ≥ 80 bp in length, and with expected value ≤ 10-30.
BLAST results parsing and duplication detection
Each BLAST report was sorted by chromosomal coordinates. All identical hits (same coordinate alignments), including suboptimal BLAST alignments recognized by multiple, overlapping alignments, as well as mirror hits (reverse coordinate alignments) from the BLAST results of the intrachromosomal set were removed. Contiguous alignments separated by a distance of less than 3 kb, then 5 kb, and subsequently 9 kb were joined (stepwise) into modules in order to traverse masked repetitive sequences and to overcome breaks in the BLAST alignments caused by insertions/deletions and sequence gaps. Such contiguous sequence alignment modules represent sequence similarity between the subject and query chromosome sequence in question (at their respective positional coordinates). This pairwise sequence comparison procedure serves as a rapid and robust way to detect duplication relationships. However, because of the use of masked sequences, our method would only yield a poor (on average 0.1-0.5 kb) resolution for the determination of the precise duplication alignment boundaries. Results were classified as either duplications or 'questionable' results based on sequencing status of the region and the percent sequence similarity between the detected alignments. Questionable duplications are results that fall within regions containing draft sequences with > 99.5 % detected sequence identity with another region. We consider these questionable duplications to be involved in potential sequence misassignment errors in the human genome assembly and would require further effort to achieve resolution.
Fine mapping of segmental duplications
Detailed information regarding segmental duplications as well as potential sequence misassignment errors identified by our analysis were presented using the Generic Genome Browser [11, 12]. We have also summarized our results in table formats  that include information on size of duplications, chromosomal band locations, level of identity between duplicated copies, sequenced clones (accession numbers) and their sequencing status, as well as genes mapped to these regions. In addition, we have plotted the size of each intrachromosomal duplication (y-axis) against its chromosome position (x-axis) along each chromosome to indicate the intrachromosomal segmental duplication content of each chromosome (Figure 1) using the publicly available visualization tool GenomePixelizer [29, 30]. Results generated from the detection of segmental duplications were subsequently converted into coordinate files as input for display using GenomePixelizer.
Paralogous sequence variants (PSV) density map
SNP mapping data from dbSNP were obtained through the NCBI ftp site . Each chromosome SNP table, containing annotation regarding ambSNPs that have appeared twice in a particular chromosome, were extracted and sorted along with their corresponding chromosomal positions. The number of ambSNPs was tabulated along a 10-kb window to produce density plots of ambSNPs along the length of each chromosome (Figure 1).
Additional data file
The following additional data are availiable as a single file with this article: a table describing the relative frequency increase of repeat types at duplicon junctions (additional data file 1); a table giving the number of ambiguousgly mapped SNPs within segmental duplications (number of PSVs) (additional data file 1); a figure (additional data file 1) showing functional profiling of genes involved in recent whole-gene duplication vs human genome average.
Reiter LT, Murakami T, Koeuth T, Pentao L, Muzny DM, Gibbs RA, Lupski JR: A recombination hotspot responsible for two inherited peripheral neuropathies is located near a mariner transposon-like element. Nat Genet. 1996, 12: 288-297.
Giglio S, Broman KW, Matsumoto N, Calvari V, Gimelli G, Neumann T, Ohashi H, Voullaire L, Larizza D, Giorda R, et al: Olfactory receptor-gene clusters, genomic-inversion polymorphisms, and common chromosome rearrangements. Am J Hum Genet. 2001, 68: 874-883. 10.1086/319506.
Giglio S, Calvari V, Gregato G, Gimelli G, Camanini S, Giorda R, Ragusa A, Guerneri S, Selicorni A, Stumm M, et al: Heterozygous submicroscopic inversions involving olfactory receptor-gene clusters mediate the recurrent t(4;8) (p16; p23) translocation. Am J Hum Genet. 2002, 71: 276-285. 10.1086/341610.
Saglio G, Storlazzi CT, Giugliano E, Surace C, Anelli L, Rege-Cambrin G, Zagaria A, Jimenez Velasco A, Heiniger A, Scaravaglio P, et al: A 76-kb duplicon maps close to the BCR gene on chromosome 22 and the ABL gene on chromosome 9: possible involvement in the genesis of the Philadelphia chromosome translocation. Proc Natl Acad Sci USA. 2002, 99: 9882-9887. 10.1073/pnas.152171299.
Lupski JR: Genomic disorders: structural features of the genome can lead to DNA rearrangements and human disease traits. Trends Genet. 1998, 14: 417-422. 10.1016/S0168-9525(98)01555-8.
Bailey JA, Yavor AM, Massa HF, Trask BJ, Eichler EE: Segmental duplications: organization and impact within the current human genome project assembly. Genome Res. 2001, 11: 1005-1017. 10.1101/gr.GR-1871R.
International Human Genome Sequencing Consortium: Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921. 10.1038/35057062.
Bailey JA, Gu Z, Clark RA, Reinert K, Samonte RV, Schwartz S, Adams MD, Myers EW, Li PW, Eichler EE: Recent segmental duplications in the human genome. Science. 2002, 297: 1003-1007. 10.1126/science.1072047.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410. 10.1006/jmbi.1990.9999.
Estivill X, Cheung J, Pujana MA, Nakabayashi K, Scherer SW, Tsui LC: Chromosomal regions containing high-density and ambiguous-mapped single nucleotide polymorphisms (SNPs) correlate with segmental duplications in the human genome. Hum Mol Genet. 2002, 11: 1987-1995. 10.1093/hmg/11.17.1987.
Stein LD, Mungall C, Shu S, Caudy M, Mangone M, Day A, Nickerson E, Stajich JE, Harris TW, Arva A, et al: The generic genome browser: a building block for a model organism system database. Genome Res. 2002, 12: 1599-1610. 10.1101/gr.403602.
TCAG - Human Segmental Duplication Homepage. [http://chr7.ocgc.ca/humandup]
Kolomietz E, Meyn MS, Pandita A, Squire JA: The role of Alu repeat clusters as mediators of recurrent chromosomal aberrations in tumors. Genes Chromosomes Cancer. 2002, 35: 97-112. 10.1002/gcc.10111.
Samonte RV, Eichler EE: Segmental duplications and the evolution of the primate genome. Nat Rev Genet. 2002, 3: 65-72. 10.1038/nrg705.
Lynch M, Conery JS: The evolutionary fate and consequences of duplicate genes. Science. 2000, 290: 1151-1155. 10.1126/science.290.5494.1151.
Johnson ME, Viggiano L, Bailey JA, Abdul-Rauf M, Goodwin G, Rocchi M, Eichler EE: Positive selection of a gene family during the emergence of humans and African apes. Nature. 2001, 413: 514-519. 10.1038/35097067.
Gene Ontology Consortium. [http://www.geneontology.org]
Kuroda-Kawaguchi T, Skaletsky H, Brown LG, Minx PJ, Cordum HS, Waterston RH, Wilson RK, Silber S, Oates R, Rozen S, et al: The AZFc region of the Y chromosome features massive palindromes and uniform recurrent deletions in infertile men. Nat Genet. 2001, 29: 279-286. 10.1038/ng757.
Electronic PCR. [http://www.ncbi.nlm.nih.gov/genome/sts/epcr.cgi]
Valero MC, de Luis O, Cruces J, Perez Jurado LA: Fine-scale comparative mapping of the human 7q11.23 region and the orthologous region on mouse chromosome 5G: the low-copy repeats that flank the Williams-Beuren syndrome deletion arose at breakpoint sites of an evolutionary inversion(s). Genomics. 2000, 69: 1-13. 10.1006/geno.2000.6312.
Osborne LR, Li M, Pober B, Chitayat D, Bodurtha J, Mandel A, Costa T, Grebe T, Cox S, Tsui LC, et al: A 1.5 million-base pair inversion polymorphism in families with Williams-Beuren syndrome. Nat Genet. 2001, 29: 321-325. 10.1038/ng753.
Emanuel BS, Shaikh TH: Segmental duplications: an 'expanding' role in genomic instability and disease. Nat Rev Genet. 2001, 2: 791-800. 10.1038/35093500.
Stankiewicz P, Lupski JR: Genome architecture, rearrangements and genomic disorders. Trends Genet. 2002, 18: 74-82. 10.1016/S0168-9525(02)02592-1.
Shaikh TH, Kurahashi H, Saitta SC, O'Hare AM, Hu P, Roe BA, Driscoll DA, McDonald-McGinn DM, Zackai EH, Budarf ML, et al: Chromosome 22-specific low copy repeats and the 22q11.2 deletion syndrome: genomic organization and deletion endpoint analysis. Hum Mol Genet. 2000, 9: 489-501. 10.1093/hmg/9.4.489.
Venter JC, Adams MD, Myers EW, Li PW, Mural RJ, Sutton GG, Smith HO, Yandell M, Evans CA, Holt RA, et al: The sequence of the human genome. Science. 2001, 291: 1304-1351. 10.1126/science.1058040.
Scherer SW, Cheung J: Discovery of the human genome sequence in the public and private databases. Curr Biol. 2001, 11: R808-R811. 10.1016/S0960-9822(01)00490-0.
Gratacos M, Nadal M, Martin-Santos R, Pujana MA, Gago J, Peral B, Armengol L, Ponsa I, Miro R, Bulbena A, et al: A polymorphic genomic duplication on human chromosome 15 is a susceptibility factor for panic and phobic disorders. Cell. 2001, 106: 367-379.
University of California, Santa Cruz, Genome Informatics. [http://genome.ucsc.edu]
Kozik A, Kochetkova E, Michelmore R: GenomePixelizer - a visualization program for comparative genomics within and between species. Bioinformatics. 2002, 18: 335-336. 10.1093/bioinformatics/18.2.335.
GenomePixelizer: genome visualization tool. [http://www.atgc.org/GenomePixelizer/GenomePixelizer_Welcome.html]
NCBI ftp site. [ftp://ftp.ncbi.nih.gov/snp/human/chr_rpts]
We thank L. Zaifman, J. Zhang and G. Duggan for expert programming assistance. This work was supported by the Canadian Institutes of Health Research (CIHR) and Genome Canada to S.W.S. and L.-C.T. X.E. is a Senior Scientist at the Centre de Regulació Genòmica (CRG) and a Visiting Scientist of the Hospital for Sick Children Research Institute. L.-C.T is a Distinguished Scientist of CIHR and Sellers Chair of Cystic Fibrosis Research. S.W.S. is a Scholar of CIHR and International Scholar of the Howard Hughes Medical Institute.