Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified

The organelle genomes of Hassawi rice (Oryza sativa L.) and its hybrid in saudi arabia: genome variation, rearrangement, and origins

Hassawi rice (Oryza sativa L.) is a landrace adapted to the climate of Saudi Arabia, characterized by its strong resistance to soil salinity and drought. Using high quality sequencing reads extracted from raw data of a whole genome sequencing project, we assembled both chloroplast (cp) and mitochondrial (mt) genomes of the wild-type Hassawi rice (Hassawi-1) and its dwarf hybrid (Hassawi-2). We discovered 16 InDels (insertions and deletions) but no SNP (single nucleotide polymorphism) is present between the two Hassawi cp genomes. We identified 48 InDels and 26 SNPs in the two Hassawi mt genomes and a new type of sequence variation, termed reverse complementary variation (RCV) in the rice cp genomes. There are two and four RCVs identified in Hassawi-1 when compared to 93–11 (indica) and Nipponbare (japonica), respectively. Microsatellite sequence analysis showed there are more SSRs in the genic regions of both cp and mt genomes in the Hassawi rice than in the other rice varieties. There are also large repeats in the Hassawi mt genomes, with the longest length of 96,168 bp and 96,165 bp in Hassawi-1 and Hassawi-2, respectively. We believe that frequent DNA rearrangement in the Hassawi mt and cp genomes indicate ongoing dynamic processes to reach genetic stability under strong environmental pressures. Based on sequence variation analysis and the breeding history, we suggest that both Hassawi-1 and Hassawi-2 originated from the Indonesian variety Peta since genetic diversity between the two Hassawi cultivars is very low albeit an unknown historic origin of the wild-type Hassawi rice.

Phylogenetic relationships of Palaearctic Formica species (Hymenoptera, Formicidae) based on mitochondrial cytochrome B sequences

Ants of genus Formica demonstrate variation in social organization and represent model species for ecological, behavioral, evolutionary studies and testing theoretical implications of the kin selection theory. Subgeneric division of the Formica ants based on morphology has been questioned and remained unclear after an allozyme study on genetic differentiation between 13 species representing all subgenera was conducted. In the present study, the phylogenetic relationships within the genus were examined using mitochondrial DNA sequences of the cytochrome b and a part of the NADH dehydrogenase subunit 6. All 23 Formica species sampled in the Palaearctic clustered according to the subgeneric affiliation except F. uralensis that formed a separate phylogenetic group. Unlike Coptoformica and Formica s. str., the subgenus Serviformica did not form a tight cluster but more likely consisted of a few small clades. The genetic distances between the subgenera were around 10%, implying approximate divergence time of 5 Myr if we used the conventional insect divergence rate of 2% per Myr. Within-subgenus divergence estimates were 6.69% in Serviformica, 3.61% in Coptoformica, 1.18% in Formica s. str., which supported our previous results on relatively rapid speciation in the latter subgenus. The phylogeny inferred from DNA sequences provides a necessary framework against which the evolution of social traits can be compared. We discuss implications of inferred phylogeny for the evolution of social traits.

Characterization of pro-inflammatory flagellin proteins produced by Lactobacillus ruminis and related motile Lactobacilli

Lactobacillus ruminis is one of at least twelve motile but poorly characterized species found in the genus Lactobacillus. Of these, only L. ruminis has been isolated from mammals, and this species may be considered as an autochthonous member of the gastrointestinal microbiota of humans, pigs and cows. Nine L. ruminis strains were investigated here to elucidate the biochemistry and genetics of Lactobacillus motility. Six strains isolated from humans were non-motile while three bovine isolates were motile. A complete set of flagellum biogenesis genes was annotated in the sequenced genomes of two strains, ATCC25644 (human isolate) and ATCC27782 (bovine isolate), but only the latter strain produced flagella. Comparison of the L. ruminis and L. mali DSM20444(T) motility loci showed that their genetic content and gene-order were broadly similar, although the L. mali motility locus was interrupted by an 11.8 Kb region encoding rhamnose utilization genes that is absent from the L. ruminis motility locus. Phylogenetic analysis of 39 motile bacteria indicated that Lactobacillus motility genes were most closely related to those of motile carnobacteria and enterococci. Transcriptome analysis revealed that motility genes were transcribed at a significantly higher level in motile L. ruminis ATCC27782 than in non-motile ATCC25644. Flagellin proteins were isolated from L. ruminis ATCC27782 and from three other Lactobacillus species, while recombinant flagellin of aflagellate L. ruminis ATCC25644 was expressed and purified from E. coli. These native and recombinant Lactobacillus flagellins, and also flagellate L. ruminis cells, triggered interleukin-8 production in cultured human intestinal epithelial cells in a manner suppressed by short interfering RNA directed against Toll-Like Receptor 5. This study provides genetic, transcriptomic, phylogenetic and immunological insights into the trait of flagellum-mediated motility in the lactobacilli.

Molecular phylogeny, classification and evolution of conopeptides

Conopeptides are toxins expressed in the venom duct of cone snails (Conoidea, Conus). These are mostly well-structured peptides and mini-proteins with high potency and selectivity for a broad range of cellular targets. In view of these properties, they are widely used as pharmacological tools and many are candidates for innovative drugs. The conopeptides are primarily classified into superfamilies according to their peptide signal sequence, a classification that is thought to reflect the evolution of the multigenic system. However, this hypothesis has never been thoroughly tested. Here we present a phylogenetic analysis of 1,364 conopeptide signal sequences extracted from GenBank. The results validate the current conopeptide superfamily classification, but also reveal several important new features. The so-called "cysteine-poor" conopeptides are revealed to be closely related to "cysteine-rich" conopeptides; with some of them sharing very similar signal sequences, suggesting that a distinction based on cysteine content and configuration is not phylogenetically relevant and does not reflect the evolutionary history of conopeptides. A given cysteine pattern or pharmacological activity can be found across different superfamilies. Furthermore, a few conopeptides from GenBank do not cluster in any of the known superfamilies, and could represent yet-undefined superfamilies. A clear phylogenetically based classification should help to disentangle the diversity of conopeptides, and could also serve as a rationale to understand the evolution of the toxins in the numerous other species of conoideans and venomous animals at large.

