Circularity,evolution, systematics … and circularity

Many have argued strongly that incorporation of evolutionary theory into systematics is dangerously circular, while others have maintained that such an integrated approach increases the accuracy of phylogenetic inference. Here, it is demonstrated that such blanket statements regarding exclusion or inclusion of evolutionary principles in systematics fail to distinguish between two very different types of principles. ‘Phylogeny-neutral’ evolutionary principles are those inferred without any recourse to specific phylogenetic hypotheses (e.g. via developmental genetics, biomechanics). In contrast, ‘phylogeny-dependent’ principles are those which can only be inferred on the basis of specific phylogenetic hypotheses (e.g. character associations detected via ‘comparative methods’). Inclusion of phylogeny-neutral principles in systematic studies as a priori assumptions can be justified, since these principles have (often strong) external empirical support from other spheres of investigation. However, inclusion of phylogeny-dependent principles in systematic studies is circular, since such principles have no external empirical support but are themselves derived from systematic studies. Advocating inclusion or exclusion of all (or as many as possible) evolutionary principles from phylogenetic analysis is therefore misguided. Rather, phylogeny-neutral principles are independently supported and can be included, while phylogeny-dependent principles are unjustified assumptions and should be excluded to avoid circularity. However, integration of complex phylogeny-neutral principles in systematics can create operational problems, even though there are no methodological reasons against their inclusion.     

Application of multigene phylogenetics and site‐stripping to resolve intraordinal relationships in the Rhodymeniales (Rhodophyta)

Previous molecular assessments of the red algal order Rhodymeniales have confirmed its monophyly and distinguished the six currently recognized families (viz. Champiaceae, Faucheaceae, Fryeellaceae, Hymenocladiaceae, Lomentariaceae, and Rhodymeniaceae); however, relationships among most of these families have remained unresolved possibly as a result of substitution saturation at deeper phylogenetic nodes. The objective of the current study was to improve rhodymenialean systematics by increasing taxonomic representation and using a more robust multigene dataset of mitochondrial (COB, COI/COI‐5P), nuclear (LSU, EF2) and plastid markers (psbA, rbcL). Additionally, we aimed to prevent phylogenetic inference problems associated with substitution saturation (particularly at the interfamilial nodes) by removing fast‐evolving sites and analyzing a series of progressively more conservative alignments. The Rhodymeniales was resolved as two major lineages: (i) the Fryeellaceae as sister to the Faucheaceae and Lomentariaceae; and (ii) the Rhodymeniaceae allied to the Champiaceae and Hymenocladiaceae. Support at the interfamilial nodes was highest when 20% of variable sites were removed. Inclusion of Binghamiopsis, Chamaebotrys, and Minium, which were absent in previous phylogenetic investigations, established their phylogenetic affinities while assessment of two genera consistently polyphyletic in phylogenetic analyses, Erythrymenia and Lomentaria, resulted in the proposition of the novel genera Perbella and Fushitsunagia. The taxonomic position of Drouetia was reinvestigated with re‐examination of holotype material of D. coalescens to clarify tetrasporangial development in this genus. In addition, we added three novel Australian species to Drouetia as a result of ongoing DNA barcoding assessments—D. aggregata sp. nov., D. scutellata sp. nov., and D. viridescens sp. nov.     

Improved phylogenetic resolution of toxic and non-toxic Alexandrium strains using a concatenated rDNA approach

Dinoflagellates of the genus Alexandrium are known producers of paralytic shellfish toxins. Species within the genus have similar phenotypes making morphological identification problematical. The use of Alexandrium rDNA sequence data is therefore increasing, resulting in the improved resolution of evolutionary relationships by phylogenetic inferences. However, the true branching pattern within Alexandrium remains unresolved, with minimal support shown for the main phylogentic branch. The aim of this study is to improve phylogenetic resolution via a concatenated rDNA approach with a broad sample of taxa, allowing inference of the evolutionary pattern between species and toxins. 27 Alexandrium strains from 10 species were tested with HPLC for PSP toxin presence and additionally sequenced for 18S, ITS1, 5.8S, ITS2 and 28S rDNA before being phylogenetically inferred together with all available orthologous sequences from NCBI. The resulting alignment is the largest to date for the genus, in terms of both inferred characters and taxa, thus allowing for the improved phylogenetic resolution of evolutionary patterns there in. No phylogenetic pattern between PSP producing and non-producing strains could be established, however the terminal tamarense complex was shown to produce more PSP analogues than basal clades. Additionally, we distinguish a high number of polymorphic regions between the two copies of A. fundyense rDNA, thus allowing us to demonstrate the presence of chimeric sequences within GenBank, as well as a possible over estimation of diversification within the tamarense complex.     

