Multilocus phylogenetics of a rapid radiation in the genus Thomomys (Rodentia: Geomyidae)

Species complexes undergoing rapid radiation present a challenge in molecular systematics because of the possibility that ancestral polymorphism is retained in component gene trees. Coalescent theory has demonstrated that gene trees often fail to match lineage trees when taxon divergence times are less than the ancestral effective population sizes. Suggestions to increase the number of loci and the number of individuals per taxon have been proposed; however, phylogenetic methods to adequately analyze these data in a coalescent framework are scarce. We compare two approaches to estimating lineage (species) trees using multiple individuals and multiple loci: the commonly used partitioned Bayesian analysis of concatenated sequences and a modification of a newly developed hierarchical Bayesian method (BEST) that simultaneously estimates gene trees and species trees from multilocus data. We test these approaches on a phylogeny of rapidly radiating species wherein divergence times are likely to be smaller than effective population sizes, and incomplete lineage sorting is known, in the rodent genus, Thomomys. We use seven independent noncoding nuclear sequence loci (total approximately 4300 bp) and between 1 and 12 individuals per taxon to construct a phylogenetic hypothesis for eight Thomomys species. The majority-rule consensus tree from the partitioned concatenated analysis included 14 strongly supported bipartitions, corroborating monophyletic species status of five of the eight named species. The BEST tree strongly supported only the split between the two subgenera and showed very low support for any other clade. Comparison of both lineage trees to individual gene trees revealed that the concatenation method appears to ignore conflicting signals among gene trees, whereas the BEST tree considers conflicting signals and downweights support for those nodes. Bayes factor analysis of posterior tree distributions from both analyses strongly favor the model underlying the BEST analysis. This comparison underscores the risks of overreliance on results from concatenation, and ignoring the properties of coalescence, especially in cases of recent, rapid radiations.     

Basal actinopterygian relationships: a mitogenomic perspective on the phylogeny of the "ancient fish"

The basal actinopterygians comprise four major lineages (polypteriforms, acipenseriforms, lepisosteids, and Amia) and have been collectively called "ancient fish." We investigated the phylogeny of this group of fishes in relation to teleosts using mitochondrial genomic (mitogenomic) data, and compared this to the various alternative phylogenetic hypotheses that have been proposed previously. In addition to the previously determined complete mitochondrial DNA (mtDNA) sequences from 14 teleosts and two outgroups, we used newly determined mitogenomic sequences of 12 purposefully chosen species representing all the ancient fish lineages plus related teleosts. This data set comprised concatenated nucleotide sequences from 12 protein-coding genes (excluding the ND6 gene and third codon positions) and 22 transfer RNA (tRNA) genes (stem regions only) and these data were subjected to maximum parsimony, maximum likelihood, and Bayesian analyses. The resultant trees from the three methods were well resolved and largely congruent, with most internal branches being supported by high statistical values. Mitogenomic data strongly supported not only the monophyly of the teleosts (osteoglossomorphs and above), but also a sister-group relationship between the teleosts and a clade comprising the acipenseriforms, lepisosteids, and Amia, with the polypteriforms occupying the most basal position in the actinopterygian phylogeny. Although the tree topology differed from any of the previously proposed hypotheses based on morphology, it exhibited congruence with a recently proposed novel hypothesis based on nuclear markers.     

A molecular phylogeny of the Canidae based on six nuclear loci

We have reconstructed the phylogenetic relationships of 23 species in the dog family, Canidae, using DNA sequence data from six nuclear loci. Individual gene trees were generated with maximum parsimony (MP) and maximum likelihood (ML) analysis. In general, these individual gene trees were not well resolved, but several identical groupings were supported by more than one locus. Phylogenetic analysis with a data set combining the six nuclear loci using MP, ML, and Bayesian approaches produced a more resolved tree that agreed with previously published mitochondrial trees in finding three well-defined clades, including the red fox-like canids, the South American foxes, and the wolf-like canids. In addition, the nuclear data set provides novel indel support for several previously inferred clades. Differences between trees derived from the nuclear data and those from the mitochondrial data include the grouping of the bush dog and maned wolf into a clade with the South American foxes, the grouping of the side-striped jackal (Canis adustus) and black-backed jackal (Canis mesomelas) and the grouping of the bat-eared fox (Otocyon megalotis) with the raccoon dog (Nycteruetes procyonoides). We also analyzed the combined nuclear + mitochondrial tree. Many nodes that were strongly supported in the nuclear tree or the mitochondrial tree remained strongly supported in the nuclear + mitochondrial tree. Relationships within the clades containing the red fox-like canids and South American canids are well resolved, whereas the relationships among the wolf-like canids remain largely undetermined. The lack of resolution within the wolf-like canids may be due to their recent divergence and insufficient time for the accumulation of phylogenetically informative signal.     