Macroevolution of venom apparatus innovations in auger snails (Gastropoda; Conoidea; Terebridae)

The Terebridae are a diverse family of tropical and subtropical marine gastropods that use a complex and modular venom apparatus to produce toxins that capture polychaete and enteropneust preys. The complexity of the terebrid venom apparatus suggests that venom apparatus development in the Terebridae could be linked to the diversification of the group and can be analyzed within a molecular phylogenetic scaffold to better understand terebrid evolution. Presented here is a molecular phylogeny of 89 terebrid species belonging to 12 of the 15 currently accepted genera, based on Bayesian inference and Maximum Likelihood analyses of amplicons of 3 mitochondrial (COI, 16S and 12S) and one nuclear (28S) genes. The evolution of the anatomy of the terebrid venom apparatus was assessed by mapping traits of six related characters: proboscis, venom gland, odontophore, accessory proboscis structure, radula, and salivary glands. A novel result concerning terebrid phylogeny was the discovery of a previously unrecognized lineage, which includes species of Euterebra and Duplicaria. The non-monophyly of most terebrid genera analyzed indicates that the current genus-level classification of the group is plagued with homoplasy and requires further taxonomic investigations. Foregut anatomy in the family Terebridae reveals an inordinate diversity of features that covers the range of variability within the entire superfamily Conoidea, and that hypodermic radulae have likely evolved independently on at least three occasions. These findings illustrate that terebrid venom apparatus evolution is not perfunctory, and involves independent and numerous changes of central features in the foregut anatomy. The multiple emergence of hypodermic marginal radular teeth in terebrids are presumably associated with variable functionalities, suggesting that terebrids have adapted to dietary changes that may have resulted from predator-prey relationships. The anatomical and phylogenetic results presented serve as a starting point to advance investigations about the role of predator-prey interactions in the diversification of the Terebridae and the impact on their peptide toxins, which are promising bioactive compounds for biomedical research and therapeutic drug development.

Genomic and morphological evidence converge to resolve the enigma of Strepsiptera

The phylogeny of insects, one of the most spectacular radiations of life on earth, has received considerable attention. However, the evolutionary roots of one intriguing group of insects, the twisted-wing parasites (Strepsiptera), remain unclear despite centuries of study and debate. Strepsiptera exhibit exceptional larval developmental features, consistent with a predicted step from direct (hemimetabolous) larval development to complete metamorphosis that could have set the stage for the spectacular radiation of metamorphic (holometabolous) insects. Here we report the sequencing of a Strepsiptera genome and show that the analysis of sequence-based genomic data (comprising more than 18 million nucleotides from nearly 4,500 genes obtained from a total of 13 insect genomes), along with genomic metacharacters, clarifies the phylogenetic origin of Strepsiptera and sheds light on the evolution of holometabolous insect development. Our results provide overwhelming support for Strepsiptera as the closest living relatives of beetles (Coleoptera). They demonstrate that the larval developmental features of Strepsiptera, reminiscent of those of hemimetabolous insects, are the result of convergence. Our analyses solve the long-standing enigma of the evolutionary roots of Strepsiptera and reveal that the holometabolous mode of insect development is more malleable than previously thought.

Ancient gene transfer from algae to animals: mechanisms and evolutionary significance

Background

Horizontal gene transfer (HGT) is traditionally considered to be rare in multicellular eukaryotes such as animals. Recently, many genes of miscellaneous algal origins were discovered in choanoflagellates. Considering that choanoflagellates are the existing closest relatives of animals, we speculated that ancient HGT might have occurred in the unicellular ancestor of animals and affected the long-term evolution of animals.

Results

Through genome screening, phylogenetic and domain analyses, we identified 14 gene families, including 92 genes, in the tunicate Ciona intestinalis that are likely derived from miscellaneous photosynthetic eukaryotes. Almost all of these gene families are distributed in diverse animals, suggesting that they were mostly acquired by the common ancestor of animals. Their miscellaneous origins also suggest that these genes are not derived from a particular algal endosymbiont. In addition, most genes identified in our analyses are functionally related to molecule transport, cellular regulation and methylation signaling, suggesting that the acquisition of these genes might have facilitated the intercellular communication in the ancestral animal.

Conclusions

Our findings provide additional evidence that algal genes in aplastidic eukaryotes are not exclusively derived from historical plastids and thus important for interpreting the evolution of eukaryotic photosynthesis. Most importantly, our data represent the first evidence that more anciently acquired genes might exist in animals and that ancient HGT events have played an important role in animal evolution.

Molecular characterization of an adaptive response to alkylating agents in the opportunistic pathogen Aspergillus fumigatus

An adaptive response to alkylating agents based upon the conformational change of a methylphosphotriester (MPT) DNA repair protein to a transcriptional activator has been demonstrated in a number of bacterial species, but this mechanism appears largely absent from eukaryotes. Here, we demonstrate that the human pathogen Aspergillus fumigatus elicits an adaptive response to sub-lethal doses of the mono-functional alkylating agent N-methyl-N'-nitro-N-nitrosoguanidine (MNNG). We have identified genes that encode MPT and O(6)-alkylguanine DNA alkyltransferase (AGT) DNA repair proteins; deletions of either of these genes abolish the adaptive response and sensitize the organism to MNNG. In vitro DNA repair assays confirm the ability of MPT and AGT to repair methylphosphotriester and O(6)-methylguanine lesions respectively. In eukaryotes, the MPT protein is confined to a select group of fungal species, some of which are major mammalian and plant pathogens. The evolutionary origin of the adaptive response is bacterial and rooted within the Firmicutes phylum. Inter-kingdom horizontal gene transfer between Firmicutes and Ascomycete ancestors introduced the adaptive response into the Fungal kingdom. Our data constitute the first detailed characterization of the molecular mechanism of the adaptive response in a lower eukaryote and has applications for development of novel fungal therapeutics targeting this DNA repair system.

Mitochondrial genomes from the lichenized fungi Peltigera membranacea and Peltigera malacea: features and phylogeny

