MrBayes on a Graphics Processing Unit

Phylogenetic placement and comparative analysis of the mitochondrial genomes of Idiostoloidea (Hemiptera: Heteroptera)

The classification system and the higher level phylogenetic relationships of Pentatomomorpha, the second largest infraorder of Heteroptera (Insecta: Hemiptera), have been debated and remain controversial over decades. In particular, the placement and phylogenetic relationship of Idiostoloidea are not well resolved, which hampers a better understanding of the evolutionary history of Pentatomomorpha. In this study, for the first time, we reported the complete mitochondrial genome for two narrowly distributed families of Idiostoloidea (including Idiostolidae and Henicocoridae), respectively. The length of the mitochondrial genome of Monteithocoris hirsutus and Henicocoris sp. is 16,632 and 16,013 bp, respectively. The content of AT is ranging from 75.15% to 80.48%. The mitogenomic structure of Idiostoloidea is highly conservative and there are no gene arrangements. By using the Bayesian inference, maximum likelihood, and Bayesian site-heterogeneous mixture model, we inferred the phylogenetic relationships within Pentatomomorpha and estimated their divergence times based on concatenated mitogenomes and nuclear ribosomal genes. Our results support the classification system of six superfamilies within Pentatomomorpha and confirm the monophyletic groups of each superfamily, with the following phylogenetic relationships: (Aradoidea + (Pentatomoidea + (Idiostoloidea + (Coreoidea + (Pyrrhocoroidea + Lygaeoidea))))). Furthermore, estimated divergence times revealed that most pentatomomorphan superfamilies and families diverged during the Late Jurassic to Early Cretaceous, which coincides with the explosive radiation of angiosperms.

Two new species of Erythroneurini (Hemiptera, Cicadellidae, Typhlocybinae) from southern China based on morphology and complete mitogenomes

Erythroneurine leafhoppers (Hemiptera, Cicadellidae, Typhlocybinae, Erythroneurini) are utilized to resolve the relationship between the four erythroneurine leafhopper (Hemiptera, Cicadellidae, Typhlocybinae, Erythroneurini): Arboridia (Arboridia) rongchangensis sp. nov., Thaia (Thaia) jiulongensis sp. nov., Mitjaevia bifurcata Luo, Song & Song, 2021 and Mitjaevia diana Luo, Song & Song, 2021, the two new species are described and illustrated. The mitochondrial gene sequences of these four species were determined to update the mitochondrial genome database of Erythroneurini. The mitochondrial genomes of four species shared high parallelism in nucleotide composition, base composition and gene order, comprising 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs) and an AT control region, which was consistent with majority of species in Cicadellidae; all genes revealed common trait of a positive AT skew and negative GC skew. The mitogenomes of four species were ultra-conservative in structure, and which isanalogous to that of others in size and A + T content. Phylogenetic trees based on the mitogenome data of these species and another 24 species were built employing the maximum likelihood and Bayesian inference methods. The results indicated that the four species belong to the tribe Erythroneurini, M. diana is the sister-group relationship of M. protuberanta + M. bifurcata. The two species Arboridia (Arboridia) rongchangensis sp. nov. and Thaia (Thaia) jiulongensis sp. nov. also have a relatively close genetic relationship with the genus Mitjaevia.

An imported human case with the SARS-CoV-2 Omicron subvariant BA.2.75 in Yunnan Province, China

The Omicron variants spread rapidly worldwide after being initially detected in South Africa in November 2021. It showed increased transmissibility and immune evasion with far more amino acid mutations in the spike (S) protein than the previously circulating variants of concern (VOCs). Notably, on 15 July 2022, we monitored the first VOC / Omicron subvariant BA.2.75 in China from an imported case. Moreover, nowadays, this subvariant still is predominant in India. It has nine additional mutations in the S protein compared to BA.2, three of which (W152R, G446S, and R493Q reversion) might contribute to higher transmissibility and immune escape. This subvariant could cause wider spread and pose a threat to the global situation. Our timely reporting and continuous genomic analysis are essential to fully elucidate the characteristics of the subvariant BA.2.75 in the future.

The complete mitochondrial genomes of two erythroneurine leafhoppers (Hemiptera, Cicadellidae, Typhlocybinae, Erythroneurini) with assessment of the phylogenetic status and relationships of tribes of Typhlocybinae

