Evolutionary dynamics of Waxy and the origin of hexaploid Spartina species (Poaceae)

We investigated the evolutionary dynamics of duplicated copies of the granule-bound starch synthase I gene (GBSSI or Waxy) within polyploid Spartina species. Molecular cloning, sequencing, and phylogenetic analyses revealed incongruences between the expected species phylogeny and the inferred gene trees. Some genes within species were more divergent than expected from ploidy level alone, suggesting the existence of paralogous sets of Waxy loci in Spartina. Phylogenetic analyses indicate that this paralogy originated from a duplication that occurred prior to the divergence of Spartina from other Chloridoideae. Gene tree topologies revealed three divergent homoeologous sequences in the hexaploid S. alterniflora that are consistent with the proposal of an allopolyploid origin of the hexaploid clade. Waxy sequences differ in insertion–deletion events in introns, which may be used to diagnose gene copies. Both paralogous and homoeologous coding regions appear to evolving under selective constraints.     

Resolving tylenchid evolutionary relationships through multiple gene analysis derived from EST data

Sequence-based phylogenetic analyses typically are based on a small number of character sets and report gene trees which may not reflect the true species tree. We employed an EST mining strategy to suppress such incongruencies, and recovered the most robust phylogeny for five species of plant-parasitic nematode (Meloidogyne arenaria, M. chitwoodi, M. hapla, M. incognita, and M. javanica), three closely related tylenchid taxa (Heterodera glycines, Globodera pallida, and G. rostochiensis) and a distant taxon, Caenorhabditis elegans. Our multiple-gene approach is based on sampling more than 80,000 publicly available tylenchid EST sequences to identify phylum-wide orthologues. Bayesian inference, minimum evolution, maximum likelihood and protein distance methods were employed for phylogenetic reconstruction and hypothesis tests were constructed to elucidate differential selective pressures across the phylogeny for each gene. Our results place M. incognita and M. javanica as sister taxa, with M. arenaria as the next closely related nematode. Significant differences in selective pressure were revealed for some genes under some hypotheses, though all but one gene are exclusively under purifying selection, indicating conservation across the orthologous groups. This EST-based multi-gene analysis is a first step towards accomplishing genome-wide coverage for tylenchid evolutionary analyses.     

Diversity and evolution of 5S rRNA gene family organization in Pythium

The 5S rRNA gene family organization among 87 species and varieties of Pythium was investigated to assess evolutionary stability of the two patterns detected and to determine which pattern is likely the ancestral state in the genus. Species with filamentous sporangia (Groups A-C according to the ITS phylogenetic tree for Pythium) had 5S genes linked to the rDNA repeat that were predominantly coded for on the DNA strand opposite to the one with the other rRNA genes (‘inverted’ orientation). A small group of species with contiguous sporangia (Group D) is related to Groups A-C but had unlinked 5S genes. The main group of species with spherical zoosporangia (Groups E-J) generally had unlinked 5S genes in tandem arrays. The six species in Group K, although they also have spherical sporangia, had linked genes on the same strand as the other rRNA genes ‘non-inverted’ and most of them had pairs of tandem 5S genes. The evolutionary stability of 5S sequence organization was compared with the stability of morphological characters as interpreted from a phylogeny based on ITS sequence analysis. Features of 5S sequence organization were found to be just as consistent within groups as were the morphological characters. To determine the ancestral type of 5S family organization, a survey of Phytophthora strains was conducted to supply an outgroup reference. The most parsimonious interpretation of the data in this survey yielded the tentative conclusion that the linked condition of the 5S sequences was ancestral.     

Phylogenetic affiliation of Aeromonas culicicola MTCC 3249(T) based on gyrB gene sequence and PCR-amplicon sequence analysis of cytolytic enterotoxin gene

We determined the gyrB gene sequences of all 17 hybridizations groups of Aeromonas. Phylogenetic trees showing the evolutionary relatedness of gyrB and 16S rRNA genes in the type strains of Aeromonas were compared. Using this approach, we determined the phylogenetic position of Aeromonas culicicola MTCC 3249(T), isolated from midgut of Culex quinquefasciatus. In the gyrB based-analysis A. culicicola MTCC 3249(T) grouped with A. veronii whereas, it grouped with A. jandaei in the 16S rRNA based tree. The number of nucleotide differences in 16S rRNA sequences was less than found with the gyrB sequence data. Most of the observed nucleotide differences in the gyrB gene were synonymous. The Cophenetic Correlation Coefficient (CCC) for gyrB sequences was 0.87 indicating this gene to be a better molecular chronometer compared to 16S rRNA for delineation of Aeromonas species. This strain was found to be positive for the cytolytic enterotoxin gene. PCR-Amplicon Sequence Analysis (PCR-ASA) of this gene showed that the isolate is affiliated to type I and is potentially pathogenic. These PCR-ASA results agreed in part with the gyrB sequence results.     

