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

HIV behind bars: human immunodeficiency virus cluster analysis and drug resistance in a reference correctional unit from southern Brazil

People deprived of liberty in prisons are at higher risk of infection by the human immunodeficiency virus (HIV) due to their increased exposure through intravenous drug use, unprotected sexual activity, tattooing in prison and blood exposure in fights and rebellions. Yet, the contribution of intramural HIV transmission to the epidemic is scarcely known, especially in low- and middle-income settings. In this study, we surveyed 1,667 inmates incarcerated at Presídio Central de Porto Alegre, located in southern Brazil, for HIV infection and molecular characterization. The HIV seroprevalence was 6.6% (110/1,667). Further analyses were carried out on 40 HIV-seropositive inmates to assess HIV transmission clusters and drug resistance within the facility with the use of molecular and phylogenetic techniques. The molecular epidemiology of HIV-1 subtypes observed was similar to the one reported for the general population in southern Brazil, with the predominance of HIV-1 subtypes C, B, CRF31_BC and unique BC recombinants. In particular, the high rate (24%) of URF_BC found here may reflect multiple exposures of the population investigated to HIV infection. We failed to find HIV-infected inmates sharing transmission clusters with each other. Importantly, the analysis of HIV-1 pol genomic fragments evidenced high rates of HIV primary and secondary (acquired) drug resistance and an alarming proportion of virologic failure among patients under treatment, unveiling suboptimal access to antiretroviral therapy (ARV), low ARV adherence and dissemination of drug resistant HIV strains in primary infections. Our results call for immediate actions of public authority to implement preventive measures, serological screening and, for HIV-seropositive subjects, clinical and treatment follow-up in order to control HIV infection and limit the spread of drug resistance strains in Brazilian prisons.

Genetic characterization of EV71 isolates from 2004 to 2010 reveals predominance and persistent circulation of the newly proposed genotype D and recent emergence of a distinct lineage of subgenotype C2 in Hong Kong

Background

Enterovirus 71 (EV71) is a common etiological agent of hand, foot and mouth disease (HFMD) in children. EV71 epidemics have been reported in Hong Kong in recent years, and yet the genetic information of EV71 strains circulating in our locality is limited. The objective of this study was to investigate the genetic evolution of these EV71 isolates in Hong Kong over a 7-year period.

Methods

Twenty-two EV71 isolates from Hong Kong during 2004–2010 were included for phylogenetic analysis using partial VP2-VP3, 2C and 3D regions. Eight EV71 strains were selected for complete genome sequencing and recombination analysis.

Results

Among the 22 EV71 isolates, 20 belonged to subgenotype C4 and 2 belonged to subgenotype C2 based on the phylogenetic analysis of partial VP2-VP3, 2C and 3D gene regions. Phylogenetic, similarity plot and bootscan analyses using complete genome sequences of seven EV71 isolates of subgenotype C4 supported that the “double-recombinant” strains of subgenotype C4 persistently circulating in Hong Kong should belong to a newly proposed genotype D. Further analysis revealed two clusters, subgenotypes C4b and C4a (proposed genotypes D1a and D1b respectively), with “genotype D1b” strains being predominant in recent years in Hong Kong. A distinct lineage of EV71 subgenotype C2 has emerged in Hong Kong in 2008. The evolutionary rate of EV71 was 3.1 × 10^-3 nucleotide substitutions per site per year similar to that of other enterovirus, such as EV68, but was relatively lower than those of echovirus 30 and poliovirus. Molecular clock analysis using VP1 gene dated the time to the most recent common ancestor of all EV71 genotypes to 1900s, while the EV71 “double-recombinant” strains of “genotype D” were detected as early as 1998.

Conclusions

This study provides the molecular basis for proposing a new “genotype D” of EV71 and assigning a discrete lineage of subgenotype C2. EV71 strains of “genotype D” have been circulating in Hong Kong for over 7 years, with “genotype D1b” being predominant.

Recalibrating Equus evolution using the genome sequence of an early Middle Pleistocene horse

The rich fossil record of equids has made them a model for evolutionary processes. Here we present a 1.12-times coverage draft genome from a horse bone recovered from permafrost dated to approximately 560-780 thousand years before present (kyr BP). Our data represent the oldest full genome sequence determined so far by almost an order of magnitude. For comparison, we sequenced the genome of a Late Pleistocene horse (43 kyr BP), and modern genomes of five domestic horse breeds (Equus ferus caballus), a Przewalski's horse (E. f. przewalskii) and a donkey (E. asinus). Our analyses suggest that the Equus lineage giving rise to all contemporary horses, zebras and donkeys originated 4.0-4.5 million years before present (Myr BP), twice the conventionally accepted time to the most recent common ancestor of the genus Equus. We also find that horse population size fluctuated multiple times over the past 2 Myr, particularly during periods of severe climatic changes. We estimate that the Przewalski's and domestic horse populations diverged 38-72 kyr BP, and find no evidence of recent admixture between the domestic horse breeds and the Przewalski's horse investigated. This supports the contention that Przewalski's horses represent the last surviving wild horse population. We find similar levels of genetic variation among Przewalski's and domestic populations, indicating that the former are genetically viable and worthy of conservation efforts. We also find evidence for continuous selection on the immune system and olfaction throughout horse evolution. Finally, we identify 29 genomic regions among horse breeds that deviate from neutrality and show low levels of genetic variation compared to the Przewalski's horse. Such regions could correspond to loci selected early during domestication.

Molecular evolutionary analysis of the Alfin-like protein family in Arabidopsis lyrata, Arabidopsis thaliana, and Thellungiella halophila

In previous studies, the Alfin1 gene, a transcription factor, enhanced salt tolerance in alfalfa, primarily through altering gene expression levels in the root. Here, we examined the molecular evolution of the Alfin-like (AL) proteins in two Arabidopsis species (A. lyrata and A. thaliana) and a salt-tolerant close relative Thellungiella halophila. These AL-like proteins could be divided into four groups and the two known DUF3594 and PHD-finger domains had co-evolved within each group of genes, irrespective of species, due to gene duplication events in the common ancestor of all three species while gene loss was observed only in T. halophila. To detect whether natural selection acted in the evolution of AL genes, we calculated synonymous substitution ratios (dn/ds) and codon usage statistics, finding positive selection operated on four branches and significant differences in biased codon usage in the AL family between T. halophila and A. lyrata or A. thaliana. Distinctively, only the AL7 branch was under positive selection on the PHD-finger domain and the three members on the branch showed the smallest difference when codon bias was evaluated among the seven clusters. Functional analysis based on transgenic overexpression lines and T-DNA insertion mutants indicated that salt-stress-induced AtAL7 could play a negative role in salt tolerance of A. thaliana, suggesting that adaptive evolution occurred in the members of AL gene family.

Heterogeneous models place the root of the placental mammal phylogeny

Heterogeneity among life traits in mammals has resulted in considerable phylogenetic conflict, particularly concerning the position of the placental root. Layered upon this are gene- and lineage-specific variation in amino acid substitution rates and compositional biases. Life trait variations that may impact upon mutational rates are longevity, metabolic rate, body size, and germ line generation time. Over the past 12 years, three main conflicting hypotheses have emerged for the placement of the placental root. These hypotheses place the Atlantogenata (common ancestor of Xenarthra plus Afrotheria), the Afrotheria, or the Xenarthra as the sister group to all other placental mammals. Model adequacy is critical for accurate tree reconstruction and by failing to account for these compositional and character exchange heterogeneities across the tree and data set, previous studies have not provided a strongly supported hypothesis for the placental root. For the first time, models that accommodate both tree and data set heterogeneity have been applied to mammal data. Here, we show the impact of accurate model assignment and the importance of data sets in accommodating model parameters while maintaining the power to reject competing hypotheses. Through these sophisticated methods, we demonstrate the importance of model adequacy, data set power and provide strong support for the Atlantogenata over other competing hypotheses for the position of the placental root.

