BMC Evolutionary Biology

Chaperonin genes on the rise: new divergent classes and intense duplication in human and other vertebrate genomes.

Publication Date: 2010 Mar 1 PMID: 20193073
Authors: Mukherjee, K. - Conway de Macario, E. - Macario, A. J. - Brocchieri, L.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Chaperonin proteins are well known for the critical role they play in protein folding and in disease. However, the recent identification of three diverged chaperonin paralogs associated with the human Bardet-Biedl and McKusick-Kaufman Syndromes (BBS and MKKS, respectively) indicates that the eukaryotic chaperonin-gene family is larger and more differentiated than previously thought. The availability of complete genome sequences makes possible a definitive characterization of the complete set of chaperonin sequences in human and other species. RESULTS: We identified fifty-four chaperonin-like sequences in the human genome and similar numbers in the genomes of the model organisms mouse and rat. In mammal genomes we identified, besides the well-known CCT chaperonin genes and the three genes associated with the MKKS and BBS pathological conditions, a newly-defined class of chaperonin genes named CCT8L, represented in human by the two sequences CCT8L1 and CCT8L2. Comparative analyses from several vertebrate genomes established the monophyletic origin of chaperonin-like MKKS and BBS genes from the CCT8 lineage. The CCT8L gene originated from a later duplication also in the CCT8 lineage at the onset of mammal evolution and duplicated in primate genomes. The functionality of CCT8L genes in different species was confirmed by evolutionary analyses and in human by expression data. Detailed sequence analysis and structural predictions of MKKS, BBS and CCT8L proteins strongly suggested that they conserve a typical chaperonin-like core structure but that they are unlikely to form a CCT-like oligomeric complex. The characterization of many newly-discovered chaperonin pseudogenes uncovered the intense duplication activity of eukaryotic chaperonin genes. CONCLUSIONS: In vertebrates, chaperonin genes, driven by intense duplication processes, have diversified into multiple classes and functionalities that extend beyond their well-known protein-folding role as part of the typical oligomeric chaperonin complex, emphasizing previous observations on the involvement of individual CCT monomers in microtubule elongation. The functional characterization of newly identified chaperonin genes will be a challenge for future experimental analyses.

post to: CiteULike

Molecular differentiation of commercial varieties and feral populations of oilseed rape (Brassica napus L.).

Publication Date: 2010 Mar 1 PMID: 20193061
Authors: Pascher, K. - Macalka, S. - Rau, D. - Gollmann, G. - Reiner, H. - Glossl, J. - Grabherr, G.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: For assessing the risk of escape of transgenes from cultivation, the persistence of feral populations of crop plants is an important aspect. Feral populations of oilseed rape, Brassica napus, are well known, but only scarce information is available on their population dynamics, particularly in Central Europe. To investigate genetic diversity, origin and persistence of feral oilseed rape in Austria, we compared variation at nine polymorphic microsatellite loci in eight feral populations with 19 commercial varieties. RESULTS: Overall, commercial varieties and feral populations showed a similar pattern of genetic variation and a similar level of observed heterozygosity. The two groups, however, shared less than 50% of the alleles and no multilocus genotype. A significant among-group (commercial varieties versus feral populations) component of genetic variation was observed (AMOVA: FCT = 0.132). Pairwise comparisons between varieties and feral populations showed moderate to very high genetic differentiation (FST = 0.209 - 0.900). The software STRUCTURE also demonstrated a clear separation between commercial varieties and feral samples: out of 17 identified genetic clusters, only one comprised plants from both a commercial variety and feral sites. CONCLUSIONS: The results suggest that feral oilseed rape is able to maintain persistent populations. The feral populations may have derived from older cultivars that were not included in our analyses or perhaps have already hybridised with related crops or wild relatives. Feral populations therefore have to be considered in ecological risk assessment and future coexistence measures as a potential hybridisation partner of transgenic oilseed rape.

post to: CiteULike

Molecular evolution of the crustacean hyperglycemic hormone family in ecdysozoans.

Publication Date: 2010 Feb 25 PMID: 20184761
Authors: Montagne, N. - Desdevises, Y. - Soyez, D. - Toullec, J. Y.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Crustacean Hyperglycemic Hormone (CHH) family peptides are neurohormones known to regulate several important functions in decapod crustaceans such as ionic and energetic metabolism, molting and reproduction. The structural conservation of these peptides, together with the variety of functions they display, led us to investigate their evolutionary history. CHH family peptides exist in insects (Ion Transport Peptides) and may be present in all ecdysozoans as well. In order to extend the evolutionary study to the entire family, CHH family peptides were thus searched in taxa outside decapods, where they have been, to date, poorly investigated. RESULTS: CHH family peptides were characterized by molecular cloning in a branchiopod crustacean, Daphnia magna, and in a collembolan, Folsomia candida. Genes encoding such peptides were also rebuilt in silico from genomic sequences of another branchiopod, a chelicerate and two nematodes. These sequences were included in updated datasets to build phylogenies of the CHH family in pancrustaceans. These phylogenies suggest that peptides found in Branchiopoda and Collembola are more closely related to insect ITPs than to crustacean CHHs. Datasets were also used to support a phylogenetic hypothesis about pancrustacean relationships, which, in addition to gene structures, allowed us to propose two evolutionary scenarios of this multigenic family in ecdysozoans. CONCLUSIONS: Evolutionary scenarios suggest that CHH family genes of ecdysozoans originate from an ancestral two-exon gene, and genes of arthropods from a three-exon one. In malacostracans, the evolution of the CHH family has involved several duplication, insertion or deletion events, leading to neuropeptides with a wide variety of functions, as observed in decapods. This family could thus constitute a promising model to investigate the links between gene duplications and functional divergence.

post to: CiteULike

Identification of shared single copy nuclear genes in Arabidopsis, Populus, Vitis and Oryza and their phylogenetic utility across various taxonomic levels.

Publication Date: 2010 Feb 24 PMID: 20181251
Authors: Duarte, J. M. - Wall, P. K. - Edger, P. P. - Landherr, L. L. - Ma, H. - Pires, J. C. - Leebens-Mack, J. - Depamphilis, C. W.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Although the overwhelming majority of genes found in angiosperms are members of gene families, and both gene- and genome-duplication are pervasive forces in plant genomes, some genes are sufficiently distinct from all other genes in a genome that they can be operationally defined as 'single copy'. Using the gene clustering algorithm MCL-tribe, we have identified a set of 959 single copy genes that are common among the genomes of Arabidopsis thaliana, Populus trichocarpa, Vitis vinifera and Oryza sativa. To characterize these genes, we have performed a number of analyses examining GO annotations, coding sequence length, number of exons, number of domains, presence in distant lineages, such as Selaginella and Physcomitrella, and phylogenetic analysis to estimate copy number throughout seed plants and to demonstrate their phylogenetic utility. We then provide examples of how these genes may be used in phylogenetic analyses to reconstruct organismal history, both by using extant coverage in EST databases for seed plants and de novo amplification via RT-PCR in the family Brassicaceae. RESULTS: There are 959 single copy nuclear genes shared in Arabidopsis, Populus, Vitis and Oryza ["APVO SC genes"]. The majority of these genes are also present in the Selaginella and Physcomitrella genomes. Public EST sets for 197 species suggest that most of these genes are present across a diverse collection of seed plants, and appear to exist as single or very low copy genes, though exceptions are seen in recently polyploid taxa and in lineages where there is significant evidence for a shared large-scale duplication event. Genes encoding proteins localized in organelles are more commonly single copy than expected by chance, but the evolutionary forces responsible for this bias are unknown. Regardless of the evolutionay mechanisms responsible for the large number of shared single copy genes in diverse flowering plant lineages, these genes are valuable for phylogenetic and comparative analyses. Eighteen of the APVO SC single copy genes were amplified in the Brassicaceae using RT-PCR and directly sequenced. Alignments of these sequences provide improved resolution of Brassicaceae phylogeny compared to recent studies using plastid and ITS sequences. An analysis of sequences from 13 APVO SC genes from 69 species of seed plants, derived mainly from public EST databases yielded a phylogeny that was largely congruent with prior hypotheses based on multiple plastid sequences. Whereas single gene phylogenies that rely on EST sequences have limited bootstrap support as the result of limited sequence information, concatenated alignments result in phylogenetic trees with strong bootstrap support for already established relationships. Overall, these single copy nuclear genes are promising markers for phylogenetics, and contain a greater proportion of phylogenetically-informative sites than commonly used protein-coding sequences from the plastid or mitochondrial genomes. CONCLUSIONS: Putatively orthologous, shared single copy nuclear genes provide a vast source of new evidence for plant phylogenetics, genome mapping, and other applications, as well as a substantial class of genes for which functional characterization is needed. Preliminary evidence indicates that many of the shared single copy nuclear genes identified in this study may be well suited as markers for addressing phylogenetic hypotheses at a variety of taxonomic levels.