Mitochondrial genomes from the fungal partners of two terricolous foliose lichen symbioses, Peltigera membranacea and Peltigera malacea, have been determined using metagenomic approaches, including RNA-seq. The roughly 63 kb genomes show all the major features found in other Pezizomycotina, such as unidirectional transcription, 14 conserved protein genes, genes for the two subunit rRNAs and for a set of 26 tRNAs used in translating the 62 amino acid codons. In one of the tRNAs a CAU anticodon is proposed to be modified, via the action of the nuclear-encoded enzyme, tRNA Ile lysidine synthase, so that it recognizes the codon AUA (Ile) instead of AUG (Met). The overall arrangements and sequences of the two circular genomes are similar, the major difference being the inversion and deterioration of a gene encoding a type B DNA polymerase. Both genomes encode the RNA component of RNAse P, a feature seldom found in ascomycetes. The difference in genome size from the minimal ascomycete mitochondrial genomes is largely due to 17 and 20 group I introns, respectively, most associated with homing endonucleases and all found within protein-coding genes and the gene encoding the large subunit rRNA. One new intron insertion point was found, and an unusually small exon of seven nucleotides (nt) was identified and verified by RNA sequencing. Comparative analysis of mitochondrion-encoded proteins places the Peltigera spp., representatives of the class Lecanoromycetes, close to Leotiomycetes, Dothidiomycetes, and Sordariomycetes, in contrast to phylogenies found using nuclear genes.

Diversity at the major histocompatibility complex Class II in the platypus, Ornithorhynchus anatinus

The platypus (Ornithorhynchus anatinus) is the sole survivor of a previously widely distributed and diverse lineage of ornithorhynchid monotremes. Its dependence on healthy water systems imposes an inherent sensitivity to habitat degradation and climate change. Here, we compare genetic diversity at the major histocompatibility complex (MHC) Class II-DZB gene and 3 MHC-associated microsatellite markers with diversity at 6 neutral microsatellite markers in 70 platypuses from across their range, including the mainland of Australia and the isolated populations of Tasmania, King Island, and Kangaroo Island. Overall, high DZB diversity was observed in the platypus, with 57 DZB β1 alleles characterized. Significant positive selection was detected within the DZB peptide-binding region, promoting variation in this domain. Low levels of genetic diversity were detected at all markers in the 2 island populations, King Island (endemic) and Kangaroo Island (introduced), with the King Island platypuses monomorphic at the DZB locus. Loss of MHC diversity on King Island is of concern, as the population may have compromised immunological fitness and reduced ability to resist changing environmental conditions.

Phylogenetic analysis of ancient DNA using BEAST

Under exceptional circumstances, it is possible to obtain DNA sequences from samples that are up to hundreds of thousands of years old. These data provide an opportunity to look directly at past genetic diversity, to trace the evolutionary process through time, and to infer demographic and phylogeographic trends. Ancient DNA (aDNA) data sets have some degree of intrinsic temporal structure because the sequences have been obtained from samples of different ages. When analyzing these data sets, it is usually necessary to take the sampling times into account. A number of phylogenetic methods have been designed with this purpose in mind. Here I describe the steps involved in Bayesian phylogenetic analysis of aDNA data. I outline a procedure that can be used to co-estimate the genealogical relationships, mutation rate, evolutionary timescale, and demographic history of the study species in a single analytical framework. A number of modifications to the methodology can be made in order to deal with complicating factors such as postmortem damage, sequences from undated samples, and data sets with low information content.

Identification and characterization of bocaviruses in cats and dogs reveals a novel feline bocavirus and a novel genetic group of canine bocavirus

We report the identification and genome characterization of a novel bocavirus, feline bocavirus (FBoV), and novel bocaviruses closely related to canine bocavirus (CBoV) strain Con-161 in stray cats and dogs in Hong Kong, respectively. FBoV was detected by PCR in 7.2, 0.3, 1.6, 2.0 and 0.8% of faecal, nasal, urine, kidney and blood samples, respectively, from 364 cats, while CBoV was detected in 4.6, 5.1, 6.3 and 0.3% of faecal, nasal, urine and blood samples, respectively, from 351 dogs. Three FBoV genomes sequenced revealed the presence of three ORFs characteristic of bocaviruses. Phylogenetic analysis showed that FBoVs were related only distantly to other bocaviruses, forming a distinct cluster within the genus, with ≤ 5.7% nucleotide identities to the genome of minute virus of canines. The four CBoV genomes sequenced shared 87.4-89.2% nucleotide identities with that of CBoV strain Con-161. In addition to the three bocavirus ORFs, they encoded an additional ORF, ORF4, immediately downstream of the ORF for non-structural protein 1 (NS1), which was not found in other bocaviruses including CBoV strain Con-161. They also possessed a putative second exon encoding the C-terminal region of NS1 and conserved RNA-splicing signals, previously described in human bocaviruses. Partial VP1/VP2 sequence analysis of 23 FBoV and 25 CBoV strains demonstrated inter-host genetic diversity, with two potential genetic groups of FBoV and a novel CBoV group, CBoV-HK, distinct from the three groups, CBoV-A to -C, found in the USA. Although the pathogenicity of FBoV and CBoV remains to be determined, their presence in different host tissues suggested wide tissue tropism.

Large-scale species delimitation method for hyperdiverse groups

Accelerating the description of biodiversity is a major challenge as extinction rates increase. Integrative taxonomy combining molecular, morphological, ecological and geographical data is seen as the best route to reliably identify species. Classic molluscan taxonomic methodology proposes primary species hypotheses (PSHs) based on shell morphology. However, in hyperdiverse groups, such as the molluscan family Turridae, where most of the species remain unknown and for which homoplasy and plasticity of morphological characters is common, shell-based PSHs can be arduous. A four-pronged approach was employed to generate robust species hypotheses of a 1000 specimen South-West Pacific Turridae data set in which: (i) analysis of COI DNA Barcode gene is coupled with (ii) species delimitation tools GMYC (General Mixed Yule Coalescence Method) and ABGD (Automatic Barcode Gap Discovery) to propose PSHs that are then (iii) visualized using Klee diagrams and (iv) evaluated with additional evidence, such as nuclear gene rRNA 28S, morphological characters, geographical and bathymetrical distribution to determine conclusive secondary species hypotheses (SSHs). The integrative taxonomy approach applied identified 87 Turridae species, more than doubling the amount previously known in the Gemmula genus. In contrast to a predominantly shell-based morphological approach, which over the last 30 years proposed only 13 new species names for the Turridae genus Gemmula, the integrative approach described here identified 27 novel species hypotheses not linked to available species names in the literature. The formalized strategy applied here outlines an effective and reproducible protocol for large-scale species delimitation of hyperdiverse groups.

