GC composition of the human genome: in search of isochores

Characterization of in vitro and in vivo hypomethylating effects of decitabine in acute myeloid leukemia by a rapid, specific and sensitive LC-MS/MS method

DNA hypermethylation is a common finding in malignant cells and has been explored as a therapeutic target for hypomethylating agents (e.g., decitabine). Detection of changes in DNA methylation might serve as a pharmacodynamic endpoint to establish the biological activity of these agents and predict clinical response. We developed and validated a rapid, sensitive and specific LC-MS/MS method for determination of global DNA methylation (GDM) in vitro and in vivo. Ratios of 5-methyl-2'-deoxycytidine (5mdC) to the internal standard 2-deoxyguanosine (2dG) in mass signal were used to quantify GDM levels. The assay was validated in a linear range from 40 fmol to 200 pmol 5mdC. The intra-day precision values ranged from 2.8 to 9.9% and the inter-day values from 1.1 to 15.0%. The accuracy of the assay varied between 96.7 and 109.5%. This method was initially applied for characterization of decitabine-induced GDM changes in in-vitro-treated leukemia cells. Following exposure to 2.5 microM decitabine, GDM decreased to approximately 50% of the baseline value. The clinical applicability of this method was then demonstrated in bone marrow samples from patients with acute myeloid leukemia treated with decitabine. Our data support the use of our LC-MS/MS method for clinical pharmacodynamic determination of changes in GDM in vivo.

Insights into social insects from the genome of the honeybee Apis mellifera

Here we report the genome sequence of the honeybee Apis mellifera, a key model for social behaviour and essential to global ecology through pollination. Compared with other sequenced insect genomes, the A. mellifera genome has high A+T and CpG contents, lacks major transposon families, evolves more slowly, and is more similar to vertebrates for circadian rhythm, RNA interference and DNA methylation genes, among others. Furthermore, A. mellifera has fewer genes for innate immunity, detoxification enzymes, cuticle-forming proteins and gustatory receptors, more genes for odorant receptors, and novel genes for nectar and pollen utilization, consistent with its ecology and social organization. Compared to Drosophila, genes in early developmental pathways differ in Apis, whereas similarities exist for functions that differ markedly, such as sex determination, brain function and behaviour. Population genetics suggests a novel African origin for the species A. mellifera and insights into whether Africanized bees spread throughout the New World via hybridization or displacement.

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.

Repetitive sequence environment distinguishes housekeeping genes

Housekeeping genes are expressed across a wide variety of tissues. Since repetitive sequences have been reported to influence the expression of individual genes, we employed a novel approach to determine whether housekeeping genes can be distinguished from tissue-specific genes by their repetitive sequence context. We show that Alu elements are more highly concentrated around housekeeping genes while various longer (>400-bp) repetitive sequences ("repeats"), including Long Interspersed Nuclear Element-1 (LINE-1) elements, are excluded from these regions. We further show that isochore membership does not distinguish housekeeping genes from tissue-specific genes and that repetitive sequence environment distinguishes housekeeping genes from tissue-specific genes in every isochore. The distinct repetitive sequence environment, in combination with other previously published sequence properties of housekeeping genes, was used to develop a method of predicting housekeeping genes on the basis of DNA sequence alone. Using expression across tissue types as a measure of success, we demonstrate that repetitive sequence environment is by far the most important sequence feature identified to date for distinguishing housekeeping genes.

Isochores exhibit evidence of genes interacting with the large-scale genomic environment

The genomes of mammals and birds can be partitioned into megabase-long regions, termed isochores, with consistently high, or low, average C + G content. Isochores with high CG contain a mixture of CG-rich and AT-rich genes, while high-AT isochores contain predominantly AT-rich genes. The two gene populations in the high-CG isochores are functionally distinguishable by statistical analysis of their gene ontology categories. However, the aggregate of the two populations in CG isochores is not statistically distinct from AT-rich genes in AT isochores. Genes tend to be located at local extrema of composition within the isochores, indicating that the CG-enriching mechanism acted differently when near to genes. On the other hand, maximum-likelihood reconstruction of molecular phylogenetic trees shows that branch lengths (evolutionary distances) for third codon positions in CG-rich genes are not substantially larger than those for AT-rich genes. In the context of neutral mutation theory this argues against any strong positive selection. Disparate features of isochores might be explained by a model in which about half of all genes functionally require AT richness, while, in warm-blooded organisms, about half the genome (in large coherent blocks) acquired a strong bias for mutations to CG. Using mutations in CG-rich genes as convenient indicators, we show that approximately 20% of amino acids in proteins are broadly substitutable, without regard to chemical similarity.

