Genomic Conflict Settled in Favour of the Species Rather Than the Gene at Extreme GC Percentage Values

Wada and colleagues have shown that, whether prokaryotic or eukaryotic, each gene has a "homostabilising propensity" to adopt a relatively uniform GC percentage (GC%). Accordingly, each gene can be viewed as a "microisochore" occupying a discrete GC% niche of relatively uniform base composition amongst its fellow genes. Although first, second and third codon positions usually differ in GC%, each position tends to maintain a uniform, gene-specific GC% value. Thus, within a genome, genic GC% values can cover a wide range. This is most evident at third codon positions, which are least constrained by amino acid encoding needs. In 1991, Wada and colleagues further noted that, within a phylogenetic group, genomic GC% values can also cover a wide range. This is again most evident at third codon positions. Thus, the dispersion of GC% values among genes within a genome matches the dispersion of GC% values among genomes within a phylogenetic group. Wada described the context-independence of plots of different codon position GC% values against total GC% as a "universal" characteristic. Several studies relate this to recombination. We have confirmed that third codon positions usually relate more to the genes that contain them than to the species. However, in genomes with extreme GC% values (low or high), third codon positions tend to maintain a constant GC%, thus relating more to the species than to the genes that contain them. Genes in an extreme-GC% genome collectively span a smaller GC% range, and mainly rely on first and second codon positions for differentiation as "microisochores". Our results are consistent with the view that differences in GC% serve to recombinationally isolate both genome sectors (facilitating gene duplication) and genomes (facilitating genome duplication, e.g. speciation). In intermediate-GC% genomes, conflict between the needs of the species and the needs of individual genes within that species is minimal. However, in extreme-GC% genomes there is a conflict, which is settled in favour of the species (i.e. group selection) rather than in favour of the gene (genic selection).     

Metabolic rate and genomic GC: what we can learn from teleost fish

Teleosts are a highly diverse group of animals occupying all kind of aquatic environment. Data on routine mass specific metabolic rate were re-examined correcting them for the Boltzmann's factor. Teleostean fish were grouped in five broad groups, corresponding to major environmental classifications: polar, temperate, sub-tropical, tropical and deep-water. The specific routine metabolic rate, temperature-corrected using the Boltzmann's factor (MR), and the average base composition of genomes (GC%) were calculated in each group. Fish of the polar habitat showed the highest MR. Temperate fish displayed a significantly higher MR than tropical fish, which had the lowest average value. These results were apparently in agreement with the cold adaptation hypothesis. In contrast with this hypothesis, however, the MR of fish living in deep-water environment turned out to be not significantly different from that of fish living in tropical habitats. Most probably, the amount of oxygen dissolved in the water directly affects MR adaptation. Regarding the different habitats, the genomic GC levels showed a decreasing trend similar to that of MR. Indeed, both polar and temperate fish showed a GC level significantly higher than that of both sub-tropical and tropical fish. Plotting the genomic GC levels versus the MR a significant positive correlation was found, supporting the hypothesis that metabolic rate can explain not only the compositional transition mode (e.g. amphibian/mammals), but also the compositional shifting mode (e.g. fish/fish) of evolution observed for vertebrate genomes.     

Retrotransposon Sequence Variation in Four Asexual Plant Species

Transposable elements (TEs) can be viewed as genetic parasites that persist in populations due to their capacity for increase in copy number and the inefficacy of selection against them. A corollary of this hypothesis is that TEs are more likely to spread within sexual populations and be eliminated or inactivated within asexual populations. While previous work with animals has shown that asexual taxa may contain less TE diversity than sexual taxa, comparable work with plants has been lacking. Here we report the results of a study of Ty1/copia, Ty3/gypsy, and LINE-like retroelement diversity in four asexual plant species. Retroelement-like sequences, with a high degree of conservation both within and between species, were isolated from all four species. The sequences correspond to several previously annotated retroelement subfamilies. They also exhibit a pattern of nucleotide substitution characterized by an excess of synonymous substitutions, suggestive of a history of purifying selection. These findings were compared with retroelement sequence evolution in sexual plant taxa. One likely explanation for the discovery of conserved TE sequences in the genomes of these asexual taxa is simply that asexuality within these taxa evolved relatively recently, such that the loss and breakdown of TEs is not yet detectable through analysis of sequence diversity. This explanation is examined by conducting stochastic simulation of TE evolution and by using published information to infer rough estimates of the ages of asexual taxa. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Dmitri Petrov     

Molecular hyperdiversity and evolution in very large populations

