Chih-Horng Kuo

Evolutionary Relationships of Wild Hominids Recapitulated by Gut Microbial Communities

Although bacteria are continually acquired over the lifetime of an individual, the phylogenetic relationships of great ape species is mirrored in the compositions of their gut microbial communities.

Genetic Effects of a Persistent Bottleneck on a Natural Population of Ornate Box Turtles (Terrapene ornata)

Human activities in the past few hundred years have caused enormous impacts on many ecosystems, greatly accelerating the rate of population decline and extinction. In addition to habitat alteration and destruction, the loss of genetic diversity due to reduced population size has become a major conservation issue for many imperiled species. However, the genetic effects of persistent population bottlenecks can be very different for long-lived and short-lived species when considering the time scale of centuries. To investigate the genetic effects of persistent population bottlenecks on long-lived species, we use microsatellite markers to assess the level of genetic diversity of a small ornate box turtle population that has experienced a persistent bottleneck in the past century, and compare it to a large relatively undisturbed population. The genetic signature of a recent bottleneck is detected by examining the deviation from mutation-drift equilibrium in the small population, but the bottleneck had little effect on its level of genetic diversity. Computer simulations combined with information on population structure suggest that an effective population size of 300, which results in a census population size of 700, would be required for the small population to maintain 90% of the average number of alleles per locus in the next 200 years. The life history of long-lived species could mask the accelerated rate of genetic drift, making population recovery a relatively slow process. Statistical analysis of genetic data and empirical-based computer simulations can be important tools to facilitate conservation planning.

Deletional Bias across the Three Domains of Life

Elevated levels of genetic drift are hypothesized to be a dominant factor that influences genome size evolution across all life-forms. However, increased levels of drift appear to be correlated with genome expansion in eukaryotes but with genome contraction in bacteria, suggesting that these two groups of organisms experience vastly different mutational inputs and selective constraints. To determine the contribution of small insertion and deletion events to the differences in genome organization between eukaryotes and prokaryotes, we systematically surveyed 17 taxonomic groups across the three domains of life. Based on over 5,000 indel events in noncoding regions, we found that deletional events outnumbered insertions in all groups examined. The extent of deletional bias, when measured by the total length of insertions to deletions, revealed a marked disparity between eukaryotes and prokaryotes, whereas the ratio was close to one in the three eukaryotic groups examined, deletions outweighed insertions by at least a factor of 10 in most prokaryotes. Moreover, the strength of deletional bias is associated with the proportion of coding regions in prokaryotic genomes. Considering that genetic drift is a stochastic process and does not discriminate the exact nature of mutations, the degree of bias toward deletions provides an explanation to the differential responses of eukaryotes and prokaryotes to elevated levels of drift. Furthermore, deletional bias, rather than natural selection, is the primary mechanism by which the compact gene packing within most prokaryotic genomes is maintained.

The Apicomplexan Whole-Genome Phylogeny: An Analysis of Incongruence among Gene Trees

The protistan phylum Apicomplexa contains many important pathogens and is the subject of intense genome sequencing efforts. Based upon the genome sequences from seven apicomplexan species and a ciliate outgroup, we identified 268 single-copy genes suitable for phylogenetic inference. Both concatenation and consensus approaches inferred the same species tree topology. This topology is consistent with most prior conceptions of apicomplexan evolution based upon ultrastructural and developmental characters, that is, the piroplasm genera Theileria and Babesia form the sister group to the Plasmodium species, the coccidian genera Eimeria and Toxoplasma are monophyletic and are the sister group to the Plasmodium species and piroplasm genera, and Cryptosporidium forms the sister group to the above mentioned with the ciliate Tetrahymena as the outgroup. The level of incongruence among gene trees appears to be high at first glance; only 19% of the genes support the species tree, and a total of 48 different gene-tree topologies are observed. Detailed investigations suggest that the low signal-to-noise ratio in many genes may be the main source of incongruence. The probability of being consistent with the species tree increases as a function of the minimum bootstrap support observed at tree nodes for a given gene tree. Moreover, gene sequences that generate high bootstrap support are robust to the changes in alignment parameters or phylogenetic method used. However, caution should be taken in that some genes can infer a “wrong” tree with strong support because of paralogy, model violations, or other causes. The importance of examining multiple, unlinked genes that possess a strong phylogenetic signal cannot be overstated.