A new perspective on isochore evolution

The genomes of mammals and birds show dramatic variation in base composition over large scales, the so called isochore structure of the genome. The origin of isochores is still controversial: various neutral and selectionist models have been proposed--and criticized--since the discovery of isochores in the 1970s. The availability of complete mammalian genomes has yielded new opportunities for addressing this question. In particular, it was recently proposed that (i) the isochore structure is declining in many mammalian groups, and that (ii) GC-content is influenced by local recombination rate, possibly via the mechanism of GC-biased gene conversion. In this article we review the existing support for these two hypotheses, and discuss how they can be combined to provide a new perspective on isochore evolution.

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.

The Platypus Is in Its Place: Nuclear Genes and Indels Confirm the Sister Group Relation of Monotremes and Therians

Morphological data supports monotremes as the sister group of Theria (extant marsupials + eutherians), but phylogenetic analyses of 12 mitochondrial protein-coding genes have strongly supported the grouping of monotremes with marsupials: the Marsupionta hypothesis. Various nuclear genes tend to support Theria, but a comprehensive study of long concatenated sequences and broad taxon sampling is lacking. We therefore determined sequences from six nuclear genes and obtained additional sequences from the databases to create two large and independent nuclear data sets. One (data set I) emphasized taxon sampling and comprised five genes, with a concatenated length of 2,793 bp, from 21 species (two monotremes, six marsupials, nine placentals, and four outgroups). The other (data set II) emphasized gene sampling and comprised eight genes and three proteins, with a concatenated length of 10,773 bp or 3,669 amino acids, from five taxa (a monotreme, a marsupial, a rodent, human, and chicken). Both data sets were analyzed by parsimony, minimum evolution, maximum likelihood, and Bayesian methods using various models and data partitions. Data set I gave bootstrap support values for Theria between 55% and 100%, while support for Marsupionta was at most 12.3%. Taking base compositional bias into account generally increased the support for Theria. Data set II exclusively supported Theria, with the highest possible values and significantly rejected Marsupionta. Independent phylogenetic evidence in support of Theria was obtained from two single amino acid deletions and one insertion, while no supporting insertions and deletions were found for Marsupionta. On the basis of our data sets, the time of divergence between Monotremata and Theria was estimated at 231-217 MYA and between Marsupialia and Eutheria at 193-186 MYA. The morphological evidence for a basal position of Monotremata, well separated from Theria, is thus fully supported by the available molecular data from nuclear genes.

Human microRNAs target a functionally distinct population of genes with AT-rich 3' UTRs

While investigating microRNA targets, we have found that human genes divide into two roughly equal populations, based on the fraction of A plus T bases in their 3' UTRs. Using the Gene Ontology database, we find significant functional differences between the two gene populations, with AT-rich genes implicated in transcription and translation processes, and GC-rich genes implicated in signal transduction and posttranslational protein modification. Better understanding of the background distribution of nucleotides in 3' UTRs may allow improved prediction of microRNA-targeted genes in humans. We predict at least 1,200 KnownGene transcripts to be regulated by microRNAs. The large majority of these microRNA targets are in the AT-rich 3' UTR population. However, notwithstanding this preference for AT-rich targets, microRNA targets are found preferentially to be regulatory genes themselves, including both transcription factors and posttranslational modifiers. These results suggest that some processes involving mRNA, of which microRNA regulation may be just one, require AT-richness of 3' UTRs for functionality. A relationship, not simply one-to-one, between these 3' UTR populations and large-scale genomic isochores is described.

Variation in mutation dynamics across the maize genome as a function of regional and flanking base composition

We examine variation in mutation dynamics across a single genome (Zea mays ssp. mays) in relation to regional and flanking base composition using a data set of 10,472 SNPs generated by resequencing 1776 transcribed regions. We report several relationships between flanking base composition and mutation pattern. The A + T content of the two sites immediately flanking the mutation site is correlated with rate, transition bias, and GC --> AT pressure. We also observe a significant CpG effect, or increase in transition rate at CpG sites. At the regional level we find that the strength of the CpG effect is correlated with regional A + T content, ranging from a 1.7-fold increase in transition rate in relatively G + C-rich regions to a 2.6-fold increase in A + T-rich regions. We also observe a relationship between locus A + T content and GC --> AT pressure. This regional effect is in opposition to the influence of the two immediate neighbors in that GC --> AT pressure increases with increasing locus A + T content but decreases with increasing flanking base A + T content and may represent a relationship between genome location and mutation bias. The data indicate multiple context effects on mutations, resulting in significant variation in mutation dynamics across the genome.