How do insect nuclear and mitochondrial gene substitution patterns differ? Insights from Bayesian analyses of combined datasets

The impact of molecular data on our understanding of bee phylogeny and evolution

Our understanding of bee phylogeny has improved over the past fifteen years as a result of new data, primarily nucleotide sequence data, and new methods, primarily model-based methods of phylogeny reconstruction. Phylogenetic studies based on single or, more commonly, multilocus data sets have helped resolve the placement of bees within the superfamily Apoidea; the relationships among the seven families of bees; and the relationships among bee subfamilies, tribes, genera, and species. In addition, molecular phylogenies have played an important role in inferring evolutionary patterns and processes in bees. Phylogenies have provided the comparative framework for understanding the evolution of host-plant associations and pollen specialization, the evolution of social behavior, and the evolution of parasitism. In this paper, we present an overview of significant discoveries in bee phylogeny based primarily on the application of molecular data. We review the phylogenetic hypotheses family-by-family and then describe how the new phylogenetic insights have altered our understanding of bee biology.

Molecular-based identification and phylogeny of Oligonychus species (Acari: Tetranychidae)

The genus Oligonychus has been morphologically divided into two groups based on the direction of curvature of the aedeagus and includes some morphologically similar species that are difficult to distinguish. To develop DNA-based methods for identifying Oligonychus species and to determine the phylogenetic relationships among them, we examined the cytochrome c oxidase subunit I gene of mitochondrial DNA and the internal transcribed spacer and 28S regions of nuclear ribosomal RNA gene for 17 species. Based on the genetic distances (p-distances) of the three DNA regions, the range of intraspecific divergence was found to be below (and not overlap) the range of interspecific divergence, which allowed the 17 species to be discriminated correctly, consistent with their classification based on morphology. Phylogenetic trees constructed by neighbor-joining and Bayesian methods clearly showed two clades, consisting of species whose aedeagi curve ventrally and dorsally, respectively. Three Oligonychus species inhabiting gramineous plants formed clearly defined subclades.

Molecular phylogenetics of a South Pacific sap beetle species complex (Carpophilus spp., Coleoptera: Nitidulidae)

Several species of sap beetles in the genus Carpophilus are minor pests of fresh produce and stored products, and are frequently intercepted in biosecurity operations. In the South Pacific region, the superficially similar species C. maculatus and C. oculatus are frequently encountered in these situations. Three subspecies of C. oculatus have been described, and the complex of these four taxa has led to inaccurate identification and questions regarding the validity of these taxa. A molecular phylogenetic study using the mitochondrial gene cytochrome c oxidase I (COI) and two nuclear markers comprising the rDNA internal transcribed spacer 2 (ITS2) and the D1-D2 region of the large (28S) ribosomal RNA subunit showed that C. maculatus, and C. o. cheesmani were easily differentiated from the two other subspecies of C. oculatus. COI also showed differentiation between C. o. gilloglyi and C. o. oculatus, but this was not shown when third codon positions were removed and when RY-coding analyses were conducted. Generalised mixed Yule-Coalescent (GMYC) models were fitted to trees estimated from the COI data and were analysed using a multimodel approach to consider the evidence for three taxonomic groupings of the C. oculatus group. While the arrangement with the highest cumulative weight was not the arrangement ultimately accepted, the accepted taxonomy also had an acceptable level of support. ITS2 showed structure within C. oculatus, however C. o. oculatus was resolved as paraphyletic with respect to C. o. gilloglyi. COI showed evidence of sequence saturation and did not adequately resolve higher relationships between species represented in the dataset. 28S resolved higher relationships, but did not perform well at the species level. This study supports the validity of C. maculatus as a separate species, and provides sufficient evidence to raise C. o. cheesmani to the level of species. This study also shows significant structure within and between C. o. gilloglyi and C. o. oculatus, giving an indication of recent speciation events occurring. To highlight the interesting biology between these two taxa, C. o. gilloglyi is retained as a subspecies of C. oculatus. These results give clarity regarding the taxonomic status of C. maculatus and the subspecies of C. oculatus and provide a platform for future systematic research on Carpophilus.

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.

The complete mitochondrial genome and novel gene arrangement of the unique-headed bug Stenopirates sp. (Hemiptera: Enicocephalidae)

Many of true bugs are important insect pests to cultivated crops and some are important vectors of human diseases, but few cladistic analyses have addressed relationships among the seven infraorders of Heteroptera. The Enicocephalomorpha and Nepomorpha are consider the basal groups of Heteroptera, but the basal-most lineage remains unresolved. Here we report the mitochondrial genome of the unique-headed bug Stenopirates sp., the first mitochondrial genome sequenced from Enicocephalomorpha. The Stenopirates sp. mitochondrial genome is a typical circular DNA molecule of 15, 384 bp in length, and contains 37 genes and a large non-coding fragment. The gene order differs substantially from other known insect mitochondrial genomes, with rearrangements of both tRNA genes and protein-coding genes. The overall AT content (82.5%) of Stenopirates sp. is the highest among all the known heteropteran mitochondrial genomes. The strand bias is consistent with other true bugs with negative GC-skew and positive AT-skew for the J-strand. The heteropteran mitochondrial atp8 exhibits the highest evolutionary rate, whereas cox1 appears to have the lowest rate. Furthermore, a negative correlation was observed between the variation of nucleotide substitutions and the GC content of each protein-coding gene. A microsatellite was identified in the putative control region. Finally, phylogenetic reconstruction suggests that Enicocephalomorpha is the sister group to all the remaining Heteroptera.

What is the phylogenetic signal limit from mitogenomes? The reconciliation between mitochondrial and nuclear data in the Insecta class phylogeny

BACKGROUND

Efforts to solve higher-level evolutionary relationships within the class Insecta by using mitochondrial genomic data are hindered due to fast sequence evolution of several groups, most notably Hymenoptera, Strepsiptera, Phthiraptera, Hemiptera and Thysanoptera. Accelerated rates of substitution on their sequences have been shown to have negative consequences in phylogenetic inference. In this study, we tested several methodological approaches to recover phylogenetic signal from whole mitochondrial genomes. As a model, we used two classical problems in insect phylogenetics: The relationships within Paraneoptera and within Holometabola. Moreover, we assessed the mitochondrial phylogenetic signal limits in the deeper Eumetabola dataset, and we studied the contribution of individual genes.

RESULTS

Long-branch attraction (LBA) artefacts were detected in all the datasets. Methods using Bayesian inference outperformed maximum likelihood approaches, and LBA was avoided in Paraneoptera and Holometabola when using protein sequences and the site-heterogeneous mixture model CAT. The better performance of this method was evidenced by resulting topologies matching generally accepted hypotheses based on nuclear and/or morphological data, and was confirmed by cross-validation and simulation analyses. Using the CAT model, the order Strepsiptera was recovered as sister to Coleoptera for the first time using mitochondrial sequences, in agreement with recent results based on large nuclear and morphological datasets. Also the Hymenoptera-Mecopterida association was obtained, leaving Coleoptera and Strepsiptera as the basal groups of the holometabolan insects, which coincides with one of the two main competing hypotheses. For the Paraneroptera, the currently accepted non-monophyly of Homoptera was documented as a phylogenetic novelty for mitochondrial data. However, results were not satisfactory when exploring the entire Eumetabola, revealing the limits of the phylogenetic signal that can be extracted from Insecta mitogenomes. Based on the combined use of the five best topology-performing genes we obtained comparable results to whole mitogenomes, highlighting the important role of data quality.

