Molecular Biology and Evolution - current issue

jModelTest: Phylogenetic Model Averaging

Posada, D. — 2008-06-13

jModelTest is a new program for the statistical selection of models of nucleotide substitution based on "Phyml" (Guindon and Gascuel 2003. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 52:696–704.). It implements 5 different selection strategies, including "hierarchical and dynamical likelihood ratio tests," the "Akaike information criterion," the "Bayesian information criterion," and a "decision-theoretic performance-based" approach. This program also calculates the relative importance and model-averaged estimates of substitution parameters, including a model-averaged estimate of the phylogeny. jModelTest is written in Java and runs under Mac OSX, Windows, and Unix systems with a Java Runtime Environment installed. The program, including documentation, can be freely downloaded from the software section at http://darwin.uvigo.es.

The Timing of Selection at the Human FOXP2 Gene

Coop, G., Bullaughey, K., Luca, F., Przeworski, M. — 2008-06-13

Krause J, Lalueza-Fox C, Orlando L, et al. recently examined patterns of genetic variation at FOXP2 in 2 Neanderthals. This gene is of particular interest because it is involved in speech and language and was previously shown to harbor the signature of recent positive selection. The authors found the same 2 amino acid substitutions in Neanderthals as in modern humans. Assuming that these sites were the targets of selection and no interbreeding between the 2 groups, they concluded that selection at FOXP2 occurred before the populations split, over 300 thousand years ago. Here, we show that the data are unlikely under this scenario but may instead be consistent with low rates of gene flow between modern humans and Neanderthals. We also collect additional data and introduce a modeling framework to estimate levels of modern human contamination of the Neanderthal samples. We find that, depending on the assumptions, additional control experiments may be needed to rule out contamination at FOXP2.

Positive Selection and Expression Divergence Following Gene Duplication in the Sunflower CYCLOIDEA Gene Family

Chapman, M. A., Leebens-Mack, J. H., Burke, J. M. — 2008-06-13

Members of the CYCLOIDEA (CYC)/TEOSINTE-BRANCHED1 (TB1) group of transcription factors have been implicated in the evolution of zygomorphic (i.e., bilaterally symmetric) flowers in Antirrhinum and Lotus and the loss of branching phenotype during the domestication of maize. The composite inflorescences of sunflower (Helianthus annuus L. Asteraceae) contain both zygomorphic and actinomorphic (i.e., radially symmetric) florets (rays and disks, respectively), and the cultivated sunflower has evolved an unbranched phenotype in response to domestication from its highly branched wild progenitor; hence, genes related to CYC/TB1 are of great interest in this study system. We identified 10 members of the CYC/TB1 gene family in sunflower, which is more than found in any other species investigated to date. Phylogenetic analysis indicates that these genes occur in 3 distinct clades, consistent with previous research in other eudicot species. A combination of dating the duplication events and linkage mapping indicates that only some of the duplications were associated with polyploidization. Cosegregation between CYC-like genes and branching-related quantitative trait loci suggest a minor, if any, role for these genes in conferring differences in branching. However, the expression patterns of one gene suggest a possible role in the development of ray versus disk florets. Molecular evolutionary analyses reveal that residues in the conserved domains were the targets of positive selection following gene duplication. Taken together, these results indicate that gene duplication and functional divergence have played a major role in diversification of the sunflower CYC gene family.

Novel Transcriptome Patterns Accompany Evolutionary Restoration of Defective Social Development in the Bacterium Myxococcus xanthus

Kadam, S. V., Wegener-Feldbrugge, S., Sogaard-Andersen, L., Velicer, G. J. — 2008-06-13

Evolutionary trait losses can be restored by direct reversion or by compensatory pathways. Upon starvation, the bacterium Myxococcus xanthus develops into spore-bearing fruiting bodies, but this ability can be rapidly lost during evolution. Some developmentally defective strains of M. xanthus "cheat" on proficient strains during development by superior sporulation in mixed cultures. Here, we examine transcriptomic patterns accompanying the evolution of a cheater (obligate cheater [OC]) to a developmentally competent strain (PX) by a single mutation. Using quantitative real-time–polymerase chain reaction analysis of 5 genes essential for development, we initially show that restoration of development in strain PX was associated with increased expression of 4 of these genes, not only relative to OC but also relative to the developmentally proficient ancestor of both OC and PX (wild type [WT]). Global transcriptome analyses showed further that developmental expression of well more than 100 genes differ significantly between PX and the proficient WT ancestor. Moreover, the expression profile of PX was found to differ from that of WT more than does that of the defective intermediate strain OC. These results show that the restoration of a complex trait is accompanied by novel expression patterns across a large number and wide variety of genes, rather than by a large-scale return to ancestral expression patterns.

Parallel Rate Heterogeneity in Chloroplast and Mitochondrial Genomes of Brazil Nut Trees (Lecythidaceae) Is Consistent with Lineage Effects

Soria-Hernanz, D. F., Braverman, J. M., Hamilton, M. B. — 2008-06-13

