Molecular evidence for the monophyly of East Asian groups of Cyprinidae (Teleostei: Cypriniformes) derived from the nuclear recombination activating gene 2 sequences

The family Cyprinidae is one of the largest families of fishes in the world and a well-known component of the East Asian freshwater fish fauna. However, the phylogenetic relationships among cyprinids are still poorly understood despite much effort paid on the cyprinid molecular phylogenetics. Original nucleotide sequence data of the nuclear recombination activating gene 2 were collected from 109 cyprinid species and four non-cyprinid cypriniform outgroup taxa and used to infer the cyprinid phylogenetic relationships and to estimate node divergence times. Phylogenetic reconstructions using maximum parsimony, maximum likelihood, and Bayesian analysis retrieved the same clades, only branching order within these clades varied slightly between trees. Although the morphological diversity is remarkable, the endemic cyprinid taxa in East Asia emerged as a monophyletic clade referred to as Xenocypridini. The monophyly for the subfamilies including Cyprininae and Leuciscinae, as well as the tribes including Labeonini, Gobionini, Acheilognathini, and Leuciscini, was also well resolved with high nodal support. Analysis of the RAG2 gene supported the following cyprinid molecular phylogeny: the Danioninae is the most basal subfamily within the family Cyprinidae and the Cyprininae is the sister group of the Leuciscinae. The divergence times were estimated for the nodes corresponding to the principal clades within the Cyprinidae. The family Cyprinidae appears to have originated in the mid-Eocene in Asia, with the cladogenic event of the key basal group Danioninae occurring in the early Oligocene (about 31-30 MYA), and the origins of the two subfamilies, Cyprininae and Leuciscinae, occurring in the mid-Oligocene (around 26 MYA).     

Complete mitochondrial genome of the grass carp (Ctenopharyngodon idella, Teleostei): insight into its phylogenic position within Cyprinidae

Genome comparison has shed light on many fields of both basic and applied research, including the study of species phylogeny. Grass carp (Ctenopharyngodon idella) belongs to Cyprinidae, the largest freshwater fish family; but which subfamily it belongs to remains a controversial issue. In this study, the complete mitochondrial genome (mitogenome) sequence of grass carp was determined and phylogenetic analyses of all mitochondrial protein-coding genes and a nuclear gene (RAG 2) were conducted to explore the evolutionary relationship of grass carp with other cyprinid species. The mitogenome of grass carp is 16,609 bp in length. As with most other vertebrates, it contains the same gene order and an identical number of genes or regions, including 13 protein-coding genes, two rRNA genes, 22 tRNA genes and one putative control region. Phylogenetic analyses using two different datasets (mitochondrial and nuclear) and three different computational algorithms (Bayesian, MP and ML) all revealed two distinct groups with high statistical support, indicating that Cyprininae and Leuciscinae are two separate, valid subfamilies. Importantly, our phylogenetic result provides strong molecular evidence in support of the placement of Ctenopharyngodon in Leuciscinae rather than in Cyprininae.     

Molecular phylogenetics of cixiid planthoppers (Hemiptera: Fulgoromorpha): new insights from combined analyses of mitochondrial and nuclear genes

The planthopper family Cixiidae (Hemiptera: Fulgoromorpha) comprises approximately 160 genera and 2000 species divided in three subfamilies: Borystheninae, Bothriocerinae and Cixiinae, the later with 16 tribes. The current paper represents the first attempt to estimate phylogenetic relationships within Cixiidae based on molecular data. We use a total of 3652 bp sequence alignment of four genes: the mitochondrial coding genes Cytochrome c Oxidase subunit 1 (Cox1) and Cytochrome b (Cytb), a portion of the nuclear 18S rDNA and two non-contiguous portions of the nuclear 28S rDNA. The phylogenetic relationships of 72 terminal specimens were reconstructed using both maximum parsimony and Bayesian inference methods. Through the analysis of this empirical dataset, we also provide comparisons among different a priori partitioning strategies and the use of mixture models in a Bayesian framework. Our comparisons suggest that mixture models overcome the benefits obtained by partitioning the data according to codon position and gene identity, as they provide better accuracy in phylogenetic reconstructions. The recovered maximum parsimony and Bayesian inference phylogenies suggest that the family Cixiidae is paraphyletic in respect with Delphacidae. The paraphyly of the subfamily Cixiinae is also recovered by both approaches. In contrast to a morphological phylogeny recently proposed for cixiids, subfamilies Borystheninae and Bothriocerinae form a monophyletic group.     

