Etsuko N. Moriyama

Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors

In vertebrates, including humans, individuals harbor gut microbial communities whose species composition and relative proportions of dominant microbial groups are tremendously varied. Although external and stochastic factors clearly contribute to the individuality of the microbiota, the fundamental principles dictating how environmental factors and host genetic factors combine to shape this complex ecosystem are largely unknown and require systematic study. Here we examined factors that affect microbiota composition in a large (n = 645) mouse advanced intercross line originating from a cross between C57BL/6J and an ICR-derived outbred line (HR). Quantitative pyrosequencing of the microbiota defined a core measurable microbiota (CMM) of 64 conserved taxonomic groups that varied quantitatively across most animals in the population. Although some of this variation can be explained by litter and cohort effects, individual host genotype had a measurable contribution. Testing of the CMM abundances for cosegregation with 530 fully informative SNP markers identified 18 host quantitative trait loci (QTL) that show significant or suggestive genome-wide linkage with relative abundances of specific microbial taxa. These QTL affect microbiota composition in three ways; some loci control individual microbial species, some control groups of related taxa, and some have putative pleiotropic effects on groups of distantly related organisms. These data provide clear evidence for the importance of host genetic control in shaping individual microbiome diversity in mammals, a key step toward understanding the factors that govern the assemblages of gut microbiota associated with complex diseases.

Diversification of the gut symbiont Lactobacillus reuteri as a result of host-driven evolution

The vertebrate digestive tract, including that of humans, is the habitat to trillions of bacteria that are of significant importance to host biology and health. Although these communities are often postulated to have coevolved with their hosts, evidence is lacking, yet critical for our understanding of microbial symbiosis in vertebrates. To gain insight into the evolution of a gut symbiont, we have characterized the population genetic structure and phylogeny of Lactobacillus reuteri strains isolated from six different host species (human, mouse, rat, pig, chicken and turkey) using Amplified-Fragment Length Polymorphism (AFLP) and Multi-Locus Sequence Analysis (MLSA). The results revealed considerable genetic heterogeneity within the L. reuteri population and distinct monophyletic clades reflecting host origin but not provenance. The evolutionary patterns detected indicate a long-term association of L. reuteri lineages with particular vertebrate species and host-driven diversification. Results from a competition experiment in a gnotobiotic mouse model revealed that rodent isolates showed elevated ecological performance, indicating that evolution of L. reuteri lineages was adaptive. These findings provide evidence that some vertebrate gut microbes are not promiscuous, but have diversified into host-adapted lineages by a long-term evolutionary process, allowing the development of a highly specialized symbiosis.

Codon usage in twelve species of Drosophila

BackgroundCodon usage bias (CUB), the uneven use of synonymous codons, is a ubiquitous observation in virtually all organisms examined. The pattern of codon usage is generally similar among closely related species, but differs significantly among distantly related organisms, e.g., bacteria, yeast, and Drosophila. Several explanations for CUB have been offered and some have been supported by observations and experiments, although a thorough understanding of the evolutionary forces (random drift, mutation bias, and selection) and their relative importance remains to be determined. The recently available complete genome DNA sequences of twelve phylogenetically defined species of Drosophila offer a hitherto unprecedented opportunity to examine these problems. We report here the patterns of codon usage in the twelve species and offer insights on possible evolutionary forces involved.Results(1) Codon usage is quite stable across 11/12 of the species: G- and especially C-ending codons are used most frequently, thus defining the preferred codons. (2) The only amino acid that changes in preferred codon is Serine with six species of the melanogaster group favoring TCC while the other species, particularly subgenus Drosophila species, favor AGC. (3) D. willistoni is an exception to these generalizations in having a shifted codon usage for seven amino acids toward A/T in the wobble position. (4) Amino acids differ in their contribution to overall CUB, Leu having the greatest and Asp the least. (5) Among two-fold degenerate amino acids, A/G ending amino acids have more selection on codon usage than T/C ending amino acids. (6) Among the different chromosome arms or elements, genes on the non-recombining element F (dot chromosome) have the least CUB, while genes on the element A (X chromosome) have the most. (7) Introns indicate that mutation bias in all species is approximately 2:1, AT:GC, the opposite of codon usage bias. (8) There is also evidence for some overall regional bias in base composition that may influence codon usage.ConclusionOverall, these results suggest that natural selection has acted on codon usage in the genus Drosophila, at least often enough to leave a footprint of selection in modern genomes. However, there is evidence in the data that random forces (drift and mutation) have also left patterns in the data, especially in genes under weak selection for codon usage for example genes in regions of low recombination. The documentation of codon usage patterns in each of these twelve genomes also aids in ongoing annotation efforts.

