Flying lemurs – The 'flying tree shrews'? Molecular cytogenetic evidence for a Scandentia-Dermoptera sister clade

Background  

Flying lemurs or Colugos (order Dermoptera) represent an ancient mammalian lineage that contains only two extant species. Although molecular evidence strongly supports that the orders Dermoptera, Scandentia, Lagomorpha, Rodentia and Primates form a superordinal clade called Supraprimates (or Euarchontoglires), the phylogenetic placement of Dermoptera within Supraprimates remains ambiguous.     

Evaluating the Phylogenetic Position of Chinese Tree Shrew （Tupaia belangeri chinensis） Based on Complete Mitochondrial Genome： Implication for Using Tree Shrew as an Alternative Experimental Animal to Primates in Biomedical Research

Tree shrew(Tupaia belangeri) is currently placed in Order Scandentia and has a wide distribution in Southeast Asia and Southwest China.Due to its unique characteristics,such as small body size,high brain-to-body mass ratio,short reproductive cycle and life span,and low-cost of maintenance,tree shrew has been proposed to be an alternative experimental animal to primates in biomedical research.However,there are some debates regarding the exact phylogenetic affinity of tree shrew to primates.In this study,we determined the mtDNA entire genomes of three Chinese tree shrews(T.belangeri chinensis) and one Malayan flying lemur(Galeopterus variegatus).Combined with the published data for species in Euarchonta,we intended to discern the phylogenetic relationship among representative species of Dermoptera,Scandentia and Primates.The mtDNA genomes of Chinese tree shrews and Malayan flying lemur shared similar gene organization and structure with those of other mammals.Phylogenetic analysis based on 12 concatenated mitochondrial proteinencoding genes revealed a closer relationship between species of Scandentia and Glires,whereas species of Dermoptera were clustered with Primates.This pattern was consistent with previously reported phylogeny based on mtDNA data,but differed from the one reconstructed on the basis of nuclear genes.Our result suggested that the matrilineal affinity of tree shrew to primates may not be as close as we had thought.The ongoing project for sequencing the entire genome of Chinese tree shrew will provide more information to clarify this important issue.     

Ultraconserved elements show utility in phylogenetic inference of Adephaga (Coleoptera) and suggest paraphyly of ‘Hydradephaga’

The beetle suborder Adephaga has been the subject of many phylogenetic reconstructions utilizing a variety of data sources and inference methods. However, no strong consensus has yet emerged on the relationships among major adephagan lineages. Ultraconserved elements (UCEs) have proved useful for inferring difficult or unresolved phylogenies at varying timescales in vertebrates, arachnids and Hymenoptera. Recently, a UCE bait set was developed for Coleoptera using polyphagan genomes and a member of the order Strepsiptera as an outgroup. Here, we examine the utility of UCEs for reconstructing the phylogeny of adephagan families, in the first in vitro application a UCE bait set in Coleoptera. Our final dataset included 305 UCE loci for 18 representatives of all adephagan families except Aspidytidae, and two polyphagan outgroups, with a total concatenated length of 83 547 bp. We inferred trees using maximum likelihood analyses of the concatenated UCE alignment and coalescent species tree methods (astral ii , ASTRID, svdquartets ). Although the coalescent species tree methods had poor resolution and weak support, concatenated analyses produced well‐resolved, highly supported trees. Hydradephaga was recovered as paraphyletic, with Gyrinidae sister to Geadephaga and all other adephagans. Haliplidae was recovered as sister to Dytiscoidea, with Hygrobiidae and Amphizoidae successive sisters to Dytiscidae. Finally, Noteridae was recovered as monophyletic and sister to Meruidae. Given the success of UCE data for resolving phylogenetic relationships within Adephaga, we suggest the potential for further resolution of relationships within Adephaga using UCEs with improved taxon sampling, and by developing Adephaga‐specific probes.     

Coalescent methods for estimating species trees from phylogenomic data

