Adding unaligned sequences into an existing alignment using MAFFT and LAST

Environment and host identity structure communities of green algal symbionts in lichens

An understanding of how biotic interactions shape species' distributions is central to predicting host-symbiont responses under climate change. Switches to locally adapted algae have been proposed to be an adaptive strategy of lichen-forming fungi to cope with environmental change. However, it is unclear how lichen photobionts respond to environmental gradients, and whether they play a role in determining the fungal host's upper and lower elevational limits. Deep-coverage Illumina DNA metabarcoding was used to track changes in the community composition of Trebouxia algae associated with two phylogenetically closely related, but ecologically divergent fungal hosts along a steep altitudinal gradient in the Mediterranean region. We detected the presence of multiple Trebouxia species in the majority of thalli. Both altitude and host genetic identity were strong predictors of photobiont community assembly in these two species. The predominantly clonally dispersing fungus showed stronger altitudinal structuring of photobiont communities than the sexually reproducing host. Elevation ranges of the host were not limited by the lack of compatible photobionts. Our study sheds light on the processes guiding the formation and distribution of specific fungal-algal combinations in the lichen symbiosis. The effect of environmental filtering acting on both symbiotic partners appears to shape the distribution of lichens.

Divergent HIV-1 strains (CRF92_C2U and CRF93_cpx) co-circulating in the Democratic Republic of the Congo: Phylogenetic insights on the early evolutionary history of subtype C

Molecular epidemiological studies revealed that the epicenter of the HIV pandemic was Kinshasa, the capital city of the Democratic Republic of the Congo (DRC) in Central Africa. All known subtypes and numerous complex recombinant strains co-circulate in the DRC. Moreover, high intra-subtype diversity has been also documented. During two previous surveys on HIV-1 antiretroviral drug resistance in the DRC, we identified two divergent subtype C lineages in the protease and partial reverse transcriptase gene regions. We sequenced eight near full-length genomes and classified them using bootscanning and likelihood-based phylogenetic analyses. Four strains are more closely related to subtype C although within the range of inter sub-subtype distances. However, these strains also have small unclassified fragments and thus were named CRF92_C2U. Another strain is a unique recombinant of CRF92_C2U with an additional small unclassified fragment and a small divergent subtype A fragment. The three remaining strains represent a complex mosaic named CRF93_cpx. CRF93_cpx have two fragments of divergent subtype C sequences, which are not conventional subtype C nor the above described C2, and multiple divergent subtype A-like fragments. We then inferred the time-scaled evolutionary history of subtype C following a Bayesian approach and a partitioned analysis using major genomic regions. CRF92_C2U and CRF93_cpx had the most recent common ancestor with conventional subtype C around 1932 and 1928, respectively. A Bayesian demographic reconstruction corroborated that the subtype C transition to a faster phase of exponential growth occurred during the 1950s. Our analysis showed considerable differences between the newly discovered early-divergent strains and the conventional subtype C and therefore suggested that this virus has been diverging in humans for several decades before the HIV/M diversity boom in the 1950s.

A comparison of chloroplast genome sequences in Aconitum (Ranunculaceae): a traditional herbal medicinal genus

The herbal medicinal genus Aconitum L., belonging to the Ranunculaceae family, represents the earliest diverging lineage within the eudicots. It currently comprises of two subgenera, A. subgenus Lycoctonum and A. subg. Aconitum. The complete chloroplast (cp) genome sequences were characterized in three species: A. angustius, A. finetianum, and A. sinomontanum in subg. Lycoctonum and compared to other Aconitum species to clarify their phylogenetic relationship and provide molecular information for utilization of Aconitum species particularly in Eastern Asia. The length of the chloroplast genome sequences were 156,109 bp in A. angustius, 155,625 bp in A. finetianum and 157,215 bp in A. sinomontanum, with each species possessing 126 genes with 84 protein coding genes (PCGs). While genomic rearrangements were absent, structural variation was detected in the LSC/IR/SSC boundaries. Five pseudogenes were identified, among which Ψrps19 and Ψycf1 were in the LSC/IR/SSC boundaries, Ψrps16 and ΨinfA in the LSC region, and Ψycf15 in the IRb region. The nucleotide variability (Pi) of Aconitum was estimated to be 0.00549, with comparably higher variations in the LSC and SSC than the IR regions. Eight intergenic regions were revealed to be highly variable and a total of 58-62 simple sequence repeats (SSRs) were detected in all three species. More than 80% of SSRs were present in the LSC region. Altogether, 64.41% and 46.81% of SSRs are mononucleotides in subg. Lycoctonum and subg. Aconitum, respectively, while a higher percentage of di-, tri-, tetra-, and penta- SSRs were present in subg. Aconitum. Most species of subg. Aconitum in Eastern Asia were first used for phylogenetic analyses. The availability of the complete cp genome sequences of these species in subg. Lycoctonum will benefit future phylogenetic analyses and aid in germplasm utilization in Aconitum species.

Evolutionary and molecular foundations of multiple contemporary functions of the nitroreductase superfamily

Insight regarding how diverse enzymatic functions and reactions have evolved from ancestral scaffolds is fundamental to understanding chemical and evolutionary biology, and for the exploitation of enzymes for biotechnology. We undertook an extensive computational analysis using a unique and comprehensive combination of tools that include large-scale phylogenetic reconstruction to determine the sequence, structural, and functional relationships of the functionally diverse flavin mononucleotide-dependent nitroreductase (NTR) superfamily (>24,000 sequences from all domains of life, 54 structures, and >10 enzymatic functions). Our results suggest an evolutionary model in which contemporary subgroups of the superfamily have diverged in a radial manner from a minimal flavin-binding scaffold. We identified the structural design principle for this divergence: Insertions at key positions in the minimal scaffold that, combined with the fixation of key residues, have led to functional specialization. These results will aid future efforts to delineate the emergence of functional diversity in enzyme superfamilies, provide clues for functional inference for superfamily members of unknown function, and facilitate rational redesign of the NTR scaffold.

First Molecular Detection of Babesia gibsoni in Dogs from Wuhan, China

Canine piroplasmosis is a significant disease in dogs caused by Babesia and Theileria parasites. The clinical manifestations range from mild illness to serious disease depending on the parasite species and the physical condition of the infected dog. Canine piroplasmosis has been reported to be prevalent in China. However, no molecular evidence of the disease has been reported in pet dogs from Wuhan. In this study, 118 blood samples were randomly collected from pet dogs in veterinary clinics. The blood samples were subjected to both microscopic examination and reverse line blot (RLB) hybridization assays to detect piroplasm infection. Parasites were observed in 10 blood samples via microscopic examination, whereas there were 14 Babesia gibsoni-positive RLB tests. Phylogenetic analysis was performed after the 18S rRNA and ITS gene sequences from the 14 positive samples were cloned and sequenced. The results confirmed the existence of B. gibsoni in this area. This is the first molecular report of canine babesiosis in pet dogs from Wuhan, China. Pet dogs are companion animals, and the prevalence of babesiosis will be of concern in daily life. This study will help veterinarians better understand the prevalence of canine babesiosis and provide a guide for disease control in pet dogs.

Multi-gene phylogenetic analysis reveals that shochu-fermenting Saccharomyces cerevisiae strains form a distinct sub-clade of the Japanese sake cluster

Shochu is a traditional Japanese distilled spirit. The formation of the distinguishing flavour of shochu produced in individual distilleries is attributed to putative indigenous yeast strains. In this study, we performed the first (to our knowledge) phylogenetic classification of shochu strains based on nucleotide gene sequences. We performed phylogenetic classification of 21 putative indigenous shochu yeast strains isolated from 11 distilleries. All of these strains were shown or confirmed to be Saccharomyces cerevisiae, sharing species identification with 34 known S. cerevisiae strains (including commonly used shochu, sake, ale, whisky, bakery, bioethanol and laboratory yeast strains and clinical isolate) that were tested in parallel. Our analysis used five genes that reflect genome-level phylogeny for the strain-level classification. In a first step, we demonstrated that partial regions of the ZAP1, THI7, PXL1, YRR1 and GLG1 genes were sufficient to reproduce previous sub-species classifications. In a second step, these five analysed regions from each of 25 strains (four commonly used shochu strains and the 21 putative indigenous shochu strains) were concatenated and used to generate a phylogenetic tree. Further analysis revealed that the putative indigenous shochu yeast strains form a monophyletic group that includes both the shochu yeasts and a subset of the sake group strains; this cluster is a sister group to other sake yeast strains, together comprising a sake-shochu group. Differences among shochu strains were small, suggesting that it may be possible to correlate subtle phenotypic differences among shochu flavours with specific differences in genome sequences. Copyright © 2017 John Wiley & Sons, Ltd.