Modeling protein evolution with several amino acid replacement matrices depending on site rates

Most protein substitution models use a single amino acid replacement matrix summarizing the biochemical properties of amino acids. However, site evolution is highly heterogeneous and depends on many factors that influence the substitution patterns. In this paper, we investigate the use of different substitution matrices for different site evolutionary rates. Indeed, the variability of evolutionary rates corresponds to one of the most apparent heterogeneity factors among sites, and there is no reason to assume that the substitution patterns remain identical regardless of the evolutionary rate. We first introduce LG4M, which is composed of four matrices, each corresponding to one discrete gamma rate category (of four). These matrices differ in their amino acid equilibrium distributions and in their exchangeabilities, contrary to the standard gamma model where only the global rate differs from one category to another. Next, we present LG4X, which also uses four different matrices, but leaves aside the gamma distribution and follows a distribution-free scheme for the site rates. All these matrices are estimated from a very large alignment database, and our two models are tested using a large sample of independent alignments. Detailed analysis of resulting matrices and models shows the complexity of amino acid substitutions and the advantage of flexible models such as LG4M and LG4X. Both significantly outperform single-matrix models, providing gains of dozens to hundreds of log-likelihood units for most data sets. LG4X obtains substantial gains compared with LG4M, thanks to its distribution-free scheme for site rates. Since LG4M and LG4X display such advantages but require the same memory space and have comparable running times to standard models, we believe that LG4M and LG4X are relevant alternatives to single replacement matrices. Our models, data, and software are available from http://www.atgc-montpellier.fr/models/lg4x.

No evidence for infection of UK prostate cancer patients with XMRV, BK virus, Trichomonas vaginalis or human papilloma viruses

The prevalence of specific infections in UK prostate cancer patients was investigated. Serum from 84 patients and 62 controls was tested for neutralisation of xenotropic murine leukaemia virus-related virus (XMRV) Envelope. No reactivity was found in the patient samples. In addition, a further 100 prostate DNA samples were tested for XMRV, BK virus, Trichomonas vaginalis and human papilloma viruses by nucleic acid detection techniques. Despite demonstrating DNA integrity and assay sensitivity, we failed to detect the presence of any of these agents in DNA samples, bar one sample that was weakly positive for HPV16. Therefore we conclude that these infections are absent in this typical cohort of men with prostate cancer.

Extensive genetic diversity, unique population structure and evidence of genetic exchange in the sexually transmitted parasite Trichomonas vaginalis

Background

Trichomonas vaginalis is the causative agent of human trichomoniasis, the most common non-viral sexually transmitted infection world-wide. Despite its prevalence, little is known about the genetic diversity and population structure of this haploid parasite due to the lack of appropriate tools. The development of a panel of microsatellite makers and SNPs from mining the parasite's genome sequence has paved the way to a global analysis of the genetic structure of the pathogen and association with clinical phenotypes.

Methodology/Principal Findings

Here we utilize a panel of T. vaginalis-specific genetic markers to genotype 235 isolates from Mexico, Chile, India, Australia, Papua New Guinea, Italy, Africa and the United States, including 19 clinical isolates recently collected from 270 women attending New York City sexually transmitted disease clinics. Using population genetic analysis, we show that T. vaginalis is a genetically diverse parasite with a unique population structure consisting of two types present in equal proportions world-wide. Parasites belonging to the two types (type 1 and type 2) differ significantly in the rate at which they harbor the T. vaginalis virus, a dsRNA virus implicated in parasite pathogenesis, and in their sensitivity to the widely-used drug, metronidazole. We also uncover evidence of genetic exchange, indicating a sexual life-cycle of the parasite despite an absence of morphologically-distinct sexual stages.

Conclusions/Significance

Our study represents the first robust and comprehensive evaluation of global T. vaginalis genetic diversity and population structure. Our identification of a unique two-type structure, and the clinically relevant phenotypes associated with them, provides a new dimension for understanding T. vaginalis pathogenesis. In addition, our demonstration of the possibility of genetic exchange in the parasite has important implications for genetic research and control of the disease.

The human parasite Trichomonas vaginalis causes trichomoniasis, the world's most common non-viral sexually transmitted infection. Research on T. vaginalis genetic diversity has been limited by a lack of appropriate genotyping tools. To address this problem, we recently published a panel of T. vaginalis-specific genetic markers; here we use these markers to genotype isolates collected from ten regions around the globe. We detect high levels of genetic diversity, infer a two-type population structure, identify clinically relevant differences between the two types, and uncover evidence of genetic exchange in what was believed to be a clonal organism. Together, these results greatly improve our understanding of the population genetics of T. vaginalis and provide insights into the possibility of genetic exchange in the parasite, with implications for the epidemiology and control of the disease. By taking into account the existence of different types and their unique characteristics, we can improve understanding of the wide range of symptoms that patients manifest and better implement appropriate drug treatment. In addition, by recognizing the possibility of genetic exchange, we are more equipped to address the growing concern of drug resistance and the mechanisms by which it may spread within parasite populations.

Molecular evolution of a malaria resistance gene (DARC) in primates