post to: CiteULike

Microevolutionary dynamics of a macroevolutionary key innovation in a Lepidopteran herbivore.

Publication Date: 2010 Feb 24 PMID: 20181249
Authors: Heidel-Fischer, H. M. - Vogel, H. - Heckel, D. G. - Wheat, C. W.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: A molecular population genetics understanding is central to the study of ecological and evolutionary functional genomics. Population genetics identifies genetic variation and its distribution within and among populations, it reveals the demographic history of the populations studied, and can provide indirect insights into historical selection dynamics. Here we use this approach to examine the demographic and selective dynamics acting of a candidate gene involved in plant-insect interactions. Previous work documents the macroevolutionary and historical ecological importance of the nitrile-specifier protein (Nsp), which facilitated the host shift of Pieridae butterflies onto Brassicales host plants ~80 Myr ago. RESULTS: Here we assess the microevolutionary dynamics of the Nsp gene by studying the within and among-population variation at Nsp and reference genes in the butterfly Pieris rapae (Small Cabbage White). Nsp exhibits unexpectedly high amounts of amino acid polymorphism, unequally distributed across the gene. The vast majority of genetic variation exists within populations, with little to no genetic differentiation among four populations on two continents. A comparison of synonymous and nonsynonymous substitutions in 70 randomly chosen genes among P. rapae and its close relative Pieris brassicae (Large Cabbage White) finds Nsp to have a significantly relaxed functional constraint compared to housekeeping genes. We find strong evidence for a recent population expansion and no role for strong purifying or directional selection upon the Nsp gene. CONCLUSIONS: The microevolutionary dynamics of the Nsp gene in P. rapae are dominated by recent population expansion and variation in functional constraint across the repeated domains of the Nsp gene. While the high amounts of amino acid diversity suggest there may be significant functional differences among allelic variants segregating within populations, indirect tests of selection could not conclusively identify a signature of historical selection. The importance of using this information for planning future studies of potential performance and fitness consequences of the observed variation is discussed.

post to: CiteULike

Exploring the evolutionary dynamics of plasmids: the Acinetobacter pan-plasmidome.

Publication Date: 2010 Feb 24 PMID: 20181243
Authors: Fondi, M. - Bacci, G. - Brilli, M. - Papaleo, C. M. - Mengoni, A. - Vaneechoutte, M. - Dijkshoorn, L. - Fani, R.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Prokaryotic plasmids have a dual importance in the microbial world: first they have a great impact on the metabolic functions of the host cell, providing additional traits that can be accumulated in the cell without altering the gene content of the bacterial chromosome. Additionally and/or alternatively, from a genome perspective, plasmids can provide a basis for genomic rearrangements via homologous recombination and so they can facilitate the loss or acquisition of genes during these events, which eventually may lead to horizontal gene transfer (HGT). Given their importance for conferring adaptive traits to the host organisms, the interest in plasmid sequencing is growing and now many complete plasmid sequences are available online. RESULTS: By using the newly developed Blast2Network bioinformatic tool, a comparative analysis was performed on the plasmid and chromosome sequence data available for bacteria belonging to the genus Acinetobacter, an ubiquitous and clinically important group of gamma-proteobacteria. Data obtained showed that, although most of the plasmids lack mobilization and transfer functions, they have probably a long history of rearrangements with other plasmids and with chromosomes. Indeed, traces of transfers between different species can be disclosed. CONCLUSIONS: We show that, by combining plasmid and chromosome similarity, identity based, network analysis, an evolutionary scenario can be described even for highly mobile genetic elements that lack extensively shared genes. In particular we found that transposases and selective pressure for mercury resistance seem to have played a pivotal role in plasmid evolution in Acinetobacter genomes sequenced so far.

post to: CiteULike

Evolutionary history of anglerfishes (Teleostei: Lophiiformes): a mitogenomic perspective.

Publication Date: 2010 Feb 23 PMID: 20178642
Authors: Miya, M. - Pietsch, T. W. - Orr, J. W. - Arnold, R. J. - Satoh, T. P. - Shedlock, A. M. - Ho, H. C. - Shimazaki, M. - Yabe, M. - Nishida, M.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: The teleost order Lophiiformes, commonly known as the anglerfishes, contains a diverse array of marine fishes, ranging from benthic shallow-water dwellers to highly modified deep-sea midwater species. They comprise 321 living species placed in 68 genera, 18 families and 5 suborders, but approximately half of the species diversity is occupied by deep-sea ceratioids distributed among 11 families. The evolutionary origins of such remarkable habitat and species diversity, however, remain elusive because of the lack of fresh material for a majority of the deep-sea ceratioids and incompleteness of the fossil record across all of the Lophiiformes. To obtain a comprehensive picture of the phylogeny and evolutionary history of the anglerfishes, we assembled whole mitochondrial genome (mitogenome) sequences from 39 lophiiforms (33 newly determined during this study) representing all five suborders and 17 of the 18 families. Sequences of 77 higher teleosts including the 39 lophiiform sequences were unambiguously aligned and subjected to phylogenetic analysis and divergence time estimation. RESULTS: Partitioned maximum likelihood analysis confidently recovered monophyly for all of the higher taxa (including the order itself) with the exception of the Thaumatichthyidae (Lasiognathus was deeply nested within the Oneirodidae). The mitogenomic trees strongly support the most basal and an apical position of the Lophioidei and a clade comprising Chaunacoidei + Ceratioidei, respectively, although alternative phylogenetic positions of the remaining two suborders (Antennarioidei and Ogcocephaloidei) with respect to the above two lineages are statistically indistinguishable. While morphology-based intra-subordinal relationships for relatively shallow, benthic dwellers (Lophioidei, Antennarioidei, Ogcocephaloidei, Chaunacoidei) are either congruent with or statistically indistinguishable from the present mitogenomic tree, those of the principally deep-sea midwater dwellers (Ceratioidei) cannot be reconciled with the molecular phylogeny. A relaxed molecular-clock Bayesian analysis of the divergence times suggests that all of the subordinal diversifications have occurred during a relatively short time period between 100 and 130 Myr ago (early to mid Cretaceous). CONCLUSIONS: The mitogenomic analyses revealed previously unappreciated phylogenetic relationships among the lophiiform suborders and ceratioid familes. Although the latter relationships cannot be reconciled with the earlier hypotheses based on morphology, we found that simple exclusion of the reductive or simplified characters can alleviate some of the conflict. The acquisition of novel features, such as male dwarfism, bioluminescent lures, and unique reproductive modes allowed the deep-sea ceratioids to diversify rapidly in a largely unexploited, food-poor bathypelagic zone (200-2000 m depth) relative to the other lophiiforms occurring in shallow coastal areas.