We investigated whether relative rates of divergence were correlated between the mitochondrial and chloroplast genomes as expected under lineage effects or were genome specific as expected with locus-specific effects. Five mitochondrial noncoding regions (nad1B_C, nad4exon1_2, nad7exon2_3, nad7exon3_4, and rps14-cob) for 21 samples from Lecythidaceae were sequenced. Three chloroplast regions (rpl20-5'rps12, trnS-trnG, and psbA-trnH) were sequenced to expand the taxa in an existing data set. Absolute rates of nucleotide and insertion and deletion (indel) changes were 13 times faster in the chloroplast genome than in the mitochondrial genome. Similar indel length frequency distributions for both organelles suggested that common mechanisms were responsible for generating indels. Molecular clock tests applied to phylogenetic trees estimated from mitochondrial and chloroplast sequences revealed global rate heterogeneity of nucleotide substitution. Maximum likelihood and Tajima's 1D relative rate tests show that Lecythis zabucajo exhibited a rate acceleration for both the mitochondrial and chloroplast sequences. Whereas Eschweilera romeu-cardosoi showed a significant rate slowdown for chloroplast sequences, the mitochondrial sequences for 3 Eschweilera taxa showed evidence for a rate slowdown only when compared with L. zabucajo. Significant rate heterogeneity was also observed for indel changes in the mitochondrial genome but not for the chloroplast. The lack of mitochondrial nucleotide changes for some taxa as well as chloroplast indel homoplasy may have limited the power of relative rate tests to detect rate variation. Relative ratio tests consistently indicated rate proportionality among branch lengths between the mitochondrial and chloroplast phylogenetic trees. The relative ratio tests showed that taxa possessing rate heterogeneity had parallel relative divergence rates in both mitochondrial and chloroplast sequences as expected under lineage effects. A neutral replication–dependent model of rate heterogeneity for both nucleotide and indel changes provides a simple explanation for common patterns of rate heterogeneity across the 2 organelle genomes in Lecythidaceae. The lineage effects observed here were uncoupled from annual/perennial habit because all the species from this study are perennial.

Plastid-Derived Genes in the Nonphotosynthetic Alveolate Oxyrrhis marina

Slamovits, C. H., Keeling, P. J. — 2008-06-13

Reconstructing the history of plastid acquisition and loss in the alveolate protists is a difficult problem because our knowledge of the distribution of plastids in extant lineages is incomplete due to the possible presence of cryptic, nonphotosynthetic plastids in several lineages. The discovery of the apicoplast in apicomplexan parasites has drawn attention to this problem and, more specifically, to the question of whether many other nonphotosynthetic lineages also contain cryptic plastids or are derived from plastid-containing ancestors. Oxyrrhis marina is one such organism: It is a heterotrophic, early-branching member of the dinoflagellate lineage for which there is no evidence of a plastid. To investigate the possibility that O. marina is derived from a photosynthetic ancestor, we have generated and analyzed a large-scale EST database and searched for evidence of plastid-derived genes. Here, we describe 8 genes whose phylogeny shows them to be derived from plastid-targeted homologues. These genes encode proteins from several pathways known to be localized in the plastids of other algae, including synthesis of tetrapyrroles, isoprenoids, and amino acids, as well as carbon metabolism and oxygen detoxification. The 5' end of 5 cDNAs were also characterized using cap-dependent or spliced leader–mediated reverse transcriptase–polymerase chain reaction, revealing that at least 4 of these genes have retained leaders that are similar in nature to the plastid-targeting signals of other secondary plastids, suggesting that these proteins may be targeted to a cryptic organelle. At least 2 genes do not encode such leaders, and their products may presently function in the cytosol. Altogether, the presence of plastid-derived genes in O. marina shows that its ancestors contained a plastid, and the pathways represented by the genes and presence of targeting signals on at least some of the genes further suggests that a relict organelle may still exist to fulfill plastid metabolic functions.

An Improved General Amino Acid Replacement Matrix

Le, S. Q., Gascuel, O. — 2008-06-13

Amino acid replacement matrices are an essential basis of protein phylogenetics. They are used to compute substitution probabilities along phylogeny branches and thus the likelihood of the data. They are also essential in protein alignment. A number of replacement matrices and methods to estimate these matrices from protein alignments have been proposed since the seminal work of Dayhoff et al. (1972). An important advance was achieved by Whelan and Goldman (2001) and their WAG matrix, thanks to an efficient maximum likelihood estimation approach that accounts for the phylogenies of sequences within each training alignment. We further refine this method by incorporating the variability of evolutionary rates across sites in the matrix estimation and using a much larger and diverse database than BRKALN, which was used to estimate WAG. To estimate our new matrix (called LG after the authors), we use an adaptation of the XRATE software and 3,912 alignments from Pfam, comprising ~50,000 sequences and ~6.5 million residues overall. To evaluate the LG performance, we use an independent sample consisting of 59 alignments from TreeBase and randomly divide Pfam alignments into 3,412 training and 500 test alignments. The comparison with WAG and JTT shows a clear likelihood improvement. With TreeBase, we find that 1) the average Akaike information criterion gain per site is 0.25 and 0.42, when compared with WAG and JTT, respectively; 2) LG is significantly better than WAG for 38 alignments (among 59), and significantly worse with 2 alignments only; and 3) tree topologies inferred with LG, WAG, and JTT frequently differ, indicating that using LG impacts not only the likelihood value but also the output tree. Results with the test alignments from Pfam are analogous. LG and a PHYML implementation can be downloaded from http://atgc.lirmm.fr/LG.