On shortcomings of using mtDNA sequence divergence for the systematics of hares (genus Lepus): An example from cape hares

Despite the well-known fact that evolutionary patterns of single genes or sequences are not necessarily paralleled by organismic evolution, using mitochondrial sequence divergence for inferring phylogenetic relationships among hare species is not uncommon. Earlier studies reported interspecific introgression in the genus Lepus that may slur systematic conclusions drawn exclusively from mtDNA data. We examined pairwise divergence in partial mtHV-1 sequences and nuclear gene pools based on microsatellite and allozyme loci separately in a South and North African population of cape hares, Lepus capensis, and did not find significant correspondence on the individual level within either population. Also, the former population had a significantly higher level of mtHV-1 sequence divergence compared to the latter, but levels of individual nuclear gene pool divergence did not differ significantly between the two populations. Hence, the absence of correspondence between mtHV-1 sequence divergence and nuclear gene pool divergence among individuals within a population might be independent of the population-specific level of differentiation among individuals. Our results complement earlier findings in hares, and strongly recommend to include nuclear gene pool evidence for systematic inferences within this genus, even under the absence of mitochondrial introgression.     

Molecular Evolution at the Cytochrome Oxidase Subunit 2 Gene Among Divergent Populations of the Intertidal Copepod, Tigriopus californicus

The cytochrome c oxidase subunit 2 gene (COII) encodes a highly conserved protein that is directly responsible for the initial transfer of electrons from cytochrome c to cytochrome c oxidase (COX) crucial to the production of ATP during cellular respiration. Despite its integral role in electron transport, we have observed extensive intraspecific nucleotide and amino acid variation among 26 full-length COII sequences sampled from seven populations of the marine copepod, Tigriopus californicus. Although intrapopulation divergence was virtually nonexistent, interpopulation divergence at the COII locus was nearly 20% at the nucleotide level, including 38 nonsynonymous substitutions. Given the high degree of interaction between the cytochrome c oxidase subunit 2 protein (COX2) and the nuclear-encoded subunits of COX and cytochrome c (CYC), we hypothesized that some codons in the COII gene are likely to be under positive selection in order to compensate for amino acid substitutions in other subunits. Estimates of the ratio of nonsynonymous to synonymous substitution (ω), obtained using a series of maximum likelihood models of codon substitution, indicated that the majority of codons in T. californicus COII are under strong purifying selection (ω << 1), while approximately 4% of the sites in this gene appear to evolve under relaxed selective constraint (ω = 1). A branch-site maximum likelihood model identified three sites that may have experienced positive selection within the central California sequence clade in our COII phylogeny; these results are consistent with previous studies showing functional and fitness consequences among interpopulation hybrids between central and northern California populations. [Reviewing Editor: Dr. Willie Swanson]     

Interspecific divergence, intrachromosomal recombination, and phylogenetic utility of Y-chromosomal genes in Drosophila

Reconstruction of phylogenetic relationships among recently diverged species is complicated by three general problems: segregation of polymorphisms that pre-date species divergence, gene flow during and after speciation, and intra-locus recombination. In light of these difficulties, the Y chromosome offers several important advantages over other genomic regions as a source of phylogenetic information. These advantages include the absence of recombination, rapid coalescence, and reduced opportunity for interspecific introgression due to hybrid male sterility. In this report, we test the phylogenetic utility of Y-chromosomal sequences in two groups of closely related and partially inter-fertile Drosophila species. In the D. bipectinata species complex, Y-chromosomal loci unambiguously recover the phylogeny most consistent with previous multi-locus analysis and with reproductive relationships, and show no evidence of either post-speciation gene flow or persisting ancestral polymorphisms. In the D. simulans species complex, the situation is complicated by the duplication of at least one Y-linked gene region, followed by intrachromosomal recombination between the duplicate genes that scrambles their genealogy. We suggest that Y-chromosomal sequences are a useful tool for resolving phylogenetic relationships among recently diverged species, especially in male-heterogametic organisms that conform to Haldane's rule. However, duplication of Y-linked genes may not be uncommon, and special care should be taken to distinguish between orthologous and paralogous sequences.     