Phylogeny of the gudgeons (Teleostei: Cyprinidae: Gobioninae)

The members of the cyprinid subfamily Gobioninae, commonly called gudgeons, form one of the most well-established assemblages in the family Cyprinidae. The subfamily is a species-rich group of fishes, these fishes display diverse life histories, appearances, and behavior. The phylogenetic relationships of Gobioninae are examined using sequence data from four loci: cytochrome b, cytochrome c oxidase I, opsin, and recombination activating gene 1. This investigation produced a data matrix of 4114 bp for 162 taxa that was analyzed using parsimony, maximum likelihood, and Bayesian inference methods. The phylogenies our analyses recovered corroborate recent studies on the group. The subfamily Gobioninae is monophyletic and composed of three major lineages. We find evidence for a Hemibarbus-Squalidus group, and the tribes Gobionini and Sarcocheilichthyini, with the Hemibarbus-Squalidus group sister to a clade of Gobionini-Sarcocheilichthyini. The Hemibarbus-Squalidus group includes those two genera; the tribe Sarcocheilichthyini includes Coreius, Coreoleuciscus, Gnathopogon, Gobiocypris, Ladislavia, Paracanthobrama, Pseudorasbora, Pseudopungtungia, Pungtungia, Rhinogobio, and Sarcocheilichthys; the tribe Gobionini includes Abbottina, Biwia, Gobio, Gobiobotia, Huigobio, Microphysogobio, Platysmacheilus, Pseudogobio, Romanogobio, Saurogobio, and Xenophysogobio. The monotypic Acanthogobio is placed into the synonymy of Gobio. We tentatively assign Belligobio to the Hemibarbus-Squalidus group and Mesogobio to Gobionini; Paraleucogobio and Parasqualidus remain incertae sedis. Based on the topologies presented, the evolution of swim bladder specializations, a distinctive feature among cyprinids, has occurred more than once within the subfamily.     

Ancestral sequence reconstruction in primate mitochondrial DNA: compositional bias and effect on functional inference

Reconstruction of ancestral DNA and amino acid sequences is an important means of inferring information about past evolutionary events. Such reconstructions suggest changes in molecular function and evolutionary processes over the course of evolution and are used to infer adaptation and convergence. Maximum likelihood (ML) is generally thought to provide relatively accurate reconstructed sequences compared to parsimony, but both methods lead to the inference of multiple directional changes in nucleotide frequencies in primate mitochondrial DNA (mtDNA). To better understand this surprising result, as well as to better understand how parsimony and ML differ, we constructed a series of computationally simple "conditional pathway" methods that differed in the number of substitutions allowed per site along each branch, and we also evaluated the entire Bayesian posterior frequency distribution of reconstructed ancestral states. We analyzed primate mitochondrial cytochrome b (Cyt-b) and cytochrome oxidase subunit I (COI) genes and found that ML reconstructs ancestral frequencies that are often more different from tip sequences than are parsimony reconstructions. In contrast, frequency reconstructions based on the posterior ensemble more closely resemble extant nucleotide frequencies. Simulations indicate that these differences in ancestral sequence inference are probably due to deterministic bias caused by high uncertainty in the optimization-based ancestral reconstruction methods (parsimony, ML, Bayesian maximum a posteriori). In contrast, ancestral nucleotide frequencies based on an average of the Bayesian set of credible ancestral sequences are much less biased. The methods involving simpler conditional pathway calculations have slightly reduced likelihood values compared to full likelihood calculations, but they can provide fairly unbiased nucleotide reconstructions and may be useful in more complex phylogenetic analyses than considered here due to their speed and flexibility. To determine whether biased reconstructions using optimization methods might affect inferences of functional properties, ancestral primate mitochondrial tRNA sequences were inferred and helix-forming propensities for conserved pairs were evaluated in silico. For ambiguously reconstructed nucleotides at sites with high base composition variability, ancestral tRNA sequences from Bayesian analyses were more compatible with canonical base pairing than were those inferred by other methods. Thus, nucleotide bias in reconstructed sequences apparently can lead to serious bias and inaccuracies in functional predictions.     

