Evolution and phylogenetic utility of alignment gaps within intron sequences of three nuclear genes in bumble bees (Bombus)

Using the combined gene approach and multiple analytical methods to improve the phylogeny and classification of Bombus (Hymenoptera, Apidae) in China

Bumble bees are vital to our agro-ecological system, with approximately 250 species reported around the world in the single genus Bombus. However, the health of bumble bees is threatened by multiple factors: habitat loss, climate change, pesticide use, and disease caused by pathogens and parasites. It is therefore vitally important to have a fully developed phylogeny for bumble bee species as part of our conservation efforts. The purpose of this study was to explore the phylogenetic relationships of the dominant bumble bees on the Tibetan plateau and in northern China as well as their placement and classification within the genus Bombus. The study used combined gene analysis consisting of sequence fragments from six genes, 16S rRNA, COI, EF-1α, Argk, Opsin and PEPCK, and the phylogenetic relationships of 209 Bombus species were explored. Twenty-six species, including 152 gene sequences, were collected from different regions throughout China, and 1037 gene sequences representing 183 species were obtained from GenBank or BOLD. The results suggest that the 209 analyzed species belong to fifteen subgenera and that most of the subgenera in Bombus are monophyletic, which is in accordance with conventional morphology-based classifications. The phylogenetic trees also show that nearly all subgenera easily fall into two distinct clades: short-faced and long-faced. The study is the first to investigate the phylogenetic placement of Bombus turneri (Richards), Bombus opulentus Smith, Bombus pyrosoma Morawitz, Bombus longipennis Friese, Bombus minshanensis Bischoff, and Bombus lantschouensis Vogt, all of which are widely distributed throughout different regions of China. The knowledge and understanding gained from the findings can provide a molecular basis to accurately classify Bombus in China and to define strategies to conserve biodiversity and promote pollinator populations.

The evolution of haploid chromosome numbers in Meliponini

It is thought that two evolutionary mechanisms gave rise to chromosomal variation in bees: the first one points to polyploidy as the main cause of chromosomal evolution, while the second, Minimum Interaction Theory (MIT), is more frequently used to explain chromosomal changes in Meliponini and suggests that centric fission is responsible for variations in karyotype. However, differences in chromosome number between Meliponini and its sister taxa and in the karyotype patterns of the Melipona genus cannot be explained by MIT, suggesting that other events were involved in chromosomal evolution. Thus, we assembled cytogenetical and molecular information to reconstruct an ancestral chromosome number for Meliponini and its sister group, Bombini, and propose a hypothesis to explain the evolutionary pathways underpinning chromosomal changes in Meliponini. We hypothesize that the common ancestor shared by the Meliponini and Bombini tribes possessed a chromosome number of n = 18. The karyotype with n = 17 chromosomes was maintained in Meliponini, and variations of haploid numbers possibly originated through additional Robertsonian fissions and fusions. Thus, the low chromosome number would not be an ancestral condition, as predicted by MIT. We then conclude that Robertsonian fission and fusions are unlikely to be the cause of chromosomal rearrangements that originated the current karyotypes in Meliponini.

Phylogeography of the Assassin Bug Sphedanolestes impressicollis in East Asia Inferred From Mitochondrial and Nuclear Gene Sequences

The assassin bug, Sphedanolestes impressicollis (Hemiptera: Reduviidae), is widely distributed in East Asia. It is an ideal model for evaluating the effects of climatic fluctuation and geographical events on the distribution patterns of East Asian reduviids. Here, we used two mitochondrial genes and one nuclear gene to investigate the phylogeographic pattern of the assassin bug based on comprehensive sampling in China, Japan, South Korea, Vietnam, and Laos. High levels of genetic differentiation were detected among the geographic populations classified into the northern and southern groups. A significant correlation was detected between genetic and geographical distances. The East China Sea land bridge served as a “dispersal corridor” during Pleistocene glaciation. The estimated divergence time indicated that the northern group may have separated from the eastern Chinese populations when the sea level rapidly rose during the “Ryukyu Coral Sea Stage” and the East China Sea land bridge was completely submerged. Demographic history and ecological niche modeling suggested that appropriate climatic conditions may have accounted for the rapid spread across the Korean Peninsula and Japan during the late Pleistocene. Our study underscores the pivotal roles of the Pleistocene sea level changes and climatic fluctuations in determining the distribution patterns of East Asian reduviids.

Invasion of the assassin bug Agriosphodrus dohrni (Hemiptera: Reduviidae) to Japan: Source estimation inferred from mitochondrial and nuclear gene sequences

A large-sized assassin bug Agriosphodrus dohrni (Signoret), has been recorded from India, Vietnam, China and Japan. It is one of the potential biological control agents against some important agricultural and forest pests. This species is speculated to have invaded Japan from its native range in China about 60 years ago. We used three mitochondrial gene fragments (COI, Cytb, and ND5) and one nuclear gene fragment (EF-1α) to clarify the invasion history of A. dohrni and assess the effects of geographic events and associated ecological adaptation on the distribution pattern. The native populations of A. dohrni in China are divided into three distinct groups, which might be molded by the Early Pleistocene glaciation event and diverged during the Calabrian Stage. However, consistent with the hypothesis of a recent invasion, extremely low level of genetic variation was detected in the Japanese populations, with only two haplotypes for the combined mitochondrial genes. Both the splits network and the ML/BI phylogenetic trees revealed that haplotypes of Japan were more closely-related to those from eastern China. Therefore, we postulate that there has been only one introduction event, probably from somewhere around the Nanjing (NJ) and Lin'an (LA) populations of eastern China.

Proteomic Characterization of the Venom of Five Bombus (Thoracobombus) Species

Venomous animals use venom, a complex biofluid composed of unique mixtures of proteins and peptides, to act on vital systems of the prey or predator. In bees, venom is solely used for defense against predators. However, the venom composition of bumble bees (Bombus sp.) is largely unknown. The Thoracobombus subgenus of Bombus sp. is a diverse subgenus represented by 14 members across Turkey. In this study, we sought out to proteomically characterize the venom of five Thoracobombus species by using bottom-up proteomic techniques. We have obtained two-dimensional polyacrylamide gel (2D-PAGE) images of each species’ venom sample. We have subsequently identified the protein spots by using matrix assisted laser desorption ionization/time-of-flight mass spectrometry (MALDI-TOF MS). We have identified 47 proteins for Bombus humilis, 32 for B. pascuorum, 60 for B. ruderarius, 39 for B. sylvarum, and 35 for B. zonatus. Moreover, we illustrated that intensities of 2DE protein spots corresponding to putative venom toxins vary in a species-specific manner. Our analyses provide the primary proteomic characterization of five bumble bee species’ venom composition.

Choosing and Using Introns in Molecular Phylogenetics

Introns are now commonly used in molecular phylogenetics in an attempt to recover gene trees that are concordant with species trees, but there are a range of genomic, logistical and analytical considerations that are infrequently discussed in empirical studies that utilize intron data. This review outlines expedient approaches for locus selection, overcoming paralogy problems, recombination detection methods and the identification and incorporation of LVHs in molecular systematics. A range of parsimony and Bayesian analytical approaches are also described in order to highlight the methods that can currently be employed to align sequences and treat indels in subsequent analyses. By covering the main points associated with the generation and analysis of intron data, this review aims to provide a comprehensive introduction to using introns (or any non-coding nuclear data partition) in contemporary phylogenetics.

Comparative phylogeography in the Atlantic forest and Brazilian savannas: pleistocene fluctuations and dispersal shape spatial patterns in two bumblebees

Background

Bombus morio and B. pauloensis are sympatric widespread bumblebee species that occupy two major Brazilian biomes, the Atlantic forest and the savannas of the Cerrado. Differences in dispersion capacity, which is greater in B. morio, likely influence their phylogeographic patterns. This study asks which processes best explain the patterns of genetic variation observed in B. morio and B. pauloensis, shedding light on the phenomena that shaped the range of local populations and the spatial distribution of intra-specific lineages.

Results

Results suggest that Pleistocene climatic oscillations directly influenced the population structure of both species. Correlative species distribution models predict that the warmer conditions of the Last Interglacial contributed to population contraction, while demographic expansion happened during the Last Glacial Maximum. These results are consistent with physiological data suggesting that bumblebees are well adapted to colder conditions. Intra-specific mitochondrial genealogies are not congruent between the two species, which may be explained by their documented differences in dispersal ability.

Conclusions

While populations of the high-dispersal B. morio are morphologically and genetically homogeneous across the species range, B. pauloensis encompasses multiple (three) mitochondrial lineages, and show clear genetic, geographic, and morphological differences. Because the lineages of B. pauloensis are currently exposed to distinct climatic conditions (and elevations), parapatric diversification may occur within this taxon. The eastern portion of the state of São Paulo, the most urbanized area in Brazil, represents the center of genetic diversity for B. pauloensis.

Electronic supplementary material

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

Distinguishing between weedy Amaranthus species based on intron 1 sequences from the 5-enolpyruvylshikimate-3-phosphate synthase gene

BACKGROUND

Hybridization between Amaranthus species and the potential for herbicide resistance to be transferred by hybridization are of growing concern in the weed science community. Early detection of evolved herbicide resistance and hybrids expressing resistance to single or multiple herbicides is important to develop an effective control strategy.

RESULTS

A PCR test was developed for quick identification of weedy amaranths and any hybrids. The sequences of intron 1 for the 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS; EC 2.5.1.19) gene were determined for Amaranthus palmeri, A. spinosus, A. retroflexus, A. blitoides, A. viridis, A. tuberculatus and A. hybridus. These sequences were aligned and primers were developed in areas where the sequence differed between species. Species-specific primers and cycle conditions were successfully developed. These primers produce a single robust band only for the species for which they were designed.

CONCLUSION

The PCR techniques described here allow identification of a weedy amaranth or suspect hybrid in a few hours. Using a similar target, it may be possible to design simple PCR tests to identify even more difficult to distinguish weed species or weeds prone to interspecific hybridization. Published 2016. This article is a U.S. Government work and is in the public domain in the USA.

The converse to Bergmann's rule in bumblebees, a phylogenetic approach

Two patterns commonly emerge when animal body size is analyzed as a function of latitudinal distribution. First, body size increases with latitude, a temperature effect known as Bergmann's rule, and second, the converse to Bergmann's rule, a pattern in which body size decreases with latitude. However, other geographic patterns can emerge when the mechanisms that generate Bergmann's and the converse to Bergmann's clines operate together. Here, we use phylogenetic comparative analysis in order to control for phylogenetic inertia, and we show that bumblebees exhibit the converse to Bergmann's rule. Bumblebee taxa are distributed worldwide in temperate and tropical regions. The largest species are found in places with high water availability during the driest time of the year. Nonetheless, large body size is constrained by extreme temperatures. Bumblebees’ body size could be related to a higher extent to the size of food rewards to be harvested than to the energetic advantages of thermoregulation. Moreover, we found that the body size of eusocial and cuckoo species responded in the same way to environmental variables, suggesting that they have not diverged due to different selective pressures.

High genetic diversity and geographic subdivision of three lance nematode species (Hoplolaimus spp.) in the United States

Lance nematodes (Hoplolaimus spp.) feed on the roots of a wide range of plants, some of which are agronomic crops. Morphometric values of amphimictic lance nematode species overlap considerably, and useful morphological characters for their discrimination require high magnification and significant diagnostic time. Given their morphological similarity, these Hoplolaimus species provide an interesting model to investigate hidden diversity in crop agroecosystems. In this scenario, H. galeatus may have been over-reported and the related species that are morphologically similar could be more widespread in the United States that has been recognized thus far. The main objectives of this study were to delimit Hoplolaimus galeatus and morphologically similar species using morphology, phylogeny, and a barcoding approach, and to estimate the genetic diversity and population structure of the species found. Molecular analyses were performed using sequences of the cytochrome c oxidase subunit 1 (Cox1) and the internal transcribed spacer (ITS1) on 23 populations. Four morphospecies were identified: H. galeatus, H. magnistylus, H. concaudajuvencus, and H. stephanus, along with a currently undescribed species. Pronounced genetic structure correlated with geographic origin was found for all species, except for H. galeatus. Hoplolaimus galeatus also exhibited low genetic diversity and the shortest genetic distances among populations. In contrast, H. stephanus, the species with the fewest reports from agricultural soils, was the most common and diverse species found. Results of this project may lead to better delimitation of lance nematode species in the United States by contributing to the understanding the diversity within this group.

Cladistic analysis of self-grooming indicates a single origin of eusociality in corbiculate bees (Hymenoptera: Apidae)

Behavioural traits have been used extensively in recent years as an important character source for making phylogenetic inferences. The phylogenetic positions of the members of the Apini subtribe are increasingly being debated, and new characters must be examined. We analysed the presence and absence of certain behavioural patterns, as well as the sequences of some of these patterns, to generate 79 characters. Eleven species comprised the ingroup, and Xylocopini comprised the outgroup. Parsimony analysis showed that the most parsimonious tree was (Euglossina(Bombina(Apina+Meliponina))). This topology is consistent with most studies that use morphological data and the few that use behavioural data, which suggests that advanced eusociality arose only once in a common ancestor of the clade Apina plus Meliponina; however, this hypothesis is inconsistent with our molecular data. Thus we considered behavioural, molecular, and morphological data and recovered the same topology, in which eusociality has a single origin in corbiculate bees.

Delphinid systematics and biogeography with a focus on the current genus Lagenorhynchus: multiple pathways for antitropical and trans-oceanic radiation

The six species currently classified within the genus Lagenorhynchus exhibit a pattern of antitropical distribution common among marine taxa. In spite of their morphological similarities they are now considered an artificial grouping, and include both recent and the oldest representatives of the Delphinidae radiation. They are, therefore, a good model for studying questions about the evolutionary processes that have driven dolphin speciation, dispersion and distribution. Here we used two different approaches. First we constructed a multigenic phylogeny with a minimum amount of missing data (based on 9 genes, 11,030bp, using the 6 species of the genus and their closest relatives) to infer their relationships. Second, we built a supermatrix phylogeny (based on 33 species and 27 genes) to test the effect of taxon sampling on the phylogeny of the genus, to provide inference on biogeographic history, and provide inference on the main events shaping the dispersion and radiation of delphinids. Our analyses suggested an early evolutionary history of marine dolphins in the North Atlantic Ocean and revealed multiple pathways of migration and radiation, probably guided by paleoceanographic changes during the Miocene and Pliocene. L. acutus and L. albirostris likely shared a common ancestor that arose in the North Atlantic around the Middle Miocene, predating the radiation of subfamilies Delphininae, Globicephalinae and Lissodelphininae.

A phylogenetic analysis of the genus Fragaria (strawberry) using intron-containing sequence from the ADH-1 gene

The genus Fragaria encompasses species at ploidy levels ranging from diploid to decaploid. The cultivated strawberry, Fragaria×ananassa, and its two immediate progenitors, F. chiloensis and F. virginiana, are octoploids. To elucidate the ancestries of these octoploid species, we performed a phylogenetic analysis using intron-containing sequences of the nuclear ADH-1 gene from 39 germplasm accessions representing nineteen Fragaria species and one outgroup species, Dasiphora fruticosa. All trees from Maximum Parsimony and Maximum Likelihood analyses showed two major clades, Clade A and Clade B. Each of the sampled octoploids contributed alleles to both major clades. All octoploid-derived alleles in Clade A clustered with alleles of diploid F. vesca, with the exception of one octoploid allele that clustered with the alleles of diploid F. mandshurica. All octoploid-derived alleles in clade B clustered with the alleles of only one diploid species, F. iinumae. When gaps encoded as binary characters were included in the Maximum Parsimony analysis, tree resolution was improved with the addition of six nodes, and the bootstrap support was generally higher, rising above the 50% threshold for an additional nine branches. These results, coupled with the congruence of the sequence data and the coded gap data, validate and encourage the employment of sequence sets containing gaps for phylogenetic analysis. Our phylogenetic conclusions, based upon sequence data from the ADH-1 gene located on F. vesca linkage group II, complement and generally agree with those obtained from analyses of protein-encoding genes GBSSI-2 and DHAR located on F. vesca linkage groups V and VII, respectively, but differ from a previous study that utilized rDNA sequences and did not detect the ancestral role of F. iinumae.

Phylogeny of suckermouth catfishes (Mochokidae: Chiloglanis) from Kenya: the utility of Growth Hormone introns in species level phylogenies

African suckermouth catfishes (Mochokidae: Chiloglanis) occur in freshwater throughout tropical Africa. Specimens from all major drainages across Kenya were collected over three field seasons. Here we present a phylogeny inferred from both mitochondrial cytochrome b (cyt b) and introns of the nuclear Growth Hormone gene (GH). The phylogeny inferred from introns is largely congruent with the results from an analysis of cyt b. The length and variability of GH introns make them ideal species level nuclear markers without the problem of introgression commonly encountered with mitochondrial genes. This analysis confirmed the presence of two previously known undescribed Chiloglanis species and also suggests the presence of previously unknown diversity within the Athi River system. The resulting phylogeny also indicates the presence of two separate lineages within C. brevibarbis. The historical biogeography of Chiloglanis within Kenya is discussed. The utility of GH intron for species level phylogenies of Siluriformes is compared to that in other groups.

A note on probabilistic models over strings: the linear algebra approach

Probabilistic models over strings have played a key role in developing methods that take into consideration indels as phylogenetically informative events. There is an extensive literature on using automata and transducers on phylogenies to do inference on these probabilistic models, in which an important theoretical question is the complexity of computing the normalization of a class of string-valued graphical models. This question has been investigated using tools from combinatorics, dynamic programming, and graph theory, and has practical applications in Bayesian phylogenetics. In this work, we revisit this theoretical question from a different point of view, based on linear algebra. The main contribution is a set of results based on this linear algebra view that facilitate the analysis and design of inference algorithms on string-valued graphical models. As an illustration, we use this method to give a new elementary proof of a known result on the complexity of inference on the "TKF91" model, a well-known probabilistic model over strings. Compared to previous work, our proving method is easier to extend to other models, since it relies on a novel weak condition, triangular transducers, which is easy to establish in practice. The linear algebra view provides a concise way of describing transducer algorithms and their compositions, opens the possibility of transferring fast linear algebra libraries (for example, based on GPUs), as well as low rank matrix approximation methods, to string-valued inference problems.

Convergent intron gains in hymenopteran elongation factor-1α

The eukaryotic translation elongation factor-1α gene (eEF1A) has been used extensively in higher level phylogenetics of insects and other groups, despite being present in two or more copies in several taxa. Orthology assessment has relied heavily on the position of introns, but the basic assumption of low rates of intron loss and absence of convergent intron gains has not been tested thoroughly. Here, we study the evolution of eEF1A based on a broad sample of taxa in the insect order Hymenoptera. The gene is universally present in two copies - F1 and F2 - both of which apparently originated before the emergence of the order. An elevated ratio of non-synonymous versus synonymous substitutions and differences in rates of amino acid replacements between the copies suggest that they evolve independently, and phylogenetic methods clearly cluster the copies separately. The F2 copy appears to be ancient; it is orthologous with the copy known as F1 in Diptera, and is likely present in most insect orders. The hymenopteran F1 copy, which may or may not be unique to this order, apparently originated through retroposition and was originally intron free. During the evolution of the Hymenoptera, it has successively accumulated introns, at least three of which have appeared at the same position as introns in the F2 copy or in eEF1A copies in other insects. The sites of convergent intron gain are characterized by highly conserved nucleotides that strongly resemble specific intron-associated sequence motifs, so-called proto-splice sites. The significant rate of convergent intron gain renders intron-exon structure unreliable as an indicator of orthology in eEF1A, and probably also in other protein-coding genes.

Re-mind the gap! Insertion - deletion data reveal neglected phylogenetic potential of the nuclear ribosomal internal transcribed spacer (ITS) of fungi

Rapidly evolving, indel-rich phylogenetic markers play a pivotal role in our understanding of the relationships at multiple levels of the tree of life. There is extensive evidence that indels provide conserved phylogenetic signal, however, the range of phylogenetic depths for which gaps retain tree signal has not been investigated in detail. Here we address this question using the fungal internal transcribed spacer (ITS), which is central in many phylogenetic studies, molecular ecology, detection and identification of pathogenic and non-pathogenic species. ITS is repeatedly criticized for indel-induced alignment problems and the lack of phylogenetic resolution above species level, although these have not been critically investigated. In this study, we examined whether the inclusion of gap characters in the analyses shifts the phylogenetic utility of ITS alignments towards earlier divergences. By re-analyzing 115 published fungal ITS alignments, we found that indels are slightly more conserved than nucleotide substitutions, and when included in phylogenetic analyses, improved the resolution and branch support of phylogenies across an array of taxonomic ranges and extended the resolving power of ITS towards earlier nodes of phylogenetic trees. Our results reconcile previous contradicting evidence for the effects of data exclusion: in the case of more sophisticated indel placement, the exclusion of indel-rich regions from the analyses results in a loss of tree resolution, whereas in the case of simpler alignment methods, the exclusion of gapped sites improves it. Although the empirical datasets do not provide to measure alignment accuracy objectively, our results for the ITS region are consistent with previous simulations studies alignment algorithms. We suggest that sophisticated alignment algorithms and the inclusion of indels make the ITS region and potentially other rapidly evolving indel-rich loci valuable sources of phylogenetic information, which can be exploited at multiple taxonomic levels.

Accelerated evolution of mitochondrial but not nuclear genomes of Hymenoptera: new evidence from crabronid wasps

Mitochondrial genes in animals are especially useful as molecular markers for the reconstruction of phylogenies among closely related taxa, due to the generally high substitution rates. Several insect orders, notably Hymenoptera and Phthiraptera, show exceptionally high rates of mitochondrial molecular evolution, which has been attributed to the parasitic lifestyle of current or ancestral members of these taxa. Parasitism has been hypothesized to entail frequent population bottlenecks that increase rates of molecular evolution by reducing the efficiency of purifying selection. This effect should result in elevated substitution rates of both nuclear and mitochondrial genes, but to date no extensive comparative study has tested this hypothesis in insects. Here we report the mitochondrial genome of a crabronid wasp, the European beewolf (Philanthus triangulum, Hymenoptera, Crabronidae), and we use it to compare evolutionary rates among the four largest holometabolous insect orders (Coleoptera, Diptera, Hymenoptera, Lepidoptera) based on phylogenies reconstructed with whole mitochondrial genomes as well as four single-copy nuclear genes (18S rRNA, arginine kinase, wingless, phosphoenolpyruvate carboxykinase). The mt-genome of P. triangulum is 16,029 bp in size with a mean A+T content of 83.6%, and it encodes the 37 genes typically found in arthropod mt genomes (13 protein-coding, 22 tRNA, and two rRNA genes). Five translocations of tRNA genes were discovered relative to the putative ancestral genome arrangement in insects, and the unusual start codon TTG was predicted for cox2. Phylogenetic analyses revealed significantly longer branches leading to the apocritan Hymenoptera as well as the Orussoidea, to a lesser extent the Cephoidea, and, possibly, the Tenthredinoidea than any of the other holometabolous insect orders for all mitochondrial but none of the four nuclear genes tested. Thus, our results suggest that the ancestral parasitic lifestyle of Apocrita is unlikely to be the major cause for the elevated substitution rates observed in hymenopteran mitochondrial genomes.

Host tracking or cryptic adaptation? Phylogeography of Pediobius saulius (Hymenoptera, Eulophidae), a parasitoid of the highly invasive horse-chestnut leafminer

Classical biological control is often advocated as a tool for managing invasive species. However, accurate evaluations of parasitoid species complexes and assessment of host specificity are impeded by the lack of morphological variation. Here, we study the possibility of host races/species within the eulophid wasp Pediobius saulius, a pupal generalist parasitoid that parasitize the highly invasive horse-chestnut leaf-mining moth Cameraria ohridella. We analysed the population genetic structure, host associations and phylogeographic patterns of P. saulius in Europe using the COI mitochondrial gene. This marker strongly supports a division into at least five highly differentiated parasitoid complexes, within two of which clades with differing degrees of host specialization were found: a Balkan clade that mainly (but not only) attacks C. ohridella and a more generalist European group that attacks many hosts, including C. ohridella. The divergence in COI (up to 7.6%) suggests the existence of cryptic species, although this is neither confirmed by nuclear divergence nor morphology. We do not find evidence of host tracking. The higher parasitism rates observed in the Balkans and the scarcity of the Balkan–Cameraria haplotypes out of the Balkans open the possibility of using these Balkan haplotypes as biological control agents of C. ohridella elsewhere in Europe.

Recent trends in molecular phylogenetic analysis: where to next?

The acquisition of large multilocus sequence data is providing researchers with an unprecedented amount of information to resolve difficult phylogenetic problems. With these large quantities of data comes the increasing challenge regarding the best methods of analysis. We review the current trends in molecular phylogenetic analysis, focusing specifically on the topics of multiple sequence alignment and methods of tree reconstruction. We suggest that traditional methods are inadequate for these highly heterogeneous data sets and that researchers employ newer more sophisticated search algorithms in their analyses. If we are to best extract the information present in these data sets, a sound understanding of basic phylogenetic principles combined with modern methodological techniques are necessary.

Diversity of phylogenetic information according to the locus and the taxonomic level: an example from a parasitic mesostigmatid mite genus

Molecular markers for cladistic analyses may perform differently according to the taxonomic group considered and the historical level under investigation. Here we evaluate the phylogenetic potential of five different markers for resolving evolutionary relationships within the ectoparasitic genus Dermanyssus at the species level, and their ability to address questions about the evolution of specialization. COI provided 9–18% divergence between species (up to 9% within species), 16S rRNA 10–16% (up to 4% within species), ITS1 and 2 2–9% (up to 1% within species) and Tropomyosin intron n 8–20% (up to 6% within species). EF-1α revealed different non-orthologous copies within individuals of Dermanyssus and Ornithonyssus. Tropomyosin intron n was shown containing consistent phylogenetic signal at the specific level within Dermanyssus and represents a promising marker for future prospects in phylogenetics of Acari. Phylogenetic analyses revealed that the generalist condition is apomorphic and D. gallinae might represent a complex of hybridized lineages. The split into hirsutus-group and gallinae-group in Dermanyssus does not seem to be appropriate based upon these results and D. longipes appears to be composed of two different entities.

Conservation phylogenetics of helodermatid lizards using multiple molecular markers and a supertree approach

We analyzed both mitochondrial (mt-) and nuclear (n) DNAs in a conservation phylogenetic framework to examine deep and shallow histories of the Beaded Lizard (Heloderma horridum) and Gila Monster (H. suspectum) throughout their geographic ranges in North and Central America. Both mtDNA and intron markers clearly partitioned each species. One intron and mtDNA further subdivided H. horridum into its four recognized subspecies (H. n. alvarezi, charlesbogerti,exasperatum, and horridum). However, the two subspecies of H. suspectum (H. s. suspectum and H. s. cinctum) were undefined. A supertree approach sustained these relationships. Overall, the Helodermatidae is reaffirmed as an ancient and conserved group. Its most recent common ancestor (MRCA) was Lower Eocene [35.4 million years ago (mya)], with a approximately 25 my period of stasis before the MRCA of H. horridum diversified in Lower Miocene. Another approximately 5 my passed before H. h. exasperatum and H. h. horridum diverged, followed by approximately 1.5 my before H. h. alvarezi and H. h. charlesbogerti separated. Heloderma suspectum reflects an even longer period of stasis (approximately 30 my) before diversifying from its MRCA. Both H. suspectum (México) and H. h. alvarezi also revealed evidence of historic range expansion following a recent bottleneck. Our conservation phylogenetic approach emphasizes the origin and diversification of this group, yields information on the manner by which past environmental variance may have impacted its populations and, in turn, allows us to disentangle historic from contemporary impacts that might threaten its long-term persistence. The value of helodermatid conservation resides in natural services and medicinal products, particularly venom constituents, and these are only now being realized.

A two-locus DNA sequence database for typing plant and human pathogens within the Fusarium oxysporum species complex