The Capsid of the T4 Phage Superfamily: The Evolution, Diversity, and Structure of Some of the Most Prevalent Proteins in the Biosphere

Comeau, A. M., Krisch, H. M. — 2008-06-13

The Escherichia coli bacteriophage T4 has served as a classic system in phage biology for more than 60 years. Only recently have phylogenetic analyses and genomic comparisons demonstrated the existence of a large, diverse, and widespread superfamily of T4-like phages in the environment. We report here on the T4-like major capsid protein (MCP) sequences that were obtained by targeted polymerase chain reaction (PCR) of marine environmental samples. This analysis was then expanded to include 1,000s of new sequences of T4-like capsid genes from the metagenomic data obtained during the Sorcerer II Global Ocean Sampling (GOS) expedition. This data compilation reveals that the diversity of the major and minor capsid proteins from the GOS metagenome follows the same general patterns as the sequences from cultured phage genomes. Interestingly, the new MCP sequences obtained by PCR targeted to MCP sequences in environmental samples are more divergent (deeper branching) than the vast majority of the MCP sequences coming from the other sources. The marine T4-like phage population appears to be largely dominated by the T4-like cyanophages. Using ~1,400 T4-like MCP sequences from various sources, we mapped the degree of sequence conservation on a structural model of the T4-like MCP. The results indicate that within the T4 superfamily there are some clear phylogenetic groups with regard to the more conserved and more variable domains of the MCP. Such differences can be correlated with variations in capsid morphology, the arrangement of the MCP lattice, and the presence of different capsid accessory proteins between the subgroups of the T4 superfamily.

Differential Evolution of the 13 Atlantic Salmon Hox Clusters

Mungpakdee, S., Seo, H.-C., Angotzi, A. R., Dong, X., Akalin, A., Chourrout, D. — 2008-06-13

Hox cluster organization represents a valuable marker to study the effects of recent genome duplication in salmonid fish (25–100 Mya). Using polymerase chain reaction amplification of cDNAs, BAC library screening, and genome walking, we reconstructed 13 Hox clusters in the Atlantic salmon containing 118 Hox genes including 8 pseudogenes. Hox paralogs resulting from the genome duplication preceding the radiation of ray-finned fish have been much better preserved in salmon than in other model teleosts. The last genome duplication in the salmon lineage has been followed by the loss of 1 of the 4 HoxA clusters. Four rounds of genome duplication after the vertebrate ancestor salmon Hox clusters display the main organizational features of vertebrate Hox clusters, with Hox genes exclusively that are densely packed in the same orientation. Recently, duplicated Hox clusters have engaged a process of divergence, with several cases of pseudogenization or asymmetrical evolution of Hox gene duplicates, and a marked erosion of identity in noncoding sequences. Strikingly, the level of divergence attained strongly depends on the Hox cluster pairs rather than on the Hox genes within each cluster. It is particularly high between both HoxBb clusters and both HoxDa clusters, whereas both HoxBa clusters remained virtually identical. Positive selection on the Hox protein–coding sequences could not be detected.

Unequal Rates of Y Chromosome Gene Divergence during Speciation of the Family Ursidae

Nakagome, S., Pecon-Slattery, J., Masuda, R. — 2008-06-13

Evolution of the bear family Ursidae is well investigated in terms of morphological, paleontological, and genetic features. However, several phylogenetic ambiguities occur within the subfamily Ursinae (the family Ursidae excluding the giant panda and spectacled bear), which may correlate with behavioral traits of female philopatry and male-biased dispersal which form the basis of the observed matriarchal population structure in these species. In the process of bear evolution, we investigate the premise that such behavioral traits may be reflected in patterns of variation among genes with different modes of inheritance: matrilineal mitochondrial DNA (mtDNA), patrilineal Y chromosome, biparentally inherited autosomes, and the X chromosome. In the present study, we sequenced 3 Y-linked genes (3,453 bp) and 4 X-linked genes (4,960 bp) and reanalyzed previously published sequences from autosome genes (2,347 bp) in ursid species to investigate differences in evolutionary rates associated with patterns of inheritance. The results describe topological incongruence between sex-linked genes and autosome genes and between nuclear DNA and mtDNA. In more ancestral branches within the bear phylogeny, Y-linked genes evolved faster than autosome and X-linked genes, consistent with expectations based on male-driven evolution. However, this pattern changes among branches leading to each species within the lineage of Ursinae whereby the evolutionary rates of Y-linked genes have fewer than expected substitutions. This inconsistency between more recent nodes of the bear phylogeny with more ancestral nodes may reflect the influences of sex-biased dispersal as well as molecular evolutionary characteristics of the Y chromosome, and stochastic events in species natural history, and phylogeography unique to ursine bears.

Unicellular Ca2+ Signaling 'Toolkit' at the Origin of Metazoa

Cai, X. — 2008-06-13