Molecular phylogeny of Lymantriinae (Lepidoptera,Noctuoidea, Erebidae) inferred from eight gene regions

To understand the evolutionary history of Lymantriinae and test the present higher‐level classification, we performed the first broad‐scale molecular phylogenetic analysis of the subfamily, based on 154 exemplars representing all recognized tribes and drawn from all major biogeographical regions. We used two mitochondrial genes (cytochrome c oxidase subunit I and 16S ribosomal RNA) and six nuclear genes (elongation factor‐1α, carbamoylphosphate synthase domain protein, ribosomal protein S5, cytosolic malate dehydrogenase, glyceraldehyde‐3‐phosphate dehydrogenase and wingless). Data matrices (in total 5424 bp) were analysed by parsimony and model‐based evolutionary methods (maximum likelihood and Bayesian inference). Based on the results of the analyses, we present a new phylogenetic classification for Lymantriinae composed of seven well‐supported tribes, two of which are proposed here as new: Arctornithini, Leucomini, Lymantriini, Orgyiini, Nygmiini, Daplasini trib. nov. and Locharnini trib. nov. We discuss the internal structure of each of these tribes and address some of the more complex problems with the genus‐level classification, particularly within Orgyiini and Nygmiini.     

Phylogenetic signals from point mutations and polymorphic Alu

Allelic frequency data derived from five polymorphic Alu insertion loci and five point mutation polymorphic loci were compared to determine their ability to infer phylogenetic relationships among human populations. While point mutation polymorphisms inferred a monophyletic Caucasian clade that is corroborated by other studies, these data failed to support the generally accepted monophyly of Orientals with native Americans. In addition, there is less statistical bootstrap support for the maximum-likelihood tree derived from the point mutation polymorphisms as compared to those generated from either the Alu insertion data or the combined Alu insertion+point mutation data. The Alu data and the combined Alu insertion+point mutation data inferred a monophyletic relationship among the Oriental and native American populations. The Alu insertion data and the combined Alu insertion+point mutation data also displayed two separate, well defined, tight clusters of the Caucasian and the Oriental+native American populations which was not inferred from the point mutation data. These findings indicate greater phylogenetic information contained in Alu insertion frequencies than in allelic frequencies derived from point-mutations. This revised version was published online in July 2006 with corrections to the Cover Date.     

Inclusion of rare taxa from Blattidae and Anaplectidae improves phylogenetic resolution in the cockroach superfamily Blattoidea

Cockroaches are an ecologically and economically important insect group, but some fundamental aspects of their evolutionary history remain unresolved. In particular, there are outstanding questions about some of the deeper relationships among cockroach families. As a group transferred from Blaberoidea Saussure to Blattoidea Latreille, the evolutionary history of the family Anaplectidae Walker requires re-evaluation. In our study, we infer the phylogeny of Blattoidea based on the mitochondrial genomes of 28 outgroup taxa and 67 ingroup taxa, including 25 newly sequenced blattoid species mainly from the families Anaplectidae and Blattidae Latreille. Our results indicate that Blattoidea is the sister group of the remaining Blattodea Brunner von Wattenwyl and that Blattoidea can be divided into three main clades: Blattidae + Tryonicidae McKittrick & Mackerras, Lamproblattidae McKittrick + Anaplectidae and Termitoidae Latreille + Cryptocercidae Handlirsch. Our analyses provide robust support for previously uncertain hypotheses. The sister group of Termitoidae + Cryptocercidae (Xylophagodea Engel) is inferred to constitute the rest of Blattoidea, for the first time. Within Blattidae, Hebardina Bey-Bienko is placed as the sister lineage to the rest of Blattidae. The subfamily Archiblattinae is polyphyletic, Blattinae is paraphyletic and Polyzosteriinae is monophyletic (Macrocercinae Roth not included); the genus Periplaneta Burmrister is polyphyletic. Based on the results of our phylogenetic analyses, we have revised these taxa. A new subfamily, Hebardininae subfam.nov. , is proposed in Blattidae. Archiblattinae and Shelfordella Adelung are synonymized with Blattinae and Periplaneta, respectively: Archiblattinae Kirby syn.nov. and Shelfordella Adelung syn.nov. Our inferred divergence times indicate that Blattoidea emerged in the Late Triassic, with six families in Blattoidea diverging in the Middle and Late Jurassic. We suggest that the divergences among lineages of Asian Blattidae and Anaplectidae were driven by the uplift of the Himalayas and deglaciation during the Quaternary, leading to the present-day distributions of these taxa.     

