Molecular Biology and Evolution - recent issues

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.

The Preferential Retention of Starch Synthesis Genes Reveals the Impact of Whole-Genome Duplication on Grass Evolution

Wu, Y., Zhu, Z., Ma, L., Chen, M. — 2008-05-09

Gene duplication is a major force in evolution. Here, we analyzed the fate of duplicated genes following the ancient whole-genome duplication (WGD) in rice. Polyploidy-derived duplicated genes were found to be preferentially lost from one of each pair of duplicated chromosomal segments, suggesting that the asymmetric gene loss may result from transcriptome dominance of the ancestral allotetraploid genome. Genes involved in synthesis and catabolism of saccharides were found to be preferentially retained in rice, reflecting different trajectories of duplicated genes formed by polyploidy between rice and Arabidopsis. Further studies demonstrated all 3 catalyzing steps in the starch biosynthesis pathway have polyploidy-derived duplicated genes and one copy in each step forms a dominant pathway in the grain filling–stage rice. The new starch biosynthesis pathway reflects one aspect of the impact of WGD on grass evolution.

The McDonald-Kreitman Test and Slightly Deleterious Mutations

Charlesworth, J., Eyre-Walker, A. — 2008-05-09

It is possible to estimate the proportion of substitutions that are due to adaptive evolution using the numbers of silent and nonsilent polymorphisms and substitutions in a McDonald and Kreitman-type analysis. Unfortunately, this estimate of adaptive evolution is biased downward by the segregation of slightly deleterious mutations. It has been suggested that 1 way to cope with the effects of these slightly deleterious mutations is to remove low-frequency polymorphisms from the analysis. We investigate the performance of this method theoretically. We show that although removing low-frequency polymorphisms does indeed reduce the bias in the estimate of adaptive evolution, the estimate is always downwardly biased, often to the extent that one would not be able to detect adaptive evolution, even if it existed. The method is reasonably satisfactory, only if the rate of adaptive evolution is high and the distribution of fitness effects for slightly deleterious mutations is very leptokurtic. Our analysis suggests that adaptive evolution could be quite prevalent in humans (>8%) and still not be detectable using current methodologies. Our analysis also suggests that the level of adaptive evolution has probably been underestimated, possibly substantially, in both bacteria and Drosophila.

Toxin-Resistant Sodium Channels: Parallel Adaptive Evolution across a Complete Gene Family

Jost, M. C., Hillis, D. M., Lu, Y., Kyle, J. W., Fozzard, H. A., Zakon, H. H. — 2008-05-09

Approximately 75% of vertebrate proteins belong to protein families encoded by multiple evolutionarily related genes, a pattern that emerged as a result of gene and genome duplications over the course of vertebrate evolution. In families of genes with similar or related functions, adaptation to a strong selective agent should involve multiple adaptive changes across the entire gene family. However, we know of no evolutionary studies that have explicitly addressed this point. Here, we show how 4 taxonomically diverse species of pufferfishes (Tetraodontidae) each evolved resistance to the guanidinium toxins tetrodotoxin (TTX) and saxitoxin (STX) via parallel amino acid replacements across all 8 sodium channels present in teleost fish genomes. This resulted in diverse suites of coexisting sodium channel types that all confer varying degrees of toxin resistance, yet show remarkable convergence among genes and phylogenetically diverse species. Using site-directed mutagenesis and expression of a vertebrate sodium channel, we also demonstrate that resistance to TTX/STX is enhanced up to 15-fold by single, frequently observed replacements at 2 sites that have not previously been implicated in toxin binding but show similar or identical replacements in pufferfishes and in distantly related vertebrate and nonvertebrate animals. This study presents an example of natural selection acting upon a complete gene family, repeatedly arriving at a diverse but limited number of adaptive changes within the same genome. To be maximally informative, we suggest that future studies of molecular adaptation should consider all functionally similar paralogs of the affected gene family.

Nonadaptive Explanations for Signatures of Partial Selective Sweeps in Drosophila

2008-05-09