Ca²⁺ signaling pathways control many physiological processes in almost all types of animal cells such as fertilization, muscle contraction, hormone release, and learning and memory. Each animal cell type expresses a unique group of molecules from the Ca²⁺ signaling ‘toolkit’ to control spatiotemporal patterns of Ca²⁺ signaling. It is generally believed that the complex Ca²⁺ signaling ‘toolkit’ has arisen from the ancestral multicellular organisms to fit unique physiological roles of specialized cell types. Here, we demonstrate for the first time the presence of an extensive Ca²⁺ signaling ‘toolkit’ in the unicellular choanoflagellate Monosiga brevicollis. Choanoflagellates possess homologues of various types of animal plasma membrane Ca²⁺ channels including the store-operated channel, ligand-operated channels, voltage-operated channels, second messenger-operated channels, and 5 out of 6 animal transient receptor potential channel families. Choanoflagellates also contain homologues of inositol 1,4,5-trisphosphate receptors. Furthermore, choanoflagellates master a complete set of Ca²⁺ removal systems including plasma membrane and sarco/endoplasmic reticulum Ca²⁺ ATPases and homologues of 3 animal cation/Ca²⁺ exchanger families. Therefore, a complex Ca²⁺ signaling ‘toolkit’ might have evolved before the emergence of multicellular animals.

The Impact of the Austronesian Expansion: Evidence from mtDNA and Y Chromosome Diversity in the Admiralty Islands of Melanesia

2008-06-13

The genetic ancestry of Polynesians can be traced to both Asia and Melanesia, which presumably reflects admixture occurring between incoming Austronesians and resident non-Austronesians in Melanesia before the subsequent occupation of the greater Pacific; however, the genetic impact of the Austronesian expansion to Melanesia remains largely unknown. We therefore studied the diversity of nonrecombining Y chromosomal (NRY) and mitochondrial (mt) DNA in the Admiralty Islands, located north of mainland Papua New Guinea, and updated our previous data from Asia, Melanesia, and Polynesia with new NRY markers. The Admiralties are occupied today solely by Austronesian-speaking groups, but their human settlement history goes back 20,000 years prior to the arrival of Austronesians about 3,400 years ago. On the Admiralties, we found substantial mtDNA and NRY variation of both Austronesian and non-Austronesian origins, with higher frequencies of Asian mtDNA and Melanesian NRY haplogroups, similar to previous findings in Polynesia and perhaps as a consequence of Austronesian matrilocality. Thus, the Austronesian language replacement on the Admiralties (and elsewhere in Island Melanesia and coastal New Guinea) was accompanied by an incomplete genetic replacement that is more associated with mtDNA than with NRY diversity. These results provide further support for the "Slow Boat" model of Polynesian origins, according to which Polynesian ancestors originated from East Asia but genetically mixed with Melanesians before colonizing the Pacific. We also observed that non-Austronesian groups of coastal New Guinea and Island Melanesia had significantly higher frequencies of Asian mtDNA haplogroups than of Asian NRY haplogroups, suggesting sex-biased admixture perhaps as a consequence of non-Austronesian patrilocality. We additionally found that the predominant NRY haplogroup of Asian origin in the Admiralties (O-M110) likely originated in Taiwan, thus providing the first direct Y chromosome evidence for a Taiwanese origin of the Austronesian expansion. Furthermore, we identified a NRY haplogroup (K-P79, also found on the Admiralties) in Polynesians that most likely arose in the Bismarck Archipelago, providing the first direct link between northern Island Melanesia and Polynesia. These results significantly advance our understanding of the impact of the Austronesian expansion and human history in the Pacific region.

Selection on Amino Acid Substitutions in Arabidopsis

Foxe, J. P., Dar, V.-u.-N., Zheng, H., Nordborg, M., Gaut, B. S., Wright, S. I. — 2008-06-13

Studies of nucleotide diversity have found an excess of low-frequency amino acid polymorphisms segregating in Arabidopsis thaliana, suggesting a predominance of weak purifying selection acting on amino acid polymorphism in this inbreeding species. Here, we investigate levels of diversity and divergence at synonymous and nonsynonymous sites in 6 circumpolar populations of the outbreeding Arabidopsis lyrata and compare these results with A. thaliana, to test for differences in mutation and selection parameters across genes, populations, and species. We find that A. lyrata shows an excess of low-frequency nonsynonymous polymorphisms both within populations and species wide, consistent with weak purifying selection similar to the patterns observed in A. thaliana. Furthermore, nonsynonymous polymorphisms tend to be more restricted in their population distribution in A. lyrata, consistent with purifying selection preventing their geographic spread. Highly expressed genes show a reduced ratio of amino acid to synonymous change for both polymorphism and fixed differences, suggesting a general pattern of stronger purifying selection on high-expression proteins.

Coyotes Demonstrate How Habitat Specialization by Individuals of a Generalist Species Can Diversify Populations in a Heterogeneous Ecoregion

Sacks, B. N., Bannasch, D. L., Chomel, B. B., Ernest, H. B. — 2008-06-13

