PAL: an object-oriented programming library for molecular evolution and phylogenetics

SubRecon: ancestral reconstruction of amino acid substitutions along a branch in a phylogeny

Summary

Existing ancestral sequence reconstruction techniques are ill-suited to investigating substitutions on a single branch of interest. We present SubRecon, an implementation of a hybrid technique integrating joint and marginal reconstruction for protein sequence data. SubRecon calculates the joint probability of states at adjacent internal nodes in a phylogeny, i.e. how the state has changed along a branch. This does not condition on states at other internal nodes and includes site rate variation. Simulation experiments show the technique to be accurate and powerful. SubRecon has a user-friendly command line interface and produces concise output that is intuitive yet suitable for subsequent parsing in an automated pipeline.

Availability and implementation

SubRecon is platform independent, requiring Java v1.8 or above. Source code, installation instructions and an example dataset are freely available under the Apache 2.0 license at https://github.com/chrismonit/SubRecon.

Characterization and evolution of porcine deltacoronavirus in the United States

Porcine deltacoronavirus (PDCoV) was identified in multiple states across the United States (US) in 2014. In this study, we investigate the presence of PDCoV in diagnostic samples, which were further categorized by case identification (ID), and the association between occurrence, age, specimen and location between March and September 2014. Approximately, 7% of the case IDs submitted from the US were positive for PDCoV. Specimens were categorized into eight groups, and the univariate analysis indicated that oral fluids had 1.89 times higher odds of detecting PDCoV compared to feces. While the 43-56 day age group had the highest percentage of PDCoV positives (8.4%), the univariate analysis indicated no significant differences between age groups. However, multivariable analysis for age adjusted by specimen indicated the >147 day age group had 59% lower odds than suckling pigs of being positive for PDCoV. The percentage of PDCoV in diagnostic samples decreased to <1% in September 2014. In addition, 19 complete PDCoV genomes were sequenced, and Bayesian analysis was conducted to estimate the emergence of the US clade. The evolutionary rate of the PDCoV genome is estimated to be 3.8×10(-4) substitutions/site/year (2.3×10(-4)-5.4×10(-4), 95% HPD). Our results indicate that oral fluids continue to be a valuable specimen to monitor swineherd health, and PDCoV has been circulating in the US prior to 2014.

Genetic diversity of Candidatus Liberibacter asiaticus based on two hypervariable effector genes in Thailand

Huanglongbing (HLB), also known as citrus greening, is one of the most destructive diseases of citrus worldwide. HLB is associated with three species of 'Candidatus Liberibacter' with 'Ca. L. asiaticus' (Las) being the most widely distributed around the world, and the only species detected in Thailand. To understand the genetic diversity of Las bacteria in Thailand, we evaluated two closely-related effector genes, lasAI and lasAII, found within the Las prophages from 239 infected citrus and 55 infected psyllid samples collected from different provinces in Thailand. The results indicated that most of the Las-infected samples collected from Thailand contained at least one prophage sequence with 48.29% containing prophage 1 (FP1), 63.26% containing prophage 2 (FP2), and 19.38% containing both prophages. Interestingly, FP2 was found to be the predominant population in Las-infected citrus samples while Las-infected psyllids contained primarily FP1. The multiple banding patterns that resulted from amplification of lasAI imply extensive variation exists within the full and partial repeat sequence while the single band from lasAII indicates a low amount of variation within the repeat sequence. Phylogenetic analysis of Las-infected samples from 22 provinces in Thailand suggested that the bacterial pathogen may have been introduced to Thailand from China and the Philippines. This is the first report evaluating the genetic variation of a large population of Ca. L. asiaticus infected samples in Thailand using the two effector genes from Las prophage regions.

Evolutionary and sequence-based relationships in bacterial AdoMet-dependent non-coding RNA methyltransferases

Background

RNA post-transcriptional modification is an exciting field of research that has evidenced this editing process as a sophisticated epigenetic mechanism to fine tune the ribosome function and to control gene expression. Although tRNA modifications seem to be more relevant for the ribosome function and cell physiology as a whole, some rRNA modifications have also been seen to play pivotal roles, essentially those located in central ribosome regions. RNA methylation at nucleobases and ribose moieties of nucleotides appear to frequently modulate its chemistry and structure. RNA methyltransferases comprise a superfamily of highly specialized enzymes that accomplish a wide variety of modifications. These enzymes exhibit a poor degree of sequence similarity in spite of using a common reaction cofactor and modifying the same substrate type.

Results

Relationships and lineages of RNA methyltransferases have been extensively discussed, but no consensus has been reached. To shed light on this topic, we performed amino acid and codon-based sequence analyses to determine phylogenetic relationships and molecular evolution. We found that most Class I RNA MTases are evolutionarily related to protein and cofactor/vitamin biosynthesis methyltransferases. Additionally, we found that at least nine lineages explain the diversity of RNA MTases. We evidenced that RNA methyltransferases have high content of polar and positively charged amino acid, which coincides with the electrochemistry of their substrates.

Conclusions

After studying almost 12,000 bacterial genomes and 2,000 patho-pangenomes, we revealed that molecular evolution of Class I methyltransferases matches the different rates of synonymous and non-synonymous substitutions along the coding region. Consequently, evolution on Class I methyltransferases selects against amino acid changes affecting the structure conformation.

MiningABs: mining associated biomarkers across multi-connected gene expression datasets

Background

Human disease often arises as a consequence of alterations in a set of associated genes rather than alterations to a set of unassociated individual genes. Most previous microarray-based meta-analyses identified disease-associated genes or biomarkers independent of genetic interactions. Therefore, in this study, we present the first meta-analysis method capable of taking gene combination effects into account to efficiently identify associated biomarkers (ABs) across different microarray platforms.

Results

We propose a new meta-analysis approach called MiningABs to mine ABs across different array-based datasets. The similarity between paired probe sequences is quantified as a bridge to connect these datasets together. The ABs can be subsequently identified from an “improved” common logit model (c-LM) by combining several sibling-like LMs in a heuristic genetic algorithm selection process. Our approach is evaluated with two sets of gene expression datasets: i) 4 esophageal squamous cell carcinoma and ii) 3 hepatocellular carcinoma datasets. Based on an unbiased reciprocal test, we demonstrate that each gene in a group of ABs is required to maintain high cancer sample classification accuracy, and we observe that ABs are not limited to genes common to all platforms. Investigating the ABs using Gene Ontology (GO) enrichment, literature survey, and network analyses indicated that our ABs are not only strongly related to cancer development but also highly connected in a diverse network of biological interactions.

Conclusions

The proposed meta-analysis method called MiningABs is able to efficiently identify ABs from different independently performed array-based datasets, and we show its validity in cancer biology via GO enrichment, literature survey and network analyses. We postulate that the ABs may facilitate novel target and drug discovery, leading to improved clinical treatment. Java source code, tutorial, example and related materials are available at “http://sourceforge.net/projects/miningabs/”.

VP6 genetic diversity, reassortment, intragenic recombination and classification of rotavirus B in American and Japanese pigs

Rotavirus B (RVB) has been identified as a causative agent of diarrhea in rats, humans, cattle, lambs, and swine. Recently, 20 RVB VP7 genotypes were determined based on an 80% nucleotide percent cut-off value. In this study, we sequenced the RVB VP6 gene segment from 80 RVB positive swine samples from the United States and Japan. Phylogenetic analyses, using the 30 available RVB VP6 sequences from GenBank and our 80 novel RVB VP6 sequences, revealed a large genetic diversity of RVB strains, mainly in pigs. For classification purposes, pairwise identity frequency analyses suggested an 81% nucleotide percent cut-off value, resulting in 13 RVB VP6 (I) genotypes. In addition, an intragenic recombinant RVB VP6 segment was identified from Japan. Furthermore, the data indicates frequent reassortment events occurred between the porcine RVB VP7 and VP6 gene segments.

Molecular evolution of calcification genes in morphologically similar but phylogenetically unrelated scleractinian corals

