Discovering isochores by least-squares optimal segmentation

The isochore structure of a genome is observable by variation in the G+C (guanine and cytosine) content within and between the chromosomes. Describing the isochore structure of vertebrate genomes is a challenging task, and many computational methods have been developed and applied to it. Here we apply a well-known least-squares optimal segmentation algorithm to isochore discovery. The algorithm finds the best division of the sequence into k pieces, such that the segments are internally as homogeneous as possible. We show how this simple segmentation method can be applied to isochore discovery using as input the G+C content of sliding windows on the sequence. To evaluate the performance of this segmentation technique on isochore detection, we present results from segmenting previously studied isochore regions of the human genome. Detailed results on the MHC locus, on parts of chromosomes 21 and 22, and on a 100 Mb region from chromosome 1 are similar to previously suggested isochore structures. We also give results on segmenting all 22 autosomal human chromosomes. An advantage of this technique is that oversegmentation of G+C rich regions can generally be avoided. This is because the technique concentrates on greater global, instead of smaller local, differences in the sequence composition. The effect is further emphasized by a log-transformation of the data that lowers the high variance that is observed in G+C rich regions. We conclude that the least-squares optimal segmentation method is computationally efficient and yields results close to previous biologically motivated isochore structures.     

Biological implications of isochore boundaries in the human genome

The human genome is composed of large sequence segments with fairly homogeneous GC content, namely isochores, which have been linked to many important functions; biological implications of most isochore boundaries, however, remain elusive, partly due to the difficulty in determining these boundaries at high resolution. Using the segmentation algorithm based on the quadratic divergence, we re-determined all 79 boundaries of previously identified human isochores at single-nucleotide resolution, and then compared the boundary coordinates with other genome features. We found that 55.7% of isochore boundaries coincide with termini of repeat elements; 45.6% of isochore boundaries coincide with termini of highly conserved sequences based on alignment of 17 vertebrate genomes, i.e., the highly conserved genome sequence switches to a less or non-conserved one at the isochore boundary; some isochore boundaries coincide with abrupt change of CpG island distribution (note that one boundary can associate with more than one genome feature). In addition, sequences around isochore boundaries are highly conserved. It seems reasonable to deduce that the boundaries of all the isochores studied here would be replication timing sites in the human genome. These results suggest possible key roles of the isochore boundaries and may further our understanding of the human genome organization.     

Biological Implications of Isochore Boundaries in the Human Genome

Abstract

The human genome is composed of large sequence segments with fairly homogeneous GC content, namely isochores, which have been linked to many important functions; biological implications of most isochore boundaries, however, remain elusive, partly due to the difficulty in determining these boundaries at high resolution. Using the segmentation algorithm based on the quadratic divergence, we re-determined all 79 boundaries of previously identified human isochores at single-nucleotide resolution, and then compared the boundary coordinates with other genome features. We found that 55.7% of isochore boundaries coincide with termini of repeat elements; 45.6% of isochore boundaries coincide with termini of highly conserved sequences based on alignment of 17 vertebrate genomes, i.e., the highly conserved genome sequence switches to a less or non-conserved one at the isochore boundary; some isochore boundaries coincide with abrupt change of CpG island distribution (note that one boundary can associate with more than one genome feature). In addition, sequences around isochore boundaries are highly conserved. It seems reasonable to deduce that the boundaries of all the isochores studied here would be replication timing sites in the human genome. These results suggest possible key roles of the isochore boundaries and may further our understanding of the human genome organization.     

Evolution of isochores in rodents