Ongoing hybridization obscures phylogenetic relationships in the Drosophila subquinaria species complex

Inferring evolutionary relationships among recently diverged lineages is necessary to understand how isolating barriers produce independent lineages. Here, we investigate the phylogenetic relationships between three incompletely isolated and closely related mushroom‐feeding Drosophila species. These species form the Drosophila subquinaria species complex and consist of one Eurasian species (D. transversa) and two widespread North American species (D. subquinaria and D. recens) that are sympatric in central Canada. Although patterns of pre‐ and post‐mating isolation among these species are well characterized, previous work on their phylogenetic relationships is limited and conflicting. In this study, we generated a multi‐locus data set of 29 loci from across the genome sequenced in a population sample from each species, and then, we inferred species relationships and patterns of introgression. We find strong statistical support that D. subquinaria is paraphyletic, showing that samples from the geographic region sympatric with D. recens are most closely related to D. recens, whereas samples from the geographic region allopatric with D. recens are most closely related to D. transversa. We present several lines of evidence that both incomplete lineage sorting and gene flow are causing phylogenetic discordance. We suggest that ongoing gene flow primarily from D. recens into D. subquinaria in the sympatric part of their ranges causes phylogenetic uncertainty in the evolutionary history of these species. Our results highlight how population genetic data can be used to disentangle the sources of phylogenetic discordance among closely related species.     

Is RAD‐seq suitable for phylogenetic inference? An in silico assessment and optimization

Inferring phylogenetic relationships between closely related taxa can be hindered by three factors: (1) the lack of informative molecular variation at short evolutionary timescale; (2) the lack of established markers in poorly studied taxa; and (3) the potential phylogenetic conflicts among different genomic regions due to incomplete lineage sorting or introgression. In this context, Restriction site Associated DNA sequencing (RAD‐seq) seems promising as this technique can generate sequence data from numerous DNA fragments scattered throughout the genome, from a large number of samples, and without preliminary knowledge on the taxa under study. However, divergence beyond the within‐species level will necessarily reduce the number of conserved and non‐duplicated restriction sites, and therefore the number of loci usable for phylogenetic inference. Here, we assess the suitability of RAD‐seq for phylogeny using a simulated experiment on the 12 Drosophila genomes, with divergence times ranging from 5 to 63 million years. These simulations show that RAD‐seq allows the recovery of the known Drosophila phylogeny with strong statistical support, even for relatively ancient nodes. Notably, this conclusion is robust to the potentially confounding effects of sequencing errors, heterozygosity, and low coverage. We further show that clustering RAD‐seq data using the BLASTN and SiLiX programs significantly improves the recovery of orthologous RAD loci compared with previously proposed approaches, especially for distantly related species. This study therefore validates the view that RAD sequencing is a powerful tool for phylogenetic inference.     

Catching the phylogenic history through the ontogenic hourglass: a phylogenomic analysis of Drosophila body segmentation genes