We constructed a two-locus database, comprising partial translation elongation factor (EF-1alpha) gene sequences and nearly full-length sequences of the nuclear ribosomal intergenic spacer region (IGS rDNA) for 850 isolates spanning the phylogenetic breadth of the Fusarium oxysporum species complex (FOSC). Of the 850 isolates typed, 101 EF-1alpha, 203 IGS rDNA, and 256 two-locus sequence types (STs) were differentiated. Analysis of the combined dataset suggests that two-thirds of the STs might be associated with a single host plant. This analysis also revealed that the 26 STs associated with human mycoses were genetically diverse, including several which appear to be nosocomial in origin. A congruence analysis, comparing partial EF-1alpha and IGS rDNA bootstrap consensus, identified a significant number of conflicting relationships dispersed throughout the bipartitions, suggesting that some of the IGS rDNA sequences may be non-orthologous. We also evaluated enniatin, fumonisin and moniliformin mycotoxin production in vitro within a phylogenetic framework.

Phylogeny and character evolution of the coprinoid mushroom genus Parasola as inferred from LSU and ITS nrDNA sequence data

Phylogenetic relationships, species concepts and morphological evolution of the coprinoid mushroom genus Parasola were studied. A combined dataset of nuclear ribosomal ITS and LSU sequences was used to infer phylogenetic relationships of Parasola species and several outgroup taxa. Clades recovered in the phylogenetic analyses corresponded well to morphologically discernable species, although in the case of P. leiocephala, P. lilatincta and P. plicatilis amended concepts proved necessary. Parasola galericuliformis and P. hemerobia are shown to be synonymous with P. leiocephala and P. plicatilis, respectively. By mapping morphological characters on the phylogeny, it is shown that the emergence of deliquescent Parasola taxa was accompanied by the development of pleurocystidia, brachybasidia and a plicate pileus. Spore shape and the position of the germ pore on the spores showed definite evolutionary trends within the group: from ellipsoid the former becomes more voluminous and heart-shaped, the latter evolves from central to eccentric in taxa referred to as 'crown' Parasola species. The results are discussed and compared to other Coprinus s.l. and Psathyrella taxa. Homoplasy and phylogenetic significance of various morphological characters, as well as indels in ITS and LSU sequences, are also evaluated.

Exon-primed intron-crossing (EPIC) PCR markers of Helicoverpa armigera (Lepidoptera: Noctuidae)

Applying microsatellite DNA markers in population genetic studies of the pest moth Helicoverpa armigera is subject to numerous technical problems, such as the high frequency of null alleles, occurrence of size homoplasy, presence of multiple copies of flanking sequence in the genome and the lack of PCR amplification robustness between populations. To overcome these difficulties, we developed exon-primed intron-crossing (EPIC) nuclear DNA markers for H. armigera based on ribosomal protein (Rp) and the Dopa Decarboxylase (DDC) genes and sequenced alleles showing length polymorphisms. Allele length polymorphisms were usually from random indels (insertions or deletions) within introns, although variation of short dinucleotide DNA repeat units was also detected. Mapping crosses demonstrated Mendelian inheritance patterns for these EPIC markers and the absence of both null alleles and allele 'dropouts'. Three examples of allele size homoplasies due to indels were detected in EPIC markers RpL3, RpS6 and DDC, while sequencing of multiple individuals across 11 randomly selected alleles did not detect indel size homoplasies. The robustness of the EPIC-PCR markers was demonstrated by PCR amplification in the related species, H. zea, H. assulta and H. punctigera.

Hosts are more important than destinations: What genetic variation in Microctonus aethiopoides (Hymenoptera: Braconidae) means for foreign exploration for natural enemies

Nucleotide sequence data were generated from the gene regions COI, 16S, and arginine kinase to assess genetic variation within the Palearctic parasitoid, Microctonus aethiopoides, reared from Sitona discoideus, S. hispidulus, and Hypera postica collected from two proximate locations in Mediterranean France. Partitioned Bayesian phylogenetic analyses of the molecular data provided strong support for the presence of at least two M. aethiopoides biotypes, one associated with Hypera species and the other with Sitona species. These new results combined with previously published data from 14 countries show that M. aethiopoides genetic variation is much more strongly correlated with host taxon than with sampling location. This contrasts with earlier perceptions that M. aethiopoides exhibits significant geographic variation, and helps to explain the widely varying biological control outcomes that have been obtained following the introductions of M. aethiopoides to Australia, New Zealand, and North America. The results strongly suggest that success rates and environmental safety in biological control would both be improved by ensuring that parasitoids collected in the native range are reared from the same host species as the one being targeted for control in the region of introduction. The results also provided insights both on the evolution of M. aethiopoides' host range, and on its evolutionary transition between solitary and gregarious larval development.

A molecular phylogenetic and fruit evolutionary analysis of the major groups of the paleotropical woody bamboos (Gramineae: Bambusoideae) based on nuclear ITS, GBSSI gene and plastid trnL-F DNA sequences

This study presented the first molecular phylogenetic analysis of the major clades of woody bamboos of the Old World tropics based on nuclear and chloroplast sequences (ITS, GBSSI and trnL-F). Sequence data from 53 species, representing 17 paleotropical woody bamboo genera, were analyzed using the maximum parsimony and Bayesian inference methods. All examined ingroup taxa were clustered into two clades, i.e., the Bambusinae+Dinochloa clade and the Melocanninae clade. The former clade included Bambusa, Bonia, Dendrocalamus, Dendrocalamopsis, Dinochloa, Gigantochloa, Molecalamus, Neomicrocalamus, Neosinocalamus, Oxytenanthera s. str. (sensu stricto), Racemobambos and Thyrsostachys. The Melocanninae clade consisted of Cephalostachyum, Leptocanna (better treated as part of Cephalostachyum), Melocanna, Pseudostachyum and Schizostachyum s. str. The subtribe Racemobambosinae and tribes Dendrocalameae and Oxytenanthereae were not supported and may be better placed in subtribe Bambusinae. The ovary characters seemed to be good criteria to distinguish these two clades. The reconstruction of ancestral fruit characters indicated that the bacoid caryopsis, namely, fleshy or berry-like fruits, was found to be scattered in three lineages of the examined paleotropical woody bamboos. Fruit characters are thus not reliable indicators of phylogeny and bacoid caryopsis likely represents a specialization for particular ecological conditions.