Pan-genome and phylogeny of Bacillus cereus sensu lato

BACKGROUND

Bacillus cereus sensu lato (s. l.) is an ecologically diverse bacterial group of medical and agricultural significance. In this study, I use publicly available genomes and novel bioinformatic workflows to characterize the B. cereus s. l. pan-genome and perform the largest phylogenetic and population genetic analyses of this group to date in terms of the number of genes and taxa included. With these fundamental data in hand, I identify genes associated with particular phenotypic traits (i.e., "pan-GWAS" analysis), and quantify the degree to which taxa sharing common attributes are phylogenetically clustered.

METHODS

A rapid k-mer based approach (Mash) was used to create reduced representations of selected Bacillus genomes, and a fast distance-based phylogenetic analysis of this data (FastME) was performed to determine which species should be included in B. cereus s. l. The complete genomes of eight B. cereus s. l. species were annotated de novo with Prokka, and these annotations were used by Roary to produce the B. cereus s. l. pan-genome. Scoary was used to associate gene presence and absence patterns with various phenotypes. The orthologous protein sequence clusters produced by Roary were filtered and used to build HaMStR databases of gene models that were used in turn to construct phylogenetic data matrices. Phylogenetic analyses used RAxML, DendroPy, ClonalFrameML, PAUP*, and SplitsTree. Bayesian model-based population genetic analysis assigned taxa to clusters using hierBAPS. The genealogical sorting index was used to quantify the phylogenetic clustering of taxa sharing common attributes.

RESULTS

The B. cereus s. l. pan-genome currently consists of ≈60,000 genes, ≈600 of which are "core" (common to at least 99% of taxa sampled). Pan-GWAS analysis revealed genes associated with phenotypes such as isolation source, oxygen requirement, and ability to cause diseases such as anthrax or food poisoning. Extensive phylogenetic analyses using an unprecedented amount of data produced phylogenies that were largely concordant with each other and with previous studies. Phylogenetic support as measured by bootstrap probabilities increased markedly when all suitable pan-genome data was included in phylogenetic analyses, as opposed to when only core genes were used. Bayesian population genetic analysis recommended subdividing the three major clades of B. cereus s. l. into nine clusters. Taxa sharing common traits and species designations exhibited varying degrees of phylogenetic clustering.

CONCLUSIONS

All phylogenetic analyses recapitulated two previously used classification systems, and taxa were consistently assigned to the same major clade and group. By including accessory genes from the pan-genome in the phylogenetic analyses, I produced an exceptionally well-supported phylogeny of 114 complete B. cereus s. l. genomes. The best-performing methods were used to produce a phylogeny of all 498 publicly available B. cereus s. l. genomes, which was in turn used to compare three different classification systems and to test the monophyly status of various B. cereus s. l. species. The majority of the methodology used in this study is generic and could be leveraged to produce pan-genome estimates and similarly robust phylogenetic hypotheses for other bacterial groups.

Diversity and oceanic distribution of prasinophytes clade VII, the dominant group of green algae in oceanic waters

Prasinophytes clade VII is a group of pico/nano-planktonic green algae (division Chlorophyta) for which numerous ribosomal RNA (rRNA) sequences have been retrieved from the marine environment in the last 15 years. A large number of strains have also been isolated but have not yet received a formal taxonomic description. A phylogenetic analysis of available strains using both the nuclear 18S and plastidial 16S rRNA genes demonstrates that this group composes at least 10 different clades: A1-A7 and B1-B3. Analysis of sequences from the variable V9 region of the 18S rRNA gene collected during the Tara Oceans expedition and in the frame of the Ocean Sampling Day consortium reveal that clade VII is the dominant Chlorophyta group in oceanic waters, replacing Mamiellophyceae, which have this role in coastal waters. At some location, prasinophytes clade VII can even be the dominant photosynthetic eukaryote representing up to 80% of photosynthetic metabarcodes overall. B1 and A4 are the overall dominant clades and different clades seem to occupy distinct niches, for example, A6 is dominant in surface Mediterranean Sea waters, whereas A4 extend to high temperate latitudes. Our work demonstrates that prasinophytes clade VII constitute a highly diversified group, which is a key component of phytoplankton in open oceanic waters but has been neglected in the conceptualization of marine microbial diversity and carbon cycle.

Oligotrophic lagoons of the South Pacific Ocean are home to a surprising number of novel eukaryotic microorganisms

The diversity of microbial eukaryotes was surveyed by environmental sequencing from tropical lagoon sites of the South Pacific, collected through the American Museum of Natural History (AMNH)'s Explore21 expedition to the Solomon Islands in September 2013. The sampled lagoons presented low nutrient concentrations typical of oligotrophic waters, but contained levels of chlorophyll a, a proxy for phytoplankton biomass, characteristic of meso- to eutrophic waters. Two 18S rDNA hypervariable sites, the V4 and V8-V9 regions, were amplified from the total of eight lagoon samples and sequenced on the MiSeq system. After assembly, clustering at 97% similarity, and removal of singletons and chimeras, a total of 2741 (V4) and 2606 (V8-V9) operational taxonomic units (OTUs) were identified. Taxonomic annotation of these reads, including phylogeny, was based on a combination of automated pipeline and manual inspection. About 18.4% (V4) and 13.8% (V8-V9) of the OTUs could not be assigned to any of the known eukaryotic groups. Of these, we focused on OTUs that were not divergent and possessed multiple sources of evidence for their existence. Phylogenetic analyses of these sequences revealed more than ten branches that might represent new deeply-branching lineages of microbial eukaryotes, currently without any cultured representatives or morphological information.

Identification of Complete Repertoire of Apis florea Odorant Receptors Reveals Complex Orthologous Relationships with Apis mellifera

We developed a computational pipeline for homology based identification of the complete repertoire of olfactory receptor (OR) genes in the Asian honey bee species, Apis florea. Apis florea is phylogenetically the most basal honey bee species and also the most distant sister species to the Western honey bee Apis mellifera, for which all OR genes had been identified before. Using our pipeline, we identified 180 OR genes in A. florea, which is very similar to the number of ORs identified in A. mellifera (177 ORs). Many characteristics of the ORs including gene structure, synteny of tandemly repeated ORs and basic phylogenetic clustering are highly conserved. The composite phylogenetic tree of A. florea and A. mellifera ORs could be divided into 21 clades which are in harmony with the existing Hymenopteran tree. However, we found a few nonorthologous OR relationships between both species as well as independent pseudogenization of ORs suggesting separate evolutionary changes. Particularly, a subgroup of the OR gene clade XI, which had been hypothesized to code cuticular hydrocarbon receptors showed a high number of species-specific ORs. RNAseq analysis detected a total number of 145 OR transcripts in male and 162 in female antennae. Most of the OR genes were highly expressed on the female antennae. However, we detected five distinct male-biased OR genes, out of which three genes (AfOr11, AfOr18, AfOr170P) were shown to be male-biased in A. mellifera, too, thus corroborating a behavioral function in sex-pheromone communication.

Substrate prediction of Ixodes ricinus salivary lipocalins differentially expressed during Borrelia afzelii infection