SUMMARY The phylogenetic information content of different developmental stages is a long‐standing issue in the study of development and evolution. We performed phylogenetic analyses of 51 body segmentation genes in 12 species of Drosophila in order to investigate the impact of the mode of evolution of development on phylogeny inference. Previous studies of these genes in Drosophila using pairwise phenetic comparisons at the species group level revealed the presence of an “hourglass model” (HG), wherein mid‐embryonic stages are the most evolutionarily constrained. We utilized two character‐based approaches: taxonomic congruence using the relative consensus fork index (RCFI), in which phylogenies are inferred from each gene separately and compared with a total evidence tree (TET), and partitioned simultaneous analysis using several indices such as branch support (BS) and localized incongruence length difference (LILD) test. We also proposed a new index, the recapitulatory index (R), which divides the number of synapomorphies on the total number of informative characters in a data set. Polynomial adjustment of both BS and R indices showed strong support for the hourglass model regardless of the taxonomic level (species subgroup vs. subgenera), showing less phylogenetic information content for mid‐developmental stages (mainly the zygotic segment polarity stage). Significant LILD scores were randomly distributed among developmental stages revealing the absence of differential selective constraints, but were significantly related to chromosomal location showing physical (linkage) impact on phylogenetic incongruence. RCFI was the most sensitive measure to taxonomic level, having a convex parabola at the species subgroup level in support of the hourglass model and a concave parabola at the subgeneric level in support of the adaptive penetrance model. This time‐dependent discrepancy of best fit developmental model parallels previous conflicting results from the vertebrates. Because of the quasi‐phenetic nature of this index, we argue that the discrepancy is due to the evolutionary rate heterogeneity of developmental genes rather than to fundamental differences among organisms. We suggest that simultaneous character‐based analyses give better macroevolutionary support to the hourglass model of the developmental constraints on genome evolution than pairwise phenetic comparisons.     

ClipKIT: A multiple sequence alignment trimming software for accurate phylogenomic inference

Highly divergent sites in multiple sequence alignments (MSAs), which can stem from erroneous inference of homology and saturation of substitutions, are thought to negatively impact phylogenetic inference. Thus, several different trimming strategies have been developed for identifying and removing these sites prior to phylogenetic inference. However, a recent study reported that doing so can worsen inference, underscoring the need for alternative alignment trimming strategies. Here, we introduce ClipKIT, an alignment trimming software that, rather than identifying and removing putatively phylogenetically uninformative sites, instead aims to identify and retain parsimony-informative sites, which are known to be phylogenetically informative. To test the efficacy of ClipKIT, we examined the accuracy and support of phylogenies inferred from 14 different alignment trimming strategies, including those implemented in ClipKIT, across nearly 140,000 alignments from a broad sampling of evolutionary histories. Phylogenies inferred from ClipKIT-trimmed alignments are accurate, robust, and time saving. Furthermore, ClipKIT consistently outperformed other trimming methods across diverse datasets, suggesting that strategies based on identifying and retaining parsimony-informative sites provide a robust framework for alignment trimming.

Highly divergent sites in multiple sequence alignments are thought to negatively impact phylogenetic inference; trimming methods aim to remove these sites, but recent analysis suggests that doing so can worsen inference. This study introduces ClipKIT, a trimming method that instead aims to retain parsimony-informative sites; phylogenetic inference using ClipKIT-trimmed alignments is accurate, robust and time-saving.     

A molecular phylogenetic analysis of Bulinus (Gastropoda: Planorbidae) with conserved nuclear genes

Jørgensen, A., Madsen, H., Nalugwa, A., Nyakaana, S., Rollinson, D., Stothard, J. R. & Kristensen, T. K. A molecular phylogenetic analysis of Bulinus (Gastropoda: Planorbidae) with conserved nuclear genes. —Zoologica Scripta, 40, 126–136. Mutational saturation of inspected DNA loci and topological incongruence in the phylogenetic inferences have previously confounded attempts to resolve the evolutionary relationships within the freshwater snail genus Bulinus. Traditionally, the 37 species of Bulinus are placed within the four species groups and the evolutionary divergence between groups is substantial. With an intention to shed new light on species group relationships, the present study was designed to investigate the basal divergences in the phylogeny of Bulinus using highly conserved nuclear genes. The resolved phylogeny inferred that the four species groups of Bulinus were monophyletic and Shimodaira‐Hasegawa topology tests found them to be significantly supported. The Bulinus truncatus/tropicus species complex and Bulinus wrighti (Bulinus reticulatus group) formed a well‐supported sister‐group relationship. The Bulinus africanus species group was the sister‐group to the clade (Bulinus truncatus/tropicus + B. wrighti) with the Bulinus forskalii species group as the sister‐group to these taxa. The sister‐group relationship between Indoplanorbis and Bulinus was non‐significant and the basal clade support of Bulinus improved upon exclusion of Indoplanorbis. The finding of basal long branches of Bulinus species originating from Madagascar strongly suggests the presence of additional cryptic species and an evolutionary scenario influenced by this island’s geological vicariance from the African mainland. Speciation by polyploidy was inferred to have evolved within a clade in the Bulinus truncatus/tropicus species complex. Although the monophyletic status of each species group was firmly supported, it was difficult to establish species group concepts equally across the variations and place this precisely in a specific temporal framework.     

