Molecular Phylogenetics at the Felsenstein Zone: Approaching the Strepsiptera Problem Using 5.8S and 28S rDNA Sequences

Recent efforts to reconstruct the phylogenetic position of the insect order Strepsiptera have elicited a major controversy in molecular phylogenetics. We sequenced the 5.8S rDNA and major parts of the 28S rDNA 5′ region of the strepsipteran speciesStylops melittae.Their evolutionary dynamics were analyzed together with previously published insect rDNA sequences to identify tree estimation bias risks and to explore additional sources of phylogenetic information. Several major secondary structure changes were found as being autapomorphic for the Diptera, the Strepsiptera, or the Archaeognatha. Besides elevated substitution rates a significant AT bias was present in dipteran and strepsipteran 28S rDNA which, however, was restricted to stem sites in the Diptera while also affecting single-stranded sites in the Strepsiptera. When dipteran taxa were excluded from tree estimation all methods consistently supported the placement of Strepsiptera to within the Holometabola. When dipteran taxa were included maximum likelihood continued to favor a sister-group relationship of Strepsiptera with Mecoptera while remaining methods strongly supported a sister-group relationship with Diptera. Parametric bootstrap analysis revealed maximum likelihood as a consistent estimator if rate heterogeneity across sites was taken into account. Though the position of Strepsiptera within Holometabola remains elusive, we conclude that the evolution of dipteran and strepsipteran rDNA involved similar yet independent changes of substitution parameters.     

The phylogenetic relationships of flies in the superfamily Empidoidea (Insecta: Diptera)

We conducted a molecular phylogenetic study of the Empidoidea, a diverse group of 10,000 species of true flies, with two major goals: to reconstruct a taxonomically complete and robustly supported phylogeny for the group and to use this information to assess several competing classifications for the clade. We amassed 3900+ nucleotides of coding data from the carbamoylphosphate synthase domain of the rudimentary locus (CAD) and 1200+ nucleotides from the large nuclear ribosomal subunit (28S) from 72 and 71 species, respectively, representing several orthorrhaphan and cyclorrhaphan families and all previously recognized empidoidean subfamilies. Independent and combined phylogenetic analyses of these data were conducted using parsimony, maximum likelihood, and Bayesian criteria. The combined matrix included 61 taxa for which both CAD and 28S sequences were obtained. Analyses of CAD first and second codon positions alone and when concatenated with 28S sequences yielded trees with similar and largely stable topologies. Analyses of 28S data alone supported many clades although resolution is limited by low sequence divergence. The following major empidoid clades were recovered with convincing support in a majority of analyses: Atelestidae, Empidoidea exclusive of Atelestidae, Hybotidae sensu lato, Dolichopodidae+Microphorinae (including Parathallassius), and Empididae sensu lato (including Brachystomatinae, Ceratomerinae, Clinocerinae, Empidinae, Hemerodromiinae, Oreogetoninae, and Trichopezinae). The branching arrangement among these four major clades was Atelestidae, Hybotidae, Dolichopodidae/Microphorinae, Empididae. Previously recognized subclades recovered with robust support included Hybotinae, Brachystomatinae, Tachydromiinae, Clinocerinae (in part), Hemerodromiinae, Empidinae, and Empidiini.     

CHS基因起源初探及其在被子植物中的进化分析

利用PCR与TAIL-PCR方法,从半月苔(Lunulariacruciata(L.)Dum.exLindb.)中获得了一段长约1000bp的基因片段,它与已知的CHS基因在核苷酸水平上的相似性大于56%,在氨基酸水平上的相似性大于60%,所推断的氨基酸序列中酶反应的4个催化位点与已知晶体结构的紫花苜蓿MCHS2A上的催化位点相同,首次证明了苔类植物中可能存在类CHS基因,将CHS基因的起源时间推到苔藓类植物出现之前。以该序列和两种蕨类植物(Psilotumnudum(L.)Griseb.和EquisetumarvenseL.)的CHS序列作为外类群,应用邻接法、最大简约法和最大似然法分别构建了被子植物的CHS的分子系统树。结果表明,大部分科中的CHS分布在不同的分支上,而十字花科、豆科和禾本科各自聚成一个单系类群。以邻接树为依据,对茄科、旋花科和菊科的CHS基因进行了相对碱基替换速率的检测,发现这三个科内或科间序列的替换速率不一致。被子植物的CHS基因在基因拷贝数目、碱基替换速率以及重复/丢失事件的发生上都存在较大的差异,这种差异可能与被子植物的生活史、生活环境、花的特性以及对外界的防御系统等的多样性相关。     

