Correlations between isochores and chromosomal bands in the human genome

Chromosomal Rearrangements and Altered Nuclear Organization: Recent Mechanistic Models in Cancer

Simple Summary

New methodologies and technologies developed in the last few decades have highlighted the precise spatial organization of the genome into the cell nucleus, with chromatin architecture playing a central role in controlling several genome functions. Genes are expressed in a well-defined way and at a well-defined time during cell differentiation, and alterations in genome organization can lead to genetic diseases, such as cancers. Here we review how the genome is organized in the cell nucleus and the evidence of genome misorganization leading to cancer diseases.

Abstract

The last decade has seen significant progress in understanding how the genome is organized spatially within interphase nuclei. Recent analyses have confirmed earlier molecular cytogenetic studies on chromosome positioning within interphase nuclei and provided new information about the topologically associated domains (TADs). Examining the nuances of how genomes are organized within interphase nuclei will provide information fundamental to understanding gene regulation and expression in health and disease. Indeed, the radial spatial positioning of individual gene loci within nuclei has been associated with up- and down-regulation of specific genes, and disruption of normal genome organization within nuclei will result in compromised cellular health. In cancer cells, where reorganization of the nuclear architecture may occur in the presence of chromosomal rearrangements such as translocations, inversions, or deletions, gene repositioning can change their expression. To date, very few studies have focused on radial gene positioning and the correlation to gene expression in cancers. Further investigations would improve our understanding of the biological mechanisms at the basis of cancer and, in particular, in leukemia initiation and progression, especially in those cases where the molecular consequences of chromosomal rearrangements are still unclear. In this review, we summarize the main milestones in the field of genome organization in the nucleus and the alterations to this organization that can lead to cancer diseases.

The "Genomic Code": DNA Pervasively Moulds Chromatin Structures Leaving no Room for "Junk"

The chromatin of the human genome was analyzed at three DNA size levels. At the first, compartment level, two "gene spaces" were found many years ago: A GC-rich, gene-rich "genome core" and a GC-poor, gene-poor "genome desert", the former corresponding to open chromatin centrally located in the interphase nucleus, the latter to closed chromatin located peripherally. This bimodality was later confirmed and extended by the discoveries (1) of LADs, the Lamina-Associated Domains, and InterLADs; (2) of two "spatial compartments", A and B, identified on the basis of chromatin interactions; and (3) of "forests and prairies" characterized by high and low CpG islands densities. Chromatin compartments were shown to be associated with the compositionally different, flat and single- or multi-peak DNA structures of the two, GC-poor and GC-rich, "super-families" of isochores. At the second, sub-compartment, level, chromatin corresponds to flat isochores and to isochore loops (due to compositional DNA gradients) that are susceptible to extrusion. Finally, at the short-sequence level, two sets of sequences, GC-poor and GC-rich, define two different nucleosome spacings, a short one and a long one. In conclusion, chromatin structures are moulded according to a "genomic code" by DNA sequences that pervade the genome and leave no room for "junk".

The Genomic Code: A Pervasive Encoding/Molding of Chromatin Structures and a Solution of the "Non-Coding DNA" Mystery

Recent investigations have revealed 1) that the isochores of the human genome group into two super-families characterized by two different long-range 3D structures, and 2) that these structures, essentially based on the distribution and topology of short sequences, mold primary chromatin domains (and define nucleosome binding). More specifically, GC-poor, gene-poor isochores are low-heterogeneity sequences with oligo-A spikes that mold the lamina-associated domains (LADs), whereas GC-rich, gene-rich isochores are characterized by single or multiple GC peaks that mold the topologically associating domains (TADs). The formation of these "primary TADs" may be followed by extrusion under the action of cohesin and CTCF. Finally, the genomic code, which is responsible for the pervasive encoding and molding of primary chromatin domains (LADs and primary TADs, namely the "gene spaces"/"spatial compartments") resolves the longstanding problems of "non-coding DNA," "junk DNA," and "selfish DNA" leading to a new vision of the genome as shaped by DNA sequences.

R/G-band boundaries: genomic instability and human disease

The human genome is composed of large-scale compartmentalized structures resulting from variations in the amount of guanine and cytosine residues (GC%) and in the timing of DNA replication. These compartmentalized structures are related to the light- and dark-staining bands along chromosomes after the appropriate staining. Here we describe our current understanding of the biological importance of the boundaries between these light and dark bands (the so-called R/G boundaries). These R/G boundaries were identified following integration of information obtained from analyses of chromosome bands and genome sequences. This review also discusses the potential medical significance of these chromosomal regions for conditions related to genomic instability, such as cancer and neural disease. We propose that R/G-chromosomal boundaries, which correspond to regions showing a switch in replication timing from early to late S phase (early/late-switch regions) and of transition in GC%, have an extremely low number of replication origins and more non-B-form DNA structures than other genomic regions. Further, we suggest that genes located at R/G boundaries and which contain such DNA sequences have an increased risk of genetic instability and of being associated with human diseases. Finally, we propose strategies for genome and epigenome analyses based on R/G boundaries.