Evaluating phylogenetic informativeness as a predictor of phylogenetic signal for metazoan, fungal, and mammalian phylogenomic data sets

Phylogenetic research is often stymied by selection of a marker that leads to poor phylogenetic resolution despite considerable cost and effort. Profiles of phylogenetic informativeness provide a quantitative measure for prioritizing gene sampling to resolve branching order in a particular epoch. To evaluate the utility of these profiles, we analyzed phylogenomic data sets from metazoans, fungi, and mammals, thus encompassing diverse time scales and taxonomic groups. We also evaluated the utility of profiles created based on simulated data sets. We found that genes selected via their informativeness dramatically outperformed haphazard sampling of markers. Furthermore, our analyses demonstrate that the original phylogenetic informativeness method can be extended to trees with more than four taxa. Thus, although the method currently predicts phylogenetic signal without specifically accounting for the misleading effects of stochastic noise, it is robust to the effects of homoplasy. The phylogenetic informativeness rankings obtained will allow other researchers to select advantageous genes for future studies within these clades, maximizing return on effort and investment. Genes identified might also yield efficient experimental designs for phylogenetic inference for many sister clades and outgroup taxa that are closely related to the diverse groups of organisms analyzed.

The origins of the enigmatic Falkland Islands wolf

The origins of the extinct Falkland Islands wolf (FIW), Dusicyon australis, have remained a mystery since it was first recorded by Europeans in the seventeenth century. It is the only terrestrial mammal on the Falkland Islands (also known as the Malvinas Islands), which lie ~460 km from Argentina, leading to suggestions of either human-mediated transport or overwater dispersal. Previous studies used ancient DNA from museum specimens to suggest that the FIW diverged from its closest living relative, the South American maned wolf (Chrysocyon brachyurus) around 7 Ma, and colonized the islands ~330 ka by unknown means. Here we retrieve ancient DNA from subfossils of an extinct mainland relative, Dusicyon avus, and reveal the FIW lineage became isolated only 16 ka (8-31 ka), during the last glacial phase. Submarine terraces, formed on the Argentine coastal shelf by low sea-stands during this period, suggest that the FIW colonized via a narrow, shallow marine strait, potentially while it was frozen over.

Positive selection along the evolution of primate mitogenomes

The mitochondrial genomes of four neotropical primates, Aotus infulatus, Chiropotes israelita, Callimico goeldii and Callicebus lugens were sequenced and annotated. Phylogenetic reconstructions with mitochondrial genes of other 66 primates showed a similar arrangement to a topology based on nuclear genes. Screening for positive selection identified 15 codons in 7 genes along 9 independent lineages, three with two or more genes and five in internal nodes, ruling out false positive estimates. Mitochondrial genes of the electron transport chain (ETC.) complexes evolved with high substitution rates. A study of nuclear ETC. genes might elucidate whether they co-evolved with their mitochondrial counterparts.

Mitogenomic sequences and evidence from unique gene rearrangements corroborate evolutionary relationships of myctophiformes (Neoteleostei)

Background

A skewed assemblage of two epi-, meso- and bathypelagic fish families makes up the order Myctophiformes – the blackchins Neoscopelidae and the lanternfishes Myctophidae. The six rare neoscopelids show few morphological specializations whereas the divergent myctophids have evolved into about 250 species, of which many show massive abundances and wide distributions. In fact, Myctophidae is by far the most abundant fish family in the world, with plausible estimates of more than half of the oceans combined fish biomass. Myctophids possess a unique communication system of species-specific photophore patterns and traditional intrafamilial classification has been established to reflect arrangements of photophores. Myctophids present the most diverse array of larval body forms found in fishes although this attribute has both corroborated and confounded phylogenetic hypotheses based on adult morphology. No molecular phylogeny is available for Myctophiformes, despite their importance within all ocean trophic cycles, open-ocean speciation and as an important part of neoteleost divergence. This study attempts to resolve major myctophiform phylogenies from both mitogenomic sequences and corroborating evidence in the form of unique mitochondrial gene order rearrangements.

Results

Mitogenomic evidence from DNA sequences and unique gene orders are highly congruent concerning phylogenetic resolution on several myctophiform classification levels, corroborating evidence from osteology, larval ontogeny and photophore patterns, although the lack of larval morphological characters within the subfamily Lampanyctinae stands out. Neoscopelidae is resolved as the sister family to myctophids with Solivomer arenidens positioned as a sister taxon to the remaining neoscopelids. The enigmatic Notolychnus valdiviae is placed as a sister taxon to all other myctophids and exhibits an unusual second copy of the tRNA-Met gene – a gene order rearrangement reminiscent of that found in the tribe Diaphini although our analyses show it to be independently derived. Most tribes are resolved in accordance with adult morphology although Gonichthyini is found within a subclade of the tribe Myctophini consisting of ctenoid scaled species. Mitogenomic sequence data from this study recognize 10 reciprocally monophyletic lineages within Myctophidae, with five of these clades delimited from additional rearranged gene orders or intergenic non-coding sequences.

Conclusions

Mitogenomic results from DNA sequences and unique gene orders corroborate morphology in phylogeny reconstruction and provide a likely scenario for the phylogenetic history of Myctophiformes. The extent of gene order rearrangements found within the mitochondrial genomes of myctophids is unique for phylogenetic purposes.

Genomic and proteomic dissection of the ubiquitous plant pathogen, Armillaria mellea: toward a new infection model system

Armillaria mellea is a major plant pathogen. Yet, no large-scale “-omics” data are available to enable new studies, and limited experimental models are available to investigate basidiomycete pathogenicity. Here we reveal that the A. mellea genome comprises 58.35 Mb, contains 14473 gene models, of average length 1575 bp (4.72 introns/gene). Tandem mass spectrometry identified 921 mycelial (n = 629 unique) and secreted (n = 183 unique) proteins. Almost 100 mycelial proteins were either species-specific or previously unidentified at the protein level. A number of proteins (n = 111) was detected in both mycelia and culture supernatant extracts. Signal sequence occurrence was 4-fold greater for secreted (50.2%) compared to mycelial (12%) proteins. Analyses revealed a rich reservoir of carbohydrate degrading enzymes, laccases, and lignin peroxidases in the A. mellea proteome, reminiscent of both basidiomycete and ascomycete glycodegradative arsenals. We discovered that A. mellea exhibits a specific killing effect against Candida albicans during coculture. Proteomic investigation of this interaction revealed the unique expression of defensive and potentially offensive A. mellea proteins (n = 30). Overall, our data reveal new insights into the origin of basidiomycete virulence and we present a new model system for further studies aimed at deciphering fungal pathogenic mechanisms.

New insights into the evolution of Entomopoxvirinae from the complete genome sequences of four entomopoxviruses infecting Adoxophyes honmai, Choristoneura biennis, Choristoneura rosaceana, and Mythimna separata