A beneficial mutation that has nearly but not yet fixed in a population produces a characteristic haplotype configuration, called a partial selective sweep. Whether nonadaptive processes might generate similar haplotype configurations has not been extensively explored. Here, we consider 5 population genetic data sets taken from regions flanking high-frequency transposable elements in North American strains of Drosophila melanogaster, each of which appears to be consistent with the expectations of a partial selective sweep. We use coalescent simulations to explore whether incorporation of the species' demographic history, purifying selection against the element, or suppression of recombination caused by the element could generate putatively adaptive haplotype configurations. Whereas most of the data sets would be rejected as nonneutral under the standard neutral null model, only the data set for which there is strong external evidence in support of an adaptive transposition appears to be nonneutral under the more complex null model and in particular when demography is taken into account. High-frequency, derived mutations from a recently bottlenecked population, such as we study here, are of great interest to evolutionary genetics in the context of scans for adaptive events; we discuss the broader implications of our findings in this context.

Evolution of the Insulin Receptor Family and Receptor Isoform Expression in Vertebrates

Hernandez-Sanchez, C., Mansilla, A., de Pablo, F., Zardoya, R. — 2008-05-09

The molecular phylogeny of the vertebrate insulin receptor (IR) family was reconstructed under maximum likelihood (ML) to establish homologous relationships among its members. A sister group relationship between the orphan insulin–related receptor (IRR) and the insulin-like growth factor 1 receptor (IGF1R) to the exclusion of the IR obtained maximal bootstrap support. Although both IR and IGF1R were identified in all vertebrates, IRR could not be found in any teleost fish. The ancestral character states at each position of the receptor molecule were inferred for IR, IRR + IGF1R, and all 3 paralogous groups based on the recovered phylogeny using ML in order to determine those residues that could be important for the specific function of IR. For 18 residues, ancestral character state of IR was significantly distinct (probability >0.95) with respect to the corresponding inferred ancestral character states both of IRR + IGF1R and of all 3 vertebrate paralogs. Most of these IR distinct (shared derived) residues were located on the extracellular portion of the receptor (because this portion is larger and the rate of generation of IR shared derived sites is uniform along the receptor), suggesting that functional diversification during the evolutionary history of the family was largely generated modifying ligand affinity rather than signal transduction at the tyrosine kinase domain. In addition, 2 residues at positions 436 and 1095 of the human IR sequence were identified as radical cluster-specific sites in IRR + IGF1R. Both Ir and Irr have an extra exon (namely exon 11) with respect to Igf1r. We used the molecular phylogeny to infer the evolution of this additional exon. The Irr exon 11 can be traced back to amphibians, whereas we show that presence and alternative splicing of Ir exon 11 seems to be restricted exclusively to mammals. The highly divergent sequence of both exons and the reconstructed phylogeny of the vertebrate IR family strongly indicate that both exons were acquired independently by each paralog.

A Mixed Branch Length Model of Heterotachy Improves Phylogenetic Accuracy

Kolaczkowski, B., Thornton, J. W. — 2008-05-09

Evolutionary relationships are typically inferred from molecular sequence data using a statistical model of the evolutionary process. When the model accurately reflects the underlying process, probabilistic phylogenetic methods recover the correct relationships with high accuracy. There is ample evidence, however, that models commonly used today do not adequately reflect real-world evolutionary dynamics. Virtually all contemporary models assume that relatively fast-evolving sites are fast across the entire tree, whereas slower sites always evolve at relatively slower rates. Many molecular sequences, however, exhibit site-specific changes in evolutionary rates, called "heterotachy." Here we examine the accuracy of 2 phylogenetic methods for incorporating heterotachy, the mixed branch length model—which incorporates site-specific rate changes by summing likelihoods over multiple sets of branch lengths on the same tree—and the covarion model, which uses a hidden Markov process to allow sites to switch between variable and invariable as they evolve. Under a variety of simple heterogeneous simulation conditions, the mixed model was dramatically more accurate than homotachous models, which were subject to topological biases as well as biases in branch length estimates. When data were simulated with strong versions of the types of heterotachy observed in real molecular sequences, the mixed branch length model was more accurate than homotachous techniques. Analyses of empirical data sets confirmed that the mixed branch length model can improve phylogenetic accuracy under conditions that cause homotachous models to fail. In contrast, the covarion model did not improve phylogenetic accuracy compared with homotachous models and was sometimes substantially less accurate. We conclude that a mixed branch length approach, although not the solution to all phylogenetic errors, is a valuable strategy for improving the accuracy of inferred trees.