Genome-scale sequence data have become increasingly available in the phylogenetic studies for understanding the evolutionary histories of species. However, it is challenging to develop probabilistic models to account for heterogeneity of phylogenomic data. The multispecies coalescent model describes gene trees as independent random variables generated from a coalescence process occurring along the lineages of the species tree. Since the multispecies coalescent model allows gene trees to vary across genes, coalescent-based methods have been popularly used to account for heterogeneous gene trees in phylogenomic data analysis. In this paper, we summarize and evaluate the performance of coalescent-based methods for estimating species trees from genome-scale sequence data. We investigate the effects of deep coalescence and mutation on the performance of species tree estimation methods. We found that the coalescent-based methods perform well in estimating species trees for a large number of genes, regardless of the degree of deep coalescence and mutation. The performance of the coalescent methods is negatively correlated with the lengths of internal branches of the species tree.     

Reconciling gene trees of eastern North American minnows

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

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.     

Phylogenetic analyses of complete mitochondrial genome sequences suggest a basal divergence of the enigmatic rodent <Emphasis Type="Italic">Anomalurus</Emphasis>

Background  

Phylogenetic relationships between Lagomorpha, Rodentia and Primates and their allies (Euarchontoglires) have long been debated. While it is now generally agreed that Rodentia constitutes a monophyletic sister-group of Lagomorpha and that this clade (Glires) is sister to Primates and Dermoptera, higher-level relationships within Rodentia remain contentious.     

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

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

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

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

The sequence of the largest subunit of RNA polymerase II is a useful marker for inferring seed plant phylogeny

We used RT-PCR to sequence approximately 3 kb of the gene coding for the largest subunit of RNA polymerase II (rpb1) from nine land plants. Our results show that plant rpb1 genes all have a similar GC-content and that their amino acid sequences evolve at a similar rate in most species we examined, except for the Arabidopsis thaliana and rice sequences which evolve faster. This gene also exists as a single copy in most species and contains enough phylogenetically informative sites to resolve the evolutionary relationships among seed plants. Protein maximum parsimony, as well as neighbor-joining and maximum likelihood analyses of DNA and protein sequences, all generated identical tree topologies with similar strong support values at each node. The angiosperms are a clade comprising Amborella as a sister group to all other angiosperms, followed by Nymphaea, Magnolia, Arabidopsis, and a monocot clade containing maize and rice. The gymnosperms also form a monophyletic clade with Welwitschia and pine grouped together and sister to a Cycas and Zamia clade. These findings concur with recent studies that refute the Anthophyte Hypothesis and place Amborella at the base of the angiosperm tree. These rpb1 sequences also give a more consistent picture of seed plant relationships than similar analyses performed on data sets made of 18S rDNA, atpB, and rbcL sequences from the same species. These sequences therefore show great promise to help further resolve the phylogenetic relationships of seed plants.     

Phylogenetic Analysis of Eukaryotes Using Heat-Shock Protein Hsp90

Most eukaryote molecular phylogenies have been based on small-subunit ribosomal RNA as its database includes the most species, but serious problems have been encountered that can make these trees misleading. More recent studies using concatenated protein sequences have increased the data per organism, reducing misleading signals from a single sequence, but taxon sampling is limited. To increase the database of protein-coding genes we sequenced the cytosolic form of heat-shock protein Hsp90 from a broad variety of previously unsampled eukaryote groups: protozoan flagellates (phyla Choanozoa, Apusozoa, Cercozoa) and all three groups of chromists (Cryptophyta, Heterokonta, Haptophyta). Gamma-corrected distance trees robustly show three groups: bacterial sequences are sister to all eukaryote sequences, which are cleanly subdivided into the cytosolic sequences and a clade comprising the chloroplast and endoplasmic reticulum (ER) Hsp90 sequences. The eukaryote cytosolic sequences comprise a robust opisthokont clade (animals/Choanozoa/fungi), a bikont clade, and an amoebozoan branch. However their topology is not robust. When the cytosolic sequences are rooted using only the ER/chloroplast clade as outgroup the amoebozoan Dictyostelium is sister to the opisthokonts forming a unikont clade in the distance tree. Congruence of this tree with that for concatenated mitochondrial proteins suggests that the root of the eukaryote tree is between unikonts and bikonts. Gamma-corrected maximum likelihood analyses of cytosolic sequences alone (519 unambiguously aligned amino acid positions) show bikonts as a clade, as do least-squares distance trees, but with other distance methods and parsimony the sole amoebozoan species branches weakly within bikonts. Choanozoa are clearly sisters to animals. Some major bikont groups (e.g. green plants, alveolates, Euglenozoa) are consistently recovered, but others (e.g. discicristates, chromalveolates) appear only in some trees; the backbone of the bikont subtree is not resolved, the position of groups represented only by single sequences being particularly unclear. Although single-gene trees will probably never resolve these uncertainties, the congruence of Hsp90 trees with other data is greater than for most other molecules and further taxon sampling of this molecule is recommended.     