CHS基因起源初探及其在被子植物中的进化分析

利用PCR与TAlL-PCR方法,从半月苔(Lunularia cructata(L.)Dum.ex Lindb)中获得了一段长约l 000 bp的基因片段,它与已知的CHS基因在核苷酸水平上的相似性大于56%,在氨基酸水平上的相似性大于60%,所推断的氨基酸序列中酶反应的4个催化位点与已知晶体结构的紫花苜蓿MCHS2A上的催化位点相同,首次证明了苔类植物中可能存在类CHS基因,将CHS基因的起源时间推到苔藓类植物出现之前.以该序列和两种蕨类植物(Psilotumnudum(L.)Griseb.和Equisetum arvense L.)的CHS序列作为外类群,应用邻接法、最大简约法和最大似然法分别构建了被子植物的CHS的分子系统树.结果表明,大部分科中的CHS分布在不同的分支上,而十字花科、可科和禾本科各自聚成一个单系类群.以邻接树为依据,对茄科、旋花科和菊科的CHS基因进行了相对碱基替换速率的检测,发现这三个科内或科间序列的替换速率不一致.被子植物的CHS基因在基因拷贝数目、碱基替换速率以及重复/丢失事件的发生上都存在较大的差异,这种差异可能与被子植物的生活史、生活环境、花的特性以及对外界的防御系统等的多样性相关.     

Spatial and temporal heterogeneity in nucleotide sequence evolution

Models of nucleotide substitution make many simplifying assumptions about the evolutionary process, including that the same process acts on all sites in an alignment and on all branches on the phylogenetic tree. Many studies have shown that in reality the substitution process is heterogeneous and that this variability can introduce systematic errors into many forms of phylogenetic analyses. I propose a new rigorous approach for describing heterogeneity called a temporal hidden Markov model (THMM), which can distinguish between among site (spatial) heterogeneity and among lineage (temporal) heterogeneity. Several versions of the THMM are applied to 16 sets of aligned sequences to quantitatively assess the different forms of heterogeneity acting within them. The most general THMM provides the best fit in all the data sets examined, providing strong evidence of pervasive heterogeneity during evolution. Investigating individual forms of heterogeneity provides further insights. In agreement with previous studies, spatial rate heterogeneity (rates across sites [RAS]) is inferred to be the single most prevalent form of heterogeneity. Interestingly, RAS appears so dominant that failure to independently include it in the THMM masks other forms of heterogeneity, particularly temporal heterogeneity. Incorporating RAS into the THMM reveals substantial temporal and spatial heterogeneity in nucleotide composition and bias toward transition substitution in all alignments examined, although the relative importance of different forms of heterogeneity varies between data sets. Furthermore, the improvements in model fit observed by adding complexity to the model suggest that the THMMs used in this study do not capture all the evolutionary heterogeneity occurring in the data. These observations all indicate that current tests may consistently underestimate the degree of temporal heterogeneity occurring in data. Finally, there is a weak link between the amount of heterogeneity detected and the level of divergence between the sequences, suggesting that variability in the evolutionary process will be a particular problem for deep phylogeny.     

Molecular phylogeny, historical biogeography, and divergence time estimates for swallowtail butterflies of the genus Papilio (Lepidoptera: Papilionidae)

Swallowtail butterflies are recognized as model organisms in ecology, evolutionary biology, genetics, and conservation biology but present numerous unresolved phylogenetic problems. We inferred phylogenetic relationships for 51 of about 205 species of the genus Papilio (sensu lato) from 3.3-Kilobase (kb) sequences of mitochondrial and nuclear DNA (2.3 kb of cytochrome oxidases I and II and 1.0 kb of elongation factor 1 alpha). Congruent phylogenetic trees were recovered within Papilio from analyses of combined data using maximum likelihood, Bayesian analysis, and maximum parsimony bootstrap consensus. Several disagreements with the traditional classification of Papilio were found. Five major previously hypothesized subdivisions within Papilio were well supported: Heraclides, Pterourus, Chilasa, Papilio (sensu stricto), and Eleppone. Further studies are required to clarify relationships within traditional "Princeps," which was paraphyletic. Several biologically interesting characteristics of Papilio appear to have polyphyletic origins, including mimetic adults, larval host associations, and larval morphology. Early diversification within Papilio is estimated at 55-65 million years ago based on a combination of biogeographic time constraints rather than fossils. This divergence time suggests that Papilio has slower apparent substitution rates than do Drosophila and fig-pollinating wasps and/or divergences corrected using best-fit substitution models are still being consistently underestimated. The amount of sequence divergence between Papilio subdivisions is equivalent to divergences between genera in other tribes of the Papilionidae, and between genera of moths of the noctuid subfamily Heliothinae.     

