Phylogeny,landmark analysis and the use of wing venation to study the evolution of social wasps (Hymenoptera: Vespidae: Vespinae)

Wing venation provides useful characters with which to classify extant and fossil insects. Recently, quantification of its shape using landmarks has increased the potential of wing venation to distinguish taxa. However, the use of wing landmarks in phylogenetic analyses remains largely unexplored. Here, we tested landmark analysis under parsimony (LAUP) to include wing shape data in a phylogenetic analysis of hornets and yellow jackets. Using 68 morphological characters, nine genes and wing landmarks, we produced the first total‐evidence phylogeny of Vespinae. We also tested the influence of LAUP parameters using simulated landmarks. Our data confirmed that optimization parameters, alignment method, landmark number and, under low optimization parameters, the initial orientation of aligned shapes can influence LAUP results. Furthermore, single landmark configurations never accurately reflected the topology used for data simulation, but results were significantly close when compared to random topologies. Thus, wing landmark configurations were unreliable phylogenetic characters when treated independently, but provided some useful insights when combined with other data. Our phylogeny corroborated the monophyly of most groups proposed on the basis of morphology and showed the fossil Palaeovespa is distantly related to extant genera. Unstable relationships among genera suggest that rapid radiations occurred in the early history of the Vespinae.     

The use of DNA barcodes to estimate phylogenetic diversity in forest communities of southern China

To elucidate potential ecological and evolutionary processes associated with the assembly of plant communities, there is now widespread use of estimates of phylogenetic diversity that are based on a variety of DNA barcode regions and phylogenetic construction methods. However, relatively few studies consider how estimates of phylogenetic diversity may be influenced by single DNA barcodes incorporated into a sequence matrix (conservative regions vs. hypervariable regions) and the use of a backbone family‐level phylogeny. Here, we use general linear mixed‐effects models to examine the influence of different combinations of core DNA barcodes (rbcL, matK, ITS, and ITS2) and phylogeny construction methods on a series of estimates of community phylogenetic diversity for two subtropical forest plots in Guangdong, southern China. We ask: (a) What are the relative influences of single DNA barcodes on estimates phylogenetic diversity metrics? and (b) What is the effect of using a backbone family‐level phylogeny to estimate topology‐based phylogenetic diversity metrics? The combination of more than one barcode (i.e., rbcL + matK + ITS) and the use of a backbone family‐level phylogeny provided the most parsimonious explanation of variation in estimates of phylogenetic diversity. The use of a backbone family‐level phylogeny showed a stronger effect on phylogenetic diversity metrics that are based on tree topology compared to those that are based on branch lengths. In addition, the variation in the estimates of phylogenetic diversity that was explained by the top‐rank models ranged from 0.1% to 31% and was dependent on the type of phylogenetic community structure metric. Our study underscores the importance of incorporating a multilocus DNA barcode and the use of a backbone family‐level phylogeny to infer phylogenetic diversity, where the type of DNA barcode employed and the phylogenetic construction method used can serve as a significant source of variation in estimates of phylogenetic community structure.     

Molecular phylogeny of the Oriental butterfly genus Arhopala (Lycaenidae,Theclinae) inferred from mitochondrial and nuclear genes

Abstract. We present a phylogeny for a selection of species of the butterfly genus Arhopala Boisduval, 1832 based on molecular characters. We sequenced 1778 bases of the mitochondrial genes Cytochrome Oxidase 1 and 2 including tRNA^Leu, and a 393‐bp fragment of the nuclear wingless gene for a total of 42 specimens of 33 species, representing all major species groups. Analyses of mtDNA and wingless genes show congruent phylogenetic signal. The phylogeny presented here confirms the monophyly of the centaurus, eumolphus, camdeo and epimuta groups and the amphimuta subgroup. It confirms close relationships between species within the agelastus group, that together with the amphimuta subgroup, centaurus and camdeo groups form a monophyletic group. However, incongruencies with previous taxonomic studies also occur; the amphimuta and silhetensis groups are not monophyletic, as is the genus Arhopala itself. One enigmatic species, A. kinabala, was evaluated further for topology and the support for basal placement of this species is due mainly to the wingless gene. However, in the Parsimony analysis, and subsequent Maximum Likelihood evaluations, certain nodes could not be resolved due to insufficient support. The mtDNA shows extreme AT bias with compositional heterogeneity at 3rd codon positions, which may result in saturation. By contrast, the wingless gene does not show compositional bias, suggesting that poor support is not due solely to saturation. The evaluation of morphological characters used in previous studies on Arhopala systematics on the molecular tree indicates that the macular pattern and the absence of tails at the hind wings show extensive homoplasy. A significant phylogenetic signal (as indicated by T‐PTP tests) is present in several of these morphological characters, which are nevertheless of limited use in phylogenetic studies due to their labile nature.     

