Nonallelic gene conversion is not GC-biased in Drosophila or primates

Gene conversion is the unidirectional transfer of genetic information between allelic (orthologous) or nonallelic (paralogous) DNA segments. Recently, there has been much interest in understanding how gene conversion shapes the nucleotide composition of the genomic landscape. A widely held hypothesis is that gene conversion is universally GC-biased. However, direct experimental evidence of this hypothesis is limited to a single study of meiotic crossovers in yeast. Although there have been a number of indirect studies of gene conversion, evidence of GC-biased replacements gathered from such studies can also be attributed to positive selection, which has the same evolutionary dynamics as biased gene conversion. Here, we apply a direct phylogenetic approach to examine nucleotide replacements produced by nonallelic gene conversion in Drosophila and primate genomes. We find no evidence for GC-biased gene conversion in either lineage, suggesting that previously observed GC biases may be due to positive selection rather than to biased gene conversion.     

Intraphylum diversity and complex evolution of cyanobacterial aminoacyl-tRNA synthetases

A comparative genomic analysis of 35 cyanobacterial strains has revealed that the gene complement of aminoacyl-tRNA synthetases (AARSs) and routes for aminoacyl-tRNA synthesis may differ among the species of this phylum. Several genes encoding AARS paralogues were identified in some genomes. In-depth phylogenetic analysis was done for each of these proteins to gain insight into their evolutionary history. GluRS, HisRS, ArgRS, ThrRS, CysRS, and Glu-Q-RS showed evidence of a complex evolutionary course as indicated by a number of inconsistencies with our reference tree for cyanobacterial phylogeny. In addition to sequence data, support for evolutionary hypotheses involving horizontal gene transfer or gene duplication events was obtained from other observations including biased sequence conservation, the presence of indels (insertions or deletions), or vestigial traces of ancestral redundant genes. We present evidences for a novel protein domain with two putative transmembrane helices recruited independently by distinct AARS in particular cyanobacteria.     

芸薹属A，B和C基因组之间关系研究进展

芸薹属A,B和C基因组之间的亲缘关系近年来取得了很大进展,大量细胞遗传学和分子生物学的研究结果表明A和C基因组之间的亲缘关系较A和B基因组以及B和C基因组之间更为接近。A,B和C基因组之间的比较基因组结果表明,这3个基因组是由更加原始物种进化而来的。在芸薹属基因组演化过程中发生了大量的染色体变异,如重复、缺失、重排等,从而造成了现在不同基因组之间的差别。最后,文章对芸薹属不同基因组和拟南芥基因组之间的亲缘关系进行了综述。     

Context-dependent mutation rates may cause spurious signatures of a fixation bias favoring higher GC-content in humans

Understanding the proximate and ultimate causes underlying the evolution of nucleotide composition in mammalian genomes is of fundamental interest to the study of molecular evolution. Comparative genomics studies have revealed that many more substitutions occur from G and C nucleotides to A and T nucleotides than the reverse, suggesting that mammalian genomes are not at equilibrium for base composition. Analysis of human polymorphism data suggests that mutations that increase GC-content tend to be at much higher frequencies than those that decrease or preserve GC-content when the ancestral allele is inferred via parsimony using the chimpanzee genome. These observations have been interpreted as evidence for a fixation bias in favor of G and C alleles due to either positive natural selection or biased gene conversion. Here, we test the robustness of this interpretation to violations of the parsimony assumption using a data set of 21,488 noncoding single nucleotide polymorphisms (SNPs) discovered by the National Institute of Environmental Health Sciences (NIEHS) SNPs project via direct resequencing of n = 95 individuals. Applying standard nonparametric and parametric population genetic approaches, we replicate the signatures of a fixation bias in favor of G and C alleles when the ancestral base is assumed to be the base found in the chimpanzee outgroup. However, upon taking into account the probability of misidentifying the ancestral state of each SNP using a context-dependent mutation model, the corrected distribution of SNP frequencies for GC-content increasing SNPs are nearly indistinguishable from the patterns observed for other types of mutations, suggesting that the signature of fixation bias is a spurious artifact of the parsimony assumption.     

GC-biased gene conversion impacts ribosomal DNA evolution in vertebrates, angiosperms, and other eukaryotes

