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.     

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.     

The length of chromatin loops in meiotic prophase I of warm-blooded vertebrates depends on the DNA compositional organization

In meiotic prophase I, chromatin fibrils attached to the lateral elements of the synaptonemal complexes (SC) form loops. SCAR DNA (synaptonemal complex associated regions of DNA) is a family of genomic DNA tightly associated with the SC and located at the chromatin loop basements. Using the hybridization technique, it was demonstrated that localization of SCAR DNA was evolutionarily conserved in the isochore compositional fractions of the three examined genomes of warm-blooded vertebrates—human, chicken, and golden hamster. The introduction of the concept of the comparative loops (CL) of DNA that form of chromatin attach to SC in the isochore compositional fractions provided the calculation of their length. An inverse proportional relationship between the length of CL DNA and the GC level in the isochore compartments of the studied warm-blooded vertebrate genomes was revealed. An exception was the GCpoorest L1 isochore family. For different compositional isochores of the human and chicken genomes, the number of genes in the CL DNA was evaluated. A model of the formation of GC-rich isochores in vertebrate genomes, according to which there was not only an increase in the GC level but also the elimination of functionally insignificant noncoding DNA regions, as well as joining of isochores decreasing in size, was suggested.     

A compositional map of the cen-q21 region of human chromosome 21

A compositional map of the centromere and of the subcentromeric region of the long arm of human chromosome 21 was established by determining the GC levels (GC is the molar fraction of guanine+cytosine in DNA) of 11 YACs (yeast artificial chromosomes) covering this 13–14 Mb region which extends from the α-satellite sequences of the C(entromeric) band qll.1, through R(everse) band q11.2, to the proximal part of G(iemsa) band q21. The entire region is made up of GC-poor, or L, isochores with only one GC-rich H1 isochore, at least 2 Mb in size, located in band q21. The almost identical GC levels of the centromeric α-satellite repeats (38.5%), of R band q11.2 (39%), and of G bands (38–40%) provide a direct demonstration that base composition cannot be the only cause of the cytogenetic differences between C, G, and the majority of R bands, namely the H3^- R bands (which do not contain the GC-richest H3 isochores). The results obtained also show that isochores may be as long as 6 Mb, at least in the GC-poor regions of the genome, and support previous observations suggesting that YACs from isochore borders are unstable and/or difficult to clone. Genes and CpG islands are very rare in the GC-poor region investigated, as expected from the fact that their concentration is proportional to the GC levels of the isochores in which they are contained.     

A compositional map of human chromosome 21.

GC-poor and GC-rich isochores, the long (greater than 300 kb) compositionally homogeneous DNA segments that form the genome of warm-blooded vertebrates, are located in G- and R-bands respectively of metaphase chromosomes. The precise correspondence between GC-rich isochores and R-band structure is still, however, an open problem, because GC-rich isochores are compositionally heterogeneous and only represent one-third of the genome, with the GC-richest family (which is by far the highest in gene concentration) corresponding to less than 5% of the genome. In order to clarify this issue and, more generally, to correlate DNA composition and chromosomal structure in an unequivocal way, we have developed a new approach, compositional mapping. This consists of assessing the base composition over 0.2-0.3 Mb (megabase) regions surrounding landmarks that were previously localized on the physical map. Compositional mapping was applied here to the long arm of human chromosome 21, using 53 probes that had already been used in physical mapping. The results obtained provide a direct demonstration that the DNA stretches of G-bands essentially correspond to GC-poor isochores, and that R-band DNA is characterized by a compositional heterogeneity that is much more striking than expected, in that it comprises isochores covering the full spectrum of GC levels. GC-poor isochores of R-bands may, however, correspond to 'thin' G-bands, as visualized at high resolution, leaving GC-rich and very GC-rich isochores as the real components of (high-resolution) R-band DNA.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

Turtle isochore structure is intermediate between amphibians and other amniotes