CONCLUSION

We show for the first time that mitogenomic data agrees with nuclear and morphological data for several of the most controversial insect evolutionary relationships, adding a new independent source of evidence to study relationships among insect orders. We propose that deeper divergences cannot be inferred with the current available methods due to sequence saturation and compositional bias inconsistencies. Our exploratory analysis indicates that the CAT model is the best dealing with LBA and it could be useful for other groups and datasets with similar phylogenetic difficulties.

Pitfalls in comparisons of genetic distances: a case study of the avian family Acrocephalidae

Genetic distances are increasingly being used for identification and species delimitation, especially since the introduction of "barcoding". While for phylogenetic inferences great care is generally taken to choose the best-fit evolutionary model, this is usually neglected in calculating genetic distances. Moreover, distances obtained from others than best-fit models, different lengths of sequences, and even different loci are often freely compared. We examined the influence of different methods on calculating genetic distances using mitochondrial cytochrome b sequences for the passerine family Acrocephalidae. We found substantial differences between: (1) corrected distances based on the best-fit model (TrN+Γ) vs. uncorrected p-distances; (2) distances calculated based on different parts of the same gene; and (3) distances calculated using the methods of "complete deletion" vs. "pairwise deletion" for sequences that included uncertain nucleotides. All these methodological differences affected comparisons between species and potential taxonomical conclusions. We suggest that (1) different loci are incomparable. (2) Only perfectly homologous regions (same length, same part of locus) should be compared. (3) In the case of sequences with some uncertain nucleotides, only distances calculated by the method of "complete deletion" are fully comparable. (4) Only distances based on the optimal substitution model should be used. (5) Even within the same locus, corrected genetic distances are unique to the study in which they are calculated, as they are conditional on the particular dataset and model selected for that dataset.

Lymnaea schirazensis, an overlooked snail distorting fascioliasis data: genotype, phenotype, ecology, worldwide spread, susceptibility, applicability

BACKGROUND

Lymnaeid snails transmit medical and veterinary important trematodiases, mainly fascioliasis. Vector specificity of fasciolid parasites defines disease distribution and characteristics. Different lymnaeid species appear linked to different transmission and epidemiological patterns. Pronounced susceptibility differences to absolute resistance have been described among lymnaeid populations. When assessing disease characteristics in different endemic areas, unexpected results were obtained in studies on lymnaeid susceptibility to Fasciola. We undertook studies to understand this disease transmission heterogeneity.

METHODOLOGY/PRINCIPAL FINDINGS

A ten-year study in Iran, Egypt, Spain, the Dominican Republic, Mexico, Venezuela, Ecuador and Peru, demonstrated that such heterogeneity is not due to susceptibility differences, but to a hitherto overlooked cryptic species, Lymnaea schirazensis, confused with the main vector Galba truncatula and/or other Galba/Fossaria vectors. Nuclear rDNA and mtDNA sequences and phylogenetic reconstruction highlighted an old evolutionary divergence from other Galba/Fossaria species, and a low intraspecific variability suggesting a recent spread from one geographical source. Morphometry, anatomy and egg cluster analyses allowed for phenotypic differentiation. Selfing, egg laying, and habitat characteristics indicated a migration capacity by passive transport. Studies showed that it is not a vector species (n = 8572 field collected, 20 populations): snail finding and penetration by F. hepatica miracidium occur but never lead to cercarial production (n = 338 experimentally infected).

CONCLUSIONS/SIGNIFICANCE

This species has been distorting fasciolid specificity/susceptibility and fascioliasis geographical distribution data. Hence, a large body of literature on G. truncatula should be revised. Its existence has henceforth to be considered in research. Genetic data on livestock, archeology and history along the 10,000-year post-domestication period explain its wide spread from the Neolithic Fertile Crescent. It is an efficient biomarker for the follow-up of livestock movements, a crucial aspect in fascioliasis emergence. It offers an outstanding laboratory model for genetic studies on susceptibility/resistance in F. hepatica/lymnaeid interaction, a field of applied research with disease control perspectives.

The molecular systematics of blowflies and screwworm flies (Diptera: Calliphoridae) using 28S rRNA, COX1 and EF-1α: insights into the evolution of dipteran parasitism

The Calliphoridae include some of the most economically significant myiasis-causing flies in the world - blowflies and screwworm flies - with many being notorious for their parasitism of livestock. However, despite more than 50 years of research, key taxonomic relationships within the family remain unresolved. This study utilizes nucleotide sequence data from the protein-coding genes COX1 (mitochondrial) and EF1α (nuclear), and the 28S rRNA (nuclear) gene, from 57 blowfly taxa to improve resolution of key evolutionary relationships within the family Calliphoridae. Bayesian phylogenetic inference was carried out for each single-gene data set, demonstrating significant topological difference between the three gene trees. Nevertheless, all gene trees supported a Calliphorinae-Luciliinae subfamily sister-lineage, with respect to Chrysomyinae. In addition, this study also elucidates the taxonomic and evolutionary status of several less well-studied groups, including the genus Bengalia (either within Calliphoridae or as a separate sister-family), genus Onesia (as a sister-genera to, or sub-genera within, Calliphora), genus Dyscritomyia and Lucilia bufonivora, a specialised parasite of frogs and toads. The occurrence of cross-species hybridisation within Calliphoridae is also further explored, focusing on the two economically significant species Lucilia cuprina and Lucilia sericata. In summary, this study represents the most comprehensive molecular phylogenetic analysis of family Calliphoridae undertaken to date.

DNA sequence characterisation and phylogeography of Lymnaea cousini and related species, vectors of fascioliasis in northern Andean countries, with description of L. meridensis n. sp. (Gastropoda: Lymnaeidae)

Background

Livestock fascioliasis is a problem throughout Ecuador, Colombia and Venezuela, mainly in Andean areas where the disease also appears to affect humans. Transmission patterns and epidemiological scenarios of liver fluke infection have shown to differ according to the lymnaeid vector snail species involved. These Andean countries present the vectors Lymnaea cousini, L. bogotensis and L. ubaquensis, unknown in the rest of Latin America. An exhaustive combined haplotype study of these species is performed by means of DNA sequencing of the nuclear ribosomal 18S RNA gene, ITS-2 and ITS-1, and mitochondrial DNA cox1 gene.

Results