Ribosomal DNA (rDNA) is one of the most conserved genes in eukaryotes. The multiples copies of rDNA in the genome evolve in a concerted manner, through unequal crossing over and/or gene conversion, two mechanisms related to homologous recombination. Recombination increases local GC content in several organisms through a process known as GC-biased gene conversion (gBGC). gBGC has been well characterized in mammals, birds, and grasses, but its phylogenetic distribution across the tree of life is poorly understood. Here, we test the hypothesis that recombination affects the evolution of base composition in 18S rDNA and examine the reliability of this thoroughly studied molecule as a marker of gBGC in eukaryotes. Phylogenetic analyses of 18S rDNA in vertebrates and angiosperms reveal significant heterogeneity in the evolution of base composition across both groups. Mammals, birds, and grasses experience increases in the GC content of the 18S rDNA, consistent with previous genome-wide analyses. In addition, we observe increased GC contents in Ostariophysi ray-finned fishes and commelinid monocots (i.e., the clade including grasses), suggesting that the genomes of these two groups have been affected by gBGC. Polymorphism analyses in rDNA confirm that gBGC, not mutation bias, is the most plausible explanation for these patterns. We also find that helix and loop sites of the secondary structure of ribosomal RNA do not evolve at the same pace: loops evolve faster than helices, whereas helices are GC richer than loops. We extend analyses to major lineages of eukaryotes and suggest that gBGC might have also affected base composition in Giardia (Diplomonadina), nudibranch gastropods (Mollusca), and Asterozoa (Echinodermata).     

Can mutation or fixation biases explain the allele frequency distribution of human single nucleotide polymorphisms (SNPs)?

One of the most abiding controversies in evolutionary biology concerns the role of neutral processes in molecular evolution. A main focus of the debate has been the evolution of isochores, the strong and systematic variation of base composition in mammalian genomes. One set of hypotheses argue that regions of similar GC are owing to localised mutational biases coupled with neutral evolution. The alternatives point to either selection or biased gene conversion as mechanisms to preferentially remove A or T bases, favouring G and C instead. Using a novel method, we compare models including such fixation biases to models based on mutation bias alone, under the assumption that non-coding, non-repetitive human DNA is at compositional equilibrium. While failing to fully explain the allele frequency distributions of recent single nucleotide polymorphism data, we show that the data are best fitted if the mutation bias is assumed to be constant across the genome, while fixation bias varies with GC content. We also attempt to estimate the strength of fixation bias, which increases linearly with increasing GC. Our approximation suggests that this force exists within the necessary parameter range: it is not so weak as to be drowned by random drift, but not so strong as to lead to exclusive use of G and C alone. Together these results demonstrate that mutation bias fails to explain the evolution of isochores, and suggest that either selection or biased gene conversion are involved.     

The evolution of chloroplast genes and genomes in ferns

Most of the publicly available data on chloroplast (plastid) genes and genomes come from seed plants, with relatively little information from their sister group, the ferns. Here we describe several broad evolutionary patterns and processes in fern plastid genomes (plastomes), and we include some new plastome sequence data. We review what we know about the evolutionary history of plastome structure across the fern phylogeny and we compare plastome organization and patterns of evolution in ferns to those in seed plants. A large clade of ferns is characterized by a plastome that has been reorganized with respect to the ancestral gene order (a similar order that is ancestral in seed plants). We review the sequence of inversions that gave rise to this organization. We also explore global nucleotide substitution patterns in ferns versus those found in seed plants across plastid genes, and we review the high levels of RNA editing observed in fern plastomes.     

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.     

Physical mapping of rDNA loci in Brassica species.

The number of major rDNA loci (the genes coding for 18S-5.8S-26S rRNA) was investigated in the economically important Brassica species and their wild relatives by in situ hybridization of an rDNA probe to metaphase chromosomes and interphase nuclei. The diploid species B. nigra (B genome) has two major pairs of rDNA loci, B. oleracea (C genome) has two major pairs and one minor pair of loci, while B. campestris (A genome) has five pairs of loci. Among the three tetraploid species arising from these three diploid ancestors, B. carinata (BBCC genomes) has four loci, B. juncea (AABB genomes) has five major pairs and one minor pair of loci, and B. napus (AACC genomes) has six pairs of loci, indicating that the number of loci has been reduced during evolution. The complexity of the known rDNA restriction fragment length polymorphism patterns gave little indication of number of rDNA loci. It is probable that chromosome rearrangements have occurred during evolution of the amphidiploid species. The data will be useful for physical mapping of genes relative to rDNA loci, micro- and macro-evolutionary studies and analysis of aneuploids including addition and substitution lines used in Brassica breeding programs.     