Dynamically heterogenous partitions and phylogenetic inference: an evaluation of analytical strategies with cytochrome b and ND6 gene sequences in cranes

ki ctes over whether molecular sequence data should be partitioned for phylogenetic analysis often confound two types of heterogeneity among partitions. We distinguish historical heterogeneity (i.e., different partitions have different evolutionary relationships) from dynamic heterogeneity (i.e., different partitions show different patterns of sequence evolution) and explore the impact of the latter on phylogenetic accuracy and precision with a two-gene, mitochondrial data set for cranes. The well-established phylogeny of cranes allows us to contrast tree-based estimates of relevant parameter values with estimates based on pairwise comparisons and to ascertain the effects of incorporating different amounts of process information into phylogenetic estimates. We show that codon positions in the cytochrome b and NADH dehydrogenase subunit 6 genes are dynamically heterogenous under both Poisson and invariable-sites + gamma-rates versions of the F84 model and that heterogeneity includes variation in base composition and transition bias as well as substitution rate. Estimates of transition-bias and relative-rate parameters from pairwise sequence comparisons were comparable to those obtained as tree-based maximum likelihood estimates. Neither rate-category nor mixed-model partitioning strategies resulted in a loss of phylogenetic precision relative to unpartitioned analyses. We suggest that weighted-average distances provide a computationally feasible alternative to direct maximum likelihood estimates of phylogeny for mixed-model analyses of large, dynamically heterogenous data sets.     

Evolutionary stability of the R1 retrotransposable element in the genus Drosophila

R1 is a non-long terminal repeat (non-LTR) retrotransposable element that inserts into a specific sequence of insect 28S ribosomal RNA genes. We have previously shown that this element has been maintained through vertical transmission in the melanogaster species subgroup of Drosophila. To address whether R1 elements have been vertically transmitted for longer periods of evolutionary time, the analysis has been extended to 11 other species from four species groups of the genus Drosophila (melanogaster, obscura, testecea, and repleta). All sequenced elements appeared functional on the basis of the preservation of their open-reading frames and consistently higher rate of substitution at synonymous sites relative to replacement sites. The phylogenetic relationships of the R1 elements from all species analyzed were congruent with the species phylogenies, suggesting that the R1 elements have been vertically transmitted since the inception of the Drosophila genus, an estimated 50-70 Mya. The stable maintenance of R1 through the germ line appears to be the major mechanism for the widespread distribution of these elements in Drosophila. In two species, D. neotestecea of the testecea group and D. takahashii of the melanogaster group, a second family of R1 elements was also present that differed in sequence by 46% and 31%, respectively, from the family that was congruent with the species phylogeny. These second families may represent occasional horizontal transfers or, alternatively, they could reflect the ability of R1 elements to diverge into new families within a species and evolve independently.      

Heterogeneity of nucleotide frequencies among evolutionary lineages and phylogenetic inference

A major assumption of many molecular phylogenetic methods is the homogeneity of nucleotide frequencies among taxa, which refers to the equality of the nucleotide frequency bias among species. Changes in nucleotide frequency among different lineages in a data set are thought to lead to erroneous phylogenetic inference because unrelated clades may appear similar because of evolutionarily unrelated similarities in nucleotide frequencies. We tested the effects of the heterogeneity of nucleotide frequency bias on phylogenetic inference, along with the interaction between this heterogeneity and stratified taxon sampling, by means of computer simulations using evolutionary parameters derived from genomic databases. We found that the phylogenetic trees inferred from data sets simulated under realistic, observed levels of heterogeneity for mammalian genes were reconstructed with accuracy comparable to those simulated with homogeneous nucleotide frequencies; the results hold for Neighbor-Joining, minimum evolution, maximum parsimony, and maximum-likelihood methods. The LogDet distance method, specifically designed to deal with heterogeneous nucleotide frequencies, does not perform better than distance methods that assume substitution pattern homogeneity among sequences. In these specific simulation conditions, we did not find a significant interaction between phylogenetic accuracy and substitution pattern heterogeneity among lineages, even when the taxon sampling is increased.     

Phylogenetics and temporal diversification of the earliest true flies (Insecta: Diptera) based on multiple nuclear genes

Abstract Relationships among families of the lower Diptera (formerly suborder ‘Nematocera’) have been exceptionally difficult to resolve. Multiple hypotheses based on morphology have been proposed to identify the earliest lineages of flies and place the phylogenetic origin of the higher flies (Brachycera), but convincing support is limited. Here we resolve relationships among the major groups of lower Diptera using sequence data from four nuclear markers, including both ribosomal (28S rDNA) and protein‐coding (CAD, TPI and PGD) genes. Our results support both novel and traditional arrangements. Most unexpectedly, the small, highly‐specialized family Deuterophlebiidae appears to be sister to all remaining Diptera. Other results include the resolution of the traditional infra‐orders Culicomorpha (including a novel superfamily Simulioidea = Thaumaleidae + Simuliidae), Tipulomorpha (Tipulidae sensu lato + Trichoceridae) and Bibionomorpha sensu lato. We find support for a limited Psychodomorpha (Blephariceridae, Tanyderidae and Psychodidae) and Ptychopteromorpha (Ptychopteridae), whereas the placement of several enigmatic families (Nymphomyiidae, Axymyiidae and Perissommatidae) remains ambiguous. According to genetic data, the infra‐order Bibionomorpha is sister to the Brachycera. Much of the phylogenetic signal for major lineages was found in the 28S rDNA gene, whereas protein‐coding genes performed variably at different levels. In addition to elucidating relationships, we also estimate the age of major lower dipteran clades, based on molecular divergence time estimates using relaxed‐clock Bayesian methods and fossil calibration points.     

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.     

Molecular phylogeny of the carnivora (mammalia): assessing the impact of increased sampling on resolving enigmatic relationships

This study analyzed 76 species of Carnivora using a concatenated sequence of 6243 bp from six genes (nuclear TR-i-I, TBG, and IRBP; mitochondrial ND2, CYTB, and 12S rRNA), representing the most comprehensive sampling yet undertaken for reconstructing the phylogeny of this clade. Maximum parsimony and Bayesian methods were remarkably congruent in topologies observed and in nodal support measures. We recovered all of the higher level carnivoran clades that had been robustly supported in previous analyses (by analyses of morphological and molecular data), including the monophyly of Caniformia, Feliformia, Arctoidea, Pinnipedia, Musteloidea, Procyonidae + Mustelidae sensu stricto, and a clade of (Hyaenidae + (Herpestidae + Malagasy carnivorans)). All of the traditional "families," with the exception of Viverridae and Mustelidae, were robustly supported as monophyletic groups. We further have determined the relative positions of the major lineages within the Caniformia, which previous studies could not resolve, including the first robust support for the phylogenetic position of marine carnivorans (Pinnipedia) within the Arctoidea (as the sister-group to musteloids [sensu lato], with ursids as their sister group). Within the pinnipeds, Odobenidae (walrus) was more closely allied with otariids (sea lions/fur seals) than with phocids ("true" seals). In addition, we recovered a monophyletic clade of skunks and stink badgers (Mephitidae) and resolved the topology of musteloid interrelationships as: Ailurus (Mephitidae (Procyonidae, Mustelidae [sensu stricto])). This pattern of interrelationships of living caniforms suggests a novel inference that large body size may have been the primitive condition for Arctoidea, with secondary size reduction evolving later in some musteloids. Within Mustelidae, Bayesian analyses are unambiguous in supporting otter monophyly (Lutrinae), and in both MP and Bayesian analyses Martes is paraphyletic with respect to Gulo and Eira, as has been observed in some previous molecular studies. Within Feliformia, we have confirmed that Nandinia is the outgroup to all other extant feliforms, and that the Malagasy Carnivora are a monophyletic clade closely allied with the mongooses (Herpestidae [sensu stricto]). Although the monophyly of each of the three major feliform clades (Viverridae sensu stricto, Felidae, and the clade of Hyaenidae + (Herpestidae + Malagasy carnivorans)) is robust in all of our analyses, the relative phylogenetic positions of these three lineages is not resolvable at present. Our analyses document the monophyly of the "social mongooses," strengthening evidence for a single origin of eusociality within the Herpestidae. For a single caniform node, the position of pinnipeds relative to Ursidae and Musteloidea, parsimony analyses of data for the entire Carnivora did not replicate the robust support observed for both parsimony and Bayesian analyses of the caniform ingroup alone. More detailed analyses and these results demonstrate that outgroup choice can have a considerable effect on the strength of support for a particular topology. Therefore, the use of exemplar taxa as proxies for entire clades with diverse evolutionary histories should be approached with caution.The Bayesian analysis likelihood functions generally were better able to reconstruct phylogenetic relationships (increased resolution and more robust support for various nodes) than parsimony analyses when incompletely sampled taxa were included. Bayesian analyses were not immune, however, to the effects of missing data; lower resolution and support in those analyses likely arise from non-overlap of gene sequence data among less well-sampled taxa. These issues are a concern for similar studies, in which different gene sequences are concatenated in an effort to increase resolving power.     