Molecular phylogenies of scleractinian corals often fail to agree with traditional phylogenies derived from morphological characters. These discrepancies are generally attributed to non-homologous or morphologically plastic characters used in taxonomic descriptions. Consequently, morphological convergence of coral skeletons among phylogenetically unrelated groups is considered to be the major evolutionary process confounding molecular and morphological hypotheses. A strategy that may help identify cases of convergence and/or diversification in coral morphology is to compare phylogenies of existing "neutral" genetic markers used to estimate genealogic phylogenetic history with phylogenies generated from non-neutral genes involved in calcification (biomineralization). We tested the hypothesis that differences among calcification gene phylogenies with respect to the "neutral" trees may represent convergent or divergent functional strategies among calcification gene proteins that may correlate to aspects of coral skeletal morphology. Partial sequences of two nuclear genes previously determined to be involved in the calcification process in corals, "Cnidaria-III" membrane-bound/secreted α-carbonic anhydrase (CIII-MBSα-CA) and bone morphogenic protein (BMP) 2/4, were PCR-amplified, cloned and sequenced from 31 scleractinian coral species in 26 genera and 9 families. For comparison, "neutral" gene phylogenies were generated from sequences from two protein-coding "non-calcification" genes, one nuclear (β-tubulin) and one mitochondrial (cytochrome b), from the same individuals. Cloned CIII-MBSα-CA sequences were found to be non-neutral, and phylogenetic analyses revealed CIII-MBSα-CAs to exhibit a complex evolutionary history with clones distributed between at least 2 putative gene copies. However, for several coral taxa only one gene copy was recovered. With CIII-MBSα-CA, several recovered clades grouped taxa that differed from the "non-calcification" loci. In some cases, these taxa shared aspects of their skeletal morphology (i.e., convergence or diversification relative to the "non-calcification" loci), but in other cases they did not. For example, the "non-calcification" loci recovered Atlantic and Pacific mussids as separate evolutionary lineages, whereas with CIII-MBSα-CA, clones of two species of Atlantic mussids (Isophyllia sinuosa and Mycetophyllia sp.) and two species of Pacific mussids (Acanthastrea echinata and Lobophyllia hemprichii) were united in a distinct clade (except for one individual of Mycetophyllia). However, this clade also contained other taxa which were not unambiguously correlated with morphological features. BMP2/4 also contained clones that likely represent different gene copies. However, many of the sequences showed no significant deviation from neutrality, and reconstructed phylogenies were similar to the "non-calcification" tree topologies with a few exceptions. Although individual calcification genes are unlikely to precisely explain the diverse morphological features exhibited by scleractinian corals, this study demonstrates an approach for identifying cases where morphological taxonomy may have been misled by convergent and/or divergent molecular evolutionary processes in corals. Studies such as this may help illuminate our understanding of the likely complex evolution of genes involved in the calcification process, and enhance our knowledge of the natural history and biodiversity within this central ecological group.

DoMosaics: software for domain arrangement visualization and domain-centric analysis of proteins

DoMosaics is an application that unifies protein domain annotation, domain arrangement analysis and visualization in a single tool. It simplifies the analysis of protein families by consolidating disjunct procedures based on often inconvenient command-line applications and complex analysis tools. It provides a simple user interface with access to domain annotation services such as InterProScan or a local HMMER installation, and can be used to compare, analyze and visualize the evolution of domain architectures.

AVAILABILITY AND IMPLEMENTATION

DoMosaics is licensed under theApache License, Version 2.0, and binaries can be freely obtained from www.domosaics.net.

Evolutionary Modeling of Genotype-Phenotype Associations, and Application to Primate Coding and Non-coding mtDNA Rate Variation

Variation in substitution rates across a phylogeny can be indicative of shifts in the evolutionary dynamics of a protein or non-protein coding regions. One way to understand these signals is to seek the phenotypic correlates of rate variation. Here, we extended a previously published likelihood method designed to detect evolutionary associations between genotypic evolutionary rate and phenotype over a phylogeny. In simulation with two discrete categories of phenotype, the method has a low false-positive rate and detects greater than 80% of true-positives with a tree length of three or greater and a three-fold or greater change in substitution rate given the phenotype. In addition, we successfully extend the test from two to four phenotype categories and evaluated its performance. We then applied the method to two major hypotheses for rate variation in the mitochondrial genome of primates-longevity and generation time as well as body mass which is correlated with many aspects of life history-using three categories of phenotype through discretization of continuous values. Similar to previous results for mammals, we find that the majority of mitochondrial protein-coding genes show associations consistent with the longevity and body mass predictions and that the predominant signal of association comes from the third codon position. We also found a significant association between maximum lifespan and the evolutionary rate of the control region of the mtDNA. In contrast, 24 protein-coding genes from the nuclear genome do not show a consistent pattern of association, which is inconsistent with the generation time hypothesis. These results show the extended method can robustly identify genotype-phenotype associations up to at least four phenotypic categories, and demonstrate the successful application of the method to study factors affecting neutral evolutionary rate in protein-coding and non-coding loci.

More than one way to produce protein diversity: duplication and limited alternative splicing of an adhesion molecule gene in basal arthropods

Exon duplication and alternative splicing evolved multiple times in metazoa and are of overall importance in shaping genomes and allowing organisms to produce many fold more proteins than there are genes in the genome. No other example is as striking as the one of the Down syndrome cell adhesion molecule (Dscam) of insects and crustaceans (pancrustaceans) involved in the nervous system differentiation and in the immune system. To elucidate the evolutionary history of this extraordinary gene, we investigated Dscam homologs in two basal arthropods, the myriapod Strigamia maritima and the chelicerate Ixodes scapularis. In both, Dscam diversified extensively by whole gene duplications resulting in multigene expansions. Within some of the S. maritima genes, exons coding for one of the immunoglobulin domains (Ig7) duplicated and are mutually exclusively alternatively spliced. Our results suggest that Dscam diversification was selected independently in chelicerates, myriapods, and pancrustaceans and that the usage of Dscam diversity by immune cells evolved for the first time in basal arthropods. We propose an evolutionary scenario for the appearance of the highly variable Dscam gene of pancrustaceans, adding to the understanding of how alternative splicing, exon, and gene duplication contribute to create molecular diversity associated with potentially new cellular functions.

Analysis of the WUSCHEL-RELATED HOMEOBOX gene family in the conifer picea abies reveals extensive conservation as well as dynamic patterns

BACKGROUND

Members of the WUSCHEL-RELATED HOMEOBOX (WOX) gene family have important functions during all stages of plant development and have been implicated in the development of morphological novelties during evolution. Most studies have examined the function of these genes in angiosperms and very little is known from other plant species.

RESULTS

In this study we examine the presence and expression of WOX genes in the conifer Picea abies. We have cloned 11 WOX genes from both mRNA and genomic DNA and examined their phylogenetic relationship to WOX genes from other species as well as their expression during somatic embryogenesis and in adult tissues.

CONCLUSIONS

Our study shows that all major radiations within the WOX gene family took place before the angiosperm-gymnosperm split and that there has been a recent expansion within the intermediate clade in the Pinaceae family. Furthermore, we show that the genes from the intermediate clade are preferentially expressed during embryo development in Picea abies. Our data also indicates that there are clear orthologs of both WUS and WOX5 present in the P. abies genome.

The effects of model choice and mitigating bias on the ribosomal tree of life