Intra-genomic GC heterogeneity in sauropsids: evolutionary insights from cDNA mapping and GC(3) profiling in snake

BACKGROUND

Extant sauropsids (reptiles and birds) are divided into two major lineages, the lineage of Testudines (turtles) and Archosauria (crocodilians and birds) and the lineage of Lepidosauria (tuatara, lizards, worm lizards and snakes). Karyotypes of these sauropsidan groups generally consist of macrochromosomes and microchromosomes. In chicken, microchromosomes exhibit a higher GC-content than macrochromosomes. To examine the pattern of intra-genomic GC heterogeneity in lepidosaurian genomes, we constructed a cytogenetic map of the Japanese four-striped rat snake (Elaphe quadrivirgata) with 183 cDNA clones by fluorescence in situ hybridization, and examined the correlation between the GC-content of exonic third codon positions (GC3) of the genes and the size of chromosomes on which the genes were localized.

RESULTS

Although GC3 distribution of snake genes was relatively homogeneous compared with those of the other amniotes, microchromosomal genes showed significantly higher GC3 than macrochromosomal genes as in chicken. Our snake cytogenetic map also identified several conserved segments between the snake macrochromosomes and the chicken microchromosomes. Cross-species comparisons revealed that GC3 of most snake orthologs in such macrochromosomal segments were GC-poor (GC3 < 50%) whereas those of chicken orthologs in microchromosomes were relatively GC-rich (GC3 ≥ 50%).

CONCLUSION

Our results suggest that the chromosome size-dependent GC heterogeneity had already occurred before the lepidosaur-archosaur split, 275 million years ago. This character was probably present in the common ancestor of lepidosaurs and but lost in the lineage leading to Anolis during the diversification of lepidosaurs. We also identified several genes whose GC-content might have been influenced by the size of the chromosomes on which they were harbored over the course of sauropsid evolution.

Spatial analyses of mono, di and trinucleotide trends in plant genes

Genomic DNA sequences display compositional heterogeneity on many scales. In this paper we analyzed tendencies and anomalies in the occurence of mono, di and trinucleotides in structural regions of plant genes. Representation of these trends as a function of position along genic sequences highlighted compositional features peculiar of either monocots or eudicots that were remarkably uniform within these two evolutionary clades. The most evident of these features appeared in the form of gradient of base content along the direction of transcription. The robustness of such a representation was validated in sequences sub-datasets generated considering structural and compositional features such as total length of cds, overall GC content and genic orientation in the genome. Piecewise regression analyses indicated that the gradients could be conveniently approximated to a two segmented model where a first region featuring a steep slope is followed by a second segment fitting a milder variation. In general, monocots species showed steeper segments than eudicots. The guanine gradient was the most distinctive feature between the two evolutionary clades, being moderately increasing in eudicots and firmly decreasing in monocots. Single gene investigation revealed that a high proportion of genes show compositional trends compatible with a segmented model suggesting that these features are essential attributes of gene organization. Dinucleotide and trinucleotide biases were referred to expectation based on a random union of the component elements. The average bias at dinucleotide level identified a significant undererpresentation of some dinucleotide and the overrepresention of others. The bias at trinucleotide level was on average low. Finally, the analysis of bryophyte coding sequences showed mononucleotide, dinucleotide and trinucleotide compositional trends resembling those of higher plants. This finding suggested that the emergenge of compositional bias is an ancient event in evolution which was already present at the time of land conquest by green plants.

Disorders caused by chromosome abnormalities

Many human genetic disorders result from unbalanced chromosome abnormalities, in which there is a net gain or loss of genetic material. Such imbalances often disrupt large numbers of dosage-sensitive, developmentally important genes and result in specific and complex phenotypes. Alternately, some chromosomal syndromes may be caused by a deletion or duplication of a single gene with pleiotropic effects. Traditionally, chromosome abnormalities were identified by visual inspection of the chromosomes under a microscope. The use of molecular cytogenetic technologies, such as fluorescence in situ hybridization and microarrays, has allowed for the identification of cryptic or submicroscopic imbalances, which are not visible under the light microscope. Microarrays have allowed for the identification of numerous new syndromes through a genotype-first approach in which patients with the same or overlapping genomic alterations are identified and then the phenotypes are described. Because many chromosomal alterations are large and encompass numerous genes, the ascertainment of individuals with overlapping deletions and varying clinical features may allow researchers to narrow the region in which to search for candidate genes.

Capture of AT-rich chromatin by ELYS recruits POM121 and NDC1 to initiate nuclear pore assembly