Evolution has provided ticks with an arsenal of bioactive saliva molecules that counteract host defense mechanisms. This salivary pharmacopoeia enables blood-feeding while enabling pathogen transmission. High-throughput sequencing of tick salivary glands has thus become a major focus, revealing large expansion within protein encoding gene families. Among these are lipocalins, ubiquitous barrel-shaped proteins that sequester small, typically hydrophobic molecules. This study was initiated by mining the Ixodes ricinus salivary gland transcriptome for specific, uncharacterized lipocalins: three were identified. Differential expression of these I. ricinus lipocalins during feeding at distinct developmental stages and in response to Borrelia afzelii infection suggests a role in transmission of this Lyme disease spirochete. A phylogenetic analysis using 803 sequences places the three I. ricinus lipocalins with tick lipocalins that sequester monoamines, leukotrienes and fatty acids. Both structural analysis and biophysical simulations generated robust predictions showing these I. ricinus lipocalins have the potential to bind monoamines similar to other tick species previously reported. The multidisciplinary approach employed in this study characterized unique lipocalins that play a role in tick blood-feeding and transmission of the most important tick-borne pathogen in North America and Eurasia.

The utility of protein structure as a predictor of site-wise dN/dS varies widely among HIV-1 proteins

Protein structure acts as a general constraint on the evolution of viral proteins. One widely recognized structural constraint explaining evolutionary variation among sites is the relative solvent accessibility (RSA) of residues in the folded protein. In influenza virus, the distance from functional sites has been found to explain an additional portion of the evolutionary variation in the external antigenic proteins. However, to what extent RSA and distance from a reference site in the protein can be used more generally to explain protein adaptation in other viruses and in the different proteins of any given virus remains an open question. To address this question, we have carried out an analysis of the distribution and structural predictors of site-wise dN/dS in HIV-1. Our results indicate that the distribution of dN/dS in HIV follows a smooth gamma distribution, with no special enrichment or depletion of sites with dN/dS at or above one. The variation in dN/dS can be partially explained by RSA and distance from a reference site in the protein, but these structural constraints do not act uniformly among the different HIV-1 proteins. Structural constraints are highly predictive in just one of the three enzymes and one of three structural proteins in HIV-1. For these two proteins, the protease enzyme and the gp120 structural protein, structure explains between 30 and 40% of the variation in dN/dS. Finally, for the gp120 protein of the receptor-binding complex, we also find that glycosylation sites explain just 2% of the variation in dN/dS and do not explain gp120 evolution independently of either RSA or distance from the apical surface.

Phytophthora infestans Argonaute 1 binds microRNA and small RNAs from effector genes and transposable elements

Phytophthora spp. encode large sets of effector proteins and distinct populations of small RNAs (sRNAs). Recent evidence has suggested that pathogen-derived sRNAs can modulate the expression of plant defense genes. Here, we studied the sRNA classes and functions associated with Phytophthora infestans Argonaute (Ago) proteins. sRNAs were co-immunoprecipitated with three PiAgo proteins and deep sequenced. Twenty- to twenty-two-nucleotide (nt) sRNAs were identified as the main interaction partners of PiAgo1 and high enrichment of 24-26-nt sRNAs was seen in the PiAgo4-bound sample. The frequencies and sizes of transposable element (TE)-derived sRNAs in the different PiAgo libraries suggested diversified roles of the PiAgo proteins in the control of different TE classes. We further provide evidence for the involvement of PiAgo1 in the P. infestans microRNA (miRNA) pathway. Protein-coding genes are probably regulated by the shared action of PiAgo1 and PiAgo5, as demonstrated by analysis of differential expression. An abundance of sRNAs from genes encoding host cell death-inducing Crinkler (CRN) effectors was bound to PiAgo1, implicating this protein in the regulation of the expanded CRN gene family. The data suggest that PiAgo1 plays an essential role in gene regulation and that at least two RNA silencing pathways regulate TEs in the plant-pathogenic oomycete P. infestans.

Phylogenetically Structured Differences in rRNA Gene Sequence Variation among Species of Arbuscular Mycorrhizal Fungi and Their Implications for Sequence Clustering

Arbuscular mycorrhizal (AM) fungi form mutualisms with plant roots that increase plant growth and shape plant communities. Each AM fungal cell contains a large amount of genetic diversity, but it is unclear if this diversity varies across evolutionary lineages. We found that sequence variation in the nuclear large-subunit (LSU) rRNA gene from 29 isolates representing 21 AM fungal species generally assorted into genus- and species-level clades, with the exception of species of the genera Claroideoglomus and Entrophospora However, there were significant differences in the levels of sequence variation across the phylogeny and between genera, indicating that it is an evolutionarily constrained trait in AM fungi. These consistent patterns of sequence variation across both phylogenetic and taxonomic groups pose challenges to interpreting operational taxonomic units (OTUs) as approximations of species-level groups of AM fungi. We demonstrate that the OTUs produced by five sequence clustering methods using 97% or equivalent sequence similarity thresholds failed to match the expected species of AM fungi, although OTUs from AbundantOTU, CD-HIT-OTU, and CROP corresponded better to species than did OTUs from mothur or UPARSE. This lack of OTU-to-species correspondence resulted both from sequences of one species being split into multiple OTUs and from sequences of multiple species being lumped into the same OTU. The OTU richness therefore will not reliably correspond to the AM fungal species richness in environmental samples. Conservatively, this error can overestimate species richness by 4-fold or underestimate richness by one-half, and the direction of this error will depend on the genera represented in the sample.

IMPORTANCE

Arbuscular mycorrhizal (AM) fungi form important mutualisms with the roots of most plant species. Individual AM fungi are genetically diverse, but it is unclear whether the level of this diversity differs among evolutionary lineages. We found that the amount of sequence variation in an rRNA gene that is commonly used to identify AM fungal species varied significantly between evolutionary groups that correspond to different genera, with the exception of two genera that are genetically indistinguishable from each other. When we clustered groups of similar sequences into operational taxonomic units (OTUs) using five different clustering methods, these patterns of sequence variation caused the number of OTUs to either over- or underestimate the actual number of AM fungal species, depending on the genus. Our results indicate that OTU-based inferences about AM fungal species composition from environmental sequences can be improved if they take these taxonomically structured patterns of sequence variation into account.

Application of the MAFFT sequence alignment program to large data-reexamination of the usefulness of chained guide trees

Motivation: Large multiple sequence alignments (MSAs), consisting of thousands of sequences, are becoming more and more common, due to advances in sequencing technologies. The MAFFT MSA program has several options for building large MSAs, but their performances have not been sufficiently assessed yet, because realistic benchmarking of large MSAs has been difficult. Recently, such assessments have been made possible through the HomFam and ContTest benchmark protein datasets. Along with the development of these datasets, an interesting theory was proposed: chained guide trees increase the accuracy of MSAs of structurally conserved regions. This theory challenges the basis of progressive alignment methods and needs to be examined by being compared with other known methods including computationally intensive ones.

Results: We used HomFam, ContTest and OXFam (an extended version of OXBench) to evaluate several methods enabled in MAFFT: (1) a progressive method with approximate guide trees, (2) a progressive method with chained guide trees, (3) a combination of an iterative refinement method and a progressive method and (4) a less approximate progressive method that uses a rigorous guide tree and consistency score. Other programs, Clustal Omega and UPP, available for large MSAs, were also included into the comparison. The effect of method 2 (chained guide trees) was positive in ContTest but negative in HomFam and OXFam. Methods 3 and 4 increased the benchmark scores more consistently than method 2 for the three datasets, suggesting that they are safer to use.

Availability and Implementation:http://mafft.cbrc.jp/alignment/software/

Contact:katoh@ifrec.osaka-u.ac.jp

Supplementary information:Supplementary data are available at Bioinformatics online.

Metagenomics Reveals Pervasive Bacterial Populations and Reduced Community Diversity across the Alaska Tundra Ecosystem