Deep-level relationships within Bacteria, Archaea, and Eukarya as well as the relationships of these three domains to each other require resolution. The ribosomal machinery, universal to all cellular life, represents a protein repertoire resistant to horizontal gene transfer, which provides a largely congruent signal necessary for reconstructing a tree suitable as a backbone for life's reticulate history. Here, we generate a ribosomal tree of life from a robust taxonomic sampling of Bacteria, Archaea, and Eukarya to elucidate deep-level intra-domain and inter-domain relationships. Lack of phylogenetic information and systematic errors caused by inadequate models (that cannot account for substitution rate or compositional heterogeneities) or improper model selection compound conflicting phylogenetic signals from HGT and/or paralogy. Thus, we tested several models of varying sophistication on three different datasets, performed removal of fast-evolving or long-branched Archaea and Eukarya, and employed three different strategies to remove compositional heterogeneity to examine their effects on the topological outcome. Our results support a two-domain topology for the tree of life, where Eukarya emerges from within Archaea as sister to a Korarchaeota/Thaumarchaeota (KT) or Crenarchaeota/KT clade for all models under all or at least one of the strategies employed. Taxonomic manipulation allows single-matrix and certain mixture models to vacillate between two-domain and three-domain phylogenies. We find that models vary in their ability to resolve different areas of the tree of life, which does not necessarily correlate with model complexity. For example, both single-matrix and some mixture models recover monophyletic Crenarchaeota and Euryarchaeota archaeal phyla. In contrast, the most sophisticated model recovers a paraphyletic Euryarchaeota but detects two large clades that comprise the Bacteria, which were recovered separately but never together in the other models. Overall, models recovered consistent topologies despite dataset modifications due to the removal of compositional bias, which reflects either ineffective bias reduction or robust datasets that allow models to overcome reconstruction artifacts. We recommend a comparative approach for evolutionary models to identify model weaknesses as well as consensus relationships.

Molecular characterization of three GnRH receptor paralogs in the European eel, Anguilla anguilla: tissue-distribution and changes in transcript abundance during artificially induced sexual development

Gonadotropin-releasing hormone receptor (GnRH-R) activation stimulates synthesis and release of gonadotropins in the vertebrate pituitary and also mediates other processes both in the brain and in peripheral tissues. To better understand the differential function of multiple GnRH-R paralogs, three GnRH-R genes (gnrhr1a, 1b, and 2) were isolated and characterized in the European eel. All three gnrhr genes were expressed in the brain and pituitary of pre-pubertal eels, and also in several peripheral tissues, notably gills and kidneys. During hormonally induced sexual maturation, pituitary expression of gnrhr1a (female) and gnrhr2 (male and female) was up-regulated in parallel with gonad development. In the brain, a clear regulation during maturation was seen only for gnrhr2 in the midbrain, with highest levels recorded during early vitellogenesis. These data suggest that GnRH-R2 is the likely hypophysiotropic GnRH-R in male eel, while both GnRH-R1a and GnRH-R2 seems to play this role in female eels.

Examining new phylogenetic markers to uncover the evolutionary history of early-diverging fungi: comparing MCM7, TSR1 and rRNA genes for single- and multi-gene analyses of the Kickxellomycotina

The recently recognised protein-coding genes MCM7 and TSR1 have shown significant promise for phylogenetic resolution within the Ascomycota and Basidiomycota, but have remained unexamined within other fungal groups (except for Mucorales). We designed and tested primers to amplify these genes across early-diverging fungal clades, with emphasis on the Kickxellomycotina, zygomycetous fungi with characteristic flared septal walls forming pores with lenticular plugs. Phylogenetic tree resolution and congruence with MCM7 and TSR1 were compared against those inferred with nuclear small (SSU) and large subunit (LSU) rRNA genes. We also combined MCM7 and TSR1 data with the rDNA data to create 3- and 4-gene trees of the Kickxellomycotina that help to resolve evolutionary relationships among and within the core clades of this subphylum. Phylogenetic inference suggests that Barbatospora, Orphella, Ramicandelaber and Spiromyces may represent unique lineages. It is suggested that these markers may be more broadly useful for phylogenetic studies among other groups of early-diverging fungi.

The genetic architecture of degenerin/epithelial sodium channels in Drosophila

Degenerin/epithelial sodium channels (DEG/ENaC) represent a large family of animal-specific membrane proteins. Although the physiological functions of most family members are not known, some have been shown to act as nonvoltage gated, amiloride-sensitive sodium channels. The DEG/ENaC family is exceptionally large in genomes of Drosophila species relative to vertebrates and other insects. To elucidate the evolutionary history of the DEG/ENaC family in Drosophila, we took advantage of the genomic and genetic information available for 12 Drosophila species that represent all the major species groups in the Drosophila clade. We have identified 31 family members (termed pickpocket genes) in Drosophila melanogaster, which can be divided into six subfamilies, which are represented in all 12 species. Structure prediction analyses suggested that some subunits evolved unique structural features in the large extracellular domain, possibly supporting mechanosensory functions. This finding is further supported by experimental data that show that both ppk1 and ppk26 are expressed in multidendritic neurons, which can sense mechanical nociceptive stimuli in larvae. We also identified representative genes from five of the six DEG/ENaC subfamilies in a mosquito genome, suggesting that the core DEG/ENaC subfamilies were already present early in the dipteran radiation. Spatial and temporal analyses of expression patterns of the various pickpocket genes indicated that paralogous genes often show very different expression patterns, possibly indicating that gene duplication events have led to new physiological or cellular functions rather than redundancy. In summary, our analyses support a rapid early diversification of the DEG/ENaC family in Diptera followed by physiological and/or cellular specialization. Some members of the family may have diversified to support the physiological functions of a yet unknown class of ligands.

The evolutionary origins of chordate hematopoiesis and vertebrate endothelia

The vertebrate circulatory system is the most complex vascular system among those of metazoans, with key innovations including a multi-chambered heart and highly specialized blood cells. Invertebrate vessels, on the other hand, consist of hemal spaces between the basal laminae of epithelia. How the evolutionary transition from an invertebrate-type system to the complex vertebrate one occurred is, however, poorly understood. We investigate here the development of the cardiovascular system of the cephalochordate amphioxus Branchiostoma lanceolatum in order to gain insight into the origin of the vertebrate cardiovascular system. The cardiac markers Hand, Csx (Nkx2-5) and Tbx4/5 reveal a broad cardiac-like domain in amphioxus; such a decentralized organization during development parallels that seen in the adult anatomy. Our data therefore support the hypothesis that amphioxus never possessed a proper heart, even transiently during development. We also define a putative hematopoietic domain, supported by the expression of the hematopoietic markers Scl and Pdvegfr. We show that this area is closed to the dorsal aorta anlages, partially linked to excretory tissues, and that its development is regulated by retinoic acid, thus recalling the aorta-gonads-mesonephros (AGM) area of vertebrates. This region probably produces Pdvegfr+ hemal cells, with an important role in amphioxus vessel formation, since treatments with an inhibitor of PDGFR/VEGFR lead to a decrease of Laminin in the basal laminae of developing vessels. Our results point to a chordate origin of hematopoiesis in an AGM-like area from where hemal Pdvegfr+ cells are produced. These Pdvegfr+ cells probably resemble the ancestral chordate blood cells from which the vertebrate endothelium later originated.

Evolutionary study of duplications of the miRNA machinery in aphids associated with striking rate acceleration and changes in expression profiles

Background

The sequencing of the genome of the pea aphid Acyrthosiphon pisum revealed an unusual expansion of the miRNA machinery, with two argonaute-1, two dicer-1 and four pasha gene copies. In this report, we have undertaken a deeper evolutionary analysis of the phylogenetic timing of these gene duplications and of the associated selective pressures by sequencing the two copies of ago-1 and dcr-1 in different aphid species of the subfamily Aphidinae. We have also carried out an analysis of the expression of both copies of ago-1 and dcr-1 by semi-quantitative PCR in different morphs of the pea aphid life cycle.

Results

The analysis has shown that the duplication of ago-1 occurred in an ancestor of the subfamily Aphidinae while the duplication of dcr-1 appears to be more recent. Besides, it has confirmed a pattern of one conserved copy and one accelerated copy for both genes, and has revealed the action of positive selection on several regions of the fast-evolving ago-1b. On the other hand, the semi-quantitative PCR experiments have revealed a differential expression of these genes between the morphs of the parthenogenetic and the sexual phases of Acyrthosiphon pisum.

Conclusions