Assembly of the nuclear pore, gateway to the genome, from its component subunits is a complex process. In higher eukaryotes, nuclear pore assembly begins with the binding of ELYS/MEL-28 to chromatin and recruitment of the large critical Nup107-160 pore subunit. The choreography of steps that follow is largely speculative. Here, we set out to molecularly define early steps in nuclear pore assembly, beginning with chromatin binding. Point mutation analysis indicates that pore assembly is exquisitely sensitive to the change of only two amino acids in the AT-hook motif of ELYS. The dependence on AT-rich chromatin for ELYS binding is borne out by the use of two DNA-binding antibiotics. AT-binding Distamycin A largely blocks nuclear pore assembly, whereas GC-binding Chromomycin A(3) does not. Next, we find that recruitment of vesicles containing the key integral membrane pore proteins POM121 and NDC1 to the forming nucleus is dependent on chromatin-bound ELYS/Nup107-160 complex, whereas recruitment of gp210 vesicles is not. Indeed, we reveal an interaction between the cytoplasmic domain of POM121 and the Nup107-160 complex. Our data thus suggest an order for nuclear pore assembly of 1) AT-rich chromatin sites, 2) ELYS, 3) the Nup107-160 complex, and 4) POM121- and NDC1-containing membrane vesicles and/or sheets, followed by (5) assembly of the bulk of the remaining soluble pore subunits.

Assignment of isochores for all completely sequenced vertebrate genomes using a consensus

A new consensus isochore assignment method and a database of isochore maps for all completely sequenced vertebrate genomes are presented.

We show that although the currently available isochore mapping methods agree on the isochore classification of about two-thirds of the human DNA, they produce significantly different results with regard to the location of isochore boundaries and isochore length distribution. We present a new consensus isochore assignment method based on majority voting and provide IsoBase, a comprehensive on-line database of isochore maps for all completely sequenced vertebrate genomes.

Prediction of replication time zones at single nucleotide resolution in the human genome

The human genome is structured at multiple levels: it is organized into a series of replication time zones, and meanwhile it is composed of isochores. Accumulating evidence suggests a match between these two genome features. Based on newly developed software GC-Profile, we obtained a complete coverage of the human genome by 3198 isochores with boundaries at single nucleotide resolution. Interestingly, the experimentally confirmed replication timing sites in the regions of 1p36.1, 6p21.32, 17q11.2 and 22q12.1 nearly all coincide with the determined isochore boundaries. The precise boundaries of the 3198 isochores are available via the website: http://tubic.tju.edu.cn/isomap/.

Reverse polarization in amino acid and nucleotide substitution patterns between human-mouse orthologs of two compositional extrema

Genome-wide analysis of sequence divergence patterns in 12,024 human-mouse orthologous pairs reveals, for the first time, that the trends in nucleotide and amino acid substitutions in orthologs of high and low GC composition are highly asymmetric and polarized to opposite directions. The entire dataset has been divided into three groups on the basis of the GC content at third codon sites of human genes: high, medium, and low. High-GC orthologs exhibit significant bias in favor of the replacements, Thr --> Ala, Ser --> Ala, Val --> Ala, Lys --> Arg, Asn --> Ser, Ile --> Val etc., from mouse to human, whereas in low-GC orthologs, the reverse trends prevail. In general, in the high-GC group, residues encoded by A/U-rich codons of mouse proteins tend to be replaced by the residues encoded by relatively G/C-rich codons in their human orthologs, whereas the opposite trend is observed among the low-GC orthologous pairs. The medium-GC group shares some trends with high-GC group and some with low-GC group. The only significant trend common in all groups of orthologs, irrespective of their GC bias, is (Asp)(Mouse) --> (Glu)(Human) replacement. At the nucleotide level, high-GC orthologs have undergone a large excess of (A/T)(Mouse) --> (G/C)(Human) substitutions over (G/C)(Mouse) --> (A/T)(Human) at each codon position, whereas for low-GC orthologs, the reverse is true.

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.

Pan-S replication patterns and chromosomal domains defined by genome-tiling arrays of ENCODE genomic areas

In eukaryotes, accurate control of replication time is required for the efficient completion of S phase and maintenance of genome stability. We present a high-resolution genome-tiling array-based profile of replication timing for approximately 1% of the human genome studied by The ENCODE Project Consortium. Twenty percent of the investigated segments replicate asynchronously (pan-S). These areas are rich in genes and CpG islands, features they share with early-replicating loci. Interphase FISH showed that pan-S replication is a consequence of interallelic variation in replication time and is not an artifact derived from a specific cell cycle synchronization method or from aneuploidy. The interallelic variation in replication time is likely due to interallelic variation in chromatin environment, because while the early- or late-replicating areas were exclusively enriched in activating or repressing histone modifications, respectively, the pan-S areas had both types of histone modification. The replication profile of the chromosomes identified contiguous chromosomal segments of hundreds of kilobases separated by smaller segments where the replication time underwent an acute transition. Close examination of one such segment demonstrated that the delay of replication time was accompanied by a decrease in level of gene expression and appearance of repressive chromatin marks, suggesting that the transition segments are boundary elements separating chromosomal domains with different chromatin environments.