Genes involved in host-pathogen interactions are often strongly affected by positive natural selection. The Duffy antigen, coded by the Duffy antigen receptor for chemokines (DARC) gene, serves as a receptor for Plasmodium vivax in humans and for Plasmodium knowlesi in some nonhuman primates. In the majority of sub-Saharan Africans, a nucleic acid variant in GATA-1 of the gene promoter is responsible for the nonexpression of the Duffy antigen on red blood cells and consequently resistance to invasion by P. vivax. The Duffy antigen also acts as a receptor for chemokines and is expressed in red blood cells and many other tissues of the body. Because of this dual role, we sequenced a ~3,000-bp region encompassing the entire DARC gene as well as part of its 5' and 3' flanking regions in a phylogenetic sample of primates and used statistical methods to evaluate the nature of selection pressures acting on the gene during its evolution. We analyzed both coding and regulatory regions of the DARC gene. The regulatory analysis showed accelerated rates of substitution at several sites near known motifs. Our tests of positive selection in the coding region using maximum likelihood by branch sites and maximum likelihood by codon sites did not yield statistically significant evidence for the action of positive selection. However, the maximum likelihood test in which the gene was subdivided into different structural regions showed that the known binding region for P. vivax/P. knowlesi is under very different selective pressures than the remainder of the gene. In fact, most of the gene appears to be under strong purifying selection, but this is not evident in the binding region. We suggest that the binding region is under the influence of two opposing selective pressures, positive selection possibly exerted by the parasite and purifying selection exerted by chemokines.

Evolution and architecture of the inner membrane complex in asexual and sexual stages of the malaria parasite

The inner membrane complex (IMC) is a unifying morphological feature of all alveolate organisms. It consists of flattened vesicles underlying the plasma membrane and is interconnected with the cytoskeleton. Depending on the ecological niche of the organisms, the function of the IMC ranges from a fundamental role as reinforcement system to more specialized roles in motility and cytokinesis. In this article, we present a comprehensive evolutionary analysis of IMC components, which exemplifies the adaptive nature of the IMCs' protein composition. Focusing on eight structurally distinct proteins in the most prominent "genus" of the Alveolata-the malaria parasite Plasmodium-we demonstrate that the level of conservation is reflected in phenotypic characteristics, accentuated in differential spatial-temporal patterns of these proteins in the motile stages of the parasite's life cycle. Colocalization studies with the centromere and the spindle apparatus reveal their discriminative biogenesis. We also reveal that the IMC is an essential structural compartment for the development of the sexual stages of Plasmodium, as it seems to drive the morphological changes of the parasite during the long and multistaged process of sexual differentiation. We further found a Plasmodium-specific IMC membrane matrix protein that highlights transversal structures in gametocytes, which could represent a genus-specific structural innovation required by Plasmodium. We conclude that the IMC has an additional role during sexual development supporting morphogenesis of the cell, which in addition to its functions in the asexual stages highlights the multifunctional nature of the IMC in the Plasmodium life cycle.

Identification and characterization of a novel calcium-activated apyrase from Cryptosporidium parasites and its potential role in pathogenesis

Herein, we report the biochemical and functional characterization of a novel Ca²⁺-activated nucleoside diphosphatase (apyrase), CApy, of the intracellular gut pathogen Cryptosporidium. The purified recombinant CApy protein displayed activity, substrate specificity and calcium dependency strikingly similar to the previously described human apyrase, SCAN-1 (soluble calcium-activated nucleotidase 1). CApy was found to be expressed in both Cryptosporidium parvum oocysts and sporozoites, and displayed a polar localization in the latter, suggesting a possible co-localization with the apical complex of the parasite. In vitro binding experiments revealed that CApy interacts with the host cell in a dose-dependent fashion, implying the presence of an interacting partner on the surface of the host cell. Antibodies directed against CApy block Cryptosporidium parvum sporozoite invasion of HCT-8 cells, suggesting that CApy may play an active role during the early stages of parasite invasion. Sequence analyses revealed that the capy gene shares a high degree of homology with apyrases identified in other organisms, including parasites, insects and humans. Phylogenetic analysis argues that the capy gene is most likely an ancestral feature that has been lost from most apicomplexan genomes except Cryptosporidium, Neospora and Toxoplasma.

The complete chloroplast and mitochondrial genome sequences of Boea hygrometrica: insights into the evolution of plant organellar genomes

The complete nucleotide sequences of the chloroplast (cp) and mitochondrial (mt) genomes of resurrection plant Boea hygrometrica (Bh, Gesneriaceae) have been determined with the lengths of 153,493 bp and 510,519 bp, respectively. The smaller chloroplast genome contains more genes (147) with a 72% coding sequence, and the larger mitochondrial genome have less genes (65) with a coding faction of 12%. Similar to other seed plants, the Bh cp genome has a typical quadripartite organization with a conserved gene in each region. The Bh mt genome has three recombinant sequence repeats of 222 bp, 843 bp, and 1474 bp in length, which divide the genome into a single master circle (MC) and four isomeric molecules. Compared to other angiosperms, one remarkable feature of the Bh mt genome is the frequent transfer of genetic material from the cp genome during recent Bh evolution. We also analyzed organellar genome evolution in general regarding genome features as well as compositional dynamics of sequence and gene structure/organization, providing clues for the understanding of the evolution of organellar genomes in plants. The cp-derived sequences including tRNAs found in angiosperm mt genomes support the conclusion that frequent gene transfer events may have begun early in the land plant lineage.

Spatially structured populations with a low level of cryptic diversity in European marine Gastrotricha

Species of the marine meiofauna such as Gastrotricha are known to lack dispersal stages and are thus assumed to have low dispersal ability and low levels of gene flow between populations. Yet, most species are widely distributed, and this creates a paradox. To shed light on this apparent paradox, we test (i) whether such wide distribution may be due to misidentification and lumping of cryptic species with restricted distributions and (ii) whether spatial structures exist for the phylogeography of gastrotrichs. As a model, we used the genus Turbanella in NW Europe. DNA taxonomy using a mitochondrial and a nuclear marker supports distinctness of four traditional species (Turbanella ambronensis, T. bocqueti, T. mustela and T. cornuta) and provides evidence for two cryptic species within T. hyalina. An effect of geography on the within-species genetic structure is indeed present, with the potential for understanding colonization processes and for performing phylogeographic inference from microscopic animals. On the other hand, the occurrence of widely distributed haplotypes indicates long-distance dispersal as well, despite the assumed low dispersal ability of gastrotrichs.

Interkingdom gene transfer may contribute to the evolution of phytopathogenicity in botrytis cinerea