Resolving deep phylogenetic relationships in salamanders: analyses of mitochondrial and nuclear genomic data

Phylogenetic relationships among salamander families illustrate analytical challenges inherent to inferring phylogenies in which terminal branches are temporally very long relative to internal branches. We present new mitochondrial DNA sequences, approximately 2,100 base pairs from the genes encoding ND1, ND2, COI, and the intervening tRNA genes for 34 species representing all 10 salamander families, to examine these relationships. Parsimony analysis of these mtDNA sequences supports monophyly of all families except Proteidae, but yields a tree largely unresolved with respect to interfamilial relationships and the phylogenetic positions of the proteid genera Necturus and Proteus. In contrast, Bayesian and maximum-likelihood analyses of the mtDNA data produce a topology concordant with phylogenetic results from nuclear-encoded rRNA sequences, and they statistically reject monophyly of the internally fertilizing salamanders, suborder Salamandroidea. Phylogenetic simulations based on our mitochondrial DNA sequences reveal that Bayesian analyses outperform parsimony in reconstructing short branches located deep in the phylogenetic history of a taxon. However, phylogenetic conflicts between our results and a recent analysis of nuclear RAG-1 gene sequences suggest that statistical rejection of a monophyletic Salamandroidea by Bayesian analyses of our mitochondrial genomic data is probably erroneous. Bayesian and likelihood-based analyses may overestimate phylogenetic precision when estimating short branches located deep in a phylogeny from data showing substitutional saturation; an analysis of nucleotide substitutions indicates that these methods may be overly sensitive to a relatively small number of sites that show substitutions judged uncommon by the favored evolutionary model.     

Cyprinid phylogeny based on Bayesian and maximum likelihood analyses of partitioned data: implications for Cyprinidae systematics

Cyprinidae is the biggest family of freshwater fish, but the phylogenetic relationships among its higher-level taxa are not yet fully resolved. In this study, we used the nuclear recombination activating gene 2 and the mitochondrial 16S ribosomal RNA and cytochrome b genes to reconstruct cyprinid phylogeny. Our aims were to (i) demonstrate the effects of partitioned phylogenetic analyses on phylogeny reconstruction of cyprinid fishes; (ii) provide new insights into the phylogeny of cyprinids. Our study indicated that unpartitioned strategy was optimal for our analyses; partitioned analyses did not provide better-resolved or -supported estimates of cyprinid phylogeny. Bayesian analyses support the following relationships among the major monophyletic groups within Cyprinidae: (Cyprininae, Labeoninae), ((Acheilognathinae, ((Leuciscinae, Tincinae), Gobioninae)), Xenocyprininae). The placement of Danioninae was poorly resolved. Estimates of divergence dates within the family showed that radiation of the major cyprinid groups occurred during the Late Oligocene through the Late Miocene. Our phylogenetic analyses improved our understanding of the evolutionary history of this important fish family.     

Bias and conflict in phylogenetic inference of myco-heterotrophic plants: a case study in Thismiaceae

Due to morphological reduction and absence of amplifiable plastid genes, the identification of photosynthetic relatives of heterotrophic plants is problematic. Although nuclear and mitochondrial gene sequences may offer a welcome alternative source of phylogenetic markers, the presence of rate heterogeneity in these genes may introduce bias/systematic error in phylogenetic analyses. We examine the phylogenetic position of Thismiaceae based on nuclear 18S rDNA and mitochondrial atpA DNA sequence data, as well as using parsimony, likelihood and Bayesian inference methods. Significant differences in evolutionary rates of these genes between closely related taxa lead to conflicting results: while parsimony analyses of 18S rDNA and combined data strongly support the monophyly of Thismiaceae, Bayesian inference, with and without a relaxed molecular clock, as well as the Swofford–Olsen–Waddell–Hillis (SOWH) test confidently reject this hypothesis. We show that rate heterogeneity in our data leads to long-branch attraction artifacts in parsimony analysis. However, using model-based inference methods the question of whether Thismiaceae are monophyletic remains elusive. On the one hand maximum likelihood nonparametric bootstrapping and parametric hypothesis tests fail to support a paraphyletic Thismiaceae, on the other hand Bayesian inference methods (both without and with a relaxed clock) significantly reject a monophyletic Thismiaceae. These results show that an adequate sampling, the use of rate homogeneous data, and the application of different inference methods are important factors for developing phylogenetic hypotheses of myco-heterotrophic plants. © The Willi Hennig Society 2009.     

