Representing GC variation along eukaryotic chromosomes

Nanodosimetric Calculations of Radiation-Induced DNA Damage in a New Nucleus Geometrical Model Based on the Isochore Theory

Double-strand breaks (DSBs) in nuclear DNA represents radiation-induced damage that has been identified as particularly deleterious. Calculating this damage using Monte Carlo track structure modeling could be a suitable indicator to better assess and anticipate the side-effects of radiation therapy. However, as already demonstrated in previous work, the geometrical description of the nucleus and the DNA content used in the simulation significantly influence damage calculations. Therefore, in order to obtain accurate results, this geometry must be as realistic as possible. In this study, a new geometrical model of an endothelial cell nucleus and DNA distribution according to the isochore theory are presented and used in a Monte Carlo simulation chain based on the Geant4-DNA toolkit. In this theory, heterochromatin and euchromatin compaction are distributed along the genome according to five different families (L1, L2, H1, H2, and H3). Each of these families is associated with a different hetero/euchromatin rate related to its compaction level. In order to compare the results with those obtained using a previous nuclear geometry, simulations were performed for protons with linear energy transfers (LETs) of 4.29 keV/µm, 19.51 keV/µm, and 43.25 keV/µm. The organization of the chromatin fibers at different compaction levels linked to isochore families increased the DSB yield by 6-10%, and it allowed the most affected part of the genome to be identified. These new results indicate that the genome core is more radiosensitive than the genome desert, with a 3-8% increase in damage depending on the LET. This work highlights the importance of using realistic distributions of chromatin compaction levels to calculate radio-induced damage using Monte Carlo simulation methods.

A common genomic code for chromatin architecture and recombination landscape

Recent findings established a link between DNA sequence composition and interphase chromatin architecture and explained the evolutionary conservation of TADs (Topologically Associated Domains) and LADs (Lamina Associated Domains) in mammals. This prompted us to analyse conformation capture and recombination rate data to study the relationship between chromatin architecture and recombination landscape of human and mouse genomes. The results reveal that: (1) low recombination domains and blocks of elevated linkage disequilibrium tend to coincide with TADs and isochores, indicating co-evolving regulatory elements and genes in insulated neighbourhoods; (2) double strand break (DSB) and recombination frequencies increase in the short loops of GC-rich TADs, whereas recombination cold spots are typical of LADs and (3) the binding and loading of proteins, which are critical for DSB and meiotic recombination (SPO11, DMC1, H3K4me3 and PRMD9) are higher in GC-rich TADs. One explanation for these observations is that the occurrence of DSB and recombination in meiotic cells are associated with compositional and epigenetic features (genomic code) that influence DNA stiffness/flexibility and appear to be similar to those guiding the chromatin architecture in the interphase nucleus of pre-leptotene cells.

An Isochore Framework Underlies Chromatin Architecture

A recent investigation showed the existence of correlations between the architectural features of mammalian interphase chromosomes and the compositional properties of isochores. This result prompted us to compare maps of the Topologically Associating Domains (TADs) and of the Lamina Associated Domains (LADs) with the corresponding isochore maps of mouse and human chromosomes. This approach revealed that: 1) TADs and LADs correspond to isochores, i.e., isochores are the genomic units that underlie chromatin domains; 2) the conservation of TADs and LADs in mammalian genomes is explained by the evolutionary conservation of isochores; 3) chromatin domains corresponding to GC-poor isochores (e.g., LADs) show not only self-interactions but also intrachromosomal interactions with other domains also corresponding to GC-poor isochores even if located far away; in contrast, chromatin domains corresponding to GC-rich isochores (e.g., TADs) show more localized chromosomal interactions, many of which are inter-chromosomal. In conclusion, this investigation establishes a link between DNA sequences and chromatin architecture, explains the evolutionary conservation of TADs and LADs and provides new information on the spatial distribution of GC-poor/gene-poor and GC-rich/gene-rich chromosomal regions in the interphase nucleus.