A new platform linking chromosomal and sequence information

We have tested whether a direct correlation of sequence information and staining properties of chromosomes is possible and whether this combined information can be used to precisely map any position on the chromosome. Despite huge differences of compaction between the naked DNA and the DNA packed in chromosomes we found a striking correlation when visualizing the GGCC density on both levels. Software was developed that allows one to superimpose chromosomal fluorescence intensity profiles generated by chromolysin A3 (CMA3) staining with GGCC density extracted from the Ensembl database. Thus, any position along the chromosome can be defined in megabase pairs (Mb) besides the cytoband information, enabling direct alignment of chromosomal information with the sequence data. The mapping tool was validated using 13 different BAC clones, resulting in a mean difference from Ensembl data of 2 Mb (ranging from 0.79 to 3.57 Mb). Our results indicate that the sequence density information and information gained with sequence-specific fluorochromes are superimposable. Thus, the visualized GGCC motif density along the chromosome (sequence bands) provides a unique platform for comparing different types of genomic information.

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.

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.

The rate, not the spectrum, of base pair substitutions changes at a GC-content transition in the human NF1 gene region: implications for the evolution of the mammalian genome structure

The human genome is composed of long stretches of DNA with distinct GC contents, called isochores or GC-content domains. A boundary between two GC-content domains in the human NF1 gene region is also a boundary between domains of early- and late-replicating sequences and of regions with high and low recombination frequencies. The perfect conservation of the GC-content distribution in this region between human and mouse demonstrates that GC-content stabilizing forces must act regionally on a fine scale at this locus. To further elucidate the nature of these forces, we report here on the spectrum of human SNPs and base pair substitutions between human and chimpanzee. The results show that the mutation rate changes exactly at the GC-content transition zone from low values in the GC-poor sequences to high values in GC-rich ones. The GC content of the GC-poor sequences can be explained by a bias in favor of GC > AT mutations, whereas the GC content of the GC-rich segment may result from a fixation bias in favor of AT > GC substitutions. This fixation bias may be explained by direct selection by the GC content or by biased gene conversion.

Chromosome organization and chromatin modification: influence on genome function and evolution

Histone modifications of nucleosomes distinguish euchromatic from heterochromatic chromatin states, distinguish gene regulation in eukaryotes from that of prokaryotes, and appear to allow eukaryotes to focus recombination events on regions of highest gene concentrations. Four additional epigenetic mechanisms that regulate commitment of cell lineages to their differentiated states are involved in the inheritance of differentiated states, e.g., DNA methylation, RNA interference, gene repositioning between interphase compartments, and gene replication time. The number of additional mechanisms used increases with the taxon's somatic complexity. The ability of siRNA transcribed from one locus to target, in trans, RNAi-associated nucleation of heterochromatin in distal, but complementary, loci seems central to orchestration of chromatin states along chromosomes. Most genes are inactive when heterochromatic. However, genes within beta-heterochromatin actually require the heterochromatic state for their activity, a property that uniquely positions such genes as sources of siRNA to target heterochromatinization of both the source locus and distal loci. Vertebrate chromosomes are organized into permanent structures that, during S-phase, regulate simultaneous firing of replicon clusters. The late replicating clusters, seen as G-bands during metaphase and as meiotic chromomeres during meiosis, epitomize an ontological utilization of all five self-reinforcing epigenetic mechanisms to regulate the reversible chromatin state called facultative (conditional) heterochromatin. Alternating euchromatin/heterochromatin domains separated by band boundaries, and interphase repositioning of G-band genes during ontological commitment can impose constraints on both meiotic interactions and mammalian karyotype evolution.

Isochore structures in the chicken genome

The availability of the complete chicken genome sequence provides an unprecedented opportunity to study the global genome organization at the sequence level. Delineating compositionally homogeneous G + C domains in DNA sequences can provide much insight into the understanding of the organization and biological functions of the chicken genome. A new segmentation algorithm, which is simple and fast, has been proposed to partition a given genome or DNA sequence into compositionally distinct domains. By applying the new segmentation algorithm to the draft chicken genome sequence, the mosaic organization of the chicken genome can be confirmed at the sequence level. It is shown herein that the chicken genome is also characterized by a mosaic structure of isochores, long DNA segments that are fairly homogeneous in the G + C content. Consequently, 25 isochores longer than 2 Mb (megabases) have been identified in the chicken genome. These isochores have a fairly homogeneous G + C content and often correspond to meaningful biological units. With the aid of the technique of cumulative GC profile, we proposed an intuitive picture to display the distribution of segmentation points. The relationships between G + C content and the distributions of genes (CpG islands, and other genomic elements) were analyzed in a perceivable manner. The cumulative GC profile, equipped with the new segmentation algorithm, would be an appropriate starting point for analyzing the isochore structures of higher eukaryotic genomes.