The systematic utility of floral and vegetative fragrance in two genera of nyctaginaceae

We examined relationships between fragrance and phylogeny using a number of approaches to coding fragrance data and comparing the hierarchical information in fragrance data with the phylogenetic signal in a DNA sequence data set. We first used distance analyses to determine which coding method(s) best distinguishes species while grouping conspecifics. Results suggest that interspecific differences in fragrance composition were maximized by coding as presence/absence of fragrance compounds and biosynthetic pathways rather than when quantitative information was also included. Useful systematic information came from both compounds and pathways and from fragrance emitted by both floral and vegetative tissues. The coding methods that emerged from the distance analyses as best distinguishing species were then adapted for use in phylogenetic analysis. Although hierarchical signal among fragrance data sets was congruent, this signal was highly incongruent with the phylogenetic signal in the DNA sequence data. Notably, topologies inferred from fragrance data sets were congruent with the DNA topology only in the most distal portions (e.g., sister group pairs or closely related species that had similar fragrance profiles were often recovered by analyses of fragrance). Examination of consistency and retention indices for individual fragrance compounds and pathways as optimized onto one of the most-parsimonious trees inferred from DNA data revealed that although most compounds were homoplastic, some compounds were perfectly congruent with the DNA phylogeny. In particular, compounds and pathways found in a few taxa were less homoplastic than those found in many taxa. Pathways that synthesize few volatiles also seem to have lower homoplasy than those that produce many. Although fragrance data as a whole may not be useful in phylogeny reconstruction, these data can provide additional support for clades reconstructed with other types of characters. Factors other than phylogeny, including pollinator interactions, also likely influence fragrance composition.     

Phylogenetic signal in bone histology of amniotes revisited

Legendre, L, Le Roy, N, Martinez‐Maza, C, Montes, L, Laurin, M & Cubo, J. (2012). Phylogenetic signal in bone histology of amniotes revisited. —Zoologica Scripta, 42, 44–53. There is currently a debate about the presence of a phylogenetic signal in bone histological data, but very few rigorous tests have fuelled the discussions on this topic. Here, we performed new analyses using a larger set of seven histological traits and including 25 taxa (nine extinct and 16 extant taxa), using three methods: the phylogenetic eigenvector regression, the tree length distribution and the regressions on distance matrices. Our results clearly show that the phylogenetic signal in our sample of bone histological characters is strong, even after correcting for multiple testing. Most characters exhibit a significant phylogenetic signal according to at least one of our three tests, with the phylogeny often explaining 20–60% of the variation in the histological characters. Thus, we conclude that the phylogenetic comparative method should be systematically used in interspecific analyses of bone histodiversity to avoid problems of non‐independence among observations.     

Phylogeny of prickly poppies,Argemone (Papaveraceae), and the evolution of morphological and alkaloid characters based on ITS nrDNA sequence variation

Evolutionary relationships withinArgemone (Papaveraceae) were inferred from nucleotide variation in nuclear ribosomal DNA. A complete representation of the genus was achieved by using herbarium material to a large extent (74%). Four distinct clades can be recognized based on the molecular results. The support for the different clades varies greatly due to a very uneven distribution of characters. Although some clades are largely unresolved, some unexpected relationships, for example the sister group relationship of the relict speciesA. subintegrifolia to the rest ofArgemone, were found. The evolution of morphological and alkaloid characters against the background of the obtained phylogeny is discussed. Whereas the morphological characters provide some support for the phylogeny found, the alkaloid characters turned out not to contain much phylogenetic information. Past difficulties in recognizing phylogenetic relationships withinArgemone result from largely continuous morphological variation and from parallelisms correlated with and probably caused by polyploidization.     

Adult and early‐stage characters of Brassolini contain conflicting phylogenetic signal (Lepidoptera,Nymphalidae)