The discovery of these gene duplications in the miRNA machinery of aphids opens new perspectives of research about the regulation of gene expression in these insects. Accelerated evolution, positive selection and differential expression affecting some of the copies of these genes suggests the possibility of a neofunctionalization of these duplicates, which might play a role in the display of the striking phenotypic plasticity of aphids.

A mitogenomic re-evaluation of the bdelloid phylogeny and relationships among the Syndermata

Molecular and morphological data regarding the relationships among the three classes of Rotifera (Bdelloidea, Seisonidea, and Monogononta) and the phylum Acanthocephala are inconclusive. In particular, Bdelloidea lacks molecular-based phylogenetic appraisal. I obtained coding sequences from the mitochondrial genomes of twelve bdelloids and two monogononts to explore the molecular phylogeny of Bdelloidea and provide insight into the relationships among lineages of Syndermata (Rotifera + Acanthocephala). With additional sequences taken from previously published mitochondrial genomes, the total dataset included nine species of bdelloids, three species of monogononts, and two species of acanthocephalans. A supermatrix of these 10-12 mitochondrial proteins consistently recovered a bdelloid phylogeny that questions the validity of a generally accepted classification scheme despite different methods of inference and various parameter adjustments. Specifically, results showed that neither the family Philodinidae nor the order Philodinida are monophyletic as currently defined. The application of a similar analytical strategy to assess syndermate relationships recovered either a tree with Bdelloidea and Monogononta as sister taxa (Eurotatoria) or Bdelloidea and Acanthocephala as sister taxa (Lemniscea). Both outgroup choice and method of inference affected the topological outcome emphasizing the need for sequences from more closely related outgroups and more sophisticated methods of analysis that can account for the complexity of the data.

TreeCmp: Comparison of Trees in Polynomial Time

When a phylogenetic reconstruction does not result in one tree but in several, tree metrics permit finding out how far the reconstructed trees are from one another. They also permit to assess the accuracy of a reconstruction if a true tree is known. TreeCmp implements eight metrics that can be calculated in polynomial time for arbitrary (not only bifurcating) trees: four for unrooted (Matching Split metric, which we have recently proposed, Robinson-Foulds, Path Difference, Quartet) and four for rooted trees (Matching Cluster, Robinson-Foulds cluster, Nodal Splitted and Triple). TreeCmp is the first implementation of Matching Split/Cluster metrics and the first efficient and convenient implementation of Nodal Splitted. It allows to compare relatively large trees. We provide an example of the application of TreeCmp to compare the accuracy of ten approaches to phylogenetic reconstruction with trees up to 5000 external nodes, using a measure of accuracy based on normalized similarity between trees.

A genome-wide study of recombination rate variation in Bartonella henselae

Background

Rates of recombination vary by three orders of magnitude in bacteria but the reasons for this variation is unclear. We performed a genome-wide study of recombination rate variation among genes in the intracellular bacterium Bartonella henselae, which has among the lowest estimated ratio of recombination relative to mutation in prokaryotes.

Results

The 1.9 Mb genomes of B. henselae strains IC11, UGA10 and Houston-1 genomes showed only minor gene content variation. Nucleotide sequence divergence levels were less than 1% and the relative rate of recombination to mutation was estimated to 1.1 for the genome overall. Four to eight segments per genome presented significantly enhanced divergences, the most pronounced of which were the virB and trw gene clusters for type IV secretion systems that play essential roles in the infection process. Consistently, multiple recombination events were identified inside these gene clusters. High recombination frequencies were also observed for a gene putatively involved in iron metabolism. A phylogenetic study of this gene in 80 strains of Bartonella quintana, B. henselae and B. grahamii indicated different population structures for each species and revealed horizontal gene transfers across Bartonella species with different host preferences.

Conclusions

Our analysis has shown little novel gene acquisition in B. henselae, indicative of a closed pan-genome, but higher recombination frequencies within the population than previously estimated. We propose that the dramatically increased fixation rate for recombination events at gene clusters for type IV secretion systems is driven by selection for sequence variability.

Basal jawed vertebrate phylogenomics using transcriptomic data from Solexa sequencing

The traditionally accepted relationships among basal jawed vertebrates have been challenged by some molecular phylogenetic analyses based on mitochondrial sequences. Those studies split extant gnathostomes into two monophyletic groups: tetrapods and piscine branch, including Chondrichthyes, Actinopterygii and sarcopterygian fishes. Lungfish and bichir are found in a basal position on the piscine branch. Based on transcriptomes of an armored bichir (Polypterus delhezi) and an African lungfish (Protopterus sp.) we generated, expressed sequences and whole genome sequences available from public databases, we obtained 111 genes to reconstruct the phylogenetic tree of basal jawed vertebrates and estimated their times of divergence. Our phylogenomic study supports the traditional relationship. We found that gnathostomes are divided into Chondrichthyes and the Osteichthyes, both with 100% support values (posterior probabilities and bootstrap values). Chimaeras were found to have a basal position among cartilaginous fishes with a 100% support value. Osteichthyes were divided into Actinopterygii and Sarcopterygii with 100% support value. Lungfish and tetrapods form a monophyletic group with 100% posterior probability. Bichir and two teleost species form a monophyletic group with 100% support value. The previous tree, based on mitochondrial data, was significantly rejected by an approximately unbiased test (AU test, p = 0). The time of divergence between lungfish and tetrapods was estimated to be 391.8 Ma and the divergence of bichir from pufferfish and medaka was estimated to be 330.6 Ma. These estimates closely match the fossil record. In conclusion, our phylogenomic study successfully resolved the relationship of basal jawed vertebrates based on transtriptomes, EST and whole genome sequences.

GeneclusterViz: a tool for conserved gene cluster visualization, exploration and analysis

MOTIVATION

Gene clusters are arrangements of functionally related genes on a chromosome. In bacteria, it is expected that evolutionary pressures would conserve these arrangements due to the functional advantages they provide. Visualization of conserved gene clusters across multiple genomes provides key insights into their evolutionary histories. Therefore, a software tool that enables visualization and functional analyses of gene clusters would be a great asset to the biological research community.

RESULTS

We have developed GeneclusterViz, a Java-based tool that allows for the visualization, exploration and downstream analyses of conserved gene clusters across multiple genomes. GeneclusterViz combines an easy-to-use exploration interface for gene clusters with a host of other analysis features such as multiple sequence alignments, phylogenetic analyses and integration with the KEGG pathway database.

AVAILABILITY

http://biohealth.snu.ac.kr/GeneclusterViz/; http://microbial.informatics.indiana.edu/GeneclusterViz/

Horizontal gene transfer confers fermentative metabolism in the respiratory-deficient plant trypanosomatid Phytomonas serpens

Among trypanosomatids, the genus Phytomonas is the only one specifically adapted to infect plants. These hosts provide a particular habitat with a plentiful supply of carbohydrates. Phytomonas sp. lacks a cytochrome-mediated respiratory chain and Krebs cycle, and ATP production relies predominantly on glycolysis. We have characterised the complete gene encoding a putative pyruvate/indolepyruvate decarboxylase (PDC/IPDC) (548 amino acids) of P. serpens, that displays high amino acid sequence similarity with phytobacteria and Leishmania enzymes. No orthologous PDC/IPDC genes were found in Trypanosoma cruzi or T. brucei. Conservation of the PDC/IPDC gene sequence was verified in 14 Phytomonas isolates. A phylogenetic analysis shows that Phytomonas protein is robustly monophyletic with Leishmania spp. and C. fasciculata enzymes. In the trees this clade appears as a sister group of indolepyruvate decarboxylases of γ-proteobacteria. This supports the proposition that a horizontal gene transfer event from a donor phytobacteria to a recipient ancestral trypanosome has occurred prior to the separation between Phytomonas, Leishmania and Crithidia. We have measured the PDC activity in P. serpens cell extracts. The enzyme has a Km value for pyruvate of 1.4mM. The acquisition of a PDC, a key enzyme in alcoholic fermentation, explains earlier observations that ethanol is one of the major end-products of glucose catabolism under aerobic and anaerobic conditions. This represents an alternative and necessary route to reoxidise part of the NADH produced in the highly demanding glycolytic pathway and highlights the importance of this type of event in metabolic adaptation.

