The compositional distribution of coding sequences and DNA molecules in humans and murids

Summary The compositional distributions of coding sequences and DNA molecules (in the 50-100-kb range) are remarkably narrower in murids (rat and mouse) compared to humans (as well as to all other mammals explored so far). In murids, both distributions begin at higher and end at lower GC values. A comparison of homologous coding sequences from murids and humans revealed that their different compositional distributions are due to differences in GC levels in all three codon positions, particularly of genes located at both ends of the distribution. In turn, these differences are responsible for differences in both codon usage and amino acids. When GC levels at first+second codon positions and third codon positions, respectively, of murid genes are plotted against corresponding GC levels of homologous human genes, linear relationships (with very high correlation coefficients and slopes of about 0.78 and 0.60, respectively) are found. This indicates a conservation of the order of GC levels in homologous genes from humans and murids. (The same comparison for mouse and rat genes indicates a conservation of GC levels of homologous genes.) A similar linear relationship was observed when plotting GC levels of corresponding DNA fractions (as obtained by density gradient centrifugation in the presence of a sequence-specific ligand) from mouse and human. These findings indicate that orderly compositional changes affecting not only coding sequences but also noncoding sequences took place since the divergence of murids. Such directional fixations of mutations point to the existence of selective pressures affecting the genome as a whole.     

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

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

Relationship between base composition in non-coding DNA of genes and codon composition

The C + G percentage in third position of codons is linearly dependent on the C + G composition of flanking regions and introns. A similar relationship is shown for the first and second position which significantly influence the nature of amino acid sequence. If mutations would be oriented according to the local base composition, this will imply that genes of the same multigenic family would evolve at different rate.     

Suppressor of cytokine signaling 1 interacts with the macrophage colony-stimulating factor receptor and negatively regulates its proliferation signal

Macrophage colony-stimulating factor receptor (M-CSF-R) is a tyrosine kinase that regulates proliferation, differentiation, and cell survival during monocytic lineage development. Upon activation, M-CSF-R dimerizes and autophosphorylates on specific tyrosines, creating binding sites for several cytoplasmic SH2-containing signaling molecules that relay and modulate the M-CSF signal. Here we show that M-CSF-R interacts with suppressor of cytokine signaling 1 (Socs1), a negative regulator of various cytokine and growth factor signaling pathways. Using the yeast two-hybrid system, in vitro glutathione S-transferase-M-CSF-R pull-down, and in vivo coimmunoprecipitation experiments, we demonstrated a direct interaction between the SH2 domain of Socs1 and phosphorylated tyrosines 697 or 721 of the M-CSF-R kinase insert region. Moreover, Socs1 is tyrosine-phosphorylated in response to M-CSF. Ectopic expression of Socs1 in FDC-P1/MAC and EML hematopoietic cell lines decreased their growth rates in the presence of limiting concentrations of M-CSF. However, Socs1 expression did not totally suppress long term cell growth in the presence of saturating M-CSF concentrations, in contrast to other cytokines such as stem cell factor and interleukin 3. Taken together, these results suggest that Socs1 is an M-CSF-R-binding partner involved in negative regulation of proliferation signaling and that it differentially affects cytokine receptor signals.     

Neutral effect of recombination on base composition in Drosophila

Recombination is thought to have various evolutionary effects on genome evolution. In this study, we investigated the relationship between the base composition and recombination rate in the Drosophila melanogaster genome. Because of a current debate about the accuracy of the estimates of recombination rate in Drosophila, we used eight different measures of recombination rate from recent work. We confirmed that the G + C content of large introns and flanking regions is positively correlated with recombination rate, suggesting that recombination has a neutral effect on base composition in Drosophila. We also confirmed that this neutral effect of recombination is the main determinant of the correlation between synonymous codon usage bias and recombination rate in Drosophila.     

Characterization of the promoter controlling Mona/Gads expression in the megakaryocytic lineage

The deletion of a 260-kb segment containing all the coding DNA sequences (CDS) of chromosome 1 of Leishmania major Friedlin strain was performed through homologous recombination during a transfection experiment. This allowed the selection of a mutant clone containing a linear extra chromosome sizing 155 kb (XC155). The structure of XC155 was determined by restriction analysis and DNA cloning and sequencing of the gel-purified chromosome: it is made of a 'mirror' inverted duplication of the 'right' end of chromosome 1a (approximately 25 kb at each end), and in its central part of a complex tandem amplification of the linearized transfection vector containing the hygromycin resistance gene (over approximately 105 kb). No sequence of the coding region of chromosome 1 (including the 1.6-kb 'switch' region) was found. By contrast, XC155 contains two large (approximately 13 kb) clusters of tandemly repeated subtelomeric sequences (272-bp 'satellite' DNA) as well as telomeric hexamer repeats. This extra chromosome was found to be mitotically stable after >150 generations without selective pressure in vitro. Two sequence elements are considered which may have an effect on mitotic stability and participate to centromeric function in this extra chromosome: the amplification of the input vector and the 272-bp 'satellite' DNA bound by telomeric repeats.     