The number and classification of tribes in the leafhopper subfamily Typhlocybinae are not yet fully clear, and molecular data has recently been used to help resolve the problem. In this study, the mitochondrial genomes of Mitjaeviashibingensis Chen, Song & Webb, 2020 and M.dworakowskae Chen, Song & Webb, 2020 of the tribe Erythroneurini (Cicadellidae, Typhlocybinae) were sequenced. Most protein-coding genes (PCGs) start with ATN and end with TAA or TAG, and the AT content of these three codons were found differ from previous results that show that the first codon has the highest incidence. Two rRNA genes are highly conserved, and the AT content in 16S is higher than that of 12S. The nucleotide diversity and genetic distance among 13 PCGs of the four tribes from Typhlocybinae show that Empoascini nucleotide diversity is significantly less than in the other three tribes, and have the largest distance from the others, while Typhlocybini and Zyginellini have the smallest distance, indicating that the relationship between the two is the closest. The nad2, nad4, nad4L, and nad5 genes have greater nucleotide diversity, showing potential for use as the main markers for species identification. The phylogenetic analysis yielded a well-supported topology with most branches receiving maximum support and a few branches pertaining to relationships within Zyginellini and Typhlocybini receiving lower support. The species of these two tribes are intertwined, and it was impossible to resolve them into separate branches. In addition, the tribes Empoascini and Erythroneurini were recovered as monophyletic, and Alebrini was placed at the base of the tree as the most primitive. These results are broadly in line with other molecular phylogenetical studies which differ from traditional morphological classification.

Structural features and phylogenetic implications of Cicadellidae subfamily and two new mitogenomes leafhoppers

Complete mitochondrial genome sequences facilitate species identification and analyses of phylogenetic relationships. However, the available data are limited to the diverse and widespread insect family Cicadellidae. This study analyzes and summarizes the complete mitochondrial genome structure characteristics of 11 leafhopper subfamilies and two newly sequenced Typhlocybinae species, Empoascanara wengangensis and E. gracilis. Moreover, using 13PCGs and rRNA data to analyze the nucleotide diversity, evolution rate, and the phylogenetic relationship between the subfamilies of 56 species, verifying the taxonomic status analysis of E. wengangensis and E. gracilis. The analysis results show that the genome structures of the subfamilies and the newly sequenced two species are very similar, and the size of the CR region is significantly related to the repeat unit. However, in the entire AT-skews and CG-skews, the AT-skews of other subfamilies are all positive, and CG-skews are negative, while Empoascini of Typhlocybinae and Ledrinae are the opposite. Furthermore, among 13PCGs, the AT-skews of 13 species are all negative while CG-skews are positive, which from Empoascini in Typhlocybinae, Idiocerinae, Cicadellinae, Ledrinae, and Evacanthinae. Phylogenetic analysis shows that ML and PB analysis produce almost consistent topologies between different data sets and models, and some relationships are highly supported and remain unchanged. Mileewinae is a monophyletic group and is a sister group with Typhlocybinae, and the sister group of Evacanthinae is Ledrinae + Cicadellinae. Phylogenetic analysis grouped the two newly sequenced species with other species of Typhlocybinae, which was separated from other subfamilies, and all Erythroneurini insects gathered together. However, E. gracilis grouped into a single group, not grouped with species of the same genus (Empoascanara). This result does not match the traditional classification, and other nuclear genes or transcriptome genes may be needed to verify the result. Nucleotide diversity analysis shows that nad4 and nad5 may be evaluated as potential DNA markers defining the Cicadellidae insect species.

Mitogenome and phylogenetic analysis of typhlocybine leafhoppers (Hemiptera: Cicadellidae)

Mitogenomes have been widely used to estimate phylogenetic relationships among insects and provide data useful for augmenting traditional morphological characters in delimiting species. Here, complete mitogenome sequences of two closely related typhlocybine leafhoppers, Cassianeuracassiae (Ahmed, 1970) and C. bimaculata Dworakowska, 1984, were obtained and found to be 15,423 bp and 14,597 bp in length, respectively. The gene order was found to be similar to other published leafhopper mitogenomes, but the control region of C.bimaculata is the shortest among known leafhoppers and lacks tandem repeats. Phylogenetic analysis of 13 protein-coding genes (PCGs), the first and second codons of 13 PCGs, 13 PCGs and two rRNAs formed three well-supported tree topologies. The topologies of phylogenetic trees inferred from three datasets were almost identical, which was consistent with previous molecular phylogenies of this group. Comparative morphological study of the ovipositors revealed several characters potentially useful for diagnosing genera and resolving their phylogenetic relationships. Phylogenetic analysis of these and other morphological characters yielded a tree that is mostly consistent with the tree obtained from analysis of mitogenome sequences. In both molecular and morphological phylogenenies, Typhlocybini and Zyginellini clustered into one clade, but neither was recovered as monophyletic.

Phylogenetic relationships of Grapsoidea and insights into the higher phylogeny of Brachyuran