Accounting for alignment uncertainty in phylogenomics

Uncertainty in multiple sequence alignments has a large impact on phylogenetic analyses. Little has been done to evaluate the quality of individual positions in protein sequence alignments, which directly impact the accuracy of phylogenetic trees. Here we describe ZORRO, a probabilistic masking program that accounts for alignment uncertainty by assigning confidence scores to each alignment position. Using the BALIBASE database and in simulation studies, we demonstrate that masking by ZORRO significantly reduces the alignment uncertainty and improves the tree accuracy.

An ancient genomic regulatory block conserved across bilaterians and its dismantling in tetrapods by retrogene replacement

Developmental genes are regulated by complex, distantly located cis-regulatory modules (CRMs), often forming genomic regulatory blocks (GRBs) that are conserved among vertebrates and among insects. We have investigated GRBs associated with Iroquois homeobox genes in 39 metazoans. Despite 600 million years of independent evolution, Iroquois genes are linked to ankyrin-repeat-containing Sowah genes in nearly all studied bilaterians. We show that Iroquois-specific CRMs populate the Sowah locus, suggesting that regulatory constraints underlie the maintenance of the Iroquois-Sowah syntenic block. Surprisingly, tetrapod Sowah orthologs are intronless and not associated with Iroquois; however, teleost and elephant shark data demonstrate that this is a derived feature, and that many Iroquois-CRMs were ancestrally located within Sowah introns. Retroposition, gene, and genome duplication have allowed selective elimination of Sowah exons from the Iroquois regulatory landscape while keeping associated CRMs, resulting in large associated gene deserts. These results highlight the importance of CRMs in imposing constraints to genome architecture, even across large phylogenetic distances, and of gene duplication-mediated genetic redundancy to disentangle these constraints, increasing genomic plasticity.

Characterizing the phylogenetic tree-search problem

Phylogenetic trees are important in many areas of biological research, ranging from systematic studies to the methods used for genome annotation. Finding the best scoring tree under any optimality criterion is an NP-hard problem, which necessitates the use of heuristics for tree-search. Although tree-search plays a major role in obtaining a tree estimate, there remains a limited understanding of its characteristics and how the elements of the statistical inferential procedure interact with the algorithms used. This study begins to answer some of these questions through a detailed examination of maximum likelihood tree-search on a wide range of real genome-scale data sets. We examine all 10,395 trees for each of the 106 genes of an eight-taxa yeast phylogenomic data set, then apply different tree-search algorithms to investigate their performance. We extend our findings by examining two larger genome-scale data sets and a large disparate data set that has been previously used to benchmark the performance of tree-search programs. We identify several broad trends occurring during tree-search that provide an insight into the performance of heuristics and may, in the future, aid their development. These trends include a tendency for the true maximum likelihood (best) tree to also be the shortest tree in terms of branch lengths, a weak tendency for tree-search to recover the best tree, and a tendency for tree-search to encounter fewer local optima in genes that have a high information content. When examining current heuristics for tree-search, we find that nearest-neighbor-interchange performs poorly, and frequently finds trees that are significantly different from the best tree. In contrast, subtree-pruning-and-regrafting tends to perform well, nearly always finding trees that are not significantly different to the best tree. Finally, we demonstrate that the precise implementation of a tree-search strategy, including when and where parameters are optimized, can change the character of tree-search, and that good strategies for tree-search may combine existing tree-search programs.

Study of Sylvilagus rabbit TRIM5α species-specific domain: how ancient endoviruses could have shaped the antiviral repertoire in Lagomorpha

Background

Since the first report of the antiretroviral restriction factor TRIM5α in primates, several orthologs in other mammals have been described. Recent studies suggest that leporid retroviruses like RELIK, the first reported endogenous lentivirus ever, may have imposed positive selection in TRIM5α orthologs of the European rabbit and European brown hare. Considering that RELIK must already have been present in a common ancestor of the leporid genera Lepus, Sylvilagus and Oryctolagus, we extended the study of evolutionary patterns of TRIM5α to other members of the Leporidae family, particularly to the genus Sylvilagus. Therefore, we obtained the TRIM5α nucleotide sequences of additional subspecies and species of the three leporid genera. We also compared lagomorph TRIM5α deduced protein sequences and established TRIM5α gene and TRIM5α protein phylogenies.

Results

The deduced protein sequence of Iberian hare TRIM5α was 89% identical to European rabbit TRIM5α, although high divergence was observed at the PRYSPRY v1 region between rabbit and the identified alleles from this hare species (allele 1: 50% divergence; allele 2: 53% divergence). A high identity was expected between the Sylvilagus and Oryctolagus TRIM5α proteins and, in fact, the Sylvilagus TRIM5α was 91% identical to the Oryctolagus protein. Nevertheless, the PRYSPRY v1 region was only 50% similar between these genera. Selection analysis of Lagomorpha TRIM5α proteins identified 25 positively-selected codons, 11 of which are located in the PRYSPRY v1 region, responsible for species specific differences in viral capsid recognition.

Conclusions

By extending Lagomorpha TRIM5α studies to an additional genus known to bear RELIK, we verified that the divergent species-specific pattern observed between the Oryctolagus and Lepus PRYSPRY-domains is also present in Sylvilagus TRIM5α. This work is one of the first known studies that compare the evolution of the antiretroviral restriction factor TRIM5α in different mammalian groups, Lagomorpha and Primates.

Substrate specificity of haloalkane dehalogenases

An enzyme's substrate specificity is one of its most important characteristics. The quantitative comparison of broad-specificity enzymes requires the selection of a homogenous set of substrates for experimental testing, determination of substrate-specificity data and analysis using multivariate statistics. We describe a systematic analysis of the substrate specificities of nine wild-type and four engineered haloalkane dehalogenases. The enzymes were characterized experimentally using a set of 30 substrates selected using statistical experimental design from a set of nearly 200 halogenated compounds. Analysis of the activity data showed that the most universally useful substrates in the assessment of haloalkane dehalogenase activity are 1-bromobutane, 1-iodopropane, 1-iodobutane, 1,2-dibromoethane and 4-bromobutanenitrile. Functional relationships among the enzymes were explored using principal component analysis. Analysis of the untransformed specific activity data revealed that the overall activity of wild-type haloalkane dehalogenases decreases in the following order: LinB~DbjA>DhlA~DhaA~DbeA~DmbA>DatA~DmbC~DrbA. After transforming the data, we were able to classify haloalkane dehalogenases into four SSGs (substrate-specificity groups). These functional groups are clearly distinct from the evolutionary subfamilies, suggesting that phylogenetic analysis cannot be used to predict the substrate specificity of individual haloalkane dehalogenases. Structural and functional comparisons of wild-type and mutant enzymes revealed that the architecture of the active site and the main access tunnel significantly influences the substrate specificity of these enzymes, but is not its only determinant. The identification of other structural determinants of the substrate specificity remains a challenge for further research on haloalkane dehalogenases.

Four gonadotropin releasing hormone receptor genes in Atlantic cod are differentially expressed in the brain and pituitary during puberty