Reconciling gene trees of eastern North American minnows

Most eastern North American cyprinid fishes belong to a clade known as the "open posterior myodome" (OPM) minnows, but phylogenetic relationships within this clade have been difficult to ascertain. Previous attempts to resolve relationships among the generally benthic "chubs" and the more pelagic "shiners", that constitute the majority of OPM minnows, have led to highly discordant phylogenetic hypotheses. To further examine relationships among the OPM minnows, we utilized both a concatenated Bayesian approach and a coalescent-based species tree method to analyze data from six protein coding nuclear loci (Enc1, Ptr, Ryr3, Sh3px3, Tbr1, and Zic1), as well as the mitochondrial locus (Cytb). We focused our analyses on the chub-like genus Phenacobius, a group that has drifted topologically between other benthic chubs and the more pelagic shiners, and also included exemplar taxa from 11 other OPM lineages. Individual gene trees were highly discordant regarding relationships within Phenacobius and across the OPM clade. The concatenated Bayesian analysis and coalescent-based species tree reconstruction recovered slightly different phylogenetic topologies. Additionally, the posterior support values for clades using the coalescent-based approach were consistently lower than the concatenated analysis. However, Phenacobius was resolved as monophyletic and as the sister lineage to Erimystax regardless of the combined data approach taken. Furthermore, Phenacobius+Erimystax was recovered as more closely related to the shiners we examined than to other chubs. Relationships within Phenacobius varied depending on the combined phylogenetic method utilized. Our results highlight the importance of multi-locus, coalescent-based approaches for resolving the phylogeny of diverse clades like the eastern North American OPM minnows.     

Species tree estimation and the historical biogeography of heroine cichlids

Heroine cichlids are major components of the fish faunas in both Central America and the Caribbean. To examine the evolutionary patterns of how cichlids colonized both of these regions, we reconstructed the phylogenetic relationships among 23 cichlid lineages. We used three phylogenetically novel nuclear markers (Dystropin b, Myomesin1, and Wnt7b) in combination with sequence data from seven other gene regions (Nd2, Rag1, Enc1, Sreb2, Ptr, Plagl2, and Zic1) to elucidate the species tree of these cichlids. The species examined represent major heroine lineages in South America, Central America, and the Greater Antilles. The individual gene trees of these groups were topologically quite discordant. Therefore, we combined the genetic partitions and inferred the species tree using both concatenation and a coalescent-based Bayesian method. The two resulting phylogenetic topologies were largely concordant but differed in two fundamental ways. First, more nodes in the concatenated tree were supported with substantial or 100% Bayesian posterior support than in the coalescent-based tree. Second, there was a minor, but biogeographically critical, topological difference between the concatenated and coalescent-based trees. Nevertheless, both analyses recovered topologies consistent with the Greater Antillean heroines being phylogenetically nested within the largely Central American heroine radiation. This study suggests that reconstructions of cichlid phylogeny and historical biogeography should account for the vagaries of individual gene histories.     

Toward the resolution of an explosive radiation--a multilocus phylogeny of oceanic dolphins (Delphinidae)

Oceanic dolphins (Delphinidae) are the product of a rapid radiation that yielded ～36 extant species of small to medium-sized cetaceans that first emerged in the Late Miocene. Although they are a charismatic group of organisms that have become poster children for marine conservation, many phylogenetic relationships within Delphinidae remain elusive due to the slow molecular evolution of the group and the difficulty of resolving short branches from successive cladogenic events. Here I combine existing and newly generated sequences from four mitochondrial (mt) genes and 20 nuclear (nu) genes to reconstruct a well-supported phylogenetic hypothesis for Delphinidae. This study compares maximum-likelihood and Bayesian inference methods of several data sets including mtDNA, combined nuDNA, gene trees of individual nuDNA loci, and concatenated mtDNA+nuDNA. In addition, I contrast these standard phylogenetic analyses with the species tree reconstruction method of Bayesian concordance analysis (BCA). Despite finding discordance between mtDNA and individual nuDNA loci, the concatenated matrix recovers a completely resolved and robustly supported phylogeny that is also broadly congruent with BCA trees. This study strongly supports groupings such as Delphininae, Lissodelphininae, Globicephalinae, Sotalia+Delphininae, Steno+Orcaella+Globicephalinae, and Leucopleurus acutus, Lagenorhynchus albirostris, and Orcinus orca as basal delphinid taxa.     