How soil microbial communities contrast with respect to taxonomic and functional composition within and between ecosystems remains an unresolved question that is central to predicting how global anthropogenic change will affect soil functioning and services. In particular, it remains unclear how small-scale observations of soil communities based on the typical volume sampled (1-2 g) are generalizable to ecosystem-scale responses and processes. This is especially relevant for remote, northern latitude soils, which are challenging to sample and are also thought to be more vulnerable to climate change compared to temperate soils. Here, we employed well-replicated shotgun metagenome and 16S rRNA gene amplicon sequencing to characterize community composition and metabolic potential in Alaskan tundra soils, combining our own datasets with those publically available from distant tundra and temperate grassland and agriculture habitats. We found that the abundance of many taxa and metabolic functions differed substantially between tundra soil metagenomes relative to those from temperate soils, and that a high degree of OTU-sharing exists between tundra locations. Tundra soils were an order of magnitude less complex than their temperate counterparts, allowing for near-complete coverage of microbial community richness (~92% breadth) by sequencing, and the recovery of 27 high-quality, almost complete (>80% completeness) population bins. These population bins, collectively, made up to ~10% of the metagenomic datasets, and represented diverse taxonomic groups and metabolic lifestyles tuned toward sulfur cycling, hydrogen metabolism, methanotrophy, and organic matter oxidation. Several population bins, including members of Acidobacteria, Actinobacteria, and Proteobacteria, were also present in geographically distant (~100-530 km apart) tundra habitats (full genome representation and up to 99.6% genome-derived average nucleotide identity). Collectively, our results revealed that Alaska tundra microbial communities are less diverse and more homogenous across spatial scales than previously anticipated, and provided DNA sequences of abundant populations and genes that would be relevant for future studies of the effects of environmental change on tundra ecosystems.

Xerotolerant Cladosporium sphaerospermum Are Predominant on Indoor Surfaces Compared to Other Cladosporium Species

Indoor fungi are a major cause of cosmetic and structural damage of buildings worldwide and prolonged exposure of these fungi poses a health risk. Aspergillus, Penicillium and Cladosporium species are the most predominant fungi in indoor environments. Cladosporium species predominate under ambient conditions. A total of 123 Cladosporium isolates originating from indoor air and indoor surfaces of archives, industrial factories, laboratories, and other buildings from four continents were identified by sequencing the internal transcribed spacer (ITS), and a part of the translation elongation factor 1α gene (TEF) and actin gene (ACT). Species from the Cladosporium sphaerospermum species complex were most predominant representing 44.7% of all isolates, while the Cladosporium cladosporioides and Cladosporium herbarum species complexes represented 33.3% and 22.0%, respectively. The contribution of the C. sphaerospermum species complex was 23.1% and 58.2% in the indoor air and isolates from indoor surfaces, respectively. Isolates from this species complex showed growth at lower water activity (≥ 0.82) when compared to species from the C. cladosporioides and C. herbarum species complexes (≥ 0.85). Together, these data indicate that xerotolerance provide the C. sphaerospermum species complex advantage in colonizing indoor surfaces. As a consequence, C. sphaerospermum are proposed to be the most predominant fungus at these locations under ambient conditions. Findings are discussed in relation to the specificity of allergy test, as the current species of Cladosporium used to develop these tests are not the predominant indoor species.

Halophilic and thermotolerant Gymnoascus species from several special environments, China

This study introduces three new Gymnoascus species (Gymnoascaceae, Onygenales), G. halophilus, G. stercorarius and G. thermotolerans, isolated from sediments in Chaka Salt Lake, compost and cornfield soil, respectively, in China, based on a polyphasic characterization including morphology, physiology and molecular phylogeny. Phylogenetic relationships were assessed based on the nuclear internal transcribed spacer (ITS = ITS1 + 5.8S + ITS2) region and a combined multilocus alignment of the ITS, 18S subunit rRNA gene and 28S subunit rRNA genes. Our study identified phylogenetic and phenotypic characters that differentiated the three new species from known species in the genus. Salinity and temperature tolerance tests revealed that G. halophilus was an obligate halophile while G. stercorarius and G. thermotolerans were halotolerant and thermotolerant. A key to accepted species of Gymnoascus is provided.

Identification and characterization of a novel 34 kDa merozoite protein in Babesia orientalis

A novel Babesia orientalis 34 kDa protein (designated BoP34) was obtained by immunoscreening of a cDNA expression library using B. orientalis infected water buffalo serum. The complete nucleotide sequence of the BoP34 was 1088 bp, which contained one open reading frame (ORF), two untranslated regions (UTRs) and a poly (A) tail. The length of ORF was 933 bp, encoding a polypeptide of 310 aa with a predicted size of 34 kDa. BLAST analysis showed that the nucleotide sequence of BoP34 had 71% similarity with that of the Babesia bovis gene XM_001611335, which encodes a nuclear movement family protein. This suggested that BoP34 is a homologous of the movement family protein. Structural analysis of the BoP34 protein indicated a CS domain which may interact with the ATPase domain of the heat shock protein 90. A truncated version of BoP34 was cloned into the expression vector pET-32a and subsequently expressed and purified from the Escherichia coli Rosetta™ (DE3) pLysS stain as a Trx-fusion (designated rBoP34-T). Antibodies in the serum of a B. orientalis-infected water buffalo were able to recognize this protein in immune-bloting analysis. Rabbit antibodies raised against rBoP34-T could detecte native BoP34 (34 kDa) in B. orientalis-infected water buffalo erythrocytes. These results suggested that BoP34 might be a good diagnostic antigen for the specific detection of anti-B. orientalis antibody in water buffalo. Further research is required to explore the biological function and diagnostic potential of this molecule.

PlantOrDB: a genome-wide ortholog database for land plants and green algae

BACKGROUND

Genes with different functions are originally generated from some ancestral genes by gene duplication, mutation and functional recombination. It is widely accepted that orthologs are homologous genes evolved from speciation events while paralogs are homologous genes resulted from gene duplication events.With the rapid increase of genomic data, identifying and distinguishing these genes among different species is becoming an important part of functional genomics research.

DESCRIPTION

Using 35 plant and 6 green algal genomes from Phytozome v9, we clustered 1,291,670 peptide sequences into 49,355 homologous gene families in terms of sequence similarity. For each gene family, we have generated a peptide sequence alignment and phylogenetic tree, and identified the speciation/duplication events for every node within the tree. For each node, we also identified and highlighted diagnostic characters that facilitate appropriate addition of a new query sequence into the existing phylogenetic tree and sequence alignment of its best matched gene family. Based on a desired species or subgroup of all species, users can view the phylogenetic tree, sequence alignment and diagnostic characters for a given gene family selectively. PlantOrDB not only allows users to identify orthologs or paralogs from phylogenetic trees, but also provides all orthologs that are built using Reciprocal Best Hit (RBH) pairwise alignment method. Users can upload their own sequences to find the best matched gene families, and visualize their query sequences within the relevant phylogenetic trees and sequence alignments.

CONCLUSION

PlantOrDB ( http://bioinfolab.miamioh.edu/plantordb ) is a genome-wide ortholog database for land plants and green algae. PlantOrDB offers highly interactive visualization, accurate query classification and powerful search functions useful for functional genomic research.

Population structure of mitochondrial genomes in Saccharomyces cerevisiae

Background

Rigorous study of mitochondrial functions and cell biology in the budding yeast, Saccharomyces cerevisiae has advanced our understanding of mitochondrial genetics. This yeast is now a powerful model for population genetics, owing to large genetic diversity and highly structured populations among wild isolates. Comparative mitochondrial genomic analyses between yeast species have revealed broad evolutionary changes in genome organization and architecture. A fine-scale view of recent evolutionary changes within S. cerevisiae has not been possible due to low numbers of complete mitochondrial sequences.

Results

To address challenges of sequencing AT-rich and repetitive mitochondrial DNAs (mtDNAs), we sequenced two divergent S. cerevisiae mtDNAs using a single-molecule sequencing platform (PacBio RS). Using de novo assemblies, we generated highly accurate complete mtDNA sequences. These mtDNA sequences were compared with 98 additional mtDNA sequences gathered from various published collections. Phylogenies based on mitochondrial coding sequences and intron profiles revealed that intraspecific diversity in mitochondrial genomes generally recapitulated the population structure of nuclear genomes. Analysis of intergenic sequence indicated a recent expansion of mobile elements in certain populations. Additionally, our analyses revealed that certain populations lacked introns previously believed conserved throughout the species, as well as the presence of introns never before reported in S. cerevisiae.

Conclusions

Our results revealed that the extensive variation in S. cerevisiae mtDNAs is often population specific, thus offering a window into the recent evolutionary processes shaping these genomes. In addition, we offer an effective strategy for sequencing these challenging AT-rich mitochondrial genomes for small scale projects.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1664-4) contains supplementary material, which is available to authorized users.