Comparative Sequence Analysis of the <Emphasis Type="Italic">Phytochrome C</Emphasis> Gene and its Upstream Region in Allohexaploid Wheat Reveals New Data on the Evolution of its Three Constituent Genomes

Bread wheat is an allohexaploid with genome composition AABBDD. Phytochrome C is a gene involved in photomorphogenesis that has been used extensively for phylogenetic analyses. In wheat, the PhyC genes are single copy in each of the three homoeologous genomes and map to orthologous positions on the long arms of the group 5 chromosomes. Comparative sequence analysis of the three homoeologous copies of the wheat PhyC gene and of some 5 kb of upstream region has demonstrated a high level of conservation of PhyC, but frequent interruption of the upstream regions by the insertion of retroelements and other repeats. One of the repeats in the region under investigation appeared to have inserted before the divergence of the diploid wheat genomes, but was degraded to the extent that similarity between the A and D copies could only be observed at the amino acid level. Evidence was found for the differential presence of a foldback element and a miniature inverted-repeat transposable element (MITE) 5′ to PhyC in different wheat cultivars. The latter may represent the first example of an active MITE family in the wheat genome. Several conserved non-coding sequences were also identified that may represent functional regulatory elements. The level of sequence divergence (K_s) between the three wheat PhyC homoeologs suggests that the divergence of the diploid wheat ancestors occurred some 6.9 Mya, which is considerably earlier than the previously estimated 2.5–4.5 Mya. K_a/K_s ratios were <0.15 indicating that all three homoeologs are under purifying selection and presumably represent functional PhyC genes. RT-PCR confirmed expression of the A, B and D copies. The discrepancy in evolutionary age of the wheat genomes estimated using sequences from different parts of the genome may reflect a mosaic origin of some of the Triticeae genomes.     

Molecular evolution of the nontandemly repeated genes of the histone 3 multigene family.

In some species, histone gene clusters consist of tandem arrays of each type of histone gene, whereas in other species the genes may be clustered but not arranged in tandem. In certain species, however, histone genes are found scattered across several different chromosomes. This study examines the evolution of histone 3 (H3) genes that are not arranged in large clusters of tandem repeats. Although H3 amino acid sequences are highly conserved both within and between species, we found that the nucleotide sequence divergence at synonymous sites is high, indicating that purifying selection is the major force for maintaining H3 amino acid sequence homogeneity over long-term evolution. In cases where synonymous-site divergence was low, recent gene duplication appeared to be a better explanation than gene conversion. These results, and other observations on gene inactivation, organization, and phylogeny, indicated that these H3 genes evolve according to a birth-and-death process under strong purifying selection. Thus, we found little evidence to support previous claims that all H3 proteins, regardless of their genome organization, undergo concerted evolution. Further analyses of the structure of H3 proteins revealed that the histones of higher eukaryotes might have evolved from a replication-independent-like H3 gene.     

The genome of <Emphasis Type="Italic">Streptomyces rimosus</Emphasis> subsp. <Emphasis Type="Italic">rimosus</Emphasis> shows a novel structure compared to other <Emphasis Type="Italic">Streptomyces</Emphasis> using DNA/DNA microarray analysis

DNA/DNA genome microarray analysis together with genome sequencing suggests that the genome of members of the genus Streptomyces would seem to have a common structure including a linear genomic structure, a core of common syntenous Actinomycete genes, the presence of species specific terminal regions and two intermediate group of syntenous genes that seem to be genus specific. We analyzed Streptomyces species using DNA/DNA microarray comparative genome analysis. Only Streptomyces rimosus failed to give a congruent genome pattern for the genes found in Streptomyces coelicolor. We expanded the analysis to include a number of strains related to the type strain of S. rimosus and obtained a similar divergence from the main body of Streptomyces species. These strains showed very close identity to the original strain with no gene deletion or duplication detected. The 16S rRNA sequences of these S. rimosus strains were confirmed as very similar to the S. rimosus sequences available from the Ribosomal Database Project. When the SSU ribosomal RNA phylogeny of S. rimosus is analyzed, the species is positioned at the edge of the Streptomyces clade. We conclude that S. rimosus represents a distinct evolutionary lineage making the species a worthy possibility for genome sequencing.     