This study examines the contribution of early‐stages and adult characters to the reconstruction of the phylogeny of Brassolini butterflies. Parsimony analyses used both equal weights and implied weights, and a series of analyses were performed. First, we analysed adult and early‐stages partitions independently and in combination for a subset of 27 species; in these cases the matrices were mostly complete. Whereas the adult partition alone produced a topology that was well resolved and congruent with previous studies, the early‐stages partition produced a poorly resolved tree under equal weights. Furthermore, implied weights produced a well‐resolved early‐stages topology that differed significantly from the adult topology. When both partitions were combined for 27 species, implied weights yielded a topology that resembled the adult tree except for the positions of Bia and Penetes, but statistical node support was generally lower. This suggests that stochastic noise increased when early‐stage characters were added to the adult partition, but the combined partitions topology was not statistically different from that based on adult characters alone. Second, given that preserved early stages are not as readily available as adults, we analysed a matrix including 45 species in which early‐stage data were missing for 18 species, and compared the topology to that produced by the adult partition alone. Results were similar to the analyses including fewer species; the combined partitions tree was similar to that from the adult partition except for the position of Bia and Penetes. We compare our findings to other genus‐level phylogenetic studies within Lepidoptera that have also used early‐stages and adult characters.     

Developmental characters in phylogenetic inference and their absolute timing information

Despite the recent surge of interest in studying the evolution of development, surprisingly little work has been done to investigate the phylogenetic signal in developmental characters. Yet, both the potential usefulness of developmental characters in phylogenetic reconstruction and the validity of inferences on the evolution of developmental characters depend on the presence of such a phylogenetic signal and on the ability of our coding scheme to capture it. In a recent study, we showed, using simulations, that a new method (called the continuous analysis) using standardized time or ontogenetic sequence data and squared-change parsimony outperformed event pairing and event cracking in analyzing developmental data on a reference phylogeny. Using the same simulated data, we demonstrate that all these coding methods (event pairing and standardized time or ontogenetic sequence data) can be used to produce phylogenetically informative data. Despite some dependence between characters (the position of an event in an ontogenetic sequence is not independent of the position of other events in the same sequence), parsimony analysis of such characters converges on the correct phylogeny as the amount of data increases. In this context, the new coding method (developed for the continuous analysis) outperforms event pairing; it recovers a lower proportion of incorrect clades. This study thus validates the use of ontogenetic data in phylogenetic inference and presents a simple coding scheme that can extract a reliable phylogenetic signal from these data.     

Biogeography and evolution of Dolichandra (Bignonieae,Bignoniaceae)

Little information on evolutionary relationships of Neotropical organisms or on the factors that have shaped the diversity currently encountered in this region is available. However, it is clear that biotic interactions and abiotic aspects have played important roles for species diversification in the region. This study focuses on Dolichandra (Bignonieae, Bignoniaceae), a clade of Neotropical lianas that is distributed broadly across different habitats and with diverse pollination and dispersal systems. We used sequences from two plastid DNA markers (ndhF and rpl32‐trnL) and one nuclear gene (PepC) to infer phylogenetic relationships in Dolichandra using parsimony and Bayesian approaches. We then used this phylogenetic framework as basis to study the biogeographic history, reconstruct the evolution of morphological characters and test the impact of morphology and environment on the diversification of the genus. More specifically, we: (1) time‐calibrate the phylogenetic tree of Dolichandra; (2) estimate the ancestral areas of the various lineages; (3) estimate the ancestral states of discrete and continuous morphological traits; (4) test for phylogenetic signal in environmental and phenotypic data; and (5) test whether morphological characters and/or niche evolution are correlated with cladogenesis. All Dolichandra spp. are monophyletic in the combined molecular phylogeny; relationships among species are generally well resolved, although poorly supported in some instances. The genus is inferred to have originated 36.43–26.23 Mya, possibly in eastern South America. Ancestral state reconstructions of continuous and discrete floral characters inferred a mixed morphology as the ancestral condition for the group. Phylogenetic signal differed between perianth and sexual whorls and gradual evolution was recovered for all traits except style length and anther length. Environmental variables showed no phylogenetic signal and a pattern of variation that was not correlated with branch length, suggesting that environmental transitions were concomitant with speciation. Dispersal is inferred to be the main driver of the differential distribution observed among species. In addition, climatic preferences and floral characters seem to have been important reproductive barriers in Dolichandra. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 403–420.     

Recovering phylogenetic signal from frog mating calls