Phylogeographic breaks in low-dispersal species: the emergence of concordance across gene trees

Computer simulations were used to investigate population conditions under which phylogeographic breaks in gene genealogies can be interpreted with confidence to infer the existence and location of historical barriers to gene flow in continuously distributed, low-dispersal species. We generated collections of haplotypic gene trees under a variety of demographic scenarios and analyzed them with regard to salient genealogical breaks in their spatial patterns. In the first part of the analysis, we estimated the frequency in which the spatial location of the deepest phylogeographic break between successive pairs of populations along a linear habitat coincided with a spatial physical barrier to dispersal. Results confirm previous reports that individual gene trees can show ‘haphazard’ phylogeographic discontinuities even in the absence of historical barriers to gene flow. In the second part of the analysis, we assessed the probability that pairs of gene genealogies from a set of population samples agree upon the location of a geographical barrier. Our findings extend earlier reports by demonstrating that spatially concordant phylogeographic breaks across independent neutral loci normally emerge only in the presence of longstanding historical barriers to gene flow. Genealogical concordance across multiple loci thus becomes a deciding criterion by which to distinguish between stochastic and deterministic causation in accounting for phylogeographic discontinuities in continuously distributed species.

The Extinction Dynamics of Bacterial Pseudogenes

Pseudogenes are usually considered to be completely neutral sequences whose evolution is shaped by random mutations and chance events. It is possible, however, for disrupted genes to generate products that are deleterious due either to the energetic costs of their transcription and translation or to the formation of toxic proteins. We found that after their initial formation, the youngest pseudogenes in Salmonella genomes have a very high likelihood of being removed by deletional processes and are eliminated too rapidly to be governed by a strictly neutral model of stochastic loss. Those few highly degraded pseudogenes that have persisted in Salmonella genomes correspond to genes with low expression levels and low connectivity in gene networks, such that their inactivation and any initial deleterious effects associated with their inactivation are buffered. Although pseudogenes have long been considered the paradigm of neutral evolution, the distribution of pseudogenes among Salmonella strains indicates that removal of many of these apparently functionless regions is attributable to positive selection.

Inferring clocks when lacking rocks: the variable rates of molecular evolution in bacteria

BackgroundBecause bacteria do not have a robust fossil record, attempts to infer the timing of events in their evolutionary history requires comparisons of molecular sequences. This use of molecular clocks is based on the assumptions that substitution rates for homologous genes or sites are fairly constant through time and across taxa. Violation of these conditions can lead to erroneous inferences and result in estimates that are off by orders of magnitude. In this study, we examine the consistency of substitution rates among a set of conserved genes in diverse bacterial lineages, and address the questions regarding the validity of molecular dating.ResultsBy examining the evolution of 16S rRNA gene in obligate endosymbionts, which can be calibrated by the fossil record of their hosts, we found that the rates are consistent within a clade but varied widely across different bacterial lineages. Genome-wide estimates of nonsynonymous and synonymous substitutions suggest that these two measures are highly variable in their rates across bacterial taxa. Genetic drift plays a fundamental role in determining the accumulation of substitutions in 16S rRNA genes and at nonsynonymous sites. Moreover, divergence estimates based on a set of universally conserved protein-coding genes also exhibit low correspondence to those based on 16S rRNA genes.ConclusionOur results document a wide range of substitution rates across genes and bacterial taxa. This high level of variation cautions against the assumption of a universal molecular clock for inferring divergence times in bacteria. However, by applying relative-rate tests to homologous genes, it is possible to derive reliable local clocks that can be used to calibrate bacterial evolution.ReviewersThis article was reviewed by Adam Eyre-Walker, Simonetta Gribaldo and Tal Pupko (nominated by Dan Graur).

The fate of new bacterial genes