post to: CiteULike

Evolution of spliceosomal introns following endosymbiotic gene transfer.

Publication Date: 2010 Feb 23 PMID: 20178587
Authors: Ahmadinejad, N. - Dagan, T. - Gruenheit, N. - Martin, W. - Gabaldon, T.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Spliceosomal introns are an ancient, widespread hallmark of eukaryotic genomes. Despite much research, many questions regarding the origin and evolution of spliceosomal introns remain unsolved, partly due to the difficulty of inferring ancestral gene structures. We circumvent this problem by using genes originated by endosymbiotic gene transfer, in which an intron-less structure at the time of the transfer can be assumed. RESULTS: By comparing the exon-intron structures of 64 mitochondrial-derived genes that were transferred to the nucleus at different evolutionary periods, we can trace the history of intron gains in different eukaryotic lineages. Our results show that the intron density of genes transferred relatively recently to the nuclear genome is similar to that of genes originated by more ancient transfers, indicating that gene structure can be rapidly shaped by intron gain after the integration of the gene into the genome and that this process is mainly determined by forces acting specifically on each lineage. We analyze 12 cases of mitochondrial-derived genes that have been transferred to the nucleus independently in more than one lineage. CONCLUSIONS: Remarkably, the proportion of shared intron positions that were gained independently in homologous genes is similar to that proportion observed in genes that were transferred prior to the speciation event and whose shared intron positions might be due to vertical inheritance. A particular case of parallel intron gain in the nad7 gene is discussed in more detail.

post to: CiteULike

Animal HECT ubiquitin ligases: evolution and functional implications.

Publication Date: 2010 Feb 22 PMID: 20175895
Authors: Marin, I.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: HECT ubiquitin ligases (HECT E3s) are key components of the eukaryotic ubiquitin-proteasome system and are involved in the genesis of several human diseases. In this study, I analyze the patterns of diversification of HECT E3s since animals emerged in order to provide the right framework to understand the functional data available for proteins of this family. RESULTS: I show that the current classification of HECT E3s into three groups (NEDD4-like E3s, HERCs and single-HECT E3s) is fundamentally incorrect. First, the existence of a "Single-HECT E3s" group is not supported by phylogenetic analyses. Second, the HERC proteins must be divided into two subfamilies (Large HERCs, Small HERCs) that are evolutionarily very distant, their structural similarity being due to convergence and not to a common origin. Sequence and structural analyses show that animal HECT E3s can be naturally classified into 16 subfamilies. Almost all of them appeared either before animals originated or in early animal evolution. More recently, multiple gene losses have occurred independently in some lineages (nematodes, insects, urochordates), the same groups that have also lost genes of another type of E3s (RBR family). Interestingly, the emergence of some animal HECT E3s precedes the origin of key cellular systems that they regulate (TGF-beta and EGF signal transduction pathways; p53 family of transcription factors) and it can be deduced that distantly related HECT proteins have been independently co-opted to perform similar roles. This may contribute to explain why distantly related HECT E3s are involved in the genesis of multiple types of cancer. CONCLUSIONS: The complex evolutionary history of HECT ubiquitin ligases in animals has been deciphered. The most appropriate model animals to study them and new theoretical and experimental lines of research are suggested by these results.

post to: CiteULike

Retrotransposition and mutation events yield Rap1 GTPases with differential signalling capacity.

Publication Date: 2010 PMID: 20170508
Authors: Zemojtel, T. - Duchniewicz, M. - Zhang, Z. - Paluch, T. - Luz, H. - Penzkofer, T. - Scheele, J. S. - Zwartkruis, F. J.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Retrotransposition of mRNA transcripts gives occasionally rise to functional retrogenes. Through acquiring tempero-spatial expression patterns distinct from their parental genes and/or functional mutations in their coding sequences, such retrogenes may in principle reshape signalling networks. RESULTS: Here we present evidence for such a scenario, involving retrogenes of Rap1 belonging to the Ras family of small GTPases. We identified two murine and one human-specific retrogene of Rap1A and Rap1B, which encode proteins that differ by only a few amino acids from their parental Rap1 proteins. Markedly, human hRap1B-retro and mouse mRap1A-retro1 acquired mutations in the 12th and 59th amino acids, respectively, corresponding to residues mutated in constitutively active oncogenic Ras proteins. Statistical and structural analyses support a functional evolution scenario, where Rap1 isoforms of retrogenic origin are functionally distinct from their parental proteins. Indeed, all retrogene-encoded GTPases have an increased GTP/GDP binding ratio in vivo, indicating that their conformations resemble that of active GTP-bound Rap1. We furthermore demonstrate that these three Rap1 isoforms exhibit distinct affinities for the Ras-binding domain of RalGDS. Finally, when tested for their capacity to induce key cellular processes like integrin-mediated cell adhesion or cell spreading, marked differences are seen. CONCLUSIONS: Together, these data lend strong support for an evolution scenario, where retrotransposition and subsequent mutation events generated species-specific Rap1 isoforms with differential signaling potential. Expression of the constitutively active human Rap1B-retro in cells like those derived from Ramos Burkitt's lymphoma and bone marrow from a patient with myelodysplastic syndrome (MDS) warrants further investigation into its role in disease development.

post to: CiteULike

Two Frequenins in Drosophila: unveiling the evolutionary history of an unusual Neuronal Calcium Sensor (NCS) duplication.

Publication Date: 2010 Feb 19 PMID: 20170488
Authors: Sanchez-Gracia, A. - Romero-Pozuelo, J. - Ferrus, A.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Drosophila Frequenin (Frq), the homolog of the mammalian Neuronal Calcium Sensor-1 (NCS-1), is a high affinity calcium-binding protein with ubiquitous expression in the nervous system. This protein has an important role in the regulation of neurotransmitter release per synapse, axonal growth and button formation. In D. melanogaster, Frequenin is encoded by two genes (frq1 and frq2), a very unexpected feature in the Frq/NCS-1 subfamily. These genes are located in tandem in the same genomic region, and their products are 95% identical in their amino acid sequence, clearly indicating their recent origin by gene duplication. Here, we have investigated the factors involved in this unusual feature by examining the molecular evolution of the two frq genes in Drosophila and the evolutionary dynamics of NCS family in a large set of bilaterian species. RESULTS: Surprisingly, we have found no amino acid replacements fixed across the twelve Drosophila species surveyed. In contrast, synonymous substitutions have been prevalent in the evolution of the coding region of frq1 and frq2, indicating the presence of strong functional constraints following gene duplication. Despite that, we have detected that significant evolutionary rate acceleration had occurred in Frq1 in early times from the duplication, in which positive selection (likely promoting functional diversification) had probably an important role. The analysis of sequence conservation and DNA topology at the non-coding regions of both genes has allowed the identification of DNA regions candidates to be cis-regulatory elements. The results reveal a possible mechanism of regulatory diversification between frq1 and frq2. CONCLUSIONS: The presence of two Frequenins in Drosophila and the rapid accumulation of amino acid substitutions after gene duplication are very unusual features in the evolution of the Frq/NCS-1 subfamily. Here we show that the action of positive selection in concordance with some extent of regulatory diversification might explain these findings. Selected amino acid substitutions in Frq1 likely contributed to the functional divergence between the two duplicates, which, in turn, should have diverged in their regulation by Ecdysone-induced early genes.