Goicoechea, N., De La Riva, I. & Padial, J. M. (2010). Recovering phylogenetic signal from frog mating calls. —Zoologica Scripta, 39, 411–154. Few studies have tried to analyse the phylogenetic information contained in frog mating calls. While some of those studies suggest that sexual selection deletes any phylogenetic signal, others indicate that frog calls do retain phylogenetic informative characters. Discordant results can be the outcome of disparate rates of character evolution and evolutionary plasticity of call characters in different groups of frogs, but also the result of applying different coding methods. No study to date has compared the relative performance of different coding methods in detecting phylogenetic signal in calls, hampering thus potential consilience between previous results. In this study, we analyse the strength of phylogenetic signal in 10 mating call characters of 11 related species of frogs belonging to three genera of Andean and Amazonian frogs (Anura: Terrarana: Strabomantidae). We use six quantitative characters (number of notes per call, note length, call length, number of pulses per note, fundamental frequency and dominant frequency) and four qualitative ones (presence/absence of: pseudopulses, frequency modulation in notes, amplitude modulation in notes and amplitude modulation in pulses). We code quantitative characters using four different coding and scaling methods: (i) gap‐coding, (ii) fixed‐scale, (iii) step‐matrix gap‐weighting with between‐characters scaling, and (iv) step‐matrix gap‐weighting with between‐states scaling. All four coding methods indicate that frog calls contain phylogenetic information. These results suggest that divergent selection on frog mating calls may not always be strong enough to eliminate phylogenetic signal. However, coding methods strongly affect the amount of recoverable information. Step‐matrix gap‐weighting with between‐characters scaling and gap‐coding are suggested as the best methods available for coding quantitative characters of frog calls. Also, our results indicate that the arbitrariness in selecting character states and the method for scaling transitions weights, rather than the number of character states, is what potentially biases phylogenetic analyses with quantitative characters.     

Phylogeny and reclassification of the Caucasigenini radiation from the Caucasus region (Gastropoda,Hygromiidae)

The Caucasigenini is an endemic radiation of hygromiid land snails from the Caucasus region. A phylogenetic analysis of morphological characters of the genitalia and the shell showed that the morphological characters are insufficient for resolving the relationships within the Caucasigenini. Convergences of the few parsimony informative characters in other groups of the Hygromiidae demonstrate that these characters are not reliable indicators of phylogenetic relationships. Phylogenetic analyses of sequences of cox1, 16S rDNA, 5.8S rDNA, ITS2 and 28S rDNA revealed several well‐supported groups. The relationships among these groups could not be resolved. It is likely that these groups originated in a rapid radiation during the uplift of the Caucasus. Based on the molecular phylogeny, we propose a new classification of the species of the Caucasigenini and establish a new genus, Lazicana gen. n.     

A phylogenetic evaluation of Leandra (Miconieae, Melastomataceae): a polyphyletic genus where the seeds tell the story, not the petals

Melastomataceae is a tropical family of 4500–5000 species divided into nine tribes. The largest tribe, Miconieae, is composed of approximately 2200 species in 30 genera and is found exclusively in the Neotropics. Previous phylogenetic analyses of the Miconieae have suggested that many of the genera are derived from a paraphyletic Miconia. However, these analyses only included six species of the large genus Leandra, so its phylogenetic affinities remained unclear. As currently defined Leandra is characterized by acute petals and terminal inflorescences, but some species of Miconia, Clidemia and Ossaea also have these characters. In this study, we present an analysis of nrITS sequence data for a sample of 63 species of Leandra. The genus is clearly resolved as polyphyletic, but some distinct and well‐supported clades exist. Some of these partially correspond to sections recognized in the nineteenth century by Cogniaux, or to geographic distribution. The distribution of seed structure characters is better correlated with the phylogeny than traditional characters, such as petal morphology. Seed appendages in Leandra have evolved independently at least four times. © The Willi Hennig Society 2007.     

Relationships within the family Actiniidae (Cnidaria,Actiniaria) based on molecular characters

Current sea anemone systematics is based on relatively few morphological characters, and potentially could benefit from the use of molecular characters. In this paper, the phylogenetic relationships of 12 species from 6 genera in the family Actiniidae have been investigated using electrophoretically separated isozymes. A numerical cladistic analysis has produced an estimated phylogeny. The implications of this phylogeny for the taxonomic use of certain morphological characters are discussed.     

Can phylogenetic signal,character displacement,or random phenotypic drift explain the morphological variation in the genus Geonoma (Arecaceae)?