Poxviruses are nucleocytoplasmic large DNA viruses encompassing two subfamilies, the Chordopoxvirinae and the Entomopoxvirinae, infecting vertebrates and insects, respectively. While chordopoxvirus genomics have been widely studied, only two entomopoxvirus (EPV) genomes have been entirely sequenced. We report the genome sequences of four EPVs of the Betaentomopoxvirus genus infecting the Lepidoptera: Adoxophyes honmai EPV (AHEV), Choristoneura biennis EPV (CBEV), Choristoneura rosaceana EPV (CREV), and Mythimna separata EPV (MySEV). The genomes are 80% AT rich, are 228 to 307 kbp long, and contain 247 to 334 open reading frames (ORFs). Most genes are homologous to those of Amsacta moorei entomopoxvirus and encode several protein families repeated in tandem in terminal regions. Some genomes also encode proteins of unknown functions with similarity to those of other insect viruses. Comparative genomic analyses highlight a high colinearity among the lepidopteran EPV genomes and little gene order conservation with other poxvirus genomes. As with previously sequenced EPVs, the genomes include a relatively conserved central region flanked by inverted terminal repeats. Protein clustering identified 104 core EPV genes. Among betaentomopoxviruses, 148 core genes were found in relatively high synteny, pointing to low genomic diversity. Whole-genome and spheroidin gene phylogenetic analyses showed that the lepidopteran EPVs group closely in a monophyletic lineage, corroborating their affiliation with the Betaentomopoxvirus genus as well as a clear division of the EPVs according to the orders of insect hosts (Lepidoptera, Coleoptera, and Orthoptera). This suggests an ancient coevolution of EPVs with their insect hosts and the need to revise the current EPV taxonomy to separate orthopteran EPVs from the lepidopteran-specific betaentomopoxviruses so as to form a new genus.

The role of historical barriers in the diversification processes in open vegetation formations during the Miocene/Pliocene using an ancient rodent lineage as a model

The Neotropics harbors a high diversity of species and several hypotheses have been proposed to account for this pattern. However, while species of forested domains are frequently studied, less is known of species from open vegetation formations occupying, altogether, a larger area than the Amazon Forest. Here we evaluate the role of historical barriers and the riverine hypothesis in the speciation patterns of small mammals by analyzing an ancient rodent lineage (Thrichomys, Hystricomorpha). Phylogenetic and biogeographic analyses were carried out with mitochondrial and nuclear DNA markers to analyze the evolutionary relationships between Thrichomys lineages occurring in dry domains along both banks of the Rio São Francisco. This river is one of the longest of South America whose course and water flow have been modified by inland tectonic activities and climate changes. Molecular data showed a higher number of lineages than previously described. The T. inermis species complex with 2n = 26, FN = 48 was observed in both banks of the river showing a paraphyletic arrangement, suggesting that river crossing had occurred, from east to west. A similar pattern was also observed for the T. apereoides complex. Thrichomys speciation occurred in Late Miocene when the river followed a different course. The current geographic distribution of Thrichomys species and their phylogenetic relationships suggested the existence of frequent past connections between both banks in the middle section of the Rio São Francisco. The extensive palaeodune region found in this area has been identified as a centre of endemism of several vertebrate species and is likely to be a center of Thrichomys diversification.

Inferring the evolutionary histories of divergences in Hylobates and Nomascus gibbons through multilocus sequence data

BACKGROUND

Gibbons (Hylobatidae) are the most diverse group of living apes. They exist as geographically-contiguous species which diverged more rapidly than did their close relatives, the great apes (Hominidae). Of the four extant gibbon genera, the evolutionary histories of two polyspecific genera, Hylobates and Nomascus, have been the particular focus of research but the DNA sequence data used was largely derived from the maternally inherited mitochondrial DNA (mtDNA) locus.

RESULTS

To investigate the evolutionary relationships and divergence processes of gibbon species, particularly those of the Hylobates genus, we produced and analyzed a total of 11.5 kb DNA of sequence at 14 biparentally inherited autosomal loci. We find that on average gibbon genera have a high average sequence diversity but a lower degree of genetic differentiation as compared to great ape genera. Our multilocus species tree features H. pileatus in a basal position and a grouping of the four Sundaic island species (H. agilis, H. klossii, H. moloch and H. muelleri). We conducted pairwise comparisons based on an isolation-with-migration (IM) model and detect signals of asymmetric gene flow between H. lar and H. moloch, between H. agilis and H. muelleri, and between N. leucogenys and N. siki.

CONCLUSIONS

Our multilocus analyses provide inferences of gibbon evolutionary histories complementary to those based on single gene data. The results of IM analyses suggest that the divergence processes of gibbons may be accompanied by gene flow. Future studies using analyses of multi-population model with samples of known provenance for Hylobates and Nomascus species would expand the understanding of histories of gene flow during divergences for these two gibbon genera.

Genetic diversity analysis of the Uruguayan Creole cattle breed using microsatellites and mtDNA markers

The Uruguayan Creole cattle population (N = 600) is located in a native habitat in south-east Uruguay. We analyzed its genetic diversity and compared it to other populations of American Creole cattle. A random sample of 64 animals was genotyped for a set of 17 microsatellite loci, and the D-loop hyper-variable region of mtDNA was sequenced for 28 calves of the same generation. We identified an average of 5.59 alleles per locus, with expected heterozygosities between 0.466 and 0.850 and an expected mean heterozygosity of 0.664. The polymorphic information content ranged from 0.360 to 0.820, and the global FIS index was 0.037. The D-loop analysis revealed three haplotypes (UY1, UY2 and UY3), belonging to the European matriline group, with a haplotype diversity of 0.532. The history of the population, changes in the effective population size, bottlenecks, and genetic drift are possible causes of the genetic variability patterns that we detected.

Genomics of Loa loa, a Wolbachia-free filarial parasite of humans

Loa loa, the African eyeworm, is a major filarial pathogen of humans. Unlike most filariae, Loa loa does not contain the obligate intracellular Wolbachia endosymbiont. We describe the 91.4 Mb genome of Loa loa, and the genome of the related filarial parasite Wuchereria bancrofti, and predict 14,907 Loa loa genes based on microfilarial RNA sequencing. By comparing these genomes to that of another filarial parasite, Brugia malayi, and to several other nematode genomes, we demonstrate synteny among filariae but not with non-parasitic nematodes. The Loa loa genome encodes many immunologically relevant genes, as well as protein kinases targeted by drugs currently approved for humans. Despite lacking Wolbachia, Loa loa shows no new metabolic synthesis or transport capabilities compared to other filariae. These results suggest that the role played by Wolbachia in filarial biology is more subtle than previously thought and reveal marked differences between parasitic and non-parasitic nematodes.

Horizontal transfer of a subtilisin gene from plants into an ancestor of the plant pathogenic fungal genus Colletotrichum

The genus Colletotrichum contains a large number of phytopathogenic fungi that produce enormous economic losses around the world. The effect of horizontal gene transfer (HGT) has not been studied yet in these organisms. Inter-Kingdom HGT into fungal genomes has been reported in the past but knowledge about the HGT between plants and fungi is particularly limited. We describe a gene in the genome of several species of the genus Colletotrichum with a strong resemblance to subtilisins typically found in plant genomes. Subtilisins are an important group of serine proteases, widely distributed in all of the kingdoms of life. Our hypothesis is that the gene was acquired by Colletotrichum spp. through (HGT) from plants to a Colletotrichum ancestor. We provide evidence to support this hypothesis in the form of phylogenetic analyses as well as a characterization of the similarity of the subtilisin at the primary, secondary and tertiary structural levels. The remarkable level of structural conservation of Colletotrichum plant-like subtilisin (CPLS) with plant subtilisins and the differences with the rest of Colletotrichum subtilisins suggests the possibility of molecular mimicry. Our phylogenetic analysis indicates that the HGT event would have occurred approximately 150–155 million years ago, after the divergence of the Colletotrichum lineage from other fungi. Gene expression analysis shows that the gene is modulated during the infection of maize by C. graminicola suggesting that it has a role in plant disease. Furthermore, the upregulation of the CPLS coincides with the downregulation of several plant genes encoding subtilisins. Based on the known roles of subtilisins in plant pathogenic fungi and the gene expression pattern that we observed, we postulate that the CPLSs have an important role in plant infection.

Nosema ceranae escapes fumagillin control in honey bees