Altered miRNA Repertoire in the Simplified Chordate, Oikopleura dioica

Fu, X., Adamski, M., Thompson, E. M. — 2008-05-09

Recent studies reveal correlation between microRNA (miRNA) innovation and increased developmental complexity. This is exemplified by dramatic expansion of the miRNA inventory in vertebrates, a lineage where genome duplication has played a significant evolutionary role. Urochordates, the closest extant group to the vertebrates, exhibit an opposite trend to genome and morphological simplification. We show that the urochordate, larvacean, Oikopleura dioica, possesses the requisite miRNA biogenic machinery. The miRNAs isolated by small RNA cloning were expressed throughout the short life cycle, a number of which were stocked as maternal determinants prior to rapid embryonic development. We identify sex-specific miRNAs that appeared as male/female gonad differentiation became apparent and were maintained throughout spermatogenesis. Whereas 80% of mammalian miRNAs are hosted in introns of protein-coding genes, the majority of O. dioica miRNA loci were located in antisense orientations to such genes. Including sister group ascidians in analysis of the urochordate miRNA repertoire, we find that 11 highly conserved bilaterian miRNA families have been lost or derived to the point they are not recognizable in urochordates and a further 4 of these families are absent in larvaceans. Subsequent to this loss/derivation, at least 29 novel miRNA families have been acquired in larvaceans. This suggests a profound reorganization of the miRNA repertoire integral to evolution in the urochordate lineage.

Evolution of Gene Expression in the Drosophila Olfactory System

2008-05-09

Host plant shifts by phytophagous insects play a key role in insect evolution and plant ecology. Such shifts often involve major behavioral changes as the insects must acquire an attraction and/or lose the repulsion to the new host plant's odor and taste. The evolution of chemotactic behavior may be due, in part, to gene expression changes in the peripheral sensory system. To test this hypothesis, we compared gene expression in the olfactory organs of Drosophila sechellia, a narrow ecological specialist that feeds on the fruit of Morinda citrifolia, with its close relatives Drosophila simulans and Drosophila melanogaster, which feed on a wide variety of decaying plant matter. Using whole-genome microarrays and quantitative polymerase chain reaction, we surveyed the entire repertoire of Drosophila odorant receptors (ORs) and odorant-binding proteins (OBPs) expressed in the antennae. We found that the evolution of OR and OBP expression was accelerated in D. sechellia compared both with the genome average in that species and with the rate of OR and OBP evolution in the other species. However, some of the gene expression changes that correlate with D. sechellia’s increased sensitivity to Morinda odorants may predate its divergence from D. simulans. Interspecific divergence of olfactory gene expression cannot be fully explained by changes in the relative abundance of different sensilla as some ORs and OBPs have evolved independently of other genes expressed in the same sensilla. A number of OR and OBP genes are upregulated in D. sechellia compared with its generalist relatives. These genes include Or22a, which likely responds to a key odorant of M. citrifolia, and several genes that are yet to be characterized in detail. Increased expression of these genes in D. sechellia may have contributed to the evolution of its unique chemotactic behavior.

Ancient DNA Identification of Early 20th Century Simian T-Cell Leukemia Virus Type 1

2008-05-09

The molecular identification of proviruses from ancient tissues (and particularly from bones) remains a contentious issue. It can be expected that the copy number of proviruses will be low, which magnifies the risk of contamination with retroviruses from exogenous sources. To assess the feasibility of paleoretrovirological studies, we attempted to identify proviruses from early 20th century bones of museum specimens while following a strict ancient DNA methodology. Simian T-cell leukemia virus type 1 sequences were successfully obtained and authenticated from a Chlorocebus pygerythrus specimen. This represents the first clear evidence that it will be possible to use museum specimens to better characterize simian and human T-tropic retrovirus genetic diversity and analyze their origin and evolution, in greater detail.