post to: CiteULike

Characterization of the tandem CWCH2 sequence motif: a hallmark of inter-zinc finger interactions.

Publication Date: 2010 Feb 19 PMID: 20167128
Authors: Hatayama, M. - Aruga, J.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: The C2H2 zinc finger (ZF) domain is widely conserved among eukaryotic proteins. In Zic/Gli/Zap1 C2H2 ZF proteins, the two N-terminal ZFs form a single structural unit by sharing a hydrophobic core. This structural unit defines a new motif comprised of two tryptophan side chains at the center of the hydrophobic core. Because each tryptophan residue is located between the two cysteine residues of the C2H2 motif, we have named this structure the tandem CWCH2 (tCWCH2) motif. Results: Here, we characterized 587 tCWCH2-containing genes using data derived from public databases. We categorized genes into 11 classes including Zic/Gli/Glis, Arid2/Rsc9, PacC, Mizf, Aebp2, Zap1/ZafA, Fungl, Zfp106, Twincl, Clr1, and Fungl-4ZF, based on sequence similarity, domain organization, and functional similarities. tCWCH2 motifs are mostly found in organisms belonging to the Opisthokonta (metazoa, fungi, and choanoflagellates) and Amoebozoa (amoeba, Dictyostelium discoideum). By comparison, the C2H2 ZF motif is distributed widely among the eukaryotes. The structure and organization of the tCWCH2 motif, its phylogenetic distribution, and molecular phylogenetic analysis suggest that prototypical tCWCH2 genes existed in the Opisthokonta ancestor. Within-group or between-group comparisons of the tCWCH2 amino acid sequence identified three additional sequence features (site-specific amino acid frequencies, longer linker sequence between two C2H2 ZFs, and frequent extra-sequences within C2H2 ZF motifs). Conclusion: These features suggest that the tCWCH2 motif is a specialized motif involved in inter-zinc finger interactions.

post to: CiteULike

Lineage-specific positive selection at the merozoite surface protein 1 (msp1) locus of Plasmodium vivax and related simian malaria parasites.

Publication Date: 2010 PMID: 20167126
Authors: Sawai, H. - Otani, H. - Arisue, N. - Palacpac, N. - de Oliveira Martins, L. - Pathirana, S. - Handunnetti, S. - Kawai, S. - Kishino, H. - Horii, T. - Tanabe, K.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: The 200 kDa merozoite surface protein 1 (MSP-1) of malaria parasites, a strong vaccine candidate, plays a key role during erythrocyte invasion and is a target of host protective immune response. Plasmodium vivax, the most widespread human malaria parasite, is closely related to parasites that infect Asian Old World monkeys, and has been considered to have become a parasite of man by host switch from a macaque malaria parasite. Several Asian monkey parasites have a range of natural hosts. The same parasite species shows different disease manifestations among host species. This suggests that host immune responses to P. vivax-related malaria parasites greatly differ among host species (albeit other factors). It is thus tempting to invoke that a major immune target parasite protein such as MSP-1 underwent unique evolution, depending on parasite species that exhibit difference in host range and host specificity. RESULTS: We performed comparative phylogenetic and population genetic analyses of the gene encoding MSP-1 (msp1) from P. vivax and nine P. vivax-related simian malaria parasites. The inferred phylogenetic tree of msp1 significantly differed from that of the mitochondrial genome, with a striking displacement of P. vivax from a position close to P. cynomolgi in the mitochondrial genome tree to an outlier of Asian monkey parasites. Importantly, positive selection was inferred for two ancestral branches, one leading to P. inui and P. hylobati and the other leading to P. vivax, P. fieldi and P. cynomolgi. This ancestral positive selection was estimated to have occurred three to six million years ago, coinciding with the period of radiation of Asian macaques. Comparisons of msp1 polymorphisms between P. vivax, P. inui and P. cynomolgi revealed that while some positively selected amino acid sites or regions are shared by these parasites, amino acid changes greatly differ, suggesting that diversifying selection is acting species-specifically on msp1. CONCLUSIONS: The present results indicate that the msp1 locus of P. vivax and related parasite species has lineage-specific unique evolutionary history with positive selection. P. vivax and related simian malaria parasites offer an interesting system toward understanding host species-dependent adaptive evolution of immune-target surface antigen genes such as msp1.

post to: CiteULike

Functional evolution of the trace amine associated receptors in mammals and the loss of TAAR1 in dogs.

Publication Date: 2010 Feb 18 PMID: 20167089
Authors: Vallender, E. J. - Xie, Z. - Westmoreland, S. V. - Miller, G. M.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: The trace amine associated receptor family is a diverse array of GPCRs that arose before the first vertebrates walked on land. Trace amine associated receptor 1 (TAAR1) is a wide spectrum aminergic receptor that acts as a modulator in brain monoaminergic systems. Other trace amine associated receptors appear to relate to environmental perception and show a birth-and-death pattern in mammals similar to olfactory receptors. RESULTS: Across mammals, avians, and amphibians, the TAAR1 gene is intact and appears to be under strong purifying selection based on rates of amino acid fixation compared to neutral mutations. We have found that in dogs it has become a pseudogene. Our analyses using a comparative genetics approach revealed that the pseudogenization event predated the emergence of the Canini tribe rather than being coincident with canine domestication. By assessing the effects of the TAAR1 agonist beta-phenylethylamine on [3H]dopamine uptake in canine striatal synaptosomes and comparing the degree and pattern of uptake inhibition to that seen in other mammals, including TAAR1 knockout mice, wild type mice and rhesus monkey, we found that the TAAR1 pseudogenization event resulted in an uncompensated loss of function. CONCLUSION: The gene family has seen expansions among certain mammals, notably rodents, and reductions in others, including primates. By placing the trace amine associated receptors in an evolutionary context we can better understand their function and their potential associations with behavior and neurological disease.

post to: CiteULike

Mitochondrial phylogenomics of the Bivalvia (Mollusca): searching for the origin and mitogenomic correlates of doubly uniparental inheritance of mtDNA.

Publication Date: 2010 Feb 18 PMID: 20167078
Authors: Doucet-Beaupre, H. - Breton, S. - Chapman, E. G. - Blier, P. U. - Bogan, A. E. - Stewart, D. T. - Hoeh, W. R.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Doubly uniparental inheritance (DUI) is an atypical system of animal mtDNA inheritance found only in some bivalves. Under DUI, maternally (F genome) and paternally (M genome) transmitted mtDNAs yield two distinct gender-associated mtDNA lineages. The oldest distinct M and F genomes are found in freshwater mussels (order Unionoida). Comparative analyses of unionoid mitochondrial genomes and a robust phylogenetic framework are necessary to elucidate the origin, function and molecular evolutionary consequences of DUI. Herein, F and M genomes from three unionoid species, Venustaconcha ellipsiformis, Pyganodon grandis and Quadrula quadrula have been sequenced. Comparative genomic analyses were carried out on these six genomes along with two F and one M unionoid genomes from GenBank (F and M genomes of Inversidens japanensis and F genome of Lampsilis ornata). RESULTS: Compared to their unionoid F counterparts, the M genomes contain some unique features including a novel localization of the trnH gene, an inversion of the atp8-trnD genes and a unique 3'coding extension of the cytochrome c oxidase subunit II gene. One or more of these unique M genome features could be causally associated with paternal transmission. Unionoid bivalves are characterized by extreme intraspecific sequence divergences between gender-associated mtDNAs with an average of 50% for V. ellipsiformis, 50% for I. japanensis, 51% for P. grandis and 52% for Q. quadrula (uncorrected amino acid p-distances). Phylogenetic analyses of 12 protein-coding genes from 29 bivalve and five outgroup mt genomes robustly indicate bivalve monophyly and the following branching order within the autolamellibranch bivalves: ((Pteriomorphia, Veneroida) Unionoida). CONCLUSION: The basal nature of the Unionoida within the autolamellibranch bivalves and the previously hypothesized single origin of DUI suggest that (1) DUI arose in the ancestral autolamellibranch bivalve lineage and was subsequently lost in multiple descendant lineages and (2) the mitochondrial genome characteristics observed in unionoid bivalves could more closely resemble the DUI ancestral condition. Descriptions and comparisons presented in this paper are fundamental to a more complete understanding regarding the origins and consequences of DUI.