Gonadotropin releasing hormones (GnRH) are an important part of the brain-pituitary-gonad axis in vertebrates. GnRH binding to its receptors (GnRH-R) stimulates synthesis and release of gonadotropins in the pituitary. GnRH-Rs also mediate other processes in the central nervous system such as reproductive behavior and neuromodulation. As many as five GnRH-R genes have been identified in two teleost fish species, but the function and phylogenetic relationship of these receptors is not fully understood. To gain a better understanding of the functional relationship between multiple GnRH-Rs in an important aquaculture species, the Atlantic cod (Gadus morhua), we identified four GnRH-Rs (gmGnRH-R) by RT-PCR, followed by full-length cloning and sequencing. The deduced amino acid sequences were used for phylogenetic analysis to identify conserved functional motifs and to clarify the relationship of gmGnRH-Rs with other vertebrate GnRH-Rs. The function of GnRH-R variants was investigated by quantitative PCR gene expression analysis in the brain and pituitary of female cod during a full reproductive cycle and in various peripheral tissues in sexually mature fish. Phylogenetic analysis revealed two types of teleost GnRH-Rs: Type I including gmGnRH-R1b and Type II including gmGnRH-R2a, gmGnRH-R2b and gmGnRH-R2c. All four gmGnRH-Rs are expressed in the brain, and gmGnRH-R1b, gmGnRH-R2a and gmGnRH-R2c are expressed in the pituitary. The only GnRH-R differentially expressed in the pituitary during the reproductive cycle is gmGnRH-R2a such that its expression is significantly increased during spawning. These data suggest that gmGnRH-R2a is the most likely candidate to mediate the hypophysiotropic function of GnRH in Atlantic cod.

An efficient algorithm for approximating geodesic distances in tree space

The increasing use of phylogeny in biological studies is limited by the need to make available more efficient tools for computing distances between trees. The geodesic tree distance-introduced by Billera, Holmes, and Vogtmann-combines both the tree topology and edge lengths into a single metric. Despite the conceptual simplicity of the geodesic tree distance, algorithms to compute it don't scale well to large, real-world phylogenetic trees composed of hundred or even thousand leaves. In this paper, we propose the geodesic distance as an effective tool for exploring the likelihood profile in the space of phylogenetic trees, and we give a cubic time algorithm, GeoHeuristic, in order to compute an approximation of the distance. We compare it with the GTP algorithm, which calculates the exact distance, and the cone path length, which is another approximation, showing that GeoHeuristic achieves a quite good trade-off between accuracy (relative error always lower than 0.0001) and efficiency. We also prove the equivalence among GeoHeuristic, cone path, and Robinson-Foulds distances when assuming branch lengths equal to unity and we show empirically that, under this restriction, these distances are almost always equal to the actual geodesic.

Diversity and plasticity of the intracellular plant pathogen and insect symbiont "Candidatus Liberibacter asiaticus" as revealed by hypervariable prophage genes with intragenic tandem repeats

"Candidatus Liberibacter asiaticus" is a psyllid-transmitted, phloem-limited alphaproteobacterium and the most prevalent species of "Ca. Liberibacter" associated with a devastating worldwide citrus disease known as huanglongbing (HLB). Two related and hypervariable genes (hyv(I) and hyv(II)) were identified in the prophage regions of the Psy62 "Ca. Liberibacter asiaticus" genome. Sequence analyses of the hyv(I) and hyv(II) genes in 35 "Ca. Liberibacter asiaticus" DNA isolates collected globally revealed that the hyv(I) gene contains up to 12 nearly identical tandem repeats (NITRs, 132 bp) and 4 partial repeats, while hyv(II) contains up to 2 NITRs and 4 partial repeats and shares homology with hyv(I). Frequent deletions or insertions of these repeats within the hyv(I) and hyv(II) genes were observed, none of which disrupted the open reading frames. Sequence conservation within the individual repeats but an extensive variation in repeat numbers, rearrangement, and the sequences flanking the repeat region indicate the diversity and plasticity of "Ca. Liberibacter asiaticus" bacterial populations in the world. These differences were found not only in samples of distinct geographical origins but also in samples from a single origin and even from a single "Ca. Liberibacter asiaticus"-infected sample. This is the first evidence of different "Ca. Liberibacter asiaticus" populations coexisting in a single HLB-affected sample. The Florida "Ca. Liberibacter asiaticus" isolates contain both hyv(I) and hyv(II), while all other global "Ca. Liberibacter asiaticus" isolates contain either one or the other. Interclade assignments of the putative Hyv(I) and Hyv(II) proteins from Florida isolates with other global isolates in phylogenetic trees imply multiple "Ca. Liberibacter asiaticus" populations in the world and a multisource introduction of the "Ca. Liberibacter asiaticus" bacterium into Florida.

Nephele: genotyping via complete composition vectors and MapReduce

BACKGROUND

Current sequencing technology makes it practical to sequence many samples of a given organism, raising new challenges for the processing and interpretation of large genomics data sets with associated metadata. Traditional computational phylogenetic methods are ideal for studying the evolution of gene/protein families and using those to infer the evolution of an organism, but are less than ideal for the study of the whole organism mainly due to the presence of insertions/deletions/rearrangements. These methods provide the researcher with the ability to group a set of samples into distinct genotypic groups based on sequence similarity, which can then be associated with metadata, such as host information, pathogenicity, and time or location of occurrence. Genotyping is critical to understanding, at a genomic level, the origin and spread of infectious diseases. Increasingly, genotyping is coming into use for disease surveillance activities, as well as for microbial forensics. The classic genotyping approach has been based on phylogenetic analysis, starting with a multiple sequence alignment. Genotypes are then established by expert examination of phylogenetic trees. However, these traditional single-processor methods are suboptimal for rapidly growing sequence datasets being generated by next-generation DNA sequencing machines, because they increase in computational complexity quickly with the number of sequences.

RESULTS

Nephele is a suite of tools that uses the complete composition vector algorithm to represent each sequence in the dataset as a vector derived from its constituent k-mers by passing the need for multiple sequence alignment, and affinity propagation clustering to group the sequences into genotypes based on a distance measure over the vectors. Our methods produce results that correlate well with expert-defined clades or genotypes, at a fraction of the computational cost of traditional phylogenetic methods run on traditional hardware. Nephele can use the open-source Hadoop implementation of MapReduce to parallelize execution using multiple compute nodes. We were able to generate a neighbour-joined tree of over 10,000 16S samples in less than 2 hours.

CONCLUSIONS

We conclude that using Nephele can substantially decrease the processing time required for generating genotype trees of tens to hundreds of organisms at genome scale sequence coverage.

Identification and gene expression analysis of three GnRH genes in female Atlantic cod during puberty provides insight into GnRH variant gene loss in fish

Gonadotropin releasing hormone (GnRH) is a key regulator of sexual development and reproduction in vertebrates. Fish have either two or three pre-pro-GnRH genes, encoding structurally distinct peptides. We identified three pre-pro-GnRH genes in Atlantic cod (Gadus morhua, gmGnRH) using RT-PCR, RACE-PCR and BAC DNA library clone sequencing based on synteny searching. Gene identity was confirmed by sequence alignment and subsequent phylogenetic analysis. The expression of these genes was measured by quantitative PCR in the brain and pituitary of female cod throughout their reproductive cycle and in peripheral tissues. All three gmGnRH genes have highly conserved deduced decapeptide sequences, but sequence and phylogenetic data for gmGnRH1 suggest that this is a pseudogene. gmGnRH1 shares low identity with all fish GnRH variants and grouped with the GnRH3 clade. Although gmGnRH1 is a putative pseudogene, it is transcribed in multiple tissues but at low levels in the brain, indicating the loss of conserved hypophysiotrophic function. Phylogenetic analysis reveals that gmGnRH2 and gmGnRH3 variants are located in variant-specific clades. Both gmGnRH2 and gmGnRH3 transcripts are most abundant in the brain, with lower expression in pituitaries and ovaries. Brain gmGnRH3 gene expression increases in spawning fish and is expressed in the pituitary during puberty. Brain gmGnRH2 transcripts are highly expressed relative to gmGnRH3 before and during spawning. Sequence and expression data suggest that gmGnRH1 is a pseudogene and that gmGnRH3 is likely the hypophysiotrophic form of GnRH in Atlantic cod.

The first molecular phylogeny of Strepsiptera (Insecta) reveals an early burst of molecular evolution correlated with the transition to endoparasitism