The conserved 5.8S rDNA sequence corroborated that no pseudogenes are involved in the numerous non-microsatellite/minisatellite-related indels appearing between the ITS-2 and ITS-1 sequences when comparing different L. cousini - L. bogotensis populations. Sequence analyses and phylogenetic reconstruction methods including other lymnaeid vector species show that (i) L. bogotensis is a synonym of L. cousini, (ii) L. ubaquensis is a synonym of Pseudosuccinea columella, and (iii) populations of L. cousini hitherto known from Venezuelan highlands indeed belong to a new species for which the name L. meridensis n. sp. is proposed. This new species is described and a complete phenotypic differentiation provided.

Conclusions

ITS-2, ITS-1 and cox1 prove to be good markers for specimen classification and haplotype characterisation of these morphologically similar lymnaeids in endemic areas. Analysis of the 18S gene and phylogenetic reconstructions indicate that L. cousini and L. meridensis n. sp. cluster in an evolutionary line different from the one of P. columella, despite their external resemblance. This suggests an evolutionary phenotypic convergence related to similar environments and which has given rise to frequent specimen misclassification. Body size and phylogenetic relationships of L. meridensis n. sp. with well-known vectors as Lymnaea cousini and P. columella, as well as with Galba/Fossaria species, suggest that the new species may participate in disease transmission to both animals and humans in altitude areas during the yearly window in which temperatures are higher than the F. hepatica minimum development threshold. The involvement of L. cousini and P. columella in the transmission and geographical/altitudinal distribution of fascioliasis in these Andean countries is analysed.

The architecture and complete sequence of mitochondrial genome of an assassin bug Agriosphodrus dohrni (Hemiptera: Reduviidae)

The 16, 470 bp nucleotide sequence of the mitochondrial genome (mitogenome) of an assassin bug from the reduviid subfamily Harpactorinae, Agriosphodrus dohrni, has been revealed. The entire genome encodes for two ribosomal RNA genes (rrnL and rrnS), 22 transfer RNA (tRNA) genes, 13 protein-coding genes, and a control region. The nucleotide composition is biased toward adenine and thymine (A+T = 72.2%). Comparative analysis with two other reduviid species Triatoma dimidiata and Valentia hoffmanni, exhibited highly conserved genome architectures including genome contents, gene order, nucleotide composition, codon usage, amino acid composition, as well as genome asymmetry. All protein-coding genes use standard mitochondrial initiation codons (methionine and isoleucine), except that nad1 starts with GTG. All tRNAs have the classic clover-leaf structure, except that the dihydrouridine (DHU) arm of tRNA(Ser(AGN)) forms a simple loop. Secondary structure comparisons of the two mitochondrial ribosomal subunits among sequenced assassin bugs show that the sequence and structure of rrnL is more conservative than that of rrnS. The presence of structural elements in the control region is also discussed, with emphasis on their implications in the regulation of replication and/or transcription of the reduviid mitogenome. The phylogenetic analyses indicated that within Reduviidae, Harpactorinae is a sister group to the Salyavatinae + Triatominae clade.

Molecular systematics and species limits in the Philippine fantails (Aves: Rhipidura)

Islands have long-attracted scientists because of their relatively simple biotas and stark geographic boundaries. However, for many islands and archipelagos, this simplicity may be overstated because of methodological and conceptual limitations when these biotas were described. One archipelago that has received relatively little recent attention is the Philippine islands. Although much of its biota was documented long ago, taxonomic revision and evolutionary study has been surprisingly scarce, and only a few molecular phylogenetic studies are beginning to appear. We present a molecular phylogeny and taxonomic revision for the Philippine fantails (Aves: Rhipidura) using nuclear and mitochondrial DNA sequences. Our results suggest that current taxonomy underestimates diversity in the group. Some morphologically distinct subspecies warrant species status, whereas one was indistinguishable genetically and morphologically and should not be retained. A few taxa require additional sampling for thorough taxonomic assessment. Patterns of diversity within Philippine Rhipidura mostly corroborate predictions of the Pleistocene aggregate island complex (PAIC) hypothesis, in which diversity is expected to be partitioned by deep water channels separating Pleistocene aggregate islands rather than by current islands. Substantial structure within PAIC clades indicates that additional drivers of diversification should be considered.

Taxonomic and phylogenetic relationships between species of the genus Culex (Diptera: culicidae) from Brazil inferred from the cytochrome c oxidase I mitochondrial gene

Species of the genus Culex Linnaeus have been incriminated as the main vectors of lymphatic filariases and are important vectors of arboviruses, including West Nile virus. Sequences corresponding to a fragment of 478 bp of the cytochrome c oxidase subunit I gene, which includes part of the barcode region, of 37 individuals of 17 species of genus Culex were generated to establish relationships among five subgenera, Culex, Phenacomyia, Melanoconion, Microculex, and Carrollia, and one species of the genus Lutzia that occurs in Brazil. Bayesian methods were employed for the phylogenetic analyses. Results of sequence comparisons showed that individuals identified as Culex dolosus, Culex mollis, and Culex imitator possess high intraspecific divergence (3.1, 2.3, and 3.5%, respectively) when using the Kimura two parameters model. These differences were associated either with distinct morphological characteristics of the male genitalia or larval and pupal stages, suggesting that these may represent species complexes. The Bayesian topology suggested that the genus and subgenus Culex are paraphyletic relative to Lutzia and Phenacomyia, respectively. The cytochrome c oxidase subunit I sequences may be a useful tool to both estimate phylogenetic relationships and identify morphologically similar species of the genus Culex.

Phylogenetic relationships of true butterflies (Lepidoptera: Papilionoidea) inferred from COI, 16S rRNA and EF-1α sequences

The molecular phylogenetic relationships among true butterfly families (superfamily Papilionoidea) have been a matter of substantial controversy; this debate has led to several competing hypotheses. Two of the most compelling of those hypotheses involve the relationships of (Nymphalidae + Lycaenidae) + (Pieridae + Papilionidae) and (((Nymphalidae + Lycaenidae) + Pieridae) + Papilionidae). In this study, approximately 3,500 nucleotide sequences from cytochrome oxidase subunit I (COI), 16S ribosomal RNA (16S rRNA), and elongation factor-1 alpha (EF-1α) were sequenced from 83 species belonging to four true butterfly families, along with those of three outgroup species belonging to three lepidopteran superfamilies. These sequences were subjected to phylogenetic reconstruction via Bayesian Inference (BI), Maximum Likelihood (ML), and Maximum Parsimony (MP) algorithms. The monophyletic Pieridae and monophyletic Papilionidae evidenced good recovery in all analyses, but in some analyses, the monophylies of the Lycaenidae and Nymphalidae were hampered by the inclusion of single species of the lycaenid subfamily Miletinae and the nymphalid subfamily Danainae. Excluding those singletons, all phylogenetic analyses among the four true butterfly families clearly identified the Nymphalidae as the sister to the Lycaenidae and identified this group as a sister to the Pieridae, with the Papilionidae identified as the most basal linage to the true butterfly, thus supporting the hypothesis: (Papilionidae + (Pieridae + (Nymphalidae + Lycaenidae))).

Avtzis D, Riegler M, Stauffer C. Mitochondrial phylogenies in the light of pseudogenes and Wolbachia: re-assessment of a bark beetle dataset

