Phylogeny of shorebirds, gulls, and alcids (Aves: Charadrii) from the cytochrome-b gene: parsimony, Bayesian inference, minimum evolution, and quartet puzzling

Charadrii (shorebirds, gulls, and alcids) have exceptional diversity in ecological, behavioral, and life-history traits. A phylogenetic framework is necessary to fully understand the relationships among these traits. Despite several attempts to resolve the phylogeny of the Charadrii, none have comprehensively utilized molecular sequence data. Complete and partial cytochrome-b gene sequences for 86 Charadrii and five Falconides species (as outgroup taxa) were obtained from GenBank and aligned. We analyzed the resulting matrices using parsimony, Bayesian inference, minimum evolution, and quartet puzzling methods. Posterior probabilities, decay indices, and bootstrapping provide strong support for four major lineages consisting of gulls, alcids, plovers, and sandpipers, respectively. The broad structure of the trees differ significantly from all previous hypotheses of Charadrii phylogeny in placing the plovers at the base of the tree below the sandpipers in a pectinate sequence towards a large clade of gulls and alcids. The parsimony, Bayesian, and minimum evolution models provide strong evidence for this phylogenetic hypothesis. This is further corroborated by non-tree based measures of support and conflict (Lento plots). The quartet puzzling trees are poorly resolved and inconclusive.     

Phylogenetic relationships of the New World monkeys (Primates, platyrrhini) based on nuclear G6PD DNA sequences

In order to test hypotheses about the phylogenetic relationships among living genera of New World monkeys, 1.3 kb of DNA sequence information was collected for two introns of the glucose-6-phosphate dehydrogenase (G6PD) locus, encoded on the X chromosome, for 24 species of New World monkeys. These data were analyzed using a maximum parsimony algorithm. The strict consensus of the three most-parsimonious gene trees that result shows support for the following clades: a pitheciine clade including Callicebus within which Chiropotes and Cacajao are sister taxa, an Alouatta-atelin clade within which Brachyteles is the sister taxon of Lagothrix and which is sister to another clade containing the callitrichines, and a callitrichine/Aotus/Cebus/Saimiri clade. Within the callitrichines, Callimico is the sister taxon of Callithrix. Cebus and Saimiri form a clade. These results are broadly consistent with previously published DNA sequence analyses of platyrrhine phylogeny and provide additional support for groupings provisionally proposed in those earlier studies. Nevertheless, questions remain as to the relative phylogenetic placement of Leontopithecus and Saguinus, the branching order within the Aotus/Cebus/Saimiri/callitrichine clade, and the placement of the pitheciine clade relative to the atelines and the callitrichines.     

The complete nucleotide sequence of the mitochondrial DNA of the agnathan Lampetra fluviatilis: bearings on the phylogeny of cyclostomes

There are two competing theories about the interrelationships of craniates: the cyclostome theory assumes that lampreys and hagfishes are a clade, the cyclostomes, whose sister group is the jawed vertebrates (gnathostomes); the vertebrate theory assumes that lampreys and gnathostomes are a clade, the vertebrates, whose sister group is hagfishes. The vertebrate theory is best supported by a number of unique anatomical and physiological characters. Molecular sequence data from 18S and 28S rRNA genes rather support the cyclostome theory, but mtDNA sequence of Myxine glutinosa rather supports the vertebrate theory. Additional molecular data are thus needed to elucidate this three-taxon problem. We determined the complete nucleotide sequence of the mtDNA of the lamprey Lampetra fluviatilis. The mtDNA of L. fluviatilis possesses the same genomic organization as Petromyzon marinus, which validates this gene order as a synapomorphy of lampreys. The mtDNA sequence of L. fluviatilis was used in combination with relevant mtDNA sequences for an approach to the hagfish/lamprey relationships using the maximum-parsimony, neighbor-joining, and maximum-likelihood methods. Although trees compatible with our present knowledge of the phylogeny of craniates can be reconstructed by using the three methods, the data collected do not support the vertebrate or the cyclostome hypothesis. The present data set does not allow the resolution of this three-taxon problem, and new kinds of data, such as nuclear DNA sequences, need to be collected.     

Phylogenetic relationships among loliginid squids (Cephalopoda: Myopsida) based on analyses of multiple data sets