Fumagillin is the only antibiotic approved for control of nosema disease in honey bees and has been extensively used in United States apiculture for more than 50 years for control of Nosema apis. It is toxic to mammals and must be applied seasonally and with caution to avoid residues in honey. Fumagillin degrades or is diluted in hives over the foraging season, exposing bees and the microsporidia to declining concentrations of the drug. We showed that spore production by Nosema ceranae, an emerging microsporidian pathogen in honey bees, increased in response to declining fumagillin concentrations, up to 100% higher than that of infected bees that have not been exposed to fumagillin. N. apis spore production was also higher, although not significantly so. Fumagillin inhibits the enzyme methionine aminopeptidase2 (MetAP2) in eukaryotic cells and interferes with protein modifications necessary for normal cell function. We sequenced the MetAP2 gene for apid Nosema species and determined that, although susceptibility to fumagillin differs among species, there are no apparent differences in fumagillin binding sites. Protein assays of uninfected bees showed that fumagillin altered structural and metabolic proteins in honey bee midgut tissues at concentrations that do not suppress microsporidia reproduction. The microsporidia, particularly N. ceranae, are apparently released from the suppressive effects of fumagillin at concentrations that continue to impact honey bee physiology. The current application protocol for fumagillin may exacerbate N. ceranae infection rather than suppress it.

Evolution of yolk protein genes in the Echinodermata

Vitellogenin genes (vtg) encode large lipid transfer proteins (LLTPs) that are typically female-specific, functioning as precursors to major yolk proteins (MYPs). Within the phylum Echinodermata, however, the MYP of the Echinozoa (Echinoidea + Holothuroidea) is expressed by an unrelated transferrin-like gene that has a reproductive function in both sexes. We investigated egg proteins in the Asterozoa (Asteroidea + Ophiuroidea), a sister clade to the Echinozoa, showing that eggs of the asteroid Parvulastra exigua contain a vitellogenin protein (Vtg). vtg is expressed by P. exigua, a species with large eggs and nonfeeding larvae, and by the related asterinid Patiriella regularis which has small eggs and feeding larvae. In the Asteroidea, therefore, the reproductive function of vtg is conserved despite significant life history evolution. Like the echinozoan MYP gene, asteroid vtg is expressed in both sexes and may play a role in the development of both ovaries and testes. Phylogenetic analysis indicated that a putative Vtg from the sea urchin genome, a likely pseudogene, does not clade with asteroid Vtg. We propose the following sequence as a potential pathway for the evolution of YP genes in the Echinodermata: (1) the ancestral echinoderm produced YPs derived from Vtg, (2) bisexual vtg expression subsequently evolved in the echinoderm lineage, (3) the reproductive function of vtg was assumed by a transferrin-like gene in the ancestral echinozoan, and (4) redundant echinozoan vtg was released from stabilizing selection.

Three independent determinants of protein evolutionary rate

One of the most widely accepted ideas related to the evolutionary rates of proteins is that functionally important residues or regions evolve slower than other regions, a reasonable outcome of which should be a slower evolutionary rate of the proteins with a higher density of functionally important sites. Oddly, the role of functional importance, mainly measured by essentiality, in determining evolutionary rate has been challenged in recent studies. Several variables other than protein essentiality, such as expression level, gene compactness, protein-protein interactions, etc., have been suggested to affect protein evolutionary rate. In the present review, we try to refine the concept of functional importance of a gene, and consider three factors-functional importance, expression level, and gene compactness, as independent determinants of evolutionary rate of a protein, based not only on their known correlation with evolutionary rate but also on a reasonable mechanistic model. We suggest a framework based on these mechanistic models to correctly interpret the correlations between evolutionary rates and the various variables as well as the interrelationships among the variables.

Has living on islands been so simple? Insights from the insular endemic frog Discoglossus montalentii

Island populations have been extensively used as model systems in ecology, biogeography, conservation and evolutionary biology, owing to the several simplifying assumptions that they allow. Nevertheless, recent findings from intra-island phylogeographic studies are casting doubts on the generality of some of these underlying assumptions. We investigated the phylogeography, historical demography, and population genetic structure of the Corsican endemic frog, Discoglossus montalentii. In contrast with expectations based on its insular, restricted and continuous distribution, we found evidence of 3 phylogroups, whose rather ancient divergence (Early-Middle Pleistocene) was likely primed by climatic changes that occurred during the ‘middle Pleistocene revolution’. Furthermore, their differentiation explained most (68%) of the overall genetic diversity that was observed. These results and the growing evidence from intra-island phylogeographies, suggest that island populations frequently may not conform to some long-standing assumptions, including long-term stability, range-wide panmixia and the correlation of effective population size to the island size. As a consequence, both for theoretical and for applied purposes, the extensive use of these assumptions in the study of island populations warrants a careful re-examination.

Functional divergence of the glutathione S-transferase supergene family in Physcomitrella patens reveals complex patterns of large gene family evolution in land plants

Plant glutathione S-transferases (GSTs) are multifunctional proteins encoded by a large gene family that play major roles in the detoxification of xenobiotics and oxidative stress metabolism. To date, studies on the GST gene family have focused mainly on vascular plants (particularly agricultural plants). In contrast, little information is available on the molecular characteristics of this large gene family in nonvascular plants. In addition, the evolutionary patterns of this family in land plants remain unclear. In this study, we identified 37 GST genes from the whole genome of the moss Physcomitrella patens, a nonvascular representative of early land plants. The 37 P. patens GSTs were divided into 10 classes, including two new classes (hemerythrin and iota). However, no tau GSTs were identified, which represent the largest class among vascular plants. P. patens GST gene family members showed extensive functional divergence in their gene structures, gene expression responses to abiotic stressors, enzymatic characteristics, and the subcellular locations of the encoded proteins. A joint phylogenetic analysis of GSTs from P. patens and other higher vascular plants showed that different class GSTs had distinct duplication patterns during the evolution of land plants. By examining multiple characteristics, this study revealed complex patterns of evolutionary divergence among the GST gene family in land plants.

Tetraflagellochloris mauritanica gen. et sp. nov. (Chlorophyceae), a New Flagellated Alga from the Mauritanian Desert: Morphology, Ultrastructure, and Phylogenetic Framing

Morphological, ultrastructural, and molecular-sequence data were used to assess the phylogenetic position of a tetraflagellate green alga isolated from soil samples of a saline dry basin near F'derick, Mauritania. This alga can grow as individual cells or form non-coenobial colonies of up to 12 individuals. It has a parietal chloroplast with an embedded pyrenoid covered by a starch sheath and traversed by single parallel thylakoids, and an eyespot located in a parietal position opposite to the flagellar insertion. Lipid vacuoles are present in the cytoplasm. Microspectroscopy indicated the presence of chlorophylls a and b, with lutein as the major carotenoid in the chloroplast, while the eyespot spectrum has a shape typical of green-algal eyespots. The cell has four flagella, two of them long and two considerably shorter. Sequence data from the 18S rRNA gene and ITS2 were obtained and compared with published sequences for green algae. Results from morphological and ultrastructural examinations and sequence analysis support the placement of this alga in the Chlorophyceae, as Tetraflagellochloris mauritanica L. Barsanti et A. Barsanti, gen. et sp. nov.

Discovery of a divergent HPIV4 from respiratory secretions using second and third generation metagenomic sequencing

Molecular detection of viruses has been aided by high-throughput sequencing, permitting the genomic characterization of emerging strains. In this study, we comprehensively screened 500 respiratory secretions from children with upper and/or lower respiratory tract infections for viral pathogens. The viruses detected are described, including a divergent human parainfluenza virus type 4 from GS FLX pyrosequencing of 92 specimens. Complete full-genome characterization of the virus followed, using Single Molecule, Real-Time (SMRT) sequencing. Subsequent "primer walking" combined with Sanger sequencing validated the RS platform's utility in viral sequencing from complex clinical samples. Comparative genomics reveals the divergent strain clusters with the only completely sequenced HPIV4a subtype. However, it also exhibits various structural features present in one of the HPIV4b reference strains, opening questions regarding their lifecycle and evolutionary relationships among these viruses. Clinical data from patients infected with the strain, as well as viral prevalence estimates using real-time PCR, is also described.