The temporal program of DNA replication: new insights into old questions

During the last decades it has been shown that the replication timing program in metazoans is related to chromosome structure, the nuclear positioning and AT/GC content of chromosomal loci, their patterns of histone modifications, and their transcriptional regulation. Here, the current state of knowledge concerning these relationships is reviewed. An integrated view on structure-function relationships in the nucleus is provided and the determination and functional role of the replication timing program is discussed in this context. A corresponding comprehensive model is developed and a key aspect of this model is the suggestion that mammalian chromosomes are organized into stable units equivalent to replicon clusters. It is proposed that the nuclear positions of these units would depend on their histone modifications and determine the replication timing of the whole unit. It is furthermore predicted that replication timing is only indirectly linked to transcriptional regulation and contributes to the maintenance of gene expression patterns. These clear predictions, and the fact that the tools are at hand now to further test them, open an avenue towards solving the long standing problem on how replication timing is determined in metazoan cells.

Strong Regional Biases in Nucleotide Substitution in the Chicken Genome

Interspersed repeats have emerged as a valuable tool for studying neutral patterns of molecular evolution. Here we analyze variation in the rate and pattern of nucleotide substitution across all autosomes in the chicken genome by comparing the present-day CR1 repeat sequences with their ancestral copies and reconstructing nucleotide substitutions with a maximum likelihood model. The results shed light on the origin and evolution of large-scale heterogeneity in GC content found in the genomes of birds and mammals--the isochore structure. In contrast to mammals, where GC content is becoming homogenized, heterogeneity in GC content is being reinforced in the chicken genome. This is also supported by patterns of substitution inferred from alignments of introns in chicken, turkey, and quail. Analysis of individual substitution frequencies is consistent with the biased gene conversion (BGC) model of isochore evolution, and it is likely that patterns of evolution in the chicken genome closely resemble those in the ancestral amniote genome, when it is inferred that isochores originated. Microchromosomes and distal regions of macrochromosomes are found to have elevated substitution rates and a more GC-biased pattern of nucleotide substitution. This can largely be accounted for by a strong correlation between GC content and the rate and pattern of substitution. The results suggest that an interaction between increased mutability at CpG motifs and fixation biases due to BGC could explain increased levels of divergence in GC-rich regions.

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.

Chimpanzee chromosomes: retrotransposable compound repeat DNA organization (RCRO) and its influence on meiotic prophase and crossing-over

The terminal C-bands that are a specific feature of chimpanzee chromosomes were dissected using a molecular cytogenetic technique, PRINS, with primers for telomeric sequences, subterminal satellite, and retrotransposable elements (HERV-K and -W). These DNA elements jointly formed a large block of retrotransposable compound repeat DNA organization (RCRO) at the terminal C-band regions of 30 chromosomes, and are also located at the centromeric regions of some chromosomes. Additionally, a block consisting of all members of the RCRO has transposed to the middle (q31.1) of the long arm of chromosome 6, and three members, the subterminal satellite and the two HERVs, have integrated into the proximal region (q14.4) of the long arm of chromosome 14. Terminal RCROs seem to induce and prolong the bouquet stage in meiotic prophase, and to affect chiasma formation, together with interstitial RCROs. It is also postulated that RCROs may cause a position effect to gene expression, resulting in gene silencing and/or late replication.

Substantial regional variation in substitution rates in the human genome: importance of GC content, gene density, and telomere-specific effects

This study presents the first global, 1-Mbp-level analysis of patterns of nucleotide substitutions along the human lineage. The study is based on the analysis of a large amount of repetitive elements deposited into the human genome since the mammalian radiation, yielding a number of results that would have been difficult to obtain using the more conventional comparative method of analysis. This analysis revealed substantial and consistent variability of rates of substitution, with the variability ranging up to twofold among different regions. The rates of substitutions of C or G nucleotides with A or T nucleotides vary much more sharply than the reverse rates, suggesting that much of that variation is due to differences in mutation rates rather than in the probabilities of fixation of C/G vs. A/T nucleotides across the genome. For all types of substitution we observe substantially more hotspots than coldspots, with hotspots showing substantial clustering over tens of Mbp's. Our analysis revealed that GC-content of surrounding sequences is the best predictor of the rates of substitution. The pattern of substitution appears very different near telomeres compared to the rest of the genome and cannot be explained by the genome-wide correlations of the substitution rates with GC content or exon density. The telomere pattern of substitution is consistent with natural selection or biased gene conversion acting to increase the GC-content of the sequences that are within 10-15 Mbp away from the telomere.