Molecular phylogenetic relationships and phenotypic diversity in miniaturized toadlets, genus Brachycephalus (Amphibia: Anura: Brachycephalidae)

Toadlets of the genus Brachycephalus are endemic to the Atlantic rainforests of southeastern and southern Brazil. The 14 species currently described have snout-vent lengths less than 18 mm and are thought to have evolved through miniaturization: an evolutionary process leading to an extremely small adult body size. Here, we present the first comprehensive phylogenetic analysis for Brachycephalus, using a multilocus approach based on two nuclear (Rag-1 and Tyr) and three mitochondrial (Cyt b, 12S, and 16S rRNA) gene regions. Phylogenetic relationships were inferred using a partitioned Bayesian analysis of concatenated sequences and the hierarchical Bayesian method (BEST) that estimates species trees based on the multispecies coalescent model. Individual gene trees showed conflict and also varied in resolution. With the exception of the mitochondrial gene tree, no gene tree was completely resolved. The concatenated gene tree was completely resolved and is identical in topology and degree of statistical support to the individual mtDNA gene tree. On the other hand, the BEST species tree showed reduced significant node support relative to the concatenate tree and recovered a basal trichotomy, although some bipartitions were significantly supported at the tips of the species tree. Comparison of the log likelihoods for the concatenated and BEST trees suggests that the method implemented in BEST explains the multilocus data for Brachycephalus better than the Bayesian analysis of concatenated data. Landmark-based geometric morphometrics revealed marked variation in cranial shape between the species of Brachycephalus. In addition, a statistically significant association was demonstrated between variation in cranial shape and genetic distances estimated from the mtDNA and nuclear loci. Notably, B. ephippium and B. garbeana that are predicted to be sister-species in the individual and concatenated gene trees and the BEST species tree share an evolutionary novelty, the hyperossified dorsal plate.     

Basal euteleostean relationships: a mitogenomic perspective on the phylogenetic reality of the "Protacanthopterygii"

Higher-level relationships of the basal Euteleostei (=Protacanthopterygii) are so complex and controversial that at least nine different morphology-based phylogenetic hypotheses have been proposed during the last 30 years. Relationships of the Protacanthopterygii were investigated using mitochondrial genomic (mitogenomic) data from 34 purposefully chosen species (data for 12 species being newly determined during the study) that fully represented major basal euteleostean lineages and some basal teleosts plus neoteleosts as outgroups. Unweighted and weighted maximum parsimony (MP) and maximum likelihood (ML) analyses were conducted with the data set that comprised concatenated nucleotide sequences from 12 protein-coding genes (excluding the ND6 gene and 3rd codon positions) and 22 transfer RNA (tRNA) genes (stem regions only) from the 34 species. The resultant trees were well resolved and largely congruent, with most internal branches being supported by high statistical values. Monophyly of the protacanthopterygians was confidently rejected by the mitogenomic data. Of the five major monophyletic groups that received high statistical support within the protacanthopterygians, a clade comprising members of the alepocephaloids was unexpectedly nested within the Otocephala, sister-group of the euteleosts. The remaining four major monophyletic groups, on the other hand, occupied phylogenetic positions intermediate between the otocephalans and neoteleosts, with a clade comprising esociforms + salmoniforms being more basal to the argentinoids and osmeroids. Although interrelationships of the latter two clades (argentinoids and osmeroids) with the neoteleosts remained ambiguous, the present results indicated explicitly that the protacanthopterygians as currently defined merely represent a collective, polyphyletic group of the basal euteleosts, located between the basal teleosts (elopomorphs and below) and neoteleosts (stomiiforms and above).     

Phylogenetics of Anthyllis (Leguminosae: Papilionoideae: Loteae): Partial incongruence between nuclear and plastid markers, a long branch problem and implications for morphological evolution

Phylogenetic relationships in the genus Anthyllis (Leguminosae: Papilionoideae: Loteae) were investigated using data from the nuclear ribosomal internal transcribed spacer regions (ITS) and three plastid regions (psbA-trnH intergenic spacer, petB-petD region and rps16 intron). Bayesian and maximum parsimony (MP) analysis of a concatenated plastid dataset recovered well-resolved trees that are topologically similar, with many clades supported by unique indels. MP and Bayesian analyses of the ITS sequence data recovered trees that have several well-supported topological differences, both among analyses, and to trees inferred from the plastid data. The most substantial of these concerns A. vulneraria and A. lemanniana, whose placement in the parsimony analysis of the ITS data appears to be due to a strong long-branch effect. Analysis of the secondary structure of the ITS1 spacer showed a strong bias towards transitions in A. vulneraria and A. lemanniana, many of which were also characteristic of certain outgroup taxa. This may contribute to the conflicting placement of this clade in the MP tree for the ITS data. Additional conflicts between the plastid and ITS trees were more taxonomically focused. These differences may reflect the occurrence of reticulate evolution between closely related species, including a possible hybrid origin for A. hystrix. The patterns of incongruence between the plastid and the ITS data seem to correlate with taxon ranks. All of our phylogenetic analyses supported the monophyly of Anthyllis (incl. Hymenocarpos). Although they are often taxonomically associated with Anthyllis, the genera Dorycnopsis and Tripodion are shown here to be more closely related to other genera of Loteae. We infer up to six major clades in Anthyllis that are morphologically well-characterized, and which could be recognized as sections. Four of these agree with various morphology-based classifications, while the other two are novel. We reconstruct the evolution of several morphological characteristics found only in Anthyllis or tribe Loteae. Some of these characters support major clades, while others show evidence of homoplasy within Anthyllis.     

The phylogenetic systematics of blue‐tailed skinks (Plestiodon) and the family Scincidae

Blue‐tailed skinks (genus Plestiodon) are a common component of the terrestrial herpetofauna throughout their range in eastern Eurasia and North and Middle America. Plestiodon species are also frequent subjects of ecological and evolutionary research, yet a comprehensive, well‐supported phylogenetic framework does not yet exist for this genus. We construct a comprehensive molecular phylogeny of Plestiodon using Bayesian phylogenetic analyses of a nine‐locus data set comprising 8308 base pairs of DNA, sampled from 38 of the 43 species in the genus. We evaluate potential gene tree/species tree discordance by conducting phylogenetic analyses of the concatenated and individual locus data sets, as well as employing coalescent‐based methods. Specifically, we address the placement of Plestiodon within the evolutionary tree of Scincidae, as well as the phylogenetic relationships between Plestiodon species, and their taxonomy. Given our sampling of major Scincidae lineages, we also re‐evaluate ‘deep’ relationships within the family, with the goal of resolving relationships that have been ambiguous in recent molecular phylogenetic analyses. We infer strong support for several scincid relationships, including a major clade of ‘scincines’ and the inter‐relationships of major Mediterranean and southern African genera. Although we could not estimate the precise phylogenetic affinities of Plestiodon with statistically significant support, we nonetheless infer significant support for its inclusion in a large ‘scincine’ clade exclusive of Acontinae, Lygosominae, Brachymeles, and Ophiomorus. Plestiodon comprises three major geographically cohesive clades. One of these clades is composed of mostly large‐bodied species inhabiting northern Indochina, south‐eastern China (including Taiwan), and the southern Ryukyu Islands of Japan. The second clade comprises species inhabiting central China (including Taiwan) and the entire Japanese archipelago. The third clade exclusively inhabits North and Middle America and the island of Bermuda. A vast majority of interspecific relationships are strongly supported in the concatenated data analysis, but there is nonetheless significant conflict amongst the individual gene trees. Coalescent‐based gene tree/species tree analyses indicate that incongruence amongst the nuclear loci may severely obscure the phylogenetic inter‐relationships of the primarily small‐bodied Plestiodon species that inhabit the central Mexican highlands. These same analyses do support the sister relationship between Plestiodon marginatus Hallowell, 1861 and Plestiodon stimpsonii (Thompson, 1912), and differ with the mitochondrial DNA analysis that supports Plestiodon elegans (Boulenger, 1887) + P. stimpsonii. Finally, because the existing Plestiodon taxonomy is a poor representation of evolutionary relationships, we replace the existing supraspecific taxonomy with one congruent with our phylogenetic results. © 2012 The Linnean Society of London, Zoological Journal of the Linnean Society, 2012, 165 , 163–189.     