Insights into the phylogenetic relationships within Cixiidae (Hemiptera: Fulgoromorpha): cladistic analysis of a morphological dataset

Abstract.  According to the most recent classifications proposed, the planthopper family Cixiidae comprises three subfamilies, namely Borystheninae, Bothriocerinae and Cixiinae, the latter with 16 tribes. Here we examine morphological characters to present the first phylogenetic reconstructions within Cixiidae derived from a cladistic analysis. We scored 85 characters of the head, thorax, and male and female genitalia for 50 taxa representative of all cixiid subfamilies and tribes and for six outgroup taxa. Analyses were based on maximum parsimony – using both equally weighted and successive weighting procedures – and Bayesian inferences. The monophyly of most currently accepted tribes and subfamilies was investigated through Templeton statistical tests of alternative phylogenetic hypotheses. The cladistic analyses recover the monophyly of Cixiidae, the subfamily Bothriocerinae, and the tribes Pentastirini, Mnemosynini, and Eucarpiini. Successive weighting and Bayesian inference recover the monophyly of the tribe Gelastocephalini, but only Bayesian inference supports the monophyly of Semoniini. The relationships recovered support the groups [Stenophlepsini (Borystheninae + Bothriocerinae)] arising from the tribe Oecleini, and [Andini + Brixiidini + Brixiini (polyphyletic) + Bennini]. Templeton tests reject the alternative hypothesis of a monophyletic condition for the tribe Pintaliini as presently defined.     

Application of the secondary structure model of rRNA for phylogeny: D2-D3 expansion segments of the LSU gene of plant-parasitic nematodes from the family Hoplolaimidae Filipjev, 1934

Knowledge of rRNA structure is increasingly important to assist phylogenetic analysis through reconstructing optimal alignment, utilizing molecule features as an additional source of data and refining appropriate models of evolution of the molecule. We describe a procedure of optimization for alignment and a new coding method for nucleotide sequence data using secondary structure models of the D2 and D3 expansion fragments of the LSU-rRNA gene reconstructed for fifteen nematode species of the agriculturally important and diverse family Hoplolaimidae, order Tylenchida. Using secondary structure information we converted the original sequence data into twenty-eight symbol codes and submitted the transformed data to maximum parsimony analysis. We also applied the original sequence data set for Bayesian inference. This used the doublet model with sixteen states of nucleotide doublets for the stem region and the standard model of DNA substitution with four nucleotide states for loops and bulges. By this approach, we demonstrate that using structural information for phylogenetic analyses led to trees with lower resolved relationships between clades and likely eliminated some artefactual support for misinterpreted relationships, such as paraphyly of Helicotylenchus or Rotylenchus. This study as well as future phylogenetic analyses is herein supported by the development of an on-line database, NEMrRNA, for rRNA molecules in a structural format for nematodes. We also have developed a new computer program, RNAstat, for calculation of nucleotide statistics designed and proposed for phylogenetic studies.     

Evolution and trans-species polymorphism of MHC class IIbeta genes in cyprinid fish

The polymorphism of DAB genes encoding MHC IIbeta was investigated in 11 cyprinid species from central Europe. The species belonged to four subfamilies: Cyprininae, Tincinae, Gobioninae and Leuciscinae. Two paralogous groups of sequences, DAB1 and DAB3, were recognised according to the similarity of their nucleotide and amino-acid sequences and from phylogenetic analyses using either partial exon 2 or partial exon 3 sequences. A high allelic variability among species was found for exon 2, indicating extensive MHC polymorphism. Time divergence estimation supports the separation of DAB1 and DAB3 groups predating the separation into fish subfamilies, and a cyprinid origin of the DAB genes. Phylogenetic trees using exon 2 support the hypothesis of trans-species polymorphism, which appears to be limited to the subfamily level, i.e. the presence of sequences from different species in the same allelic group was more often recognised within subfamilies Cyprininae and Leuciscinae than between them. Phylogenetic trees using exon 3 reflect the phylogenetic patterns previously found for Cyprinidae systematics. Specific nucleotides and amino-acids in exon 3 that separate both subfamilies, as well as the species within the Cyprininae subfamily were observed. A lack of segregation in leuciscin species was recognised and the alleles of different leuciscin species tend to share similar motifs in exon 3. This could be explained by the ancient and complicated dispersion history of Cyprininae and the radiation of Leuciscinae. The effects of selective pressures were investigated: (1) within species, (2) among lineages, and (3) among sites. From intraspecific analyses, exon 2 sequences were identified as the targets of diversifying selection, whilst the evolution of exon 3 seems to be under the influence of purifying selection. The analyses among lineages indicate positive selection in many branches when using exon 2, therefore confirming trans-species polymorphism, whilst the DAB lineages of exon 3 are potentially submitted to purifying selection to some extent. Moreover, our results suggest the secondary acquisition of function of DAB1 group after duplication. The analyses among sites reveal that exon 2 exhibits sites under positive selection mostly corresponding to the putative PBR sites involved in the alpha-helix structure of the protein.     