Evaluating nuclear protein-coding genes for phylogenetic utility in beetles

Although nuclear protein-coding genes have proven broadly useful for phylogenetic inference, relatively few such genes are regularly employed in studies of Coleoptera, the most diverse insect order. We increase the number of loci available for beetle systematics by developing protocols for three genes previously unused in beetles (alpha-spectrin, RNA polymerase II and topoisomerase I) and by refining protocols for five genes already in use (arginine kinase, CAD, enolase, PEPCK and wingless). We evaluate the phylogenetic performance of each gene in a Bayesian framework against a presumably known test phylogeny. The test phylogeny covers 31 beetle specimens and two outgroup taxa of varying age, including three of the four extant beetle suborders and a denser sampling in Adephaga and in the carabid genus Bembidion. All eight genes perform well for Cenozoic divergences and accurately separate closely related species within Bembidion, but individual genes differ markedly in accuracy over the older Mesozoic and Permian divergences. The concatenated data reconstruct the test phylogeny with high support in both Bayesian and parsimony analyses, indicating that combining data from multiple nuclear loci will be a fruitful approach for assembling the beetle tree of life.

Molecular phylogenetics and reproductive incompatibility in a complex of cryptic species of aphid parasitoids

We infer the phylogeny of a complex of cryptic species and populations of parasitic wasps and examine how reproductive incompatibility maps onto the molecular phylogeny. We used four nuclear (28S-D2, ITS1, ITS2, ArgK) and two mitochondrial (COI, COII) gene regions to analyze relationships among populations in the Aphelinus varipes species complex (Hymenoptera: Aphelinidae) from throughout Eurasia (France, Georgia, Israel, China, Korea and Japan) and from three aphid hosts (Aphis glycines, Diuraphis noxia and Rhopalosiphum padi; Hemiptera: Aphididae). A combined analysis of 21 genotypes of Aphelinus resulted in six most-parsimonious trees, and successive approximations character-weighting selected two of these as best supported by the data. All six gene regions were necessary to fully resolve the relationships among taxa. Four clades within the A. varipes complex were distinguished: (1) Aphelinus kurdjumovi, (2) Aphelinus hordei, (3) Aphelinus atriplicis, Aphelinus varipes, and Aphelinus albipodus, and (4) Aphelinus certus (populations from China, Korea, and Japan). Based on rates of nucleotide substitutions, these clades diverged between 78 and 526 thousand years ago during a period of repeated glaciations in Eurasia. In laboratory crosses, A. kurdjumovi, A. hordei, and A. varipes were reproductively incompatible with one another and all other populations. A. atriplicis was incompatible with these three species, but not with A. certus. The populations of A. certus from China, Japan, and Korea were reproductively compatible with one another but not with the other populations. Thus, with one exception, entities that were phylogenetically distinct were also reproductively incompatible with one another. Our evidence on molecular differentiation and reproductive incompatibility supports recognition of at least five cryptic species in the A. varipes complex. We discuss likely reasons for the high rate of speciation in this complex.

Low nuclear DNA variation supports a recent origin of Hawaiian Hylaeus bees (Hymenoptera: Colletidae)

Previous phylogenetic work on the Hawaiian bees of the genus Hylaeus, based on mitochondrial DNA and morphology, appeared to support a recent origin for the group, but support for the resulting tree was weak. Four nuclear genes with varying evolutionary rates -- arginine kinase, EF-1alpha, opsin, and wingless -- were sequenced for a reduced taxon set in an attempt to find one or more data set that would provide better support. All showed very low variation (<2%) in the ingroup. Comparison among genes revealed a much higher than expected rate of evolution in mtDNA, especially at first and second positions. While the data from the nuclear genes showed insufficient variation for phylogenetic analysis, the strong sequence similarity among the Hawaiian species supports the previous hypothesis of a recent origin for the group.

The genomic landscape of short insertion and deletion polymorphisms in the chicken (Gallus gallus) Genome: a high frequency of deletions in tandem duplicates

It is increasingly recognized that insertions and deletions (indels) are an important source of genetic as well as phenotypic divergence and diversity. We analyzed length polymorphisms identified through partial (0.25x) shotgun sequencing of three breeds of domestic chicken made by the International Chicken Polymorphism Map Consortium. A data set of 140,484 short indel polymorphisms in unique DNA was identified after filtering for microsatellite structures. There was a significant excess of tandem duplicates at indel sites, with deletions of a duplicate motif outnumbering the generation of duplicates through insertion. Indel density was lower in microchromosomes than in macrochromosomes, in the Z chromosome than in autosomes, and in 100 bp of upstream sequence, 5'-UTR, and first introns than in intergenic DNA and in other introns. Indel density was highly correlated with single nucleotide polymorphism (SNP) density. The mean density of indels in pairwise sequence comparisons was 1.9 x 10(-4) indel events/bp, approximately 5% the density of SNPs segregating in the chicken genome. The great majority of indels involved a limited number of nucleotides (median 1 bp), with A-rich motifs being overrepresented at indel sites. The overrepresentation of deletions at tandem duplicates indicates that replication slippage in duplicate sequences is a common mechanism behind indel mutation. The correlation between indel and SNP density indicates common effects of mutation and/or selection on the occurrence of indels and point mutations.

The relative performance of indel-coding methods in simulations

We used simulations to compare the performance of 10 approaches that have been used for treating unambiguously aligned gaps in phylogenetic analyses. We examined how these approaches perform under the ideal conditions of correct alignments, as well as how robust they are to errors caused by use of inferred alignments. Our results indicate that 5th-state coding dramatically outperformed all other coding methods, which in turn all outperformed treating gaps as missing data or excluding gapped positions. Simple indel coding (SIC) and modified complex indel coding (MCIC) performed about the same, and generally outperformed the other indel-coding methods. The high performance of 5th-state coding was largely found to be a weighting artifact. We suggest that MCIC-coded gap characters be scored for all unambiguously aligned gaps in parsimony-based molecular phylogenetic analyses. When the number of terminals sampled precludes the use of MCIC, SIC may be used as an effective substitute.