Multilocus patterns of polymorphism and selection across the X chromosome of Caenorhabditis remanei

Natural selection and neutral processes such as demography, mutation, and gene conversion all contribute to patterns of polymorphism within genomes. Identifying the relative importance of these varied components in evolution provides the principal challenge for population genetics. To address this issue in the nematode Caenorhabditis remanei, I sampled nucleotide polymorphism at 40 loci across the X chromosome. The site-frequency spectrum for these loci provides no evidence for population size change, and one locus presents a candidate for linkage to a target of balancing selection. Selection for codon usage bias leads to the non-neutrality of synonymous sites, and despite its weak magnitude of effect (N(e)s approximately 0.1), is responsible for profound patterns of diversity and divergence in the C. remanei genome. Although gene conversion is evident for many loci, biased gene conversion is not identified as a significant evolutionary process in this sample. No consistent association is observed between synonymous-site diversity and linkage-disequilibrium-based estimators of the population recombination parameter, despite theoretical predictions about background selection or widespread genetic hitchhiking, but genetic map-based estimates of recombination are needed to rigorously test for a diversity-recombination relationship. Coalescent simulations also illustrate how a spurious correlation between diversity and linkage-disequilibrium-based estimators of recombination can occur, due in part to the presence of unbiased gene conversion. These results illustrate the influence that subtle natural selection can exert on polymorphism and divergence, in the form of codon usage bias, and demonstrate the potential of C. remanei for detecting natural selection from genomic scans of polymorphism.     

Plastid genome sequences of legumes reveal parallel inversions and multiple losses of rps16 in papilionoids

To date, publicly available plastid genomes of legumes have for the most part been limited to the subfamily Papilionoideae. Here we report 13 new plastid genomes of legumes spanning all three subfamilies. The genomes representing Caesalpinioideae and Mimosoideae are highly conserved in gene content and gene order, similar to the ancestral angiosperm genome organization. Genomes within the Papilionoideae, however, have reduced sizes due to deletions in nine intergenic spacers primarily in the large single copy region. Our study also indicates that rps16 has been independently lost at least five times in legumes, with additional gene and intron losses scattered among the papilionoids. Additionally, genera from two distinct lineages within the papilionoids, Lupinus and Robinia, have a parallel inversion of 36 and 39 kb, respectively. This parallel inversion is novel as it appears to be caused by a 29 bp repeat within two trnS genes. This repeat is present in all available legume plastid genomes indicating that there is the potential for this inversion to be present in more species. This case of a homoplasious inversion is also evidence that some inversion events may not be reliable phylogenetic markers.     

SINEs, evolution and genome structure in the opossum

Short INterspersed Elements (SINEs) are non-autonomous retrotransposons, usually between 100 and 500 base pairs (bp) in length, which are ubiquitous components of eukaryotic genomes. Their activity, distribution, and evolution can be highly informative on genomic structure and evolutionary processes. To determine recent activity, we amplified more than one hundred SINE1 loci in a panel of 43 M. domestica individuals derived from five diverse geographic locations. The SINE1 family has expanded recently enough that many loci were polymorphic, and the SINE1 insertion-based genetic distances among populations reflected geographic distance. Genome-wide comparisons of SINE1 densities and GC content revealed that high SINE1 density is associated with high GC content in a few long and many short spans. Young SINE1s, whether fixed or polymorphic, showed an unbiased GC content preference for insertion, indicating that the GC preference accumulates over long time periods, possibly in periodic bursts. SINE1 evolution is thus broadly similar to human Alu evolution, although it has an independent origin. High GC content adjacent to SINE1s is strongly correlated with bias towards higher AT to GC substitutions and lower GC to AT substitutions. This is consistent with biased gene conversion, and also indicates that like chickens, but unlike eutherian mammals, GC content heterogeneity (isochore structure) is reinforced by substitution processes in the M. domestica genome. Nevertheless, both high and low GC content regions are apparently headed towards lower GC content equilibria, possibly due to a relative shift to lower recombination rates in the recent Monodelphis ancestral lineage. Like eutherians, metatherian (marsupial) mammals have evolved high CpG substitution rates, but this is apparently a convergence in process rather than a shared ancestral state.     