Decapoda is one of the most diverse crustacean orders, and has become an important research subject. However, the phylogenetic relationships among the main lineages of Decapoda remain uncertain, especially in the order Brachyura. Herein, we sequenced the whole mitochondrial genome of V. litterata and constructed a phylogenetic tree to understand its phylogenetic relationships with other species. The results showed that the mitochondrial genome of V. litterata was generally similar to mitogenomes of Metazoa reported in the literature, with a size of 16,247 bp, 37 genes, and a control region. Both AT-skew and GC-skew were negative, indicating more abundant Cs and Ts than Gs and As. The gene arrangement of V. litterata is identical to those of Eriocheir hepuensis, Cyclograpsus granulosus, Hemigrapsus sanguineus, Helicana wuana, and Helice tientsinensis but differs from the pancrustacean ground pattern and typical arrangement of Brachyuran crabs. Phylogenetic reconstruction showed that V. litterata belongs to the Varunidae.

Whole-genome sequencing reveals origin and evolution of influenza A(H1N1)pdm09 viruses in Lincang, China, from 2014 to 2018

The continuous variation of the seasonal influenza viruses, particularly A(H1N1)pdm09, persistently threatens human life and health around the world. In local areas of southwest china, the large time-scale genomic research on A(H1N1)pdm09 is still insufficient. Here, we sequenced 45 whole-genome sequences of influenza A(H1N1)pdm09 viruses in Lincang, China, from 2014 to 2018, by next-generation sequencing technology to characterize molecular mechanisms of their origin and evolution. Our phylogenetic analyses suggest that the A(H1N1)pdm09 strains circulating in Lincang belong to clade 6B and the subclade 6B.1A predominates in 2018. Further, the strains in 2018 possess elevated evolutionary rate as compared to strains in other years. Several newly emerged mutations for HA (hemagglutinin) in 2018 are revealed (i.e., S183P and R221K). Intriguingly, the substitution R221K falls into the RBS (receptor binding site) of HA protein, which could affect antigenic properties of influenza A(H1N1)pdm09 viruses, and another substitution S183P near to RBS with a high covering frequency (11/14 strains) in 2018 is exactly located at the epitope B. Notably, the NA (neuraminidase) protein harbors a new mutation I23T, potentially involved in N-glycosylation. Based on the background with a higher evolutionary rate in 2018 strains, we deeply evaluate the potential vaccine efficacy against Lincang strains and discover a substantive decline of the vaccine efficacy in 2018. Our analyses reaffirm that the real-time molecular surveillance and timely updated vaccine strains for prevention and control of influenza A(H1N1)pdm09 are crucial in the future.

Phylogeny and diversification of the true water bugs (Insecta: Hemiptera: Heteroptera: Nepomorpha)

Climate fluctuations and tectonic reconfigurations associated with environmental changes play large roles in determining patterns of adaptation and diversification, but studies documenting how such drivers have shaped the evolutionary history and diversification dynamics of limnic organisms during the Mesozoic are scarce. Members of the heteropteran infraorder Nepomorpha, or aquatic bugs, are ideal for testing the effects of these determinants on their diversification pulses because most species are confined to aquatic environments during their entire life. The group has a relatively mature taxonomy and is well represented in the fossil record. We investigated the evolution of Nepomorpha based on phylogenetic analyses of morphological and molecular characters sampled from 115 taxa representing all 13 families and approximately 40% of recognized genera. Our results were largely congruent with the phylogenetic relationships inferred from morphology. A divergence dating analysis indicated that Nepomorpha began to diversify in the late Permian (approximately 263 Ma), and diversification analyses suggested that palaeoecological opportunities probably promoted lineage diversification in this group.

Compositional and mutational rate heterogeneity in mitochondrial genomes and its effect on the phylogenetic inferences of Cimicomorpha (Hemiptera: Heteroptera)

Background

Mitochondrial genome (mt-genome) data can potentially return artefactual relationships in the higher-level phylogenetic inference of insects due to the biases of accelerated substitution rates and compositional heterogeneity. Previous studies based on mt-genome data alone showed a paraphyly of Cimicomorpha (Insecta, Hemiptera) due to the positions of the families Tingidae and Reduviidae rather than the monophyly that was supported based on morphological characters, morphological and molecular combined data and large scale molecular datasets. Various strategies have been proposed to ameliorate the effects of potential mt-genome biases, including dense taxon sampling, removal of third codon positions or purine-pyrimidine coding and the use of site-heterogeneous models. In this study, we sequenced the mt-genomes of five additional Tingidae species and discussed the compositional and mutational rate heterogeneity in mt-genomes and its effect on the phylogenetic inferences of Cimicomorpha by implementing the bias-reduction strategies mentioned above.