Plant clades may exhibit little or wide morphological variation as a result of (1) the retention of ancestral characteristics or phylogenetic signal, (2) character displacement, or (3) random phenotypic drift or convergence. Understanding the taxonomy and systematics of many plant lineages has been challenging due to continuous intra‐ and interspecific morphological variation. To assess which evolutionary hypothesis could explain the morphological diversity in the genus Geonoma (Arecaceae), we performed a Mantel test between phylogenetic and morphological distances of 54 taxa, and tested for phylogenetic signal using Blomberg's K‐statistic on continuous variables, and a randomization of character states. To obtain a phylogenetic (patristic) distance matrix for Geonoma, we constructed a molecular phylogeny of tribe Geonomateae using three nuclear DNA regions. A positive relationship between the patristic and a 26‐discrete‐character distance matrix (R² = 0.55, P < 0.001) supported the phylogenetic signal hypothesis. The randomization test showed that signal was present in 16 characters. No relationship was evident using a 17‐quantitative‐variable distance matrix (R² = 0.07, P = 0.13), supporting the random drift hypothesis or convergence, and all 17 K‐values were close to 0, suggesting less phylogenetic signal than under the Brownian model. If most morphological variables traditionally used to classify Geonoma evolved randomly or convergently, it might explain Geonoma's challenging taxonomy. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106 , 528–539.     

Tooth development and histology patterns in lamniform sharks (Elasmobranchii,Lamniformes) revisited

The dentition of lamniforme sharks exhibits several characters that have been used extensively to resolve the phylogenetic relationships of extant taxa, yet some uncertainties remain. Also, the development of different teeth of a tooth file within the jaws of most extant lamniforms has not been documented to date. High‐resolution micro‐computed tomography is used here to re‐evaluate the importance of two dental characters within the order Lamniformes, which were considered not to be phylogenetically informative, the histotype and the number of teeth per tooth file. Additionally, the development and mineralization patterns of the teeth of the two osteodont lamniforms Lamna nasus and Alopias superciliosus were compared. We discuss the importance of these dental characters for phylogenetic interpretations to assess the quality of these characters in resolving lamniform relationships. The dental characters suggest that (1) Lamniformes are the only modern‐level sharks exhibiting the osteodont histotype, (2) the osteodont histotype in lamniform sharks is a derived state in modern‐level sharks (Elasmobranchii), (3) the osteodont type, conversely is convergently achieved when the clade Chondrichthyes is considered and thus might comprise a functional rather than a phylogenetic signal, and (4) there is an increase in the number of teeth per file throughout lamniform phylogeny. Structural development of the teeth of L. nasus and A. superciliosus is congruent with a previous investigation of the lamniform shark Carcharodon carcharias. J. Morphol. 277:1584–1598, 2016. © 2016 Wiley Periodicals, Inc.     

Reconstructing the evolution of complex bird song in the oropendolas

The elaborate songs of songbirds are frequent models for investigating the evolution of animal signals. However, few previous studies have attempted to reconstruct historical changes in song evolution using a phylogenetic comparative approach. In particular, no comparative studies of bird song have used a large number of vocal characters and a well-supported, independently derived phylogeny. We identified 32 features in the complex vocal displays of male oropendolas (genera Psarocolius, Gymnostinops, and Ocyalus) that are relatively invariant within taxa and mapped these characters onto a robust molecular phylogeny of the group. Our analysis revealed that many aspects of oropendola song are surprisingly evolutionarily conservative and thus are potentially useful characters for reconstructing historical patterns. Of the characters that varied among taxa, nearly two thirds (19 of 29) showed no evidence of evolutionary convergence or reversal when mapped onto the tree, which was reflected in a high overall consistency index (CI = 0.78) and retention index (RI = 0.88). Some reconstructed patterns provided evidence of selection on these signals. For example, rapid divergence of the songs of the Montezuma oropendola, Gymnostinops montezuma, from those of closely related taxa suggests the recent influence of strong sexual selection. In general, our results provide insights into the mode of vocal evolution in songbirds and suggest that complex vocalizations can provide information about phylogeny. Based on this evidence, we use song characters to estimate the phylogenetic affinities of three oropendola taxa for which molecular data are not yet available.     

Phylogenetic relationships of Ruteae (Rutaceae): new evidence from the chloroplast genome and comparisons with non-molecular data