The tendency for individuals to disperse into habitat similar to their natal habitat has been observed in a wide range of species, although its population genetic consequences have received little study. Such behavior could lead to discrete habitat-specific population subdivisions even in the absence of physical dispersal barriers or habitat gaps. Previous studies of coyotes have supported this hypothesis in a small region of California, but its evolutionary significance ultimately depends on the extent and magnitude of habitat-specific subdivision. Here, we investigated these questions using autosomal, Y chromosome, and mitochondrial markers and >2,000 coyotes from a broad region, including 2 adjacent ecoregions with contrasting levels of habitat heterogeneity—the California Floristic Province (CFP) (heterogeneous landscape) and the Desert–Prairie ecoregion (DPE) (homogeneous landscape). Consistent with predictions, we found a close correspondence between population genetic structure and habitat subdivisions throughout the CFP and virtual panmixia over the larger DPE. Conversely, although genetic diversity was similar in these 2 ecoregions overall, it was lower within sites of the CFP, as would be the expected consequence of greater genetic drift within subregions. The magnitude of habitat-specific genetic subdivisions (i.e., genetic distance) in the CFP varied considerably, indicating complexity (e.g., asymmetric gene flow or extinction/recolonization), but, in general, was higher than that due to geographic distance or recent human-related barriers. Because habitat-specific structure can enhance a species' adaptive potential and resilience to changing environments, these findings suggest the CFP may constitute an evolutionarily important portion of the range for coyotes and sympatric species exhibiting habitat-specific population structure.

Systematic Survey for Novel Types of Prokaryotic Retroelements Based on Gene Neighborhood and Protein Architecture

Kojima, K. K., Kanehisa, M. — 2008-06-13

Retroelements, elements encoding reverse transcriptase (RT), are ubiquitous in eukaryotes and have a great influence on the evolution of our genome. Detailed information is available on eukaryotic retroelements; however, prokaryotic retroelements are poorly understood. Recently, new types of eukaryotic retroelements were characterized on the basis of their gene composition and their phylogenetic positions. Here we performed a systematic survey to identify novel types of prokaryotic retroelements by analyzing gene neighborhood and protein architecture. We found novel types of gene combination and examined whether they represent actual retroelements. Five monophyletic groups were identified that were distinct from characterized prokaryotic retroelements, showed specific gene combination, were distributed patchily, and included at least 1 example of recent integration. These results strongly indicated the frequent horizontal transfer of these elements. One group encoded DNA polymerase A. A possible function of DNA polymerase A in the life cycle of retroelements is catalyzing second-strand cDNA synthesis, which is DNA polymerization performed using a DNA template not an RNA template. Another group encoded both bacterial primase and carbon–nitrogen hydrolase. Primase is likely to synthesize primers to initiate reverse transcription. Two other groups also encoded carbon–nitrogen hydrolase as a fusion protein with RT. It is difficult to speculate on the function of hydrolase in the life cycle of retroelements. The last group encoded dual RT proteins, which are likely to form heterodimers during replication. The protein sets of these 5 groups of prokaryotic retroelements were completely different from those of eukaryotic retroelements, indicating that the survival constraints of prokaryotic elements were distinct from those of eukaryotic elements. It is likely that these prokaryotic retroelements are maintained as extrachromosomal DNA or RNA or are accidentally integrated into genomes. Our findings presented the possibility that many types of extrachromosomal prokaryotic retroelements remain to be characterized. In addition, we found 8 RT genes were associated with clustered regularly interspaced short palindrome repeats (CRISPRs) of the CRISPR–Cas system. These RT genes are likely to work in immunity against RNA phages via cDNA synthesis.

Organellar RNA Editing and Plant-Specific Extensions of Pentatricopeptide Repeat Proteins in Jungermanniid but not in Marchantiid Liverworts

Rudinger, M., Polsakiewicz, M., Knoop, V. — 2008-06-13

The pyrimidine exchange type of RNA editing in land plant (embryophyte) organelles has largely remained an enigma with respect to its biochemical mechanisms, the underlying specificities, and its raison d’être. Apparently arising with the earliest embryophytes, RNA editing is conspicuously absent in one clade of liverworts, the complex thalloid Marchantiidae. Several lines of evidence suggest that the large gene family of organelle-targeted RNA–binding pentatricopeptide repeat (PPR) proteins plays a fundamental role in the sequence-specific editing of organelle transcripts. We here describe the identification of PPR protein genes with plant-specific carboxyterminal (C-terminal) sequence signatures (E, E+, and DYW domains) in ferns, lycopodiophytes, mosses, hornworts, and jungermanniid liverworts, one subclass of the basal most clade of embryophytes, on DNA and cDNA level. In contrast, we were unable to identify these genes in a wide sampling of marchantiid liverworts (including the phylogenetic basal genus Blasia)—taxa for which no RNA editing is observed in the organelle transcripts. On the other hand, we found significant diversity of this type of PPR proteins also in Haplomitrium, a genus with an extremely high rate of RNA editing and a phylogenetic placement basal to all other liverworts. Although the presence of modularly extended PPR proteins correlates well with organelle RNA editing, the now apparent complete loss of an entire gene family from one clade of embryophytes, the marchantiid liverworts, remains puzzling.

A Phylogenetic Analysis of Indel Dynamics in the Cotton Genus

Grover, C. E., Yu, Y., Wing, R. A., Paterson, A. H., Wendel, J. F. — 2008-06-13