Probabilistic methods surpass parsimony when assessing clade support in phylogenetic analyses of discrete morphological data

Fossil taxa are critical to inferences of historical diversity and the origins of modern biodiversity, but realizing their evolutionary significance is contingent on restoring fossil species to their correct position within the tree of life. For most fossil species, morphology is the only source of data for phylogenetic inference; this has traditionally been analysed using parsimony, the predominance of which is currently challenged by the development of probabilistic models that achieve greater phylogenetic accuracy. Here, based on simulated and empirical datasets, we explore the relative efficacy of competing phylogenetic methods in terms of clade support. We characterize clade support using bootstrapping for parsimony and Maximum Likelihood, and intrinsic Bayesian posterior probabilities, collapsing branches that exhibit less than 50% support. Ignoring node support, Bayesian inference is the most accurate method in estimating the tree used to simulate the data. After assessing clade support, Bayesian and Maximum Likelihood exhibit comparable levels of accuracy, and parsimony remains the least accurate method. However, Maximum Likelihood is less precise than Bayesian phylogeny estimation, and Bayesian inference recaptures more correct nodes with higher support compared to all other methods, including Maximum Likelihood. We assess the effects of these findings on empirical phylogenies. Our results indicate probabilistic methods should be favoured over parsimony.     

Phylogenetic relationships of tachinid flies in subfamily Exoristinae (Tachinidae: Diptera) based on 28S rDNA and elongation factor-1α

Abstract The phylogenetic relationships within the largest subfamily of Tachinidae, Exoristinae, were explored using nucleotide sequences of two genes (EF-1α and 28S rDNA). A total of fifty-five and forty-three taxa were represented in the analyses for each gene, respectively, representing forty-three genera. Neighbour joining, parsimony and maximum likelihood inference methods were employed to reconstruct phylogenetic relationships in separate analyses of each gene, and parsimony was used to analyse the combined dataset. Although certain taxa were highly mobile, phylogenetic reconstructions generally supported recent classification schemes based on reproductive habits and genitalia. Generally, the monophyly of Tachinidae and Exoristinae was supported. Tribes Winthemiini, Exoristini and Blondeliini were repeatedly constructed as monophyletic groups, with the former two clades often occupying a basal position among Exoristinae. Goniini and Eryciini generally clustered together as a derived clade within Exoristinae; however, they were never reconstructed as two distinct clades. These results suggest that the possession of unembryonated eggs is plesiomorphic within the subfamily and that there may have been multiple transitions between microtype and macrotype egg forms.     

Phylogenetic relationships of Iberian cyprinids: systematic and biogeographical implications.

The phylogenetic relationships among all Iberian endemic cyprinids were inferred using the complete nucleotide sequence of the cytochrome b gene. The inferred molecular phylogeny included representatives from Central European, Asian and North African species, and is highly congruent with previous phylogenies based on osteological characters. Iberian cyprinids were grouped into only five, very speciose lineages (with the exception of the monotypic Anaecypris): Barbus, Luciobarbus, Chondrostoma, Leuciscus and Anaecypris. The existence of such a relatively small number of Iberian cyprinid lineages can be explained by the historical isolation of the Iberian Peninsula. North African and Asian barbels are the sister group of Iberian Luciobarbus, supporting a south-eastern route of colonization of the Iberian Peninsula for this subgenus. With leuciscins, Anaecypris hispanica was considered a relict species as it could not be related to any other Iberian cyprinid. The phylogenetic relationships among the main lineages of Iberian cyprinids based on cytochrome b sequence data supported the traditional division of the Cyprinidae into two subfamilies: Cyprininae and Leuciscinae.     