Adaptive Evolution of Hepcidin Genes in Antarctic Notothenioid Fishes

2008-05-09

Hepcidin is a small bioactive peptide with dual roles as an antimicrobial peptide and as the principal hormonal regulator of iron homeostasis in human and mouse. Hepcidin homologs of very similar structures are found in lower vertebrates, all comprise ~20–25 amino acids with 8 highly conserved cysteines forming 4 intramolecular disulfide bonds, giving hepcidin a hairpin structure. Hepcidins are particularly diverse in teleost fishes, which may be related to the diversity of aquatic environments with varying degree of pathogen challenge, oxygenation, and iron concentration, factors known to alter hepcidin expression in mammals. We characterized the diversity of hepcidin genes of the Antarctic notothenioid fishes that are endemic to the world's coldest and most oxygen-rich marine water. Notothenioid fishes have at least 4 hepcidin variants, in 2 distinctive structural types. Type I hepcidins comprise 3 distinct variants that are homologs of the widespread 8-cysteine hepcidins. Type II is a novel 4-cysteine variant and therefore only 2 possible disulfide bonds, highly expressed in hematopoietic tissues. Analyses of d_N/d_S substitution rate ratios and likelihood ratio test under site-specific models detected significant signal of positive Darwinian selection on the mature hepcidin–coding sequence, suggesting adaptive evolution of notothenioid hepcidins. Genomic polymerase chain reaction and Southern hybridization showed that the novel type II hepcidin occurs exclusively in lineages of the Antarctic notothenioid radiation but not in the basal non-Antarctic taxa, and lineage-specific positive selection was detected on the branch leading to the type II hepcidin clade under branch-site models, suggesting adaptive evolution of the reduced cysteine variant in response to the polar environment. We also isolated a structurally distinct 4-cysteine (4cys) hepcidin from an Antarctic eelpout that is unrelated to the notothenioids but inhabits the same freezing water. Neighbor-Joining (NJ) analyses of teleost hepcidins showed that the eelpout 4cys variant arose independently from the notothenioid version, which lends support to adaptive evolution of reduced cysteine hepcidin variants on cold selection. The NJ tree also showed taxonomic-specific expansions of hepcidin variants, indicating that duplication and diversification of hepcidin genes play important roles in evolutionary response to diverse ecological conditions.

Emergence of Polyproline II-Like Structure at Early Stages of Experimental Evolution from Random Polypeptides

Toyota, H., Hosokawa, M., Urabe, I., Yomo, T. — 2008-05-09

To examine whether a primordial functional protein at the early stages of evolution has structural features, we carried out experimental evolution consisting of 25 cycles (generations) of mutation and selection toward DNA-binding function using a random-sequence polypeptide of 139 amino acid residues with no secondary structure as the initial sequence. In each generation, 16 clones were sampled arbitrarily for sequence analysis, and a phylogenetic tree was constructed. Polypeptide evolution proceeded from the initial point on branch I in 2 main directions of branches II and III. The initial and 2 evolved polypeptides (one at the 24th generation on branch III and the other at the 25th generation on branch II) were purified to examine their functional and structural properties. Although binding of the initial polypeptide to the target DNA was not detected by surface plasmon resonance measurements, the 2 evolved polypeptides bound to the DNA with dissociation constants of 1.6 and 1.0 µM, respectively, indicating an increase in affinity during the experimental evolution. Circular dichroism spectra of the evolved polypeptides, but not of the initial polypeptide, showed features characteristic of the polyproline II (PPII)–like structure, a left-handed helical structure commonly found in natural proteins, suggesting that the structure emerged through the experimental evolution. The same structural feature was found in another experimental evolution toward catalytic activity. These results demonstrate that the PPII-like structure is one of the common features that could have appeared in the early evolutionary stages of primordial functional protein.