Genome size evolution is a dynamic process involving counterbalancing mechanisms whose actions vary across lineages and over time. Whereas the primary mechanism of expansion, transposable element (TE) amplification, has been widely documented, the evolutionary dynamics of genome contraction have been less thoroughly explored. To evaluate the relative impact and evolutionary stability of the mechanisms that affect genome size, we conducted a phylogenetic analysis of indel rates for 2 genomic regions in 4 Gossypium genomes: the 2 coresident genomes (A_T and D_T) of tetraploid cotton and its model diploid progenitors, Gossypium arboreum (A) and Gossypium raimondii (D). We determined the rates of sequence gain or loss along each branch, partitioned by mechanism, and how these changed during species divergence. In general, there has been a propensity toward growth of the diploid genomes and contraction in the polyploid. Most of the size difference between the diploid species occurred prior to polyploid divergence and was largely attributable to TE amplification in the A/A_T genome. After separating from the true parents of the polyploid genomes, both diploid genomes experienced slower sequence gain than in the ancestor, due to fewer TE insertions in the A genome and a combination of increased deletions and decreased TE insertions in the D genome. Both genomes of the polyploid displayed increased rates of deletion and decreased rates of insertion, leading to a rate of near stasis in D_T and overall contraction in A_T resulting in polyploid genome contraction. As expected, TE insertions contributed significantly to the genome size differences; however, intrastrand homologous recombination, although rare, had the most significant impact on the rate of deletion. Small indel data for the diploids suggest the possibility of a bias as the smaller genomes add less or delete more sequence through small indels than do the larger genomes, whereas data for the polyploid suggest increased sequence turnover in general (both as small deletions and small insertions). Illegitimate recombination, although not demonstrated to be a dominant mechanism of genome size change, was biased in the polyploid toward deletions, which may provide a partial explanation of polyploid genomic downsizing.

The Dscam Homologue of the Crustacean Daphnia Is Diversified by Alternative Splicing Like in Insects

2008-06-13

In insects, the homologue of the Down syndrome cell adhesion molecule (Dscam) is a unique case of a single-locus gene whose expression has extensive somatic diversification in both the nervous and immune systems. How this situation evolved is best understood through comparative studies. We describe structural, expression, and evolutionary aspects of a Dscam homolog in 2 species of the crustacean Daphnia. The Dscam of Daphnia generates up to 13,000 different transcripts by the alternative splicing of variable exons. This extends the taxonomic range of a highly diversified Dscam beyond the insects. Additionally, we have identified 4 alternative forms of the cytoplasmic tail that generate isoforms with or without inhibitory or activating immunoreceptor tyrosine–based motifs (ITIM and ITAM respectively), something not previously reported in insect's Dscam. In Daphnia, we detected exon usage variability in both the brain and hemocytes (the effector cells of immunity), suggesting that Dscam plays a role in the nervous and immune systems of crustaceans, as it does in insects. Phylogenetic analysis shows a high degree of amino acid conservation between Daphnia and insects except in the alternative exons, which diverge greatly between these taxa. Our analysis shows that the variable exons diverged before the split of the 2 Daphnia species and is in agreement with the nearest-neighbor model for the evolution of the alternative exons. The genealogy of the Dscam gene family from vertebrates and invertebrates confirmed that the highly diversified form of the gene evolved from a nondiversified form before the split of insects and crustaceans.

Mitochondrial Genome Evolution in the Social Amoebae

Heidel, A. J., Glockner, G. — 2008-06-13

Most mitochondria contain a core set of genes required for mitochondrial function, but beyond this base there are variable genomic features. The mitochondrial genome of the model species Dictyostelium discoideum demonstrated that the social amoebae mitochondrial genomes have a size between those of metazoans and plants, but no comparative study of social amoebae mitochondria has been performed. Here, we present a comparative analysis of social amoebae mitochondrial genomes using D. discoideum, Dictyostelium citrinum, Dictyostelium fasciculatum, and Polysphondylium pallidum. The social amoebae mitochondria have similar sizes, AT content, gene content and have a high level of synteny except for one segmental rearrangement and extensive displacement of tRNAs. From the species that contain the rearrangement, it can be concluded that the event occurred late in the evolution of social amoebas. A phylogeny using 36 mitochondrial genes produced a well-supported tree suggesting that the pairs of D. discoideum/D. citrinum and D. fasciculatum/P. pallidum are sister species although the position of the root is not certain. Group I introns and endonucleases are variable in number and location in the social amoebae. Phylogenies of the introns and endonucleases suggest that there have been multiple recent duplications or extinctions and confirm that endonucleases have the ability to insert into new areas. An analysis of dN/dS ratios in mitochondrial genes revealed that among groups of genes, adenosine triphosphate synthase complex genes have the highest ratio, whereas cytochrome oxidase and nicotinamide adenine dinucleotide (NADH) dehydrogenase genes had the lowest ratio. The genetic codes of D. citrinum, P. pallidum, and D. fasciculatum are the universal code although D. fasciculatum does not use the TGA stop codon. In D. fasciculatum, we demonstrate for the first time that a mitochondrial genome without the TGA stop codon still uses the release factor RF2 that recognizes TGA. Theories of how the genetic code can change and why RF2 may be a constraint against switching codes are discussed.

Recurrent Tandem Gene Duplication Gave Rise to Functionally Divergent Genes in Drosophila

Fan, C., Chen, Y., Long, M. — 2008-06-13