Exploration of data partitioning in an eight‐gene data set: phylogeny of metalmark moths (Lepidoptera,Choreutidae)

Rota, J. & Wahlberg, N. (2012). Exploration of data partitioning in an eight‐gene data set: phylogeny of metalmark moths (Lepidoptera, Choreutidae). —Zoologica Scripta, 41, 536–546. Molecular data sets for phylogenetic inference continue to increase in size, especially with respect to the number of genes sampled. As more and more genes are included in analyses, the importance of partitioning the data to avoid problems that can arise from underparameterization becomes more apparent. With an eight‐gene data set from 38 metalmark moth species (12 genera represented) and three outgroups, we explored different data partitioning strategies and their influence on convergence and mixing of Markov Chains Monte Carlo in a Bayesian setting. We found that in larger data sets, with an increase in the number of partitions that are made a priori (e.g. by gene and codon position), convergence and mixing become poor. This problem can be overcome by using a recently published algorithm in which homologous sites are grouped into blocks with similar evolutionary rates that can then be modelled as separate data subsets. Using this novel approach to data partitioning, our analyses resolve with strong support relationships among the genera of metalmark moths. Support for the monophyly of the family, the two subfamilies and all genera except Hemerophila is strong. Hemerophila is broken into two separate clades, Hemerophila sensu stricto and another well‐supported clade. To render Hemerophila monophyletic, we describe a new genus, Ornarantia Rota, gen. nov., and transfer 18 species from Hemerophila to it. The type species of Ornarantia is Hemerophila laciniosella Busck, 1914.     

Legumes select symbiosis island sequence variants in Bradyrhizobium

Bradyrhizobium strains sampled from 14 legume genera native to eastern North America showed substantial host‐related phylogenetic clustering at three loci in the symbiotic island (SI) region (nodC, nifD, nifH), indicating selection of distinct suites of SI lineages by different legumes. Bacteria assorted consistently with particular legumes across two regions separated by 800 km, implying recurrent assembly of the same symbiotic combinations. High genetic polymorphism of all three SI loci relative to four nonsymbiotic loci supported the inference that a form of multiple‐niche balancing selection has acted on the SI region, arising from differential symbiont utilization by different legume taxa. Extensive discordance between the tree for SI variants and a phylogenetic tree inferred for four housekeeping loci implied that lateral transfer of the symbiosis island region has been common (at least 26 transfer events among 85 Bradyrhizobium strains analysed). Patterns of linkage disequilibrium also supported the conclusion that recombination has impacted symbiotic and nonsymbiotic regions unequally. The high prevalence of lateral transfer suggests that acquisition of a novel SI variant may often confer a strong selective advantage for recipient cells.     

Phylogenetic signal and the utility of 12S and 16S mtDNA in frog phylogeny

Genes selected for a phylogenetic study need to contain conserved information that reflects the phylogenetic history at the specific taxonomic level of interest. Mitochondrial ribosomal genes have been used for a wide range of phylogenetic questions in general and in anuran systematics in particular. We checked the plausibility of phylogenetic reconstructions in anurans that were built from commonly used 12S and 16S rRNA gene sequences. For up to 27 species arranged in taxon sets of graded inclusiveness, we inferred phylogenetic hypotheses based on different a priori decisions, i.e. choice of alignment method and alignment parameters, including/excluding variable sites, choice of reconstruction algorithm and models of evolution. Alignment methods and parameters, as well as taxon sampling all had notable effects on the results leading to a large number of conflicting topologies. Very few nodes were supported in all of the analyses. Data sets in which fast evolving and ambiguously aligned sites had been excluded performed worse than the complete data sets. There was moderate support for the monophyly of the Discoglossidae, Pelobatoidea, Pelobatidae and Pipidae. The clade Neobatrachia was robustly supported and the intrageneric relationships within Bombina and Discoglossus were well resolved indicating the usefulness of the genes for relatively recent phylogenetic events. Although 12S and 16S rRNA genes seem to carry some phylogenetic signal of deep (Mesozoic) splitting events the signal was not strong enough to resolve consistently the inter‐relationships of major clades within the Anura under varied methods and parameter settings.     