Male-driven biased gene conversion governs the evolution of base composition in human alu repeats

Regional biases in substitution pattern are likely to be responsible for the large-scale variation in base composition observed in vertebrate genomes. However, the evolutionary forces responsible for these biases are still not clearly defined. In order to study the processes of mutation and fixation across the entire human genome, we analyzed patterns of substitution in Alu repeats since their insertion. We also studied patterns of human polymorphism within the repeats. There is a highly significant effect of recombination rate on the pattern of substitution, whereas no such effect is seen on the pattern of polymorphism. These results suggest that regional biases in substitution are caused by biased gene conversion, a process that increases the probability of fixation of mutations that increase GC content. Furthermore, the strongest correlate of substitution patterns is found to be male recombination rates rather than female or sex-averaged recombination rates. This indicates that in addition to sexual dimorphism in recombination rates, the sexes also differ in the relative rates of crossover and gene conversion.

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.

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.

Chromatin architecture of the human genome: gene-rich domains are enriched in open chromatin fibers

We present an analysis of chromatin fiber structure across the human genome. Compact and open chromatin fiber structures were separated by sucrose sedimentation and their distributions analyzed by hybridization to metaphase chromosomes and genomic microarrays. We show that compact chromatin fibers originate from some sites of heterochromatin (C-bands), and G-bands (euchromatin). Open chromatin fibers correlate with regions of highest gene density, but not with gene expression since inactive genes can be in domains of open chromatin, and active genes in regions of low gene density can be embedded in compact chromatin fibers. Moreover, we show that chromatin fiber structure impacts on further levels of chromatin condensation. Regions of open chromatin fibers are cytologically decondensed and have a distinctive nuclear organization. We suggest that domains of open chromatin may create an environment that facilitates transcriptional activation and could provide an evolutionary constraint to maintain clusters of genes together along chromosomes.

Representing GC variation along eukaryotic chromosomes

Genome sequencing now permits direct visual representation, at any scale, of GC heterogeneity along the chromosomes of several higher eukaryotes. Plots can be easily obtained from the chromosomal sequences, yet sequence releases of mammalian or plant chromosomes still tend to use small scales or window sizes that obscure important large-scale compositional features. To faithfully reveal, at one glance, the compositional variation at a given scale, we have devised a simple scheme that combines line plots with color-coded shading of the regions underneath the plots. The scheme can be applied to different eukaryotic genomes to facilitate their comparison, as illustrated here for a sample of chromosomes chosen from seven selected species. As a complement to a previously published compact view of isochores in the human genome sequence, we include here an analogous map for the recently sequenced mouse genome, and discuss the contribution of repetitive DNA to the GC variation along the plots. Supplementary information, including a database of color-coded GC profiles for all recently sequenced eukaryotes and the program draw_chromosomes_gc.pl used to obtain them, are available at.

Isochore Structures in the Genome of the Plant Arabidopsis thaliana

Arabidopsis thaliana is an important model system for the study of plant biology. We have analyzed the complete genome sequences of Arabidopsis by using a newly developed windowless method for the GC content computation, the cumulative GC profile. It is shown that the Arabidopsis genome is organized into a mosaic structure of isochores. All the centromeric regions are located in GC-rich isochores, called centromere-isochores, which are characterized by a high GC content but low gene and T-DNA insertion densities. This characteristic distinguishes centromere-isochores from the other class of GC-rich isochores, called GC-isochores, which have high gene and T-DNA insertion densities. Consequently, 15 isochores have been identified, i.e., 7 AT-isochores, 3 GC-isochores, and 5 centromere-isochores. The genes in centromere-isochores, which have the highest GC content, have much shorter intron lengths and lower intron numbers, compared to those of the other two types. There is also considerable difference in the numbers and lengths of transposable elements (TEs) between AT and GC-isochores, i.e., the TE number (length) of AT-isochores is 6.3 (7.3) times that of GC-isochores. It is generally believed that TEs are accumulated in the regions surrounding the centromeres. However, within these TE-rich regions, there are regions of extremely low TE numbers (TE deserts), which correspond to the positions of centromere-isochores. In addition, a heterochromatic knob is located at the boundary of an AT-isochore. Furthermore, we show that the differences in GC content among isochores are mainly due to the GC content variation of introns, the third codon positions and intergenic regions.

The DNA sequence and analysis of human chromosome 13