Vertebrate genomes are comprised of isochores that are relatively long (>100 kb) regions with a relatively homogenous (either GC-rich or AT-rich) base composition and with rather sharp boundaries with neighboring isochores. Mammals and living archosaurs (birds and crocodilians) have heterogeneous genomes that include very GC-rich isochores. In sharp contrast, the genomes of amphibians and fishes are more homogeneous and they have a lower overall GC content. Because DNA with higher GC content is more thermostable, the elevated GC content of mammalian and archosaurian DNA has been hypothesized to be an adaptation to higher body temperatures. This hypothesis can be tested by examining structure of isochores across the reptilian clade, which includes the archosaurs, testudines (turtles), and lepidosaurs (lizards and snakes), because reptiles exhibit diverse body sizes, metabolic rates, and patterns of thermoregulation. This study focuses on a comparative analysis of a new set of expressed genes of the red-eared slider turtle and orthologs of the turtle genes in mammalian (human, mouse, dog, and opossum), archosaurian (chicken and alligator), and amphibian (western clawed frog) genomes. EST (expressed sequence tag) data from a turtle cDNA library enriched for genes that have specialized functions (developmental genes) revealed using the GC content of the third-codon-position to examine isochore structure requires careful consideration of the types of genes examined. The more highly expressed genes (e.g., housekeeping genes) are more likely to be GC-rich than are genes with specialized functions. However, the set of highly expressed turtle genes demonstrated that the turtle genome has a GC content that is intermediate between the GC-poor amphibians and the GC-rich mammals and archosaurs. There was a strong correlation between the GC content of all turtle genes and the GC content of other vertebrate genes, with the slope of the line describing this relationship also indicating that the isochore structure of turtles is intermediate between that of amphibians and other amniotes. These data are consistent with some thermal hypotheses of isochore evolution, but we believe that the credible set of models for isochore evolution still includes a variety of models. These data expand the amount of genomic data available from reptiles upon which future studies of reptilian genomics can build.     

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.     

Isochores and the evolutionary genomics of vertebrates

The nuclear genomes of vertebrates are mosaics of isochores, very long stretches (>300kb) of DNA that are homogeneous in base composition and are compositionally correlated with the coding sequences that they embed. Isochores can be partitioned in a small number of families that cover a range of GC levels (GC is the molar ratio of guanine+cytosine in DNA), which is narrow in cold-blooded vertebrates, but broad in warm-blooded vertebrates. This difference is essentially due to the fact that the GC-richest 10-15% of the genomes of the ancestors of mammals and birds underwent two independent compositional transitions characterized by strong increases in GC levels. The similarity of isochore patterns across mammalian orders, on the one hand, and across avian orders, on the other, indicates that these higher GC levels were then maintained, at least since the appearance of ancestors of warm-blooded vertebrates. After a brief review of our current knowledge on the organization of the vertebrate genome, evidence will be presented here in favor of the idea that the generation and maintenance of the GC-richest isochores in the genomes of warm-blooded vertebrates were due to natural selection.     

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.     

Human chromosomal bands: nested structure,high-definition map and molecular basis

In this paper, we report investigations on the nested structure, the high-definition mapping, and the molecular basis of the classical Giemsa and Reverse bands in human chromosomes. We found the rules according to which the approximately 3,200 isochores of the human genome are assembled in high (850-band) resolution bands, and the latter in low (400-band) resolution bands, so forming the nested mosaic structure of chromosomes. Moreover, we identified the borders of both sets of chromosomal bands at the DNA sequence level on the basis of our recent map of isochores, which represent the highest-resolution, ultimate bands. Indeed, beyond the 100-kb resolution of the isochore map, the guanine and cytosine (GC) profile of DNA becomes turbulent owing to the contribution of specific sequences such as exons, introns, interspersed repeats, CpG islands, etc. The isochore-based level of definition (100 kb) of chromosomal bands is much higher than the cytogenetic definition level (2-3 Mb). The major conclusions of this work concern the high degree of order found in the structure of chromosomal bands, their mapping at a high definition, and the solution of the long-standing problem of the molecular basis of chromosomal bands, as these could be defined on the basis of compositional DNA properties alone.     

No isochores in the human chromosomes 21 and 22?

