This Article Abstract Full Text (PDF) Supplemental Research Data Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Warburton, P. E. Articles by Benson, G. Search for Related Content PubMed PubMed Citation Articles by Warburton, P. E. Articles by Benson, G. Social Bookmarking                   What's this?

Letter

Inverted Repeat Structure of the Human Genome: The X-Chromosome Contains a Preponderance of Large, Highly Homologous Inverted Repeats That Contain Testes Genes

¹ Department of Human Genetics, Mount Sinai School of Medicine, New York, New York 10029, USA ² Department of Computer Science, Department of Biology, Boston University, Boston, Massachusetts 02215, USA ³ Laboratory for Biocomputing and Informatics, Boston University, Boston, Massachusetts 02215, USA

	ABSTRACT

Top ABSTRACT RESULTS DISCUSSION METHODS REFERENCES WEB SITE REFERENCES

 
We have performed the first genome-wide analysis of the Inverted Repeat (IR) structure in the human genome, using a novel and efficient software package called Inverted Repeats Finder (IRF). After masking of known repetitive elements, IRF detected 22,624 human IRs characterized by arm size from 25 bp to >100 kb with at least 75% identity, and spacer length up to 100 kb. This analysis required 6 h on a desktop PC. In all, 166 IRs had arm lengths >8 kb. From this set, IRs were excluded if they were in unfinished/unassembled regions of the genome, or clustered with other closely related IRs, yielding a set of 96 large IRs. Of these, 24 (25%) occurred on the X-chromosome, although it represents only 5% of the genome. Of the X-chromosome IRs, 83.3% were 99% identical, compared with 28.8% of autosomal IRs. Eleven IRs from Chromosome X, one from Chromosome 11, and seven already described from Chromosome Y contain genes predominantly expressed in testis. PCR analysis of eight of these IRs correctly amplified the corresponding region in the human genome, and six were also confirmed in gorilla or chimpanzee genomes. Similarity dot-plots revealed that 22 IRs contained further secondary homologous structures partially categorized into three distinct patterns. The prevalence of large highly homologous IRs containing testes genes on the X- and Y-chromosomes suggests a possible role in male germ-line gene expression and/or maintaining sequence integrity by gene conversion.

Inverted Repeats (IRs) make up one class of human duplications, which consist of two arms of similar DNA—with one inverted and complemented relative to the other—around a central, usually nonhomologous spacer. Large IRs have been observed in physical maps of the X-chromosome, and have been associated with chromosomal rearrangements and gene deletions (Lafreniere et al. 1993; Small et al. 1997; McDonell et al. 2000; Aradhya et al. 2001). Recently, the finished sequence of the human Y-chromosome has revealed the presence of several remarkably large and highly homologous IRs, up to 1.4 Mb in size and 99.97% identity, which contain Y-specific genes expressed in testes and thought to be required for spermatogenesis (Skaletsky et al. 2003). Gene conversion is evidently maintaining the homology between the arms of these palindromes, and thus the sequence integrity and function of the genes in the absence of meiotic recombination between homologs (Rozen et al. 2003).

IRs are widespread in both prokaryotic and eukaryotic genomes, and have been associated with a myriad of possible functions, reviewed in Pearson et al. (1996). Some IRs are capable of extruding into DNA cruciforms, structures in which the normal double-stranded DNA denatures, and complementary arms in the same strand form intrastrand double helices, or stems. The spacer regions become unpaired loops at the top of each stem, and the four-way junction where the bases of the stems meet is indistinguishable from a Holliday structure. The ability of particular IRs to extrude into cruciforms depends on the size and sequence composition of both the arms and spacer region. The energy barrier required to denature the DNA and extrude into cruciforms is reduced by the unwinding torsional stress induced by local negative supercoiling (Shlyakhtenko et al. 2000; Benham et al. 2002).