The ascomycete Botrytis cinerea is a phytopathogenic fungus infecting and causing significant yield losses in a number of crops. The genome of B. cinerea has been fully sequenced while the importance of horizontal gene transfer (HGT) to extend the host range in plant pathogenic fungi has been recently appreciated. However, recent data confirm that the B. cinerea fungus shares conserved virulence factors with other fungal plant pathogens with narrow host range. Therefore, interkingdom HGT may contribute to the evolution of phytopathogenicity in B. cinerea. In this study, a stringent genome comparison pipeline was used to identify potential genes that have been obtained by B. cinerea but not by other fungi through interkingdom HGT. This search led to the identification of four genes: a UDP-glucosyltransferase (UGT), a lipoprotein and two alpha/beta hydrolase fold proteins. Phylogenetic analysis of the four genes suggests that B. cinerea acquired UGT from plants and the other 3 genes from bacteria. Based on the known gene functions and literature searching, a correlation between gene acquision and the evolution of pathogenicity in B. cinerea can be postulated.

Widespread impact of horizontal gene transfer on plant colonization of land

In complex multicellular eukaryotes such as animals and plants, horizontal gene transfer is commonly considered rare with very limited evolutionary significance. Here we show that horizontal gene transfer is a dynamic process occurring frequently in the early evolution of land plants. Our genome analyses of the moss Physcomitrella patens identified 57 families of nuclear genes that were acquired from prokaryotes, fungi or viruses. Many of these gene families were transferred to the ancestors of green or land plants. Available experimental evidence shows that these anciently acquired genes are involved in some essential or plant-specific activities such as xylem formation, plant defence, nitrogen recycling as well as the biosynthesis of starch, polyamines, hormones and glutathione. These findings suggest that horizontal gene transfer had a critical role in the transition of plants from aquatic to terrestrial environments. On the basis of these findings, we propose a model of horizontal gene transfer mechanism in nonvascular and seedless vascular plants.

Widespread impact of horizontal gene transfer on plant colonization of land

In complex multicellular eukaryotes such as animals and plants, horizontal gene transfer is commonly considered rare with very limited evolutionary significance. Here we show that horizontal gene transfer is a dynamic process occurring frequently in the early evolution of land plants. Our genome analyses of the moss Physcomitrella patens identified 57 families of nuclear genes that were acquired from prokaryotes, fungi or viruses. Many of these gene families were transferred to the ancestors of green or land plants. Available experimental evidence shows that these anciently acquired genes are involved in some essential or plant-specific activities such as xylem formation, plant defence, nitrogen recycling as well as the biosynthesis of starch, polyamines, hormones and glutathione. These findings suggest that horizontal gene transfer had a critical role in the transition of plants from aquatic to terrestrial environments. On the basis of these findings, we propose a model of horizontal gene transfer mechanism in nonvascular and seedless vascular plants.

Assessment of inactivating stop codon mutations in forty Saccharomyces cerevisiae strains: implications for [PSI] prion- mediated phenotypes

The yeast prion [PSI⁺] has been implicated in the generation of novel phenotypes by a mechanism involving a reduction in translation fidelity causing readthrough of naturally occurring stop codons. Some [PSI⁺] associated phenotypes may also be generated due to readthrough of inactivating stop codon mutations (ISCMs). Using next generation sequencing we have sequenced the genomes of two Saccharomyces cerevisiae strains that are commonly used for the study of the yeast [PSI⁺] prion. We have identified approximately 26,000 and 6,500 single nucleotide polymorphisms (SNPs) in strains 74-D694 and G600 respectively, compared to reference strain S288C. In addition to SNPs that produce non-synonymous amino acid changes we have also identified a number of SNPs that cause potential ISCMs in these strains, one of which we show is associated with a [PSI⁺]-dependent stress resistance phenotype in strain G600. We identified twenty-two potential ISCMs in strain 74-D694, present in genes involved in a variety of cellular processes including nitrogen metabolism, signal transduction and oxidative stress response. The presence of ISCMs in a subset of these genes provides possible explanations for previously identified [PSI⁺]-associated phenotypes in this strain. A comparison of ISCMs in strains G600 and 74-D694 with S. cerevisiae strains sequenced as part of the Saccharomyces Genome Resequencing Project (SGRP) shows much variation in the generation of strain-specific ISCMs and suggests this process is possible under complex genetic control. Additionally we have identified a major difference in the abilities of strains G600 and 74-D694 to grow at elevated temperatures. However, this difference appears unrelated to novel SNPs identified in strain 74-D694 present in proteins involved in the heat shock response, but may be attributed to other SNP differences in genes previously identified as playing a role in high temperature growth.

Molecular and phenotypic evidence of a new species of genus Esox (Esocidae, Esociformes, Actinopterygii): the southern pike, Esox flaviae

We address the taxonomic position of the southern European individuals of pike, performing a series of tests and comparisons from morphology, DNA taxonomy and population genetics parameters, in order to support the hypothesis that two species of pike, and not only one, exist in Europe. A strong relationship emerged between a northern genotype supported by COI, Cytb, AFLP and specific fragments, and a phenotype with round spot skin colour pattern and a large number of scales in the lateral line, clearly separated from a southern genotype with other skin colour pattern and a low number of scales in the lateral line. DNA taxonomy, based on a coalescent approach (GMYC) from phylogenetic reconstructions on COI and Cytb together with AFLP admixture analysis, supported the existence of two independently evolving entities. Such differences are not simply due to geographic distances, as northern European samples are more similar to Canadian and Chinese samples than the southern Europe ones. Thus, given that the differences between the two groups of European pike are significant at the phenotypic, genotypic and geographical levels, we propose the identification of two pike species: the already known northern pike (Esox lucius) and the southern pike (E. flaviae n.sp.). The correct identification of these two lineages as independent species should give rise to a ban on the introduction of northern pikes in southern Europe for recreational fishing, due to potential problems of hybridisation.

Comparative analysis of GT14/GT14-like gene family in Arabidopsis, Oryza, Populus, Sorghum and Vitis