Results

Heterogeneity in nucleotide composition and mutational biases were found in mt protein-coding genes, and the third codon exhibited high levels of saturation. Dense taxon sampling of Tingidae and Reduviidae and the other common strategies mentioned above were insufficient to recover the monophyly of the well-established group Cimicomorpha. When the sites with weak phylogenetic signals in the dataset were removed, the remaining dataset of mt-genomes can support the monophyly of Cimicomorpha; this support demonstrates that mt-genomes possess strong phylogenetic signals for the inference of higher-level phylogeny of this group. Comparison of the ratio of the removal of amino acids for each PCG showed that ATP8 has the highest ratio while CO1 has the lowest. This pattern is largely congruent with the evolutionary rate of 13 PCGs that ATP8 represents the highest evolutionary rate, whereas CO1 appears to be the lowest. Notably, the value of Ka/Ks ratios of all PCGs is less than 1, indicating that these genes are likely evolving under purifying selection.

Conclusions

Our results demonstrate that mt-genomes have sites with strong phylogenetic signals for the inference of higher-level phylogeny of Cimicomorpha. Consequently, bioinformatic approaches to removing sites with weak phylogenetic signals in mt-genome without relying on an a priori tree topology would greatly improve this field.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4650-9) contains supplementary material, which is available to authorized users.

The first complete mitochondrial genome from the subfamily Phyllocephalinae (Heteroptera: Pentatomidae) and its phylogenetic analysis

The true bug, Gonopsis affinis, is an important pest in China. Here, we determined the complete mitogenome of G. affinis, which is the first for the subfamily Phyllocephalinae. This 16,011-basepair (bp) mitogenome comprises of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes, and an A + T-rich region. The gene order and the orientation are similar to those of other sequenced Hemipteran species. All PCGs start with ATN codons, except COI, and end with TAA. The 22 tRNAs have a typical cloverleaf secondary structure except trnS-Ser (AGN), which lacks a dihydrouridine (DHU) arm. Phylogenetic analyses highly supported the monophyly of each family, and confirmed the reasonable placement of G. affinis.

The first mitochondrial genome for the subfamily Asopinae (Heteroptera: Pentatomidae) and its phylogenetic implications

Here, we determined the complete mitogenome of Picromerus griseus, the first for the subfamily Asopinae. This mitochondrial genome contains 16,338 bp, with an A + T content of 71.69%, and contains 37 typical mitochondrial genes (13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and a control region). The genome size, gene arrangement, A + T content, codon usage and secondary structures of 22 transfer RNA genes of the P. griseus mitogenome were similar to those of other sequenced pentatomoids. Bayesian analyses performed using the mitogenomes of P. griseus and its relatives, including 14 taxa, confirmed the reasonable placement of P. griseus.

A greedy alignment-free distance estimator for phylogenetic inference

BACKGROUND

Alignment-free sequence comparison approaches have been garnering increasing interest in various data- and compute-intensive applications such as phylogenetic inference for large-scale sequences. While k-mer based methods are predominantly used in real applications, the average common substring (ACS) approach is emerging as one of the prominent alignment-free approaches. This ACS approach has been further generalized by some recent work, either greedily or exactly, by allowing a bounded number of mismatches in the common substrings.

RESULTS

We present ALFRED-G, a greedy alignment-free distance estimator for phylogenetic tree reconstruction based on the concept of the generalized ACS approach. In this algorithm, we have investigated a new heuristic to efficiently compute the lengths of common strings with mismatches allowed, and have further applied this heuristic to phylogeny reconstruction. Performance evaluation using real sequence datasets shows that our heuristic is able to reconstruct comparable, or even more accurate, phylogenetic tree topologies than the kmacs heuristic algorithm at highly competitive speed.

CONCLUSIONS