The Anolis Lizard Genome: An Amniote Genome without Isochores?

Two articles published 5 years ago concluded that the genome of the lizard Anolis carolinensis is an amniote genome without isochores. This claim was apparently contradicting previous results on the general presence of an isochore organization in all vertebrate genomes tested (including Anolis). In this investigation, we demonstrate that the Anolis genome is indeed heterogeneous in base composition, since its macrochromosomes comprise isochores mainly from the L2 and H1 families (a moderately GC-poor and a moderately GC-rich family, respectively), and since the majority of the sequenced microchromosomes consists of H1 isochores. These families are associated with different features of genome structure, including gene density and compositional correlations (e.g., GC3 vs flanking sequence GC and intron GC), as in the case of mammalian and avian genomes. Moreover, the assembled Anolis chromosomes have an enormous number of gaps, which could be due to sequencing problems in GC-rich regions of the genome. In conclusion, the Anolis genome is no exception to the general rule of an isochore organization in the genomes of vertebrates (and other eukaryotes).

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.

Best practices for mapping replication origins in eukaryotic chromosomes

Understanding the regulatory principles ensuring complete DNA replication in each cell division is critical for deciphering the mechanisms that maintain genomic stability. Recent advances in genome sequencing technology facilitated complete mapping of DNA replication sites and helped move the field from observing replication patterns at a handful of single loci to analyzing replication patterns genome-wide. These advances address issues, such as the relationship between replication initiation events, transcription, and chromatin modifications, and identify potential replication origin consensus sequences. This unit summarizes the technological and fundamental aspects of replication profiling and briefly discusses novel insights emerging from mining large datasets, published in the last 3 years, and also describes DNA replication dynamics on a whole-genome scale.

Compositional patterns in the genomes of unicellular eukaryotes

Background

The genomes of multicellular eukaryotes are compartmentalized in mosaics of isochores, large and fairly homogeneous stretches of DNA that belong to a small number of families characterized by different average GC levels, by different gene concentration (that increase with GC), different chromatin structures, different replication timing in the cell cycle, and other different properties. A question raised by these basic results concerns how far back in evolution the compartmentalized organization of the eukaryotic genomes arose.

Results

In the present work we approached this problem by studying the compositional organization of the genomes from the unicellular eukaryotes for which full sequences are available, the sample used being representative. The average GC levels of the genomes from unicellular eukaryotes cover an extremely wide range (19%-60% GC) and the compositional patterns of individual genomes are extremely different but all genomes tested show a compositional compartmentalization.

Conclusions

The average GC range of the genomes of unicellular eukaryotes is very broad (as broad as that of prokaryotes) and individual compositional patterns cover a very broad range from very narrow to very complex. Both features are not surprising for organisms that are very far from each other both in terms of phylogenetic distances and of environmental life conditions. Most importantly, all genomes tested, a representative sample of all supergroups of unicellular eukaryotes, are compositionally compartmentalized, a major difference with prokaryotes.

Two common profiles exist for genomic oligonucleotide frequencies

Background

It was reported that there is a majority profile for trinucleotide frequencies among genomes. And further study has revealed that two common profiles, rather than one majority profile, exist for genomic trinucleotide frequencies. However, the origins of the common/majority profile remain elusive. Moreover, it is not clear whether the features of common profile may be extended to oligonucleotides other than trinucleotides.

Findings

We analyzed 571 prokaryotic genomes (chromosomes) and some selected eukaryotic nuclear genomes as well as other genetic systems to study their compositional features. We found that there are also two common profiles for genomic oligonucleotide frequencies: one is from low-GC content genomes, and the other is from high-GC content genomes. Furthermore, each common profile is highly correlated to the average profile of random sequences with corresponding GC content and generated according to first-order symmetry.

Conclusions

The causes for the existence of two common profiles would mainly be GC content variations and strand symmetry of genomic sequences. Therefore, both GC content and strand symmetry would play important roles in genome evolution.

Unraveling cell type-specific and reprogrammable human replication origin signatures associated with G-quadruplex consensus motifs