Microbiology, genomics, and clinical significance of the Pseudomonas fluorescens species complex, an unappreciated colonizer of humans

Pseudomonas fluorescens is not generally considered a bacterial pathogen in humans; however, multiple culture-based and culture-independent studies have identified it at low levels in the indigenous microbiota of various body sites. With recent advances in comparative genomics, many isolates originally identified as the "species" P. fluorescens are now being reclassified as novel Pseudomonas species within the P. fluorescens "species complex." Although most widely studied for its role in the soil and the rhizosphere, P. fluorescens possesses a number of functional traits that provide it with the capability to grow and thrive in mammalian hosts. While significantly less virulent than P. aeruginosa, P. fluorescens can cause bacteremia in humans, with most reported cases being attributable either to transfusion of contaminated blood products or to use of contaminated equipment associated with intravenous infusions. Although not suspected of being an etiologic agent of pulmonary disease, there are a number of reports identifying it in respiratory samples. There is also an intriguing association between P. fluorescens and human disease, in that approximately 50% of Crohn's disease patients develop serum antibodies to P. fluorescens. Altogether, these reports are beginning to highlight a far more common, intriguing, and potentially complex association between humans and P. fluorescens during health and disease.

Improved genome inference in the MHC using a population reference graph

Although much is known about human genetic variation, such information is typically ignored in assembling new genomes. Instead, reads are mapped to a single reference, which can lead to poor characterization of regions of high sequence or structural diversity. We introduce a population reference graph, which combines multiple reference sequences and catalogs of variation. The genomes of new samples are reconstructed as paths through the graph using an efficient hidden Markov model, allowing for recombination between different haplotypes and additional variants. By applying the method to the 4.5-Mb extended MHC region on human chromosome 6, combining 8 assembled haplotypes, the sequences of known classical HLA alleles and 87,640 SNP variants from the 1000 Genomes Project, we demonstrate using simulations, SNP genotyping, and short-read and long-read data how the method improves the accuracy of genome inference and identified regions where the current set of reference sequences is substantially incomplete.

Seasonal diversity and dynamics of haptophytes in the Skagerrak, Norway, explored by high-throughput sequencing

Microalgae in the division Haptophyta play key roles in the marine ecosystem and in global biogeochemical processes. Despite their ecological importance, knowledge on seasonal dynamics, community composition and abundance at the species level is limited due to their small cell size and few morphological features visible under the light microscope. Here, we present unique data on haptophyte seasonal diversity and dynamics from two annual cycles, with the taxonomic resolution and sampling depth obtained with high-throughput sequencing. From outer Oslofjorden, S Norway, nano- and picoplanktonic samples were collected monthly for 2 years, and the haptophytes targeted by amplification of RNA/cDNA with Haptophyta-specific 18S rDNA V4 primers. We obtained 156 operational taxonomic units (OTUs), from c. 400.000 454 pyrosequencing reads, after rigorous bioinformatic filtering and clustering at 99.5%. Most OTUs represented uncultured and/or not yet 18S rDNA-sequenced species. Haptophyte OTU richness and community composition exhibited high temporal variation and significant yearly periodicity. Richness was highest in September–October (autumn) and lowest in April–May (spring). Some taxa were detected all year, such as Chrysochromulina simplex, Emiliania huxleyi and Phaeocystis cordata, whereas most calcifying coccolithophores only appeared from summer to early winter. We also revealed the seasonal dynamics of OTUs representing putative novel classes (clades HAP-3–5) or orders (clades D, E, F). Season, light and temperature accounted for 29% of the variation in OTU composition. Residual variation may be related to biotic factors, such as competition and viral infection. This study provides new, in-depth knowledge on seasonal diversity and dynamics of haptophytes in North Atlantic coastal waters.

Sequence analysis reveals a conserved extension in the capping enzyme of the alphavirus supergroup, and a homologous domain in nodaviruses

Background

Members of the alphavirus supergroup include human pathogens such as chikungunya virus, hepatitis E virus and rubella virus. They encode a capping enzyme with methyltransferase-guanylyltransferase (MTase-GTase) activity, which is an attractive drug target owing to its unique mechanism. However, its experimental study has proven very difficult.

Results

We examined over 50 genera of viruses by sequence analyses. Earlier studies showed that the MTase-GTase contains a “Core” region conserved in sequence. We show that it is followed by a long extension, which we termed “Iceberg” region, whose secondary structure, but not sequence, is strikingly conserved throughout the alphavirus supergroup. Sequence analyses strongly suggest that the minimal capping domain corresponds to the Core and Iceberg regions combined, which is supported by earlier experimental data. The Iceberg region contains all known membrane association sites that contribute to the assembly of viral replication factories. We predict that it may also contain an overlooked, widely conserved membrane-binding amphipathic helix. Unexpectedly, we detected a sequence homolog of the alphavirus MTase-GTase in taxa related to nodaviruses and to chronic bee paralysis virus. The presence of a capping enzyme in nodaviruses is biologically consistent, since they have capped genomes but replicate in the cytoplasm, where no cellular capping enzyme is present. The putative MTase-GTase domain of nodaviruses also contains membrane-binding sites that may drive the assembly of viral replication factories, revealing an unsuspected parallel with the alphavirus supergroup.

Conclusions

Our work will guide the functional analysis of the alphaviral MTase-GTase and the production of domains for structure determination. The identification of a homologous domain in a simple model system, nodaviruses, which replicate in numerous eukaryotic cell systems (yeast, flies, worms, mammals, and plants), can further help crack the function and structure of the enzyme.

Reviewers

This article was reviewed by Valerian Dolja, Eugene Koonin and Sebastian Maurer-Stroh.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0050-0) contains supplementary material, which is available to authorized users.

Analysis of new isolates reveals new genome organization and a hypervariable region in infectious myonecrosis virus (IMNV)

Infectious myonecrosis virus (IMNV) has been the cause of many losses in shrimp farming since 2002, when the first myonecrosis outbreak was reported at Brazilian's northeast coast. Two additional genomes of Brazilian IMNV isolates collected in 2009 and 2013 were sequenced and analyzed in the present study. The sequencing revealed extra 643 bp and 22 bp, at 5' and 3' ends of IMNV genome respectively, confirming that its actual size is at least 8226 bp long. Considering these additional sequences in genome extremities, ORF1 can starts at nt 470, encoding a 1708 aa polyprotein. Computational predictions reveal two stem loops and two pseudoknots in the 5' end and a putative stem loop and a slippery motif located at 3' end, indicating that these regions can be involved in the start and termination of translation. Through a careful phylogenetic analysis, a higher genetic variability among Brazilian isolates could be observed, comparing with Indonesian IMNV isolates. It was also observed that the most variable region of IMNV genome is located in the first half of ORF1, coinciding with a region which probably encodes the capsid protrusions. The results presented here are a starting point to elucidate the viral's translational regulation and the mechanisms involved in virulence.

MarkerMiner 1.0: A new application for phylogenetic marker development using angiosperm transcriptomes

PREMISE OF THE STUDY

Targeted sequencing using next-generation sequencing (NGS) platforms offers enormous potential for plant systematics by enabling economical acquisition of multilocus data sets that can resolve difficult phylogenetic problems. However, because discovery of single-copy nuclear (SCN) loci from NGS data requires both bioinformatics skills and access to high-performance computing resources, the application of NGS data has been limited.

METHODS AND RESULTS

We developed MarkerMiner 1.0, a fully automated, open-access bioinformatic workflow and application for discovery of SCN loci in angiosperms. Our new tool identified as many as 1993 SCN loci from transcriptomic data sampled as part of four independent test cases representing marker development projects at different phylogenetic scales.

CONCLUSIONS

MarkerMiner is an easy-to-use and effective tool for discovery of putative SCN loci. It can be run locally or via the Web, and its tabular and alignment outputs facilitate efficient downstream assessments of phylogenetic utility, locus selection, intron-exon boundary prediction, and primer or probe development.

Oceanographic structure drives the assembly processes of microbial eukaryotic communities

