Genes, isochores and bands in human chromosomes 21 and 22

The formation of chromatin domains involves a primary step based on the 3-D structure of DNA

The general model presented here for the formation of chromatin domains, LADs and TADs, is primarily based on the 3-D structures of the corresponding DNA sequences, the GC-poor and GC-rich isochores. Indeed, the low-heterogeneity GC-poor isochores locally are intrinsically stiff and curved because of the presence of interspersed oligo-Adenines. In contrast, the high-heterogeneity GC-rich isochores are in the shape of peaks characterized by increasing levels of GC and of interspersed oligo-Guanines. In LADs, oligo-Adenines induce local nucleosome depletions leading to structures that are well suited for the attachment to (and embedding in) the lamina. In TADs, the gradients of GC and of oligo-Guanines are responsible for a decreasing nucleosome density, decreasing supercoiling and increasing accessibility. This "moulding step" shapes the "primary TADs" into loops that lack self-interactions, being CTCF/cohesin-free structures. The cohesin complex then binds to the tips of "primary TADs" and slides down the loops, thanks to Nipbl, an essential factor for loading cohesin and for stimulating its ATPase activity and its translocation. This "extruding step" leads to closer contacts and to self-interactions in the loops and stops at the CTCF binding sites located at the base of the loops that are thus closed and insulated.

Genome Organization and Chromosome Architecture

How the same DNA sequences can function in the three-dimensional architecture of interphase nucleus, fold in the very compact structure of metaphase chromosomes, and go precisely back to the original interphase architecture in the following cell cycle remains an unresolved question to this day. The solution to this question presented here rests on the correlations that were found to hold between the isochore organization of the genome and the architecture of chromosomes from interphase to metaphase. The key points are the following: (1) The transition from the looped domains and subdomains of interphase chromatin to the 30-nm fiber loops of early prophase chromosomes goes through their unfolding into an extended chromatin structure (probably a 10-nm "beads-on-a-string" structure); (2) the architectural proteins of interphase chromatin, such as CTCF and cohesin subunits, are retained in mitosis and are part of the discontinuous protein scaffold of mitotic chromosomes; and (3) the conservation of the link between architectural proteins and their binding sites on DNA through the cell cycle explains the reversibility of the interphase to mitosis process and the "mitotic memory" of interphase architecture.

Chromosome Architecture and Genome Organization

How the same DNA sequences can function in the three-dimensional architecture of interphase nucleus, fold in the very compact structure of metaphase chromosomes and go precisely back to the original interphase architecture in the following cell cycle remains an unresolved question to this day. The strategy used to address this issue was to analyze the correlations between chromosome architecture and the compositional patterns of DNA sequences spanning a size range from a few hundreds to a few thousands Kilobases. This is a critical range that encompasses isochores, interphase chromatin domains and boundaries, and chromosomal bands. The solution rests on the following key points: 1) the transition from the looped domains and sub-domains of interphase chromatin to the 30-nm fiber loops of early prophase chromosomes goes through the unfolding into an extended chromatin structure (probably a 10-nm "beads-on-a-string" structure); 2) the architectural proteins of interphase chromatin, such as CTCF and cohesin sub-units, are retained in mitosis and are part of the discontinuous protein scaffold of mitotic chromosomes; 3) the conservation of the link between architectural proteins and their binding sites on DNA through the cell cycle explains the "mitotic memory" of interphase architecture and the reversibility of the interphase to mitosis process. The results presented here also lead to a general conclusion which concerns the existence of correlations between the isochore organization of the genome and the architecture of chromosomes from interphase to metaphase.

Segmenting the Human Genome into Isochores

The human genome is a mosaic of isochores, which are long (>200 kb) DNA sequences that are fairly homogeneous in base composition and can be assigned to five families comprising 33%–59% of GC composition. Although the compartmentalized organization of the mammalian genome has been investigated for more than 40 years, no satisfactory automatic procedure for segmenting the genome into isochores is available so far. We present a critical discussion of the currently available methods and a new approach called isoSegmenter which allows segmenting the genome into isochores in a fast and completely automatic manner. This approach relies on two types of experimentally defined parameters, the compositional boundaries of isochore families and an optimal window size of 100 kb. The approach represents an improvement over the existing methods, is ideally suited for investigating long-range features of sequenced and assembled genomes, and is publicly available at https://github.com/bunop/isoSegmenter.

The radial arrangement of the human chromosome 7 in the lymphocyte cell nucleus is associated with chromosomal band gene density

In the nuclei of human lymphocytes, chromosome territories are distributed according to the average gene density of each chromosome. However, chromosomes are very heterogeneous in size and base composition, and can contain both very gene-dense and very gene-poor regions. Thus, a precise analysis of chromosome organisation in the nuclei should consider also the distribution of DNA belonging to the chromosomal bands in each chromosome. To improve our understanding of the chromatin organisation, we localised chromosome 7 DNA regions, endowed with different gene densities, in the nuclei of human lymphocytes. Our results showed that this chromosome in cell nuclei is arranged radially with the gene-dense/GC-richest regions exposed towards the nuclear interior and the gene-poorest/GC-poorest ones located at the nuclear periphery. Moreover, we found that chromatin fibres from the 7p22.3 and the 7q22.1 bands are not confined to the territory of the bulk of this chromosome, protruding towards the inner part of the nucleus. Overall, our work demonstrates the radial arrangement of the territory of chromosome 7 in the lymphocyte nucleus and confirms that human genes occupy specific radial positions, presumably to enhance intra- and inter-chromosomal interaction among loci displaying a similar expression pattern, and/or similar replication timing.