A comprehensive model of evolution requires an understanding of the relationship between selection at the molecular and phenotypic level. We investigate this in Strepsiptera, an order of endoparasitic insects whose evolutionary biology is poorly studied. We present the first molecular phylogeny of Strepsiptera, and use this as a framework to investigate the association between parasitism and molecular evolution. We find evidence of a significant burst in the rate of molecular evolution in the early history of Strepsiptera. The evolution of morphological traits linked to parasitism is significantly correlated with the pattern in molecular rate. The correlated burst in genotypic-phenotypic evolution precedes the main phase of strepsipteran diversification, which is characterised by the return to a low and even molecular rate, and a period of relative morphological stability. These findings suggest that the transition to endoparasitism led to relaxation of selective constraint in the strepsipteran genome. Our results indicate that a parasitic lifestyle can affect the rate of molecular evolution, although other causal life-history traits correlated with parasitism may also play an important role.

Contrasting 5' and 3' evolutionary histories and frequent evolutionary convergence in Meis/hth gene structures

Organisms show striking differences in genome structure; however, the functional implications and fundamental forces that govern these differences remain obscure. The intron–exon organization of nuclear genes is involved in a particularly large variety of structures and functional roles. We performed a 22-species study of Meis/hth genes, intron-rich homeodomain-containing transcription factors involved in a wide range of developmental processes. Our study revealed three surprising results that suggest important and very different functions for Meis intron–exon structures. First, we find unexpected conservation across species of intron positions and lengths along most of the Meis locus. This contrasts with the high degree of structural divergence found in genome-wide studies and may attest to conserved regulatory elements residing within these conserved introns. Second, we find very different evolutionary histories for the 5′ and 3′ regions of the gene. The 5′-most 10 exons, which encode the highly conserved Meis domain and homeodomain, show striking conservation. By contrast, the 3′ of the gene, which encodes several domains implicated in transcriptional activation and response to cell signaling, shows a remarkably active evolutionary history, with diverse isoforms and frequent creation and loss of new exons and splice sites. This region-specific diversity suggests evolutionary “tinkering,” with alternative splicing allowing for more subtle regulation of protein function. Third, we find a large number of cases of convergent evolution in the 3′ region, including 1) parallel losses of ancestral coding sequence, 2) parallel gains of external and internal splice sites, and 3) recurrent truncation of C-terminal coding regions. These results attest to the importance of locus-specific splicing functions in differences in structural evolution across genes, as well as to commonalities of forces shaping the evolution of individual genes along different lineages.

SuiteMSA: visual tools for multiple sequence alignment comparison and molecular sequence simulation

Background

Multiple sequence alignment (MSA) plays a central role in nearly all bioinformatics and molecular evolutionary applications. MSA reconstruction is thus one of the most heavily scrutinized bioinformatics fields. Evaluating the quality of MSA reconstruction is often hindered by the lack of good reference MSAs. The use of sequence evolution simulation can provide such reference MSAs. Furthermore, none of the MSA viewing/editing programs currently available allows the user to make direct comparisons between two or more MSAs. Considering the importance of MSA quality in a wide range of research, it is desirable if MSA assessment can be performed more easily.

Results

We have developed SuiteMSA, a java-based application that provides unique MSA viewers. Users can directly compare multiple MSAs and evaluate where the MSAs agree (are consistent) or disagree (are inconsistent). Several alignment statistics are provided to assist such comparisons. SuiteMSA also includes a graphical phylogeny editor/viewer as well as a graphical user interface for a sequence evolution simulator that can be used to construct reference MSAs.

Conclusions

SuiteMSA provides researchers easy access to a sequence evolution simulator, reference alignments generated by the simulator, and a series of tools to evaluate the performance of the MSA reconstruction programs. It will help us improve the quality of MSAs, often the most important first steps of bioinformatics and other biological research.

Widespread occurrence of secondary lipid biosynthesis potential in microbial lineages

Bacterial production of long-chain omega-3 polyunsaturated fatty acids (PUFAs), such as eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3), is constrained to a narrow subset of marine γ-proteobacteria. The genes responsible for de novo bacterial PUFA biosynthesis, designated pfaEABCD, encode large, multi-domain protein complexes akin to type I iterative fatty acid and polyketide synthases, herein referred to as "Pfa synthases". In addition to the archetypal Pfa synthase gene products from marine bacteria, we have identified homologous type I FAS/PKS gene clusters in diverse microbial lineages spanning 45 genera representing 10 phyla, presumed to be involved in long-chain fatty acid biosynthesis. In total, 20 distinct types of gene clusters were identified. Collectively, we propose the designation of "secondary lipids" to describe these biosynthetic pathways and products, a proposition consistent with the "secondary metabolite" vernacular. Phylogenomic analysis reveals a high degree of functional conservation within distinct biosynthetic pathways. Incongruence between secondary lipid synthase functional clades and taxonomic group membership combined with the lack of orthologous gene clusters in closely related strains suggests horizontal gene transfer has contributed to the dissemination of specialized lipid biosynthetic activities across disparate microbial lineages.

ProtTest 3: fast selection of best-fit models of protein evolution

We have implemented a high-performance computing (HPC) version of ProtTest that can be executed in parallel in multicore desktops and clusters. This version, called ProtTest 3, includes new features and extended capabilities.

AVAILABILITY

ProtTest 3 source code and binaries are freely available under GNU license for download from http://darwin.uvigo.es/software/prottest3, linked to a Mercurial repository at Bitbucket (https://bitbucket.org/).

CONTACT

dposada@uvigo.es

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Stepwise assembly of the Nova-regulated alternative splicing network in the vertebrate brain

Novel organismal structures in metazoans are often undergirded by complex gene regulatory networks; as such, understanding the emergence of new structures through evolution requires reconstructing the series of evolutionary steps leading to these underlying networks. Here, we reconstruct the step-by-step assembly of the vertebrate splicing network regulated by Nova, a splicing factor that modulates alternative splicing in the vertebrate central nervous system by binding to clusters of YCAY motifs on pre-RNA transcripts. Transfection of human HEK293T cells with Nova orthologs indicated vertebrate-like splicing regulatory activity in bilaterian invertebrates, thus Nova acquired the ability to bind YCAY clusters and perform vertebrate-like splicing modulation at least before the last common ancestor of bilaterians. In situ hybridization studies in several species showed that Nova expression became restricted to CNS later on, during chordate evolution. Finally, comparative genomics studies revealed a diverse history for Nova-regulated exons, with target exons arising through both de novo exon creation and acquisition of YCAY motifs by preexisting exons throughout chordate and vertebrate history. In addition, we find that tissue-specific Nova expression patterns emerged independently in other lineages, suggesting independent assembly of tissue-specific regulatory networks.

The putative AKH receptor of the tobacco hornworm, Manduca sexta, and its expression

Adipokinetic hormones are peptide hormones that mobilize lipids and/or carbohydrates for flight in adult insects and activate glycogen Phosphorylase in larvae during starvation and during molt. We previously examined the functional roles of adipokinetic hormone in Manduca sexta L. (Lepidoptera: Sphingidae). Here we report the cloning of the full-length cDNA encoding the putative adipokinetic hormone receptor from the fat body of M. sexta. The sequence analysis shows that the deduced amino acid sequence shares common motifs of G protein-coupled receptors, by having seven hydrophobic transmembrane segments. We examined the mRNA expression pattern of the adipokinetic hormone receptor by quantitative Real-Time PCR in fat body during development and in different tissues and found the strongest expression in fat body of larvae two days after molt to the fifth instar. We discuss these results in relation to some of our earlier results. We also compare the M. sexta adipokinetic hormone receptor with the known adipokinetic hormone receptors of other insects and with gonadotropin releasing hormone-like receptors of invertebrates.

Diversity and distribution of single-stranded DNA phages in the North Atlantic Ocean