post to: CiteULike

Gene duplication in the genome of parasitic Giardia lamblia.

Publication Date: 2010 PMID: 20163721
Authors: Sun, J. - Jiang, H. - Flores, R. - Wen, J.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Giardia are a group of widespread intestinal protozoan parasites in a number of vertebrates. Much evidence from G. lamblia indicated they might be the most primitive extant eukaryotes. When and how such a group of the earliest branching unicellular eukaryotes developed the ability to successfully parasitize the latest branching higher eukaryotes (vertebrates) is an intriguing question. Gene duplication has long been thought to be the most common mechanism in the production of primary resources for the origin of evolutionary novelties. In order to parse the evolutionary trajectory of Giardia parasitic lifestyle, here we carried out a genome-wide analysis about gene duplication patterns in G. lamblia. RESULTS: Although genomic comparison showed that in G. lamblia the contents of many fundamental biologic pathways are simplified and the whole genome is very compact, in our study 40% of its genes were identified as duplicated genes. Evolutionary distance analyses of these duplicated genes indicated two rounds of large scale duplication events had occurred in G. lamblia genome. Functional annotation of them further showed that the majority of recent duplicated genes are VSPs (Variant-specific Surface Proteins), which are essential for the successful parasitic life of Giardia in hosts. Based on evolutionary comparison with their hosts, it was found that the rapid expansion of VSPs in G. lamblia is consistent with the evolutionary radiation of placental mammals. CONCLUSIONS: Based on the genome-wide analysis of duplicated genes in G. lamblia, we found that gene duplication was essential for the origin and evolution of Giardia parasitic lifestyle. The recent expansion of VSPs uniquely occurring in G. lamblia is consistent with the increment of its hosts. Therefore we proposed a hypothesis that the increment of Giradia hosts might be the driving force for the rapid expansion of VSPs.

post to: CiteULike

Accelerated microevolution in an outer membrane protein (OMP) of the intracellular bacteria Wolbachia.

Publication Date: 2010 Feb 17 PMID: 20163713
Authors: Baldo, L. - Desjardins, C. A. - Russell, J. A. - Stahlhut, J. K. - Werren, J. H.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Outer membrane proteins (OMPs) of Gram-negative bacteria are key players in the biology of bacterial-host interactions. However, while considerable attention has been given to OMPs of vertebrate pathogens, relatively little is known about the role of these proteins in bacteria that primarily infect invertebrates. One such OMP is found in the intracellular bacteria Wolbachia, which are widespread symbionts of arthropods and filarial nematodes. Recent experimental studies have shown that the Wolbachia surface protein (WSP) can trigger host immune responses and control cell death programming in humans, suggesting a key role of WSP for establishment and persistence of the symbiosis in arthropods. RESULTS: Here we performed an analysis of 515 unique alleles found in 831 Wolbachia isolates, to investigate WSP structure, microevolution and population genetics. WSP shows an eight-strand transmembrane beta-barrel structure with four extracellular loops containing hypervariable regions (HVRs). A clustering approach based upon patterns of HVR haplotype diversity was used to group similar WSP sequences and to estimate the relative contribution of mutation and recombination during early stages of protein divergence. Results indicate that although point mutations generate most of the new protein haplotypes, recombination is a predominant force triggering diversity since the very first steps of protein evolution, causing at least 50% of the total amino acid variation observed in recently diverged proteins. Analysis of synonymous variants indicates that individual WSP protein types are subject to a very rapid turnover and that HVRs can accommodate a virtually unlimited repertoire of peptides. Overall distribution of WSP across hosts supports a non-random association of WSP with the host genus, although extensive horizontal transfer has occurred also in recent times. CONCLUSIONS: In OMPs of vertebrate pathogens, large recombination impact, positive selection, reduced structural and compositional constraints, and extensive lateral gene transfer are considered hallmarks of evolution in response to the adaptive immune system. However, Wolbachia do not infect vertebrates. Here we predict that the rapid turnover of WSP loop motifs could aid in evading or inhibiting the invertebrate innate immune response. Overall, these features identify WSP as a strong candidate for future studies of host-Wolbachia interactions that affect establishment and persistence of this widespread endosymbiosis.

post to: CiteULike

Evolutionary dynamics of U12-type spliceosomal introns.

Publication Date: 2010 PMID: 20163699
Authors: Lin, C. F. - Mount, S. M. - Jarmolowski, A. - Makalowski, W.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Many multicellular eukaryotes have two types of spliceosomes for the removal of introns from messenger RNA precursors. The major (U2) spliceosome processes the vast majority of introns, referred to as U2-type introns, while the minor (U12) spliceosome removes a small fraction (less than 0.5%) of introns, referred to as U12-type introns. U12-type introns have distinct sequence elements and usually occur together in genes with U2-type introns. A phylogenetic distribution of U12-type introns shows that the minor splicing pathway appeared very early in eukaryotic evolution and has been lost repeatedly. RESULTS: We have investigated the evolution of U12-type introns among eighteen metazoan genomes by analyzing orthologous U12-type intron clusters. Examination of gain, loss, and type switching shows that intron type is remarkably conserved among vertebrates. Among 180 intron clusters, only eight show intron loss in any vertebrate species and only five show conversion between the U12 and the U2-type. Although there are only nineteen U12-type introns in Drosophila melanogaster, we found one case of U2 to U12-type conversion, apparently mediated by the activation of cryptic U12 splice sites early in the dipteran lineage. Overall, loss of U12-type introns is more common than conversion to U2-type and the U12 to U2 conversion occurs more frequently among introns of the GT-AG subtype than among introns of the AT-AC subtype. We also found support for natural U12-type introns with non-canonical terminal dinucleotides (CT-AC, GG-AG, and GA-AG) that have not been previously reported. CONCLUSIONS: Although complete loss of the U12-type spliceosome has occurred repeatedly, U12 introns are extremely stable in some taxa, including eutheria. Loss of U12 introns or the genes containing them is more common than conversion to the U2-type. The degeneracy of U12-type terminal dinucleotides among natural U12-type introns is higher than previously thought.

post to: CiteULike

Phenotypic effect of mutations in evolving populations of RNA molecules.