Elimination of deleterious mutations in plastid genomes by gene conversion

Asexual reproduction is believed to be detrimental, mainly because of the accumulation of deleterious mutations over time, a hypothesis known as Muller's ratchet. In seed plants, most asexually reproducing genetic systems are polyploid, with apomictic species (plants forming seeds without fertilization) as well as plastids and mitochondria providing prominent examples. Whether or not polyploidy helps asexual genetic systems to escape Muller's ratchet is unknown. Gene conversion, particularly when slightly biased, represents a potential mechanism that could allow asexual genetic systems to reduce their mutation load in a genome copy number-dependent manner. However, direct experimental evidence for the operation of gene conversion between genome molecules to correct mutations is largely lacking. Here we describe an experimental system based on transgenic tobacco chloroplasts that allows us to analyze gene conversion events in higher plant plastid genomes. We provide evidence for gene conversion acting as a highly efficient mechanism by which the polyploid plastid genetic system can correct deleterious mutations and make one good genome out of two bad ones. Our finding that gene conversion can be biased may provide a molecular link between asexual reproduction, high genome copy numbers and low mutation rates.     

Integrating genomics,bioinformatics, and classical genetics to study the effects of recombination on genome evolution

This study presents compelling evidence that recombination significantly increases the silent GC content of a genome in a selectively neutral manner, resulting in a highly significant positive correlation between recombination and "GC3s" in the yeast Saccharomyces cerevisiae. Neither selection nor mutation can explain this relationship. A highly significant GC-biased mismatch repair system is documented for the first time in any member of the Kingdom Fungi. Much of the variation in the GC3s within yeast appears to result from GC-biased gene conversion. Evidence suggests that GC-biased mismatch repair exists in numerous organisms spanning six kingdoms. This transkingdom GC mismatch repair bias may have evolved in response to a ubiquitous AT mutational bias. A significant positive correlation between recombination and GC content is found in many of these same organisms, suggesting that the processes influencing the evolution of the yeast genome may be a general phenomenon. Nonrecombining regions of the genome and nonrecombining genomes would not be subject to this type of molecular drive. It is suggested that the low GC content characteristic of many nonrecombining genomes may be the result of three processes (1) a prevailing AT mutational bias, (2) random fixation of the most common types of mutation, and (3) the absence of the GC-biased gene conversion which, in recombining organisms, permits the reversal of the most common types of mutation. A model is proposed to explain the observation that introns, intergenic regions, and pseudogenes typically have lower GC content than the silent sites of corresponding open reading frames. This model is based on the observation that the greater the heterology between two sequences, the less likely it is that recombination will occur between them. According to this "Constraint" hypothesis, the formation and propagation of heteroduplex DNA is expected to occur, on average, more frequently within conserved coding and regulatory regions of the genome. In organisms possessing GC-biased mismatch repair, this would enhance the GC content of these regions through biased gene conversion. These findings have a number of important implications for the way we view genome evolution and suggest a new model for the evolution of sex.     

GC-Content Evolution in Bacterial Genomes: The Biased Gene Conversion Hypothesis Expands