To facilitate the study of genomic IRs, a novel and efficient computer program called Inverted Repeats Finder (IRF) was developed, and the first genome-wide analysis of IRs in the human genome was performed. The largest and most homologous IRs found in the human genome, described in detail in this report, showed a disproportionately high representation on the X-chromosome. The program IRF, as well as updated descriptions of the IR structure of future assemblies of the human genome, will be made publicly available at the Inverted Repeat Data Base (IRDB) (http://tandem.bu.edu).

	RESULTS

Top ABSTRACT RESULTS DISCUSSION METHODS REFERENCES WEB SITE REFERENCES

 
Inverted Repeats Finder (IRF) Reveals a Preponderance of Large, Highly Homologous IRs on the X-Chromosome
IRF was run against each human chromosome from the latest available version of the human genome sequence (hg16). Repeat-Masker (Jurka 2000; A.F.A. Smit and P. Green, unpubl.; http://ftp.genome.washington.edu/RM/RepeatMasker.html) was used to exclude repetitive elements during identification of candidate IRs, thereby blocking detection of hundreds of thousands of biologically uninteresting IRs that consist of two nearby homologous, oppositely oriented interspersed repeats. However, during alignment and extension of IRs, repetitive elements were included. Short and low-identity IRs initially detected by IRF were filtered out (see Methods), leaving a set of 22,624 IRs with spacer lengths up to 100 kb that were 75% identical between arms (Fig. 1A,B). These IRs had arm lengths from 25 bp up to 500 kb, albeit skewed toward lengths under 100 bp (Fig. 1A). The spacer lengths were skewed toward lengths <500 bp (mean of 9358 bp, median of 327 bp; data not shown).

View larger version (34K): [in this window] [in a new window]   Figure 1 Results from Inverted Repeat Finder. (A) Distribution of arm lengths of human IRs detected by IRF, mean of 551 bp, median of 85 bp. Note the log scale on the x-axis. (B) Distribution of arm length (x-axis) by percent identity between arms (y-axis). The x-axes in A and B correspond for direct comparison. The "birdcrest" pattern observed resulted from limited ranges of percent identity values for the shortest most similar repeats, for example, 1/25 bp, 1/26 bp,..., 2/25 bp, 2/26 bp,.... As IR length increases, the pattern becomes less constrained. In all, 166 IRs had arm lengths 8000 bp (vertical line), of which 37 were on the X-chromosome and 18 on the Y-chromosome. (C) IRs detected by IRF that are 8000 bp, 95% identical, of which 27 are from the X-chromosome and 10 from the Y-chromosome. (D) Comparison on each human chromosome of percent total IRs (22,624), percent IRs 8 kb (166) detected by IRF, and percent IRs 8 kb after exclusion (96; see Supplemental S1). Values are compared to the percent of total assembled genome (3.07 x 109 bp) for each chromosome.

View this table: [in this window] [in a new window]   Table 1. Inverted Repeats in the Human Genome, >8000 bp, >99% Homology

To produce the most robust list of large IRs possible, the 166 largest IRs detected by IRF were evaluated by visualizing them on the UCSC Genome Browser using the custom track file provided by IRF (see Methods). Note that when viewing these custom tracks, a striking mirror symmetrical pattern for the arms of each IR is seen in the RepeatMasker Tracks of the UCSC Genome Browser. We took a conservative approach and removed all IRs that were possible false positives due to assembly errors, such as IRs that span gaps (12) or abut gaps (9) (Supplemental data S1). However, several of the IRs excluded as potential false positives may be confirmed as true IRs upon further refinement of the sequence assembly. Conversely, it is possible that the arms of highly homologous IRs may have been inadvertently combined into a single sequence, resulting in false negatives undetected by IRF. It will be important to compare the results of IRF from different builds of the human genome available at IRDB (http://tandem.bu.edu) to re-evaluate IRs that were excluded as potential false positives in the current set.

Overall, 70 of the 166 large IRs initially detected by IRF in hg16 build 34 were excluded based on several criteria (see Methods; Supplemental data S1), yielding a final set of 96 IRs (Supplemental data S2). This exclusion process did not significantly alter the distribution of IRs across the chromosomes and specifically did not affect the significance of the high proportion of IRs found on the X-chromosome (24 IRs, 25%) and Y-chromosome (13 IRs, 13.6%; Fig. 1D). The high proportion of X-chromosome IRs did not appear to be caused by any bias in the quality of the DNA sequence, as the X-chromosome has a relatively high proportion of gaps (Kent et al. 2002). Table 1 lists the IRs that have 99% arm identity. Remarkably, 20 of the 24 (83.3%) X-chromosome IRs were 99% identical, whereas only 17 of the 59 (28.8%) autosomal IRs were 99% identical (Table 2).

View this table: [in this window] [in a new window]   Table 2. Total Human IRs >8000 bp on X- and Y-chromosomes and Autosomes

Large IRs on the X-Chromosome Predominantly Contain Testes Genes
Analysis of the genes contained within the set of 96 large IRs revealed a striking preference for genes that are predominantly or exclusively expressed in testes (Table 1), organized in opposite orientation on either arm of the IR, similar to the Y-chromosome (Skaletsky et al. 2003). Eleven large IRs from Chromosome X and one from Chromosome 11 contain a gene for which either an RT-PCR or a Northern blot directly demonstrating testes expression has been published (Table 1). Eight of these genes are expressed exclusively in the testes (from normal tissues), and four were expressed in at least one other tissue (Table 1). Eight of the Chromosome X genes were identified as Cancer-testes antigen (CTA) genes (Table 1), which are expressed predominantly or exclusively from testes, and in certain cancers (Zendman et al. 2003a; Scanlan et al. 2004). These results demonstrate that the human X-chromosome contains a preponderance of large, highly homologous IRs that contain testes genes (Table 2). Approximately 20% of the 52 genes on the human X-chromosome that are expressed predominantly in testes (from the GNF atlas 2 database; see Methods) are represented in these IRs. Many of the remaining IRs listed in Table 1 contained mRNAs cloned from various cDNA libraries. A representative mRNA and tissue was listed for each, with testes mRNAs preferentially identified when present, although these were not counted as known testes genes in Tables 1 and 2.

Analysis of IRs in the Mouse Genome
IRF will be useful for the analysis of any sequenced genome for which RepeatMasking is possible. For comparison to a non-primate mammal, we ran IRF on the current version of the mouse genome (NCBI build 32; Waterston et al. 2002). The large proportion of unfinished sequence and many small and large gaps in the current mouse assembly reduce confidence in the reliability of the set of IRs detected by IRF. For example, IRF detected 303 large IRs (spacer 100 kb, arm 8 kb, 75% arm identity) of which 196 (64.6%) span gaps, as compared with 12/166 (7.2%) of large IRs that span gaps in the human analysis. Furthermore, many of the internal IR arm/spacer boundaries were found exactly at the end of an assembled BAC sequence, which further suggests the possibility of an assembly error. Nevertheless, of the 107 large IRs detected in the mouse that do not span gaps, 46 (42.9%) are found on the X-chromosome. Of the nine large IRs on the mouse X-chromosome that contain genes, six contain testes genes, including mouse SSX genes (SSX4 and 5 and SSXB7), testes-expressed homeobox genes (TgiFx1 and TEX2), a testes-specific ferritin-like gene (FH17), and BC061169 (Xmr), which contains a Cor1 domain, a component of the chromosome core in the meiotic prophase chromosomes. Thus, this preliminary analysis of the IR structure of the mouse genome strongly suggests that the mouse X-chromosome also contains a preponderance of large IRs that contain testes genes. However, the analysis must be revisited as the mouse genome sequence and assembly are improved. IRF analyses of the current and future assemblies of mouse genome will be included in the IRDB.

Analysis of IR Structure in Human Chromosome Xp11.22
Similarity dot-plot analysis was performed on each large IR and surrounding genomic DNA, to reveal details about repetitive DNA organization and potential for secondary structure formation (Kuroda-Kawaguchi et al. 2001). Figure 2A shows this analysis for the 2-Mb region that contained the highest density of large IRs in the genome, located on Chromosome Xp11.22. IRX51.17 directly abuts the distal edge of a 100-kb gap that remains in Xp11.22 in hg16, and this IR was not included in the final set of 96 large IRs (Supplemental data S2). However, it is the longest and most homologous IR on the X-chromosome, clearly seen as a long vertical line (Fig. 2B). It is also highly homologous to IRX-51.4958, both of which contain GAGE D2 genes (Fig. 2B), and one of these IRs may be eliminated when the intervening gap is closed. Notably, the Y-chromosome palindromes P1 and P2 containing the DAZ genes are also highly homologous, and appear to share a similar structure (although there is no intervening gap; Kuroda-Kawaguchi et al. 2001; Skaletsky et al. 2003).

View larger version (35K): [in this window] [in a new window]   Figure 2 IRs in Xp11.2. (A) Similarity dot-plot of 2 Mb from Xp11.2. Homologous regions indicated by horizontal lines (direct repeats) or vertical lines (inverted repeats). Six large IRs were detected in this region as vertical lines (gray triangles). No other significant sequence similarities were seen in the top portion of the dot-plot not shown. The window size and percent identity for the dot-plot are indicated. (B) Higher-resolution view of the region containing IRX-51.17 and IRX-51.4958, which both contain the GAGE-D genes. Arrows on the same line represent homologous regions as indicated by the dot-plot; percent identity indicated when calculated. Internal IRs detected by IRF are indicated, for example, IRX-51.489 (Supplemental data S1), and the midpoint of the spacer indicated by a dot. IRX-51.17 was not included in the final set of 96 large IRs. (C) Higher-resolution view of IRX-51.73 and surrounding region, which contains an SSX gene cluster. Both arms of the IRX-51.73 contain inverted copies of SSX2 and SSX pseudo5. As in B, internal IRs detected by IRF are indicated, for example, IRX-51.635, and the midpoint of the spacer of the IR for which details are provided is indicated by a dot (Supplemental data S1).

Classification of Large IRs Into Three Distinct Patterns of Genomic Organization
Similarity dot-plot analysis of the IRs in Table 1 showed that 22 displayed a complex pattern of inverted and/or tandem regions of similarity, whereas the remainder simply showed a single vertical line. Of these IRs, 16 could be classified into three distinct patterns based on the genomic organization of inverted and tandem repeats within the large IRs, which each suggest different potential secondary structures (Fig. 3). The first pattern of genomic organization consists of IRs that contain mirror image IRs within both arms of the larger IR (Fig. 3A), seen in three large IRs (Table 1). This type of pattern was also present in the largest Y-chromosome palindrome P1, which contains within each arm the 9.9-kb P1.1 and P1.2 (Kuroda-Kawaguchi et al. 2001; Skaletsky et al. 2003). Notably, these types of IRs would have the potential to extrude into unusual double cruciform structures (Fig. 3D).

View larger version (42K): [in this window] [in a new window]   Figure 3 Secondary structure patterns of IRs. (A) Dot-plot similarity analysis of IRX-69.83, which contains smaller, less homologous IRs in each arm. Internal IRs detected by IRF and percent identity are indicated (Supplemental data S1 and S2). Below each dot-plot are the RepeatMasker and TRF (Benson 1999) tracks from the UCSC Genome Browser, where the IR and internal secondary structures can be observed visually as mirror symmetries. (B) Dot-plot similarity analysis of IRX-150.52. The spacer of this IR is homologous to a region in the arms. (C) Analysis of IR7-143.36 and IR7-143.41. One arm of both of these IRs is found in the spacer of the other. IRs with similar secondary structure to A, B, and C are indicated in Table 1. Internal IRs detected in these regions were eliminated from the final data set (Supplemental data S1 and S2). The dot-plot similarity analysis was performed with a window size of 50 bp and mimimum identity of 85%. (D) Possible double cruciform structure suggested by IRs that contain internal IRs, for example, in A and Figure 2C. (Black and gray DNA strands) outside of IR region; (blue and red DNA strands) large, high homology IR; (purple and orange DNA strands) internal, lower homology IRs.

Conservation of IR Arm Boundaries in Great Apes
To confirm these IRs in the human genome, PCR primers were designed from the human genome sequence (hg16) to define arm-specific STSs for several examples. One primer (primer A, Fig. 4) was designed to hybridize to sequences present in both arms, on either side of the spacer. Individual primers were designed within the spacer to be specific for either the left or right arm/spacer boundary (primers L and R, respectively). In each case (Fig. 4), PCR amplification of human genomic DNA using primer pairs A-L and A-R resulted in PCR amplification of the predicted size fragments (300-900 bp), and thus confirmed the presence of both arms of the IR in the human genome. Sequence analysis of these PCR products confirmed their correspondence to the appropriate arm-spacer boundaries in the human genome sequence (Supplemental data S3).

View larger version (20K): [in this window] [in a new window]   Figure 4 PCR analysis of internal arm/spacer boundary. For each IR indicated, a single primer (A) was designed to hybridize to both arms. Individual primers (L and R) were designed to hybridize to the spacer region. PCR amplification was performed using primer pairs A-L and A-R with human, gorilla, and chimpanzee genomic DNA. For each IR, a + indicates that PCR products were amplified using both primer pairs. All amplified IRs are >99% homologous (Table 1), except IRX-51.924, which is 97% homologous (Fig. 2A; Supplemental data S1 and S2). The PCR primers and DNA sequence of PCR products are shown in Supplemental data S3. The high homology between the arm and spacer of IRX-51.17 and IRX-51.924 (see Fig. 2B) does not allow for specific STSs, although two distinct sets of primers pairs were used (Supplemental data S3).

	DISCUSSION

Top ABSTRACT RESULTS DISCUSSION METHODS REFERENCES WEB SITE REFERENCES

 
We have performed the first genome-wide analysis of the IR structure of the current version (build 34) of the human genome DNA sequence, using a novel efficient software package called IRF, which is available for ongoing future analyses at IRDB (http://tandem.bu.edu). We found that the human X-chromosome contained a disproportionately high number of large, highly homologous IRs that contained testes genes. This is highly analogous to the IRs found on the human Y-chromosome, which are evidently undergoing conversion to preserve gene integrity and the function of male fertility genes in the absence of meiotic pairing and crossing over (Rozen et al. 2003; Skaletsky et al. 2003). Arm-to-arm gene conversion would not be as critical for maintaining gene integrity on autosomes and the X-chromosome, because they do undergo meiotic crossing over. However, the X-chromosome does so only in females, and thus at a 50% reduced frequency relative to autosomes. Although our arm-specific STSs did not provide enough DNA sequence to assess gene conversion between arms in this report, gene conversion has been previously described between the arms of at least two X-chromosome IRs in Xq28, IRX-152.06 and IRX-152.30 (Small et al. 1997; Aradhya et al. 2001).

The highly homologous IRs on the X-chromosome predominantly contain genes expressed in testes, suggesting a possible role in male germ-line gene expression. The accumulation of sex-linked genes on the X-chromosome appears to be dependent on their timing of expression in meiosis. The mammalian X- and Y-chromosomes undergo male germ-line sex chromosome inactivation (MSCI), which prevents expression of X- and Y-linked genes during meiotic pachytene. The X-chromosome appears to accumulate spermatogenesis genes that are expressed prior to MSCI (Wang et al. 2001; Wu and Xu 2003; Khil et al. 2004). These genes are consistent with a model of sexually antagonistic alleles in which recessive genes beneficial to XY males but detrimental to XX females would accumulate on the X chromosomes because the detrimental effects would initially be masked in females due to heterozygosity (Wang et al. 2001; Lercher et al. 2003). However, as these genes accumulate on the X-chromosomes in the population, modifiers would be expected to arise that limit the genes' expression to the male. The formation of IR-based structures could provide one mechanism to limit expression to the male germ line, preventing deleterious expression of these X-linked alleles in the female (Wu and Xu 2003). The IRs on the Y-chromosome may play a similar role in male germ-line expression.

Most genes expressed during later stages of spermatogenesis have been found on autosomes (Eddy and O'Brien 1998; Emerson et al. 2004; Khil et al. 2004). Furthermore, MSCI of many essential X-linked genes appears to be compensated by the testes-specific expression of autosomal retrotransposed copies (Wu and Xu 2003; Emerson et al. 2004; Wang 2004). However, none of the X-linked testes genes found in IRs (Table 1) have retrotransposed autosomal counterparts (Emerson et al. 2004; Wang 2004), although expression of at least some of them, for example, MAGE-A, GAGE-D, and SPANX, may be required during or after MSCI (Zendman et al. 2003a). Thus, the IRs on which these genes are found may, through formation of cruciforms or other unusual chromatin structures, permit escape from meiotic X inactivation and permit expression of critical spermatogenesis genes that remain exclusively on the X-chromosome (Skaletsky et al. 2003).

The largest autosomal IR observed (IR11-89.34), which was present in great apes by our PCR assay (Fig. 4), also contains the testes gene RNF18 (Table 1). IR1-147.04, also present in great apes (Fig. 4), contains the only human histone gene cluster organized in an IR, which, however, does not contain histone H2B, for which a testes-specific homolog has been identified (Zalensky et al. 2002). Thus, if autosomal IRs facilitate male germ-cell expression, this IR may represent a subset of replication-dependent histone genes used during germ cell development or meiosis (Marzluff et al. 2002).

The highly homologous IRs described here suggest the formation of large DNA cruciform structures, the arms of which would be indistinguishable from normal double-stranded DNA. Such large cruciforms could both replicate and be transcribed essentially normally. They would be exquisite structures for regulating the topological state of chromosomal regions, especially during chromatin remodeling and/or nucleosome replacement. Removal of nucleosomes from DNA creates negative superhelical twist, which could be relaxed by extrusion into a cruciform. Some IRs with complex secondary structures (Fig. 3; Table 1) could form distinct cruciform structures, such as the double cruciform shown in Figure 3D. In IRs in which the spacer is homologous to the arm (Fig. 3B; Table 1), multiple cruciforms are possible, each pairing different sets of homologous genes and permitting conversion between them. And when one IR is surrounded by another IR (Fig. 3C; Table 1), two mutually exclusive cruciforms are possible.

The large and highly homologous IRs described here could potentially lead to aberrant sister-chromatid exchange and chromosome rearrangements. Notably, several human isodicentric Xq chromosome breakpoints have been mapped to Xp11.22 (Wolff et al. 1996), the region of densest large IR occurrence (Fig. 2), and at least three isodicentric Xp chromosome breakpoints have been precisely mapped to the 119-kb IRX-70.95 (Table 1; McDonell et al. 2000). Furthermore, rearrangements between the arms of at least two IRs in Xq28, IRX152.06 and IRX-152.30, have been implicated in deletions of the emerin gene and the NEMO gene, respectively (Small et al. 1997; Aradhya et al. 2001). Thus, the IRs described in this report may represent regions that are particularly prone to genomic rearrangements.

To summarize, we have examined the IR structure of the human genome, and revealed a remarkable preponderance of large, highly homologous IRs on the human X-chromosome, in regions containing testes genes. These IRs may play an important evolutionary or regulatory role in controlling sex-specific gene expression critical during germ-cell development or meiosis.

	METHODS

Top ABSTRACT RESULTS DISCUSSION METHODS REFERENCES WEB SITE REFERENCES

 
IRF (http://tandem.bu.edu/cgi-bin/irdb/irdb.exe) is a prototype tool for identifying approximate inverted repeats in nucleotide sequences that is similar in concept to the Tandem Repeats Finder (Benson 1999). Candidate IRs are detected by finding short, exact, reverse-complement matches of 4-7 nt (k-tuples) between nonoverlapping fragments of a sequence. A "center" position is defined for each k-tuple match. Short k-tuples are used to detect short IRs with short spacers, and longer k-tuples are used to detect longer IRs with potentially larger spacers, typically 10-100 kb. The program detects "clusters" of k-tuple matches having the same or nearly the same center and falling within a small interval of sequence. Several interval sizes are prespecified, typically between 30 and 2000 nt long.

Candidate IRs are confirmed (aligned and extended) or rejected by computing Smith-Waterman style similarity alignment. An efficient "narrowband" technique is used (Benson 1999), which computes alignment scores in a band of specified width around the presumed correct alignment, shifting the band as the location of the alignment shifts. When an alignment exceeds a prespecified minimum alignment score, the IR pair is reported. The alignment will terminate at the border of an insertion/deletion in one of the repeat copies when it is longer than the bandwidth. The remaining matching parts may be detected as an independent IR pair with a significantly shifted center. In the IRF version used here, the maximum alignment was 500,000 bp, with a bandwidth of 200 bp.

IRF was run against human genome sequences of each chromosome (hg16, obtained from NCBI), using parameters 2,3,5,40::match, mismatch, indel, minimum score. Chromosomes were run simultaneously on a computing cluster consisting of four nodes, each with a single 2.4 GHz Pentium 4 processor, 2 GB RAM, and 80 GB hard disk storage. Analysis took 1.5 h, the equivalent of 6 h on a single computer. The data were then inserted into a Microsoft SQL Server database.

At these settings of IRF, the shortest IR detected has 20 nt in each arm with a identity of 100% between the arms. For all runs, repetitive elements masked by RepeatMasker (Jurka 2000; A.F.A. Smit and P. Green, unpubl.; http://ftp.genome.washington.edu/RM/RepeatMasker.html) were excluded during initial candidate detection but were included in subsequent alignment and extension of IRs. The initial set of IRs detected by IRF was filtered to remove redundant IRs, defined as those that shared a positional identity 60%, removing the IRs with the lower "homology times length" score. These IRs were then filtered to retain those with arm identity 75% and a homology-times-length score of 25 bp, yielding the final data set of 22,624 IRs (spacer 100 kb; Fig. 1A,B) or 31,163 IRs (spacer 500 kb). Analysis of these IRs was performed using the data download feature of IRF and Microsoft Excel.

The complete set of IRs detected by IRF or any subsequently filtered subset can be displayed on the assembled human genome sequence using the UCSC Genome Browser (Kent et al. 2002) by downloading a custom track gff file generated by IRF onto your computer and uploading this file into the UCSC custom track option. IRs were confirmed by visual analysis of the mirror symmetrical pattern of the RepeatMasker track, and by BLAT searches with DNA sequences from one arm, which always identified the other arm. Each IR is named by the chromosome on which it is found and the approximate genomic coordinates (in megabase pairs) of the center of the spacer, for example, IRX-47.303 (Table 1). Exclusion of IRs from the data set was performed as described in Supplemental data S1. Those that were located in unfinished or unassembled regions of the current genome assembly were excluded, for example, spanning or abutting gaps. IRs that share a common center were collapsed into single IRs. In cases in which multiple IRs occurred in complex clusters containing multigene families, for example, the UGT genes on Chromosome 4, or in regions with patterns of secondary structure (Figs. 2 and 3), only the largest most similar IR in the cluster was included. Approximately the same proportion of IRs were excluded from the X- (13) and Y-chromosomes (five) relative to the autosomes (52) as were initially detected by IRF, and thus did not contribute any bias to the overall conclusions of this study (Fig. 1D).

To estimate the number of X-chromosome genes expressed predominantly from testes, we queried the human transcriptome represented in the GNF atlas 2 data base (Su et al. 2004) using the "gene sorter" feature of the UCSC Genome Browser for genes on the X-chromosome with a minimum testes expression ratio of 1.0 and a maximum expression ratio of 1.5 for all other tissues.

Two specific IRs with spacer >100 kb are listed in Table 1. The 283-kb IRY-23.2 has a 169-kb spacer (P3; Skaletsky et al. 2003) that contained the 9.8-kb IRY-23.30274. IRY-23.2 has replaced IRY-23.30274 in Table 1 (Supplemental data S1 and S2). The 29.1-kb IRX-147.36 has a 164.9-kb spacer that contained IRX-147.49. The MAGE-A9 genes were actually found in the arms of IRX-147.36 (Table 1; Supplemental data S2). Both of these IRs were correctly detected when the spacer length was set at 500 kb. Similarity dot-plots were performed using MacVector 7.0 and analyzed with the help of Canvas 5.0. Additional alignments were performed using CLUSTAL 1.81 (Thompson et al. 1994).

PCR was performed using standard protocols. Gorilla and chimpanzee genomic DNA were obtained from Coriell Cell Repositories. DNA sequencing was performed using standard protocols.

	Acknowledgements

 
We thank Laura Carrel (Penn State, Hershey, PA) for comments on the manuscript and Alfredo Rodriquez (Boston University) for helpful advice. This work was supported in part by grants from the NSF DBI-0413462 (to G.B.) and the NIH R21 HG002919 (to P.E.W.).

	Footnotes

 
⁴ Corresponding author.
E-MAIL peter.warburton{at}mssm.edu; FAX (212) 849-2508.

[Supplemental material is available online at www.genome.org.]

Article and publication are at http://www.genome.org/cgi/doi/10.1101/gr.2542904.

	REFERENCES

Top ABSTRACT RESULTS DISCUSSION METHODS REFERENCES WEB SITE REFERENCES

 

Aradhya, S., Bardaro, T., Galgoczy, P., Yamagata, T., Esposito, T., Patlan, H., Ciccodicola, A., Munnich, A., Kenwrick, S., Platzer, M., et al. 2001. Multiple pathogenic and benign genomic rearrangements occur at a 35 kb duplication involving the NEMO and LAGE2 genes. Hum. Mol. Genet. 10: 2557-2567.[Abstract/Free Full Text]

Bailey, J.A., Gu, Z., Clark, R.A., Reinert, K., Samonte, R.V., Schwartz, S., Adams, M.D., Myers, E.W., Li, P.W., and Eichler, E.E. 2002. Recent segmental duplications in the human genome. Science 297: 1003-1007.[Abstract/Free Full Text]

Benham, C.J., Savitt, A.G., and Bauer, W.R. 2002. Extrusion of an imperfect palindrome to a cruciform in superhelical DNA: Complete determination of energetics using a statistical mechanical model. J. Mol. Biol. 316: 563-581.[CrossRef][Medline]

Benson, G. 1999. Tandem Repeats Finder: A program to analyze DNA sequences. Nucleic Acids Res. 27: 573-580.[Abstract/Free Full Text]

Chen, Y.T., Scanlan, M.J., Sahin, U., Tureci, O., Gure, A.O., Tsang, S., Williamson, B., Stockert, E., Pfreundschuh, M., and Old, L.J. 1997. A testicular antigen aberrantly expressed in human cancers detected by autologous antibody screening. Proc. Natl. Acad. Sci. 94: 1914-1918.[Abstract/Free Full Text]

Chomez, P., De Backer, O., Bertrand, M., De Plaen, E., Boon, T., and Lucas, S. 2001. An overview of the MAGE gene family with the identification of all human members of the family. Cancer Res. 61: 5544-5551.[Abstract/Free Full Text]

Eddy, E.M. and O'Brien, D.A. 1998. Gene expression during mammalian meiosis. Curr. Top. Dev. Biol. 37: 141-200.[Medline]

Emerson, J.J., Kaessmann, H., Betran, E., and Long, M. 2004. Extensive gene traffic on the mammalian X chromosome. Science 303: 537-540.[Abstract/Free Full Text]

Gure, A.O., Wei, I.J., Old, L.J., and Chen, Y.T. 2002. The SSX gene family: Characterization of 9 complete genes. Int. J. Cancer 101: 448-453.[CrossRef][Medline]

Jurka, J. 2000. Repbase update: A database and an electronic journal of repetitive elements. Trends Genet. 16: 418-420.[CrossRef][Medline]

Kent, W.J., Sugnet, C.W., Furey, T.S., Roskin, K.M., Pringle, T.H., Zahler, A.M., and Haussler, D. 2002. The Human Genome Browser at UCSC. Genome Res. 12: 996-1006.[Abstract/Free Full Text]

Khil, P.P., Smirnova, N.A., Romanienko, P.J., and Camerini-Otero, R.D. 2004. The mouse X chromosome is enriched for sex-biased genes not subject to selection by meiotic sex chromosome inactivation. Nat. Genet. 36: 642-646.[CrossRef][Medline]

Kuroda-Kawaguchi, T., Skaletsky, H., Brown, L.G., Minx, P.J., Cordum, H.S., Waterston, R.H., Wilson, R.K., Silber, S., Oates, R., Rozen, S., et al. 2001. The AZFc region of the Y chromosome features massive palindromes and uniform recurrent deletions in infertile men. Nat. Genet. 29: 279-286.[CrossRef][Medline]

Lafreniere, R.G., Brown, C.J., Rider, S., Chelly, J., Taillon-Miller, P., Chinault, A.C., Monaco, A.P., and Willard, H.F. 1993. 2.6 Mb YAC contig of the human X inactivation center region in Xq13: Physical linkage of the RPS4X, PHKA1, XIST and DXS128E genes. Hum. Mol. Genet. 2: 1105-1115.[Abstract/Free Full Text]

Lander, E.S., Linton, L.M., Birren, B., Nusbaum, C., Zody, M.C., Baldwin, J., Devon, K., Dewar, K., Doyle, M., FitzHugh, W., et al. 2001. Initial sequencing and analysis of the human genome. Nature 409: 860-921.[CrossRef][Medline]

Lercher, M.J., Urrutia, A.O., and Hurst, L.D. 2003. Evidence that the human X chromosome is enriched for male-specific but not female-specific genes. Mol. Biol. Evol. 20: 1113-1116.[Abstract/Free Full Text]

Marzluff, W.F., Gongidi, P., Woods, K.R., Jin, J., and Maltais, L.J. 2002. The human and mouse replication-dependent histone genes. Genomics 80: 487-498.[CrossRef][Medline]

McDonell, N., Ramser, J., Francis, F., Vinet, M.C., Rider, S., Sudbrak, R., Riesselman, L., Yaspo, M.L., Reinhardt, R., Monaco, A.P., et al. 2000. Characterization of a highly complex region in Xq13 and mapping of three isodicentric breakpoints associated with preleukemia. Genomics 64: 221-229.[CrossRef][Medline]

Ottolenghi, C., Fellous, M., Barbieri, M., and McElreavey, K. 2002. Novel paralogy relations among human chromosomes support a link between the phylogeny of doublesex-related genes and the evolution of sex determination. Genomics 79: 333-343.[CrossRef][Medline]

Pearson, C.E., Zorbas, H., Price, G.B., and Zannis-Hadjopoulos, M. 1996. Inverted repeats, stem-loops, and cruciforms: Significance for initiation of DNA replication. J. Cell Biochem. 63: 1-22.[Medline]

Rozen, S., Skaletsky, H., Marszalek, J.D., Minx, P.J., Cordum, H.S., Waterston, R.H., Wilson, R.K., and Page, D.C. 2003. Abundant gene conversion between arms of palindromes in human and ape Y chromosomes. Nature 423: 873-876.[CrossRef][Medline]

Samonte, R.V. and Eichler, E.E. 2002. Segmental duplications and the evolution of the primate genome. Nat. Rev. Genet. 3: 65-72.[Medline]

Scanlan, M.J., Simpson, A.J.G., and Old, L.J. 2004. The cancer/testis genes: Review, standardization, and commentary. Cancer Imm. 4: 1-15.

Shlyakhtenko, L.S., Hsieh, P., Grigoriev, M., Potaman, V.N., Sinden, R.R., and Lyubchenko, Y.L. 2000. A cruciform structural transition provides a molecular switch for chromosome structure and dynamics. J. Mol. Biol. 296: 1169-1173.[CrossRef][Medline]

Skaletsky, H., Kuroda-Kawaguchi, T., Minx, P.J., Cordum, H.S., Hillier, L., Brown, L.G., Repping, S., Pyntikova, T., Ali, J., Bieri, T., et al. 2003. The male-specific region of the human Y chromosome is a mosaic of discrete sequence classes. Nature 423: 825-837.[CrossRef][Medline]

Small, K., Iber, J., and Warren, S.T. 1997. Emerin deletion reveals a common X-chromosome inversion mediated by inverted repeats. Nat. Genet. 16: 96-99.[CrossRef][Medline]

Su, A.I., Wiltshire, T., Batalov, S., Lapp, H., Ching, K.A., Block, D., Zhang, J., Soden, R., Hayakawa, M., Kreiman, G., et al. 2004. A gene atlas of the mouse and human protein-encoding transcriptomes. Proc. Natl. Acad. Sci. 101: 6062-6067.[Abstract/Free Full Text]

Thompson, J.D., Higgins, D.G., and Gibson, T.J. 1994. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22: 4673-4680.[Abstract/Free Full Text]

Wang, P.J. 2004. X chromosomes, retrogenes and their role in male reproduction. Trends Endocrin. Metab. 15: 79-83.[CrossRef][Medline]

Wang, P.J., McCarrey, J.R., Yang, F., and Page, D.C. 2001. An abundance of X-linked genes expressed in spermatogonia. Nat. Genet. 27: 422-426.[CrossRef][Medline]

Waterston, R.H., Lindblad-Toh, K., Birney, E., Rogers, J., Abril, J.F., Agarwal, P., Agarwala, R., Ainscough, R., Alexandersson, M., An, P., et al. 2002. Initial sequencing and comparative analysis of the mouse genome. Nature 420: 520-562.[CrossRef][Medline]

Wayne, C.M., MacLean, J.A., Cornwall, G., and Wilkinson, M.F. 2002. Two novel human X-linked homeobox genes, hPEPP1 and hPEPP2, selectively expressed in the testis. Gene 301: 1-11.[CrossRef][Medline]

Wolff, D.J., Miller, A.P., Van Dyke, D.L., Schwartz, S., and Willard, H.F. 1996. Molecular definition of breakpoints associated with human Xq isochromosomes: Implications for mechanisms of formation. Am. J. Hum. Genet. 58: 154-160.[Medline]

Wu, C.I. and Xu, E.Y. 2003. Sexual antagonism and X inactivation—The SAXI hypothesis. Trends Genet. 19: 243-247.[CrossRef][Medline]

Yoshikawa, T., Seki, N., Azuma, T., Masuho, Y., Muramatsu, M., Miyajima, N., and Saito, T. 2000. Isolation of a cDNA for a novel human RING finger protein gene, RNF18, by the virtual transcribed sequence (VTS) approach(1). Biochim. Biophys. Acta 1493: 349-355.[Medline]

Zalensky, A.O., Siino, J.S., Gineitis, A.A., Zalenskaya, I.A., Tomilin, N.V., Yau, P., and Bradbury, E.M. 2002. Human testis/sperm-specific histone H2B (hTSH2B). Molecular cloning and characterization. J. Biol. Chem. 277: 43474-43480.[Abstract/Free Full Text]

Zendman, A.J., Ruiter, D.J., and Van Muijen, G.N. 2003a. Cancer/testis-associated genes: Identification, expression profile, and putative function. J. Cell Physiol. 194: 272-288.[CrossRef][Medline]

Zendman, A.J., Zschocke, J., van Kraats, A.A., de Wit, N.J., Kurpisz, M., Weidle, U.H., Ruiter, D.J., Weiss, E.H., and van Muijen, G.N. 2003b. The human SPANX multigene family: Genomic organization, alignment and expression in male germ cells and tumor cell lines. Gene 309: 125-133.[CrossRef][Medline]

	WEB SITE REFERENCES

Top ABSTRACT RESULTS DISCUSSION METHODS REFERENCES WEB SITE REFERENCES

 

http://ftp.genome.washington.edu/RM/RepeatMasker.html; RepeatMasker.

http://tandem.bu.edu/cgi-bin/irdb/irdb.exe; Inverted Repeat Finder.

http://tandem.bu.edu; Inverted Repeat Data Base (IRDB).

Received March 4, 2004; accepted in revised format July 27, 2004.

CiteULike

Connotea

Del.icio.us

Digg

Reddit

Technorati

This article has been cited by other articles:

				M. Caceres, National Institutes of Health Intramural Sequencin, R. T. Sullivan, and J. W. Thomas A recurrent inversion on the eutherian X chromosome PNAS, November 20, 2007; 104(47): 18571 - 18576. [Abstract] [Full Text] [PDF]

				H. Tanaka, Y. Cao, D. A. Bergstrom, C. Kooperberg, S. J. Tapscott, and M.-C. Yao Intrastrand Annealing Leads to the Formation of a Large DNA Palindrome and Determines the Boundaries of Genomic Amplification in Human Cancer Mol. Cell. Biol., March 15, 2007; 27(6): 1993 - 2002. [Abstract] [Full Text] [PDF]

				S. Rajender, K. Thangaraj, N. J. Gupta, N. Leelavathy, D. S. Rani, R. G. Nambiar, V. Kalavathy, S. T. Santhiya, S. Rajangam, P. M. Gopinath, et al. A Novel Human Sex-Determining Gene Linked to Xp11.21-11.23 J. Clin. Endocrinol. Metab., October 1, 2006; 91(10): 4028 - 4036. [Abstract] [Full Text] [PDF]

				A. Bacolla, J. R. Collins, B. Gold, N. Chuzhanova, M. Yi, R. M. Stephens, S. Stefanov, A. Olsh, J. P. Jakupciak, M. Dean, et al. Long homopurine*homopyrimidine sequences are characteristic of genes expressed in brain and the pseudoautosomal region. Nucleic Acids Res., January 1, 2006; 34(9): 2663 - 2675. [Abstract] [Full Text] [PDF]

				N. Kouprina, A. Pavlicek, V. N. Noskov, G. Solomon, J. Otstot, W. Isaacs, J. D. Carpten, J. M. Trent, J. Schleutker, J. C. Barrett, et al. Dynamic structure of the SPANX gene cluster mapped to the prostate cancer susceptibility locus HPCX at Xq27 Genome Res., November 1, 2005; 15(11): 1477 - 1486. [Abstract] [Full Text] [PDF]

				P. J.I. Ellis, E. J. Clemente, P. Ball, A. Toure, L. Ferguson, J. M.A. Turner, K. L. Loveland, N. A. Affara, and P. S. Burgoyne Deletions on mouse Yq lead to upregulation of multiple X- and Y-linked transcripts in spermatids Hum. Mol. Genet., September 15, 2005; 14(18): 2705 - 2715. [Abstract] [Full Text] [PDF]

This Article Abstract Full Text (PDF) Supplemental Research Data Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend