A simple method to control over-alignment in the MAFFT multiple sequence alignment program

QuickProbs 2: Towards rapid construction of high-quality alignments of large protein families

The ever-increasing size of sequence databases caused by the development of high throughput sequencing, poses to multiple alignment algorithms one of the greatest challenges yet. As we show, well-established techniques employed for increasing alignment quality, i.e., refinement and consistency, are ineffective when large protein families are investigated. We present QuickProbs 2, an algorithm for multiple sequence alignment. Based on probabilistic models, equipped with novel column-oriented refinement and selective consistency, it offers outstanding accuracy. When analysing hundreds of sequences, Quick-Probs 2 is noticeably better than ClustalΩ and MAFFT, the previous leaders for processing numerous protein families. In the case of smaller sets, for which consistency-based methods are the best performing, QuickProbs 2 is also superior to the competitors. Due to low computational requirements of selective consistency and utilization of massively parallel architectures, presented algorithm has similar execution times to ClustalΩ, and is orders of magnitude faster than full consistency approaches, like MSAProbs or PicXAA. All these make QuickProbs 2 an excellent tool for aligning families ranging from few, to hundreds of proteins.

A new subterranean species of Pseudocrangonyx from China with an identification key to all species of the genus (Crustacea, Amphipoda, Pseudocrangonyctidae)

A subterranean species of Pseudocrangonyxelegantulus Hou, sp. n. is described from caves of Wulongdong National Forest Park in Henan Province, China. Pseudocrangonyxelegantulus is characterized by both male and female with calceoli on antenna II; urosomite III dorsal margin without armature; uropod III with peduncle 0.30 times as long as outer ramus and terminal article of the outer ramus a little shorter than adjacent spines; telson cleft 0.27 of its length. Phylogenetic analysis based on 28S and COI sequences supported the species distinctness. A key to the genus Pseudocrangonyx with 22 species and a map of their distributions are provided.

Influenza A Viruses of Swine (IAV-S) in Vietnam from 2010 to 2015: Multiple Introductions of A(H1N1)pdm09 Viruses into the Pig Population and Diversifying Genetic Constellations of Enzootic IAV-S

Active surveillance of influenza A viruses of swine (IAV-S) involving 262 farms and 10 slaughterhouses in seven provinces in northern and southern Vietnam from 2010 to 2015 yielded 388 isolates from 32 farms; these viruses were classified into H1N1, H1N2, and H3N2 subtypes. Whole-genome sequencing followed by phylogenetic analysis revealed that the isolates represented 15 genotypes, according to the genetic constellation of the eight segments. All of the H1N1 viruses were entirely A(H1N1)pdm09 viruses, whereas all of the H1N2 and H3N2 viruses were reassortants among 5 distinct ancestral viruses: H1 and H3 triple-reassortant (TR) IAV-S that originated from North American pre-2009 human seasonal H1, human seasonal H3N2, and A(H1N1)pdm09 viruses. Notably, 93% of the reassortant IAV-S retained M genes that were derived from A(H1N1)pdm09, suggesting some advantage in terms of their host adaptation. Bayesian Markov chain Monte Carlo analysis revealed that multiple introductions of A(H1N1)pdm09 and TR IAV-S into the Vietnamese pig population have driven the genetic diversity of currently circulating Vietnamese IAV-S. In addition, our results indicate that a reassortant IAV-S with human-like H3 and N2 genes and an A(H1N1)pdm09 origin M gene likely caused a human case in Ho Chi Minh City in 2010. Our current findings indicate that human-to-pig transmission as well as cocirculation of different IAV-S have contributed to diversifying the gene constellations of IAV-S in Vietnam.

IMPORTANCE

This comprehensive genetic characterization of 388 influenza A viruses of swine (IAV-S) isolated through active surveillance of Vietnamese pig farms from 2010 through 2015 provides molecular epidemiological insight into the genetic diversification of IAV-S in Vietnam after the emergence of A(H1N1)pdm09 viruses. Multiple reassortments among A(H1N1)pdm09 viruses and enzootic IAV-S yielded 14 genotypes, 9 of which carried novel gene combinations. The reassortants that carried M genes derived from A(H1N1)pdm09 viruses became predominant, replacing those of the IAV-S that had been endemic in Vietnam since 2011. Notably, one of the novel reassortants likely caused a human case in Vietnam. Given that Vietnam is the second-largest pig-producing country in Asia, continued monitoring of IAV-S is highly important from the viewpoints of both the swine industry and human public health.

MSAProbs-MPI: parallel multiple sequence aligner for distributed-memory systems

MSAProbs is a state-of-the-art protein multiple sequence alignment tool based on hidden Markov models. It can achieve high alignment accuracy at the expense of relatively long runtimes for large-scale input datasets. In this work we present MSAProbs-MPI, a distributed-memory parallel version of the multithreaded MSAProbs tool that is able to reduce runtimes by exploiting the compute capabilities of common multicore CPU clusters. Our performance evaluation on a cluster with 32 nodes (each containing two Intel Haswell processors) shows reductions in execution time of over one order of magnitude for typical input datasets. Furthermore, MSAProbs-MPI using eight nodes is faster than the GPU-accelerated QuickProbs running on a Tesla K20. Another strong point is that MSAProbs-MPI can deal with large datasets for which MSAProbs and QuickProbs might fail due to time and memory constraints, respectively.

AVAILABILITY AND IMPLEMENTATION

Source code in C ++ and MPI running on Linux systems as well as a reference manual are available at http://msaprobs.sourceforge.net CONTACT: jgonzalezd@udc.esSupplementary information: Supplementary data are available at Bioinformatics online.

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.