Evolution and phylogenetic utility of CAD (rudimentary) among Mesozoic-aged Eremoneuran Diptera (Insecta)

We sequenced nearly the entire carbomoylphosphate synthase (CPS) domain of CAD, or rudimentary, (ca. 4 kb) from 29 species of flies representing all major clades within Eremoneura, or higher flies, and several orthorrhaphous brachyceran outgroups. We compared these sequences with orthologs from Anopheles gambiae and Drosophila melanogaster to assess structure, compositional bias, and phylogenetic utility. CAD is large (6.6+ kb), complex (comprised of three major and myriad minor functional domains) and relatively free of introns, extreme nucleotide bias (except third codon positions), and large hypervariable regions. The CPS domain possesses moderate levels of nonsynonymous divergence among taxa of intermediate evolutionary age and conveys considerable phylogenetic signal. Phylogenetic analysis of CPS sequences under varying methods and assumptions resulted in well-resolved, strongly supported trees concordant with many traditional ideas about higher dipteran phylogeny and with prior inferences from 28S rDNA. The most robustly supported major eremoneuran clades were Cyclorrhapha, Platypezoidea, Eumuscomorpha, Empidoidea, Atelestidae, Empidoidea exclusive of Atelestidae, Hybotidae s.l., Microphoridae+Dolichopodidae, and Empididae s. str. Because CAD is ubiquitous, apparently single copy (at least within holometabolous insects), readily obtained from several insect orders using primers described herein, and exhibits considerable phylogenetic utility, it should have wide applicability in insect molecular systematics.     

The tortoise and the hare II: relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis

Chloroplast DNA sequences are a primary source of data for plant molecular systematic studies. A few key papers have provided the molecular systematics community with universal primer pairs for noncoding regions that have dominated the field, namely trnL-trnF and trnK/matK. These two regions have provided adequate information to resolve species relationships in some taxa, but often provide little resolution at low taxonomic levels. To obtain better phylogenetic resolution, sequence data from these regions are often coupled with other sequence data. Choosing an appropriate cpDNA region for phylogenetic investigation is difficult because of the scarcity of information about the tempo of evolutionary rates among different noncoding cpDNA regions. The focus of this investigation was to determine whether there is any predictable rate heterogeneity among 21 noncoding cpDNA regions identified as phylogenetically useful at low levels. To test for rate heterogeneity among the different cpDNA regions, we used three species from each of 10 groups representing eight major phylogenetic lineages of phanerogams. The results of this study clearly show that a survey using as few as three representative taxa can be predictive of the amount of phylogenetic information offered by a cpDNA region and that rate heterogeneity exists among noncoding cpDNA regions.     

Hierarchical analysis of variation in the mitochondrial 16S rRNA gene among Hymenoptera