Chromosome 13 is the largest acrocentric human chromosome. It carries genes involved in cancer including the breast cancer type 2 (BRCA2) and retinoblastoma (RB1) genes, is frequently rearranged in B-cell chronic lymphocytic leukaemia, and contains the DAOA locus associated with bipolar disorder and schizophrenia. We describe completion and analysis of 95.5 megabases (Mb) of sequence from chromosome 13, which contains 633 genes and 296 pseudogenes. We estimate that more than 95.4% of the protein-coding genes of this chromosome have been identified, on the basis of comparison with other vertebrate genome sequences. Additionally, 105 putative non-coding RNA genes were found. Chromosome 13 has one of the lowest gene densities (6.5 genes per Mb) among human chromosomes, and contains a central region of 38 Mb where the gene density drops to only 3.1 genes per Mb.

Mitochondrial DNA-like sequences in the nucleus (NUMTs): insights into our African origins and the mechanism of foreign DNA integration

Nuclear mitochondrial DNA sequences (NUMTs) are common in eukaryotes. However, the mechanism by which they integrate into the nuclear genome remains a riddle. We analyzed 247 NUMTs in the human nuclear DNA (nDNA), along with their flanking regions. This analysis revealed that some NUMTs have accumulated many changes, and thus have resided in the nucleus a long time, while others are >94% similar to the reference human mitochondrial DNA (mtDNA), and thus must be recent. Among the latter, two NUMTs, encompassing the COI gene, carry a set of transitions characteristic of the extant African-specific L macrohaplogroup mtDNAs and are more homologous to human mtDNA than to chimp. Screening for one of these NUMTs revealed its presence in all human samples tested, confirming that the African macrohaplogroup L mtDNAs were present in the earliest modern humans and thus were the first human mtDNAs. An analysis of flanking sequences of the NUMTs revealed that 59% were within 150 bp of repetitive elements, with 26% being within 15 bp of and 33% being within 15-150 bp of repetitive elements. Only 14% were integrated into a repetitive element. This association of NUMTs with repetitive elements is highly nonrandom (p<0.001). These data suggest that the vicinity of transposable elements influences the ongoing integration of mtDNA sequences and their subsequent duplication within the nDNA. Finally, NUMTs appear to preferentially integrate into DNA with different GC content than the surrounding chromosomal band. Our results suggest that chromosomal structure might influence integration of NUMTs.

Mapping and initial analysis of human subtelomeric sequence assemblies

Physical mapping data were combined with public draft and finished sequences to derive subtelomeric sequence assemblies for each of the 41 genetically distinct human telomere regions. Sequence gaps that remain on the reference telomeres are generally small,well-defined,and for the most part,restricted to regions directly adjacent to the terminal (TTAGGG)n tract. Of the 20.66 Mb of subtelomeric DNA analyzed, 3.01 Mb are subtelomeric repeat sequences (Srpt),and an additional 2.11 Mb are segmental duplications. The subtelomeric sequence assemblies are enriched >25-fold in short,internal (TTAGGG)n-like sequences relative to the rest of the genome; a total of 114 (TTAGGG)n-like islands were found,55 within Srpt regions,35 within one-copy regions,11 at one-copy/Srpt or Srpt/segmental duplication boundaries,and 13 at the telomeric ends of assemblies. Transcripts were annotated in each assembly,noting their mapping coordinates relative to their respective telomere and whether they originate in duplicated DNA or single-copy DNA. A total of 697 transcripts were found in 15.53 Mb of one-copy DNA,76 transcripts in 2.11 Mb of segmentally duplicated DNA,and 168 transcripts in 3.01 Mb of Srpt sequence. This overall transcript density is similar (within approximately 10%) to that found genome-wide. Zinc finger-containing genes and olfactory receptor genes are duplicated within and between multiple telomere regions.

Identification of isochore boundaries in the human genome using the technique of wavelet multiresolution analysis

Incorporated with the Z curve method, the technique of wavelet multiresolution (also known as multiscale) analysis has been proposed to identify the boundaries of isochores in the human genome. The human MHC sequence and the longest contigs of human chromosomes 21 and 22 are used as examples. The boundary between the isochores of Class III and Class II in the MHC sequence has been detected and found to be situated at the position 2,490,368bp. This result is in good agreement with the experimental evidence. An isochore with a length of about 7Mb in chromosome 21 has been identified and found to be gene- and Alu-poor. We have also found that the G+C content of chromosome 21 is more homogeneous than that of chromosome 22. Compared with the window-based methods, the present method has the highest resolution for identifying the boundaries of isochores, even at a scale of single base. Compared with the entropic segmentation method, the present method has the merits of more intuitiveness and less calculations. The important conclusion drawn in this study is that the segmentation points, at which the G+C content undergoes relatively dramatic changes, do exist in the human genome. These 'singularity' points may be considered to be candidates of isochore boundaries in the human genome. The method presented is a general one and can be used to analyze any other genomes.

An isochore map of the human genome based on the Z curve method