Publication Date: 2010 Feb 17 PMID: 20163698
Authors: Stich, M. - Lazaro, E. - Manrubia, S. C.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: The secondary structure of folded RNA sequences is a good model to map phenotype onto genotype, as represented by the RNA sequence. Computational studies of the evolution of ensembles of RNA molecules towards target secondary structures yield valuable clues to the mechanisms behind adaptation of complex populations. The relationship between the space of sequences and structures, the organization of RNA ensembles at mutation-selection equilibrium, the time of adaptation as a function of the population parameters, the presence of collective effects in quasispecies, or the optimal mutation rates to promote adaptation all are issues that can be explored within this framework. RESULTS: We investigate the effect of microscopic mutations on the phenotype of RNA molecules during their in silico evolution and adaptation. We calculate the distribution of the effects of mutations on fitness, the relative fractions of beneficial and deleterious mutations and the corresponding selection coefficients for populations evolving under different mutation rates. Three different situations are explored: the mutation-selection equilibrium (optimized population) in three different fitness landscapes, the dynamics during adaptation towards a goal structure (adapting population), and the behavior under periodic population bottlenecks (perturbed population). CONCLUSIONS: The ratio between the number of beneficial and deleterious mutations experienced by a population of RNA sequences increases with the value of the mutation rate $\mu$ at which evolution proceeds. In contrast, the selective value of mutations remains almost constant, independent of $\mu$, indicating that adaptation occurs through an increase in the amount of beneficial mutations, with little variations in the average effect they have on fitness. Statistical analyses of the distribution of fitness effects reveal that small effects, either beneficial or deleterious, are well described by a Pareto distribution. These results are robust under changes in the fitness landscape, remarkably when, in addition to selecting a target secondary structure, specific subsequences or low-energy folds are required. A population perturbed by bottlenecks behaves similarly to an adapting population, struggling to return to the optimized state. Whether it can survive in the long run or whether it goes extinct depends critically on the length of the time interval between bottlenecks.

post to: CiteULike

Evolution of plant RNA polymerase IV/V genes: evidence of subneofunctionalization of duplicated NRPD2/NRPE2-like paralogs in Viola (Violaceae).

Publication Date: 2010 Feb 16 PMID: 20158916
Authors: Marcussen, T. - Oxelman, B. - Skog, A. - Jakobsen, K. S.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: DNA-dependent RNA polymerase IV and V (Pol IV and V) are multi-subunit enzymes occurring in plants. The origin of Pol V, specific to angiosperms, from Pol IV, which is present in all land plants, is linked to the duplication of the gene encoding the largest subunit and the subsequent subneofunctionalization of the two paralogs (NRPD1 and NRPE1). Additional duplication of the second-largest subunit, NRPD2/NRPE2, has happened independently in at least some eudicot lineages, but its paralogs are often subject to concerted evolution and gene death and little is known about their evolution nor their affinity with Pol IV and Pol V. RESULTS: We sequenced a ~1500 bp NRPD2/E2-like fragment from 18 Viola species, mostly paleopolyploids, and 6 non-Viola Violaceae species. Incongruence between the NRPD2/E2-like gene phylogeny and species phylogeny indicates a first duplication of NRPD2 relatively basally in Violaceae, with subsequent sorting of paralogs in the descendants, followed by a second duplication in the common ancestor of Viola and Allexis. In Viola, the mutation pattern suggested (sub-) neofunctionalization of the two NRPD2/E2-like paralogs, NRPD2/E2-a and NRPD2/E2-b. The dN/dS ratios indicated that a 54 bp region exerted strong positive selection for both paralogs immediately following duplication. This 54 bp region encodes a domain that is involved in the binding of the Nrpd2 subunit with other Pol IV/V subunits, and may be important for correct recognition of subunits specific to Pol IV and Pol V. Across all Viola taxa 73 NRPD2/E2-like sequences were obtained, of which 23 (32%) were putative pseudogenes - all occurring in polyploids. The NRPD2 duplication was conserved in all lineages except the diploid MELVIO clade, in which NRPD2/E2-b was lost, and its allopolyploid derivates from hybridization with the CHAM clade, section Viola and section Melanium, in which NRPD2/E2-a occurred in multiple copies while NRPD2/E2-b paralogs were either absent or pseudogenized. CONCLUSIONS: Following the relatively recent split of Pol IV and Pol V, our data indicate that these two multi-subunit enzymes are still in the process of specialization and each acquiring fully subfunctionalized copies of their subunit genes. Even after specialization, the NRPD2/E2-like paralogs are prone to pseudogenization and gene conversion and NRPD2 and NRPE2 copy number is a highly dynamic process modulated by allopolyploidy and gene death.

post to: CiteULike

The genetic diversity and evolution of field pea (Pisum) studied by high throughput retrotransposon based insertion polymorphism (RBIP) marker analysis.

Publication Date: 2010 Feb 15 PMID: 20156342
Authors: Jing, R. - Vershinin, A. - Grzebyta, J. - Shaw, P. - Smykal, P. - Marshall, D. - Ambrose, M. J. - Ellis, T. N. - Flavell, A. J.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: The genetic diversity of crop species is the result of natural selection on the wild progenitor and human intervention by ancient and modern farmers and breeders. The genomes of modern cultivars, old cultivated landraces, ecotypes and wild relatives reflect the effects of these forces and provide insights into germplasm structural diversity, the geographical dimension to species diversity and the process of domestication of wild organisms. This issue is also of great practical importance for crop improvement because wild germplasm represents a rich potential source of useful under-exploited alleles or allele combinations. The aim of the present study was to analyse a major Pisum germplasm collection to gain a broad understanding of the diversity and evolution of Pisum and provide a new rational framework for designing germplasm core collections of the genus. RESULTS: 3020 Pisum germplasm samples from the John Innes Pisum germplasm collection were genotyped for 45 retrotransposon based insertion polymorphism (RBIP) markers by the Tagged Array Marker (TAM) method. The data set was stored in a purpose-built Germinate relational database and analysed by both principal coordinate analysis and a nested application of the Structure program which yielded substantially similar but complementary views of the diversity of the genus Pisum. Structure revealed three Groups (1-3) corresponding approximately to landrace, cultivar and wild Pisum respectively, which were resolved by nested Structure analysis into 14 Sub-Groups, many of which correlate with taxonomic sub-divisions of Pisum, domestication related phenotypic traits and/or restricted geographical locations. Genetic distances calculated between these Sub-Groups are broadly supported by principal coordinate analysis and these, together with the trait and geographical data, were used to infer a detailed model for the domestication of Pisum. CONCLUSIONS: These data provide a clear picture of the major distinct gene pools into which the genus Pisum is partitioned and their geographical distribution. The data strongly support the model of independent domestications for P.sativum ssp abyssinicum and P.sativum. The relationships between these two cultivated germplasms and the various sub-divisions of wild Pisum have been clarified and the most likely ancestral wild gene pools for domesticated P. sativum identified. Lastly, this study provides a framework for defining global Pisum germplasm which will be useful for designing core collections.

post to: CiteULike

Protein engineering of conger eel galectins by tracing of molecular evolution using probable ancestral mutants.

Publication Date: 2010 Feb 14 PMID: 20152053
Authors: Konno, A. - Yonemaru, S. - Kitagawa, A. - Muramoto, K. - Shirai, T. - Ogawa, T.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Conger eel galectins, congerin I (ConI) and congerin II (ConII), show the different molecular characteristics resulting from accelerating evolution. We recently reconstructed a probable ancestral form of congerins, Con-anc. It showed properties similar to those of ConII in terms of thermostability and carbohydrate recognition specificity, although it shares a higher sequence similarity with ConI than ConII. RESULTS: In this study, we have focused on the different amino acid residues between Con-anc and ConI, and have performed the protein engineering of Con-anc through site-directed mutagenesis, followed by the molecular evolution analysis of the mutants. This approach revealed the functional importance of loop structures of congerins: (1) N- and C-terminal and loop 5 regions that are involved in conferring a high thermostability to ConI; (2) loops 3, 5, and 6 that are responsible for stronger binding of ConI to most sugars; and (3) loops 5 and 6, and Thr38 residue in loop 3 contribute the specificity of ConI toward lacto-N-fucopentaose-containing sugars. CONCLUSIONS: Thus, this methodology, with tracing of the molecular evolution using ancestral mutants, is a powerful tool for the analysis of not only the molecular evolutionary process, but also the structural elements of a protein responsible for its various functions.