Glycosyltransferase family14 (GT14) belongs to the glycosyltransferase (GT) superfamily that plays important roles in the biosynthesis of cell walls, the most abundant source of cellulosic biomass for bioethanol production. It has been hypothesized that DUF266 proteins are a new class of GTs related to GT14. In this study, we identified 62 GT14 and 106 DUF266 genes (named GT14-like herein) in Arabidopsis, Oryza, Populus, Sorghum and Vitis. Our phylogenetic analysis separated GT14 and GT14-like genes into two distinct clades, which were further divided into eight and five groups, respectively. Similarities in protein domain, 3D structure and gene expression were uncovered between the two phylogenetic clades, supporting the hypothesis that GT14 and GT14-like genes belong to one family. Therefore, we proposed a new family name, GT14/GT14-like family that combines both subfamilies. Variation in gene expression and protein subcellular localization within the GT14-like subfamily were greater than those within the GT14 subfamily. One-half of the Arabidopsis and Populus GT14/GT14-like genes were found to be preferentially expressed in stem/xylem, indicating that they are likely involved in cell wall biosynthesis. This study provided new insights into the evolution and functional diversification of the GT14/GT14-like family genes.

Adaptive evolution of the matrix extracellular phosphoglycoprotein in mammals

Background

Matrix extracellular phosphoglycoprotein (MEPE) belongs to a family of small integrin-binding ligand N-linked glycoproteins (SIBLINGs) that play a key role in skeleton development, particularly in mineralization, phosphate regulation and osteogenesis. MEPE associated disorders cause various physiological effects, such as loss of bone mass, tumors and disruption of renal function (hypophosphatemia). The study of this developmental gene from an evolutionary perspective could provide valuable insights on the adaptive diversification of morphological phenotypes in vertebrates.

Results

Here we studied the adaptive evolution of the MEPE gene in 26 Eutherian mammals and three birds. The comparative genomic analyses revealed a high degree of evolutionary conservation of some coding and non-coding regions of the MEPE gene across mammals indicating a possible regulatory or functional role likely related with mineralization and/or phosphate regulation. However, the majority of the coding region had a fast evolutionary rate, particularly within the largest exon (1467 bp). Rodentia and Scandentia had distinct substitution rates with an increased accumulation of both synonymous and non-synonymous mutations compared with other mammalian lineages. Characteristics of the gene (e.g. biochemical, evolutionary rate, and intronic conservation) differed greatly among lineages of the eight mammalian orders. We identified 20 sites with significant positive selection signatures (codon and protein level) outside the main regulatory motifs (dentonin and ASARM) suggestive of an adaptive role. Conversely, we find three sites under selection in the signal peptide and one in the ASARM motif that were supported by at least one selection model. The MEPE protein tends to accumulate amino acids promoting disorder and potential phosphorylation targets.

Conclusion

MEPE shows a high number of selection signatures, revealing the crucial role of positive selection in the evolution of this SIBLING member. The selection signatures were found mainly outside the functional motifs, reinforcing the idea that other regions outside the dentonin and the ASARM might be crucial for the function of the protein and future studies should be undertaken to understand its importance.

RNA polymerase beta subunit (rpoB) gene and the 16S-23S rRNA intergenic transcribed spacer region (ITS) as complementary molecular markers in addition to the 16S rRNA gene for phylogenetic analysis and identification of the species of the family Mycoplasmataceae

Conventional classification of the species in the family Mycoplasmataceae is mainly based on phenotypic criteria, which are complicated, can be difficult to measure, and have the potential to be hampered by phenotypic deviations among the isolates. The number of biochemical reactions suitable for phenotypic characterization of the Mycoplasmataceae is also very limited and therefore the strategy for the final identification of the Mycoplasmataceae species is based on comparative serological results. However, serological testing of the Mycoplasmataceae species requires a performance panel of hyperimmune sera which contains anti-serum to each known species of the family, a high level of technical expertise, and can only be properly performed by mycoplasma-reference laboratories. In addition, the existence of uncultivated and fastidious Mycoplasmataceae species/isolates in clinical materials significantly complicates, or even makes impossible, the application of conventional bacteriological tests. The analysis of available genetic markers is an additional approach for the primary identification and phylogenetic classification of cultivable species and uncultivable or fastidious organisms in standard microbiological laboratories. The partial nucleotide sequences of the RNA polymerase β-subunit gene (rpoB) and the 16S-23S rRNA intergenic transcribed spacer (ITS) were determined for all known type strains and the available non-type strains of the Mycoplasmataceae species. In addition to the available 16S rRNA gene data, the ITS and rpoB sequences were used to infer phylogenetic relationships among these species and to enable identification of the Mycoplasmataceae isolates to the species level. The comparison of the ITS and rpoB phylogenetic trees with the 16S rRNA reference phylogenetic tree revealed a similar clustering patterns for the Mycoplasmataceae species, with minor discrepancies for a few species that demonstrated higher divergence of their ITS and rpoB in comparison to their neighbor species. Overall, our results demonstrated that the ITS and rpoB gene could be useful complementary phylogenetic markers to infer phylogenetic relationships among the Mycoplasmataceae species and provide useful background information for the choice of appropriate metabolic and serological tests for the final classification of isolates. In summary, three-target sequence analysis, which includes the ITS, rpoB, and 16S rRNA genes, was demonstrated to be a reliable and useful taxonomic tool for the species differentiation within the family Mycoplasmataceae based on their phylogenetic relatedness and pairwise sequence similarities. We believe that this approach might also become a valuable tool for routine analysis and primary identification of new isolates in medical and veterinary microbiological laboratories.

Multiple geographic origins of commensalism and complex dispersal history of Black Rats