The most deviant isochore pattern within mammals was found in rat and mouse; most other mammals possess a different kind of isochore organization called the "general pattern." However, isochore patterns remain largely unknown in rodents other than mouse and rat. To investigate the taxonomic distribution of isochore patterns in rodents, we sequenced the nuclear gene LCAT (lecithin:cholesterol acyltransferase) from 17 rodents species (bringing the total of LCAT sequences in rodent to 19) and compared their GC contents at third codon positions and in introns. We also analyzed an extensive sequence database from rodents other than rat and mouse. All murid LCAT sequences are much poorer in GC than all nonrodent LCAT sequences, and the hamster sequence database shows exactly the same isochore pattern as rat and mouse. Thus, all murids share the same special isochore pattern--GC homogenization. LCAT sequences are GC-poor in hystricomorphs too, but the guinea pig sequence database indicates that large changes in GC content occur without an overall modification of the isochore pattern. This novel mode of isochore evolution is called GC reordering. LCAT sequences also show that the evolution of isochores in sciurids and glirids is nonconservative in comparison with that in nonrodents. Thus, at least two novel patterns of isochore evolution were found. No rodent investigated to date shared the general mammalian pattern.      

Genetics of isoprenoid biosynthesis in Paracoccus zeaxanthinifaciens

Three statistical/mathematical analyses are carried out on isochore sequences: spectral analysis, analysis of variance, and segmentation analysis. Spectral analysis shows that there are GC content fluctuations at different length scales in isochore sequences. The analysis of variance shows that the null hypothesis (the mean value of a group of GC contents remains the same along the sequence) may or may not be rejected for an isochore sequence, depending on the subwindow sizes at which GC contents are sampled, and the window size within which group members are defined. The segmentation analysis shows that there are stronger indications of GC content changes at isochore borders than within an isochore. These analyses support the notion of isochore sequences, but reject the assumption that isochore sequences are homogeneous at the base level. An isochore sequence may pass a homogeneity test when GC content fluctuations at smaller length scales are ignored or averaged out.     

Evolution of the isochore structure in the scale of chromosome: insight from the mutation bias and fixation bias

Following the development of reliable methods for inferring the direction of mutations of the single nucleotide polymorphism (SNP), and the revealing of the human isochore map, it has become possible to investigate the evolution of the isochore structure in a continuous region. In this study, the recent evolution of the isochore structure on human chromosome 18, as inferred from the SNP, was examined. A remarkable mutation bias was found, which was destroying the present isochore structure. However, a fixation bias contributed by the biased gene conversion (BGC) effect and a rising fixation probability of derived alleles with increasing GC content was extending the present isochore structure. Combining the two opposing processes, the old isochore structure was declining and a more homogenous isochore structure with higher GC content was being formed on the chromosome. During this process, both the CpG and genic sites, which were present in the isochore but were paid little attention to before, played an important role. In addition, the recombination was confirmed to promote the GC alleles fixed in the genome because of the BGC effect. For the first time, it was observed that with the occurrence of little recombination, AT alleles had the identical fixation probability with GC alleles in the recombination cold spots.     

Isochore chromosome maps of the human genome

The human genome is a mosaic of isochores, which are long DNA segments (300 kbp) relatively homogeneous in G+C. Human isochores were first identified by density-gradient ultracentrifugation of bulk DNA, and differ in important features, e.g. genes are found predominantly in the GC-richest isochores. Here, we use a reliable segmentation method to partition the longest contigs in the human genome draft sequence into long homogeneous genome regions (LHGRs), thereby revealing the isochore structure of the human genome. The advantages of the isochore maps presented here are: (1) sequence heterogeneities at different scales are shown in the same plot; (2) pair-wise compositional differences between adjacent regions are all statistically significant; (3) isochore boundaries are accurately defined to single base pair resolution; and (4) both gradual and abrupt isochore boundaries are simultaneously revealed. Taking advantage of the wide sample of genome sequence analyzed, we investigate the correspondence between LHGRs and true human isochores revealed through DNA centrifugation. LHGRs show many of the typical isochore features, mainly size distribution, G+C range, and proportions of the isochore classes. The relative density of genes, Alu and long interspersed nuclear element repeats and the different types of single nucleotide polymorphisms on LHGRs also coincide with expectations in true isochores. Potential applications of isochore maps range from the improvement of gene-finding algorithms to the prediction of linkage disequilibrium levels in association studies between marker genes and complex traits. The coordinates for the LHGRs identified in all the contigs longer than 2 Mb in the human genome sequence are available at the online resource on isochore mapping: http://bioinfo2.ugr.es/isochores.     