DNA replication is highly regulated, ensuring faithful inheritance of genetic information through each cell cycle. In metazoans, this process is initiated at many thousands of DNA replication origins whose cell type-specific distribution and usage are poorly understood. We exhaustively mapped the genome-wide location of replication origins in human cells using deep sequencing of short nascent strands and identified ten times more origin positions than we expected; most of these positions were conserved in four different human cell lines. Furthermore, we identified a consensus G-quadruplex-forming DNA motif that can predict the position of DNA replication origins in human cells, accounting for their distribution, usage efficiency and timing. Finally, we discovered a cell type-specific reprogrammable signature of cell identity that was revealed by specific efficiencies of conserved origin positions and not by the selection of cell type-specific subsets of origins.

The isochore patterns of invertebrate genomes

Background

Previous investigations from our laboratory were largely focused on the genome organization of vertebrates. We showed that these genomes are mosaics of isochores, megabase-size DNA sequences that are fairly homogeneous in base composition yet belong to a small number of families that cover a wide compositional spectrum. A question raised by these results concerned how far back in evolution an isochore organization of the eukaryotic genome arose.

Results

The present investigation deals with the compositional patterns of the invertebrates for which full genome sequences, or at least scaffolds, are available. We found that (i) a mosaic of isochores is the long-range organization of all the genomes that we investigated; (ii) the isochore families from the invertebrate genomes matched the corresponding families of vertebrates in GC levels; (iii) the relative amounts of isochore families were remarkably different for different genomes, except for those from phylogenetically close species, such as the Drosophilids.

Conclusion

This work demonstrates not only that an isochore organization is present in all metazoan genomes analyzed that included Nematodes, Arthropods among Protostomia, Echinoderms and Chordates among Deuterostomia, but also that the isochore families of invertebrates share GC levels with the corresponding families of vertebrates.

Classifying coding DNA with nucleotide statistics

In this report, we compared the success rate of classification of coding sequences (CDS) vs. introns by Codon Structure Factor (CSF) and by a method that we called Universal Feature Method (UFM). UFM is based on the scoring of purine bias (Rrr) and stop codon frequency. We show that the success rate of CDS/intron classification by UFM is higher than by CSF. UFM classifies ORFs as coding or non-coding through a score based on (i) the stop codon distribution, (ii) the product of purine probabilities in the three positions of nucleotide triplets, (iii) the product of Cytosine (C), Guanine (G), and Adenine (A) probabilities in the 1st, 2nd, and 3rd positions of triplets, respectively, (iv) the probabilities of G in 1st and 2nd position of triplets and (v) the distance of their GC3 vs. GC2 levels to the regression line of the universal correlation. More than 80% of CDSs (true positives) of Homo sapiens (>250 bp), Drosophila melanogaster (>250 bp) and Arabidopsis thaliana (>200 bp) are successfully classified with a false positive rate lower or equal to 5%. The method releases coding sequences in their coding strand and coding frame, which allows their automatic translation into protein sequences with 95% confidence. The method is a natural consequence of the compositional bias of nucleotides in coding sequences.

Mapping insertions, deletions and SNPs on Venter's chromosomes

BACKGROUND

The very recent availability of fully sequenced individual human genomes is a major revolution in biology which is certainly going to provide new insights into genetic diseases and genomic rearrangements.

RESULTS

We mapped the insertions, deletions and SNPs (single nucleotide polymorphisms) that are present in Craig Venter's genome, more precisely on chromosomes 17 to 22, and compared them with the human reference genome hg17. Our results show that insertions and deletions are almost absent in L1 and generally scarce in L2 isochore families (GC-poor L1+L2 isochores represent slightly over half of the human genome), whereas they increase in GC-rich isochores, largely paralleling the densities of genes, retroviral integrations and Alu sequences. The distributions of insertions/deletions are in striking contrast with those of SNPs which exhibit almost the same density across all isochore families with, however, a trend for lower concentrations in gene-rich regions.