The human genome is described in the literature as being composed of the isochores, i.e., long (hundreds of kilobases) segments with a homogeneous (G + C) content. We calculated the (G + C) content variations along the DNA molecules of the human chromosomes 21 and 22 and found the variations to be higher everywhere compared to the randomized sequences. Hence the (G + C) content is certainly not homogeneous on the isochore scale in the two human chromosomes. In addition, we found no significant difference between the two human molecules and the genome of E. coli regarding the (G + C) content variations. Hence no isochores are either present in the DNA molecules of the human chromosomes 21 and 22, or the isochores are also present in the genome of Escherichia coli. In any case, the present communication demonstrates that the isochores should be defined in unambiguous molecular terms if they are to be used for an up-to-date genome structure characterization.     

The pig genome: compositional analysis and identification of the gene-richest regions in chromosomes and nuclei

The isochore organization of the mammalian genome comprises a general pattern and some special patterns, the former being characterized by a wider compositional distribution of the DNA fragments. The large majority of the mammalian genomes belong to the former, and only some groups, such as the Myomorpha sub-order of Rodentia, belong to the latter. Here we describe the compositional organization of the pig (Sus scrofa) genome that belongs to the general mammalian pattern. We investigated (i) the compositional distribution of the genes by analysis of their GC3 levels (the GC levels at the third codon positions), and (ii) the correlation between the GC3 value of orthologous genes from pig and other vertebrates (human, calf, mouse, chicken, and Xenopus). As expected, the highest gene concentration corresponded to the H3 isochore family, and the highest GC3 correlations were observed in the pig/human and pig/calf comparisons. Then we identified, by in situ hybridization of the GC-richest H3 isochores, the pig chromosomal regions endowed by the highest gene-density that largely corresponded to the telomeric chromosomal bands. Moreover, we observed that these gene-rich bands are syntenic with the previously identified GC-richest/gene richest H3+ bands of the human chromosomes. At the cell nucleus level, we observed that the gene-dense region corresponded to the more internal compartment, as previously found in human and avian cell nuclei.     

An isochore transition in the NF1 gene region coincides with a switch in the extent of linkage disequilibrium

Whole-genome association studies will be a powerful tool to identify genes responsible for common human diseases. A crucial task for association-mapping studies is the evaluation of the relationship between linkage disequilibrium (LD) and physical distance for the genomic region under study. Since it is known that the extent of LD is nonuniformly distributed throughout the human genome, the required marker density has to be determined specifically for the region under study. These regions may be related to isochores and chromosomal bands, as indicated by earlier cytogenetic findings concerning chiasma distribution in meiosis. Therefore we analyzed the neurofibromatosis type 1 (NF1) gene region on chromosome 17q11.2, which is characterized by a nonuniform LD pattern and an L1-to-H2 isochore transition. Long-range LD within the NF1 gene was found to extend over 200 kb (D' = 0.937) in the L1 isochore, whereas, in the neighboring H2 isochore, no LD is apparent between markers spaced by 26 kb (D' = 0.144). Recombination frequencies derived from the LD are at.00019 (high LD) and.01659 (low LD) per megabase, the latter identical to the average value from segregation analysis. The boundary between these regions coincides precisely with a transition in the GC content of the sequences, with low values (37.2%) in the region with long-range LD and high values (51%) in the other. Our results suggest a correlation between the LD pattern and the isochores, at least in the NF1 region. If this correlation can be generalized, the marker densities required for association studies have to be adjusted to the regional GC content and may be chosen according to the isochores.     

Patterns of Vertebrate Isochore Evolution Revealed by Comparison of Expressed Mammalian,Avian, and Crocodilian Genes

Vertebrate genomes are mosaics of isochores, defined as long (>100 kb) regions with relatively homogeneous within-region base composition. Birds and mammals have more GC-rich isochores than amphibians and fish, and the GC-rich isochores of birds and mammals have been suggested to be an adaptation to homeothermy. If this hypothesis is correct, all poikilothermic (cold-blooded) vertebrates, including the nonavian reptiles, are expected to lack a GC-rich isochore structure. Previous studies using various methods to examine isochore structure in crocodilians, turtles, and squamates have led to different conclusions. We collected more than 6000 expressed sequence tags (ESTs) from the American alligator to overcome sample size limitations suggested to be the fundamental problem in the previous reptilian studies. The alligator ESTs were assembled and aligned with their human, mouse, chicken, and western clawed frog orthologs, resulting in 366 alignments. Analyses of third-codon-position GC content provided conclusive evidence that the poikilothermic alligator has GC-rich isochores, like homeothermic birds and mammals. We placed these results in a theoretical framework able to unify available models of isochore evolution. The data collected for this study allowed us to reject the models that explain the evolution of GC content using changes in body temperature associated with the transition from poikilothermy to homeothermy. Falsification of these models places fundamental constraints upon the plausible pathways for the evolution of isochores. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users. Reviewing Editor: Dr. Nicolas Galtier     