The characterization of functional elements in genomes relies on the identification of the footprints of natural selection. In this quest, taking into account neutral evolutionary processes such as mutation and genetic drift is crucial because these forces can generate patterns that may obscure or mimic signatures of selection. In mammals, and probably in many eukaryotes, another such confounding factor called GC-Biased Gene Conversion (gBGC) has been documented. This mechanism generates patterns identical to what is expected under selection for higher GC-content, specifically in highly recombining genomic regions. Recent results have suggested that a mysterious selective force favouring higher GC-content exists in Bacteria but the possibility that it could be gBGC has been excluded. Here, we show that gBGC is probably at work in most if not all bacterial species. First we find a consistent positive relationship between the GC-content of a gene and evidence of intra-genic recombination throughout a broad spectrum of bacterial clades. Second, we show that the evolutionary force responsible for this pattern is acting independently from selection on codon usage, and could potentially interfere with selection in favor of optimal AU-ending codons. A comparison with data from human populations shows that the intensity of gBGC in Bacteria is comparable to what has been reported in mammals. We propose that gBGC is not restricted to sexual Eukaryotes but also widespread among Bacteria and could therefore be an ancestral feature of cellular organisms. We argue that if gBGC occurs in bacteria, it can account for previously unexplained observations, such as the apparent non-equilibrium of base substitution patterns and the heterogeneity of gene composition within bacterial genomes. Because gBGC produces patterns similar to positive selection, it is essential to take this process into account when studying the evolutionary forces at work in bacterial genomes.     

The compositional evolution of the murid genome

Murid rodents show much less variation in isochore base composition than do most other mammals, a difference which has been referred to as the murid shift. We have investigated the murid shift by asking (1) whether the murid shift is ongoing and (2) whether there is any evidence of selection or biased gene conversion affecting base composition in the present-day mouse genome. By estimating the ancestral base composition of protein-coding genes in murids we can confirm that the murid shift is ongoing. Tests using nongenic polymorphism data fail to reject the hypothesis that base composition is due to mutation bias alone. However, the patterns of compositional change suggested by the polymorphism and divergence data differ, suggesting the possibility of two murid shifts.     

A hotspot of gene order rearrangement by tandem duplication and random loss in the vertebrate mitochondrial genome

Most reported examples of change in vertebrate mitochondrial (mt) gene order could be explained by a tandem duplication followed by random loss of redundant genes (tandem duplication-random loss [TDRL] model). Under this model of evolution, independent loss of genes arising from a single duplication in an ancestral species are predicted, and remnant pseudogenes expected, intermediate states that may remain in rearranged genomes. However, evidence for this is rare and largely scattered across vertebrate lineages. Here, we report new derived mt gene orders in the vertebrate "WANCY" region of four closely related caecilian amphibians. The novel arrangements found in this genomic region (one of them is convergent with the derived arrangement of marsupials), presence of pseudogenes, and positions of intergenic spacers fully satisfy predictions from the TDRL model. Our results, together with comparative data for the available vertebrate complete mt genomes, provide further evidence that the WANCY genomic region is a hotspot for gene order rearrangements and support the view that TDRL is the dominant mechanism of gene order rearrangement in vertebrate mt genomes. Convergent gene rearrangements are not unlikely in hotspots of gene order rearrangement by TDRL.     

Poor genetic differentiation but clear cytoform divergence among cryptic species in Simulium damnosum complex (Diptera: Simuliidae)

A complex cryptic species taxonomy has been developed for the flies that vector Onchocerca volvulus L euckart in A frica. This taxonomy names cytoforms of Simulium damnosum T heobald s.l. based on chromosomal inversion polymorphisms. Researchers have linked variation in cytoforms to variation in geographical distributions and biological traits, but have failed to demonstrate that many of these cytoforms are differentiated genetically or morphologically. This may be explained by recent divergence and incomplete lineage sorting of ancestral polymorphisms. Alternatively, it may be the outcome of ongoing gene flow, with inversions maintaining sets of co‐adapted alleles and selection driving inversion frequency differences across habitats. Here, our goals were to test species delimitations in the S. damnosum complex, and determine if the poor differentiation between these species is better explained by the ancestral polymorphisms or gene‐flow hypotheses. Geometric morphometric analysis and more intensive genetic sampling of populations in N igeria only reinforced the lack of clear divergence among S. damnosum complex cytoforms. By contrast, the first statistical estimates of evolutionary relationships based on chromosomal inversion polymorphism data revealed clear signs of divergence, which lends support to the ancestral polymorphisms hypothesis. But the latter analyses may be biased by non‐neutral evolution of inversion polymorphisms; more rigorous tests of cytoform‐based species hypothesis will require deeper sampling across populations and neutral genetic markers along with explicit modelling of genetic isolation and migration.