Knowledge of marine phages is highly biased toward double-stranded DNA (dsDNA) phages; however, recent metagenomic surveys have also identified single-stranded DNA (ssDNA) phages in the oceans. Here, we describe two complete ssDNA phage genomes that were reconstructed from a viral metagenome from 80 m depth at the Bermuda Atlantic Time-series Study (BATS) site in the northwestern Sargasso Sea and examine their spatial and temporal distributions. Both genomes (SARssφ1 and SARssφ2) exhibited similarity to known phages of the Microviridae family in terms of size, GC content, genome organization and protein sequence. PCR amplification of the replication initiation protein (Rep) gene revealed narrow and distinct depth distributions for the newly described ssDNA phages within the upper 200 m of the water column at the BATS site. Comparison of Rep gene sequences obtained from the BATS site over time revealed changes in the diversity of ssDNA phages over monthly time scales, although some nearly identical sequences were recovered from samples collected 4 years apart. Examination of ssDNA phage diversity along transects through the North Atlantic Ocean revealed a positive correlation between genetic distance and geographic distance between sampling sites. Together, the data suggest fundamental differences between the distribution of these ssDNA phages and the distribution of known marine dsDNA phages, possibly because of differences in host range, host distribution, virion stability, or viral evolution mechanisms and rates. Future work needs to elucidate the host ranges for oceanic ssDNA phages and determine their ecological roles in the marine ecosystem.

Diversity and distribution of microbial long-chain fatty acid biosynthetic genes in the marine environment

Bacterial production of long-chain fatty acids via a polyketide synthase-related mechanism has thus far only been investigated in isolate-based studies. Here, the genetic capacity for production of long-chain fatty acids was investigated using a culture-independent approach. PCR primers targeting the keto-acyl synthase (KS) domain of the pfaA gene involved in omega-3 polyunsaturated fatty acid (PUFA) biosynthesis were used to construct clone libraries to investigate KS sequence diversity in disparate marine habitats. Of the 446 sequences recovered, 123 (27.6%) clustered with KS sequences involved in the synthesis of eicosapentaenoic acid (EPA, C20:5n-3), docosahexaenoic acid (DHA, C22:6n-3) and arachidonic acid (AA, C20:4n-6). The remaining 72.4% of clones formed environmental-only groups or grouped with the KS domains of pfaA homologues from organisms producing unidentified products. In total, 17 groups were recovered - four known and 13 newly identified. A query of metagenomic data sets revealed sequences related to EPA KS domains, as well as sequences related to four environmental-only groups discovered in the clone libraries. The phylogenetic affiliation and end product of these environmental-only KS clusters is unknown. These findings reveal a widespread capacity for long-chain fatty acid production in marine microorganisms, including biosynthetic pathways not yet characterized.

Concept and application of a computational vaccinology workflow

BACKGROUND

The last years have seen a renaissance of the vaccine area, driven by clinical needs in infectious diseases but also chronic diseases such as cancer and autoimmune disorders. Equally important are technological improvements involving nano-scale delivery platforms as well as third generation adjuvants. In parallel immunoinformatics routines have reached essential maturity for supporting central aspects in vaccinology going beyond prediction of antigenic determinants. On this basis computational vaccinology has emerged as a discipline aimed at ab-initio rational vaccine design.Here we present a computational workflow for implementing computational vaccinology covering aspects from vaccine target identification to functional characterization and epitope selection supported by a Systems Biology assessment of central aspects in host-pathogen interaction. We exemplify the procedures for Epstein Barr Virus (EBV), a clinically relevant pathogen causing chronic infection and suspected of triggering malignancies and autoimmune disorders.

RESULTS

We introduce pBone/pView as a computational workflow supporting design and execution of immunoinformatics workflow modules, additionally involving aspects of results visualization, knowledge sharing and re-use. Specific elements of the workflow involve identification of vaccine targets in the realm of a Systems Biology assessment of host-pathogen interaction for identifying functionally relevant targets, as well as various methodologies for delineating B- and T-cell epitopes with particular emphasis on broad coverage of viral isolates as well as MHC alleles.Applying the workflow on EBV specifically proposes sequences from the viral proteins LMP2, EBNA2 and BALF4 as vaccine targets holding specific B- and T-cell epitopes promising broad strain and allele coverage.

CONCLUSION

Based on advancements in the experimental assessment of genomes, transcriptomes and proteomes for both, pathogen and (human) host, the fundaments for rational design of vaccines have been laid out. In parallel, immunoinformatics modules have been designed and successfully applied for supporting specific aspects in vaccine design. Joining these advancements, further complemented by novel vaccine formulation and delivery aspects, have paved the way for implementing computational vaccinology for rational vaccine design tackling presently unmet vaccine challenges.

Interspersed DNA repeats bcr1-bcr18 of Bacillus cereus group bacteria form three distinct groups with different evolutionary and functional patterns

Many short (<400 bp) interspersed sequence repeats exist in bacteria, yet little is known about their origins, mode of generation, or possible function. Here, we present a comprehensive analysis of 18 different previously identified repeated DNA elements, bcr1-bcr18 (Økstad OA, Hegna I, Lindback T, Rishovd AL, Kolstø AB. 1999. Genome organization is not conserved between Bacillus cereus and Bacillus subtilis. Microbiology. 145:621-631.; Tourasse NJ, Helgason E, Økstad OA, Hegna IK, Kolstø AB. 2006. The Bacillus cereus group: novel aspects of population structure and genome dynamics. J Appl Microbiol. 101:579-593.), in 36 sequenced genomes from the Bacillus cereus group of bacteria. This group consists of genetically closely related species with variable pathogenic specificity toward different hosts and includes among others B. anthracis, B. cereus, and B. thuringiensis. The B. cereus group repeat elements could be classified into three categories with different properties: Group A elements (bcr1-bcr3) exhibited highly variable copy numbers ranging from 4 to 116 copies per strain, showed a nonconserved chromosomal distribution pattern between strains, and displayed several features characteristic of mobile elements. Group B repeats (bcr4-bcr6) were present in 0-10 copies per strain and were associated with strain-specific genes and disruptions of genome synteny, implying a possible contribution to genome rearrangements and/or horizontal gene transfer events. bcr5, in particular, was associated with large gene clusters showing resemblance to integrons. In agreement with their potentially mobile nature or involvement in horizontal transfers, the sequences of the repeats from Groups A and B (bcr1-bcr6) followed a phylogeny different from that of the host strains. Conversely, repeats from Group C (bcr7-bcr18) had a conserved chromosomal location and orthologous gene neighbors in the investigated B. cereus group genomes, and their phylogeny matched that of the host chromosome. Several of the group C repeats exhibited a conserved secondary structure or had parts of the structure conserved, possibly indicating functional RNAs. Accordingly, five of the repeats in group C overlapped regions encoding previously characterized riboswitches. Similarly, other group C repeats could represent novel riboswitches, encode small RNAs, and/or constitute other types of regulatory elements with specific biological functions. The current analysis suggests that the multitude of repeat elements identified in the B. cereus group promote genome dynamics and plasticity and could contribute to the flexible and adaptive life style of these bacteria.

Conserved developmental expression of Fezf in chordates and Drosophila and the origin of the Zona Limitans Intrathalamica (ZLI) brain organizer

Background

The zona limitans intrathalamica (ZLI) and the isthmus organizer (IsO) are two major secondary organizers of vertebrate brain development. These organizers are located at the interface of the expression domains of key patterning genes (Fezf-Irx and Otx-Gbx, respectively). To gain insights into the evolutionary origin of the ZLI, we studied Fezf in bilaterians.

Results

In this paper, we identified a conserved sequence motif (Fezf box) in all bilaterians. We report the expression pattern of Fezf in amphioxus and Drosophila and compare it with those of Gbx, Otx and Irx. We found that the relative expression patterns of these genes in vertebrates are fully conserved in amphioxus and flies, indicating that the genetic subdivisions defining the location of both secondary organizers in early vertebrate brain development were probably present in the last common ancestor of extant bilaterians. However, in contrast to vertebrates, we found that Irx-defective flies do not show an affected Fezf expression pattern.

Conclusions

The absence of expression of the corresponding morphogens from cells at these conserved genetic boundaries in invertebrates suggests that the organizing properties might have evolved specifically in the vertebrate lineage by the recruitment of key morphogens to these conserved genetic locations.