Reconstructing the Fungal Tree of Life Using Phylogenomics and a Preliminary Investigation of the Distribution of Yeast Prion-Like Proteins in the Fungal Kingdom

We have used three independent phylogenomic approaches (concatenated alignments, single-, and multi-gene supertrees) to reconstruct the fungal tree of life (FTOL) using publicly available fungal genomes. This is the first time multi-gene families have been used in fungal supertree reconstruction and permits us to use up to 66% of the 1,001,217 genes in our fungal database. Our analyses show that different phylogenomic datasets derived from varying clustering criteria and alignment orientation do not have a major effect on phylogenomic supertree reconstruction. Overall the resultant phylogenomic trees are relatively congruent with one another and successfully recover the major fungal phyla, subphyla and classes. We find that where incongruences do occur, the inferences are usually poorly supported. Within the Ascomycota phylum, our phylogenies reconstruct monophyletic Saccharomycotina and Pezizomycotina subphyla clades and infer a sister group relationship between these to the exclusion of the Taphrinomycotina. Within the Pezizomycotina subphylum, all three phylogenies infer a sister group relationship between the Leotiomycetes and Sordariomycetes classes. However, there is conflict regarding the relationships with the Dothideomycetes and Eurotiomycetes classes. Within the Basidiomycota phylum, supertrees derived from single- and multi-gene families infer a sister group relationship between the Pucciniomycotina and Agaricomycotina subphyla while the concatenated phylogeny infers a poorly supported relationship between the Agaricomycotina and Ustilagomycotina. The reconstruction of a robust FTOL is important for future fungal comparative analyses. We illustrate this point by performing a preliminary investigation into the phyletic distribution of yeast prion-like proteins in the fungal kingdom.     

How well do multispecies coalescent methods perform with mitochondrial genomic data? A case study of butterflies and moths (Insecta: Lepidoptera)

Despite the broad adoption of multispecies coalescent (MSC) methods for nuclear phylogenomics, they have yet to be applied to mitochondrial (mt) genomic data. As the potential sources of phylogenomic bias that MSC methods can address, such as incomplete lineage sorting, horizontal gene transfer and gene tree heterogeneity, have been found in mt genomic data, these approaches may improve the accuracy of phylogenetic inference with these data. In the present study, we examined the behaviour of MSC methods in reconstructing the phylogeny of Lepidoptera (butterflies and moths), a group for which mt genomic data are known to have strong resolving power. Traditional concatenation methods of analysing mt genomes for Lepidoptera infer topologies highly congruent with those generated from independent nuclear datasets. Individual mt gene trees performed poorly in recovering consensus relationships at deep levels (i.e. superfamily monophyly and inter-relationships) and only moderately well for shallow relationships (i.e. within Papilionoidea). In contrast, MSC analyses with ASTRAL performed strongly with almost complete concordance to both concatenated mt genome analyses and independent nuclear analyses at both deep and shallow phylogenetic scales. Outgroup choice had a limited impact on tree accuracy, with even phylogenetically distant outgroups still resulting in topologies highly congruent with results from nuclear datasets, although MSC analyses appeared to be marginally more affected by outgroup choice than concatenation analyses. In general, discordance between concatenation and MSC analyses was found at nodes whose resolution varied between previous nuclear phylogenomic studies. The sensitivity of individual relationships to analysis with MSC vs concatenation can thus be used to test the robustness of phylogenetic hypotheses. For insect phylogenetics, MSC is a reliable inference method for mt genomic data and is thus a useful complement to the already widely used concatenation approaches.     

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

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

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

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

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.     