On the Expansion of the Pentatricopeptide Repeat Gene Family in Plants

2008-05-09

Pentatricopeptide repeat (PPR) proteins form a huge family in plants (450 members in Arabidopsis and 477 in rice) defined by tandem repetitions of characteristic sequence motifs. Some of these proteins have been shown to play a role in posttranscriptional processes within organelles, and they are thought to be sequence-specific RNA-binding proteins. The origins of this family are obscure as they are lacking from almost all prokaryotes, and the spectacular expansion of the family in land plants is equally enigmatic. In this study, we investigate the growth of the family in plants by undertaking a genome-wide identification and comparison of the PPR genes of 3 organisms: the flowering plants Arabidopsis thaliana and Oryza sativa and the moss Physcomitrella patens. A large majority of the PPR genes in each of the flowering plants are intron less. In contrast, most of the 103 PPR genes in Physcomitrella are intron rich. A phylogenetic comparison of the PPR genes in all 3 species shows similarities between the intron-rich PPR genes in Physcomitrella and the few intron-rich PPR genes in higher plants. Intron-poor PPR genes in all 3 species also display a bias toward a position of their introns at their 5' ends. These results provide compelling evidence that one or more waves of retrotransposition were responsible for the expansion of the PPR gene family in flowering plants. The differing numbers of PPR proteins are highly correlated with differences in organellar RNA editing between the 3 species.

Two Circular Chromosomes of Unequal Copy Number Make Up the Mitochondrial Genome of the Rotifer Brachionus plicatilis

Suga, K., Mark Welch, D. B., Tanaka, Y., Sakakura, Y., Hagiwara, A. — 2008-05-09

The monogonont rotifer Brachionus plicatilis is an emerging model system for a diverse array of questions in limnological ecosystem dynamics, the evolution of sexual recombination, cryptic speciation, and the phylogeny of basal metazoans. We sequenced the complete mitochondrial genome of B. plicatilis sensu strictu NH1L and found that it is composed of 2 circular chromosomes, designated mtDNA-I (11,153 bp) and mtDNA-II (12,672 bp). Hybridization to DNA isolated from mitochondria demonstrated that mtDNA-I is present at 4 times the copy number of mtDNA-II. The only nucleotide similarity between the 2 chromosomes is a 4.9-kbp region of 99.5% identity including a transfer RNA (tRNA) gene and an extensive noncoding region that contains putative D-loop and control sequence. The mtDNA-I chromosome encodes 4 proteins (ATP6, COB, NAD1, and NAD2), 13 tRNAs, and the large and small subunit ribosomal RNAs; mtDNA-II encodes 8 proteins (COX1–3, NAD3–6, and NAD4L) and 9 tRNAs. Gene order is not conserved between B. plicatilis and its closest relative with a sequenced mitochondrial genome, the acanthocephalan Leptorhynchoides thecatus, or other sequenced mitochondrial genomes. Polymerase chain reaction assays and Southern hybridization to DNA from 18 strains of Brachionus suggest that the 2-chromosome structure has been stable for millions of years. The novel organization of the B. plicatilis mitochondrial genome into 2 nearly equal chromosomes of 4-fold different copy number may provide insight into the evolution of metazoan mitochondria and the phylogenetics of rotifers and other basal animal phyla.

Evolutionary Genomics of Host Adaptation in Vesicular Stomatitis Virus

Remold, S. K., Rambaut, A., Turner, P. E. — 2008-05-09