Isochores Merit the Prefix 'Iso'

The isochore concept in human genome sequence was challenged in an analysis by the International Human Genome Sequencing Consortium (IHGSC). We argue here that a statement in IGHSC analysis concerning the existence of isochore is incorrect, because it had applied an inappropriate statistical test. To test the existence of isochores should be equivalent to a test of homogeneity of windowed GC%. The statistical test applied in the IHGSC's analysis, the binomial test, is however a test of a sequence being random on the base level. For testing the existence of isochore, or homogeneity in GC%, we propose to use another statistical test: the analysis of variance (ANOVA). It can be shown that DNA sequences that are rejected by binomial test may not be rejected by the ANOVA test.     

Isochore patterns and gene distributions in fish genomes

The compositional approach developed in our laboratory many years ago revealed a large-scale compositional heterogeneity in vertebrate genomes, in which GC-rich and GC-poor regions, the isochores, were found to be characterized by high and low gene densities, respectively. Here we mapped isochores on fish chromosomes and assessed gene densities in isochore families. Because of the availability of sequence data, we have concentrated our investigations on four species, zebrafish (Brachydanio rerio), medaka (Oryzias latipes), stickleback (Gasterosteus aculeatus), and pufferfish (Tetraodon nigroviridis), which belong to four distant orders and cover almost the entire GC range of fish genomes. These investigations produced isochore maps that were drastically different not only from those of mammals (in that only two major isochore families were essentially present in each genome vs five in the human genome) but also from each other (in that different isochore families were represented in different genomes). Gene density distributions for these fish genomes were also obtained and shown to follow the expected increase with increasing isochore GC. Finally, we discovered a remarkable conservation of the average size of the isochores (which match replicon clusters in the case of human chromosomes) and of the average GC levels of isochore families in both fish and human genomes. Moreover, in each genome the GC-poorest isochore families comprised a group of "long isochores" (2-20 Mb in size), which were the lowest in GC and varied in size distribution and relative amount from one genome to the other.     

Comparative Compositional Mapping of Chicken and Quail Chromosomes

The distribution of various isochore families on mitotic chromosomes of domestic chicken and Japanese quail was studied by the method of fluorescence in situ DNA–DNA hybridization (FISH). DNA of various isochore families was shown to be distributed irregularly and similarly on chromosomes of domestic chicken and Japanese quail. The GC-rich isochore families (H2, H3, and H4) hybridized mainly to microchromosomes and a majority of macrochromosome telomeric regions. In chicken, an intense fluorescence was also in a structural heterochromatin region of the Z chromosome long arm. In some regions of the quail macrochromosome arms, hybridization was also with isochore families H3 and H4. On macrochromosomes of both species, the pattern of hybridization with isochores of the H2 and H3 families resembled R-banding. The light isochores (L1 and L2 families) are mostly detected within macrochromosome internal regions corresponding to G bands, whereas microchromosomes lack light isochores. Although mammalian and avian karyotypes differ significantly in organization, the isochore distribution in genomes of these two lineages of the warm-blooded animals is similar in principle. On macrochromosomes of the two avian species studied, a pattern of isochore distribution resembled that of mammalian chromosomes. The main specific feature of the avian genome, a great number of microchromosomes (about 30% of the genome), determines a compositional specialization of the latter. This suggests the existence of not only structural but also functional compartmentalization of the avian genome.     

Recombination explains isochores in mammalian genomes

The mouse Fxy gene was translocated into the highly recombining pseudoautosomal region comparatively recently in evolutionary terms. This event resulted in a rapid increase of GC content. We investigated the consequences of the translocation further by sequencing exons and introns of Fxy in various rodent species. We found that the DNA fragment newly located in a highly recombining context has acquired every property of a GC-rich isochore, namely increased GC content (especially at the third codon positions of exons), shorter introns and high density of minisatellites. These results strongly suggest that recombination is the primary determinant of the isochore organization of mammalian genomes.     

Single-copy sequence homology among the GC-richest isochores of the genomes from warm-blooded vertebrates