Phylogenetic studies based on mtDNA become increasingly questioned because of potential pitfalls due to mitochondrial pseudogenes and mitochondrial selective sweeps. While the inclusion of nuclear markers should preferentially be considered for future studies, there is no need to abandon mtDNA as long as tests for the known mtDNA artefacts are performed. In this study we presentadditionaldata and test previous phylogeographical studies of Pityogenes chalcographus. We did not detect nuclear copies (numts) of the previously used mitochondrial markers by performing a combined long range/nested PCR of the COI gene and by an in silico analysis of the COI sequence data. This confirms the robustness of our previous phylogenetic study of Pityogenes chalcographus. Results of an in-situ hybridization of Wolbachia in Pityogenes chalcographus confirm the presence of this endosysmbiont in this species. However, we did not detect a correlation between infection status, geographical region and mtDNA haplotypes. The hybridisation data also support a previous hypothesis that infections do not result from parasitoids or parasitic nematodes, insect surface or laboratory contaminations and are hence a true infection of Pityogenes chalcographus. We conclude that the deep structure found in mitochondrial populations of Pityogenes chalcographus indeed represents the evolutionary history of European populations.

The complete mitochondrial genome and phylogenomics of a damselfly, Euphaea formosa support a basal Odonata within the Pterygota

This study determined the first complete mitochondrial genome of a damselfly, Euphaea formosa (Insecta: Odonata: Zygoptera), and reconstructed a phylogeny based on thirteen protein-coding genes of mitochondrial genomes in twenty-five representative hexapods to examine the relationships among the basal Pterygota. The damselfly's mitochondrial genome is a circular molecule of 15,700bp long, and contains the entire set of thirty-seven genes typically found in insects. The gene arrangement, nucleotide composition, and codon usage pattern of the mitochondrial genome are similar across the three odonate species, suggesting a conserved genome evolution within the Odonata. The presence of the intergenic spacer s5 likely represents a synapomorphy for the dragonflies (Anisoptera). Maximum parsimony, maximum likelihood, and Bayesian analyses of both nucleotide and amino acid sequences cannot support the three existing phylogenetic hypotheses of the basal Pterygota (Palaeoptera, Metapterygota, and Chiastomyaria). In contrast, the phylogenetic results indicate an alternative hypothesis of a strongly supported basal Odonata and a sister relationship of the Ephemeroptera and Plecoptera. The unexpected sister Ephemeroptera+Plecoptera clade, which contradicts with the widely accepted hypothesis of a monophyletic Neoptera, requires further analyses with additional mitochondrial genome sampling at the base of the Neoptera.

Placement of Conopidae (Diptera) within Schizophora based on mtDNA and nrDNA gene regions

The first attempt to phylogenetically place Conopidae using molecular characters, as well as the largest molecular analysis of relationships within Schizophora (Diptera) to date, is presented. Twenty-eight taxa from 11 acalyptrate families and seven acalyptrate superfamilies are represented. Nearly 12,800 bp of sequence data from 10 genes representing both mitochondrial (cytochrome oxidase I (COI), cytochrome b (cytB), and 12S) and nuclear genes (28S, the carbamoyl phosphate synthetase region of CAD (CAD), elongation factor-1alpha (EF-1alpha), white, alanyl-tRNA synthetase (AATS), triose phosphate isomerase (TPI), and phosphogluconate dehydrogenase (PGD)) are analysed. Parsimony and Bayesian analyses strongly support the monophyly of both Conopidae and Schizophora. While in the parsimony analysis, Conopidae are placed as sister to the remaining Schizophora, the Bayesian analysis recovers a Conopidae+Lauxaniidae clade. The value of nuclear, mitochondrial, ribosomal, and protein-coding gene sequence data for answering phylogenetic questions at different levels of divergence is evaluated.

Phylogenetic studies of sinipercid fish (Perciformes: Sinipercidae) based on multiple genes, with first application of an immune-related gene, the virus-induced protein (viperin) gene

The sinipercid fish represent a group of 12 species of freshwater percoid fish endemic to East Asia. To date published morphological and molecular phylogenetics hypotheses of sinipercid fish are part congruent, and there are some areas of significant disagreement with respect to species relationships. The present study used separate and combined methods to analyze 7307 bp of data from three mitochondrial genes (cyt b, CO1 and 16S rRNA; approximately 2312 bp) and three nuclear genes (viperin, the first two introns of S7 ribosomal protein gene; approximately 4995 bp) for the attempts to estimate the relationships among sinipercids and to assess the phylogenetic utility of these markers. Phylogenetic trees were reconstructed using maximum parsimony, maximum likelihood and partitioned Bayesian analyses. Despite the detection of significant heterogeneity of phylogenetic signal between the mitochondrial and nuclear partitions, the combined data analysis represented the best-supported topology of all data. The sinipercid fish form a monophyletic group with two distinct clades, one corresponding to the genus Siniperca and the other to Coreoperca. Coreoperca whiteheadi is the sister taxon to Coreoperca herzi plus Coreoperca kawamebari. In the Siniperca, Siniperca undulata is the sister taxon to the other members of Siniperca, within the subclade containing the other members of the genus, Siniperca chuatsi and Siniperca kneri are sister species, next joined by Siniperca obscura, Siniperca roulei, Sinipercascherzeri and finally by Siniperca fortis. The potential utilities of six different genes for phylogenetic resolution of closely related sinipercid species were also evaluated, with special interest in that of the novel virus-induced protein (viperin) gene.

An evaluation of phylogenetic informativeness profiles and the molecular phylogeny of diplazontinae (Hymenoptera, Ichneumonidae)

How to quantify the phylogenetic information content of a data set is a longstanding question in phylogenetics, influencing both the assessment of data quality in completed studies and the planning of future phylogenetic projects. Recently, a method has been developed that profiles the phylogenetic informativeness (PI) of a data set through time by linking its site-specific rates of change to its power to resolve relationships at different timescales. Here, we evaluate the performance of this method in the case of 2 standard genetic markers for phylogenetic reconstruction, 28S ribosomal RNA and cytochrome oxidase subunit 1 (CO1) mitochondrial DNA, with maximum parsimony, maximum likelihood, and Bayesian analyses of relationships within a group of parasitoid wasps (Hymenoptera: Ichneumonidae, Diplazontinae). Retrieving PI profiles of the 2 genes from our own and from 3 additional data sets, we find that the method repeatedly overestimates the performance of the more quickly evolving CO1 compared with 28S. We explore possible reasons for this bias, including phylogenetic uncertainty, violation of the molecular clock assumption, model misspecification, and nonstationary nucleotide composition. As none of these provides a sufficient explanation of the observed discrepancy, we use simulated data sets, based on an idealized setting, to show that the optimum evolutionary rate decreases with increasing number of taxa. We suggest that this relationship could explain why the formula derived from the 4-taxon case overrates the performance of higher versus lower rates of evolution in our case and that caution should be taken when the method is applied to data sets including more than 4 taxa.

Integrative taxonomy: a multisource approach to exploring biodiversity