Speciation in <Emphasis Type="Italic">Chlamydia</Emphasis>: Genomewide Phylogenetic Analyses Identified a Reliable Set of Acquired Genes

Horizontal gene transfer (HGT), a process through which genomes acquire sequences from distantly related organisms, is believed to be a major source of genetic diversity in bacteria. A central question concerning the impact of HGT on bacterial genome evolution is the proportion of horizontally transferred sequences within genomes. This issue, however, remains unresolved because the various methods developed to detect potential HGT events identify different sets of genes. The present-day consensus is that phylogenetic analysis of individual genes is still the most objective and accurate approach for determining the occurrence and directionality of HGT. Here we present a genome-scale phylogenetic analysis of protein-encoding genes from five closely related Chlamydia, identifying a reliable set of sequences that have arisen via HGT since the divergence of the Chlamydia lineage. According to our knowledge, this is the first systematic phylogenetic inference-based attempt to establish a reliable set of acquired genes in a bacterial genome. Although Chlamydia are obligate intracellular parasites of higher eukaryotes, and thus suspected to be isolated from HGT more than the free-living species, our results show that their diversification has involved the introduction of foreign sequences into their genome. Furthermore, we also identified a complete set of genes that have undergone deletion, duplication, or rearrangement during this evolutionary period leading to the radiation of Chlamydia species. Our analysis may provide a deeper insight into how these medically important pathogens emerged and evolved from a common ancestor.     

Evolution of ketosynthase domains of polyketide synthase genes in the Cladonia chlorophaea species complex (Cladoniaceae)

Lichen-forming fungi synthesize a diversity of polyketides, but only a few non-reducing polyketide synthase (PKS) genes from a lichen-forming fungus have been linked with a specific polyketide. While it is a challenge to link the large number of PKS paralogs in fungi with specific products, it might be expected that the PKS paralogs from closely related species would be similar because of recent evolutionary divergence. The objectives of this study were to reconstruct a PKS gene phylogeny of the Cladonia chlorophaea species complex based on the ketosynthase domain, a species phylogeny of the complex, and to explore the presence of PKS gene paralogs among members of the species complex. DNA was isolated from 51 individuals of C. chlorophaea and allies to screen for the presence of 13 PKS paralogs. A 128 sequence PKS gene phylogeny using deduced amino acid sequences estimated from the 13 PKS paralogs and sequences subjected to BLASTx comparisons showed losses of each of two PKS domains (reducing and methylation). This research provided insight into the evolution of PKS genes in the C. chlorophaea group, species evolution in the group, and it identified potential directions for further investigation of polyketide synthesis in the C. chlorophaea species complex.     

High mitochondrial diversity in geographically widespread butterflies of Madagascar: A test of the DNA barcoding approach

The standardized use of mitochondrial cytochrome c oxidase subunit I (COI) gene sequences as DNA barcodes has been widely promoted as a high-throughput method for species identification and discovery. Species delimitation has been based on the following criteria: (1) monophyletic association and less frequently (2) a minimum 10× greater divergence between than within species. Divergence estimates, however, can be inflated if sister species pairs are not included and the geographic extent of variation within any given taxon is not sampled comprehensively. This paper addresses both potential biases in DNA divergence estimation by sampling range-wide variation in several morphologically distinct, endemic butterfly species in the genus Heteropsis, some of which are sister taxa. We also explored the extent to which mitochondrial DNA from the barcode region can be used to assess the effects of historical rainforest fragmentation by comparing genetic variation across Heteropsis populations with an unrelated forest-associated taxon Saribia tepahi. Unexpectedly, generalized primers led to the inadvertent amplification of the endosymbiont Wolbachia, undermining the use of universal primers and necessitating the design of genus-specific COI primers alongside a Wolbachia-specific PCR assay. Regardless of the high intra-specific genetic variation observed, most species satisfy DNA barcoding criteria and can be differentiated in the nuclear phylogeny. Nevertheless, two morphologically distinguishable candidate species fail to satisfy the barcoding 10× genetic distance criterion, underlining the difficulties of applying a standard distance threshold to species delimitation. Phylogeographic analysis of COI data suggests that forest fragmentation may have played an important role in the recent evolutionary diversification of these butterflies. Further work on other Malagasy taxa using both mitochondrial and nuclear data will provide better insight into the role of historical habitat fragmentation in species diversification and may potentially contribute to the identification of priority areas for conservation.     