New approaches for unravelling reassortment pathways

BACKGROUND

Every year the human population encounters epidemic outbreaks of influenza, and history reveals recurring pandemics that have had devastating consequences. The current work focuses on the development of a robust algorithm for detecting influenza strains that have a composite genomic architecture. These influenza subtypes can be generated through a reassortment process, whereby a virus can inherit gene segments from two different types of influenza particles during replication. Reassortant strains are often not immediately recognised by the adaptive immune system of the hosts and hence may be the source of pandemic outbreaks. Owing to their importance in public health and their infectious ability, it is essential to identify reassortant influenza strains in order to understand the evolution of this virus and describe reassortment pathways that may be biased towards particular viral segments. Phylogenetic methods have been used traditionally to identify reassortant viruses. In many studies up to now, the assumption has been that if two phylogenetic trees differ, it is because reassortment has caused them to be different. While phylogenetic incongruence may be caused by real differences in evolutionary history, it can also be the result of phylogenetic error. Therefore, we wish to develop a method for distinguishing between topological inconsistency that is due to confounding effects and topological inconsistency that is due to reassortment.

RESULTS

The current work describes the implementation of two approaches for robustly identifying reassortment events. The algorithms rest on the idea of significance of difference between phylogenetic trees or phylogenetic tree sets, and subtree pruning and regrafting operations, which mimic the effect of reassortment on tree topologies. The first method is based on a maximum likelihood (ML) framework (MLreassort) and the second implements a Bayesian approach (Breassort) for reassortment detection. We focus on reassortment events that are found by both methods. We test both methods on a simulated dataset and on a small collection of real viral data isolated in Hong Kong in 1999.

CONCLUSIONS

The nature of segmented viral genomes present many challenges with respect to disease. The algorithms developed here can effectively identify reassortment events in small viral datasets and can be applied not only to influenza but also to other segmented viruses. Owing to computational demands of comparing tree topologies, further development in this area is necessary to allow their application to larger datasets.

Phylogenetically and spatially close marine sponges harbour divergent bacterial communities

Recent studies have unravelled the diversity of sponge-associated bacteria that may play essential roles in sponge health and metabolism. Nevertheless, our understanding of this microbiota remains limited to a few host species found in restricted geographical localities, and the extent to which the sponge host determines the composition of its own microbiome remains a matter of debate. We address bacterial abundance and diversity of two temperate marine sponges belonging to the Irciniidae family--Sarcotragus spinosulus and Ircinia variabilis--in the Northeast Atlantic. Epifluorescence microscopy revealed that S. spinosulus hosted significantly more prokaryotic cells than I. variabilis and that prokaryotic abundance in both species was about 4 orders of magnitude higher than in seawater. Polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE) profiles of S. spinosulus and I. variabilis differed markedly from each other--with higher number of ribotypes observed in S. spinosulus--and from those of seawater. Four PCR-DGGE bands, two specific to S. spinosulus, one specific to I. variabilis, and one present in both sponge species, affiliated with an uncultured sponge-specific phylogenetic cluster in the order Acidimicrobiales (Actinobacteria). Two PCR-DGGE bands present exclusively in S. spinosulus fingerprints affiliated with one sponge-specific phylogenetic cluster in the phylum Chloroflexi and with sponge-derived sequences in the order Chromatiales (Gammaproteobacteria), respectively. One Alphaproteobacteria band specific to S. spinosulus was placed in an uncultured sponge-specific phylogenetic cluster with a close relationship to the genus Rhodovulum. Our results confirm the hypothesized host-specific composition of bacterial communities between phylogenetically and spatially close sponge species in the Irciniidae family, with S. spinosulus displaying higher bacterial community diversity and distinctiveness than I. variabilis. These findings suggest a pivotal host-driven effect on the shape of the marine sponge microbiome, bearing implications to our current understanding of the distribution of microbial genetic resources in the marine realm.

Genome-scale comparative analysis of gene fusions, gene fissions, and the fungal tree of life

During the course of evolution genes undergo both fusion and fission by which ORFs are joined or separated. These processes can amend gene function and represent an important factor in the evolution of protein interaction networks. Gene fusions have been suggested to be useful characters for identifying evolutionary relationships because they constitute synapomorphies or cladistic characters. To investigate the fidelity of gene-fusion characters, we developed an approach for identifying differentially distributed gene fusions among whole-genome datasets: fdfBLAST. Applying this tool to the Fungi, we identified 63 gene fusions present in two or more genomes. Using a combination of phylogenetic and comparative genomic analyses, we then investigated the evolution of these genes across 115 fungal genomes, testing each gene fusion for evidence of homoplasy, including gene fission, convergence, and horizontal gene transfer. These analyses demonstrated 110 gene-fission events. We then identified a minimum of three mechanisms that drive gene fission: separation, degeneration, and duplication. These data suggest that gene fission plays an important and hitherto underestimated role in gene evolution. Gene fusions therefore are highly labile characters, and their use for polarizing evolutionary relationships, without reference to gene and species phylogenies, is limited. Accounting for these considerable sources of homoplasy, we identified fusion characters that provide support for multiple nodes in the phylogeny of the Fungi, including relationships within the deeply derived flagellum-forming fungi (i.e., the chytrids).

Heterochronic activation of VEGF signaling and the evolution of the skeleton in echinoderm pluteus larvae

The evolution of the echinoderm larval skeleton was examined from the aspect of interactions between skeletogenic mesenchyme cells and surrounding epithelium. We focused on vascular endothelial growth factor (VEGF) signaling, which was reported to be essential for skeletogenesis in sea urchin larvae. Here, we examined the expression patterns of vegf and vegfr in starfish and brittle stars. During starfish embryogenesis, no expression of either vegfr or vegf was detected, which contrast with previous reports on the expression of starfish homologs of sea urchin skeletogenic genes, including Ets, Tbr, and Dri. In later stages, when adult skeletogenesis commenced, vegfr and vegf expression were upregulated in skeletogenic cells and in the adjacent epidermis, respectively. These expression patterns suggest that heterochronic activation of VEGF signaling is one of the key molecular evolutionary steps in the evolution of the larval skeleton. The absence of vegf or vegfr expression during early embryogenesis in starfish suggests that the evolution of the larval skeleton requires distinct evolutionary changes, both in mesoderm cells (activation of vegfr expression) and in epidermal cells (activation of vegf expression). In brittle stars, which have well-organized skeletons like the sea urchin, vegfr and vegf were expressed in the skeletogenic mesenchyme and the overlying epidermis, respectively, in the same manner as in sea urchins. Therefore, the distinct activation of vegfr and vegf may have occurred in two lineages, sea urchins and brittle stars.

Broad distribution of TPI-GAPDH fusion proteins among eukaryotes: evidence for glycolytic reactions in the mitochondrion?

Glycolysis is a central metabolic pathway in eukaryotic and prokaryotic cells. In eukaryotes, the textbook view is that glycolysis occurs in the cytosol. However, fusion proteins comprised of two glycolytic enzymes, triosephosphate isomerase (TPI) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH), were found in members of the stramenopiles (diatoms and oomycetes) and shown to possess amino-terminal mitochondrial targeting signals. Here we show that mitochondrial TPI-GAPDH fusion protein genes are widely spread across the known diversity of stramenopiles, including non-photosynthetic species (Bicosoeca sp. and Blastocystis hominis). We also show that TPI-GAPDH fusion genes exist in three cercozoan taxa (Paulinella chromatophora, Thaumatomastix sp. and Mataza hastifera) and an apusozoan protist, Thecamonas trahens. Interestingly, subcellular localization predictions for other glycolytic enzymes in stramenopiles and a cercozoan show that a significant fraction of the glycolytic enzymes in these species have mitochondrial-targeted isoforms. These results suggest that part of the glycolytic pathway occurs inside mitochondria in these organisms, broadening our knowledge of the diversity of mitochondrial metabolism of protists.