Arabidopsis G-protein interactome reveals connections to cell wall carbohydrates and morphogenesis

The heterotrimeric G‐protein complex is minimally composed of Gα, Gβ, and Gγ subunits. In the classic scenario, the G‐protein complex is the nexus in signaling from the plasma membrane, where the heterotrimeric G‐protein associates with heptahelical G‐protein‐coupled receptors (GPCRs), to cytoplasmic target proteins called effectors. Although a number of effectors are known in metazoans and fungi, none of these are predicted to exist in their canonical forms in plants. To identify ab initio plant G‐protein effectors and scaffold proteins, we screened a set of proteins from the G‐protein complex using two‐hybrid complementation in yeast. After deep and exhaustive interrogation, we detected 544 interactions between 434 proteins, of which 68 highly interconnected proteins form the core G‐protein interactome. Within this core, over half of the interactions comprising two‐thirds of the nodes were retested and validated as genuine in planta. Co‐expression analysis in combination with phenotyping of loss‐of‐function mutations in a set of core interactome genes revealed a novel role for G‐proteins in regulating cell wall modification.

Mining the Arabidopsis thaliana genome for highly-divergent seven transmembrane receptors

To identify divergent seven-transmembrane receptor (7TMR) candidates from the Arabidopsis thaliana genome, multiple protein classification methods were combined, including both alignment-based and alignment-free classifiers. This resolved problems in optimally training individual classifiers using limited and divergent samples, and increased stringency for candidate proteins. We identified 394 proteins as 7TMR candidates and highlighted 54 with corresponding expression patterns for further investigation.

AN EVOLUTIONARY EXPRESSED SEQUENCE TAG ANALYSIS OF DROSOPHILA SPERMATHECA GENES

Abstract This study investigates genes enriched for expression in the spermatheca, the long-term sperm storage organ (SSO) of female Drosophila. SSO genes are likely to play an important role in processes of sexual selection such as sperm competition and cryptic female choice. Although there is keen interest in the mechanisms of sexual selection at the molecular level, very little is known about the female genes that are involved. In the present study, a high proportion of genes enriched for expression in the spermatheca are evolving rapidly. Most of the rapidly evolving genes are proteases and genes of unknown function that could play a specialized role in the spermatheca. A high percentage of the rapidly evolving genes have secretion signals and thus could encode proteins that directly interact with ejaculate proteins and coevolve with them. In addition to identifying rapidly evolving genes, the present study documents categories of genes that could play a role in spermatheca function such as storing, maintaining, and utilizing sperm. In general, candidate genes discovered in this study could play a key role in sperm competition, cryptic female choice of sperm, and sexually antagonistic coevolution, and ultimately speciation.

Analysis of a Shift in Codon Usage in Drosophila

In order to gain further insight into a shift in codon usage first observed in Drosophila willistoni we have analyzed seven genes in six species in the lineage leading to D. willistoni. This lineage contains the willistoni and saltans species groups. Sequences were obtained from GenBank or newly sequenced for this study. All species studied showed significant difference in codon usage compared to D. melanogaster for about one third of all amino acids. Within the willistoni/saltans lineage, codon usage is homogeneous, indicating that the shift in codon usage occurred prior to the diversification of extant species in this lineage which we estimate to date to about 20 million years ago. Thus the shift is old and has been stable. We also examined introns from these genes and the G/C composition at four-fold degenerate sites in an effort to detect a change in mutation bias. There is little or no evidence for a difference in mutation bias compared to D. melanogaster. We also considered whether relaxed selection (possibly due to reduced population sizes) or reduced recombination (due to numerous naturally occurring inversions) could account for the shift and concluded these factors alone are insufficient to explain the patterns observed. A change in the relative abundance of isoaccepting tRNAs is one of the few explanations that can account for the observations. Particularly intriguing is the fact that the greatest changes in codon usage have occurred for amino acids with two-fold C/T ending codons for which it is known that posttranscriptional modification occurs in tRNAs from a G in the wobble position to Queuosine that changes optimal binding from C to a slight preference for U. However, we do not argue that this shift was adaptive in nature, rather it may be an example of a “frozen accident.”

Biological Sequence Simulation for Testing Complex Evolutionary Hypotheses: indel-Seq-Gen Version 2.0