Populations experiencing similar selection pressures can sometimes diverge in the genetic architectures underlying evolved complex traits. We used RNA virus populations of large size and high mutation rate to study the impact of historical environment on genome evolution, thus increasing our ability to detect repeatable patterns in the evolution of genetic architecture. Experimental vesicular stomatitis virus populations were evolved on HeLa cells, on MDCK cells, or on alternating hosts. Turner and Elena (2000. Cost of host radiation in an RNA virus. Genetics. 156:1465–1470.) previously showed that virus populations evolved in single-host environments achieved high fitness on their selected hosts but failed to increase in fitness relative to their ancestor on the unselected host and that alternating-host–evolved populations had high fitness on both hosts. Here we determined the complete consensus sequence for each evolved population after 95 generations to gauge whether the parallel phenotypic changes were associated with parallel genomic changes. We also analyzed the patterns of allele substitutions to discern whether differences in fitness across hosts arose through true pleiotropy or the presence of not only a mutation that is beneficial in both hosts but also 1 or more mutations at other loci that are costly in the unselected environment (mutation accumulation [MA]). We found that ecological history may influence to what extent pleiotropy and MA contribute to fitness asymmetries across environments. We discuss the degree to which current genetic architecture is expected to constrain future evolution of complex traits, such as host use by RNA viruses.

Mobility Pathways for Vertebrate L1, L2, CR1, and RTE Clade Retrotransposons

Ichiyanagi, K., Okada, N. — 2008-05-09

Autonomous non–long terminal repeat retrotransposons (NLRs) are ubiquitous mobile genetic elements that insert their DNA copies at new locations by retrotransposition. In vertebrates, there are 4 NLR clades, L1, L2, CR1, and RTE, which diverged in the Precambrian era. It has been demonstrated that retrotransposition of L1 and L2 members proceeds via coordinated reactions of targeted DNA cleavage and reverse transcription catalyzed by the NLR-encoded proteins, which are followed by the joining of the 5' (upstream) junction. However, the study on the mobility pathways for vertebrate NLRs is so far limited to L1 and L2. In this report, using target analysis of nested transposons for genomic copies, we studied retrotransposition pathways for a variety of vertebrate NLRs, including those of the L1, L2, CR1, and RTE clades in the human, cow, opossum, chicken, and zebrafish genomes. Thus, this study constitutes the first comprehensive analysis of NLR retrotransposition products in vertebrates. Our data revealed that these elements share similar mechanisms for the cleavages of the 2 target DNA strands and for the initiation of reverse transcription. Possible endonuclease-independent insertions were also identified. Overall, our results suggest the existence of multiple retrotransposition pathways that are conserved among the diverse NLR clades in various vertebrate hosts.

The cbb3 Oxidases Are an Ancient Innovation of the Domain Bacteria

Ducluzeau, A.-L., Ouchane, S., Nitschke, W. — 2008-05-09

A survey of genomes for the presence of gene clusters related to cbb₃ oxidases detected bona fide members of the family in almost all phyla of the domain Bacteria. No archaeal representatives were found. The subunit composition was seen to vary substantially between clades observed on the phylogenetic tree of the catalytic subunit CcoN. The protein diade formed by CcoN and the monoheme cytochrome CcoO appears to constitute the functionally essential "core" of the enzyme conserved in all sampled cbb₃ gene clusters. The topology of the phylogenetic tree contradicts the scenario of a recent origin of cbb₃ oxidases and substantiates the status of this family as a phylogenetic entity on the same level as the other subgroups of the heme-copper superfamily (including nitric oxide reductase). This finding resuscitates and exacerbates the conundrum of the evolutionary origin of heme-copper oxidases.

A Cryptic Algal Group Unveiled: A Plastid Biosynthesis Pathway in the Oyster Parasite Perkinsus marinus

Matsuzaki, M., Kuroiwa, H., Kuroiwa, T., Kita, K., Nozaki, H. — 2008-05-09

