Isochore structures in the mouse genome

The distribution of the G+C content in the mouse genome has been studied using a windowless technique. We have found that: (i). Abrupt variations of the G+C content from a GC-rich region to a GC-poor region, and vice versa, occur frequently at some sites along the sequence of the mouse genome. (ii). Long domains with relatively homogeneous G+C content (isochores) exist, which usually have sharp boundaries. Consequently, 28 isochores longer than 1 Mb have been identified in the mouse genome. A homogeneity index was used to quantify the variations of the G+C content within isochores. The precise boundaries, sizes, and G+C contents of these isochores have been determined. The windowless technique for the G+C content computation was also used to analyze the DNA sequence containing the mouse MHC region, which has a GC-poor isochore. This isochore is located at the central part of the sequence with boundaries at 468459 and 812716 bp, where the sequence is extended from the centromeric end to the telomeric end. In addition, the analysis of a segment of the rat genome shows that the rat genome also has clear isochore structures.     

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.     

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.     

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)     

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.     

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.     

Compositional mapping of the human dystrophin-encoding gene.

The genomes of warm-blooded vertebrates are mosaics of long DNA segments (> 300 kb, on the average), the isochores, homogeneous in GC levels, which belong to a small number of compositional families. In the present work, the human dystrophin-encoding gene, spanning more than 2.3 Mb in Giemsa band Xp21 (on the short arm of the X chromosome), was analyzed in its isochore organization by hybridizing cDNA probes, corresponding to eight contiguous segments of the coding sequence, on compositional fractions from human DNA. Five DNA regions of uniform (+/- 0.5%) GC content, separated by compositional discontinuities of about 2% GC, were found, so providing the first high-resolution compositional map obtained for a human genome locus and the first direct estimate of isochore size (360 kb to more than 770 kb, in the locus under consideration). One of the isochores contains 71% and another one 21% of deletion breakpoints found in patients suffering from Duchenne's and Becker's muscular dystrophies.     

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

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

Statistical analysis of vertebrate sequences reveals that long genes are scarce in GC-rich isochores

We compared the exon/intron organization of vertebrate genes belonging to different isochore classes, as predicted by their GC content at third codon position. Two main features have emerged from the analysis of sequences published in GenBank: (1) genes coding for long proteins (i.e., 500 aa) are almost two times more frequent in GC-poor than in GC-rich isochores; (2) intervening sequences (=sum of introns) are on average three times longer in GC-poor than in GC-rich isochores. These patterns are observed among human, mouse, rat, cow, and even chicken genes and are therefore likely to be common to all warm-blooded vertebrates. Analysis of Xenopus sequences suggests that the same patterns exist in cold-blooded vertebrates. It could be argued that such results do not reflect the reality because sequence databases are not representative of entire genomes. However, analysis of biases in GenBank revealed that the observed discrepancies between GC-rich and GC-poor isochores are not artifactual, and are probably largely underestimated. We investigated the distribution of microsatellites and interspersed repeats in introns of human and mouse genes from different isochores. This analysis confirmed previous studies showing that Ll repeats are almost absent from GC-rich isochores. Microsatellites and SINES (Alu, B1, B2) are found at roughly equal frequencies in introns from all isochore classes. Globally, the presence of repeated sequences does not account for the increased intron length in GC-poor isochores. The relationships between gene structure and global genome organization and evolution are discussed.     

Evidence of distinct gene functional patterns in GC-poor and GC-rich isochores in Bos taurus

Vertebrate genomes are mosaics of megabase-size DNA segments with a fairly homogeneous base composition, called isochores. They are divided into five families characterized by different guanine-cytosine (GC) levels and linked to several functional and structural properties. The increased availability of fully sequenced genomes allows the investigation of isochores in several species, assessing their level of conservation across vertebrate genomes. In this work, we characterized the isochores in Bos taurus using the ARS-UCD1.2 genome version. The comparison of our results with the well-studied human isochores and those of other mammals revealed a large conservation in isochore families, in number, average GC levels and gene density. Exceptions to the established increase in gene density with the increase in isochores (GC%) were observed for the following gene biotypes: tRNA, small nuclear RNA, small nucleolar RNA and pseudogenes that have their maximum number in H2 and H1 isochores. Subsequently, we assessed the ontology of all gene biotypes looking for functional classes that are statistically over- or under-represented in each isochore. Receptor activity and sensory perception pathways were significantly over-represented in L1 and L2 (GC-poor) isochores. This was also validated for the horse genome. Our analysis of housekeeping genes confirmed a preferential localization in GC-rich isochores, as reported in other species. Finally, we assessed the SNP distribution of a bovine high-density SNP chip across the isochores, finding a higher density in the GC-rich families, reflecting a potential bias in the chip, widely used for genetic selection and biodiversity studies.     