Bacteria experience a continual influx of novel genetic material from a wide range of sources and yet their genomes remain relatively small. This aspect of bacterial evolution indicates that most newly arriving sequences are rapidly eliminated; however, numerous new genes persist, as evident from the presence of unique genes in almost all bacterial genomes. This review summarizes the methods for identifying new genes in bacterial genomes and examines the features that promote the retention and elimination of these evolutionary novelties.

Comparative genome analysis of Spiroplasma melliferum IPMB4A, a honeybee-associated bacterium

BackgroundThe genus Spiroplasma contains a group of helical, motile, and wall-less bacteria in the class Mollicutes. Similar to other members of this class, such as the animal-pathogenic Mycoplasma and the plant-pathogenic ‘Candidatus Phytoplasma’, all characterized Spiroplasma species were found to be associated with eukaryotic hosts. While most of the Spiroplasma species appeared to be harmless commensals of insects, a small number of species have evolved pathogenicity toward various arthropods and plants. In this study, we isolated a novel strain of honeybee-associated S. melliferum and investigated its genetic composition and evolutionary history by whole-genome shotgun sequencing and comparative analysis with other Mollicutes genomes.ResultsThe whole-genome shotgun sequencing of S. melliferum IPMB4A produced a draft assembly that was ~1.1 Mb in size and covered ~80% of the chromosome. Similar to other Spiroplasma genomes that have been studied to date, we found that this genome contains abundant repetitive sequences that originated from plectrovirus insertions. These phage fragments represented a major obstacle in obtaining a complete genome sequence of Spiroplasma with the current sequencing technology. Comparative analysis of S. melliferum IPMB4A with other Spiroplasma genomes revealed that these phages may have facilitated extensive genome rearrangements in these bacteria and contributed to horizontal gene transfers that led to species-specific adaptation to different eukaryotic hosts. In addition, comparison of gene content with other Mollicutes suggested that the common ancestor of the SEM (Spiroplasma, Entomoplasma, and Mycoplasma) clade may have had a relatively large genome and flexible metabolic capacity; the extremely reduced genomes of present day Mycoplasma and ‘Candidatus Phytoplasma’ species are likely to be the result of independent gene losses in these lineages.ConclusionsThe findings in this study highlighted the significance of phage insertions and horizontal gene transfer in the evolution of bacterial genomes and acquisition of pathogenicity. Furthermore, the inclusion of Spiroplasma in comparative analysis has improved our understanding of genome evolution in Mollicutes. Future improvements in the taxon sampling of available genome sequences in this group are required to provide further insights into the evolution of these important pathogens of humans, animals, and plants.

Complete Genomes of Two Dipteran-Associated Spiroplasmas Provided Insights into the Origin, Dynamics, and Impacts of Viral Invasion in Spiroplasma

Spiroplasma is a genus of wall-less, low-GC, Gram-positive bacteria with helical morphology. As commensals or pathogens of plants, insects, ticks, or crustaceans, they are closely related with mycoplasmas and form a monophyletic group (Spiroplasma–Entomoplasmataceae–Mycoides) with Mycoplasma mycoides and its relatives. In this study, we report the complete genome sequences of Spiroplasma chrysopicola and S. syrphidicola from the Chrysopicola clade. These species form the sister group to the Citri clade, which includes several well-known pathogenic spiroplasmas. Surprisingly, these two newly available genomes from the Chrysopicola clade contain no plectroviral genes, which were found to be highly repetitive in the previously sequenced genomes from the Citri clade. Based on the genome alignment and patterns of GC-skew, these two Chrysopicola genomes appear to be relatively stable, rather than being highly rearranged as those from the Citri clade. Phylogenetic analyses suggest that the susceptibility to plectroviral invasion probably originated in the common ancestor of the Citri clade or one of its subclades. This susceptibility may be attributed to the absence of antiviral systems found in the Chrysopicola clade. Using the virus-free genomes of the Chrysopicola clade as references, we inferred the putative viral integration sites in the Citri genomes. Comparisons of syntenic regions suggest that the extensive viral invasion in the Citri clade promoted genome rearrangements and expansions. More importantly, the viral invasion may have facilitated horizontal gene transfers that contributed to adaptation in the Citri clade.