The Black Rat (Rattus rattus) spread out of Asia to become one of the world's worst agricultural and urban pests, and a reservoir or vector of numerous zoonotic diseases, including the devastating plague. Despite the global scale and inestimable cost of their impacts on both human livelihoods and natural ecosystems, little is known of the global genetic diversity of Black Rats, the timing and directions of their historical dispersals, and the risks associated with contemporary movements. We surveyed mitochondrial DNA of Black Rats collected across their global range as a first step towards obtaining an historical genetic perspective on this socioeconomically important group of rodents. We found a strong phylogeographic pattern with well-differentiated lineages of Black Rats native to South Asia, the Himalayan region, southern Indochina, and northern Indochina to East Asia, and a diversification that probably commenced in the early Middle Pleistocene. We also identified two other currently recognised species of Rattus as potential derivatives of a paraphyletic R. rattus. Three of the four phylogenetic lineage units within R. rattus show clear genetic signatures of major population expansion in prehistoric times, and the distribution of particular haplogroups mirrors archaeologically and historically documented patterns of human dispersal and trade. Commensalism clearly arose multiple times in R. rattus and in widely separated geographic regions, and this may account for apparent regionalism in their associated pathogens. Our findings represent an important step towards deeper understanding the complex and influential relationship that has developed between Black Rats and humans, and invite a thorough re-examination of host-pathogen associations among Black Rats.

Supermatrix data highlight the phylogenetic relationships of photosynthetic stramenopiles

Molecular data had consistently recovered monophyletic classes for the heterokont algae, however, the relationships among the classes had remained only partially resolved. Furthermore, earlier studies did not include representatives from all taxonomic classes. We used a five-gene (nuclear encoded SSU rRNA; plastid encoded rbcL, psaA, psbA, psbC) analysis with a subset of 89 taxa representing all 16 heterokont classes to infer a phylogenetic tree. There were three major clades. The Aurearenophyceae, Chrysomerophyceae, Phaeophyceae, Phaeothamniophyceae, Raphidophyceae, Schizocladiophyceae and Xanthophyceae formed the SI clade. The Chrysophyceae, Eustigmatophyceae, Pinguiophyceae, Synchromophyceae and Synurophyceae formed the SII clade. The Bacillariophyceae, Bolidophyceae, Dictyochophyceae and Pelagophyceae formed the SIII clade. These three clades were also found in a ten-gene analysis. The approximately unbiased test rejected alternative hypotheses that forced each class into either of the other two clades. Morphological and biochemical data were not available for all 89 taxa, however, existing data were consistent with the molecular phylogenetic tree, especially for the SIII clade.

Embryonic development and expression analysis of Distal-less in Caprella scaura (Crustacea, Amphipoda, Caprellidea)

The Caprellidea generally possess rudimentary abdomens and degenerated third and fourth pereopods. Previous molecular phylogenetic studies support the concept that their unique body plan is derived from a gammarid-like body plan from which the abdomen or third and fourth pereopods have been lost in the Caprellidea. To understand the developmental and genetic mechanisms for the morphological evolution of the Caprellidea, we observed the embryonic development of Caprella scaura. Although in the early embryonic stage limb buds appeared in all of the pereonites, we found that elongation of the limb buds did not occur in the third and fourth pereonites; instead, only oval projections (possibly primordial gills) were observed. We next examined the gene expression of Distal-less (Dll) by in situ hybridization and found that Dll was not expressed in the third and fourth pereonites. This suggests that the suppression of Dll expression is responsible for the reduction of Caprellidea pereopods.

Mitosomes in trophozoites and cysts of the reptilian parasite Entamoeba invadens

Heat shock protein genes led to the discovery of mitosomes in Entamoeba histolytica, but mitosomes have not been described for any other Entamoeba species, nor have they been identified in the cyst stage. Here, we show that the distantly related reptilian pathogen Entamoeba invadens contains mitosomes, in both trophozoites and cysts, suggesting all Entamoeba species contain these organelles.

Species interactions mediate phylogenetic community structure in a hyperdiverse lizard assemblage from arid Australia

Evolutionary history can exert a profound influence on ecological communities, but few generalities have emerged concerning the relationships among phylogeny, community membership, and niche evolution. We compared phylogenetic community structure and niche evolution in three lizard clades (Ctenotus skinks, agamids, and diplodactyline geckos) from arid Australia. We surveyed lizard communities at 32 sites in the northwestern Great Victoria Desert and generated complete species-level molecular phylogenies for regional representatives of the three clades. We document a striking pattern of phylogenetic evenness within local communities for all groups: pairwise correlations in species abundance across sites are negatively related to phylogenetic similarity. By modeling site suitability on the basis of species' habitat preferences, we demonstrate that phylogenetic evenness generally persists even after controlling for habitat filtering among species. This phylogenetic evenness is coupled with evolutionary lability of habitat-associated traits, to the extent that closely related species are more divergent in habitat use than distantly related species. In contrast, lizard diets are phylogenetically conserved, and pairwise dietary overlap between species is negatively related to phylogenetic distance in two of the three clades. Our results suggest that contemporary and historical species interactions have led to similar patterns of community structure across multiple clades in one of the world's most diverse lizard communities.

Analyses of the oligopeptide transporter gene family in poplar and grape

Background

Oligopeptide transporters (OPTs) are a group of membrane-localized proteins that have a broad range of substrate transport capabilities and that are thought to contribute to many biological processes. The OPT proteins belong to a small gene family in plants, which includes about 25 members in Arabidopsis and rice. However, no comprehensive study incorporating phylogeny, chromosomal location, gene structure, expression profiling, functional divergence and selective pressure analysis has been reported thus far for Populus and Vitis.

Results

In the present study, a comprehensive analysis of the OPT gene family in Populus (P. trichocarpa) and Vitis (V. vinifera) was performed. A total of 20 and 18 full-length OPT genes have been identified in Populus and Vitis, respectively. Phylogenetic analyses indicate that these OPT genes consist of two classes that can be further subdivided into 11 groups. Gene structures are considerably conserved among the groups. The distribution of OPT genes was found to be non-random across chromosomes. A high proportion of the genes are preferentially clustered, indicating that tandem duplications may have contributed significantly to the expansion of the OPT gene family. Expression patterns based on our analyses of microarray data suggest that many OPT genes may be important in stress response and functional development of plants. Further analyses of functional divergence and adaptive evolution show that, while purifying selection may have been the main force driving the evolution of the OPTs, some of critical sites responsible for the functional divergence may have been under positive selection.

Conclusions

Overall, the data obtained from our investigation contribute to a better understanding of the complexity of the Populus and Vitis OPT gene family and of the function and evolution of the OPT gene family in higher plants.