post to: CiteULike

SPRIT: Identifying horizontal gene transfer in rooted phylogenetic trees.

Publication Date: 2010 PMID: 20152048
Authors: Hill, T. - Nordstrom, K. J. - Thollesson, M. - Safstrom, T. M. - Vernersson, A. K. - Fredriksson, R. - Schioth, H. B.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Phylogenetic trees based on sequences from a set of taxa can be incongruent due to horizontal gene transfer (HGT). By identifying the HGT events, we can reconcile the gene trees and derive a taxon tree that adequately represents the species' evolutionary history. One HGT can be represented by a rooted Subtree Prune and Regraft (RSPR) operation and the number of RSPRs separating two trees corresponds to the minimum number of HGT events. Identifying the minimum number of RSPRs separating two trees is NP-hard, but the problem can be reduced to fixed parameter tractable. A number of heuristic and two exact approaches to identifying the minimum number of RSPRs have been proposed. This is the first implementation delivering an exact solution as well as the intermediate trees connecting the input trees. RESULTS: We present the SPR Identification Tool (SPRIT), a novel algorithm that solves the fixed parameter tractable minimum RSPR problem and its GPL licensed Java implementation. The algorithm can be used in two ways, exhaustive search that guarantees the minimum RSPR distance and a heuristic approach that guarantees finding a solution, but not necessarily the minimum one. We benchmarked SPRIT against other software in two different settings, small to medium sized trees i.e. five to one hundred taxa and large trees i.e. thousands of taxa. In the small to medium tree size setting with random artificial incongruence, SPRIT's heuristic mode outperforms the other software by always delivering a solution with a low overestimation of the RSPR distance. In the large tree setting SPRIT compares well to the alternatives when benchmarked on finding a minimum solution within a reasonable time. SPRIT presents both the minimum RSPR distance and the intermediate trees. CONCLUSIONS: When used in exhaustive search mode, SPRIT identifies the minimum number of RSPRs needed to reconcile two incongruent rooted trees. SPRIT also performs quick approximations of the minimum RSPR distance, which are comparable to, and often better than, purely heuristic solutions. Put together, SPRIT is an excellent tool for identification of HGT events and pinpointing which taxa have been involved in HGT.

post to: CiteULike

Comparative analyses reveal distinct sets of lineage-specific genes within Arabidopsis thaliana.

Publication Date: 2010 PMID: 20152032
Authors: Lin, H. - Moghe, G. - Ouyang, S. - Iezzoni, A. - Shiu, S. H. - Gu, X. - Buell, C. R.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: The availability of genome and transcriptome sequences for a number of species permits the identification and characterization of conserved as well as divergent genes such as lineage-specific genes which have no detectable sequence similarity to genes from other lineages. While genes conserved among taxa provide insight into the core processes among species, lineage-specific genes provide insights into evolutionary processes and biological functions that are likely clade or species specific. RESULTS: Comparative analyses using the Arabidopsis thaliana genome and sequences from 178 other species within the Plant Kingdom enabled the identification of 24,624 A. thaliana genes (91.7%) that were termed Evolutionary Conserved (EC) as defined by sequence similarity to a database entry as well as two sets of lineage-specific genes within A. thaliana. One of the A. thaliana lineage-specific gene sets share sequence similarity only to sequences from species within the Brassicaceae family and are termed Conserved Brassicaceae-Specific Genes (914, 3.4%, CBSG). The other set of A. thaliana lineage-specific genes, the Arabidopsis Lineage-Specific Genes (1,324, 4.9%, ALSG), lack sequence similarity to any sequence outside A. thaliana. While many CBSGs (76.7%) and ALSGs (52.9%) are transcribed, the majority of the CBSGs (76.1%) and ALSGs (94.4%) have no annotated function. Co-expression analysis indicated significant enrichment of the CBSGs and ALSGs in multiple functional categories suggesting their involvement in a wide range of biological functions. Subcellular localization prediction revealed that the CBSGs were significantly enriched in proteins targeted to the secretory pathway (412, 45.1%). Among the 107 putatively secreted CBSGs with known functions, 67 encode a putative pollen coat protein or cysteine-rich protein with sequence similarity to the S-locus cysteine-rich protein that is the pollen determinant controlling allele specific pollen rejection in self-incompatible Brassicaceae species. Overall, the ALSGs and CBSGs were more highly methylated in floral tissue compared to the ECs. Single Nucleotide Polymorphism (SNP) analysis showed an elevated ratio of non-synonymous to synonymous SNPs within the ALSGs (1.99) and CBSGs (1.65) relative to the EC set (0.92), mainly caused by an elevated number of non-synonymous SNPs, indicating that they are fast-evolving at the protein sequence level. CONCLUSIONS: Our analyses suggest that while a significant fraction of the A. thaliana proteome is conserved within the Plant Kingdom, evolutionarily distinct sets of genes that may function in defining biological processes unique to these lineages have arisen within the Brassicaceae and A. thaliana.

post to: CiteULike

Structural diversity and evolution of the N-terminal isoform-specific region of ecdysone receptor-A and -B1 isoforms in insects.

Publication Date: 2010 PMID: 20152013
Authors: Watanabe, T. - Takeuchi, H. - Kubo, T.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: The ecdysone receptor (EcR) regulates various cellular responses to ecdysteroids during insect development. Insects have multiple EcR isoforms with different N-terminal A/B domains that contain the isoform-specific activation function (AF)-1 region. Although distinct physiologic functions of the EcR isoforms have been characterized in higher holometabolous insects, they remain unclear in basal direct-developing insects, in which only A isoform has been identified. To examine the structural basis of the EcR isoform-specific AF-1 regions, we performed a comprehensive structural comparison of the isoform-specific region of the EcR-A and -B1 isoforms in insects. RESULTS: The EcR isoforms were newly identified in 51 species of insects and non-insect arthropods, including direct-developing ametabolous and hemimetabolous insects. The comprehensive structural comparison revealed that the isoform-specific region of each EcR isoform contained evolutionally conserved microdomain structures and insect subgroup-specific structural modifications. The A isoform-specific region generally contained four conserved microdomains, including the SUMOylation motif and the nuclear localization signal, whereas the B1 isoform-specific region contained three conserved microdomains, including an acidic activator domain-like motif. In addition, the EcR-B1 isoform of holometabolous insects had a novel microdomain at the N-terminal end. CONCLUSIONS: Given that the nuclear receptor AF-1 is involved in cofactor recruitment and transcriptional regulation, the microdomain structures identified in the isoform-specific A/B domains might function as signature motifs and/or as targets for cofactor proteins that play essential roles in the EcR isoform-specific AF-1 regions. Moreover, the novel microdomain in the isoform-specific region of the holometabolous insect EcR-B1 isoform suggests that the holometabolous insect EcR-B1 acquired additional transcriptional regulation mechanisms.

post to: CiteULike

Positive selection neighboring functionally essential sites and disease-implicated regions of mammalian reproductive proteins.