We have hybridized a human DNA fraction corresponding to the GC-richest and gene-richest isochore family, H3, on compositional fractions of DNAs from 12 mammalian species and three avian species, representing eight and three orders, respectively. Under conditions in which repetitive sequences are competed out, the H3 isochore probe only or predominantly hybridized on the GC-richest fractions of main-band DNA from all the species investigated. These results indicate that single-copy sequences from the human H3 isochores share homology with sequences located in the compositionally corresponding compartments of the vertebrate genomes tested. These sequences are likely to be essentially formed by conserved coding sequences. The present results add to other lines of evidence indicating that isochore patterns are highly conserved in warm-blooded vertebrate genomes. Moreover, they refine recent reports (Sabeur et al., 1993; Kadi et al., 1993), and correct them in some details and also in demonstrating that the shrew genome does not exhibit the general mammalian pattern, but a special pattern.Correspondence to: G. Bernardi     

An isochore transition zone in the NF1 gene region is a conserved landmark of chromosome structure and function

The mammalian genome is organized as a mosaic of isochores, stretches of DNA with a distinct sequence composition. Isochores form the basis of the chromosomal banding pattern, which is tightly correlated with a number of structural and functional features. We have recently demonstrated that the transition from a GC-poor isochore to a GC-rich one in the NF1 gene region occurs within 5 kb and demarcates genomic regions with high and low recombination frequency. We now report that the same transition zone separates early replicating from late replicating chromatin on the molecular level. At the isochore transition the replication fork is stalled in mid-S phase and can be visualized by fiber-FISH techniques as a Y-shaped structure. The switch in GC content and in replication timing is conserved between human and mouse, emphasizing the importance of the transition zones as landmarks of chromosome organization and function.     

Localization of DNA sequences tightly associated with synaptonemal complex in compositional fractions of golden hamster

The synaptonemal complex isolated from the spermatocyte nuclei by exhaustive hydrolysis of the latter by DNase II contains tightly associated DNA sequences (SCAR DNA). Here we studied the compositional properties of a cloned family of SCAR DNA of golden hamster, namely we performed the localization of 27 SCAR DNA clones on compositionally fractionated genomic DNA from golden hamster. We observed that sequences of the SCAR DNA family are mainly localized in the GC-poor isochore families L1 and L2, that showed 63% hybridization signals. This means that 37% of signals is referred to the GC-rich isochores, indicating the presence of SCAR DNA overall the genome, even if each isochore family presents differences in density and sequence type. Moreover, the SCAR DNA sequences containing regions of homology with LINE/SINE repeats were observed in all the isochore families. The compositional localization of SCAR DNA is in agreement with the hypothesis that SC and SCAR DNA participate in the chromatin organization during the meiosis prophase I, which should result in the attachment of chromatin loops to lateral elements of SC along the whole length of the latter.     

Comparative compositional mapping of chicken and quail chromosomes

The distribution of various isochore families on mitotic chromosomes of domestic chicken and Japanese quail was studied by the method of fluorescence in situ DNA--DNA hybridization (FISH). DNA of various isochore families was shown to be distributed irregularly and similarly on chromosomes of domestic chicken and Japanese quail. The GC-rich isochore families (H2, H3, and H4) hybridized mainly to microchromosomes and a majority of macrochromosome telomeric regions. In chicken, an intense fluorescence was also in a structural heterochromatin region of the Z chromosome long arm. In some regions of the quail macrochromosome arms, hybridization was also with isochore families H3 and H4. On macrochromosomes of both species, the pattern of hybridization with isochores of the H2 and H3 families resembled R-banding. The light isochores (L1 and L2 families) are mostly detected within macrochromosome internal regions corresponding to G bands, whereas microchromosomes lack light isochores. Although mammalian and avian karyotypes differ significantly in organization, the isochore distribution in genomes of these two lineages of the warm-blooded animals is similar in principle. On macrochromosomes of the two avian species studied, a pattern of isochore distribution resembled that of mammalian chromosomes. The main specific feature of the avian genome, a great number of microchromosomes (about 30% of the genome), determines a compositional specialization of the latter. This suggests the existence of not only structural but also functional compartmentalization of the avian genome.     