Evolution of general transcription factors

Three genes GTF2IRD1, GTF2I, and GTF2IRD2, which encode members of the GTF2I (or TFII-I) family of so-called general transcription factors, were discovered and studied during the last two decades. Chromosome location and similarity of exon-intron structures suggest that the family evolved by duplications. The initial duplication of ancestral proto-GTF2IRD1 gene likely occurred in early vertebrates prior to origin of cartilaginous fish and led to formation of GTF2I (>450 MYA), which was later lost in bony fish but successfully evolved in the land vertebrates. The second duplication event, which created GTF2IRD2, occurred prior to major radiation events of eutherian mammalian evolution (>100 MYA). During recent steps of primate evolution there was another duplication which led to formation of GTF2IRD2B (<4 MYA). Two latest duplications were coupled with inversions. Genes belonging to the family have several highly conservative repeats which are implicated in DNA binding. Phylogenetic analysis of the repeats revealed a pattern of intragenic duplications, deletions and substitutions which led to diversification of the genes and proteins. Distribution of statistically rare atypical substitutions (p ≤ 0.01) sheds some light on structural differentiation of repeats and hence evolution of the genes. The atypical substitutions are often located on secondary structures joining α-helices and affect 3D arrangement of the protein globule. Such substitutions are commonly traced at the early stages of evolution in Tetrapoda, Amniota, and Mammalia.

Development of a viral biopesticide for the control of the Guatemala potato tuber moth Tecia solanivora

The Guatemala potato tuber moth Tecia solanivora (Povolny) (Lep. Gelechiidae) is an invasive species from Mesoamerica that has considerably extended its distribution area in recent decades. While this species is considered to be a major potato pest in Venezuela, Colombia, and Ecuador, currently no specific control methods are available for farmers. To address this issue we developed a biopesticide formulation to be used in integrated pest management of T. solanivora, following three steps. First, search for entomopathogenic viruses were carried out through extensive bioprospections in 12 countries worldwide. As a result, new Phthorimaea operculella granulovirus (PhopGV) isolates were found in T. solanivora and five other gelechid species. Second, twenty PhopGV isolates, including both previously known and newly found isolates, were genetically and/or biologically characterized in order to choose the best candidate for a biopesticide formulation. Sequence data were obtained for the ecdysteroid UDP-glucosyltransferase (egt) gene, a single copy gene known to play a role in pathogenicity. Three different sizes (1086, 1305 and 1353 bp) of egt were found among the virus isolates analyzed. Unexpectedly, no obvious correlation between egt size and pathogenicity was found. Bioassays on T. solanivora neonates showed a maximum of a 14-fold difference in pathogenicity among the eight PhopGV isolates tested. The most pathogenic PhopGV isolate, JLZ9f, had a medium lethal concentration (LC(50)) of 10 viral occlusion bodies per square mm of consumed tuber skin. Third, we tested biopesticide dust formulations by mixing a dry carrier (calcium carbonate) with different adjuvants (magnesium chloride or an optical brightener or soya lecithin) and different specific amounts of JLZ9f. During laboratory experiments, satisfactory control of the pest (>98% larva mortality compared to untreated control) was achieved with a formulation containing 10 macerated JLZ9f-dead T. solanivora larvae per kg of calcium carbonate mixed with 50 mL/kg of soya lecithin. The final product provides an interesting alternative to chemical pesticides for Andean farmers affected by this potato pest.

Strong purifying selection in endogenous retroviruses in the saltwater crocodile (Crocodylus porosus) in the Northern Territory of Australia

Background

Endogenous retroviruses (ERVs) are remnants of exogenous retroviruses that have integrated into the nuclear DNA of a germ-line cell. Here we present the results of a survey into the ERV complement of Crocodylus porosus, the saltwater crocodile, representing 45 individuals from 17 sampling locations in the Northern Territory of Australia. These retroelements were compared with published ERVs from other species of Crocodylia (Crocodilians; alligators, caimans, gharials and crocodiles) as well as representatives from other vertebrates. This study represents one of the first in-depth studies of ERVs within a single reptilian species shedding light on the diversity of ERVs and proliferation mechanisms in crocodilians.

Results

Analyses of the retroviral pro-pol gene region have corroborated the presence of two major clades of ERVs in C. porosus and revealed 18 potentially functional fragments out of the 227 recovered that encode intact pro-pol ORFs. Interestingly, we have identified some patterns of diversification among those ERVs as well as a novel sequence that suggests the presence of an additional retroviral genus in C. porosus. In addition, considerable diversity but low genetic divergence within one of the C. porosus ERV lineages was identified.

Conclusions

We propose that the ERV complement of C. porosus has come about through a combination of recent infections and replication of ancestral ERVs. Strong purifying selection acting on these clades suggests that this activity is recent or still occurring in the genome of this species. The discovery of potentially functional elements is an interesting development that warrants further investigation.

LEC-2, a highly variable lectin in the lichen Peltigera membranacea

Lectins are a diverse group of carbohydrate binding proteins often involved in cellular interactions. A lectin gene, lec-2, was identified in the mycobiont of the lichen Peltigera membranacea. Sequencing of lec-2 open reading frames from 21 individual samples showed an unexpectedly high level of polymorphism in the deduced protein (LEC-2), which was sorted into nine haplotypes based on amino acid sequence. Calculations showed that the rates of nonsynonymous versus synonymous nucleotide substitutions deviated significantly from the null hypothesis of neutrality, indicating strong positive selection. Molecular modeling revealed that most amino acid replacements were around the putative carbohydrate-binding pocket, indicating changes in ligand binding. Lectins have been thought to be involved in the recognition of photobiont partners in lichen symbioses, and the hypothesis that positive selection of LEC-2 is driven by variation in the Nostoc photobiont partner was tested by comparing mycobiont LEC-2 haplotypes and photobiont genotypes, as represented by the rbcLX region. It was not possible to pair up the two types of marker sequences without conflicts, suggesting that positive selection of LEC-2 was not due to variation in photobiont partners.

A method of alignment masking for refining the phylogenetic signal of multiple sequence alignments

Inaccurate inference of positional homologies in multiple sequence alignments and systematic errors introduced by alignment heuristics obfuscate phylogenetic inference. Alignment masking, the elimination of phylogenetically uninformative or misleading sites from an alignment before phylogenetic analysis, is a common practice in phylogenetic analysis. Although masking is often done manually, automated methods are necessary to handle the much larger data sets being prepared today. In this study, we introduce the concept of subsplits and demonstrate their use in extracting phylogenetic signal from alignments. We design a clustering approach for alignment masking where each cluster contains similar columns-similarity being defined on the basis of compatible subsplits; our approach then identifies noisy clusters and eliminates them. Trees inferred from the columns in the retained clusters are found to be topologically closer to the reference trees. We test our method on numerous standard benchmarks (both synthetic and biological data sets) and compare its performance with other methods of alignment masking. We find that our method can eliminate sites more accurately than other methods, particularly on divergent data, and can improve the topologies of the inferred trees in likelihood-based analyses. Software available upon request from the author.

Characterisation of four major histocompatibility complex class II genes of the koala (Phascolarctos cinereus)