ALFRED-G is an alignment-free heuristic for evolutionary distance estimation between two biological sequences. This algorithm is implemented in C++ and has been incorporated into our open-source ALFRED software package ( http://alurulab.cc.gatech.edu/phylo ).

A Mitochondrial Genome of Rhyparochromidae (Hemiptera: Heteroptera) and a Comparative Analysis of Related Mitochondrial Genomes

The Rhyparochromidae, the largest family of Lygaeoidea, encompasses more than 1,850 described species, but no mitochondrial genome has been sequenced to date. Here we describe the first mitochondrial genome for Rhyparochromidae: a complete mitochondrial genome of Panaorus albomaculatus (Scott, 1874). This mitochondrial genome is comprised of 16,345 bp, and contains the expected 37 genes and control region. The majority of the control region is made up of a large tandem-repeat region, which has a novel pattern not previously observed in other insects. The tandem-repeats region of P. albomaculatus consists of 53 tandem duplications (including one partial repeat), which is the largest number of tandem repeats among all the known insect mitochondrial genomes. Slipped-strand mispairing during replication is likely to have generated this novel pattern of tandem repeats. Comparative analysis of tRNA gene families in sequenced Pentatomomorpha and Lygaeoidea species shows that the pattern of nucleotide conservation is markedly higher on the J-strand. Phylogenetic reconstruction based on mitochondrial genomes suggests that Rhyparochromidae is not the sister group to all the remaining Lygaeoidea, and supports the monophyly of Lygaeoidea.

MrBayes tgMC3++: A High Performance and Resource-Efficient GPU-Oriented Phylogenetic Analysis Method

MrBayes is a widespread phylogenetic inference tool harnessing empirical evolutionary models and Bayesian statistics. However, the computational cost on the likelihood estimation is very expensive, resulting in undesirably long execution time. Although a number of multi-threaded optimizations have been proposed to speed up MrBayes, there are bottlenecks that severely limit the GPU thread-level parallelism of likelihood estimations. This study proposes a high performance and resource-efficient method for GPU-oriented parallelization of likelihood estimations. Instead of having to rely on empirical programming, the proposed novel decomposition storage model implements high performance data transfers implicitly. In terms of performance improvement, a speedup factor of up to 178 can be achieved on the analysis of simulated datasets by four Tesla K40 cards. In comparison to the other publicly available GPU-oriented MrBayes, the tgMC³++ method (proposed herein) outperforms the tgMC³ (v1.0), nMC³ (v2.1.1) and oMC³ (v1.00) methods by speedup factors of up to 1.6, 1.9 and 2.9, respectively. Moreover, tgMC³++ supports more evolutionary models and gamma categories, which previous GPU-oriented methods fail to take into analysis.

ALFRED: A Practical Method for Alignment-Free Distance Computation

Alignment-free approaches are gaining persistent interest in many sequence analysis applications such as phylogenetic inference and metagenomic classification/clustering, especially for large-scale sequence datasets. Besides the widely used k-mer methods, the average common substring (ACS) approach has emerged to be one of the well-known alignment-free approaches. Two recent works further generalize this ACS approach by allowing a bounded number k of mismatches in the common substrings, relying on approximation (linear time) and exact computation, respectively. Albeit having a good worst-case time complexity [Formula: see text], the exact approach is complex and unlikely to be efficient in practice. Herein, we present ALFRED, an alignment-free distance computation method, which solves the generalized common substring search problem via exact computation. Compared to the theoretical approach, our algorithm is easier to implement and more practical to use, while still providing highly competitive theoretical performances with an expected run-time of [Formula: see text]. By applying our program to phylogenetic inference as a case study, we find that our program facilitates to exactly reconstruct the topology of the reference phylogenetic tree for a set of 27 primate mitochondrial genomes, at reasonably acceptable speed. ALFRED is implemented in C++ programming language and the source code is freely available online.

The complete mitochondrial genome of Nepa hoffmanni (Hemiptera: Heteroptera: Nepidae)

The complete mitochondrial genome (mt-genome) of Nepa hoffmanni has been reported in this study. This mitochondrial genome is 15 774 bp long, with an A + T content of 72.04%, containing the typical 37 genes (13 protein-coding genes (PCGs), 22 transfer RNA genes, and two ribosomal RNA genes) and a control region. All genes are arranged in the same gene order as most other known heteropteran mt-genome. This is the second completely sequenced mt-genome from the family Nepidae of Nepomorpha. Bayesian analyses were performed using the mt-genome of Nepa hoffmanni and its relatives, including 17 taxa, showing a reasonable placement of Nepa hoffmanni.

GPU MrBayes V3.1: MrBayes on Graphics Processing Units for Protein Sequence Data

We present a modified GPU (graphics processing unit) version of MrBayes, called ta(MC)(3) (GPU MrBayes V3.1), for Bayesian phylogenetic inference on protein data sets. Our main contributions are 1) utilizing 64-bit variables, thereby enabling ta(MC)(3) to process larger data sets than MrBayes; and 2) to use Kahan summation to improve accuracy, convergence rates, and consequently runtime. Versus the current fastest software, we achieve a speedup of up to around 2.5 (and up to around 90 vs. serial MrBayes), and more on multi-GPU hardware. GPU MrBayes V3.1 is available from http://sourceforge.net/projects/mrbayes-gpu/.

The molecular symplesiomorphies shared by the stem groups of metazoan evolution: can sites as few as 1% have a significant impact on recognizing the phylogenetic position of myzostomida?