Inferring Species Trees from Gene Trees: A Phylogenetic Analysis of the Elapidae (Serpentes) Based on the Amino Acid Sequences of Venom Proteins

Toward the goal of recovering the phylogenetic relationships among elapid snakes, we separately found the shortest trees from the amino acid sequences for the venom proteins phospholipase A₂and the short neurotoxin, collectively representing 32 species in 16 genera. We then applied a method we term gene tree parsimony for inferring species trees from gene trees that works by finding the species tree which minimizes the number of deep coalescences or gene duplications plus unsampled sequences necessary to fit each gene tree to the species tree. This procedure, which is both logical and generally applicable, avoids many of the problems of previous approaches for inferring species trees from gene trees. The results support a division of the elapids examined into sister groups of the Australian and marine (laticaudines and hydrophiines) species, and the African and Asian species. Within the former clade, the sea snakes are shown to be diphyletic, with the laticaudines and hydrophiines having separate origins. This finding is corroborated by previous studies, which provide support for the usefulness of gene tree parsimony.     

Resolving the phylogeny of lizards and snakes (Squamata) with extensive sampling of genes and species

Squamate reptiles (lizards and snakes) are one of the most diverse groups of terrestrial vertebrates. Recent molecular analyses have suggested a very different squamate phylogeny relative to morphological hypotheses, but many aspects remain uncertain from molecular data. Here, we analyse higher-level squamate phylogeny with a molecular dataset of unprecedented size, including 161 squamate species for up to 44 nuclear genes each (33 717 base pairs), using both concatenated and species-tree methods for the first time. Our results strongly resolve most squamate relationships and reveal some surprising results. In contrast to most other recent studies, we find that dibamids and gekkotans are together the sister group to all other squamates. Remarkably, we find that the distinctive scolecophidians (blind snakes) are paraphyletic with respect to other snakes, suggesting that snakes were primitively burrowers and subsequently re-invaded surface habitats. Finally, we find that some clades remain poorly supported, despite our extensive data. Our analyses show that weakly supported clades are associated with relatively short branches for which individual genes often show conflicting relationships. These latter results have important implications for all studies that attempt to resolve phylogenies with large-scale phylogenomic datasets.     

Eight independent nuclear genes support monophyly of the plovers: The role of mutational variance in gene trees

Molecular phylogenies of Charadriiformes based on mtDNA genes and one to three nuclear loci do not support the traditional placement of Pluvialis in the plovers (Charadriidae), assigning it instead to oystercatchers, stilts, and avocets (Haematopodidae and Recurvirostridae). To investigate this hypothesis of plover paraphyly, the relationships among Pluvialis and closely related families were revisited by sequencing two individuals of all taxa except Peltohyas for eight independent single copy nuclear protein-coding loci selected for their informativeness at this phylogenetic depth. The species tree estimated jointly with the gene trees in the coalescent programme (*)BEAST strongly supported plover monophyly, as did Bayesian analysis of the concatenated matrix. The data sets that supported plover paraphyly in Baker et al. (2007) and Fain and Houde (2007) reflect two to four independent gene histories, and thus discordance with the plover monophyly species tree might have arisen by chance through stochastic mutational variance. For the plovers we conclude there is no conclusive evidence of coalescent variance from ancient incomplete lineage sorting across the interior branch leading to Pluvialis in the species tree. Rather, earlier studies seem have been misled by faster evolving mtDNA genes with high mutational variance, and a few nuclear genes that had low resolving power at the Pluvialis sister group level. These findings are of general relevance in avian phylogenetics, as they show that careful attention needs to be paid to the number and the phylogenetic informativeness of genes required to obtain accurate estimates of the species tree, especially where there is mutational heterogeneity in gene trees.