Major histocompatibility complex (MHC) class II molecules have an integral role in the adaptive immune response, as they bind and present antigenic peptides to T helper lymphocytes. In this study of koalas, species-specific primers were designed to amplify exon 2 of the MHC class II DA and DB genes, which contain much of the peptide-binding regions of the α and β chains. A total of two DA α1 domain variants and eight DA β1 (DAB), three DB α1 and five DB β1 variants were amplified from 20 koalas from two free-living populations from South East Queensland and the Port Macquarie region in northern New South Wales. We detected greater variation in the β1 than in the α1 domains as well as evidence of positive selection in DAB. The present study provides a springboard to future investigation of the role of MHC in disease susceptibility in koalas.

The independent prokaryotic origins of eukaryotic fructose-1, 6-bisphosphatase and sedoheptulose-1, 7-bisphosphatase and the implications of their origins for the evolution of eukaryotic Calvin cycle

BACKGROUND

In the Calvin cycle of eubacteria, the dephosphorylations of both fructose-1, 6-bisphosphate (FBP) and sedoheptulose-1, 7-bisphosphate (SBP) are catalyzed by the same bifunctional enzyme: fructose-1, 6-bisphosphatase/sedoheptulose-1, 7-bisphosphatase (F/SBPase), while in that of eukaryotic chloroplasts by two distinct enzymes: chloroplastic fructose-1, 6-bisphosphatase (FBPase) and sedoheptulose-1, 7-bisphosphatase (SBPase), respectively. It was proposed that these two eukaryotic enzymes arose from the divergence of a common ancestral eubacterial bifunctional F/SBPase of mitochondrial origin. However, no specific affinity between SBPase and eubacterial FBPase or F/SBPase can be observed in the previous phylogenetic analyses, and it is hard to explain why SBPase and/or F/SBPase are/is absent from most extant nonphotosynthetic eukaryotes according to this scenario.

RESULTS

Domain analysis indicated that eubacterial F/SBPase of two different resources contain distinct domains: proteobacterial F/SBPases contain typical FBPase domain, while cyanobacterial F/SBPases possess FBPase_glpX domain. Therefore, like prokaryotic FBPase, eubacterial F/SBPase can also be divided into two evolutionarily distant classes (Class I and II). Phylogenetic analysis based on a much larger taxonomic sampling than previous work revealed that all eukaryotic SBPase cluster together and form a close sister group to the clade of epsilon-proteobacterial Class I FBPase which are gluconeogenesis-specific enzymes, while all eukaryotic chloroplast FBPase group together with eukaryotic cytosolic FBPase and form another distinct clade which then groups with the Class I FBPase of diverse eubacteria. Motif analysis of these enzymes also supports these phylogenetic correlations.

CONCLUSIONS

There are two evolutionarily distant classes of eubacterial bifunctional F/SBPase. Eukaryotic FBPase and SBPase do not diverge from either of them but have two independent origins: SBPase share a common ancestor with the gluconeogenesis-specific Class I FBPase of epsilon-proteobacteria (or probably originated from that of the ancestor of epsilon-proteobacteria), while FBPase arise from Class I FBPase of an unknown kind of eubacteria. During the evolution of SBPase from eubacterial Class I FBPase, the SBP-dephosphorylation activity was acquired through the transition "from specialist to generalist". The evolutionary substitution of the endosymbiotic-origin cyanobacterial bifunctional F/SBPase by the two light-regulated substrate-specific enzymes made the regulation of the Calvin cycle more delicate, which contributed to the evolution of eukaryotic photosynthesis and even the entire photosynthetic eukaryotes.

Phylogenetic Estimation of Timescales Using Ancient DNA: The Effects of Temporal Sampling Scheme and Uncertainty in Sample Ages

In recent years, ancient DNA has increasingly been used for estimating molecular timescales, particularly in studies of substitution rates and demographic histories. Molecular clocks can be calibrated using temporal information from ancient DNA sequences. This information comes from the ages of the ancient samples, which can be estimated by radiocarbon dating the source material or by dating the layers in which the material was deposited. Both methods involve sources of uncertainty. The performance of Bayesian phylogenetic inference depends on the information content of the data set, which includes variation in the DNA sequences and the structure of the sample ages. Various sources of estimation error can reduce our ability to estimate rates and timescales accurately and precisely. We investigated the impact of sample-dating uncertainties on the estimation of evolutionary timescale parameters using the software BEAST. Our analyses involved 11 published data sets and focused on estimates of substitution rate and root age. We show that, provided that samples have been accurately dated and have a broad temporal span, it might be unnecessary to account for sample-dating uncertainty in Bayesian phylogenetic analyses of ancient DNA. We also investigated the sample size and temporal span of the ancient DNA sequences needed to estimate phylogenetic timescales reliably. Our results show that the range of sample ages plays a crucial role in determining the quality of the results but that accurate and precise phylogenetic estimates of timescales can be made even with only a few ancient sequences. These findings have important practical consequences for studies of molecular rates, timescales, and population dynamics.

Diversity across Seasons of Culturable Pseudomonas from a Desiccation Lagoon in Cuatro Cienegas, Mexico

Cuatro Cienegas basin (CCB) is a biodiversity reservoir within the Chihuahuan desert that includes several water systems subject to marked seasonality. While several studies have focused on biodiversity inventories, this is the first study that describes seasonal changes in diversity within the basin. We sampled Pseudomonas populations from a seasonally variable water system at four different sampling dates (August 2003, January 2004, January 2005, and August 2005). A total of 70 Pseudomonas isolates across seasons were obtained, genotyped by fingerprinting (BOX-PCR), and taxonomically characterized by 16S rDNA sequencing. We found 35 unique genotypes, and two numerically dominant lineages (16S rDNA sequences) that made up 64% of the sample: P. cuatrocienegasensis and P. otitidis. We did not recover genotypes across seasons, but lineages reoccurred across seasons; P. cuatrocienegasensis was isolated exclusively in winter, while P. otitidis was only recovered in summer. We statistically show that taxonomic identity of isolates is not independent of the sampling season, and that winter and summer populations are different. In addition to the genetic description of populations, we show exploratory measures of growth rates at different temperatures, suggesting physiological differences between populations. Altogether, the results indicate seasonal changes in diversity of free-living aquatic Pseudomonas populations from CCB.

Identification of novel glutathione transferases in Echinococcus granulosus. An evolutionary perspective

Glutathione transferase enzymes (GSTs) constitute a major detoxification system in helminth parasites and have been related to the modulation of host immune response mechanisms. At least three different GSTs classes have been described in Platyhelminthes: Mu, Sigma and Omega. Mining the genome of Echinococcus multilocularis and the ESTs databases of Taenia solium and E. granulosus identified two new GSTs from the cestode E. granulosus, named EgGST2 and EgGST3. It also revealed that the Omega class of GSTs was absent from the Taenidae family. EgGST2 and EgGST3 are actively expressed in the parasite. In order to know the origin of these new GSTs, in silico analyses were performed. While EgGST2 is classified as belonging to the Sigma class, the data obtained for EgGST3 allowed a less clear interpretation. The study of the evolutionary relatedness based on the C-terminal domain sequence, gene structure conservation and three-dimensional structure predictions, suggests that EgGST3 is derived from the Platyhelminthes' Sigma-class cluster. Interestingly, the N-terminal domain displays some characteristic Omega-class residues, including a Cys residue that is likely to be involved in the catalytic mechanism. We discuss different evolutionary scenarios that could explain the observed patterns.

High influenza a virus infection rates in Mallards bred for hunting in the Camargue, South of France