Although it is clear that taxon sampling, alignments, gene sampling, tree reconstruction methods and the total length of the sequences used are critical to the reconstruction of evolutionary history, weakly supported or misleading nodes exist in phylogenetic studies with no obvious flaw in those aspects. The phylogenetic studies focusing on the basal part of bilaterian evolution are such a case. During the past decade, Myzostomida has appeared in the basal part of Bilateria in several phylogenetic studies of Metazoa. However, most researchers have entertained only two competing hypotheses about the position of Myzostomida-an affinity with Annelida and an affinity with Platyhelminthes. In this study, dozens of symplesiomorphies were discovered by means of ancestral state reconstruction in the complete 18S and 28S rDNAs shared by the stem groups of Metazoa. By contrastive analysis on the datasets with or without such symplesiomorphic sites, we discovered that Myzostomida and other basal groups are basal lineages of Bilateria due to the corresponding symplesiomorphies shared with earlier lineages. As such, symplesiomorphies account for approximately 1-2% of the whole dataset have an essential impact on phylogenetic inference, and this study reminds molecular systematists of the importance of carrying out ancestral state reconstruction at each site in sequence-based phylogenetic studies. In addition, reasons should be explored for the low support of the hypothesis that Myzostomida belongs to Annelida in the results of phylogenomic studies. Future phylogenetic studies concerning Myzostomida should include all of the basal lineages of Bilateria to avoid directly neglecting the stand-alone basal position of Myzostomida as a potential hypothesis.

Long-branch attraction and the phylogeny of true water bugs (Hemiptera: Nepomorpha) as estimated from mitochondrial genomes

BACKGROUND

Most previous studies of morphological and molecular data have consistently supported the monophyly of the true water bugs (Hemiptera: Nepomorpha). An exception is a recent study by Hua et al. (BMC Evol Biol 9: 134, 2009) based on nine nepomorphan mitochondrial genomes. In the analysis of Hua et al. (BMC Evol Biol 9: 134, 2009), the water bugs in the group Pleoidea formed the sister group to a clade that consisted of Nepomorpha (the remaining true water bugs) + Leptopodomorpha (shore bugs) + Cimicomorpha (assassin bugs and relatives) + Pentatomomorpha (stink bugs and relatives), thereby suggesting that fully aquatic hemipterans evolved independently at least twice. Based on these results, Hua et al. (BMC Evol Biol 9: 134, 2009) elevated the Pleoidea to a new infraorder, the Plemorpha.

RESULTS

Our reanalysis suggests that the lack of support for the monophyly of the true water bugs (including Pleoidea) by Hua et al. (BMC Evol Biol 9: 134, 2009) likely resulted from inadequate taxon sampling. In particular, long-branch attraction (LBA) between the distant outgroup taxa and Pleoidea, as well as LBA among taxa in the ingroup, made Nepomorpha appear to be polyphyletic. We used three complementary strategies to test and alleviate the effects of LBA: (1) the removal of distant outgroups from the analysis; (2) the addition of closely related outgroups; and (3) the addition of a mitochondrial genome from a second family of Pleoidea. We also performed likelihood-ratio tests to examine the support for monophyly of Nepomorpha with different combinations of taxa included in the analysis. Furthermore, we found that specimens of Helotrephes sp. were misidentified as Paraplea frontalis (Fieber, 1844) by Hua et al. (BMC Evol Biol 9: 134, 2009).

CONCLUSIONS

All analyses that included the addition of more taxa significantly and consistently supported the placement of Pleoidea within the Nepomorpha (i.e., supported the monophyly of the traditional true water bugs). Our analyses further support a close relationship between Notonectoidea and Pleoidea within Nepomorpha, and the superfamilies Nepoidea, Ochteroidea, Naucoroidea, and Pleoidea are resolved as monophyletic in all trees with strong support. Our results also confirmed that monophyly of Nepomorpha clearly is not refuted by the mitochondrial genome data.

Discovering epistasis in large scale genetic association studies by exploiting graphics cards

Despite the enormous investments made in collecting DNA samples and generating germline variation data across thousands of individuals in modern genome-wide association studies (GWAS), progress has been frustratingly slow in explaining much of the heritability in common disease. Today's paradigm of testing independent hypotheses on each single nucleotide polymorphism (SNP) marker is unlikely to adequately reflect the complex biological processes in disease risk. Alternatively, modeling risk as an ensemble of SNPs that act in concert in a pathway, and/or interact non-additively on log risk for example, may be a more sensible way to approach gene mapping in modern studies. Implementing such analyzes genome-wide can quickly become intractable due to the fact that even modest size SNP panels on modern genotype arrays (500k markers) pose a combinatorial nightmare, require tens of billions of models to be tested for evidence of interaction. In this article, we provide an in-depth analysis of programs that have been developed to explicitly overcome these enormous computational barriers through the use of processors on graphics cards known as Graphics Processing Units (GPU). We include tutorials on GPU technology, which will convey why they are growing in appeal with today's numerical scientists. One obvious advantage is the impressive density of microprocessor cores that are available on only a single GPU. Whereas high end servers feature up to 24 Intel or AMD CPU cores, the latest GPU offerings from nVidia feature over 2600 cores. Each compute node may be outfitted with up to 4 GPU devices. Success on GPUs varies across problems. However, epistasis screens fare well due to the high degree of parallelism exposed in these problems. Papers that we review routinely report GPU speedups of over two orders of magnitude (>100x) over standard CPU implementations.

MrBayes tgMC³: a tight GPU implementation of MrBayes

MrBayes is model-based phylogenetic inference tool using Bayesian statistics. However, model-based assessment of phylogenetic trees adds to the computational burden of tree-searching, and so poses significant computational challenges. Graphics Processing Units (GPUs) have been proposed as high performance, low cost acceleration platforms and several parallelized versions of the Metropolis Coupled Markov Chain Mote Carlo (MC³) algorithm in MrBayes have been presented that can run on GPUs. However, some bottlenecks decrease the efficiency of these implementations. To address these bottlenecks, we propose a tight GPU MC³ (tgMC³) algorithm. tgMC³ implements a different architecture from the one-to-one acceleration architecture employed in previously proposed methods. It merges multiply discrete GPU kernels according to the data dependency and hence decreases the number of kernels launched and the complexity of data transfer. We implemented tgMC³ and made performance comparisons with an earlier proposed algorithm, nMC³, and also with MrBayes MC³ under serial and multiply concurrent CPU processes. All of the methods were benchmarked on the same computing node from DEGIMA. Experiments indicate that the tgMC³ method outstrips nMC³ (v1.0) with speedup factors from 2.1 to 2.7×. In addition, tgMC³ outperforms the serial MrBayes MC³ by a factor of 6 to 30× when using a single GTX480 card, whereas a speedup factor of around 51× can be achieved by using two GTX 480 cards on relatively long sequences. Moreover, tgMC³ was compared with MrBayes accelerated by BEAGLE, and achieved speedup factors from 3.7 to 5.7×. The reported performance improvement of tgMC³ is significant and appears to scale well with increasing dataset sizes. In addition, the strategy proposed in tgMC³ could benefit the acceleration of other Bayesian-based phylogenetic analysis methods using GPUs.

Efficient implementation of MrBayes on multi-GPU

MrBayes, using Metropolis-coupled Markov chain Monte Carlo (MCMCMC or (MC)³), is a popular program for Bayesian inference. As a leading method of using DNA data to infer phylogeny, the (MC)³ Bayesian algorithm and its improved and parallel versions are now not fast enough for biologists to analyze massive real-world DNA data. Recently, graphics processor unit (GPU) has shown its power as a coprocessor (or rather, an accelerator) in many fields. This article describes an efficient implementation a(MC)³ (aMCMCMC) for MrBayes (MC)³ on compute unified device architecture. By dynamically adjusting the task granularity to adapt to input data size and hardware configuration, it makes full use of GPU cores with different data sets. An adaptive method is also developed to split and combine DNA sequences to make full use of a large number of GPU cards. Furthermore, a new “node-by-node” task scheduling strategy is developed to improve concurrency, and several optimizing methods are used to reduce extra overhead. Experimental results show that a(MC)³ achieves up to 63× speedup over serial MrBayes on a single machine with one GPU card, and up to 170× speedup with four GPU cards, and up to 478× speedup with a 32-node GPU cluster. a(MC)³ is dramatically faster than all the previous (MC)³ algorithms and scales well to large GPU clusters.

The complete mitochondrial genome of the stalk-eyed bug Chauliops fallax Scott, and the monophyly of Malcidae (Hemiptera: Heteroptera)