Phylogenetic analysis of RAD‐seq data: examining the influence of gene genealogy conflict on analysis of concatenated data

One of the major issues in phylogenetic analysis is that gene genealogies from different gene regions may not reflect the true species tree or history of speciation. This has led to considerable debate about whether concatenation of loci is the best approach for phylogenetic analysis. The application of Next‐generation sequencing techniques such as RAD‐seq generates thousands of relatively short sequence reads from across the genomes of the sampled taxa. These data sets are typically concatenated for phylogenetic analysis leading to data sets that contain millions of base pairs per taxon. The influence of gene region conflict among so many loci in determining the phylogenetic relationships among taxa is unclear. We simulated RAD‐seq data by sampling 100 and 500 base pairs from alignments of over 6000 coding regions that each produce one of three highly supported alternative phylogenies of seven species of Drosophila. We conducted phylogenetic analyses on different sets of these regions to vary the sampling of loci with alternative gene trees to examine the effect on detecting the species tree. Irrespective of sequence length sampled per region and which subset of regions was used, phylogenetic analyses of the concatenated data always recovered the species tree. The results suggest that concatenated alignments of Next‐generation data that consist of many short sequences are robust to gene tree/species tree conflict when the goal is to determine the phylogenetic relationships among taxa.     

Improving phylogenetic inference of mushrooms with RPB1 and RPB2 nucleotide sequences (Inocybe; Agaricales)

Approximately 3000 bp across 84 taxa have been analyzed for variable regions of RPB1, RPB2, and nLSU-rDNA to infer phylogenetic relationships in the large ectomycorrhizal mushroom genus Inocybe (Agaricales; Basidiomycota). This study represents the first effort to combine variable regions of RPB1 and RPB2 with nLSU-rDNA for low-level phylogenetic studies in mushroom-forming fungi. Combination of the three loci increases non-parametric bootstrap support, Bayesian posterior probabilities, and resolution for numerous clades compared to separate gene analyses. These data suggest the evolution of at least five major lineages in Inocybe-the Inocybe clade, the Mallocybe clade, the Auritella clade, the Inosperma clade, and the Pseudosperma clade. Additionally, many clades nested within each major lineage are strongly supported. These results also suggest the family Crepiodataceae sensu stricto is sister to Inocybe. Recognition of Inocybe at the family level, the Inocybaceae, is recommended.     

Phylogenetic analyses of a combined data set suggest that the Attheya lineage is the closest living relative of the pennate diatoms (Bacillariophyceae)

A Bayesian analysis of a seven gene data set was conducted to reconstruct phylogenetic relationships among a sample of centric and pennate diatoms and to test alternative hypotheses about the closest living relative of Bacillariophyceae. A lineage, composed of two Attheya species, was inferred to share the most recent common ancestor with Bacillariophyceae--a relationship that was also corroborated by the combined parsimony analysis. All competing hypotheses about the closest living relative of Bacillariophyceae were rejected because 100% of the trees in the post-burn-in sample in the Bayesian analysis supported the Attheya-Bacillariophyceae clade. According to a partitioned Bremer support analysis, the majority of the genes in the combined data matrix supported the Attheya--Bacillariophyceae relationship. The global topology of the phylogenetic tree indicated that a monophyletic group consisting of Thalassiosirales and Toxarium undulatum formed the deepest branch followed by a node uniting a clade composed of Bacillariophyceae/Attheya species and a lineage made up of Eucampia zoodiacus, Chaetocerotales, Lithodesmiales, Triceratiales, Biddulphiales and Cymatosirales. Except for the phylogenetic positions of Lithodesmiales, Thalassiosira sp and Skeletonema costatum, the optimal tree obtained from the combined parsimony analysis showed the same branching order of taxa as those seen in the consensus tree inferred from three independent Markov chain Monte Carlo analyses. Noteworthy findings are that Toxarium undulatum shares a strongly supported node with Thalassiosirales and that the genus Attheya is not a member of the Chaetocerotales lineage.     