During the last decade, the role of wildlife in emerging pathogen transmission to domestic animals has often been pointed out. Conversely, far less attention has been paid to pathogen transmission from domestic animals to wildlife. Here, we focus on the case of game restocking, which implies the release of millions of animals worldwide each year. We conducted a 2-year study in the Camargue (Southern France) to investigate the influence of hand-reared Mallard releases on avian influenza virus dynamics in surrounding wildlife. We sampled Mallards (cloacal swabs) from several game duck facilities in 2009 and 2010 before their release. A very high (99%) infection rate caused by an H10N7 strain was detected in the game bird facility we sampled in 2009. We did not detect this strain in shot ducks we sampled, neither during the 2008/2009 nor the 2009/2010 hunting seasons. In 2010 infection rates ranged from 0 to 24% in hand-reared ducks. The 2009 H10N7 strain was fully sequenced. It results from multiple reassortment events between Eurasian low pathogenic strains. Interestingly, H10N7 strains had previously caused human infections in Egypt and Australia. The H10 and N7 segments we sequenced were clearly distinct from the Australian ones but they belonged to the same large cluster as the Egyptian ones. We did not observe any mutation linked to increased virulence, transmission to mammals, or antiviral resistance in the H10N7 strain we identified. Our results indicate that the potential role of hand-reared Mallards in influenza virus epizootics must be taken into account given the likely risk of viral exchange between game bird facilities and wild habitats, owing to duck rearing conditions. Measures implemented to limit transmission from wildlife to domestic animals as well as measures to control transmission from domestic animals to wild ones need to be equally reinforced.

Identification of a novel bat papillomavirus by metagenomics

The discovery of novel viruses in animals expands our knowledge of viral diversity and potentially emerging zoonoses. High-throughput sequencing (HTS) technology gives millions or even billions of sequence reads per run, allowing a comprehensive survey of the genetic content within a sample without prior nucleic acid amplification. In this study, we screened 156 rectal swab samples from apparently healthy bats (n = 96), pigs (n = 9), cattles (n = 9), stray dogs (n = 11), stray cats (n = 11) and monkeys (n = 20) using a HTS metagenomics approach. The complete genome of a novel papillomavirus (PV), Miniopterus schreibersii papillomavirus type 1 (MscPV1), with L1 of 60% nucleotide identity to Canine papillomavirus (CPV6), was identified in a specimen from a Common Bent-wing Bat (M. schreibersii). It is about 7.5kb in length, with a G+C content of 45.8% and a genomic organization similar to that of other PVs. Despite the higher nucleotide identity between the genomes of MscPV1 and CPV6, maximum-likelihood phylogenetic analysis of the L1 gene sequence showed that MscPV1 and Erethizon dorsatum papillomavirus (EdPV1) are most closely related. Estimated divergence time of MscPV1 from the EdPV1/MscPV1 common ancestor was approximately 60.2–91.9 millions of years ago, inferred under strict clocks using the L1 and E1 genes. The estimates were limited by the lack of reliable calibration points from co-divergence because of possible host shifts. As the nucleotide sequence of this virus only showed limited similarity with that of related animal PVs, the conventional approach of PCR using consensus primers would be unlikely to have detected the novel virus in the sample. Unlike the first bat papillomavirus RaPV1, MscPV1 was found in an asymptomatic bat with no apparent mucosal or skin lesions whereas RaPV1 was detected in the basosquamous carcinoma of a fruit bat Rousettus aegyptiacus. We propose MscPV1 as the first member of the novel Dyolambda-papillomavirus genus.

Insect phylogenomics: results, problems and the impact of matrix composition

In this study, we investigated the relationships among insect orders with a main focus on Polyneoptera (lower Neoptera: roaches, mantids, earwigs, grasshoppers, etc.), and Paraneoptera (thrips, lice, bugs in the wide sense). The relationships between and within these groups of insects are difficult to resolve because only few informative molecular and morphological characters are available. Here, we provide the first phylogenomic expressed sequence tags data ('EST': short sub-sequences from a c(opy) DNA sequence encoding for proteins) for stick insects (Phasmatodea) and webspinners (Embioptera) to complete published EST data. As recent EST datasets are characterized by a heterogeneous distribution of available genes across taxa, we use different rationales to optimize the data matrix composition. Our results suggest a monophyletic origin of Polyneoptera and Eumetabola (Paraneoptera + Holometabola). However, we identified artefacts of tree reconstruction (human louse Pediculus humanus assigned to Odonata (damselflies and dragonflies) or Holometabola (insects with a complete metamorphosis); mayfly genus Baetis nested within Neoptera), which were most probably rooted in a data matrix composition bias due to the inclusion of sequence data of entire proteomes. Until entire proteomes are available for each species in phylogenomic analyses, this potential pitfall should be carefully considered.

A comparative analysis of Y chromosome and mtDNA phylogenies of the Hylobates gibbons

Background

The evolutionary relationships of closely related species have long been of interest to biologists since these species experienced different evolutionary processes in a relatively short period of time. Comparison of phylogenies inferred from DNA sequences with differing inheritance patterns, such as mitochondrial, autosomal, and X and Y chromosomal loci, can provide more comprehensive inferences of the evolutionary histories of species. Gibbons, especially the genus Hylobates, are particularly intriguing as they consist of multiple closely related species which emerged rapidly and live in close geographic proximity. Our current understanding of relationships among Hylobates species is largely based on data from the maternally-inherited mitochondrial DNAs (mtDNAs).

Results

To infer the paternal histories of gibbon taxa, we sequenced multiple Y chromosomal loci from 26 gibbons representing 10 species. As expected, we find levels of sequence variation some five times lower than observed for the mitochondrial genome (mtgenome). Although our Y chromosome phylogenetic tree shows relatively low resolution compared to the mtgenome tree, our results are consistent with the monophyly of gibbon genera suggested by the mtgenome tree. In a comparison of the molecular dating of divergences and on the branching patterns of phylogeny trees between mtgenome and Y chromosome data, we found: 1) the inferred divergence estimates were more recent for the Y chromosome than for the mtgenome, 2) the species H. lar and H. pileatus are monophyletic in the mtgenome phylogeny, respectively, but a H. pileatus individual falls into the H. lar Y chromosome clade.

Conclusions

Based on the ~6.4 kb of Y chromosomal DNA sequence data generated for each of the 26 individuals in this study, we provide molecular inferences on gibbon and particularly on Hylobates evolution complementary to those from mtDNA data. Overall, our results illustrate the utility of comparative studies of loci with different inheritance patterns for investigating potential sex specific processes on the evolutionary histories of closely related taxa, and emphasize the need for further sampling of gibbons of known provenance.

Host selection and parasite infection in Aedes taeniorhynchus, endemic disease vector in the Galápagos Islands

Host selection in blood-sucking arthropods has important evolutionary and ecological implications for the transmission dynamics, distribution and host-specificity of the parasites they transmit. The black salt-marsh mosquito (Aedes taeniorhynchus Wiedemann) is distributed throughout tropical to temperate coastal zones in the Americas, and continental populations are primarily mammalphilic. It is the only indigenous mosquito in the Galápagos Islands, having colonised the archipelago around 200,000 years ago, potentially adapting its host selection, and in the process, altering the dynamics of vector mediated pathogen interactions in the archipelago. Here, we use blood-meal analysis and PCR-based parasite screening approach to determine the blood-feeding patterns of A. taeniorhynchus in the Galápagos Islands and identify potential parasite transmission with which this mosquito could be involved. Our results show that A. taeniorhynchus feeds equally on mammals and reptiles, and only one avian sample was observed in 190 successful PCR amplifications from blood meals. However, we detected endemic filarial worms and Haemoproteus parasites known to infect various Galápagos bird species in mosquito thoraces, suggesting that feeding on birds must occur at low frequency, and that A. taeniorhynchus may play a role in maintaining some avian vector-borne pathogens, although more work is needed to explore this possibility. We also isolated three different DNA sequences corresponding to hemogregarine parasites of the genus Hepatozoon from mosquito and iguana blood samples, suggesting that more than one species of Hepatozoon parasites are present in Galápagos. Phylogenetic analysis of Hepatozoon 18sRNA sequences indicates that A. taeniorhynchus may have facilitated a recent breakdown in host-species association of formerly isolated Hepatozoon spp. infecting the reptile populations in the Galápagos Islands.