Arctic Ocean microbial eukaryote phytoplankton form subsurface chlorophyll maximum (SCM), where much of the annual summer production occurs. This SCM is particularly persistent in the Western Arctic Ocean, which is strongly salinity stratified. The recent loss of multiyear sea ice and increased particulate-rich river discharge in the Arctic Ocean results in a greater volume of fresher water that may displace nutrient-rich saltier waters to deeper depths and decrease light penetration in areas affected by river discharge. Here, we surveyed microbial eukaryotic assemblages in the surface waters, and within and below the SCM. In most samples, we detected the pronounced SCM that usually occurs at the interface of the upper mixed layer and Pacific Summer Water (PSW). Poorly developed SCM was seen under two conditions, one above PSW and associated with a downwelling eddy, and the second in a region influenced by the Mackenzie River plume. Four phylogenetically distinct communities were identified: surface, pronounced SCM, weak SCM and a deeper community just below the SCM. Distance-decay relationships and phylogenetic structure suggested distinct ecological processes operating within these communities. In the pronounced SCM, picophytoplanktons were prevalent and community assembly was attributed to water mass history. In contrast, environmental filtering impacted the composition of the weak SCM communities, where heterotrophic Picozoa were more numerous. These results imply that displacement of Pacific waters to greater depth and increased terrigenous input may act as a control on SCM development and result in lower net summer primary production with a more heterotroph dominated eukaryotic microbial community.

Cladograms with Path to Event (ClaPTE): a novel algorithm to detect associations between genotypes or phenotypes using phylogenies

BACKGROUND

Associations between genotype and phenotype provide insight into the evolution of pathogenesis, drug resistance, and the spread of pathogens between hosts. However, common ancestry can lead to apparent associations between biologically unrelated features. The novel method Cladograms with Path to Event (ClaPTE) detects associations between character-pairs (either a pair of mutations or a mutation paired with a phenotype) while adjusting for common ancestry, using phylogenetic trees.

METHODS

ClaPTE tests for character-pairs changing close together on the phylogenetic tree, consistent with an associated character-pair. ClaPTE is compared to three existing methods (independent contrasts, mixed model, and likelihood ratio) to detect character-pair associations adjusted for common ancestry. Comparisons utilize simulations on gene trees for: HIV Env, HIV promoter, and bacterial DnaJ and GuaB; and case studies for Oseltamavir resistance in Influenza, and for DnaJ and GuaB. Simulated data include both true-positive/associated character-pairs, and true-negative/not-associated character-pairs, used to assess type I (frequency of p-values in true-negatives) and type II (sensitivity to true-positives) error control.

RESULTS AND CONCLUSIONS

ClaPTE has competitive sensitivity and better type I error control than existing methods. In the Influenza/Oseltamavir case study, ClaPTE reports no new permissive mutations but detects associations between adjacent (in primary sequence) amino acid positions which other methods miss. In the DnaJ and GuaB case study, ClaPTE reports more frequent associations between positions both from the same protein family than between positions from different families, in contrast to other methods. In both case studies, the results from ClaPTE are biologically plausible.

Characterisation of a Babesia orientalis apical membrane antigen, and comparison of its orthologues among selected apicomplexans

In the present study, we identified and characterised the complete coding sequence of Babesia orientalis apical membrane antigen 1 (designated Bo-ama1); it is 1803bp in length and encodes a polypeptide of 601 amino acids (aa). The Bo-ama-1 gene product (Bo-AMA1) is predicted to be 67kDa in size and contains a signal peptide. Mature Bo-AMA1 is predicted to have one transmembrane region and a short cytoplasmic tail (C-terminal domain). The extracellular part of Bo-AMA1 has three functional domains (DI, DII and DIII) with 14 conserved cysteine residues. A Bo-AMA1 fragment containing all three of these domains (designated Bo-AMA1-DI/II/III) was cloned into the plasmid vector pET-28a and expressed as a recombinant (His-fusion) protein of 53kDa. Antibodies in the serum from a B. orientalis-infected water buffalo specifically recognised this protein in immunoblotting analysis. Rabbit antibodies raised against the recombinant protein were able to detect native Bo-AMA1 (67kDa) from erythrocytes of B. orientalis-infected water buffalo. Bo-AMA1 is a new member of the AMA1 family and might be a good antigen for the specific detection of antibodies produced in B. orientalis infected cattle. This protein is likely to play critical roles during host cell adherence and invasion by B. orientalis, as the AMA1s reported in other organisms such as Plasmodium falciparum and Toxoplasma gondii. Further research is required to explore the biological functions of this protein and to determine whether its immunisation can induce protective effects in water buffalo against B. orientalis infection.

Parallel Implementation of MAFFT on CUDA-Enabled Graphics Hardware

Multiple sequence alignment (MSA) constitutes an extremely powerful tool for many biological applications including phylogenetic tree estimation, secondary structure prediction, and critical residue identification. However, aligning large biological sequences with popular tools such as MAFFT requires long runtimes on sequential architectures. Due to the ever increasing sizes of sequence databases, there is increasing demand to accelerate this task. In this paper, we demonstrate how graphic processing units (GPUs), powered by the compute unified device architecture (CUDA), can be used as an efficient computational platform to accelerate the MAFFT algorithm. To fully exploit the GPU's capabilities for accelerating MAFFT, we have optimized the sequence data organization to eliminate the bandwidth bottleneck of memory access, designed a memory allocation and reuse strategy to make full use of limited memory of GPUs, proposed a new modified-run-length encoding (MRLE) scheme to reduce memory consumption, and used high-performance shared memory to speed up I/O operations. Our implementation tested in three NVIDIA GPUs achieves speedup up to 11.28 on a Tesla K20m GPU compared to the sequential MAFFT 7.015.

One origin for metallo-β-lactamase activity, or two? An investigation assessing a diverse set of reconstructed ancestral sequences based on a sample of phylogenetic trees

Bacteria use metallo-β-lactamase enzymes to hydrolyse lactam rings found in many antibiotics, rendering them ineffective. Metallo-β-lactamase activity is thought to be polyphyletic, having arisen on more than one occasion within a single functionally diverse homologous superfamily. Since discovery of multiple origins of enzymatic activity conferring antibiotic resistance has broad implications for the continued clinical use of antibiotics, we test the hypothesis of polyphyly further; if lactamase function has arisen twice independently, the most recent common ancestor (MRCA) is not expected to possess lactam-hydrolysing activity. Two major problems present themselves. Firstly, even with a perfectly known phylogeny, ancestral sequence reconstruction is error prone. Secondly, the phylogeny is not known, and in fact reconstructing a single, unambiguous phylogeny for the superfamily has proven impossible. To obtain a more statistical view of the strength of evidence for or against MRCA lactamase function, we reconstructed a sample of 98 MRCAs of the metallo-β-lactamases, each based on a different tree in a bootstrap sample of reconstructed phylogenies. InterPro sequence signatures and homology modelling were then used to assess our sample of MRCAs for lactamase functionality. Only 5 % of these models conform to our criteria for metallo-β-lactamase functionality, suggesting that the ancestor was unlikely to have been a metallo-β-lactamase. On the other hand, given that ancestral proteins may have had metallo-β-lactamase functionality with variation in sequence and structural properties compared with extant enzymes, our criteria are conservative, estimating a lower bound of evidence for metallo-β-lactamase functionality but not an upper bound.

Deep-branching novel lineages and high diversity of haptophytes in the Skagerrak (Norway) uncovered by 454 pyrosequencing

Microalgae in the division Haptophyta may be difficult to identify to species by microscopy because they are small and fragile. Here, we used high-throughput sequencing to explore the diversity of haptophytes in outer Oslofjorden, Skagerrak, and supplemented this with electron microscopy. Nano- and picoplanktonic subsurface samples were collected monthly for 2 yr, and the haptophytes were targeted by amplification of RNA/cDNA with Haptophyta-specific 18S ribosomal DNA V4 primers. Pyrosequencing revealed higher species richness of haptophytes than previously observed in the Skagerrak by microscopy. From ca. 400,000 reads we obtained 156 haptophyte operational taxonomic units (OTUs) after rigorous filtering and 99.5% clustering. The majority (84%) of the OTUs matched environmental sequences not linked to a morphological species, most of which were affiliated with the order Prymnesiales. Phylogenetic analyses including Oslofjorden OTUs and available cultured and environmental haptophyte sequences showed that several of the OTUs matched sequences forming deep-branching lineages, potentially representing novel haptophyte classes. Pyrosequencing also retrieved cultured species not previously reported by microscopy in the Skagerrak. Electron microscopy revealed species not yet genetically characterised and some potentially novel taxa. This study contributes to linking genotype to phenotype within this ubiquitous and ecologically important protist group, and reveals great, unknown diversity.