Species tree estimation for a deep phylogenetic divergence in the New World monkeys (Primates: Platyrrhini)

The estimation of a robust phylogeny is a necessary first step in understanding the biological diversification of the platyrrhines. Although the most recent phylogenies are generally robust, they differ from one another in the relationship between Aotus and other genera as well as in the relationship between Pitheciidae and other families. Here, we used coding and non-coding sequences to infer the species tree and embedded gene trees of the platyrrhine genera using the Bayesian Markov chain Monte Carlo method for the multispecies coalescent (?BEAST) for the first time and to compared the results with those of a Bayesian concatenated phylogenetic analysis. Our species tree, based on all available sequences, shows a closer phylogenetic relationship between Atelidae and Cebidae and a closer relationship between Aotus and the Cebidae clade. The posterior probabilities are lower for these conflictive tree nodes compared to those in the concatenated analysis; this finding could be explained by some gene trees showing no concordant topologies between Aotus and the other genera. Moreover, the topology of our species tree also differs from the findings of previous molecular and morphological studies regarding the position of Aotus. The existence of discrepancies between morphological data, gene trees and the species tree is widely reported and can be related to processes such as incomplete lineage sorting or selection. Although these processes are common in species trees with low divergence, they can also occur in species trees with deep and rapid divergence. The sources of the inconsistency of morphological and molecular traits with the species tree could be a main focus of further research on platyrrhines.     

Phylogeny of iguanian lizards inferred from 29 nuclear loci, and a comparison of concatenated and species-tree approaches for an ancient, rapid radiation

Iguanian lizards form a diverse clade whose members have been the focus of many comparative studies of ecology, behavior, and evolution. Despite the importance of phylogeny to such studies, interrelationships among many iguanian clades remain uncertain. Within the Old World clade Acrodonta, Agamidae is sometimes found to be paraphyletic with respect to Chamaeleonidae, and recent molecular studies have produced conflicting results for many major clades. Within the largely New World clade Pleurodonta, relationships among the 12 currently recognized major subclades (mostly ranked as families) have been largely unresolved or poorly supported in previous studies. To clarify iguanian evolutionary history, we first infer phylogenies using concatenated maximum-likelihood (ML) and Bayesian analyses of DNA sequence data from 29 nuclear protein-coding genes for 47 iguanian and 29 outgroup taxa. We then estimate a relaxed-clock Bayesian chronogram for iguanians using BEAST. All three methods produce identical topologies. Within Acrodonta, we find strong support for monophyly of Agamidae with respect to Chamaeleonidae, and for almost all relationships within agamids. Within Pleurodonta, we find strong Bayesian support for almost all relationships, and strong ML support for some interfamilial relationships and for monophyly of almost all families (excepting Polychrotidae). Our phylogenetic results suggest a non-traditional biogeographic scenario in which pleurodonts originated in the Northern Hemisphere and subsequently spread southward into South America. The pleurodont portion of the tree is characterized by several very short, deep branches, raising the possibility of deep coalescences that may confound concatenated analyses. We therefore also use 27 of these genes to implement a coalescent-based species-tree approach for pleurodonts. Although this analysis strongly supports monophyly of the pleurodont families, interfamilial relationships are generally different from those in the concatenated tree, and support is uniformly poor. However, a species-tree analysis using only the seven most variable loci yields higher support and more congruence with the concatenated tree. This suggests that low support in the 27-gene species-tree analysis may be an artifact of the many loci that are uninformative for very short branches. This may be a general problem for the application of species-tree methods to rapid radiations, even with phylogenomic data sets. Finally, we correct the non-monophyly of Polychrotidae by recognizing the pleurodont genus Anolis (sensu lato) as a separate family (Dactyloidae), and we correct the non-monophyly of the agamid genus Physignathus by resurrection of the genus Istiurus for the former Physignathus lesueurii.     

Phylogenetic relationships among Phytophthora species inferred from sequence analysis of mitochondrially encoded cytochrome oxidase I and II genes