Large multi-gene phylogenetic trees of the grasses (Poaceae): progress towards complete tribal and generic level sampling

In this paper we included a very broad representation of grass family diversity (84% of tribes and 42% of genera). Phylogenetic inference was based on three plastid DNA regions rbcL, matK and trnL-F, using maximum parsimony and Bayesian methods. Our results resolved most of the subfamily relationships within the major clades (BEP and PACCMAD), which had previously been unclear, such as, among others the: (i) BEP and PACCMAD sister relationship, (ii) composition of clades and the sister-relationship of Ehrhartoideae and Bambusoideae + Pooideae, (iii) paraphyly of tribe Bambuseae, (iv) position of Gynerium as sister to Panicoideae, (v) phylogenetic position of Micrairoideae. With the presence of a relatively large amount of missing data, we were able to increase taxon sampling substantially in our analyses from 107 to 295 taxa. However, bootstrap support and to a lesser extent Bayesian inference posterior probabilities were generally lower in analyses involving missing data than those not including them. We produced a fully resolved phylogenetic summary tree for the grass family at subfamily level and indicated the most likely relationships of all included tribes in our analysis.     

Angiosperm phylogeny based on matK sequence information

Plastid matK gene sequences for 374 genera representing all angiosperm orders and 12 genera of gymnosperms were analyzed using parsimony (MP) and Bayesian inference (BI) approaches. Traditionally, slowly evolving genomic regions have been preferred for deep-level phylogenetic inference in angiosperms. The matK gene evolves approximately three times faster than the widely used plastid genes rbcL and atpB. The MP and BI trees are highly congruent. The robustness of the strict consensus tree supercedes all individual gene analyses and is comparable only to multigene-based phylogenies. Of the 385 nodes resolved, 79% are supported by high jackknife values, averaging 88%. Amborella is sister to the remaining angiosperms, followed by a grade of Nymphaeaceae and Austrobaileyales. Bayesian inference resolves Amborella + Nymphaeaceae as sister to the rest, but with weak (0.42) posterior probability. The MP analysis shows a trichotomy sister to the Austrobaileyales representing eumagnoliids, monocots + Chloranthales, and Ceratophyllum + eudicots. The matK gene produces the highest internal support yet for basal eudicots and, within core eudicots, resolves a crown group comprising Berberidopsidaceae/Aextoxicaceae, Santalales, and Caryophyllales + asterids. Moreover, matK sequences provide good resolution within many angiosperm orders. Combined analyses of matK and other rapidly evolving DNA regions with available multigene data sets have strong potential to enhance resolution and internal support in deep level angiosperm phylogenetics and provide additional insights into angiosperm evolution.     

Uninode coding vs gene tree parsimony for phylogenetic reconstruction using duplicate genes

Two different methods of using paralogous genes for phylogenetic inference have been proposed: reconciled trees (or gene tree parsimony) and uninode coding. Gene tree parsimony suffers from 10 serious problems, including differential weighting of nucleotide and gap characters, undersampling which can be misinterpreted as synapomorphy, all of the characters not being allowed to interact, and conflict between gene trees being given equal weight, regardless of branch support. These problems are largely avoided by using uninode coding. The uninode coding method is elaborated to address multiple gene duplications within a single gene tree family and handle problems caused by lack of gene tree resolution. An example of vertebrate phylogeny inferred from nine genes is reanalyzed using uninode coding. We suggest that uninode coding be used instead of gene tree parsimony for phylogenetic inference from paralogous genes.     

Partial combination applied to phylogeny of European cyprinids using the mitochondrial control region