Isochores and replication time zones: a perfect match

The mammalian genome is not a random sequence but shows a specific, evolutionarily conserved structure that becomes manifest in its isochore pattern. Isochores, i.e. stretches of DNA with a distinct sequence composition and thus a specific GC content, cause the chromosomal banding pattern. This fundamental level of genome organization is related to several functional features like the replication timing of a DNA sequence. GC richness of genomic regions generally corresponds to an early replication time during S phase. Recently, we demonstrated this interdependency on a molecular level for an abrupt transition from a GC-poor isochore to a GC-rich one in the NF1 gene region; this isochore boundary also separates late from early replicating chromatin. Now, we analyzed another genomic region containing four isochores separated by three sharp isochore transitions. Again, the GC-rich isochores were found to be replicating early, the GC-poor isochores late in S phase; one of the replication time zones was discovered to consist of one single replicon. At the boundaries between isochores, that all show no special sequence elements, the replication machinery stopped for several hours. Thus, our results emphasize the importance of isochores as functional genomic units, and of isochore transitions as genomic landmarks with a key function for chromosome organization and basic biological properties.     

Isochore chromosome maps of eukaryotic genomes

Analytical DNA ultracentrifugation revealed that eukaryotic genomes are mosaics of isochores: long DNA segments (>300 kb on average) relatively homogeneous in G+C. Important genome features are dependent on this isochore structure, e.g. genes are found predominantly in the GC-richest isochore classes. However, no reliable method is available to rigorously partition the genome sequence into relatively homogeneous regions of different composition, thereby revealing the isochore structure of chromosomes at the sequence level. Homogeneous regions are currently ascertained by plain statistics on moving windows of arbitrary length, or simply by eye on G+C plots. On the contrary, the entropic segmentation method is able to divide a DNA sequence into relatively homogeneous, statistically significant domains. An early version of this algorithm only produced domains having an average length far below the typical isochore size. Here we show that an improved segmentation method, specifically intended to determine the most statistically significant partition of the sequence at each scale, is able to identify the boundaries between long homogeneous genome regions displaying the typical features of isochores. The algorithm precisely locates classes II and III of the human major histocompatibility complex region, two well-characterized isochores at the sequence level, the boundary between them being the first isochore boundary experimentally characterized at the sequence level. The analysis is then extended to a collection of human large contigs. The relatively homogeneous regions we find show many of the features (G+C range, relative proportion of isochore classes, size distribution, and relationship with gene density) of the isochores identified through DNA centrifugation. Isochore chromosome maps, with many potential applications in genomics, are then drawn for all the completely sequenced eukaryotic genomes available.     

Isochore pattern and gene distribution in the chicken genome

We report here investigations on the isochore pattern and the distribution of genes in the chromosomes of chicken. In spite of large differences in genome size and karyotype, the compositional properties and the gene distribution of the chicken genome are very similar to those recently published for the human genome, which is a good representative of most mammalian genomes. In fact, this similarity, which extends to the relative amounts and, also, to a large extent at least, to the average base composition of isochore families, is most interesting in view of the very large distance of mammals and birds for a common ancestor, which goes back to 310–340 million years ago. This raises important questions about genome evolution in vertebrates.     

Isochore evolution in mammals: a human-like ancestral structure.

Codon usage in mammals is mainly determined by the spatial arrangement of genomic G + C-content, i.e., the isochore structure. Ancestral G + C-content at third codon positions of 27 nuclear protein-coding genes of eutherian mammals was estimated by maximum-likelihood analysis on the basis of a nonhomogeneous DNA substitution model, accounting for variable base compositions among present-day sequences. Data consistently supported a human-like ancestral pattern, i.e., highly variable G + C-content among genes. The mouse genomic structure-more narrow G + C-content distribution-would be a derived state. The circumstances of isochore evolution are discussed with respect to this result. A possible relationship between G + C-content homogenization in murid genomes and high mutation rate is proposed, consistent with the negative selection hypothesis for isochore maintenance in mammals.