Replication timing, chromosomal bands, and isochores

Chromosome replication timing is biphasic (early-late) in the cell cycle of vertebrates and of most (possibly all) eukaryotes. In the present work we have compared the extended, detailed replication timing maps that are available, namely those of human chromosomes 6, 11q, and 21q, with chromosomal bands as visualized at low (400 bands), high (850 bands), and highest (3,200 isochores) resolution. We have observed that the replicons located in a given isochore practically always show either all early or all late replication timing and that early-replicating isochores are short and GC-rich and late-replicating isochores are long and GC-poor. In the vast majority of cases, replicons are clustered in isochores, which are themselves most often clustered in early- or late-replication timing zones and may often reach the size of high-resolution bands and, very rarely, even that of low-resolution bands. Finally, we show that our results should be representative for the whole human genome and thus help to predict replication timing zones in all chromosomes.

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.

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.

Atypical mutational properties of human chromosome 21 suggested by comparative genome-scale analyses

Mutation of genetic material is a necessary component of evolutionary change. There is evidence for both intragenome and intergenome heterogeneity in terms of mutation frequencies. Reported comparisons of DNA sequence differences between human and chimpanzee (Pan troglodytes) suggest that human chromosome 21 may exhibit mutational hypervariability relative to the other autosomes. In the present study, further evidence is provided for such hypervariability based on large-scale analyses of amino acid composition of (translated) human genes and pseudogenes. A comparison of the variation in the above cases (i.e., DNA sequence differences and amino acid composition differences) yields similar ratios (1.2-1.4) for chromosome 21 relative to the other autosomes, e.g., human chromosome 22 - an autosome that is more typical in this respect and is of similar size to 21. Human chromosome 21 is also presented in this study as being atypical in terms of reported associations between mutation rates and GC content or CpG dinucleotides. In terms of GC distribution patterns, a comparison of NT_011512 and NT_011520 contigs revealed a lower heterogeneity for human chromosome 21 relative to 22. Possible hypermutability of chromosome 21 is further discussed in the context of GC patterns, reported long interspersed nuclear element content (LINE1s), and the implications of these parameters for chromatin structure.

Gene-rich and gene-poor chromosomal regions have different locations in the interphase nuclei of cold-blooded vertebrates

In situ hybridizations of single-copy GC-rich, gene-rich and GC-poor, gene-poor chicken DNA allowed us to localize the gene-rich and the gene-poor chromosomal regions in interphase nuclei of cold-blooded vertebrates. Our results showed that the gene-rich regions from amphibians (Rana esculenta) and reptiles (Podarcis sicula) occupy the more internal part of the nuclei, whereas the gene-poor regions occupy the periphery. This finding is similar to that previously reported in warm-blooded vertebrates, in spite of the lower GC levels of the gene-rich regions of cold-blooded vertebrates. This suggests that this similarity extends to chromatin structure, which is more open in the gene-rich regions of both mammals and birds and more compact in the gene-poor regions. In turn, this may explain why the compositional transition undergone by the genome at the emergence of homeothermy did not involve the entire ancestral genome but only a small part of it, and why it involved both coding and noncoding sequences. Indeed, the GC level increased only in that part of the genome that needed a thermodynamic stabilization, namely in the more open gene-rich chromatin of the nuclear interior, whereas the gene-poor chromatin of the periphery was stabilized by its own compact structure.

GC composition of the human genome: in search of isochores

The isochore theory, proposed nearly three decades ago, depicts the mammalian genome as a mosaic of long, fairly homogeneous genomic regions that are characterized by their guanine and cytosine (GC) content. The human genome, for instance, was claimed to consist of five distinct isochore families: L1, L2, H1, H2, and H3, with GC contents of <37%, 37%-42%, 42%-47%, 47%-52%, and >52%, respectively. In this paper, we address the question of the validity of the isochore theory through a rigorous sequence-based analysis of the human genome. Toward this end, we adopt a set of six attributes that are generally claimed to characterize isochores and statistically test their veracity against the available draft sequence of the complete human genome. By the selection criteria used in this study: distinctiveness, homogeneity, and minimal length of 300 kb, we identify 1,857 genomic segments that warrant the label "isochore." These putative isochores are nonuniformly scattered throughout the genome and cover about 41% of the human genome. We found that a four-family model of putative isochores is the most parsimonious multi-Gaussian model that can be fitted to the empirical data. These families, however, are GC poor, with mean GC contents of 35%, 38%, 41%, and 48% and do not resemble the five isochore families in the literature. Moreover, due to large overlaps among the families, it is impossible to classify genomic segments into isochore families reliably, according to compositional properties alone. These findings undermine the utility of the isochore theory and seem to indicate that the theory may have reached the limits of its usefulness as a description of genomic compositional 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.