Good alpha taxonomy is central to biology. On the basis of a survey of arthropod studies that used multiple disciplines for species delimitation, we evaluated the performance of single disciplines. All included disciplines had a considerable failure rate. Rigor in species delimitation can thus be increased when several disciplines chosen for complementarity are used. We present a flexible procedure and stopping rule for integrative taxonomy that uses the information from different disciplines separately. Disagreement among disciplines over the number and demarcation of species is resolved by elucidating and invoking evolutionary explanations for disagreement. With the identification of further promising study organisms and of new questions for in-depth analysis, evolutionary biology should profit from integrative taxonomy. An important rationale is clarity in researcher bias in the decision-making process. The success of integrative taxonomy will further increase through methodological progress, taxonomic training of evolutionary biologists, and balanced resource allocation.

Molecular evolution of the pDo500 satellite DNA family in Dolichopoda cave crickets (Rhaphidophoridae)

BACKGROUND

Non-coding satellite DNA (satDNA) usually has a high turn-over rate frequently leading to species specific patterns. However, some satDNA families evolve more slowly and can be found in several related species. Here, we analyzed the mode of evolution of the pDo500 satDNA family of Dolichopoda cave crickets. In addition, we discuss the potential of slowly evolving satDNAs as phylogenetic markers.

RESULTS

We sequenced 199 genomic or PCR amplified satDNA repeats of the pDo500 family from 12 Dolichopoda species. For the 38 populations under study, 39 pDo500 consensus sequences were deduced. Phylogenetic analyses using Bayesian, Maximum Parsimony, and Maximum Likelihood approaches yielded largely congruent tree topologies. The vast majority of pDo500 sequences grouped according to species designation. Scatter plots and statistical tests revealed a significant correlation between genetic distances for satDNA and mitochondrial DNA. Sliding window analyses showed species specific patterns of variable and conserved regions. The evolutionary rate of the pDo500 satDNA was estimated to be 1.63-1.78% per lineage per million years.

CONCLUSIONS

The pDo500 satDNA evolves gradually at a rate that is only slightly faster than previously published rates of insect mitochondrial COI sequences. The pDo500 phylogeny was basically congruent with the previously published mtDNA phylogenies. Accordingly, the slowly evolving pDo500 satDNA family is indeed informative as a phylogenetic marker.

ITS2 secondary structure improves phylogeny estimation in a radiation of blue butterflies of the subgenus Agrodiaetus (Lepidoptera: Lycaenidae: Polyommatus )

BACKGROUND

Current molecular phylogenetic studies of Lepidoptera and most other arthropods are predominantly based on mitochondrial genes and a limited number of nuclear genes. The nuclear genes, however, generally do not provide sufficient information for young radiations. ITS2 , which has proven to be an excellent nuclear marker for similarly aged radiations in other organisms like fungi and plants, is only rarely used for phylogeny estimation in arthropods, although universal primers exist. This is partly due to difficulties in the alignment of ITS2 sequences in more distant taxa. The present study uses ITS2 secondary structure information to elucidate the phylogeny of a species-rich young radiation of arthropods, the butterfly subgenus Agrodiaetus. One aim is to evaluate the efficiency of ITS2 to resolve the phylogeny of the subgenus in comparison with COI , the most important mitochondrial marker in arthropods. Furthermore, we assess the use of compensatory base changes in ITS2 for the delimitation of species and discuss the prospects of ITS2 as a nuclear marker for barcoding studies.

RESULTS

In the butterfly family Lycaenidae, ITS2 secondary structure enabled us to successfully align sequences of different subtribes in Polyommatini and produce a Profile Neighbour Joining tree of this tribe, the resolution of which is comparable to phylogenetic trees obtained with COI+COII . The subgenus Agrodiaetus comprises 6 major clades which are in agreement with COI analyses. A dispersal-vicariance analysis (DIVA) traced the origin of most Agrodiaetus clades to separate biogeographical areas in the region encompassing Eastern Anatolia, Transcaucasia and Iran.

CONCLUSIONS

With the inclusion of secondary structure information, ITS2 appears to be a suitable nuclear marker to infer the phylogeny of young radiations, as well as more distantly related genera within a diverse arthropod family. Its phylogenetic signal is comparable to the mitochondrial marker COI . Compensatory base changes are very rare within Polyommatini and cannot be used for species delimitation. The implementation of secondary structure information into character-based phylogenetic methods is suggested to further improve the versatility of this marker in phylogenetic studies.

Phylogenetic analysis and temporal diversification of mosquitoes (Diptera: Culicidae) based on nuclear genes and morphology

Background

Phylogenetic analyses provide a framework for examining the evolution of morphological and molecular diversity, interpreting patterns in biogeography, and achieving a stable classification. The generic and suprageneric relationships within mosquitoes (Diptera: Culicidae) are poorly resolved, making these subjects difficult to address.

Results

We carried out maximum parsimony and maximum likelihood, including Bayesian, analyses on a data set consisting of six nuclear genes and 80 morphological characters to assess their ability to resolve relationships among 25 genera. We also estimated divergence times based on sequence data and fossil calibration points, using Bayesian relaxed clock methods. Strong support was recovered for the basal position and monophyly of the subfamily Anophelinae and the tribes Aedini and Sabethini of subfamily Culicinae. Divergence times for major culicid lineages date to the early Cretaceous.

Conclusions

Deeper relationships within the family remain poorly resolved, suggesting the need for additional taxonomic sampling. Our results support the notion of rapid radiations early in the diversification of mosquitoes.

Toward reconstructing the evolution of advanced moths and butterflies (Lepidoptera: Ditrysia): an initial molecular study

BACKGROUND

In the mega-diverse insect order Lepidoptera (butterflies and moths; 165,000 described species), deeper relationships are little understood within the clade Ditrysia, to which 98% of the species belong. To begin addressing this problem, we tested the ability of five protein-coding nuclear genes (6.7 kb total), and character subsets therein, to resolve relationships among 123 species representing 27 (of 33) superfamilies and 55 (of 100) families of Ditrysia under maximum likelihood analysis.

RESULTS

Our trees show broad concordance with previous morphological hypotheses of ditrysian phylogeny, although most relationships among superfamilies are weakly supported. There are also notable surprises, such as a consistently closer relationship of Pyraloidea than of butterflies to most Macrolepidoptera. Monophyly is significantly rejected by one or more character sets for the putative clades Macrolepidoptera as currently defined (P < 0.05) and Macrolepidoptera excluding Noctuoidea and Bombycoidea sensu lato (P < or = 0.005), and nearly so for the superfamily Drepanoidea as currently defined (P < 0.08). Superfamilies are typically recovered or nearly so, but usually without strong support. Relationships within superfamilies and families, however, are often robustly resolved. We provide some of the first strong molecular evidence on deeper splits within Pyraloidea, Tortricoidea, Geometroidea, Noctuoidea and others.Separate analyses of mostly synonymous versus non-synonymous character sets revealed notable differences (though not strong conflict), including a marked influence of compositional heterogeneity on apparent signal in the third codon position (nt3). As available model partitioning methods cannot correct for this variation, we assessed overall phylogeny resolution through separate examination of trees from each character set. Exploration of "tree space" with GARLI, using grid computing, showed that hundreds of searches are typically needed to find the best-feasible phylogeny estimate for these data.