Similar integration but different stability of Alus and LINEs in the human genome

Alus and LINEs (LINE1) are widespread classes of repeats that are very unevenly distributed in the human genome. The majority of GC-poor LINEs reside in the GC-poor isochores whereas GC-rich Alus are mostly present in GC-rich isochores. The discovery that LINES and Alus share similar target site duplication and a common AT-rich insertion site specificity raised the question as to why these two families of repeats show such a different distribution in the genome. This problem was investigated here by studying the isochore distributions of subfamilies of LINES and Alus characterized by different degrees of divergence from the consensus sequences, and of Alus, LINEs and pseudogenes located on chromosomes 21 and 22. Young Alus are more frequent in the GC-poor part of the genome than old Alus. This suggests that the gradual accumulation of Alus in GC-rich isochores has occurred because of their higher stability in compositionally matching chromosomal regions. Densities of Alus and LINEs increase and decrease, respectively, with increasing GC levels, except for the telomeric regions of the analyzed chromosomes. In addition to LINEs, processed pseudogenes are also more frequent in GC-poor isochores. Finally, the present results on Alu and LINE stability/exclusion predict significant losses of Alu DNA from the GC-poor isochores during evolution, a phenomenon apparently due to negative selection against sequences that differ from the isochore composition.     

Compositional heterogeneity within and among isochores in mammalian genomes. II. Some general comments

The presence of long-range correlations and/or mosaicism in DNA sequences results in GC fluctuations, even within individual isochores that are much larger than expected correlation-free 'random' sequences. Neglecting the presence of such fluctuations can lead to incorrect conclusions regarding relative homogeneity or isochore borders. In this commentary, we address these and other methodological issues raised by the variations in GC level within human isochores. We also discuss some recent misconceptions.     

Comparative genomic sequencing reveals a strikingly similar architecture of a conserved syntenic region on human chromosome 11p15.3 (including gene ST5) and mouse chromosome 7

Comparative genomics is a superior way to identify phylogenetically conserved features like genes or regions involved in gene regulation. The comparison of extended orthologous chromosomal regions should also reveal other characteristic traits essential for chromosome or gene function. In the present study we have sequenced and compared a region of conserved synteny from human chromosome 11p15.3 and mouse chromosome 7. In human, this region is known to contain several genes involved in the development of various disorders like Beckwith-Wiedemann overgrowth syndrome and other tumor diseases. Furthermore, in the neighboring chromosome region 11p15.5 extensive imprinting of genes has been reported which might extend to region 11p15.3. The analysis of approximately 730 kb in human and 620 kb in mouse led to the identification of eleven genes. All putative genes found in the mouse DNA were also present in the same order and orientation in the human chromosome. However, in the human DNA one putative gene of unknown function could be identified which is not present in the orthologous position of the mouse chromosome. The sequence similarity between human and mouse is higher in transcribed and exon regions than in non-transcribed segments. Dot plot analysis, however, reveals a surprisingly well-conserved sequence similarity over the entire analyzed region. In particular, the positions of CpG islands, short regions of very high GC content in the 5' region of putative genes, are similar in human and mouse. With respect to base composition, two distinct segments of significantly different GC content exist as well in human as in the mouse. With a GC content of 45% the one segment would correspond to "isochore H1" and the other segment (39% GC in human, 40% GC in mouse) to "isochore L1/L2". The gene density (one gene per 66 kb) is slightly higher than the average calculated for the complete human genome (one gene per 90 kb). The comparison of the number and distribution of repetitive elements shows that the proportion of human DNA made up by interspersed repeats (43.8%) is significantly higher than in the corresponding mouse DNA (30.1%). This partly explains why the human DNA is longer between the landmark genes used to define the orthologous positions in human and mouse.