The genomic density of sequence polymorphisms critically affects the sensitivity of inferences about ongoing sequence evolution, function and demographic history. Most animal and plant genomes have relatively low densities of polymorphisms, but some species are hyperdiverse with neutral nucleotide heterozygosity exceeding 5%. Eukaryotes with extremely large populations, mimicking bacterial and viral populations, present novel opportunities for studying molecular evolution in sexually reproducing taxa with complex development. In particular, hyperdiverse species can help answer controversial questions about the evolution of genome complexity, the limits of natural selection, modes of adaptation and subtleties of the mutation process. However, such systems have some inherent complications and here we identify topics in need of theoretical developments. Close relatives of the model organisms Caenorhabditis elegans and Drosophila melanogaster provide known examples of hyperdiverse eukaryotes, encouraging functional dissection of resulting molecular evolutionary patterns. We recommend how best to exploit hyperdiverse populations for analysis, for example, in quantifying the impact of noncrossover recombination in genomes and for determining the identity and micro‐evolutionary selective pressures on noncoding regulatory elements.     

A human-specific subfamily of Alu sequences.

Of a total of 500,000 Alu family members, approximately 500 are present as a human-specific (HS) subfamily. Each of the HS subfamily members shares a high degree of nucleotide identity and is not present at orthologous positions in other primate genomes, suggesting that HS subfamily members have recently inserted within the human genome. This confirms the hypothesis that the majority of Alu family members are amplified copies of a "master" gene(s). This master gene appears to be amplifying at a rate much slower than that seen earlier in primate evolution. Some of the HS Alu subfamily members have amplified so recently that they are dimorphic in the human population, making them a potentially powerful tool for studies of human populations.     

Genome size, GC percentage and 5mC level in the Indonesian coelacanth Latimeria menadoensis

The living fossil Latimeria menadoensis is important to understand sarcopterygian evolution. To gain further insights into this fish species we studied its genome size, GC% and 5mC level. The genome size and the GC% of the Indonesian coelacanth seem to be very similar to those of the African coelacanth. Moreover the GC%, the CpG frequency and the 5mC level of L. menadoensis are more similar to those of fish and amphibians than to those of mammals, birds and reptiles and this is in line with the hypothesis that two different DNA methylation and CpG shortage equilibria arose during vertebrate evolution. Our results suggest that the genome of L. menadoensis has remained unchanged for several million years, maybe since the origin of the lineage which from lobe-finned fish led to tetrapods. These data fit a conservative evolutionary landscape and suggest that the genome of the extant crossopterygians may be a sort of evolutionarily frozen genome.     

Natural variation, functional divergence, and local adaptation of nucleotide binding site sequences in Rhododendron (Ericaceae)

To further understand natural variation and local adaptation in the evolution of plant defense, we analyzed polymorphism data of nucleotide-binding site (NBS) sequences of Rhododendron at both the species and population levels. Multiple duplication events were found in NBS sequence evolution in Rhododendron genomes, which resulted in six clades: A–F. Our results of several NBS clade pair comparisons showed significant evolutionary rate changes based on differences in substitution rates between NBS-encoding protein clades (type I functional divergence). Pairwise comparisons of NBS clades further revealed that many amino acids displayed radical biochemical property changes causing a shift in amino acid preferences between NBS-encoding protein clades (type II functional divergence). Such divergent evolution of NBSs is likely a consequence of positive selection related to differentiation of recognition signals in response to different pathogens. Primers specific to clades B and C, which differed in the number of radical amino acid changes causing type II functional divergence and levels of nucleotide diversities, were further used to amplify population clades B and C NBS sequences of Rhododendron formosanum populations. Higher levels of net nucleotide divergences (measured by D _a) between R. formosanum populations were found based on NBS sequences of population clade B compared to population clade C, suggesting local adaptation of population clade B NBS sequences. Local adaptation can be further inferred for R. formosanum population clade B NBS sequences because of significant Φ _ST based on variation in nonsynonymous substitutions. Furthermore, local adaptation was also suggested by no significant correlation of population pairwise F _ST between population clades B and C in R. formosanum.     

InSatDb: a genomic tool for insect geneticists

Microsatellites show tremendous variation between genomes in terms of their occurrence and composition. Availability of whole genome sequences allows us to study microsatellite characteristics of fully sequenced insect genomes to understand the evolution and biological significance of microsatellites. InSatDb is an insect microsatellite database that provides an interactive interface to query information on microsatellites annotated with size (in base pairs and repeat units), genomic location (exon, intron, up-stream or transposon), nature (perfect or imperfect), and sequence composition (repeat motif and GC%). Here we present a snapshot of the distribution and composition of microsatellites in introns and exons of insect genomes. The data present interesting observations regarding the microsatellite life-cycle and genome flux.     

InSatDb: A Genomic Tool for Insect Geneticists

Microsatellites show tremendous variation between genomes in terms of their occurrence and composition. Availability of whole genome sequences allows us to study microsatellite characteristics of fully sequenced insect genomes to understand the evolution and biological significance of microsatellites. InSatDb is an insect microsatellite database that provides an interactive interface to query information on microsatellites annotated with size (in base pairs and repeat units); genomic location (exon, intron, up-stream or transposon); nature (perfect or imperfect); and sequence composition (repeat motif and GC%). Here, we present a snap shot of the distribution and composition of microsatellites in introns and exons of insect genomes. The data present interesting observations regarding the microsatellite life-cycle and genome flux.     