CONCLUSION

Our results (a) corroborate the broad outlines of the current working phylogenetic hypothesis for Ditrysia, (b) demonstrate that some prominent features of that hypothesis, including the position of the butterflies, need revision, and (c) resolve the majority of family and subfamily relationships within superfamilies as thus far sampled. Much further gene and taxon sampling will be needed, however, to strongly resolve individual deeper nodes.

Cryptic failure of partitioned Bayesian phylogenetic analyses: lost in the land of long trees

Partitioned Bayesian phylogenetic analyses of routine genetic data sets, constructed using MrBayes (Ronquist and Huelsenbeck 2003), can become trapped in regions of parameter space characterized by unrealistically long trees and distorted partition rate multipliers. Such analyses commonly fail to reach stationarity during hundreds of millions of generations of sampling-many times longer than most published analyses. Some data sets are so prone to this problem that paired MrBayes runs begun from different starting trees repeatedly find the same incorrect long-tree solutions and consequently pass the most commonly employed tests of stationarity, including the average standard deviation of split frequencies (ASDSF) and the potential scale reduction factor (PSRF) statistics offered by MrBayes (Gelman and Rubin 1992). In these situations, failure to reach stationarity is recognizable only in light of prior knowledge of model parameters, such as the expectation that third-codon-position sites usually evolve fastest in protein-coding genes. The conditions that lead to the long-tree problem are frequently encountered in phylogenetic studies today, and I present 6 demonstration examples from the literature. Although the effects on tree length (TL) are often dramatic, effects on topology appear to be subtle. Susceptibility to the problem is sometimes predicted by the difference between the true TL and the starting TL. In some cases, the problems described here can be avoided or reduced by manipulation of the starting TL and/or by adjustments to the prior on branch lengths. In more difficult situations, accurate branch length estimation may not be possible with Bayesian methods because of dependence of the solution on the branch length prior.

Molecular population genetics and phylogeography of the Chagas disease vector Triatoma infestans in South America

Knowledge of the genetic variability, population structure, and evolutionary history of Triatoma infestans may be useful for developing rational vector control strategies. A 661-bp fragment of the mitochondrial gene cytochrome oxidase I (COI) was sequenced and analyzed in bugs from Argentina, Uruguay, Peru, and Bolivia, including peridomestic, domestic, Andean, and Chaco sylvatic bugs. A total of 48 polymorphic sites among 37 haplotypes were described. Nucleotide variation fluctuated among samples, with the highest nucleotide diversity observed in seven Argentinean provinces. Within this group, some populations showed patterns of variability compatible with population expansions and/or fine-scale population structure, whereas others suggested population bottlenecks and/or population admixture processes. A maximum parsimony analysis of the haplotypes showed the presence of a Bolivian/Peruvian and an Argentinean/Uruguayan clade. Bolivian sequences were further divided in Chaco sylvatic and Andean domestic and sylvatic. Two different nested clades were found within the Argentinean/Uruguayan cluster. Analysis of molecular variance (AMOVA) and K(ST)* analysis supported a strong population structure in Argentina, where genetic differentiation was correlated with geographic distance. Departures from neutrality expectations and a nested cladistic analysis suggest a recent population expansion of T. infestans in Argentina, followed by restricted gene flow and patterns of isolation by distance. This expansion could have taken place as a two-wave process, as was shown by the phylogenetic analysis and signatures of population admixture in the southern most Argentinean populations.

Populations, hybrids and the systematic concepts of species and subspecies in Chagas disease triatomine vectors inferred from nuclear ribosomal and mitochondrial DNA

In Chagas disease, triatomine vectors are the main target for control measures because of the absence of effective drugs. The broad usefulness of nuclear rDNA and mtDNA sequences explains why triatomine studies using these markers have increased so pronouncedly in recent years. This indicates the appropriateness of an updated review about these molecular markers, concentrating on aspects useful for research on Chagas disease vectors. A comparative analysis is presented on the efficiency, weight of their different characteristics, limitations and problems of each of the different DNA markers in the light of the results obtained in studies on populations, hybrids, subspecies and species of the subfamily Triatominae. The use of a standardized composite haplotype code nomenclature for both nuclear rDNA and mtDNA markers is strongly encouraged to avoid difficulties in comparative studies. Triatomine aspects related to concerted evolution, microsatellites, minisatellites and insertions/deletions in nuclear rDNA and silent/non-silent mutations, pseudogenes and weaknesses of partial sequences in mtDNA are analysed. Introgression and hybrids, nuclear and mitochondrial DNA strengths, and compared evolutionary rates of nuclear rDNA and mtDNA in triatomines are discussed. Many conclusions are obtained thanks to the availability, for the first time in triatomines, of a complete sequence of a protein-coding mtDNA gene as ND1 from very numerous triatomine species covering from different populations of a species up to members belonging to different tribes. The evolutionary rates of each nuclear rDNA marker and mtDNA marker are analysed by comparison at subspecies level (intrapopulational, interpopulational, between morphs, and between subspecies) and species level (close and distant species of the same genus, species of different genera, and species of different tribes). Weaknesses of mtDNA for systematic-taxonomic purposes detected recently and newly in insects and triatomines, respectively, are discussed in detail. Emphasis is given to taxonomic units and biological entities presenting well-known problematics, both from the systematic-taxonomic and/or epidemiological-control points of view, as well as to molecular situations which can give rise to erroneous conclusions. All these aspects constitute the background on which the key question about the systematic concepts of species and subspecies in triatomines is focused. The global purpose is to facilitate future work on triatomines by highlighting present gaps, how better choice the appropriate markers, and marker aspects which should be taken into account. Key characteristics as alpha, CI and transformation rate matrices ought to be obtained and noted to get appropriate results and allow correct interpretations. The main aim is to offer a baseline for future fundamental research on triatomines and applied research on transmission, epidemiology and control measures related to Chagas disease vectors.

Resolving genetic diversity in Australasian Culex mosquitoes: incongruence between the mitochondrial cytochrome c oxidase I and nuclear acetylcholine esterase 2

Insects that vector pathogens are under constant surveillance in Australasia although the repertoire of genetic markers to distinguish what are often cryptic mosquito species remains limited. We present a comparative assessment of the second exon-intron region of the acetylcholine esterase 2 gene (ace-2) and the mitochondrial DNA cytochrome c oxidase I (COI) using two closely related Australasia mosquitoes Culex annulirostris and Culex palpalis. The COI revealed eight divergent lineages of which four were confirmed with the ace-2. We dissect out the nuclear chromosomal haplotypes of the ace-2 as well as the exon-intron regions by assessing the protein's tertiary structure to reveal a hypervariable 5'-exon that forms part of an external protein loop and displays a higher polymorphic rate than the intron. We retrace the evolutionary history of these mosquitoes by phylogenetic inference and by testing different evolutionary hypotheses. We conclude that DNA barcoding using COI may overestimate the diversity of Culex mosquitoes in Australasia and should be applied cautiously with support from the nuclear DNA such as the ace-2. Together the COI and ace-2 provide robust evidence for distinct cryptic Culex lineages--one of which correlates exactly with the southern limit of Japanese encephalitis virus activity in Australasia.