CONCLUSIONS

Our study strongly suggests that the distribution of insertions/deletions is due to the structure of chromatin which is mostly open in gene-rich, GC-rich isochores, and largely closed in gene-poor, GC-poor isochores. The different distributions of insertions/deletions and SNPs are clearly related to the two different responsible mechanisms, namely recombination and point mutations.

The evolution of isochore patterns in vertebrate genomes

Background

Previous work from our laboratory showed that (i) vertebrate genomes are mosaics of isochores, typically megabase-size DNA segments that are fairly homogeneous in base composition; (ii) isochores belong to a small number of families (five in the human genome) characterized by different GC levels; (iii) isochore family patterns are different in fishes/amphibians and mammals/birds, the latter showing GC-rich isochore families that are absent or very scarce in the former; (iv) there are two modes of genome evolution, a conservative one in which isochore patterns basically do not change (e.g., among mammalian orders), and a transitional one, in which they do change (e.g., between amphibians and mammals); and (v) isochores are tightly linked to a number of basic biological properties, such as gene density, gene expression, replication timing and recombination.

Results

The present availability of a number of fully sequenced genomes ranging from fishes to mammals allowed us to carry out investigations that (i) more precisely quantified our previous conclusions; (ii) showed that the different isochore families of vertebrate genomes are largely conserved in GC levels and dinucleotide frequencies, as well as in isochore size; and (iii) isochore family patterns can be either conserved or change within both warm- and cold-blooded vertebrates.

Conclusion

On the basis of the results presented, we propose that (i) the large conservation of GC levels and dinucleotide frequencies may reflect the conservation of chromatin structures; (ii) the conservation of isochore size may be linked to the role played by isochores in chromosome structure and replication; (iii) the formation, the maintainance and the changes of isochore patterns are due to natural selection.

A de novo apparently balanced translocation [46,XY,t(2;9)(p13;p24)] interrupting RAB11FIP5 identifies a potential candidate gene for autism spectrum disorder

Autism spectrum disorder (ASD) is a severe developmental disorder of the central nervous system characterized by impairments in social interaction, communication, and range of interests and behaviors. The syndrome's prevalence is estimated to be as high as 1 in 150 American children yet its etiology remains largely unknown. Examination of observed cytogenetic variants in individuals with ASD may identify genes involved in its pathogenesis. As part of a multidisciplinary study, an apparently balanced de novo translocation between chromosomes 2 and 9 [46,XY,t(2;9)(p13;p24)] was identified in a subject with pervasive developmental disorder not otherwise specified (PDD-NOS), and no distinctive dysmorphic features. Molecular characterization of the rearrangement revealed direct interruption of the RAB11 family interacting protein 5 (RAB11FIP5) gene. RAB11FIP5 is a Rab effector involved in protein trafficking from apical recycling endosomes to the apical plasma membrane. It is ubiquitously expressed and reported to contribute to both neurotransmitter release and neurotransmitter uptake at the synaptic junction. Detailed analysis of the rearrangement breakpoints suggests that the reciprocal translocation may have formed secondary to incorrect repair of double strand breaks (DSBs) by nonhomologous end-joining (NHEJ).

Genomic DNA from animals shows contrasting strand bias in large and small subsequences

Background

For eukaryotes, there is almost no strand bias with regard to base composition, with exceptions for origins of replication and transcription start sites and transcribed regions. This paper revisits the question for subsequences of DNA taken at random from the genome.

Results

For a typical mammal, for example mouse or human, there is a small strand bias throughout the genomic DNA: there is a correlation between (G - C) and (A - T) on the same strand, (that is between the difference in the number of guanine and cytosine bases and the difference in the number of adenine and thymine bases). For small subsequences – up to 1 kb – this correlation is weak but positive; but for large windows – around 50 kb to 2 Mb – the correlation is strong and negative. This effect is largely independent of GC%. Transcribed and untranscribed regions give similar correlations both for small and large subsequences, but there is a difference in these regions for intermediate sized subsequences. An analysis of the human genome showed that position within the isochore structure did not affect these correlations. An analysis of available genomes of different species shows that this contrast between large and small windows is a general feature of mammals and birds. Further down the evolutionary tree, other organisms show a similar but smaller effect. Except for the nematode, all the animals analysed showed at least a small effect.