Nucleotide sequences from a 434-bp region of the 16S rRNA gene were analyzed for 65 taxa of Hymenoptera (ants, bees, wasps, parasitoid wasps, sawflies) to examine the patterns of variation within the gene fragment and the taxonomic levels for which it shows maximum utility in phylogeny estimation. A hierarchical approach was adopted in the study through comparison of levels of sequence variation among taxa at different taxonomic levels. As previously reported for many holometabolous insects, the 16S data reported here for Hymenoptera are highly AT-rich and exhibit strong site-to-site variation in substitution rate. More precise estimates of the shape parameter (alpha) of the gamma distribution and the proportion of invariant sites were obtained in this study by employing a reference phylogeny and utilizing maximum-likelihood estimation. The effectiveness of this approach to recovering expected phylogenies of selected hymenopteran taxa has been tested against the use of maximum parsimony. This study finds that the 16S gene is most informative for phylogenetic analysis at two different levels: among closely related species or populations, and among tribes, subfamilies, and families. Maximization of the phylogenetic signal extracted from the 16S gene at higher taxonomic levels may require consideration of the base composition bias and the site-to-site rate variation in a maximum-likelihood framework.      

Molecular phylogenetics of the family Cyprinidae (Actinopterygii: Cypriniformes) as evidenced by sequence variation in the first intron of S7 ribosomal protein-coding gene: further evidence from a nuclear gene of the systematic chaos in the family

The family Cyprinidae is the largest freshwater fish group in the world, including over 200 genera and 2100 species. The phylogenetic relationships of major clades within this family are simply poorly understood, largely because of the overwhelming diversity of the group; however, several investigators have advanced different hypotheses of relationships that pre- and post-date the use of shared-derived characters as advocated through phylogenetic systematics. As expected, most previous investigations used morphological characters. Recently, mitochondrial DNA (mtDNA) sequences and combined morphological and mtDNA investigations have been used to explore and advance our understanding of species relationships and test monophyletic groupings. Limitations of these studies include limited taxon sampling and a strict reliance upon maternally inherited mtDNA variation. The present study is the first endeavor to recover the phylogenetic relationships of the 12 previously recognized monophyletic subfamilies within the Cyprinidae using newly sequenced nuclear DNA (nDNA) for over 50 species representing members of the different previously hypothesized subfamily and family groupings within the Cyprinidae and from other cypriniform families as outgroup taxa. Hypothesized phylogenetic relationships are constructed using maximum parsimony and Basyesian analyses of 1042 sites, of which 971 sites were variable and 790 were phylogenetically informative. Using other appropriate cypriniform taxa of the families Catostomidae (Myxocyprinus asiaticus), Gyrinocheilidae (Gyrinocheilus aymonieri), and Balitoridae (Nemacheilus sp. and Beaufortia kweichowensis) as outgroups, the Cyprinidae is resolved as a monophyletic group. Within the family the genera Raiamas, Barilius, Danio, and Rasbora, representing many of the tropical cyprinids, represent basal members of the family. All other species can be classified into variably supported and resolved monophyletic lineages, depending upon analysis, that are consistent with or correspond to Barbini and Leuciscini. The Barbini includes taxa traditionally aligned with the subfamily Cyprininae sensu previous morphological revisionary studies by Howes (Barbinae, Labeoninae, Cyprininae and Schizothoracinae). The Leuciscini includes six other subfamilies that are mainly divided into three separate lineages. The relationships among genera and subfamilies are discussed as well as the possible origins of major lineages.     

A comparison of 18s ribosomal RNA and rubisco large subunit sequences for studying angiosperm phylogeny

Summary Partial sequences of 18s rRNA were obtained for 2 gymnosperms and 12 angiosperms from a wide range of families and these were analyzed with 5 other published sequences to form a phylogenetic tree. Using 16 published sequences of the large subunit of rubisco (rbcL), also from a wide range of angiosperm families, another phylogenetic tree was derived and the two approaches were compared. Both phylogenetic trees gave good grouping within families but in neither case was there resolution of the branching order of major taxa. Superficially the long rbcL sequences (whose base composition was homogeneous among all species) seemed very promising, but analysis showed that a large proportion of the variation did not affect the amino acid sequence. Although silent substitution contained some phylogenetic information, at the level required to order major taxa, much of it was random and obfuscating. It was concluded that neither macromolecule alone was likely to yield a solution to the problem of angiosperm phylogeny and therefore that studies of both, at least, will be required. For this reason, a method wa described for obtaining both DNA and RNA of good quality from the same preparation and which had been used successfully with a wide range of species including many with pungent leaves.     

Maximum-likelihood models for combined analyses of multiple sequence data