Publication Date: 2010 PMID: 20149245
Authors: Morgan, C. C. - Loughran, N. B. - Walsh, T. A. - Harrison, A. J. - O'Connell, M. J.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Reproductive proteins are central to the continuation of all mammalian species. The evolution of these proteins has been greatly influenced by environmental pressures induced by pathogens, rival sperm, sexual selection and sexual conflict. Positive selection has been demonstrated in many of these proteins with particular focus on primate lineages. However, the mammalia are a diverse group in terms of mating habits, population sizes and germ line generation times. We have examined the selective pressures at work on a number of novel reproductive proteins across a wide variety of mammalia. RESULTS: We show that selective pressures on reproductive proteins are highly varied. Of the 10 genes analyzed in detail, all contain signatures of positive selection either across specific sites or in specific lineages or a combination of both. Our analysis of SP56 and Col1a1 are entirely novel and the results show positively selected sites present in each gene. Our findings for the Col1a1 gene are suggestive of a link between positive selection and severe disease type. We find evidence in our dataset to suggest that interacting proteins are evolving in symphony: most likely to maintain interacting functionality. CONCLUSION: Our in silico analyses show positively selected sites are occurring near catalytically important regions suggesting selective pressure to maximize efficient fertilization. In those cases where a mechanism of protein function is not fully understood, the sites presented here represent ideal candidates for mutational study. This work has highlighted the widespread rate heterogeneity in mutational rates across the mammalia and specifically has shown that the evolution of reproductive proteins is highly varied depending on the species and interacting partners. We have shown that positive selection and disease are closely linked in the Col1a1 gene.

post to: CiteULike

The zebrafish genome encodes the largest vertebrate repertoire of functional aquaporins with dual paralogy and substrate specificities similar to mammals.

Publication Date: 2010 PMID: 20149227
Authors: Tingaud-Sequeira, A. - Calusinska, M. - Finn, R. N. - Chauvigne, F. - Lozano, J. - Cerda, J.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Aquaporins are integral membrane proteins that facilitate the transport of water and small solutes across cell membranes. These proteins are vital for maintaining water homeostasis in living organisms. In mammals, thirteen aquaporins (AQP0-12) have been characterized, but in lower vertebrates, such as fish, the diversity, structure and substrate specificity of these membrane channel proteins are largely unknown. RESULTS: The screening and isolation of transcripts from the zebrafish (Danio rerio) genome revealed eighteen sequences structurally related to the four subfamilies of tetrapod aquaporins, i.e., aquaporins (AQP0, -1 and -4), water and glycerol transporters or aquaglyceroporins (Glps; AQP3 and AQP7-10), a water and urea transporter (AQP8), and two unorthodox aquaporins (AQP11 and -12). Phylogenetic analyses of nucleotide and deduced amino acid sequences demonstrated dual paralogy between teleost and human aquaporins. Three of the duplicated zebrafish isoforms have unlinked loci, two have linked loci, while DrAqp8 was found in triplicate across two chromosomes. Genomic sequencing, structural analysis, and maximum likelihood reconstruction, further revealed the presence of a putative pseudogene that displays hybrid exons similar to tetrapod AQP5 and -1. Ectopic expression of the cloned transcripts in Xenopus laevis oocytes demonstrated that zebrafish aquaporins and Glps transport water or water, glycerol and urea, respectively, whereas DrAqp11b and -12 were not functional in oocytes. Contrary to humans and some rodents, intrachromosomal duplicates of zebrafish AQP8 were water and urea permeable, while the genomic duplicate only transported water. All aquaporin transcripts were expressed in adult tissues and found to have divergent expression patterns. In some tissues, however, redundant expression of transcripts encoding two duplicated paralogs seems to occur. CONCLUSION: The zebrafish genome encodes the largest repertoire of functional vertebrate aquaporins with dual paralogy to human isoforms. Our data reveal an early and specific diversification of these integral membrane proteins at the root of the crown-clade of Teleostei. Despite the increase in gene copy number, zebrafish aquaporins mostly retain the substrate specificity characteristic of the tetrapod counterparts. Based upon the integration of phylogenetic, genomic and functional data we propose a new classification for the piscine aquaporin superfamily.

post to: CiteULike

Discovery and identification of a male-killing agent in the Japanese ladybird Propylea japonica (Coleoptera: Coccinellidae).

Publication Date: 2010 PMID: 20149223
Authors: Majerus, T. M. - Majerus, M. E.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: Endosymbionts that manipulate the reproduction of their hosts have been reported widely in invertebrates. One such group of endosymbionts is the male-killers. To date all male-killers reported are bacterial in nature, but comprise a diverse group. Ladybirds have been described as a model system for the study of male-killing, which has been reported in multiple species from widespread geographic locations. Whilst criteria of low egg hatch-rate and female-biased progenic sex ratio have been used to identify female hosts of male-killers, variation in vertical transmission efficiency and host genetic factors may result in variation in these phenotypic indicators of male-killer presence. Molecular identification of bacteria and screening for bacterial presence provide us with a more accurate method than breeding data alone to link the presence of the bacteria to the male-killing phenotype. In addition, by identifying the bacteria responsible we may find evidence for horizontal transfer between endosymbiont hosts and can gain insight into the evolutionary origins of male-killing. Phylogenetic placement of male-killing bacteria will allow us to address the question of whether male-killing is a potential strategy for only some, or all, maternally inherited bacteria. Together, phenotypic and molecular characterisation of male-killers will allow a deeper insight into the interactions between host and endosymbiont, which ultimately may lead to an understanding of how male-killers identify and kill male-hosts. RESULTS: A male-killer was detected in the Japanese coccinellid, Propylea japonica (Thunberg) a species not previously known to harbour male-killers. Families produced by female P. japonica showed significantly female-biased sex ratios. One female produced only daughters. This male-killer trait was maternally inherited and antibiotic treatment produced a full, heritable cure. Molecular analysis identified Rickettsia to be associated with the trait in this species of ladybird. CONCLUSION: We conclude that P. japonica is host to a bacterial male-killer that is vertically inherited with variable transmission efficiency. Rickettsia presence correlates with the male-killing trait, but there is some variation in the phenotypic expression of the trait due to interaction with host factors. Phylogenetic analysis using the 16S rRNA and 17 kDa antigen genes suggests there may have been horizontal transfer of Rickettsial male-killers between different ladybird hosts.

post to: CiteULike

A worldwide correlation of lactase persistence phenotype and genotypes.

Publication Date: 2010 Feb 9 PMID: 20144208
Authors: Itan, Y. - Jones, B. L. - Ingram, C. J. - Swallow, D. M. - Thomas, M. G.
Journal: BMC Evol Biol

ABSTRACT: BACKGROUND: The ability of adult humans to digest the milk sugar lactose - lactase persistence - is a dominant Mendelian trait that has been a subject of extensive genetic, medical and evolutionary research. Lactase persistence is common in people of European ancestry as well as some African, Middle Eastern and Southern Asian groups, but is rare or absent elsewhere in the world. The recent identification of independent nucleotide changes that are strongly associated with lactase persistence in different populations worldwide has led to the possibility of genetic tests for the trait. However, it is highly unlikely that all lactase persistence-associated variants are known. Using an extensive database of lactase persistence phenotype frequencies, together with information on how those data were collected and data on the frequencies of lactase persistence variants, we present a global summary of the extent to which current genetic knowledge can explain lactase persistence phenotype frequency. RESULTS: We used surface interpolation of Old World lactase persistence genotype and phenotype frequency estimates obtained from all available literature and perform a comparison between predicted and observed trait frequencies in continuous space. By accommodating additional data on sample numbers and known false negative and false positive rates for the various lactase persistence phenotype tests (blood glucose and breath hydrogen), we also apply a Monte Carlo method to estimate the probability that known lactase persistence-associated allele frequencies can explain observed trait frequencies in different regions. CONCLUSION: Lactase persistence genotype data is currently insufficient to explain lactase persistence phenotype frequency in much of western and southern Africa, southeastern Europe, the Middle East and parts of central and southern Asia. We suggest that further studies of genetic variation in these regions should reveal additional nucleotide variants that are associated with lactase persistence.

post to: CiteULike