Altered replication timing of the HIRA/Tuple1 locus in the DiGeorge and Velocardiofacial syndromes

DiGeorge and Velocardiofacial syndromes (DGS/VCFS) are endowed by a similar complex phenotype including cardiovascular, craniofacial, and thymic malformations, and are associated with heterozygous deletions of 22q11 chromosomal band. The Typically Deleted Region in the 22q11.21 subband (here called TDR22) is very gene-dense, and the extent of the deletion has been defined precisely in several studies. However, to date there is no evidence for a mechanism of haploinsufficiency that can fully explain the DGS/VCFS phenotype. In this study, we show that the candidate gene HIRA/Tuple1 mapping on the non-deleted TDR22, in DGS/VCFS subjects presents a delayed replication timing. Moreover, we observed an increase in the cell ratio showing the HIRA/Tuple1 locus localised toward the nuclear periphery. It is known that replication timing and nuclear location are generally correlated to the transcription activity of the relative DNA region. We propose that the alteration in the replication/nuclear location pattern of the non-deleted TDR22 indicates an altered gene regulation hence an altered transcritpion in DGS/VCFS.

Gene density and transcription influence the localization of chromatin outside of chromosome territories detectable by FISH

Genes can be transcribed from within chromosome territories; however, the major histocompatibilty complex locus has been reported extending away from chromosome territories, and the incidence of this correlates with transcription from the region. A similar result has been seen for the epidermal differentiation complex region of chromosome 1. These data suggested that chromatin decondensation away from the surface of chromosome territories may result from, and/or may facilitate, transcription of densely packed genes subject to coordinate regulation.To investigate whether localization outside of the visible confines of chromosome territories can also occur for regions that are not coordinately regulated, we have examined the spatial organization of human 11p15.5 and the syntenic region on mouse chromosome 7. This region is gene rich but its genes are not coordinately expressed, rather overall high levels of transcription occur in several cell types. We found that chromatin from 11p15.5 frequently extends away from the chromosome 11 territory. Localization outside of territories was also detected for other regions of high gene density and high levels of transcription. This is shown to be partly dependent on ongoing transcription. We suggest that local gene density and transcription, rather than the activity of individual genes, influences the organization of chromosomes in the nucleus.

Initial sequencing and comparative analysis of the mouse genome

The sequence of the mouse genome is a key informational tool for understanding the contents of the human genome and a key experimental tool for biomedical research. Here, we report the results of an international collaboration to produce a high-quality draft sequence of the mouse genome. We also present an initial comparative analysis of the mouse and human genomes, describing some of the insights that can be gleaned from the two sequences. We discuss topics including the analysis of the evolutionary forces shaping the size, structure and sequence of the genomes; the conservation of large-scale synteny across most of the genomes; the much lower extent of sequence orthology covering less than half of the genomes; the proportions of the genomes under selection; the number of protein-coding genes; the expansion of gene families related to reproduction and immunity; the evolution of proteins; and the identification of intraspecies polymorphism.

Are isochore sequences homogeneous?

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

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.

Localization of the gene-richest and the gene-poorest isochores in the interphase nuclei of mammals and birds

At a resolution of 850 bands, human chromosomes comprise two subsets of bands, the GC-richest H3(+) and the GC-poorest L1(+) bands, accounting for about 17 and 26%, respectively, of all bands. The former are a subset of the R bands and the latter are a subset of the G bands. These bands showed the highest and the lowest gene densities, respectively, as well as a number of other distinct features. Here we report that human and chicken interphase nuclei are characterized by the following features. (1) The gene-richest/GC-richest chromosomal regions are predominantly distributed in internal locations, whereas the gene-poorest/GC-poorest DNA regions are close to the nuclear envelope. (2) The interphase chromosomes seem to be characterized by a polar arrangement, because the gene-richest/GC-richest bands and the gene-poorest/GC-poorest bands are predominantly located in the distal and proximal regions, respectively, of chromosomes, and because interphase chromosomes are extremely long. While this polar arrangement is evident in the larger chromosomes, it is not displayed by the chicken microchromosomes and by some small human chromosomes, namely by chromosomes that are almost only composed by GC-rich or by GC-poor DNA. (3) The gene-richest chromosomal regions display a much more spread-out conformation compared to the gene-poorest regions in human nuclei. This finding has interesting implications for the formation of GC-rich isochores of warm-blooded vertebrates.

Misunderstandings about isochores. Part 1

A few months ago the International Human Genome Sequencing Consortium (IHGSC) published a 61-page paper on the human genome (IHGSC, Nature 409 (2001) 860). Here comments will be presented on some points of the paper that were previously investigated in our laboratory, and some misunderstandings and misconceptions about the organization and the evolutionary history of the human genome will be discussed. A very recent article on the same subject (Eyre-Walker and Hurst, Nat. Rev. Genet. 2 (2001) 549) will also be addressed. The present paper is a complement to two review articles which were published last year (Bernardi, Gene 241 (2000) 3; Gene 259(1) (2000) 31).