Nematode small subunit phylogeny correlates with alignment parameters

The number of nuclear small subunit (SSU) ribosomal RNA (rRNA) sequences for Nematoda has increased dramatically in recent years, and although their use in constructing phylogenies has also increased, relatively little attention has been given to their alignment. Here we examined the sensitivity of the nematode SSU data set to different alignment parameters and to the removal of alignment ambiguous regions. Ten alignments were created with CLUSTAL W using different sets of alignment parameters (10 full alignments), and each alignment was examined by eye and alignment ambiguous regions were removed (creating 10 reduced alignments). These alignment ambiguous regions were analyzed as a third type of data set, culled alignments. Maximum parsimony, neighbor-joining, and parsimony bootstrap analyses were performed. The resulting phylogenies were compared to each other by the symmetric difference distance tree comparison metric (SymD). The correlation of the phylogenies with the alignment parameters was tested by comparing matrices from SymD with corresponding matrices of Manhattan distances representing the alignment parameters. Differences among individual parsimony trees from the full alignments were frequently correlated with the differences among alignment parameters (580/1000 tests), as were trees from the culled alignments (403/1000 tests). Differences among individual parsimony trees from the reduced alignments were less frequently correlated with the differences among alignment parameters (230/1000 tests). Differences among majority-rule consensus trees (50%) from the parsimony analysis of the full alignments were significantly correlated with the differences among alignment parameters, whereas consensus trees from the reduced and culled analyses were not correlated with the alignment parameters. These patterns of correlation confirm that choice of alignment parameters has the potential to bias the resultant phylogenies for the nematode SSU data set, and suggest that the removal of alignment ambiguous regions reduces this effect. Finally, we discuss the implications of conservative phylogenetic hypotheses for Nematoda produced by exploring alignment space and removing alignment ambiguous regions for SSU rDNA.

Phylogenetic analyses of parasites in the new millennium

Phylogenetic analysis has changed greatly in the last decade, and the most important themes in that change are reviewed here. Sequence data have become the most common source of phylogenetic information. This means that explicit models for evolutionary processes have been developed in a likelihood context, which allow more realistic data analyses. These models are becoming increasingly complex, both for nucleotides and for amino acid sequences, and so all such models need to be quantitatively assessed for each data set, to find the most appropriate one for use in any particular tree-building analysis. Bayesian analysis has been developed for tree-building and is greatly increasing in popularity. This is because a good heuristic strategy exists, which allows large data sets to be analyzed with complex evolutionary models in a practical time. Perhaps the most disappointing aspect of tree interpretation is the ongoing confusion between rooted and unrooted trees, while the effect of taxon and character sampling is often overlooked when constructing a phylogeny (especially in parasitology). The review finishes with a detailed consideration of the analysis of a multi-gene data set for several dozen taxa of Cryptosporidium (Apicomplexa), illustrating many of the theoretical and practical points highlighted in the review.

Dynamic evolution of selenocysteine utilization in bacteria: a balance between selenoprotein loss and evolution of selenocysteine from redox active cysteine residues

Comparative genomics and evolutionary analyses to examine the dynamics of selenocysteine utilization in bacteria reveal a dynamic balance between selenoprotein origin and loss.

Background

Selenocysteine (Sec) is co-translationally inserted into protein in response to UGA codons. It occurs in oxidoreductase active sites and often is catalytically superior to cysteine (Cys). However, Sec is used very selectively in proteins and organisms. The wide distribution of Sec and its restricted use have not been explained.

Results

We conducted comparative genomics and phylogenetic analyses to examine dynamics of Sec decoding in bacteria at both selenium utilization trait and selenoproteome levels. These searches revealed that 21.5% of sequenced bacteria utilize Sec, their selenoproteomes have 1 to 31 selenoproteins, and selenoprotein-rich organisms are mostly Deltaproteobacteria or Firmicutes/Clostridia. Evolutionary histories of selenoproteins suggest that Cys-to-Sec replacement is a general trend for most selenoproteins. In contrast, only a small number of Sec-to-Cys replacements were detected, and these were mostly restricted to formate dehydrogenase and selenophosphate synthetase families. In addition, specific selenoprotein gene losses were observed in many sister genomes. Thus, the Sec/Cys replacements were mostly unidirectional, and increased utilization of Sec by existing protein families was counterbalanced by loss of selenoprotein genes or entire selenoproteomes. Lateral transfers of the Sec trait were an additional factor, and we describe the first example of selenoprotein gene transfer between archaea and bacteria. Finally, oxygen requirement and optimal growth temperature were identified as environmental factors that correlate with changes in Sec utilization.

Conclusion

Our data reveal a dynamic balance between selenoprotein origin and loss, and may account for the discrepancy between catalytic advantages provided by Sec and the observed low number of selenoprotein families and Sec-utilizing organisms.

How should gaps be treated in parsimony? A comparison of approaches using simulation

Simulation with indels was used to produce alignments where true site homologies in DNA sequences were known; the gaps from these datasets were removed and the sequences were then aligned to produce hypothesized alignments. Both alignments were then analyzed under three widely used methods of treating gaps during tree reconstruction under the maximum parsimony principle. With the true alignments, for many cases (82%), there was no difference in topological accuracy for the different methods of gap coding. However, in cases where a difference was present, coding gaps as a fifth state character or as separate presence/absence characters outperformed treating gaps as unknown/missing data nearly 90% of the time. For the hypothesized alignments, on average, all gap treatment approaches performed equally well. Data sets with higher sequence divergence and more pectinate tree shapes with variable branch lengths are more affected by gap coding than datasets associated with shallower non-pectinate tree shapes.