Plant molecular diversity and applications to genomics

Surveys of nucleotide diversity are beginning to show how genomes have been shaped by evolution. Nucleotide diversity is also being used to discover the function of genes through the mapping of quantitative trait loci (QTL) in structured populations, the positional cloning of strong QTL, and association mapping.     

利用小RNA深度测序技术鉴定半夏病毒病病原

RNA沉默是真核生物体内由病毒来源的干扰小RNA（virus derived small interfering RNA, vsiRNA）沉默复合物介导目标RNA特异降解的一种保守机制,通过对vsiRNA分析可进行植物病毒病原鉴定。本文利用小RNA深度测序技术对感病半夏叶片进行鉴定,结果发现,表现典型花叶症状的半夏叶片受到大豆花叶病毒（Soybean mosaic virus, SMV）、黄瓜花叶病毒（Cucumber mosaic virus, CMV）、芋花叶病毒（Dasheen mosaic virus, DsMV）、魔芋花叶病毒（Konjac mosaic virus, KoMV）、烟草花叶病毒（Tobacco mosaic virus, TMV）等多种病毒的复合侵染。为明确SMV山西半夏分离物（SMV-SXBX）的进化关系,进行SMV-SXBX全基因组克隆与分析,获得SMV-SXBX全长为9 735 nt,编码一个由3 105个氨基酸组成的多聚蛋白质。通过核苷酸与氨基酸序列比对发现,SMV-SXBX与半夏分离物P同源性最高,分别为91.1%和94.1%,且系统发育分析表明,SMV-SXBX与半夏SMV分离物P聚为一簇。同时,也对vsiRNA进行了系统分析,研究结果有望为半夏SMV的有效防治提供一定科学依据。     

The complete chloroplast and mitochondrial genome sequences of Boea hygrometrica: insights into the evolution of plant organellar genomes

The complete nucleotide sequences of the chloroplast (cp) and mitochondrial (mt) genomes of resurrection plant Boea hygrometrica (Bh, Gesneriaceae) have been determined with the lengths of 153,493 bp and 510,519 bp, respectively. The smaller chloroplast genome contains more genes (147) with a 72% coding sequence, and the larger mitochondrial genome have less genes (65) with a coding faction of 12%. Similar to other seed plants, the Bh cp genome has a typical quadripartite organization with a conserved gene in each region. The Bh mt genome has three recombinant sequence repeats of 222 bp, 843 bp, and 1474 bp in length, which divide the genome into a single master circle (MC) and four isomeric molecules. Compared to other angiosperms, one remarkable feature of the Bh mt genome is the frequent transfer of genetic material from the cp genome during recent Bh evolution. We also analyzed organellar genome evolution in general regarding genome features as well as compositional dynamics of sequence and gene structure/organization, providing clues for the understanding of the evolution of organellar genomes in plants. The cp-derived sequences including tRNAs found in angiosperm mt genomes support the conclusion that frequent gene transfer events may have begun early in the land plant lineage.     

Avirulence gene and insertion element-based RFLP as well as RAPD markers reveal high levels of genomic polymorphism in the rice pathogen Xanthomonas oryzae pv. oryzae

Genetic polymorphism within the genomes of bacterial pathogens determines their evolutionary potential during long-term interaction with their hosts. To investigate the level of genetic variation in Xanthomonas oryzae pv. oryzae (Xoo), the causative agent of rice bacterial blight disease, three DNA marker systems, including (i) restriction fragment length polymorphism (RFLP) of the avrBs3/PthA family genes (avrXa27), (ii) RFLP of insertion (IS) elements and (iii) random amplified polymorphic DNA (RAPD) markers, were used to detect polymorphism among 32 Xoo strains that differed in their virulence patterns. All these strains contained multiple avrXa27 homologs that were variable in copy number and genomic location. RFLP of six IS elements revealed that these mobile sequences were abundant in Xoo genomes, with 150 of the total of 165 discernable markers being variable. Thirty-eight decamer primers of RAPD amplified a total of 691 bands, with 100% of them being variable. In addition, analysis of molecular variance (AMOVA) of data from RFLP analysis of IS elements and from RAPD analysis showed that most of the genetic variation residues were within Xoo populations, rather than between populations. Although all three DNA marker systems supported that substantial variation was maintained in Xoo genomes, Mantel tests did not identify significant correlation between the similarity coefficients calculated from them. The results of the present study indicated that Xoo genomes contain a high level of genetic polymorphism, which greatly facilitates the evolution of this important pathogen during interaction with its host rice plant.