Evolutionary history of alpine and subalpine Daphnia in western North America

相似文献   

Cryptic diversity in the Mexican highlands: Thousands of UCE loci help illuminate phylogenetic relationships,species limits and divergence times of montane rattlesnakes (Viperidae: Crotalus)

With the continued adoption of genome‐scale data in evolutionary biology comes the challenge of adequately harnessing the information to make accurate phylogenetic inferences. Coalescent‐based methods of species tree inference have become common, and concatenation has been shown in simulation to perform well, particularly when levels of incomplete lineage sorting are low. However, simulation conditions are often overly simplistic, leaving empiricists with uncertainty regarding analytical tools. We use a large ultraconserved element data set (>3,000 loci) from rattlesnakes of the Crotalus triseriatus group to delimit lineages and estimate species trees using concatenation and several coalescent‐based methods. Unpartitioned and partitioned maximum likelihood and Bayesian analysis of the concatenated matrix yield a topology identical to coalescent analysis of a subset of the data in bpp . ASTRAL analysis on a subset of the more variable loci also results in a tree consistent with concatenation and bpp , whereas the SVDquartets phylogeny differs at additional nodes. The size of the concatenated matrix has a strong effect on species tree inference using SVDquartets , warranting additional investigation on optimal data characteristics for this method. Species delimitation analyses suggest up to 16 unique lineages may be present within the C. triseriatus group, with divergences occurring during the Neogene and Quaternary. Network analyses suggest hybridization within the group is relatively rare. Altogether, our results reaffirm the Mexican highlands as a biodiversity hotspot and suggest that coalescent‐based species tree inference on data subsets can provide a strongly supported species tree consistent with concatenation of all loci with a large amount of missing data.     

Evolution of the RNA polymerase B' subunit gene (rpoB') in Halobacteriales: a complementary molecular marker to the SSU rRNA gene

Many prokaryotes have multiple ribosomal RNA operons. Generally, sequence differences between small subunit (SSU) rRNA genes are minor (<1%) and cause little concern for phylogenetic inference or environmental diversity studies. For Halobacteriales, an order of extremely halophilic, aerobic Archaea, within-genome SSU rRNA sequence divergence can exceed 5%, rendering phylogenetic assignment problematic. The RNA polymerase B' subunit gene (rpoB') is a single-copy conserved gene that may be an appropriate alternative phylogenetic marker for Halobacteriales. We sequenced a fragment of the rpoB' gene from 21 species, encompassing 15 genera of Halobacteriales. To examine the utility of rpoB' as a phylogenetic marker in Halobacteriales, we investigated three properties of rpoB' trees: the variation in resolution between trees inferred from the rpoB' DNA and RpoB' protein alignment, the degree of mutational saturation between taxa, and congruence with the SSU rRNA tree. The rpoB' DNA and protein trees were for the most part congruent and consistently recovered two well-supported monophyletic groups, the clade I and clade II haloarchaea, within a collection of less well resolved Halobacteriales lineages. A comparison of observed versus inferred numbers of substitution revealed mutational saturation in the rpoB' DNA data set, particularly between more distant species. Thus, the RpoB' protein sequence may be more reliable than the rpoB' DNA sequence for inferring Halobacteriales phylogeny. AU tests of tree selection indicated the trees inferred from rpoB' DNA and protein alignments were significantly incongruent with the SSU rRNA tree. We discuss possible explanations for this incongruence, including tree reconstruction artifact, differential paralog sampling, and lateral gene transfer. This is the first study of Halobacteriales evolution based on a marker other than the SSU rRNA gene. In addition, we present a valuable phylogenetic framework encompassing a broad diversity of Halobacteriales, in which novel sequences can be inserted for evolutionary, ecological, or taxonomic investigations.