Tandem gene duplication is one of the major gene duplication mechanisms in eukaryotes, as illustrated by the prevalence of gene family clusters. Tandem duplicated paralogs usually share the same regulatory element, and as a consequence, they are likely to perform similar biological functions. Here, we provide an example of a newly evolved tandem duplicate acquiring novel functions, which were driven by positive selection. CG32708, CG32706, and CG6999 are 3 clustered genes residing in the X chromosome of Drosophila melanogaster. CG6999 and CG32708 have been examined for their molecular population genetic properties (Thornton and Long 2005). We further investigated the evolutionary forces acting on these genes with greater sample sizes and a broader approach that incorporate between-species divergence, using more variety of statistical methods. We explored the possible functional implications by characterizing the tissue-specific and developmental expression patterns of these genes. Sequence comparison of species within D. melanogaster subgroup reveals that this 3-gene cluster was created by 2 rounds of tandem gene duplication in the last 5 Myr. Based on phylogenetic analysis, CG32708 is clearly the parental copy that is shared by all species. CG32706 appears to have originated in the ancestor of Drosophila simulans and D. melanogaster about 5 Mya, and CG6999 is the newest duplicate that is unique to D. melanogaster. All 3 genes have different expression profiles, and CG6999 has in addition acquired a novel transcript. Biased polymorphism frequency spectrum, linkage disequilibrium, nucleotide substitution, and McDonald–Kreitman analyses suggested that the evolution of CG6999 and CG32706 were driven by positive Darwinian selection.

Smooth Skyride through a Rough Skyline: Bayesian Coalescent-Based Inference of Population Dynamics

Minin, V. N., Bloomquist, E. W., Suchard, M. A. — 2008-06-13

Kingman's coalescent process opens the door for estimation of population genetics model parameters from molecular sequences. One paramount parameter of interest is the effective population size. Temporal variation of this quantity characterizes the demographic history of a population. Because researchers are rarely able to choose a priori a deterministic model describing effective population size dynamics for data at hand, nonparametric curve-fitting methods based on multiple change-point (MCP) models have been developed. We propose an alternative to change-point modeling that exploits Gaussian Markov random fields to achieve temporal smoothing of the effective population size in a Bayesian framework. The main advantage of our approach is that, in contrast to MCP models, the explicit temporal smoothing does not require strong prior decisions. To approximate the posterior distribution of the population dynamics, we use efficient, fast mixing Markov chain Monte Carlo algorithms designed for highly structured Gaussian models. In a simulation study, we demonstrate that the proposed temporal smoothing method, named Bayesian skyride, successfully recovers "true" population size trajectories in all simulation scenarios and competes well with the MCP approaches without evoking strong prior assumptions. We apply our Bayesian skyride method to 2 real data sets. We analyze sequences of hepatitis C virus contemporaneously sampled in Egypt, reproducing all key known aspects of the viral population dynamics. Next, we estimate the demographic histories of human influenza A hemagglutinin sequences, serially sampled throughout 3 flu seasons.

Polyploid Speciation Did Not Confer Instant Reproductive Isolation in Capsella (Brassicaceae)

Slotte, T., Huang, H., Lascoux, M., Ceplitis, A. — 2008-06-13

Polyploid formation is a major mode of sympatric speciation in flowering plants. Unlike other speciation processes, polyploidization is often assumed to confer instant reproductive isolation. Shared polymorphism across ploidy levels has therefore often been attributed to multiple polyploid origins, whereas the alternative hypothesis of introgressive hybridization has rarely been rigorously tested. Here, we sequence 12 nuclear loci representing 6 genes duplicated by polyploidy in 92 accessions of the tetraploid Capsella bursa-pastoris together with the corresponding loci in 21 accessions of its close diploid relative Capsella rubella. In C. bursa-pastoris accessions from western Eurasia, where the 2 species occur in partial sympatry, we find higher levels of nucleotide diversity than in accessions from eastern Eurasia, where C. rubella does not grow. Furthermore, haplotypes are shared across ploidy levels at 4 loci in western but not in eastern Eurasia. We test whether haplotype sharing is due to retention of ancestral polymorphism or due to hybridization and introgression using a coalescent-based isolation-with-migration model. In western but not in eastern Eurasia, there is evidence for unidirectional gene flow from C. rubella to C. bursa-pastoris. An independent estimate of the timing of dispersal of C. bursa-pastoris to eastern Eurasia indicates that it probably predated introgression. Our results show that polyploid speciation need not result in immediate and complete reproductive isolation, that postpolyploidization hybridization and introgression can contribute significantly to genetic variation in a newly formed polyploid, and that divergence population genetic analysis constitutes a powerful way of testing hypotheses on polyploid speciation.

Molecular Cloning and Characterization of a Moss (Ceratodon purpureus) Nonsymbiotic Hemoglobin Provides Insight into the Early Evolution of Plant Nonsymbiotic Hemoglobins

Garrocho-Villegas, V., Arredondo-Peter, R. — 2008-06-13