A long-range physical map of human Chromosome 21q22.1 band from the YAC continuum

The human Chromosome (Chr) 21q22.1 region contains several genes for cytokines and neurotransmitters and the gene for superoxide dismutase (mutant forms of which can cause familial amyotrophic lateral sclerosis). A region of approximately 5.8 Mb encompassing D21S82 and the glycinamide ribonucleotide transformylase (GART) loci was covered by overlapping YAC clones, which were contiguously ordered by clone walking with sequence-tagged site (STSs). A total of 76 markers, including 29 YAC end-specific STSs, were unambiguously ordered in this 5.8-Mb region, and the average interval between markers was 76 kb. Restriction maps of the YAC clones with rare-cutting enzymes were simultaneously prepared, and the restriction sites were aligned to obtain a consensus restriction map of the proximal region of the 21q22.1 band. The restriction map made from 44 overlapping YACs contains 54 physically assigned STSs. By integrating the consensus map of the adjacent 1.8-Mb region, we obtained a fine physical map spanning 6.5 Mb of human Chr 21q22.1. This map contains 24 precisely positioned end-specific STSs and 12 NotI-linking markers. More than 39 potential CpG islands were identified in this region and were found to be unevenly distributed. This physical map and the YACs should be useful as a reference map and as a resource for further structural analysis of the Giemsa-negative band (R-band) of Chr 21q22.1. Received: 1 September 1995 / Accepted: 21 November 1995     

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.     

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.     

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.     

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.     

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.      

Molecular definition of a region of chromosome 21 that causes features of the Down syndrome phenotype

Down syndrome (DS) is a major cause of mental retardation and heart disease. Although it is usually caused by the presence of an extra chromosome 21, a subset of the diagnostic features may be caused by the presence of only band 21q22. We now present evidence that significantly narrows the chromosomal region responsible for several of the phenotypic features of DS. We report a molecular and cytogenetic analysis of a three-generation family containing four individuals with clinical DS as manifested by the characteristic facial appearance, endocardial cushion defect, mental retardation, and probably dermatoglyphic changes. Autoradiograms of quantitative Southern blots of DNAs from two affected sisters, their carrier father, and a normal control were analyzed after hybridization with two to six unique DNA sequences regionally mapped on chromosome 21. These include cDNA probes for the genes for CuZn-superoxide dismutase (SOD1) mapping in 21q22.1 and for the amyloid precursor protein (APP) mapping in 21q11.2-21.05, in addition to six probes for single-copy sequences: D21S46 in 21q11.2-21.05, D21S47 and SF57 in 21q22.1-22.3, and D21S39, D21S42, and D21S43 in 21q22.3. All sequences located in 21q22.3 were present in three copies in the affected individuals, whereas those located proximal to this region were present in only two copies. In the carrier father, all DNA sequences were present in only two copies. Cytogenetic analysis of affected individuals employing R and G banding of prometaphase preparations combined with in situ hybridization revealed a translocation of the region from very distal 21q22.1 to 21qter to chromosome 4q. Except for a possible phenotypic contribution from the deletion of chromosome band 4q35, these data provide a molecular definition of the minimal region of chromosome 21 which, when duplicated, generates the facial features, heart defect, a component of the mental retardation, and probably several of the dermatoglyphic changes of DS. This region may include parts of bands 21q22.2 and 21q22.3, but it must exclude the genes S0D1 and APP and most of band 21q22.1, specifically the region defined by S0D1, SF57 and D21S47.     

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.