Chauliops fallax Scott, 1874 (Hemiptera: Heteroptera: Malcidae: Chauliopinae) is one of the most destructive insect pests of soybean and rice fields in Asia. Here we sequenced the complete mitochondrial genome of this pest. This genome is 15,739 bp long, with an A+T content of 73.7%, containing 37 typical animal mitochondrial genes and a control region. All genes were arranged in the same order as most of other Heteroptera. A remarkable strand bias was found for all nine protein coding genes (PCGs) encoded by the majority strand were positive AT-skew and negative GC-skew, whereas the reverse were found in the remaining four PCGs encoded by the minority strand and two rRNA genes. The models of secondary structures for the two rRNA genes of sequenced true bugs and Lygaeoidea were predicted. 16S rRNA consisted of six domains (domain III is absent as in other known arthropod mitochondrial genomes) and 45 helices, while three domains and 27 helices for 12S rRNA. The control region consists of five subregions: a microsatellite-like region, a tandem repeats region and other three motifs. The unusual intergenic spacer between tRNA-H and ND4 only found in the species of Lygaeoidea, not in other heteropteran species, may be the synapomorphy of this superfamily. Phylogenetic analyses were carried out based on all the 13 PCGs showed that Chauliopinae was the sister group of Malcinae and the monophyly of Lygaeoidea.

Extending the BEAGLE library to a multi-FPGA platform

Background

Maximum Likelihood (ML)-based phylogenetic inference using Felsenstein’s pruning algorithm is a standard method for estimating the evolutionary relationships amongst a set of species based on DNA sequence data, and is used in popular applications such as RAxML, PHYLIP, GARLI, BEAST, and MrBayes. The Phylogenetic Likelihood Function (PLF) and its associated scaling and normalization steps comprise the computational kernel for these tools. These computations are data intensive but contain fine grain parallelism that can be exploited by coprocessor architectures such as FPGAs and GPUs. A general purpose API called BEAGLE has recently been developed that includes optimized implementations of Felsenstein’s pruning algorithm for various data parallel architectures. In this paper, we extend the BEAGLE API to a multiple Field Programmable Gate Array (FPGA)-based platform called the Convey HC-1.

Results

The core calculation of our implementation, which includes both the phylogenetic likelihood function (PLF) and the tree likelihood calculation, has an arithmetic intensity of 130 floating-point operations per 64 bytes of I/O, or 2.03 ops/byte. Its performance can thus be calculated as a function of the host platform’s peak memory bandwidth and the implementation’s memory efficiency, as 2.03 × peak bandwidth × memory efficiency. Our FPGA-based platform has a peak bandwidth of 76.8 GB/s and our implementation achieves a memory efficiency of approximately 50%, which gives an average throughput of 78 Gflops. This represents a ~40X speedup when compared with BEAGLE’s CPU implementation on a dual Xeon 5520 and 3X speedup versus BEAGLE’s GPU implementation on a Tesla T10 GPU for very large data sizes. The power consumption is 92 W, yielding a power efficiency of 1.7 Gflops per Watt.

Conclusions

The use of data parallel architectures to achieve high performance for likelihood-based phylogenetic inference requires high memory bandwidth and a design methodology that emphasizes high memory efficiency. To achieve this objective, we integrated 32 pipelined processing elements (PEs) across four FPGAs. For the design of each PE, we developed a specialized synthesis tool to generate a floating-point pipeline with resource and throughput constraints to match the target platform. We have found that using low-latency floating-point operators can significantly reduce FPGA area and still meet timing requirement on the target platform. We found that this design methodology can achieve performance that exceeds that of a GPU-based coprocessor.

A unique box in 28S rRNA is shared by the enigmatic insect order Zoraptera and Dictyoptera

The position of the Zoraptera remains one of the most challenging and uncertain concerns in ordinal-level phylogenies of the insects. Zoraptera have been viewed as having a close relationship with five different groups of Polyneoptera, or as being allied to the Paraneoptera or even Holometabola. Although rDNAs have been widely used in phylogenetic studies of insects, the application of the complete 28S rDNA are still scattered in only a few orders. In this study, a secondary structure model of the complete 28S rRNAs of insects was reconstructed based on all orders of Insecta. It was found that one length-variable region, D3-4, is particularly distinctive. The length and/or sequence of D3-4 is conservative within each order of Polyneoptera, but it can be divided into two types between the different orders of the supercohort, of which the enigmatic order Zoraptera and Dictyoptera share one type, while the remaining orders of Polyneoptera share the other. Additionally, independent evidence from phylogenetic results support the clade (Zoraptera+Dictyoptera) as well. Thus, the similarity of D3-4 between Zoraptera and Dictyoptera can serve as potentially valuable autapomorphy or synapomorphy in phylogeny reconstruction. The clades of (Plecoptera+Dermaptera) and ((Grylloblattodea+Mantophasmatodea)+(Embiodea+Phasmatodea)) were also recovered in the phylogenetic study. In addition, considering the other studies based on rDNAs, this study reached the highest congruence with previous phylogenetic studies of Holometabola based on nuclear protein coding genes or morphology characters. Future comparative studies of secondary structures across deep divergences and additional taxa are likely to reveal conserved patterns, structures and motifs that can provide support for major phylogenetic lineages.