Phylogenetic analyses of three cpDNA markers (matK, rpl16, and trnL–trnF) were performed to evaluate previous treatments of Ruteae based on morphology and phytochemistry that contradicted each other, especially regarding the taxonomic status of Haplophyllum and Dictamnus. Trees derived from morphological, phytochemical, and molecular datasets of Ruteae were then compared to look for possible patterns of agreement among them. Furthermore, non-molecular characters were mapped on the molecular phylogeny to identify uniquely derived states and patterns of homoplasy in the morphological and phytochemical datasets. The phylogenetic analyses determined that Haplophyllum and Ruta form reciprocally exclusive monophyletic groups and that Dictamnus is not closely related to the other genera of Ruteae. The different types of datasets were partly incongruent with each other. The discordant phylogenetic patterns between the phytochemical and molecular trees might be best explained in terms of convergence in secondary chemical compounds. Finally, only a few non-molecular synapomorphies provided support for the clades of the molecular tree, while most of the morphological characters traditionally used for taxonomic purposes were found to be homoplasious. Within the context of the phylogenetic relationships supported by molecular data, Ruta, the type genus for the family, can only be diagnosed by using a combination of plesiomorphic, homoplasious, and autapomorphic morphological character states.     

Phylogenetic signal in bone microstructure of sauropsids

In spite of the fact that the potential usefulness of bone histology in systematics has been discussed for over one and a half centuries, the presence of a phylogenetic signal in the variation of histological characters has rarely been assessed. A quantitative assessment of phylogenetic signal in bone histological characters could provide a justification for performing optimizations of these traits onto independently generated phylogenetic trees (as has been done in recent years). Here we present an investigation on the quantification of the phylogenetic signal in the following bone histological, microanatomical, and morphological traits in a sample of femora of 35 species of sauropsids: vascular density, vascular orientation, index of Haversian remodeling, cortical thickness, and cross-sectional area (bone size). For this purpose, we use two methods, regressions on distance matrices tested for significance using permutations (a Mantel test) and random tree length distribution. Within sauropsids, these bone microstructural traits have an optimal systematic value in archosaurs. In this taxon, a Mantel test shows that the phylogeny explains 81.8% of the variation of bone size and 86.2% of the variation of cortical thickness. In contrast, a Mantel test suggests that the phylogenetic signal in histological traits is weak: although the phylogeny explains 18.7% of the variation of vascular density in archosaurs, the phylogenetic signal is not significant either for vascular orientation or for the index of Haversian remodeling. However, Mantel tests seem to underestimate the proportion of variance of the dependent character explained by the phylogeny, as suggested by a PVR (phylogenetic eigenvector) analysis. We also deal with some complementary questions. First, we evaluate the functional dependence of bone vascular density on bone size by using phylogenetically independent contrasts. Second, we perform a variation partitioning analysis and show that the phylogenetic signal in bone vascular density is not a by-product of phylogentic signal in bone size. Finally, we analyze the evolution of cortical thickness in diapsids by using an optimization by squared change parsimony and discuss the functional significance of this character in terms of decreased buoyancy in crocodiles and mass saving in birds. These results are placed in the framework of the constructional morphology model, according to which the variation of a character in a clade has a historical (phylogenetic) component, a functional (adaptive) component, and a structural (architectural) component.     

A MOLECULAR PHYLOGENY OF THE DROSOPHILA WILLISTONI GROUP: CONFLICTS BETWEEN SPECIES CONCEPTS?

The six sibling species of the Neotropical Drosophila willistoni group have a long history in studies of evolutionary biology, yet to date only one molecular study, which used allozymes, has been published on the phylogeny of the group. Here we present a phylogeny of the siblings based on the sequences of two nuclear genes, period (per) and Alcohol dehydrogenase (Adh), as well as the mitochondrial gene Cytochrome oxidase I (COI). Taken individually, only per has a strong phylogenetic signal supporting a well-resolved phylogeny of the group, and this phylogeny is different from that obtained using allozymes. The COI dataset by itself produces trees that disagree with per, and neither that data nor the Adh data have a strong phylogenetic signal, as indicated by low bootstrap values for all analyses. Combining the Adh and COI datasets results in the same tree as per alone. Combining all three genes results in the same topology, which is strongly supported. Two problematic taxa, D. pavlovskiana and a “Carmody strain,” which were identified as potentially separate species based on reproductive isolation, clearly cluster in the phylogenetic analyses within D. paulistorum and D. equinoxialis, respectively. Thus, there appears to be a conflict between the biological species concept and the phylogenetic species concept.