Polyphasic characterization of four new plant pathogenic Phyllosticta species from China, Japan, and the United States

The black rot disease of Vitis species and other host genera of Vitacease is caused by Phyllosticta ampelicida and allied taxa which is considered to be a species complex. In this paper, we introduce four new species of Phyllosticta, including two from the P. ampelicida complex, based on a polyphasic characterization including disease symptoms and host association, morphology, and molecular phylogeny. The phylogenetic analysis was conducted based on the ribosomal internal transcribed spacer (ITS) region and a combined multi-locus alignment of the ITS, actin (ACT), partial translation elongation factor 1-alpha (TEF-1), and glyceraldehydes 3-phosphate dehydrogenase (GPDH) gene regions. Our study confirms the phylogenetic distinctions of the four new species, as well as their phenotypic differences with known species in the genus.

Molecular cloning and characterization of Babesia orientalis rhoptry-associated protein 1

The rhoptry-associated protein 1 (RAP-1) gene of Babesia orientalis was obtained from a cDNA expression library by immunoscreening with B. orientalis-infected water buffalo sera. The nucleotide sequence of the cDNA was 1732 bp with an open reading frame (ORF) of 1434 bp, encoding a polypeptide of 478 amino acid residues with a predicted size of 52.5 kDa. The ORF was cloned into a pGEX-KG plasmid and subsequently expressed as a GST-fusion protein. The recombinant RAP-1 of B. orientalis (rBoRAP-1) was purified and evaluated as an antigen using Western blotting. The native BoRAP-1 was recognized by the antibodies raised in rabbits against rBoRAP-1. Strong immunofluorescence signals were observed in erythrocytes infected with B. orientalis. Phylogentic analysis revealed that B. orientalis fell into a Babesia clade and most closely related to Babesia bovis and Babesia ovis, which was similar to the previous reported trees based on 18S rRNA and HSP70 genes. The present study suggests that the BoRAP-1 might be a potential diagnostic antigen, and the RAP-1 genes can aid in the classification of Babesia and Theileria species.

In silico single strand melting curve: a new approach to identify nucleic acid polymorphisms in Totiviridae

Background

The PCR technique and its variations have been increasingly used in the clinical laboratory and recent advances in this field generated new higher resolution techniques based on nucleic acid denaturation dynamics. The principle of these new molecular tools is based on the comparison of melting profiles, after denaturation of a DNA double strand. Until now, the secondary structure of single-stranded nucleic acids has not been exploited to develop identification systems based on PCR. To test the potential of single-strand RNA denaturation as a new alternative to detect specific nucleic acid variations, sequences from viruses of the Totiviridae family were compared using a new in silico melting curve approach. This family comprises double-stranded RNA virus, with a genome constituted by two ORFs, ORF1 and ORF2, which encodes the capsid/RNA binding proteins and an RNA-dependent RNA polymerase (RdRp), respectively.

Results

A phylogenetic tree based on RdRp amino acid sequences was constructed, and eight monophyletic groups were defined. Alignments of RdRp RNA sequences from each group were screened to identify RNA regions with conserved secondary structure. One region in the second half of ORF2 was identified and individually modeled using the RNAfold tool. Afterwards, each DNA or RNA sequence was denatured in silico using the softwares MELTSIM and RNAheat that generate melting curves considering the denaturation of a double stranded DNA and single stranded RNA, respectively. The same groups identified in the RdRp phylogenetic tree were retrieved by a clustering analysis of the melting curves data obtained from RNAheat. Moreover, the same approach was used to successfully discriminate different variants of Trichomonas vaginalis virus, which was not possible by the visual comparison of the double stranded melting curves generated by MELTSIM.

Conclusion

In silico analysis indicate that ssRNA melting curves are more informative than dsDNA melting curves. Furthermore, conserved RNA structures may be determined from analysis of individuals that are phylogenetically related, and these regions may be used to support the reconstitution of their phylogenetic groups. These findings are a robust basis for the development of in vitro systems to ssRNA melting curves detection.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-243) contains supplementary material, which is available to authorized users.

Purification, crystal structure determination and functional characterization of type III antifreeze proteins from the European eelpout Zoarces viviparus

Antifreeze proteins (AFPs) are essential components of many organisms adaptation to cold temperatures. Fish type III AFPs are divided into two groups, SP isoforms being much less active than QAE1 isoforms. Two type III AFPs from Zoarces viviparus, a QAE1 (ZvAFP13) and an SP (ZvAFP6) isoform, are here characterized and their crystal structures determined. We conclude that the higher activity of the QAE1 isoforms cannot be attributed to single residues, but rather a combination of structural effects. Furthermore both ZvAFP6 and ZvAFP13 crystal structures have water molecules around T18 equivalent to the tetrahedral-like waters previously identified in a neutron crystal structure. Interestingly, ZvAFP6 forms dimers in the crystal, with a significant dimer interface. The presence of ZvAFP6 dimers was confirmed in solution by native electrophoresis and gel filtration. To our knowledge this is the first report of dimerization of AFP type III proteins.

Predicting the fungal CUG codon translation with Bagheera

Background

Many eukaryotes have been shown to use alternative schemes to the universal genetic code. While most Saccharomycetes, including Saccharomyces cerevisiae, use the standard genetic code translating the CUG codon as leucine, some yeasts, including many but not all of the “Candida”, translate the same codon as serine. It has been proposed that the change in codon identity was accomplished by an almost complete loss of the original CUG codons, making the CUG positions within the extant species highly discriminative for the one or other translation scheme.

Results

In order to improve the prediction of genes in yeast species by providing the correct CUG decoding scheme we implemented a web server, called Bagheera, that allows determining the most probable CUG codon translation for a given transcriptome or genome assembly based on extensive reference data. As reference data we use 2071 manually assembled and annotated sequences from 38 cytoskeletal and motor proteins belonging to 79 yeast species. The web service includes a pipeline, which starts with predicting and aligning homologous genes to the reference data. CUG codon positions within the predicted genes are analysed with respect to amino acid similarity and CUG codon conservation in related species. In addition, the tRNA_CAG gene is predicted in genomic data and compared to known leu-tRNA_CAG and ser-tRNA_CAG genes. Bagheera can also be used to evaluate any mRNA and protein sequence data with the codon usage of the respective species. The usage of the system has been demonstrated by analysing six genomes not included in the reference data.

Conclusions

Gene prediction and consecutive comparison with reference data from other Saccharomycetes are sufficient to predict the most probable decoding scheme for CUG codons. This approach has been implemented into Bagheera (http://www.motorprotein.de/bagheera).

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-411) contains supplementary material, which is available to authorized users.

Mitochondrial genome of Babesia orientalis, apicomplexan parasite of water buffalo (Bubalus babalis, Linnaeus, 1758) endemic in China

BACKGROUND

Apicomplexan parasites of the genus Babesia, Theileria and Plasmodium are very closely related organisms. Interestingly, their mitochondrial (mt) genomes are highly divergent. Among Babesia, Babesia orientalis is a new species recently identified and specifically epidemic to the southern part of China, causing severe disease to water buffalo. However, no information on the mt genome of B. orientalis was available.

METHODS

Four pairs of primers were designed based on the full genome sequence of B. orientalis (unpublished data) and by aligning reported mt genomes of B. bovis, B. bigemina, and T. parva. The entire mt genome was amplified by four sets of PCR. The obtained mt genome was annotated by aligning with published apicomplexan mt genomes and Artemis software v11. Phylogenetic analysis was performed by using cox1 and cob amino acid sequences.