Previous molecular phylogenies of European cyprinids led to some solid facts and some uncertainties. This study is based on a stretch of more than 1 kb in the mitochondrial control region newly sequenced for 35 European cyprinids and on previous cytochrome b and 16S rDNA data. The trees based on the control region are more accurate and robust than those obtained from the two other genes. Character incongruence among the three genes was tested using the incongruence length difference (ILD) test. Iterative removals of individual sequences followed by new ILD tests identified two sequences responsible for statistically significant incongruence. A partial combination was conducted, that is, a combination of the three data sets, removing the two sequences previously identified. The phylogenetic analysis of this partial combination gives a more robust and resolved picture of subfamilial interrelationships. The Rasborinae are the sister group of all other cyprinids. The monophyletic Cyprininae emerges next. Tinca tinca first and then Rhodeus are the sister groups of all the remaining nonrasborine and noncyprinine species. Gobio is the sister group of the Leuciscinae, in which the Phoxinini are the sister group of the Leuciscini. Within the Leuciscini, the genus Leuciscus and the subfamily Alburninae are both paraphyletic. The Rasborinae are the most basal cyprinid subfamily and the Tincinae are not the sister group of the Cyprininae. These two results challenge only two anatomical characters, which need to be reinterpreted or considered as homoplastic in cyprinid evolution: the modification of the first pleural rib and its parapophysis and the bony composition of the interorbital septum.     

Molecular systematics of the Eastern Fence Lizard (Sceloporus undulatus): a comparison of Parsimony, Likelihood, and Bayesian approaches

Phylogenetic analysis of large datasets using complex nucleotide substitution models under a maximum likelihood framework can be computationally infeasible, especially when attempting to infer confidence values by way of nonparametric bootstrapping. Recent developments in phylogenetics suggest the computational burden can be reduced by using Bayesian methods of phylogenetic inference. However, few empirical phylogenetic studies exist that explore the efficiency of Bayesian analysis of large datasets. To this end, we conducted an extensive phylogenetic analysis of the wide-ranging and geographically variable Eastern Fence Lizard (Sceloporus undulatus). Maximum parsimony, maximum likelihood, and Bayesian phylogenetic analyses were performed on a combined mitochondrial DNA dataset (12S and 16S rRNA, ND1 protein-coding gene, and associated tRNA; 3,688 bp total) for 56 populations of S. undulatus (78 total terminals including other S. undulatus group species and outgroups). Maximum parsimony analysis resulted in numerous equally parsimonious trees (82,646 from equally weighted parsimony and 335 from weighted parsimony). The majority rule consensus tree derived from the Bayesian analysis was topologically identical to the single best phylogeny inferred from the maximum likelihood analysis, but required approximately 80% less computational time. The mtDNA data provide strong support for the monophyly of the S. undulatus group and the paraphyly of "S. undulatus" with respect to S. belli, S. cautus, and S. woodi. Parallel evolution of ecomorphs within "S. undulatus" has masked the actual number of species within this group. This evidence, along with convincing patterns of phylogeographic differentiation suggests "S. undulatus" represents at least four lineages that should be recognized as evolutionary species.     

Phylogeny and biogeography of the family Cyprinidae in the Middle East inferred from cytochrome b DNA- evolutionary significance of this region.

The phylogenetic relationships of cyprinid species from the Middle East and neighboring biogeographical areas were investigated using cytochrome b sequence variation in order to test hypotheses that consider the Middle Eastern area as an important interchange area or a center of speciation for the freshwater fauna. A total of 62 cyprinid species were analyzed over the complete cytochrome b fragment (1140 bp); 28 belong to the Leuciscinae subfamily and 34 to the Cyprininae. All the Leuciscinae lineage fish recorded in the Middle East are also found in Europe, which was interpreted as an important Palearctic influence in the Middle Eastern ichthyofauna consistent with the Lago Mare dispersion. However, it has also been suggested that several Danube species have their origins in the Middle East. In contrast, the Cyprininae subfamily showed three highly divergent lineages, one shared with the Euro-Mediterranean area (Barbus/Luciobarbus genus) relict of the Lago Mare dispersion, one shared with Africa (Carasobarbus/Varicorhinus subgenus), and the third shared with Asia (Garra genus). Furthermore, clades observed in the phylogenetic reconstructions are not consistent with morphometric or karyological data and disagree with previous taxonomic assumptions. Lastly, the dispersion history in the Middle East of this subfamily appears much more complicated and ancient than that of the Leuciscinae. However, taking into account Cyprininae and Leuciscinae distribution, the Middle East appears more like an important interchange area for the freshwater ichtyofauna than a center of speciation.