Conclusion

The correlations on the large scale may be explained by DNA replication. Transcription may be a modifier of these effects but is not the fundamental cause. These results cast light on how DNA mutations affect the genome over evolutionary time. At least for vertebrates, there is a broad relationship between body temperature and the size of the correlation. The genome of mammals and birds has a structure marked by strand bias segments.

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.

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.

Identification of sequence motifs at the breakpoint junctions in three t(1;9)(p36.3;q34) and delineation of mechanisms involved in generating balanced translocations

Although approximately 1 in 500 individuals carries a reciprocal translocation, little is known about the mechanisms that result in their formation. We analyzed the sequences surrounding the breakpoints in three unbalanced translocations of 1p and 9q, all of which were designated t(1;9)(p36.3;q34), to investigate the presence of sequence motifs that might mediate nonhomologous end joining (NHEJ). The breakpoint regions were unique in all individuals. Two of three translocations demonstrated insertions and duplications at the junctions, suggesting NHEJ in the formation of the rearrangements. No homology was identified in the breakpoint regions, further supporting NHEJ. We found translin motifs at the breakpoint junctions, suggesting the involvement of translin in the joining of the broken chromosome ends. We propose a model for balanced translocation formation in humans similar to transposition in bacteria, in which staggered nicks are repaired resulting in duplications and insertions at the translocation breakpoints.

An isochore map of human chromosomes

Isochores are large DNA segments (>>300 kb on average) that are characterized by an internal variation in GC well below the full variation seen in the mammalian genome. Precisely defining in terms of size and composition as well as mapping the isochores on human chromosomes have, however, remained largely unsolved problems. Here we used a very simple approach to segment the human chromosomes de novo, based on assessments of GC and its variation within and between adjacent regions. We obtain a complete coverage of the human genome (neglecting the remaining gaps) by approximately 3200 isochores, which may be visualized as the ultimate chromosomal bands. Isochores visibly belong to five families characterized by different GC levels, as expected from previous investigations. Since we previously showed that isochores are tightly linked to basic biological properties such as gene density, replication timing, and recombination, the new level of detail provided by the isochore map will help the understanding of genome structure, function, and evolution.

Similar chromosomal changes in cisplatin and oxaliplatin-resistant sublines of the H69 SCLC cell line are not associated with platinum resistance

Small cell lung cancer (SCLC) initially responds well to DNA damaging drugs such as cisplatin, however this is transitory as resistance normally develops. To investigate whether changes in chromosomal copy number caused by platinum drug treatment contributes to platinum resistance, we have analyzed H69 SCLC cells and two low-level platinum-resistant sublines, H69CIS200 and H69OX400, derived by cisplatin and oxaliplatin treatment, respectively. Affymetrix 10K SNP array showed that cisplatin and oxaliplatin have independently caused similar changes including loss of segments 6q21-qter and 13pter-13q.14.11 and duplication of chromosome 21. Interestingly, despite using equally cytotoxic doses of drug in the development of the cell lines, oxaliplatin caused three times more chromosomal changes than cisplatin. The resistant cell lines lose their resistant phenotype after 3 months of drug-free culture. The revertant cell lines, denoted H69CIS200-S and H69OX400-S, were also analyzed by Affymetrix array to determine if chromosomal changes associated with resistance remain after the resistant phenotype is lost. In the H69OX400-S many of the changes observed in the resistant cells were absent suggesting that they contributed to the resistant phenotype including: loss of 1q23.3-qter, 10q11.23, and 19q13.12-q13.2 and duplication of segments 6p21.2-p12.3, 16q12.1-16q13, 16q21-q23.1, and 19q12. However, out of the similar changes induced by cisplatin and oxaliplatin, both the loss of 6q21-qter and gain of 21 were still present in the H69CIS200-S and H69OX400-S cells. This suggests that cisplatin and oxaliplatin induced similar changes due to inherent vulnerabilities in the H69 cells rather than changes associated with platinum resistance.