Nonsymbiotic hemoglobins (nsHbs) are widespread in plants including bryophytes. Bryophytes (such as mosses) are among the oldest land plants, thus an analysis of a bryophyte nsHb is of interest from an evolutionary perspective. However, very little is known about bryophyte nsHbs. Here, we report the cloning and characterization of an nshb gene (cerhb) from the moss Ceratodon purpureus. Sequence analysis showed that cerhb is interrupted by 3 introns in identical position as all known plant nshb genes, which suggests that the ancestral nshb gene was interrupted by 3 introns. Expression analysis showed that cerhb expresses in protonemas and gametophytes growing in normal conditions and that it overexpresses in protonemas subjected to osmotic (sucrose), heat-shock, cold-, and nitrate-stress conditions. Also, modeling of the Ceratodon nsHb (CerHb) tertiary structure suggests that CerHb is hexacoordinate and that it binds O₂ with high affinity. Comparative analysis of the predicted CerHb with native rice Hb1 and soybean leghemoglobin a structures revealed that the major evolutionary changes that probably occurred during the evolution of plant Hbs were 1) a hexacoordinate to pentacoordinate transition at the heme prosthetic group, 2) a length decrease at the CD-loop and N- and C-termini regions, and 3) the compaction of the protein into a globular structure.

High Rates of Molecular Evolution in Hantaviruses

2008-06-13

Hantaviruses are rodent-borne Bunyaviruses that infect the Arvicolinae, Murinae, and Sigmodontinae subfamilies of Muridae. The rate of molecular evolution in the hantaviruses has been previously estimated at approximately 10^–7 nucleotide substitutions per site, per year (substitutions/site/year), based on the assumption of codivergence and hence shared divergence times with their rodent hosts. If substantiated, this would make the hantaviruses among the slowest evolving of all RNA viruses. However, as hantaviruses replicate with an RNA-dependent RNA polymerase, with error rates in the region of one mutation per genome replication, this low rate of nucleotide substitution is anomalous. Here, we use a Bayesian coalescent approach to estimate the rate of nucleotide substitution from serially sampled gene sequence data for hantaviruses known to infect each of the 3 rodent subfamilies: Araraquara virus (Sigmodontinae), Dobrava virus (Murinae), Puumala virus (Arvicolinae), and Tula virus (Arvicolinae). Our results reveal that hantaviruses exhibit short-term substitution rates of 10^–2 to 10^–4 substitutions/site/year and so are within the range exhibited by other RNA viruses. The disparity between this substitution rate and that estimated assuming rodent–hantavirus codivergence suggests that the codivergence hypothesis may need to be reevaluated.

Molecular Evolution of a Primate-Specific microRNA Family

Zhang, R., Wang, Y.-Q., Su, B. — 2008-06-13

Lineage-specific microRNA (miRNA) families may contribute to developmental novelties during evolution. However, little is known about the origin and evolution of new miRNA families. We report evidence of an Alu-mediated rapid expansion of miRNA genes in a previously identified primate-specific miRNA family, drawn from sequencing and comparative analysis of 9 diverse primate species. Evolutionary analysis reveals similar divergence among miRNA copies whether they are within or between species, lineage-specific gain and loss of miRNAs, and gene pseudolization in multiple species. These observations support a birth-and-death process of miRNA genes in this family, implicating functional diversification during primate evolution. In addition, both secondary structure conservation and reduced single nucleotide polymorphisms density attest to functional constraint of this family in primates. Finally, we observed preferential expression of miRNAs in human placenta and fetal brain, suggesting a functional importance of this family for primate development.

Bayesian Inference of Errors in Ancient DNA Caused by Postmortem Degradation

Mateiu, L. M., Rannala, B. H. — 2008-06-13

Methods for extracting and amplifying sequences using ancient DNA (aDNA) can be prone to errors caused by postmortem modifications of the DNA strand. A new statistical method is developed for predicting errors in aDNA sequences caused by such processes. In addition to the canonical DNA substitution model parameters, a discrete Markov chain is used to describe nucleotide substitutions occurring via postmortem degradation of the aDNA sequences. A computer program, BYPASSR-degr, was developed implementing the method and was used in subsequent analyses of simulated data sets under the new model. Simulation studies show that the new method can be powerful and accurate in identifying damaged sites. The method is applied to analyze aDNA sequences of Etruscans, Adélie penguins, and horses. No significant signals of degradation were observed at any sites of the aDNA sequences we analyzed.

Difficulties in Testing for Covarion-Like Properties of Sequences under the Confounding Influence of Changing Proportions of Variable Sites

Gruenheit, N., Lockhart, P. J., Steel, M., Martin, W. — 2008-06-13

The covarion (COV)-like properties of sequences are poorly described and their impact on phylogenetic analyses poorly understood. We demonstrate using simulations that, under an evolutionary model where the proportion of variable sites changes in nonadjacent lineages, log likelihood values for rates across site (RAS) and COV models become similar, making models difficult to distinguish. Further, although COV and RAS models provide a great improvement in likelihood scores over a homogeneous model with these simulated data, reconstruction accuracy of tree building is low, suggesting caution when it is suspected that proportions of variable sites differ in different evolutionary lineages. We study the performance of a recently developed contingency test that detects the presence of COV-type evolution modified for protein data. We report that if proportions of variable sites (p_var) change in a lineage-specific manner such that their distributions in different lineages become sufficiently nonoverlapping, then the contingency test can incorrectly suggest a homogeneous model. Also of concern is the possibility of different proportions of variable sites between the groups being studied. In a study of chloroplast proteins, interpretation of the test is found to be susceptible to different partitioning of taxon groups, making the test very subjective in its implementation. Extreme intergroup differences in the extent of divergence and difference in proportions of variable sites could be contributing to this effect.