The distribution of the G+C content in the human genome has been studied by using a windowless technique derived from the Z curve method. The most important findings presented in this paper are twofold. First, abrupt variations of the G+C content along human chromosome sequences are the main variation patterns of G+C content. It is found that at some sites, the G+C content undergoes abrupt changes from a G+C-rich region to a G+C-poor region alternatively and vice versa. Second, it is shown that long domains with relatively homogeneous G+C content along each chromosome do exist. These domains are thought to be isochores, which usually have sharp boundaries. Consequently, 56 isochores longer than 3 Mb have been identified in chromosomes 1-22, X and Y. Boundaries, size and G+C content of each isochore identified are listed in detail. As an example to demonstrate the power of the method, the boundary between the Classes III and II isochores of the MHC sequence has been determined and found to be at 2,477,936, which is in good agreement with the experimental evidence. A homogeneity index is introduced to measure the homogeneity of G+C content in isochores. We emphasize that the homogeneity of G+C content is relative. The isochores in which the G+C content keeps absolutely constant do not exist. Isochore structures appear to be a basic organization of the human genome. Due to the relevance to many important biological functions, the clarification of isochore structures will provide much insight into the understanding of the human genome.

A unification of mosaic structures in the human genome

The human genome is a mosaic structure on many levels: there exist cytogenetic bands, GC composition bands (isochores) and clusters of broadly expressed genes. How might these inter-relate? It has been proposed that to optimize gene regulation, housekeeping genes should concentrate on transcriptionally competent chromosomal domains. Prior evidence suggests that regions of high GC and R bands are associated with such domains. Here we report that broadly expressed genes cluster in regions of high GC, and in R and lightest Giemsa bands. This is not only a confirmation of the adaptive hypothesis, but is also the first direct systematic evidence of a general interdependence of expression patterns with base composition and chromosome structure.

Structure of human chromosomes studied by atomic force microscopy. Part II. Relationship between structure and cytogenetic bands

In the first part of this work, human chromosomes were characterized by atomic force microscopy (AFM) in air and in aqueous solution. The analysis of the images suggests that the last level of organization consists of a radial arrangement of chromatin loops which are anchored to a fiber which is folded giving a pattern of bands which differs in volume. Here the pattern of bands observed by AFM is compared to the cytogenetic map at the 850-band level. Thus thicker and thinner bands are identified as G and R bands, respectively. Finally a model is proposed which links genome sequence, cytogenetics, and chromosome structure.

Expression patterns and gene distribution in the human genome

Genes are non-uniformly distributed in the human genome, reaching the highest concentration in GC-rich isochores. This is one of the fundamental aspects of the human genome organization (Gene 241/259 (2000a,b) 3/31, for a review). In the present paper the gene distribution was analyzed in relationship to the gene expression pattern and levels. In this study evidence is produced showing: (i) that a biased gene distribution towards GC-rich isochores applies to both tissue-specific and housekeeping genes; and (ii) that genes localized in GC-rich isochores have high transcriptional levels. Since gene density and transcriptional levels are correlated with each other and both are correlated with the GC content of the isochores, the biased gene distribution in the human genome presumably is the result of selection at the gene expression levels.

Isochore chromosome maps of the human genome

The human genome is a mosaic of isochores, which are long DNA segments (z.Gt;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.

Compositional mapping of chicken chromosomes and identification of the gene-richest regions

'Compositional chromosomal mapping', namely the assessment of the GC level of chromosomal bands, led to the identification, in the human chromosomes, of the GC-richest H3+ bands and of the GC-poorest L1+ bands, which were so called on the basis of the isochore family predominantly present in the bands. The isochore organization of the avian genome is very similar to those of most mammals, the only difference being the presence of an additional, GC-richest, H4 isochore family. In contrast, the avian karyotypes are very different from those of mammals, being characterized, in most species, by few macrochromosomes and by a large number of microchromosomes. The 'compositional mapping' of chicken mitotic and meiotic chromosomes by in-situ hybridization of isochore families showed that the chicken GC-richest isochores are localized not only on a large number of microchromosomes but also on almost all telomeric bands of macrochromosomes. On the other hand, the GC-poorest isochores are generally localized on the internal regions of macrochromosomes and are almost absent in microchromosomes. Thus, the distinct localization of the GC-richest and the GC-poorest bands observed on human chromosomes appears to be a general feature of chromosomes from warm-blooded vertebrates.

Genes, isochores and bands in human chromosomes 21 and 22

The recently available DNA sequences from chromosomes 21 and 22 enabled us to define the relationships of different band types with isochores and with gene concentration and to compare these relationships with previous results. We showed that chromosomal bands appear as Giemsa or Reverse bands depending not on their absolute GC level, but on the composition GC level relative to those of adjacent contiguous bands. We also demonstrated that the GC-richest, and gene-richest H3+ bands are characterized by a lower DNA compaction compared with the GC-poorest, gene-poorest L1+ bands. Moreover, our results indicate that the human genome contains about 30,000 genes.