Genome complexity and repetitive DNA in metazoans from extreme marine environments

As genomics converges with ecology and evolution to identify the fundamental linkages between genome structure and function, genome and transcriptome complexity will need to be measured in organisms from more diverse habitats, most often in the absence of complete sequence data. Here, we describe the complexity of ten genomes measured by a novel, high-throughput fluorescence-based kinetic hybridization assay. We applied the Shannon information index, H, and a related, fluorescence-adjusted index, H(f), as unique metrics of the hybridization kinetics to complement the conventional rate constant, k. A strong, positive relationship was present between H(f), and the repetitive DNA content of five eukaryotic genomes previously determined by Cot kinetic analyses (Onchorynchus keta, Ilyanassa obsoleta, Bos taurus, Limulus polyphemus, Saccharyomyces cerevisiae). This relationship was used to characterize the complexity of previously unstudied genomic samples in five metazoan taxa from three marine environments, including deep-sea hydrothermal vents (Alvinella pompejana), the temperate subtidal (Streblospio benedicti), and Antarctic coastal bays (Sterechinus neumayeri, Odontaster validus, Tritonia antarctica). Contrary to the predictions of nucleotypic theory, Antarctic invertebrates consistently had the lowest quantities of repetitive DNA in conjunction with low metabolic rates and highly protracted rates of cell division and larval development. Conversely, hydrothermal vent species with rapid cell division and growth do not have significantly different genome sizes or particularly low amounts of repetitive DNA as compared to non-vent, deep-sea taxa. Furthermore, there appears to be a positive correlation between the temperature at which the most abundant repetitive sequence classes anneal and habitat thermal stability. Thus, our study reveals a potential shift in repetitive sequence representation between these extreme environments that may be related to genome function in species living at these different thermal regimes.

Amplification of whole tumor genomes and gene-by-gene mapping of genomic aberrations from limited sources of fresh-frozen and paraffin-embedded DNA

Sufficient quantity of genomic DNA can be a bottleneck in genome-wide analysis of clinical tissue samples. DNA polymerase Phi29 can be used for the random-primed amplification of whole genomes, although the amplification may introduce bias in gene dosage. We have performed a detailed investigation of this technique in archival fresh-frozen and formalin-fixed/paraffin-embedded tumor DNA by using cDNA microarray-based comparative genomic hybridization. Phi29 amplified DNA from matched pairs of fresh-frozen and formalin-fixed/paraffin-embedded tumor samples with similar efficiency. The distortion in gene dosage representation in the amplified DNA was nonrandom and reproducibly involved distinct genomic loci. Regional amplification efficiency was significantly linked to regional GC content of the template genome. The biased gene representation in amplified tumor DNA could be effectively normalized by using amplified reference DNA. Our data suggest that genome-wide gene dosage alterations in clinical tumor samples can be reliably assessed from a few hundred tumor cells. Therefore, this amplification method should lend itself to high-throughput genetic analyses of limited sources of tumor, such as fine-needle biopsies, laser-microdissected tissue, and small paraffin-embedded specimens.

Detection of nucleolar organizer and mitochondrial DNA insertion regions based on the isochore map of Arabidopsis thaliana