RESULTS

The complete mt genome of B. orientalis (Wuhan strain) was sequenced and characterized. The entire mt genome is 5996 bp in length with a linear form, containing three protein-coding genes including cytochrome c oxidase I (cox1), cytochrome b (cob) and cytochrome c oxidase III (cox3) and six rRNA large subunit gene fragments. The gene arrangement in B. orientalis mt genome is similar to those of B. bovis, B. gibsoni and Theileria parva, but different from those of T. orientalis, T. equi and Plasmodium falciparum. Comparative analysis indicated that cox1 and cob genes were more conserved than cox3. Phylogenetic analysis based on amino acid sequences of cox1, cob and cox1 + cob, respectively, revealed that B. orientalis fell into Babesia clade with the closest relationship to B. bovis.

CONCLUSIONS

The availability of the entire mt genome sequences of B. orientalis provides valuable information for future phylogenetic, population genetics and molecular epidemiological studies of apicomplexan parasites.

ANIA: ANnotation and Integrated Analysis of the 14-3-3 interactome

The dimeric 14-3-3 proteins dock onto pairs of phosphorylated Ser and Thr residues on hundreds of proteins, and thereby regulate many events in mammalian cells. To facilitate global analyses of these interactions, we developed a web resource named ANIA: ANnotation and Integrated Analysis of the 14-3-3 interactome, which integrates multiple data sets on 14-3-3-binding phosphoproteins. ANIA also pinpoints candidate 14-3-3-binding phosphosites using predictor algorithms, assisted by our recent discovery that the human 14-3-3-interactome is highly enriched in 2R-ohnologues. 2R-ohnologues are proteins in families of two to four, generated by two rounds of whole genome duplication at the origin of the vertebrate animals. ANIA identifies candidate ‘lynchpins’, which are 14-3-3-binding phosphosites that are conserved across members of a given 2R-ohnologue protein family. Other features of ANIA include a link to the catalogue of somatic mutations in cancer database to find cancer polymorphisms that map to 14-3-3-binding phosphosites, which would be expected to interfere with 14-3-3 interactions. We used ANIA to map known and candidate 14-3-3-binding enzymes within the 2R-ohnologue complement of the human kinome. Our projections indicate that 14-3-3s dock onto many more human kinases than has been realized. Guided by ANIA, PAK4, 6 and 7 (p21-activated kinases 4, 6 and 7) were experimentally validated as a 2R-ohnologue family of 14-3-3-binding phosphoproteins. PAK4 binding to 14-3-3 is stimulated by phorbol ester, and involves the ‘lynchpin’ site phosphoSer99 and a major contribution from Ser181. In contrast, PAK6 and PAK7 display strong phorbol ester-independent binding to 14-3-3, with Ser113 critical for the interaction with PAK6. These data point to differential 14-3-3 regulation of PAKs in control of cell morphology.

Database URL: https://ania-1433.lifesci.dundee.ac.uk/prediction/webserver/index.py

MAFFT: iterative refinement and additional methods

This chapter outlines several methods implemented in the MAFFT package. MAFFT is a popular multiple sequence alignment (MSA) program with various options for the progressive method, the iterative refinement method and other methods. We first outline basic usage of MAFFT and then describe recent practical extensions, such as dot plot and adjustment of direction in DNA alignment. We also refer to MUSCLE, another high-performance MSA program.

Large-scale multiple sequence alignment and tree estimation using SATé

SATé is a method for estimating multiple sequence alignments and trees that has been shown to produce highly accurate results for datasets with large numbers of sequences. Running SATé using its default settings is very simple, but improved accuracy can be obtained by modifying its algorithmic parameters. We provide a detailed introduction to the algorithmic approach used by SATé, and instructions for running a SATé analysis using the GUI under default settings. We also provide a discussion of how to modify these settings to obtain improved results, and how to use SATé in a phylogenetic analysis pipeline.

Can RNA-Seq resolve the rapid radiation of advanced moths and butterflies (Hexapoda: Lepidoptera: Apoditrysia)? An exploratory study

Recent molecular phylogenetic studies of the insect order Lepidoptera have robustly resolved family-level divergences within most superfamilies, and most divergences among the relatively species-poor early-arising superfamilies. In sharp contrast, relationships among the superfamilies of more advanced moths and butterflies that comprise the mega-diverse clade Apoditrysia (ca. 145,000 spp.) remain mostly poorly supported. This uncertainty, in turn, limits our ability to discern the origins, ages and evolutionary consequences of traits hypothesized to promote the spectacular diversification of Apoditrysia. Low support along the apoditrysian "backbone" probably reflects rapid diversification. If so, it may be feasible to strengthen resolution by radically increasing the gene sample, but case studies have been few. We explored the potential of next-generation sequencing to conclusively resolve apoditrysian relationships. We used transcriptome RNA-Seq to generate 1579 putatively orthologous gene sequences across a broad sample of 40 apoditrysians plus four outgroups, to which we added two taxa from previously published data. Phylogenetic analysis of a 46-taxon, 741-gene matrix, resulting from a strict filter that eliminated ortholog groups containing any apparent paralogs, yielded dramatic overall increase in bootstrap support for deeper nodes within Apoditrysia as compared to results from previous and concurrent 19-gene analyses. High support was restricted mainly to the huge subclade Obtectomera broadly defined, in which 11 of 12 nodes subtending multiple superfamilies had bootstrap support of 100%. The strongly supported nodes showed little conflict with groupings from previous studies, and were little affected by changes in taxon sampling, suggesting that they reflect true signal rather than artifacts of massive gene sampling. In contrast, strong support was seen at only 2 of 11 deeper nodes among the "lower", non-obtectomeran apoditrysians. These represent a much harder phylogenetic problem, for which one path to resolution might include further increase in gene sampling, together with improved orthology assignments.

The hymenopteran tree of life: evidence from protein-coding genes and objectively aligned ribosomal data

Previous molecular analyses of higher hymenopteran relationships have largely been based on subjectively aligned ribosomal sequences (18S and 28S). Here, we reanalyze the 18S and 28S data (unaligned about 4.4 kb) using an objective and a semi-objective alignment approach, based on MAFFT and BAli-Phy, respectively. Furthermore, we present the first analyses of a substantial protein-coding data set (4.6 kb from one mitochondrial and four nuclear genes). Our results indicate that previous studies may have suffered from inflated support values due to subjective alignment of the ribosomal sequences, but apparently not from significant biases. The protein data provide independent confirmation of several earlier results, including the monophyly of non-xyelid hymenopterans, Pamphilioidea + Unicalcarida, Unicalcarida, Vespina, Apocrita, Proctotrupomorpha and core Proctotrupomorpha. The protein data confirm that Aculeata are nested within a paraphyletic Evaniomorpha, but cast doubt on the monophyly of Evanioidea. Combining the available morphological, ribosomal and protein-coding data, we examine the total-evidence signal as well as congruence and conflict among the three data sources. Despite an emerging consensus on many higher-level hymenopteran relationships, several problems remain unresolved or contentious, including rooting of the hymenopteran tree, relationships of the woodwasps, placement of Stephanoidea and Ceraphronoidea, and the sister group of Aculeata.

Fungal community analysis by high-throughput sequencing of amplified markers--a user's guide

Novel high-throughput sequencing methods outperform earlier approaches in terms of resolution and magnitude. They enable identification and relative quantification of community members and offer new insights into fungal community ecology. These methods are currently taking over as the primary tool to assess fungal communities of plant-associated endophytes, pathogens, and mycorrhizal symbionts, as well as free-living saprotrophs. Taking advantage of the collective experience of six research groups, we here review the different stages involved in fungal community analysis, from field sampling via laboratory procedures to bioinformatics and data interpretation. We discuss potential pitfalls, alternatives, and solutions. Highlighted topics are challenges involved in: obtaining representative DNA/RNA samples and replicates that encompass the targeted variation in community composition, selection of marker regions and primers, options for amplification and multiplexing, handling of sequencing errors, and taxonomic identification. Without awareness of methodological biases, limitations of markers, and bioinformatics challenges, large-scale sequencing projects risk yielding artificial results and misleading conclusions.