The phylogenetic relationships of 51 isolates representing 27 species of Phytophthora were assessed by sequence alignment of 568 bp of the mitochondrially encoded cytochrome oxidase II gene. A total of 1299 bp of the cytochrome oxidase I gene also were examined for a subset of 13 species. The cox II gene trees constructed by a heuristic search, based on maximum parsimony for a bootstrap 50% majority-rule consensus tree, revealed 18 species grouping into seven clades and nine species unaffiliated with a specific clade. The phylogenetic relationships among species observed on cox II gene trees did not exhibit consistent similarities in groupings for morphology, pathogenicity, host range or temperature optima. The topology of cox I gene trees, constructed by a heuristic search based on maximum parsimony for a bootstrap 50% majority-rule consensus tree for 13 species of Phytophthora, revealed 10 species grouping into three clades and three species unaffiliated with a specific clade. The groupings in general agreed with what was observed in the cox II tree. Species relationships observed for the cox II gene tree were in agreement with those based on ITS regions, with several notable exceptions. Some of these differences were noted in species in which the same isolates were used for both ITS and cox II analysis, suggesting either a differential rate of evolutionary divergence for these two regions or incorrect assumptions about alignment of ITS sequences. Analysis of combined data sets of ITS and cox II sequences generated a tree that did not differ substantially from analysis of ITS data alone, however, the results of a partition homogeneity test suggest that combining data sets may not be valid.     

Assessing phylogenetic resolution among mitochondrial, nuclear, and morphological datasets in Nothonotus darters (Teleostei: Percidae)

External morphological characters are the basis of our understanding of diversity and species relationships in many darter clades. The past decade has seen the publication of many studies utilizing mtDNA sequence data to investigate darter phylogenetics, but only recently have nuclear genes been used to investigate darter relationships. Despite a long tradition of use in darter systematics few studies have examined the phylogenetic utility of external morphological characters in estimating relationships among species in darter clades. We present DNA sequence data from the mitochondrial cytochrome b (cytb) gene, the nuclear encoded S7 intron 1, and discretely coded external morphological characters for all 20 species in the darter clade Nothonotus. Bayesian phylogenetic analyses result in phylogenies that are in broad agreement with previous studies. The cytb gene tree is well resolved, while the nuclear S7 gene tree lacks phylogenetic resolution, node support, and is characterized by a lack of reciprocal monophyly for many of the Nothonotus species. The phylogenies resulting from analysis of the morphological dataset lack resolution, but nodes present are found in the cytb and S7 gene trees. The highest resolution and node support is found in the Bayesian combined data phylogeny. Based on our results we propose continued exploration of the phylogenetic utility of external morphological characters in other darter clades. Given the extensive lack of reciprocal monophyly of species observed in the S7 gene tree we predict that nuclear gene sequences may have limited utility in intraspecific phylogeographic studies of Nothonotus darters.     

Testing the monophyly of the New Zealand and Australian endemic family Conoesucidae Ross based on combined molecular and morphological data (Insecta: Trichoptera: Sericostomatoidea)

Conoesucidae (Trichoptera, Insecta) are restricted to SE Australia, Tasmania and New Zealand. The family includes 42 described species in 12 genera, and each genus is endemic to either New Zealand or Australia. Although monophyly has been previously assumed, no morphological characters have been proposed to represent synapomorphies for the group. We collected molecular data from two mitochondrial genes (16S and cytochrome oxidase I), one nuclear gene (elongation factor 1-α) (2237–2277 bp in total), and 12 morphological characters to produce the first phylogeny of the family. We combined the molecular and morphological characters and performed both a maximum parsimony analysis and a Bayesian analysis to test the monophyly of the family, and to hypothesize the phylogeny among its genera. The parsimony analysis revealed a single most parsimonious tree with Conoesucidae being a monophyletic taxon and sistergroup to the Calocidae. The Bayesian inference produced a distribution of trees, the consensus of which is supported with posterior probabilities of 100% for 15 out of 22 possible ingroup clades including the most basal branch of the family, indicating strong support for a monophyletic Conoesucidae. The most parsimonious tree and the tree from the Bayesian analysis were identical except that the ingroup genus Pycnocentria changed position by jumping to a neighbouring clade. Based on the assumption that the ancestral conoesucid species was present on both New Zealand and Australia, a biogeographical analysis using the dispersal-vicariance criteria demonstrated that one or two (depending on which of the two phylogenetic reconstructions were applied) sympatric speciation events took place on New Zealand prior to a single, late dispersal from New Zealand to Australia.     