Eukaryotic genomes are composed of isochores, i.e. long sequences relatively homogeneous in GC content. In this paper, the isochore structure of Arabidopsis thaliana genome has been studied using a windowless technique based on the Z curve method and intuitive curves are drawn for all the five chromosomes. Using these curves, we can calculate the GC content at any resolution, even at the base level. It is observed that all the five chromosomes are composed of several GC-rich and AT-rich regions alternatively. Usually, these regions, named 'isochore-like regions', have large fluctuations in the GC content. Five isochores with little fluctuations are also observed. Detailed analyses have been performed for these isochores. A GC-rich 'isochore-like region' and a GC-isochore in chromosome II and IV, respectively, are the nucleolar organizer regions (NORs), and genes located in the two regions prefer to use GC-ending codons. Another GC-isochore located in chromosome II is a mitochondrial DNA insertion region, the position and size of this region is precisely predicted by the current method. The amino acid usage and codon preference of genes in this organellar-to-nuclear transfer region show significant difference from other regions. Moreover, the centromeres are located in GC-rich 'isochore-like regions' in all the five chromosomes. The current method can provide a useful tool for analyzing whole genomic sequences of eukaryotes.

The maize gene space is compositionally compartimentalized

Previous investigations by Southern hybridization of cDNA with compositional DNA fractions showed that the majority of maize genes are located in a narrow GC range of DNA fragments and that the corresponding gene space was GC-richer than the region of the genome where zein genes are found. Here, we revisited the maize gene space using new data from the maize genome sequencing initiative. We found that the maize gene space itself is formed of two compositional compartments, i.e., a GC-poor and a GC-rich, characterized by a different distribution of Opie and Huck retrotransposons. The GC-rich compartment tends to be richer in GC-rich genes than the GC-poor compartment. However, the gene space compartimentalization of maize is much simpler than that of human.

The Biased Distribution of Alus in Human Isochores Might Be Driven by Recombination

Alu retrotransposons do not show a homogeneous distribution over the human genome but have a higher density in GC-rich (H) than in AT-rich (L) isochores. However, since they preferentially insert into the L isochores, the question arises: What is the evolutionary mechanism that shifts the Alu density maximum from L to H isochores? To disclose the role played by each of the potential mechanisms involved in such biased distribution, we carried out a genome-wide analysis of the density of the Alus as a function of their evolutionary age, isochore membership, and intron vs. intergene location. Since Alus depend on the retrotransposase encoded by the LINE1 elements, we also studied the distribution of LINE1 to provide a complete evolutionary scenario. We consecutively check, and discard, the contributions of the Alu/LINE1 competition for retrotransposase, compositional matching pressure, and Alu overrepresentation in introns. In analyzing the role played by unequal recombination, we scan the genome for Alu trimers, a direct product of Alu-Alu recombination. Through computer simulations, we show that such trimers are much more frequent than expected, the observed/expected ratio being higher in L than in H isochores. This result, together with the known higher selective disadvantage of recombination products in H isochores, points to Alu-Alu recombination as the main agent provoking the density shift of Alus toward the GC-rich parts of the genome. Two independent pieces of evidence-the lower evolutionary divergence shown by recently inserted Alu subfamilies and the higher frequency of old stand-alone Alus in L isochores-support such a conclusion. Other evolutionary factors, such as population bottlenecks during primate speciation, may have accelerated the fast accumulation of Alus in GC-rich isochores.

Breakpoints of variant 9;22 translocations in chronic myeloid leukemia locate preferentially in the CG-richest regions of the genome

From 5% to 10% of 9;22 translocations in chronic myeloid leukemia (CML) are reported to occur in variant form, that is, with the involvement of other regions of the genome in 3-way or more rearrangements. The literature indicates that the alternative breakpoints are not distributed randomly in the genome but show hotspots. We present data on 289 unpublished cases of CML with variant 9;22 translocations having a total of 342 variant breakpoints, the largest independent series to date. We found that the distribution of breaks was in loose agreement with the literature but that some new hotspots were identified; furthermore, some published hotspots were not fully supported by our data. Moreover, when our 342 variant breakpoints were plotted against profiles of CG heterogeneity in the genome, a significant positive correlation between breakpoint locations and CG composition was observed. In an ancillary study, we compared the frequency of variant t(9;22) with that of variants of t(15;17) associated with acute promyelocytic leukemia (AML M3). We found that the frequency of the former, 9.3%, was significantly higher than that of the latter, 2.6%.