Sequence simulation is an important tool in validating biological hypotheses as well as testing various bioinformatics and molecular evolutionary methods. Hypothesis testing relies on the representational ability of the sequence simulation method. Simple hypotheses are testable through simulation of random, homogeneously evolving sequence sets. However, testing complex hypotheses, for example, local similarities, requires simulation of sequence evolution under heterogeneous models. To this end, we previously introduced indel-Seq-Gen version 1.0 (iSGv1.0; indel, insertion/deletion). iSGv1.0 allowed heterogeneous protein evolution and motif conservation as well as insertion and deletion constraints in subsequences. Despite these advances, for complex hypothesis testing, neither iSGv1.0 nor other currently available sequence simulation methods is sufficient. indel-Seq-Gen version 2.0 (iSGv2.0) aims at simulating evolution of highly divergent DNA sequences and protein superfamilies. iSGv2.0 improves upon iSGv1.0 through the addition of lineage-specific evolution, motif conservation using PROSITE-like regular expressions, indel tracking, subsequence-length constraints, as well as coding and noncoding DNA evolution. Furthermore, we formalize the sequence representation used for iSGv2.0 and uncover a flaw in the modeling of indels used in current state of the art methods, which biases simulation results for hypotheses involving indels. We fix this flaw in iSGv2.0 by using a novel discrete stepping procedure. Finally, we present an example simulation of the calycin-superfamily sequences and compare the performance of iSGv2.0 with iSGv1.0 and random model of sequence evolution.

Molecular evolution of urea amidolyase and urea carboxylase in fungi

BackgroundUrea amidolyase breaks down urea into ammonia and carbon dioxide in a two-step process, while another enzyme, urease, does this in a one step-process. Urea amidolyase has been found only in some fungal species among eukaryotes. It contains two major domains: the amidase and urea carboxylase domains. A shorter form of urea amidolyase is known as urea carboxylase and has no amidase domain. Eukaryotic urea carboxylase has been found only in several fungal species and green algae. In order to elucidate the evolutionary origin of urea amidolyase and urea carboxylase, we studied the distribution of urea amidolyase, urea carboxylase, as well as other proteins including urease, across kingdoms.ResultsAmong the 64 fungal species we examined, only those in two Ascomycota classes (Sordariomycetes and Saccharomycetes) had the urea amidolyase sequences. Urea carboxylase was found in many but not all of the species in the phylum Basidiomycota and in the subphylum Pezizomycotina (phylum Ascomycota). It was completely absent from the class Saccharomycetes (phylum Ascomycota; subphylum Saccharomycotina). Four Sordariomycetes species we examined had both the urea carboxylase and the urea amidolyase sequences. Phylogenetic analysis showed that these two enzymes appeared to have gone through independent evolution since their bacterial origin. The amidase domain and the urea carboxylase domain sequences from fungal urea amidolyases clustered strongly together with the amidase and urea carboxylase sequences, respectively, from a small number of beta- and gammaproteobacteria. On the other hand, fungal urea carboxylase proteins clustered together with another copy of urea carboxylases distributed broadly among bacteria. The urease proteins were found in all the fungal species examined except for those of the subphylum Saccharomycotina.ConclusionsWe conclude that the urea amidolyase genes currently found only in fungi are the results of a horizontal gene transfer event from beta-, gamma-, or related species of proteobacteria. The event took place before the divergence of the subphyla Pezizomycotina and Saccharomycotina but after the divergence of the subphylum Taphrinomycotina. Urea carboxylase genes currently found in fungi and other limited organisms were also likely derived from another ancestral gene in bacteria. Our study presented another important example showing plastic and opportunistic genome evolution in bacteria and fungi and their evolutionary interplay.

MITOCHONDRIAL DNA PHYLOGENIES FOR THE DROSOPHILA OBSCURA GROUP

Species belonging to the obscura group of the genus Drosophila have long held a central position in evolutionary studies, especially in experimental population genetics. Despite the considerable amount of accumulated knowledge, many of the phylogenetic relationships of the species in the group remain unclear. Here we present DNA sequence data for the mitochondrial gene cytochrome oxidase I (COI) for 13 species native to both the Old and New Worlds. We combine these data with seven other mitochondrial gene sequences from previous studies, for a total of over 3 kb per species. Strongly supported conclusions include: (1) the two North American subgroups, pseudoobscura and affinis, are each monophyletic; and (2) among Eurasian species two unambiguous clades are identified, one containing D. tristis, D. ambigua, and D. obscura and the other containing D. guanche, D. subobscura, and D. madeirensis. Constructing firm hypotheses connecting these four major clades is problematic with all datasets. Major ambiguities are the number of invasions giving rise to the North American obscura species and the relationships among the Eurasian species. The inadequacy of the mtDNA data to resolve these ambiguities does not reside in lack of changes; the transversions‐only parsimony tree has 283 informative characters. Rather, the problems are likely intrinsic to the history of the group: while radiating in temperate Eurasia, North America was colonized once or twice, followed by one or two radiations in the New World.