Target enrichment improves phylogenetic resolution in the genus Zanthoxylum (Rutaceae) and indicates both incomplete lineage sorting and hybridization events

Background and Aims Zanthoxylum is the only pantropical genus within Rutaceae, with a few species native to temperate eastern Asia and North America. Efforts using Sanger sequencing failed to resolve the backbone phylogeny of Zanthoxylum. In this study, we employed target-enrichment high-throughput sequencing to improve resolution. Gene trees were examined for concordance and sectional classifications of Zanthoxylum were evaluated. Off-target reads were investigated to identify putative single-copy markers for bait refinement, and low-copy markers for evidence of putative hybridization events.MethodsA custom bait set targeting 354 genes, with a median of 321 bp, was designed for Zanthoxylum and applied to 44 Zanthoxylum species and one Tetradium species as the outgroup. Illumina reads were processed via the HybPhyloMaker pipeline. Phylogenetic inferences were conducted using coalescent and maximum likelihood methods based on concatenated datasets. Concordance was assessed using quartet sampling. Additional phylogenetic analyses were performed on putative single and low-copy genes extracted from off-target reads.Key ResultsFour major clades are supported within Zanthoxylum: the African clade, the Z. asiaticum clade, the Asian–Pacific–Australian clade and the American–eastern Asian clade. While overall support has improved, regions of conflict are similar to those previously observed. Gene tree discordances indicate a hybridization event in the ancestor of the Hawaiian lineage, and incomplete lineage sorting in the American backbone. Off-target putative single-copy genes largely confirm on-target results, and putative low-copy genes provide additional evidence for hybridization in the Hawaiian lineage. Only two of the five sections of Zanthoxylum are resolved as monophyletic.ConclusionsTarget enrichment is suitable for assessing phylogenetic relationships in Zanthoxylum. Our phylogenetic analyses reveal that current sectional classifications need revision. Quartet tree concordance indicates several instances of reticulate evolution. Off-target reads are proven useful to identify additional phylogenetically informative regions for bait refinement or gene tree based approaches.     

Phylogenetic analysis and molecular signatures defining a monophyletic clade of heterocystous cyanobacteria and identifying its closest relatives

Detailed phylogenetic and comparative genomic analyses are reported on 140 genome sequenced cyanobacteria with the main focus on the heterocyst-differentiating cyanobacteria. In a phylogenetic tree for cyanobacteria based upon concatenated sequences for 32 conserved proteins, the available cyanobacteria formed 8–9 strongly supported clades at the highest level, which may correspond to the higher taxonomic clades of this phylum. One of these clades contained all heterocystous cyanobacteria; within this clade, the members exhibiting either true (Nostocales) or false (Stigonematales) branching of filaments were intermixed indicating that the division of the heterocysts-forming cyanobacteria into these two groups is not supported by phylogenetic considerations. However, in both the protein tree as well as in the 16S rRNA gene tree, the akinete-forming heterocystous cyanobacteria formed a distinct clade. Within this clade, the members which differentiate into hormogonia or those which lack this ability were also separated into distinct groups. A novel molecular signature identified in this work that is uniquely shared by the akinete-forming heterocystous cyanobacteria provides further evidence that the members of this group are specifically related and they shared a common ancestor exclusive of the other cyanobacteria. Detailed comparative analyses on protein sequences from the genomes of heterocystous cyanobacteria reported here have also identified eight conserved signature indels (CSIs) in proteins involved in a broad range of functions, and three conserved signature proteins, that are either uniquely or mainly found in all heterocysts-forming cyanobacteria, but generally not found in other cyanobacteria. These molecular markers provide novel means for the identification of heterocystous cyanobacteria, and they provide evidence of their monophyletic origin. Additionally, this work has also identified seven CSIs in other proteins which in addition to the heterocystous cyanobacteria are uniquely shared by two smaller clades of cyanobacteria, which form the successive outgroups of the clade comprising of the heterocystous cyanobacteria in the protein trees. Based upon their close relationship to the heterocystous cyanobacteria, the members of these clades are indicated to be the closest relatives of the heterocysts-forming cyanobacteria.