Phylogeny of colletid bees (Hymenoptera: Colletidae) inferred from four nuclear genes

Colletidae comprise approximately 2500 species of bees primarily distributed in the southern continents (only two colletid genera are widely distributed: Colletes and Hylaeus). Previously published studies have failed to resolve phylogenetic relationships on a worldwide basis and this has been a major barrier to the progress of research regarding systematics and evolution of colletid bees. For this study, data from four nuclear gene loci: elongation factor-1alpha (F2 copy), opsin, wingless, and 28S rRNA were analyzed for 122 species of colletid bees, representing all subfamilies and tribes currently recognized; 22 species belonging to three other bee families were used as outgroups. Bayesian, maximum likelihood, and parsimony methods were employed to investigate the phylogenetic relationships within Colletidae and resulted in highly congruent and well-resolved trees. The phylogenetic results show that Colletidae are monophyletic and that all traditionally recognized subfamilies (except Paracolletinae) are also strongly supported as monophyletic. Our phylogenetic hypothesis provides a framework within which broad questions related to the taxonomy, biogeography, morphology, evolution, and ecology of colletid bees can be addressed.

Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

Mitochondrial DNA (mtDNA) variants are widely used in evolutionary genetics as markers for population history and to estimate divergence times among taxa. Inferences of species history are generally based on phylogenetic comparisons, which assume that molecular evolution is clock-like. Between-species comparisons have also been used to estimate the mutation rate, using sites that are thought to evolve neutrally. We directly estimated the mtDNA mutation rate by scanning the mitochondrial genome of Drosophila melanogaster lines that had undergone approximately 200 generations of spontaneous mutation accumulation (MA). We detected a total of 28 point mutations and eight insertion-deletion (indel) mutations, yielding an estimate for the single-nucleotide mutation rate of 6.2 × 10⁻⁸ per site per fly generation. Most mutations were heteroplasmic within a line, and their frequency distribution suggests that the effective number of mitochondrial genomes transmitted per female per generation is about 30. We observed repeated occurrences of some indel mutations, suggesting that indel mutational hotspots are common. Among the point mutations, there is a large excess of G→A mutations on the major strand (the sense strand for the majority of mitochondrial genes). These mutations tend to occur at nonsynonymous sites of protein-coding genes, and they are expected to be deleterious, so do not become fixed between species. The overall mtDNA mutation rate per base pair per fly generation in Drosophila is estimated to be about 10× higher than the nuclear mutation rate, but the mitochondrial major strand G→A mutation rate is about 70× higher than the nuclear rate. Silent sites are substantially more strongly biased towards A and T than nonsynonymous sites, consistent with the extreme mutation bias towards A+T. Strand-asymmetric mutation bias, coupled with selection to maintain specific nonsynonymous bases, therefore provides an explanation for the extreme base composition of the mitochondrial genome of Drosophila.

Mitochondria are the energy-producing organelles of the cell, and they contain genetic information encoded on their own genome. Because rates of mutation for mitochondrial genomes are believed to be much higher than those in nuclear DNA, mitochondrial genetic differences between and within species are particularly useful in population genetics, for example, as markers of population movements. We have directly estimated the mutation rate in the mitochondrial genome of the fruit fly Drosophila melanogaster in lines that had been allowed to randomly accumulate mutations in the virtual absence of effective natural selection. We scanned for new mutations by comparing the DNA of different lines by a sensitive mutation detection technique. We show that the mitochondrial mutation rate is about ten times higher than the nuclear DNA mutation rate. Strikingly, however, almost all of the single–base pair mutations that we detected change G to A at an amino acid site of a protein-coding gene. The explanation for this effect seems to be that natural selection maintains the nucleotide G at amino acid sites, whereas most silent sites are under weaker selection and have previously mutated to A or T. The mutation rate for G to A changes is 70 times higher than the nuclear DNA mutation rate. This extreme mutation bias maintains the high A+T content of the Drosophila mitochondrial genome.

We show that the mitochondrial DNA mutation rate in D. melanogaster is 10X higher than the nuclear mutation rate, and that major-strand G-->A hypermutablity explains the extremely biased base composition of the mitochondrial genome.

Early Xanthochorema (Trichoptera, Insecta) radiations in New Caledonia originated on ultrabasic rocks

The toxic and nutrient poor ultrabasic rock substrate covering one-third of New Caledonia greatly influenced on the biogeography and diversity of plants in the island. Studies on the effect of ultrabasic substrate on fauna are almost entirely absent. In this paper we examine whether the diversification of Trichoptera of the New Caledonian endemic genus Xanthochorema Kimmins, 1953 was related to the presence of ultrabasic substrate. The analysis is based on data from a phylogeny derived from DNA sequences of mitochondrial COX1, COX2 and 16S, and nuclear EF1a genes. The study of the relationships between ancestral species and substrate was carried out using dispersal-vicariance analysis and tracing the history of substrate association with ultrabasic and non-ultrabasic distributions representing the terminals in the fully resolved phylogenetic tree. Our results show that (1) the ancestor of all Xanthochorema species was present on ultrabasic substrate, (2) early speciation events were restricted to ultrabasic substrate, (3) younger ancestral species dispersed into non-ultrabasic substrates, and (4) late speciation events were restricted to non-ultrabasic substrate. These results correspond to the hypothesis that New Caledonia once was more extensively covered by ultrabasic rocks than at present.

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.

Phylogeography and genetic variation of Triatoma dimidiata, the main Chagas disease vector in Central America, and its position within the genus Triatoma

Background

Among Chagas disease triatomine vectors, the largest genus, Triatoma, includes species of high public health interest. Triatoma dimidiata, the main vector throughout Central America and up to Ecuador, presents extensive phenotypic, genotypic, and behavioral diversity in sylvatic, peridomestic and domestic habitats, and non-domiciliated populations acting as reinfestation sources. DNA sequence analyses, phylogenetic reconstruction methods, and genetic variation approaches are combined to investigate the haplotype profiling, genetic polymorphism, phylogeography, and evolutionary trends of T. dimidiata and its closest relatives within Triatoma. This is the largest interpopulational analysis performed on a triatomine species so far.

Methodology and Findings

Triatomines from Mexico, Guatemala, Honduras, Nicaragua, Panama, Cuba, Colombia, Ecuador, and Brazil were used. Triatoma dimidiata populations follow different evolutionary divergences in which geographical isolation appears to have had an important influence. A southern Mexican–northern Guatemalan ancestral form gave rise to two main clades. One clade remained confined to the Yucatan peninsula and northern parts of Chiapas State, Guatemala, and Honduras, with extant descendants deserving specific status. Within the second clade, extant subspecies diversity was shaped by adaptive radiation derived from Guatemalan ancestral populations. Central American populations correspond to subspecies T. d. dimidiata. A southern spread into Panama and Colombia gave the T. d. capitata forms, and a northwestern spread rising from Guatemala into Mexico gave the T. d. maculipennis forms. Triatoma hegneri appears as a subspecific insular form.