Plastids are widespread in plant and algal lineages. They are also exploited by some nonphotosynthetic protists, including malarial parasites, to support their diverse modes of life. However, cryptic plastids may exist in other nonphotosynthetic protists, which could be important in studies on the diversity and evolution of plastids. The parasite Perkinsus marinus, which causes mass mortality in oyster farms, is a nonphotosynthetic protist that is phylogenetically related to plastid-bearing dinoflagellates and apicomplexans. In this study, we searched for P. marinus methylerythritol phosphate (MEP) pathway genes, responsible for de novo isoprenoid synthesis in plastids, and determined the full-length gene sequences for 6 of 7 of these genes. Phylogenetic analyses revealed that each P. marinus gene clusters with orthologs from plastid-bearing eukaryotes, which have MEP pathway genes with essentially the same mosaic pattern of evolutionary origin. A new analytical method called sliding-window iteration of TargetP was developed to examine the distribution of targeting preferences. This analysis revealed that the sequenced genes encode bipartite targeting peptides that are characteristic of proteins targeted to secondary plastids originating from endosymbiosis of eukaryotic algae. These results support our idea that Perkinsus is a cryptic algal group containing nonphotosynthetic secondary plastids. In fact, immunofluorescent microscopy indicated that 1 of the MEP pathway enzymes, 1-deoxy-D-xylulose 5-phosphate reductoisomerase, was localized to small compartments near mitochondrion, which are possibly plastids. This tiny organelle seems to contain very low quantities of DNA or may even lack DNA entirely. The MEP pathway genes are a useful tool for investigating plastid evolution in both of the photosynthetic and nonphotosynthetic eukaryotes and led us to propose the hypothesis that ancestral "chromalveolates" harbored plastids before a secondary endosymbiotic event.

Evolutionary Patterns of MHC Class II B in Owls and Their Implications for the Understanding of Avian MHC Evolution

Burri, R., Hirzel, H. N., Salamin, N., Roulin, A., Fumagalli, L. — 2008-05-09

Owing to its special mode of evolution and central role in the adaptive immune system, the major histocompatibility complex (MHC) has become the focus of diverse disciplines such as immunology, evolutionary ecology, and molecular evolution. MHC evolution has been studied extensively in diverse vertebrate lineages over the last few decades, and it has been suggested that birds differ from the established mammalian norm. Mammalian MHC genes evolve independently, and duplication history (i.e., orthology) can usually be traced back within lineages. In birds, this has been observed in only 3 pairs of closely related species. Here we report strong evidence for the persistence of orthology of MHC genes throughout an entire avian order. Phylogenetic reconstructions of MHC class II B genes in 14 species of owls trace back orthology over tens of thousands of years in exon 3. Moreover, exon 2 sequences from several species show closer relationships than sequences within species, resembling transspecies evolution typically observed in mammals. Thus, although previous studies suggested that long-term evolutionary dynamics of the avian MHC was characterized by high rates of concerted evolution, resulting in rapid masking of orthology, our results question the generality of this conclusion. The owl MHC thus opens new perspectives for a more comprehensive understanding of avian MHC evolution.

Origins of Human Malaria: Rare Genomic Changes and Full Mitochondrial Genomes Confirm the Relationship of Plasmodium falciparum to Other Mammalian Parasites but Complicate the Origins of Plasmodium vivax

Roy, S. W., Irimia, M. — 2008-05-09

Despite substantial work, the phylogeny of malaria parasites remains debated. The matter is complicated by concerns about patterns of evolution in potentially strongly selected genes as well as the extreme AT bias of some Plasmodium genomes. Particularly contentious has been the position of the most virulent human parasite Plasmodium falciparum, whether grouped with avian parasites or within a larger clade of mammalian parasites. Here, we study 3 classes of rare genomic changes, as well as the sequences of mitochondrial ribosomal RNA (rRNA) genes. We report 3 lines of support for a clade of mammalian parasites: 1) we find no instances of spliceosomal intron loss in a hypothetical ancestor of P. falciparum and the avian parasite Plasmodium gallinaceum, suggesting against a close relationship between those species; 2) we find 4 genomic mitochondrial indels supporting a mammalian clade, but none grouping P. falciparum with avian parasites; and 3) slowly evolving mitochondrial rRNA sequences support a mammalian parasite clade with 100% posterior probability. We further report a large deletion in the mitochondrial large subunit rRNA gene, which suggests a subclade including both African and Asian parasites within the clade of closely related primate malarias. This contrasts with previous studies that provided strong support for separate Asian and African clades, and reduces certainty about the historical and geographic origins of Plasmodium vivax. Finally, we find a lack of synapomorphic gene losses, suggesting a low rate of ancestral gene loss in Plasmodium.