Models of nucleotide substitution were constructed for combined analyses of heterogeneous sequence data (such as those of multiple genes) from the same set of species. The models account for different aspects of the heterogeneity in the evolutionary process of different genes, such as differences in nucleotide frequencies, in substitution rate bias (for example, the transition/transversion rate bias), and in the extent of rate variation across sites. Model parameters were estimated by maximum likelihood and the likelihood ratio test was used to test hypotheses concerning sequence evolution, such as rate constancy among lineages (the assumption of a molecular clock) and proportionality of branch lengths for different genes. The example data from a segment of the mitochondrial genome of six hominoid species (human, common and pygmy chimpanzees, gorilla, orangutan, and siamang) were analyzed. Nucleotides at the three codon positions in the protein-coding regions and from the tRNA-coding regions were considered heterogeneous data sets. Statistical tests showed that the amount of evolution in the sequence data reflected in the estimated branch lengths can be explained by the codon-position effect and lineage effect of substitution rates. The assumption of a molecular clock could not be rejected when the data were analyzed separately or when the rate variation among sites was ignored. However, significant differences in substitution rate among lineages were found when the data sets were combined and when the rate variation among sites was accounted for in the models. Under the assumption that the orangutan and African apes diverged 13 million years ago, the combined analysis of the sequence data estimated the times for the human-chimpanzee separation and for the separation of the gorilla as 4.3 and 6.8 million years ago, respectively.     

Phylogenetic relationships among clonal groups of extraintestinal pathogenic Escherichia coli as assessed by multi-locus sequence analysis

The evolutionary origins of extraintestinal pathogenic Escherichia coli (ExPEC) remain uncertain despite these organisms' relevance to human disease. A valid understanding of ExPEC phylogeny is needed as a framework against which the observed distribution of virulence factors and clinical associations can be analyzed. Accordingly, phylogenetic relationships were defined by multi-locus sequence analysis among 44 representatives of selected ExPEC clonal groups and the E. coli Reference (ECOR) collection. Recombination, which significantly obscured the phylogenetic signal for several strains, was dealt with by excluding strains or specific sequences. Conflicting overall phylogenies, and internal phylogenies for virulence-associated phylogenetic group B2, were inferred depending on the specific dataset (i.e., how extensively purged of recombination), outgroup (Salmonella enterica and/or Escherichia fergusonii), and analysis method (neighbor joining, maximum parsimony, maximum likelihood, or Bayesian likelihood). Nonetheless, the major E. coli phylogenetic groups A, B1, and B2 were consistently well resolved, as was a major sub-component of group D and an ECOR 37-O157:H7 clade. Moreover, nine important ExPEC clonal groups within groups B2 and D, characterized by serotypes O6:K2:H1, O18:K1:H7, O6:H31, and O4:K+:H+ (from group B2), and O1:K1:H-, O7:K1:H-, O157:K+:H (non-7), O15:K52:H1, and O11/17/77:K52:H18 ("clonal group A") (from group D), were consistently well resolved, regardless of clinical background (cystitis, pyelonephritis, neonatal meningitis, sepsis, or fecal), host group, geographical origin, and virulence profile. Among the group B2-derived clonal groups the O6:K2:H1 clade appeared basal. Within group D, "clonal group A" and the O15:K52:H1 clonal group were consistently placed with ECOR 47 and ECOR 44, respectively, as nearest neighbors. These findings clarify phylogenetic relationships among key ExPEC clonal groups but also emphasize that recombination appears to obscure the oldest evolutionary relationships, despite extensive targeted sequencing and use of a wide range of analysis techniques.     

Utility of the D1 domain of nuclear 28S rRNA for phylogenetic inference in the Digenea

The relatively variable D1 domain near the 5 end of the 28S ribosomal RNA gene (large subunit rRNA) has been sequenced for 12 species of digenean trematodes. Phylogenetic relationships among these species and Schistosoma mansoni (previously published sequence) were investigated with maximum parsimony and distance methods. A DNA sequence from the tapeworm (Cestoda) Hymenolepis diminuta was used for outgroup comparison. In all analyses, several clusters of taxa appeared repeatedly: (1) four fasciolids; (2) three didymozoids with a hemiurid; and (3) two lepocreadiids with a gyliauchenid. Bootstrap resampling of the data revealed that the first two clusters were well supported. In contrast there was less support for the cluster containing the lepocreadiids and gyliauchenid. We conclude that the D1 domain is too variable for phylogenetic inference among distantly related families, but it is well suited to phylogenetic inference within and among closely related families in the Digenea.