Conclusions

The comparison with very numerous Triatoma species allows us to reach highly supported conclusions not only about T. dimidiata, but also on different, important Triatoma species groupings and their evolution. The very large intraspecific genetic variability found in T. dimidiata sensu lato has never been detected in a triatomine species before. The distinction between the five different taxa furnishes a new frame for future analyses of the different vector transmission capacities and epidemiological characteristics of Chagas disease. Results indicate that T. dimidiata will offer problems for control, although dwelling insecticide spraying might be successful against introduced populations in Ecuador.

Chagas disease is a serious parasitic disease of Latin America. Human contamination in poor rural or periurban areas is mainly attributed to haematophagous triatomine insects. Triatoma includes important vector species, as T. dimidiata in Central and Meso-America. DNA sequences, phylogenetic methods and genetic variation analyses are combined in a large interpopulational approach to investigate T. dimidiata and its closest relatives within Triatoma. The phylogeography of Triatoma indicates two colonization lineages northward and southward of the Panama isthmus during ancient periods, with T. dimidiata presenting a large genetic variability related to evolutionary divergences from a Mexican-Guatemalan origin. One clade remained confined to Yucatan, Chiapas, Guatemala and Honduras, with extant descendants deserving species status: T. sp. aff. dimidiata. The second clade gave rise to four subspecies: T. d. dimidiata in Guatemala and Mexico (Chiapas) up to Honduras, Nicaragua, Providencia island, and introduced into Ecuador; T. d. capitata in Panama and Colombia; T. d. maculipennis in Mexico and Guatemala; and T. d. hegneri in Cozumel island. This taxa distinction may facilitate the understanding of the diversity of vectors formerly included under T. dimidiata, their different transmission capacities and the disease epidemiology. Triatoma dimidiata will offer more problems for control than T. infestans in Uruguay, Chile and Brazil, although populations in Ecuador are appropriate targets for insecticide-spraying.

Revealing pancrustacean relationships: phylogenetic analysis of ribosomal protein genes places Collembola (springtails) in a monophyletic Hexapoda and reinforces the discrepancy between mitochondrial and nuclear DNA markers

BACKGROUND

In recent years, several new hypotheses on phylogenetic relations among arthropods have been proposed on the basis of DNA sequences. One of the challenged hypotheses is the monophyly of hexapods. This discussion originated from analyses based on mitochondrial DNA datasets that, due to an unusual positioning of Collembola, suggested that the hexapod body plan evolved at least twice. Here, we re-evaluate the position of Collembola using ribosomal protein gene sequences.

RESULTS

In total 48 ribosomal proteins were obtained for the collembolan Folsomia candida. These 48 sequences were aligned with sequence data on 35 other ecdysozoans. Each ribosomal protein gene was available for 25% to 86% of the taxa. However, the total sequence information was unequally distributed over the taxa and ranged between 4% and 100%. A concatenated dataset was constructed (5034 inferred amino acids in length), of which ~66% of the positions were filled. Phylogenetic tree reconstructions, using Maximum Likelihood, Maximum Parsimony, and Bayesian methods, resulted in a topology that supports monophyly of Hexapoda.

CONCLUSION

Although ribosomal proteins in general may not evolve independently, they once more appear highly valuable for phylogenetic reconstruction. Our analyses clearly suggest that Hexapoda is monophyletic. This underpins the inconsistency between nuclear and mitochondrial datasets when analyzing pancrustacean relationships. Caution is needed when applying mitochondrial markers in deep phylogeny.

Patterns and processes underlying evolutionary significant units in the Platypleura stridula L. species complex (Hemiptera: Cicadidae) in the Cape Floristic Region, South Africa

Cicadas have been shown to be useful organisms for examining the effects of distribution, plant association and geographical barriers on gene flow between populations. The cicadas of the Platypleura stridula species complex are restricted to the biologically diverse Cape Floristic Region (CFR) of South Africa. They are thus an excellent study group for elucidating the mechanisms by which hemipteran diversity is generated and maintained in the CFR. Phylogeographical analysis of this species complex using mitochondrial DNA Cytochrome Oxidase I (COI) and ribosomal 16S sequence data, coupled with preliminary morphological and acoustic data, resolves six clades, each of which has specific host-plant associations and distinct geographical ranges. The phylogeographical structure implies simultaneous or near-simultaneous radiation events, coupled with shifts in host-plant associations. When calibrated using published COI and 16S substitution rates typical for related insects, these lineages date back to the late Pliocene - early Pleistocene, coincident with vegetation change, altered drainage patterns and accelerated erosion in response to neotectonic crustal uplift and cyclic Pleistocene climate change, and glaciation-associated changes in climate and sea level.

A phylogeny of long-tongued horse flies (Diptera:Tabanidae:Philoliche) with the first cladistic review of higher relationships within the family

Long-tongued horse flies (Diptera : Tabanidae : Pangoniinae) have proboscis lengths at least as long as their heads, the longest belonging to the Old World genus Philoliche (Wiedemann, 1820). These long proboscides are used to probe for nectar in deep-throated flowers. For some flower species, these flies are the only known pollinators. Although horse flies are both vectors of disease and important pollinators, there has been no previous study of phylogenetic relationships below the subfamily level. The first comprehensive phylogenetic hypothesis for Philoliche (and consequently, Tabanidae) based on molecular data from one nuclear (CAD) and one mitochondrial (COI) gene is presented here. Using an exemplar approach, species from eight of nine tribes in all three subfamilies in Tabanidae were sampled, giving some of the first insights into relationships among the family as a whole. All nine subgenera of Philoliche, and multiple geographic representatives of the subgenus P. (Philoliche) (Wiedemann, 1820) in southern Africa were also sampled. Within the subgenus Philoliche, molecular support was found for a previously synonymised species. In general, these analyses sustain earlier intuitive classifications, but do not support the monophyly of all currently recognised subfamilies.

Nuclear mitochondrial-like sequences in ants: evidence from Atta cephalotes (Formicidae: Attini)

Nuclear mitochondrial-like sequences (numts) are copies of mitochondrial DNA that have migrated to the genomic DNA. We present the first characterization of numts in ants, these numts being homologues to a mitochondrial DNA fragment containing loci the 3' portion of the cytochrome oxidase I gene, an intergenic spacer, the tRNA leucine gene and the 5' portion of the cytochrome oxidase II gene. All 67 specimens of Atta cephalotes (Hymenoptera: Formicidae: Attini) investigated had these homologues, which are within two monophyletic groups that we called numt1 and numt2. Numt1 and numt2 sequences are less variable than mitochondrial sequences and released from the severe purifying selection constraining the evolution of mitochondrial genes. Their formation probably involved bottlenecks related to two distinct transfer events of ancient and fast evolving mitochondrial DNA fragments to comparative slowly evolving nuclear DNA regions.