Vanishing GC-rich isochores in mammalian genomes

To understand the origin and evolution of isochores-the peculiar spatial distribution of GC content within mammalian genomes-we analyzed the synonymous substitution pattern in coding sequences from closely related species in different mammalian orders. In primate and cetartiodactyls, GC-rich genes are undergoing a large excess of GC --> AT substitutions over AT --> GC substitutions: GC-rich isochores are slowly disappearing from the genome of these two mammalian orders. In rodents, our analyses suggest both a decrease in GC content of GC-rich isochores and an increase in GC-poor isochores, but more data will be necessary to assess the significance of this pattern. These observations question the conclusions of previous works that assumed that base composition was at equilibrium. Analysis of allele frequency in human polymorphism data, however, confirmed that in the GC-rich parts of the genome, GC alleles have a higher probability of fixation than AT alleles. This fixation bias appears not strong enough to overcome the large excess of GC --> AT mutations. Thus, whatever the evolutionary force (neutral or selective) at the origin of GC-rich isochores, this force is no longer effective in mammals. We propose a model based on the biased gene conversion hypothesis that accounts for the origin of GC-rich isochores in the ancestral amniote genome and for their decline in present-day mammals.     

Determinants of substitution rates in mammalian genes: expression pattern affects selection intensity but not mutation rate

To determine whether gene expression patterns affect mutation rates and/or selection intensity in mammalian genes, we studied the relationships between substitution rates and tissue distribution of gene expression. For this purpose, we analyzed 2,400 human/rodent and 834 mouse/rat orthologous genes, and we measured (using expressed sequence tag data) their expression patterns in 19 tissues from three development states. We show that substitution rates at nonsynonymous sites are strongly negatively correlated with tissue distribution breadth: almost threefold lower in ubiquitous than in tissue-specific genes. Nonsynonymous substitution rates also vary considerably according to the tissues: the average rate is twofold lower in brain-, muscle-, retina- and neuron-specific genes than in lymphocyte-, lung-, and liver-specific genes. Interestingly, 5' and 3' untranslated regions (UTRs) show exactly the same trend. These results demonstrate that the expression pattern is an essential factor in determining the selective pressure on functional sites in both coding and noncoding regions. Conversely, silent substitution rates do not vary with expression pattern, even in ubiquitously expressed genes. This latter result thus suggests that synonymous codon usage is not constrained by selection in mammals. Furthermore, this result also indicates that there is no reduction of mutation rates in genes expressed in the germ line, contrary to what had been hypothesized based on the fact that transcribed DNA is more efficiently repaired than nontranscribed DNA.     

Type 5 adenylyl cyclase disruption leads to enhanced exercise performance

The most important physiological mechanism mediating enhanced exercise performance is increased sympathetic, beta adrenergic receptor (β‐AR), and adenylyl cyclase (AC) activity. This is the first report of decreased AC activity mediating increased exercise performance. We demonstrated that AC5 disruption, that is, knock out (KO) mice, a longevity model, increases exercise performance. Importantly for its relation to longevity, exercise was also improved in old AC5 KO. The mechanism resided in skeletal muscle rather than in the heart, as confirmed by cardiac‐ and skeletal muscle‐specific AC5 KO's, where exercise performance was no longer improved by the cardiac‐specific AC5 KO, but was by the skeletal muscle‐specific AC5 KO, and there was no difference in cardiac output during exercise in AC5 KO vs. WT. Mitochondrial biogenesis was a major mechanism mediating the enhanced exercise. SIRT1, FoxO3a, MEK, and the anti‐oxidant, MnSOD were upregulated in AC5 KO mice. The improved exercise in the AC5 KO was blocked with either a SIRT1 inhibitor, MEK inhibitor, or by mating the AC5 KO with MnSOD hetero KO mice, confirming the role of SIRT1, MEK, and oxidative stress mechanisms. The Caenorhabditis elegans worm AC5 ortholog, acy‐3 by RNAi, also improved fitness, mitochondrial function, antioxidant defense, and lifespan, attesting to the evolutionary conservation of this pathway. Thus, decreasing sympathetic signaling through loss of AC5 is not only a mechanism to improve exercise performance, but is also a mechanism to improve healthful aging, as exercise also protects against diabetes, obesity, and cardiovascular disease, which all limit healthful aging.