The phylogenetic relationships among the loliginid squids, a species-rich group of shallowwater muscular squids, have been investigated recently using several approaches, including allozyme electrophoresis and analyses of morphological and DNA sequence data, yet no consensus has been reached. This study examines the effects of combining multiple data sets (morphology, allozymes and DNA sequence data from two mitochondrial genes) on estimates of loliginid phylogeny. Various data combinations were analysed under three maximum parsimony weighting schemes: equal weights for all characters, successive approximations and implicit weights parsimony. When feasible, support for branches within trees was assessed with nonparametric bootstrapping and decay analysis. Some ingroup relationships were consistent across all analyses, but relationships among outgroup taxa and basal ingroup taxa varied. Combining data increased bootstrap support for several nodes. Methods that downweight highly variable characters (i.e. successive approximations and implicit weights parsimony) produced very similar trees which included two major clades: a clade consisting of all species sampled from American waters (except Sepioteuthis ), and a clade of several east Atlantic species ( Loligo forbesi Steenstrup, Loligo vulgaris Lamarck and Loligo reynaudi d'Orbigny) plus several Indo-West Pacific species in the genera Uroteuthis and Loliolus. The Sepioteuthis species occupied a basal position within Loliginidae, but Sepioteuthis itself was not always monophyletic. The position of a clade of a few Lolliguncula species and Loligo (Alloteuthis) also varied across analyses. A new loliginid classification is proposed based on these findings.     

Nuclear and mtDNA phylogenies of the Trimeresurus complex: implications for the gene versus species tree debate

Phylogenies based on mitochondrial DNA (mtDNA) may represent gene trees that may not be congruent with the equivalent species tree. One solution to this problem is to include additional, independent loci from the nuclear genome. Sequence data from the seventh intron of the beta-fibrinogen gene were generated for 25 specimens of vipers, including 8 nominal species of the Trimeresurus complex of Asian pit vipers. Phylogenetic trees were generated using maximum-parsimony and maximum-likelihood methods. The taxonomic level at which the intron provided significant phylogenetic information was examined and the trees were compared to those produced from previously obtained mtDNA cytochrome b sequences. A variety of different approaches (separate analyses, conditional data combination, and consensus) were used in an attempt to provide a sound organismal phylogeny based on both nuclear and mtDNA data sets. We discuss the implications for the gene tree-species tree debate and its particular relevance to medically important organisms.     

Phylogenetic relationships in Maloideae (Rosaceae): evidence from sequences of the internal transcribed spacers of nuclear ribosomal DNA and its congruence with morphology

We used sequences from both internal transcribed spacers (ITS) and a small portion of the 5.8S gene of nuclear ribosomal DNA (nrDNA) for phylogenetic reconstruction of 19 genera of Maloideae and four potential outgroups from the Rosaceae. Parsimony analyses indicate that Maloideae are not monophyletic; Vauquelinia, which is traditionally placed in Spiraeoideae, and two genera of the Maloideae, Eriobotrya and Rhaphiolepis, form a well-supported clade that is the sister to the remainder of the subfamily. Although our ITS phylogenetic hypothesis is highly resolved, there is considerable homoplasy, and support, as indicated by bootstrap values and decay indices, is relatively weak for all groups except four: Eriobotrya-Rhaphiolepis-Vauquelinia, Crataegus-Mespilus, Amelanchier-Peraphyllum-Malacomeles, and Cydonia-Pseudocydonia. Our DNA sequence data do not support a broad interpretation of Sorbus. Intergeneric hybridization, which is prevalent in Maloideae, occurs between genera that are far removed from one another on our most-parsimonious trees. We infer an overall phylogeny from separate analyses of ITS DNA sequences and recently published morphological and wood anatomical studies of Maloideae and from analyses after pooling these data sets. The four most strongly supported clades of the ITS phylogeny appear in the phylogeny based on pooled data.     

A molecular view of primate phylogeny and important systematic and evolutionary questions

Phylogenetic analysis of extensive nucleotide sequence data from primate beta-globin gene clusters elucidates the systematics and evolution of the order Primates and reveals that rates of accumulation of mutations vary by as much as a factor of seven among different primate lineages. The picture of primate phylogeny from DNA sequences clarifies many ambiguities of the morphological picture. In the molecular picture, dwarf and brown lemurs group together into superfamily Lemuroidea, Lemuroidea and Lorisoidea into suborder Strepsirhini, and Tarsius and Anthropoidea into suborder Haplorhini. The molecular picture also provides both significant evidence for a human-chimpanzee clade that narrowly excludes gorilla and overwhelming evidence for the gorilla-chimpanzee-human clade within Hominoidea. Rates of DNA sequence evolution appear to have been fastest in the early primates ancestral to Anthropoidea and next fastest on the lorisoid branch. Rates were slowest over the past 25 Myr of hominoid descent, suggesting that mechanisms lowering the mutation rate evolved in correlation with lengthened life spans.     

Comparison of Bayesian and maximum likelihood bootstrap measures of phylogenetic reliability

Owing to the exponential growth of genome databases, phylogenetic trees are now widely used to test a variety of evolutionary hypotheses. Nevertheless, computation time burden limits the application of methods such as maximum likelihood nonparametric bootstrap to assess reliability of evolutionary trees. As an alternative, the much faster Bayesian inference of phylogeny, which expresses branch support as posterior probabilities, has been introduced. However, marked discrepancies exist between nonparametric bootstrap proportions and Bayesian posterior probabilities, leading to difficulties in the interpretation of sometimes strongly conflicting results. As an attempt to reconcile these two indices of node reliability, we apply the nonparametric bootstrap resampling procedure to the Bayesian approach. The correlation between posterior probabilities, bootstrap maximum likelihood percentages, and bootstrapped posterior probabilities was studied for eight highly diverse empirical data sets and were also investigated using experimental simulation. Our results show that the relation between posterior probabilities and bootstrapped maximum likelihood percentages is highly variable but that very strong correlations always exist when Bayesian node support is estimated on bootstrapped character matrices. Moreover, simulations corroborate empirical observations in suggesting that, being more conservative, the bootstrap approach might be less prone to strongly supporting a false phylogenetic hypothesis. Thus, apparent conflicts in topology recovered by the Bayesian approach were reduced after bootstrapping. Both posterior probabilities and bootstrap supports are of great interest to phylogeny as potential upper and lower bounds of node reliability, but they are surely not interchangeable and cannot be directly compared.     

Cladistic relationships of extant and fossil hominoids

There is general agreement that the hominoid primates form a monophyletic group, that the extant great apes and humans form a second clade within that group with the gibbons as the sister group, and that the African apes and humans form a third clade. Although it has recently been proposed that humans and orang utans are sister taxa and also that the great apes form a clade to the exclusion of humans, our analysis, particularly of the molecular evidence, supports the existence of an African ape and human clade. The major problem in hominoid phylogeny at present is the relationships of the species within this clade: morphological data generally support the existence of an African ape clade which is the sister group to humans; some molecular data also support this conclusion, but most molecular evidence indicates the existence of a chimpanzee/human clade. We have cladistically re-analysed the DNA and protein sequence data for which apomorphic character states can be assessed. It is clear that there is a high degree of homoplasy whichever branching pattern is produced, with some characters supporting the existence of a chimpanzee/human clade and others supporting an African ape clade. When the cladistic analyses of morphological and molecular data are combined we believe that the most parsimonious interpretation of the data is that the African apes form a clade which is the sister taxon of the human (i.e., Australopithecus, Homo and Paranthropus) clade.This paper is not intended as a survey of all hominoid fossils but as a study of branching points in hominoid evolution and fossils are included which are relevant to this branching pattern. The analysis of fossil taxa in this study leads us to conclude that Proconsul is the sister taxon to the later Hominoidea. A number of middle Miocene forms such as Dryopithecus, Kenyapithecus, Heliopithecus and Afropithecus are shown to share derived characters with great apes and humans and provide evidence for the divergence of that clade from the gibbon lineage prior to 18 Ma. The position that Sivapithecus represents the sister group of the orang utan clade is supported here and shows that the orang utan lineage had diverged from the African ape and human lineage prior to 11·5 Ma. There is unfortunately no definitive fossil cvidence on branching sequences within the African ape and human clade, although a new specimen from Samburu, Kenya may be related to the gorilla.     

Genome-scale phylogenetics: inferring the plant tree of life from 18,896 gene trees

Phylogenetic analyses using genome-scale data sets must confront incongruence among gene trees, which in plants is exacerbated by frequent gene duplications and losses. Gene tree parsimony (GTP) is a phylogenetic optimization criterion in which a species tree that minimizes the number of gene duplications induced among a set of gene trees is selected. The run time performance of previous implementations has limited its use on large-scale data sets. We used new software that incorporates recent algorithmic advances to examine the performance of GTP on a plant data set consisting of 18,896 gene trees containing 510,922 protein sequences from 136 plant taxa (giving a combined alignment length of >2.9 million characters). The relationships inferred from the GTP analysis were largely consistent with previous large-scale studies of backbone plant phylogeny and resolved some controversial nodes. The placement of taxa that were present in few gene trees generally varied the most among GTP bootstrap replicates. Excluding these taxa either before or after the GTP analysis revealed high levels of phylogenetic support across plants. The analyses supported magnoliids sister to a eudicot + monocot clade and did not support the eurosid I and II clades. This study presents a nuclear genomic perspective on the broad-scale phylogenic relationships among plants, and it demonstrates that nuclear genes with a history of duplication and loss can be phylogenetically informative for resolving the plant tree of life.     

Free from mitochondrial DNA: Nuclear genes and the inference of species trees among closely related darter lineages (Teleostei: Percidae: Etheostomatinae)

Investigations into the phylogenetics of closely related animal species are dominated by the use of mitochondrial DNA (mtDNA) sequence data. However, the near-ubiquitous use of mtDNA to infer phylogeny among closely related animal lineages is tempered by an increasing number of studies that document high rates of transfer of mtDNA genomes among closely related species through hybridization, leading to substantial discordance between phylogenies inferred from mtDNA and nuclear gene sequences. In addition, the recent development of methods that simultaneously infer a species phylogeny and estimate divergence times, while accounting for incongruence among individual gene trees, has ushered in a new era in the investigation of phylogeny among closely related species. In this study we assess if DNA sequence data sampled from a modest number of nuclear genes can resolve relationships of a species-rich clade of North American freshwater teleost fishes, the darters. We articulate and expand on a recently introduced method to infer a time-calibrated multi-species coalescent phylogeny using the computer program ^*BEAST. Our analyses result in well-resolved and strongly supported time-calibrated darter species tree. Contrary to the expectation that mtDNA will provide greater phylogenetic resolution than nuclear gene data; the darter species tree inferred exclusively from nuclear genes exhibits a higher frequency of strongly supported nodes than the mtDNA time-calibrated gene tree.     

Branch-length prior influences Bayesian posterior probability of phylogeny

The Bayesian method for estimating species phylogenies from molecular sequence data provides an attractive alternative to maximum likelihood with nonparametric bootstrap due to the easy interpretation of posterior probabilities for trees and to availability of efficient computational algorithms. However, for many data sets it produces extremely high posterior probabilities, sometimes for apparently incorrect clades. Here we use both computer simulation and empirical data analysis to examine the effect of the prior model for internal branch lengths. We found that posterior probabilities for trees and clades are sensitive to the prior for internal branch lengths, and priors assuming long internal branches cause high posterior probabilities for trees. In particular, uniform priors with high upper bounds bias Bayesian clade probabilities in favor of extreme values. We discuss possible remedies to the problem, including empirical and full Bayesian methods and subjective procedures suggested in Bayesian hypothesis testing. Our results also suggest that the bootstrap proportion and Bayesian posterior probability are different measures of accuracy, and that the bootstrap proportion, if interpreted as the probability that the clade is true, can be either too liberal or too conservative.     

Phylogenetic relationships among seed plants: Persistent questions and the limits of molecular data

Trees inferred from DNA sequence data provide only limited insight into the phylogeny of seed plants because the living lineages (cycads, Ginkgo, conifers, gnetophytes, and angiosperms) represent fewer than half of the major lineages that have been detected in the fossil record. Nevertheless, phylogenetic trees of living seed plants inferred from sequence data can provide a test of relationships inferred in analyses that include fossils. So far, however, significant uncertainty persists because nucleotide data support several conflicting hypotheses. It is likely that improved sampling of gymnosperm diversity in nucleotide data sets will help alleviate some of the analytical issues encountered in the estimation of seed plant phylogeny, providing a more definitive test of morphological trees. Still, rigorous morphological analyses will be required to answer certain fundamental questions, such as the identity of the angiosperm sister group and the rooting of crown seed plants. Moreover, it will be important to identify approaches for incorporating insights from data that may be accurate but less likely than sequence data to generate results supported by high bootstrap values. How best to weigh evidence and distinguish among hypotheses when some types of data give high support values and others do not remains an important problem.     

Strange bayes indeed: uniform topological priors imply non-uniform clade priors

While Bayesian analysis has become common in phylogenetics, the effects of topological prior probabilities on tree inference have not been investigated. In Bayesian analyses, the prior probability of topologies is almost always considered equal for all possible trees, and clade support is calculated from the majority rule consensus of the approximated posterior distribution of topologies. These uniform priors on tree topologies imply non-uniform prior probabilities of clades, which are dependent on the number of taxa in a clade as well as the number of taxa in the analysis. As such, uniform topological priors do not model ignorance with respect to clades. Here, we demonstrate that Bayesian clade support, bootstrap support, and jackknife support from 17 empirical studies are significantly and positively correlated with non-uniform clade priors resulting from uniform topological priors. Further, we demonstrate that this effect disappears for bootstrap and jackknife when data sets are free from character conflict, but remains pronounced for Bayesian clade supports, regardless of tree shape. Finally, we propose the use of a Bayes factor to account for the fact that uniform topological priors do not model ignorance with respect to clade probability.     

Assessing horizontal transfer of nifHDK genes in eubacteria: nucleotide sequence of nifK from Frankia strain HFPCcI3

The structural genes for nitrogenase, nifK, nifD, and nifH, are crucial for nitrogen fixation. Previous phylogenetic analysis of the amino acid sequence of nifH suggested that this gene had been horizontally transferred from a proteobacterium to the gram-positive/cyanobacterial clade, although the confounding effects of paralogous comparisons made interpretation of the data difficult. An additional test of nif gene horizontal transfer using nifD was made, but the NifD phylogeny lacked resolution. Here nif gene phylogeny is addressed with a phylogenetic analysis of a third and longer nif gene, nifK. As part of the study, the nifK gene of the key taxon Frankia was sequenced. Parsimony and some distance analyses of the nifK amino acid sequences provide support for vertical descent of nifK, but other distance trees provide support for the lateral transfer of the gene. Bootstrap support was found for both hypotheses in all trees; the nifK data do not definitively favor one or the other hypothesis. A parsimony analysis of NifH provides support for horizontal transfer in accord with previous reports, although bootstrap analysis also shows some support for vertical descent of the orthologous nifH genes. A wider sampling of taxa and more sophisticated methods of phylogenetic inference are needed to understand the evolution of nif genes. The nif genes may also be powerful phylogenetic tools. If nifK evolved by vertical descent, it provides strong evidence that the cyanobacteria and proteobacteria are sister groups to the exclusion of the firmicutes, whereas 16S rRNA sequences are unable to resolve the relationships of these three major eubacterial lineages.      

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.     

New nuclear evidence for the oldest divergence among neognath birds: the phylogenetic utility of ZENK (i)

To date, there is little consensus concerning the phylogenetic relationships among neognath orders, which include all extant birds except ratites and tinamous. Different data sets, both molecular and morphologic, have yielded radically different and often unresolved ordinal topologies, especially within the neoaves clade. This lack of resolution and ongoing conflict indicates a need for additional phylogenetic characters to be applied to the question of higher-level avian phylogeny. In this study, sequences of a single-copy nuclear gene, ZENK, were used to reconstruct an ordinal-level phylogeny of neognath birds. Strong support was indicated for the oldest divergence within Neognathae; the chicken- and duck-like birds formed a clade that was sister to all other modern birds. In addition, many families of traditional taxonomic orders clustered together in the ZENK tree, indicating the gene's general phylogenetic reliability. However, within the neoaves clade, there was little support for relationships among orders, which is a result similar to all other recent molecular studies of higher-level avian phylogeny. This similarity among studies suggests the possibility of a rapid radiation of the major neoaves lineages. Despite the ongoing lack of neoaves resolution, ZENK's sequence divergence and base composition patterns indicate its general utility as a new phylogenetic marker for higher-level avian systematics.     

Measures of clade confidence do not correlate with accuracy of phylogenetic trees

Metrics of phylogenetic tree reliability, such as parametric bootstrap percentages or Bayesian posterior probabilities, represent internal measures of the topological reproducibility of a phylogenetic tree, while the recently introduced aLRT (approximate likelihood ratio test) assesses the likelihood that a branch exists on a maximum-likelihood tree. Although those values are often equated with phylogenetic tree accuracy, they do not necessarily estimate how well a reconstructed phylogeny represents cladistic relationships that actually exist in nature. The authors have therefore attempted to quantify how well bootstrap percentages, posterior probabilities, and aLRT measures reflect the probability that a deduced phylogenetic clade is present in a known phylogeny. The authors simulated the evolution of bacterial genes of varying lengths under biologically realistic conditions, and reconstructed those known phylogenies using both maximum likelihood and Bayesian methods. Then, they measured how frequently clades in the reconstructed trees exhibiting particular bootstrap percentages, aLRT values, or posterior probabilities were found in the true trees. The authors have observed that none of these values correlate with the probability that a given clade is present in the known phylogeny. The major conclusion is that none of the measures provide any information about the likelihood that an individual clade actually exists. It is also found that the mean of all clade support values on a tree closely reflects the average proportion of all clades that have been assigned correctly, and is thus a good representation of the overall accuracy of a phylogenetic tree.     

Viburnum phylogeny based on combined molecular data: implications for taxonomy and biogeography

We investigated Viburnum phylogeny using separate and combined analyses of DNA sequence data from two chloroplast and three nuclear loci. Separate analyses of nuclear and chloroplast data sets resulted in gene trees that were generally congruent with one another and with trees from two previous analyses. Our gene trees do differ in the position of section Pseudotinus, as well as in species relationships within sections Pseudotinus and Lentago. However, tests for incongruence indicate that differences between the nuclear and chloroplast data are not significant. Furthermore, gene trees from combined analyses were highly similar to those found in separate analyses, suggesting that these localized differences do not affect other parts of the tree. Our analyses provide convincing support for numerous relationships, although there is still uncertainty at the base of the tree. To facilitate future study, we propose informal names for 12 well-supported species groups, as well as for several higher-level clades. We also discuss the biogeographic implications of our phylogeny, focusing on repeated, although apparently temporally incongruent, patterns of disjunction between the Old and New Worlds.     

Molecular Phylogeny of Fulgoromorpha (Insecta, Hemiptera, Archaeorrhyncha). The Enigmatic Tettigometridae: Evolutionary Affiliations and Historical Biogeography

Previous morphologically based studies by several taxonomists disagree on whether the Tettigometridae represent a sister group to all other fulgoromorphans or whether they are a relatively derived family within the Fulgoromorpha. In this study, several parsimony-based analyses using data sets composed of nucleotide sequences of 18S rDNAs (genes encoding 18S rRNAs) support a monophyletic Fulgoromorpha. All analyses depict Tettigometridae as a relatively derived lineage in a monophyletic relationship with Tropiduchidae. Unscored and unweighted data sets position the tettigometrid + tropiduchid clade as sister to Flatidae. The tettigometrid + tropiduchid clade is supported by four synapomorphic sites, two deletions and two transversions. Constraining tettigometrids to a basal fulgoromorphan lineage significantly reduces parsimony, with several hundreds to thousand of trees being more parsimonious. An analysis employing the Barriel method for scoring potential phylogenetic information in regions containing deletions also supports a non-basal tettigometrid + tropiduchid clade. However, this method results in three equally most parsimonious trees and a sister relationship of the tettigometrid + tropiduchid clade to Flatidae becomes ambiguous.
The molecular-based results showing a derived Tettigometridae agree with previous morphological interpretations of Bourgoin. The current biogeographical distribution of tettigometrids and morphological features supporting the 18S rDNA-based phylogeny are discussed.