Adaptive Evolution and Recombination of Rickettsia Antigens

The genus Rickettsia consists of intracellular bacteria that cause a variety of arthropod vectored human diseases. I have examined the evolutionary processes that are generating variation in antigens that are potential vaccine candidates. The surface proteins rOmpA and rOmpB are subject to intense positive natural selection, causing rapid diversification of their amino acid sequences between species. The positively selected amino acids were mapped and cluster together in regions that may indicate the location of functionally important regions such as epitopes. In contrast to the rOmp antigens, there is no evidence of positive selection on the intracytoplasmic antigen PS120 despite low selective constraints on this gene. All three genes showed evidence of recombination between species, and certain sequences are clear chimeras of two parental sequences. However, recombination has been sufficiently infrequent that the phylogenies of the three genes are similar, although not identical. [Reviewing Editor: Dr. Willie J. Swanson]     

Evolution of the gibbon subgenera inferred from cytochrome b DNA sequence data

DNA sequences for the mitochondrial cytochrome b gene from the four extant gibbon subgenera are described. The data confirm that the gibbon subgenera evolved from a common hylobatid ancestor and suggest that they diverged from each other after the divergence of the extant African great ape species. The cytochrome b gene does not resolve the evolutionary relationships between the gibbon subgenera themselves.     

Gene duplication,allelic diversity,selection processes and adaptive value of MHC class II <Emphasis Type="Italic">DRB</Emphasis> genes of the bank vole, <Emphasis Type="Italic">Clethrionomys glareolus</Emphasis>

The generation and maintenance of allelic polymorphism in genes of the major histocompatibility complex (MHC) is a central issue in evolutionary genetics. Recently, the focus has changed from ex situ to in situ populations to understand the mechanisms that determine adaptive MHC polymorphism under natural selection. Birth-and-death evolution and gene conversion events are considered to generate sequence diversity in MHC genes, which subsequently is maintained by balancing selection through parasites. The ongoing arms race between the host and parasites leads to an adaptive selection pressure upon the MHC, evident in high rates of non-synonymous vs synonymous substitution rates. We characterised the MHC class II DRB exon 2 of free living bank voles, Clethrionomys glareolus by single-strand conformation polymorphism and direct sequencing. Unlike other arvicolid species, the DRB locus of the bank vole is at least quadruplicated. No evidence for gene conversion events in the Clgl-DRB sequences was observed. We found not only high allelic polymorphism with 26 alleles in 36 individuals but also high rates of silent polymorphism. Exceptional for MHC class II genes is a purifying selection pressure upon the majority of MHC-DRB sequences. Further, we analysed the association between certain DRB alleles and the parasite burden with gastrointestinal trichostrongyle nematodes Heligmosomum mixtum and Heligmosomoides glareoli and found significant quality differences between specific alleles with respect to infection intensity. Our findings suggest a snapshot in an evolutionary process of ongoing birth-and-death evolution. One allele cluster has lost its function and is already silenced, another is loosing its adaptive value in terms of gastrointestinal nematode resistance, while a third group of alleles indicates all signs of classical functional MHC alleles.     

Molecular evolution of growth hormone gene family in old world monkeys and hominoids

Growth hormone is a classic molecule in the study of the molecular clock hypothesis as it exhibits a relatively constant rate of evolution in most mammalian orders except primates and artiodactyls, where dramatically enhanced rate of evolution (25–50-fold) has been reported. The rapid evolution of primate growth hormone occurred after the divergence of tarsiers and simians, but before the separation of old world monkeys (OWM) from new world monkeys (NWM). Interestingly, this event of rapid sequence evolution coincided with multiple duplications of the growth hormone gene, suggesting gene duplication as a possible cause of the accelerated sequence evolution. Here we determined 21 different GH-like sequences from four species of OWM and hominoids. Combining with published sequences from OWM and hominoids, our analysis demonstrates that multiple gene duplications and several gene conversion events both occurred in the evolutionary history of this gene family in OWM/hominoids. The episode of recent duplications of CSH-like genes in gibbon is accompanied with rapid sequence evolution likely resulting from relaxation of purifying selection. GHN genes in both hominoids and OWM are under strong purifying selection. In contrast, CSH genes in both lineages are probably not. GHV genes in OWM and hominoids evolved at different evolutionary rates and underwent different selective constraints. Our results disclosed the complex history of the primate growth hormone gene family and raised intriguing questions on the consequences of these evolutionary events.