BIRDS IN A BUSH: FIVE GENES INDICATE EXPLOSIVE EVOLUTION OF AVIAN ORDERS

Phylogenomics Reveals Ancient Gene Tree Discordance in the Amphibian Tree of Life

Molecular phylogenies have yielded strong support for many parts of the amphibian Tree of Life, but poor support for the resolution of deeper nodes, including relationships among families and orders. To clarify these relationships, we provide a phylogenomic perspective on amphibian relationships by developing a taxon-specific Anchored Hybrid Enrichment protocol targeting hundreds of conserved exons which are effective across the class. After obtaining data from 220 loci for 286 species (representing 94% of the families and 44% of the genera), we estimate a phylogeny for extant amphibians and identify gene tree-species tree conflict across the deepest branches of the amphibian phylogeny. We perform locus-by-locus genealogical interrogation of alternative topological hypotheses for amphibian monophyly, focusing on interordinal relationships. We find that phylogenetic signal deep in the amphibian phylogeny varies greatly across loci in a manner that is consistent with incomplete lineage sorting in the ancestral lineage of extant amphibians. Our results overwhelmingly support amphibian monophyly and a sister relationship between frogs and salamanders, consistent with the Batrachia hypothesis. Species tree analyses converge on a small set of topological hypotheses for the relationships among extant amphibian families. These results clarify several contentious portions of the amphibian Tree of Life, which in conjunction with a set of vetted fossil calibrations, support a surprisingly younger timescale for crown and ordinal amphibian diversification than previously reported. More broadly, our study provides insight into the sources, magnitudes, and heterogeneity of support across loci in phylogenomic data sets.[AIC; Amphibia; Batrachia; Phylogeny; gene tree-species tree discordance; genomics; information theory.].

Deep-Time Demographic Inference Suggests Ecological Release as Driver of Neoavian Adaptive Radiation

Assessing the applicability of theory to major adaptive radiations in deep time represents an extremely difficult problem in evolutionary biology. Neoaves, which includes 95% of living birds, is believed to have undergone a period of rapid diversification roughly coincident with the Cretaceous–Paleogene (K-Pg) boundary. We investigate whether basal neoavian lineages experienced an ecological release in response to ecological opportunity, as evidenced by density compensation. We estimated effective population sizes (Ne) of basal neoavian lineages by combining coalescent branch lengths (CBLs) and the numbers of generations between successive divergences. We used a modified version of Accurate Species TRee Algorithm (ASTRAL) to estimate CBLs directly from insertion–deletion (indel) data, as well as from gene trees using DNA sequence and/or indel data. We found that some divergences near the K-Pg boundary involved unexpectedly high gene tree discordance relative to the estimated number of generations between speciation events. The simplest explanation for this result is an increase in Ne, despite the caveats discussed herein. It appears that at least some early neoavian lineages, similar to the ancestor of the clade comprising doves, mesites, and sandgrouse, experienced ecological release near the time of the K-Pg mass extinction.

Forcing external constraints on tree inference using ASTRAL

Background

To account for genome-wide discordance among gene trees, several widely-used methods seek to find a species tree with the minimum distance to input gene trees. To efficiently explore the large space of species trees, some of these methods, including ASTRAL, use dynamic programming (DP). The DP paradigm can restrict the search space, and thus, ASTRAL and similar methods use heuristic methods to define a restricted search space. However, arbitrary constraints provided by the user on the output tree cannot be trivially incorporated into such restrictions. The ability to infer trees that honor user-defined constraints is needed for many phylogenetic analyses, but no solution currently exists for constraining the output of ASTRAL.

Results

We introduce methods that enable the ASTRAL dynamic programming to infer constrained trees in an effective and scalable manner. To do so, we adopt a recently developed tree completion algorithm and extend it to allow multifurcating input and output trees. In simulation studies, we show that the approach for honoring constraints is both effective and fast. On real data, we show that constrained searches can help interrogate branches not recovered in the optimal ASTRAL tree to reveal support for alternative hypotheses.

Conclusions

The new algorithm is added ASTRAL to all user-provided constraints on the species tree.

Comparative Phylogenomics, a Stepping Stone for Bird Biodiversity Studies

Birds are a group with immense availability of genomic resources, and hundreds of forthcoming genomes at the doorstep. We review recent developments in whole genome sequencing, phylogenomics, and comparative genomics of birds. Short read based genome assemblies are common, largely due to efforts of the Bird 10K genome project (B10K). Chromosome-level assemblies are expected to increase due to improved long-read sequencing. The available genomic data has enabled the reconstruction of the bird tree of life with increasing confidence and resolution, but challenges remain in the early splits of Neoaves due to their explosive diversification after the Cretaceous-Paleogene (K-Pg) event. Continued genomic sampling of the bird tree of life will not just better reflect their evolutionary history but also shine new light onto the organization of phylogenetic signal and conflict across the genome. The comparatively simple architecture of avian genomes makes them a powerful system to study the molecular foundation of bird specific traits. Birds are on the verge of becoming an extremely resourceful system to study biodiversity from the nucleotide up.

Phylogenetic Signal of Indels and the Neoavian Radiation

The early radiation of Neoaves has been hypothesized to be an intractable “hard polytomy”. We explore the fundamental properties of insertion/deletion alleles (indels), an under-utilized form of genomic data with the potential to help solve this. We scored >5 million indels from >7000 pan-genomic intronic and ultraconserved element (UCE) loci in 48 representatives of all neoavian orders. We found that intronic and UCE indels exhibited less homoplasy than nucleotide (nt) data. Gene trees estimated using indel data were less resolved than those estimated using nt data. Nevertheless, Accurate Species TRee Algorithm (ASTRAL) species trees estimated using indels were generally similar to nt-based ASTRAL trees, albeit with lower support. However, the power of indel gene trees became clear when we combined them with nt gene trees, including a striking result for UCEs. The individual UCE indel and nt ASTRAL trees were incongruent with each other and with the intron ASTRAL trees; however, the combined indel+nt ASTRAL tree was much more congruent with the intronic trees. Finally, combining indel and nt data for both introns and UCEs provided sufficient power to reduce the scope of the polytomy that was previously proposed for several supraordinal lineages of Neoaves.

Filtering nucleotide sites by phylogenetic signal to noise ratio increases confidence in the Neoaves phylogeny generated from ultraconserved elements

Despite genome scale analyses, high-level relationships among Neoaves birds remain contentious. The placements of the Neoaves superorders are notoriously difficult to resolve because they involve deep splits followed by short internodes. Using our approach, we investigate whether filtering UCE loci on their phylogenetic signal to noise ratio helps to resolve key nodes in the Neoaves tree of life. We find that our analysis of data sets filtered for high signal to noise ratio results in topologies that are inconsistent with unfiltered results but that are congruent with whole-genome analyses. These relationships include the Columbea + Passerea sister relationship and the Phaethontimorphae + Aequornithia sister relationship. We also find increased statistical support for more recent nodes (i.e. the Pelecanidae + Ardeidae sister relationship, the Eucavitaves clade, and the Otidiformes + Musophagiformes sister relationship). We also find instances where support is reduced for well-established clades, possibly due to the removal of sites with moderate signal-to-noise ratio. Our results suggest that filtering on the basis of signal to noise ratio is a useful tool for resolving problematic splits in phylogenomic data sets.

Testing for Polytomies in Phylogenetic Species Trees Using Quartet Frequencies

Phylogenetic species trees typically represent the speciation history as a bifurcating tree. Speciation events that simultaneously create more than two descendants, thereby creating polytomies in the phylogeny, are possible. Moreover, the inability to resolve relationships is often shown as a (soft) polytomy. Both types of polytomies have been traditionally studied in the context of gene tree reconstruction from sequence data. However, polytomies in the species tree cannot be detected or ruled out without considering gene tree discordance. In this paper, we describe a statistical test based on properties of the multi-species coalescent model to test the null hypothesis that a branch in an estimated species tree should be replaced by a polytomy. On both simulated and biological datasets, we show that the null hypothesis is rejected for all but the shortest branches, and in most cases, it is retained for true polytomies. The test, available as part of the Accurate Species TRee ALgorithm (ASTRAL) package, can help systematists decide whether their datasets are sufficient to resolve specific relationships of interest.

Why Do Phylogenomic Data Sets Yield Conflicting Trees? Data Type Influences the Avian Tree of Life more than Taxon Sampling

Phylogenomics, the use of large-scale data matrices in phylogenetic analyses, has been viewed as the ultimate solution to the problem of resolving difficult nodes in the tree of life. However, it has become clear that analyses of these large genomic data sets can also result in conflicting estimates of phylogeny. Here, we use the early divergences in Neoaves, the largest clade of extant birds, as a "model system" to understand the basis for incongruence among phylogenomic trees. We were motivated by the observation that trees from two recent avian phylogenomic studies exhibit conflicts. Those studies used different strategies: 1) collecting many characters [$\sim$ 42 mega base pairs (Mbp) of sequence data] from 48 birds, sometimes including only one taxon for each major clade; and 2) collecting fewer characters ($\sim$ 0.4 Mbp) from 198 birds, selected to subdivide long branches. However, the studies also used different data types: the taxon-poor data matrix comprised 68% non-coding sequences whereas coding exons dominated the taxon-rich data matrix. This difference raises the question of whether the primary reason for incongruence is the number of sites, the number of taxa, or the data type. To test among these alternative hypotheses we assembled a novel, large-scale data matrix comprising 90% non-coding sequences from 235 bird species. Although increased taxon sampling appeared to have a positive impact on phylogenetic analyses the most important variable was data type. Indeed, by analyzing different subsets of the taxa in our data matrix we found that increased taxon sampling actually resulted in increased congruence with the tree from the previous taxon-poor study (which had a majority of non-coding data) instead of the taxon-rich study (which largely used coding data). We suggest that the observed differences in the estimates of topology for these studies reflect data-type effects due to violations of the models used in phylogenetic analyses, some of which may be difficult to detect. If incongruence among trees estimated using phylogenomic methods largely reflects problems with model fit developing more "biologically-realistic" models is likely to be critical for efforts to reconstruct the tree of life. [Birds; coding exons; GTR model; model fit; Neoaves; non-coding DNA; phylogenomics; taxon sampling.].

Tapping the woodpecker tree for evolutionary insight

Molecular phylogenetic studies of woodpeckers (Picidae) have generally focused on relationships within specific clades or have sampled sparsely across the family. We compared DNA sequences of six loci from 203 of the 217 recognized species of woodpeckers to construct a comprehensive tree of intrafamilial relationships. We recovered many known, but also numerous unknown, relationships among clades and species. We found, for example, that the three picine tribes are related as follows (Picini, (Campephilini, Melanerpini)) and that the genus Dinopium is paraphyletic. We used the tree to analyze rates of diversification and biogeographic patterns within the family. Diversification rate increased on two occasions during woodpecker history. We also tested diversification rates between temperate and tropical species but found no significant difference. Biogeographic analysis supported an Old World origin of the family and identified at least six independent cases of New World-Old World sister relationships. In light of the tree, we discuss how convergence, mimicry, and potential cases of hybridization have complicated woodpecker taxonomy.

Simultaneous speciation in the European high mountain flowering plant genus Facchinia (Minuartia s.l., Caryophyllaceae) revealed by genotyping-by-sequencing

Understanding the relative importance of different mechanisms of speciation in a given lineage requires fully resolved interspecific relationships. Using Facchinia, a genus of seven species centred in the European Alps, we explore whether the polytomy found by Sanger sequencing analyses of standard nuclear (ITS) and plastid markers (trnQ-rps16) is a hard or soft polytomy by substantially increasing the amount of DNA sequence data, generated by genotyping-by-sequencing. In comparison to 142 phylogenetically informative sites in the Sanger sequences the GBS sequences yielded 3363 phylogenetically informative sites after exclusion of apparently oversaturated SNPs. Maximum parsimony, maximum likelihood, NeighborNet, SVDquartets and Astral-II analyses all resulted in phylogenetic trees (and networks) in which interspecific relationships were largely unresolved. After excluding incomplete lineage sorting, hybridisation and oversaturation of characters as possible causes for lack of phylogenetic resolution, we conclude that the polytomy obtained most likely represents a hard polytomy. We hypothesize that diversification of Facchinia is best interpreted as the result of multiple simultaneous vicariance in response to climatic changes during the Early Quaternary.

Phylogenomics of a rapid radiation: is chromosomal evolution linked to increased diversification in north american spiny lizards (Genus Sceloporus)?

Background

Resolving the short phylogenetic branches that result from rapid evolutionary diversification often requires large numbers of loci. We collected targeted sequence capture data from 585 nuclear loci (541 ultraconserved elements and 44 protein-coding genes) to estimate the phylogenetic relationships among iguanian lizards in the North American genus Sceloporus. We tested for diversification rate shifts to determine if rapid radiation in the genus is correlated with chromosomal evolution.

Results

The phylogenomic trees that we obtained for Sceloporus using concatenation and coalescent-based species tree inference provide strong support for the monophyly and interrelationships among nearly all major groups. The diversification analysis supported one rate shift on the Sceloporus phylogeny approximately 20–25 million years ago that is associated with the doubling of the speciation rate from 0.06 species/million years (Ma) to 0.15 species/Ma. The posterior probability for this rate shift occurring on the branch leading to the Sceloporus species groups exhibiting increased chromosomal diversity is high (posterior probability = 0.997).

Conclusions

Despite high levels of gene tree discordance, we were able to estimate a phylogenomic tree for Sceloporus that solves some of the taxonomic problems caused by previous analyses of fewer loci. The taxonomic changes that we propose using this new phylogenomic tree help clarify the number and composition of the major species groups in the genus. Our study provides new evidence for a putative link between chromosomal evolution and the rapid divergence and radiation of Sceloporus across North America.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0628-x) contains supplementary material, which is available to authorized users.

Mitochondrial introgression via ancient hybridization, and systematics of the Australian endemic pygopodid gecko genus Delma

Of the more than 1500 species of geckos found across six continents, few remain as unfamiliar as the pygopodids - Family Pygopodidae (Gray, 1845). These gekkotans are limited to Australia (44 species) and New Guinea (2 species), but have diverged extensively into the most ecologically diverse limbless radiation save Serpentes. Current phylogenetic understanding of the family has relied almost exclusively on two works, which have produced and synthesized an immense amount of morphological, geographical, and molecular data. However, current interspecific relationships within the largest genus Delma Gray 1831 are based chiefly upon data from two mitochondrial loci (16s, ND2). Here, we reevaluate the interspecific relationships within the genus Delma using two mitochondrial and four nuclear loci (RAG1, MXRA5, MOS, DYNLL1), and identify points of strong conflict between nuclear and mitochondrial genomic data. We address mito-nuclear discordance, and remedy this conflict by recognizing several points of mitochondrial introgression as the result of ancient hybridization events. Owing to the legacy value and intraspecific informativeness, we suggest the continued use of ND2 as a phylogenetic marker. Results identify strong support for species groups, but relationships among these clades, and the placement of several enigmatic taxa remain uncertain. We suggest a more careful review of Delma australis and the 'northwest Australia' clade. Accurately assessing and addressing species richness and relationships within this endemic Australian Gekkotan genus is relevant for understanding patterns of squamate speciation across the region.

Genomic support for a moa-tinamou clade and adaptive morphological convergence in flightless ratites

One of the most startling discoveries in avian molecular phylogenetics is that the volant tinamous are embedded in the flightless ratites, but this topology remains controversial because recent morphological phylogenies place tinamous as the closest relative of a monophyletic ratite clade. Here, we integrate new phylogenomic sequences from 1,448 nuclear DNA loci totaling almost 1 million bp from the extinct little bush moa, Chilean tinamou, and emu with available sequences from ostrich, elegant crested tinamou, four neognaths, and the green anole. Phylogenetic analysis using standard homogeneous models and heterogeneous models robust to common topological artifacts recovered compelling support for ratite paraphyly with the little bush moa closest to tinamous within ratites. Ratite paraphyly was further corroborated by eight independent CR1 retroposon insertions. Analysis of morphological characters reinterpreted on a 27-gene paleognath topology indicates that many characters are convergent in the ratites, probably as the result of adaptation to a cursorial life style.

Tandem duplications in the C-terminal domain of the mesotocin receptor exclusively identified among East Eurasian thrushes

Mesotocin is a neurohypophyseal hormone found in some non-mammalian vertebrates, including birds, reptiles, and amphibians. In this study, we identified and characterized 18-amino acid duplications in the C-terminal domain of the mesotocin receptor (MTR), specifically found in Turdus thrushes (Aves: Passeriforms: Turdidae). These duplicated elements are located in the distal part of the C-terminal tails of MTR and consist of amino acids that are highly conserved among major vertebrates. Intraspecific polymorphisms in a variable number of tandem duplications are commonly found in East Eurasian Turdus, but not in any other genus of Turdidae. Moreover, the genus Turdus can be further classified into 2 groups according to the presence or absence of a 3-amino acid deletion just adjacent to the putative palmitoylation site in the cytoplasmic C-terminal tail. The phylogeny presented here strongly supports the conspecific group of 4 East Eurasian thrushes (Turdus pallidus, T. chrysolaus, T. obscurus, and T. celaenops). Our findings, therefore, provide a new synapomorphy that can be used for phylogenetic assumptions and shed a light on the history of diversification within Eurasian Turdus clades.

Identifying localized biases in large datasets: a case study using the avian tree of life

Large-scale multi-locus studies have become common in molecular phylogenetics, with new studies continually adding to previous datasets in an effort to fully resolve the tree of life. Total evidence analyses that combine existing data with newly collected data are expected to increase the power of phylogenetic analyses to resolve difficult relationships. However, they might be subject to localized biases, with one or a few loci having a strong and potentially misleading influence upon the results. To examine this possibility we combined a newly collected 31-locus dataset that includes representatives of all major avian lineages with a published dataset of 19 loci that has a comparable number of sites (Hackett et al., 2008. Science 320, 1763-1768). This allowed us to explore the advantages of conducting total evidence analyses, and to determine whether it was also important to analyze new datasets independent of published ones. The total evidence analysis yielded results very similar to the published results, with only slightly increased support at a few nodes. However, analyzing the 31- and 19-locus datasets separately highlighted several differences. Two clades received strong support in the published dataset and total evidence analysis, but the support appeared to reflect bias at a single locus (β-fibrinogen [FGB]). The signal in FGB that supported these relationships was sufficient to result in their recovery with bootstrap support, even when combined with 49 loci lacking that signal. FGB did not appear to have a substantial impact upon the results of species tree methods, but another locus (brain-derived neurotrophic factor [BDNF]) did have an impact upon those analyses. These results demonstrated that localized biases can influence large-scale phylogenetic analyses but they also indicated that considering independent evidence and exploring multiple analytical approaches could reveal them.

Parsimony and model-based analyses of indels in avian nuclear genes reveal congruent and incongruent phylogenetic signals

Insertion/deletion (indel) mutations, which are represented by gaps in multiple sequence alignments, have been used to examine phylogenetic hypotheses for some time. However, most analyses combine gap data with the nucleotide sequences in which they are embedded, probably because most phylogenetic datasets include few gap characters. Here, we report analyses of 12,030 gap characters from an alignment of avian nuclear genes using maximum parsimony (MP) and a simple maximum likelihood (ML) framework. Both trees were similar, and they exhibited almost all of the strongly supported relationships in the nucleotide tree, although neither gap tree supported many relationships that have proven difficult to recover in previous studies. Moreover, independent lines of evidence typically corroborated the nucleotide topology instead of the gap topology when they disagreed, although the number of conflicting nodes with high bootstrap support was limited. Filtering to remove short indels did not substantially reduce homoplasy or reduce conflict. Combined analyses of nucleotides and gaps resulted in the nucleotide topology, but with increased support, suggesting that gap data may prove most useful when analyzed in combination with nucleotide substitutions.

A phylogeny of birds based on over 1,500 loci collected by target enrichment and high-throughput sequencing

Evolutionary relationships among birds in Neoaves, the clade comprising the vast majority of avian diversity, have vexed systematists due to the ancient, rapid radiation of numerous lineages. We applied a new phylogenomic approach to resolve relationships in Neoaves using target enrichment (sequence capture) and high-throughput sequencing of ultraconserved elements (UCEs) in avian genomes. We collected sequence data from UCE loci for 32 members of Neoaves and one outgroup (chicken) and analyzed data sets that differed in their amount of missing data. An alignment of 1,541 loci that allowed missing data was 87% complete and resulted in a highly resolved phylogeny with broad agreement between the Bayesian and maximum-likelihood (ML) trees. Although results from the 100% complete matrix of 416 UCE loci were similar, the Bayesian and ML trees differed to a greater extent in this analysis, suggesting that increasing from 416 to 1,541 loci led to increased stability and resolution of the tree. Novel results of our study include surprisingly close relationships between phenotypically divergent bird families, such as tropicbirds (Phaethontidae) and the sunbittern (Eurypygidae) as well as between bustards (Otididae) and turacos (Musophagidae). This phylogeny bolsters support for monophyletic waterbird and landbird clades and also strongly supports controversial results from previous studies, including the sister relationship between passerines and parrots and the non-monophyly of raptorial birds in the hawk and falcon families. Although significant challenges remain to fully resolving some of the deep relationships in Neoaves, especially among lineages outside the waterbirds and landbirds, this study suggests that increased data will yield an increasingly resolved avian phylogeny.

An Update of Wallace’s Zoogeographic Regions of the World

Modern attempts to produce biogeographic maps focus on the distribution of species, and the maps are typically drawn without phylogenetic considerations. Here, we generate a global map of zoogeographic regions by combining data on the distributions and phylogenetic relationships of 21,037 species of amphibians, birds, and mammals. We identify 20 distinct zoogeographic regions, which are grouped into 11 larger realms. We document the lack of support for several regions previously defined based on distributional data and show that spatial turnover in the phylogenetic composition of vertebrate assemblages is higher in the Southern than in the Northern Hemisphere. We further show that the integration of phylogenetic information provides valuable insight on historical relationships among regions, permitting the identification of evolutionarily unique regions of the world.

A new genus of miniaturized and pug-nosed gecko from South America (Sphaerodactylidae: Gekkota)

Sphaerodactyl geckos comprise five genera distributed across Central and South America and the Caribbean. We estimated phylogenetic relationships among sphaerodactyl genera using both separate and combined analyses of seven nuclear genes. Relationships among genera were incongruent at different loci and phylogenies were characterized by short, in some cases zero length, internal branches and poor phylogenetic support at most nodes. We recovered a polyphyletic Coleodactylus, with Coleodactylus amazonicus being deeply divergent from the remaining Coleodactylus species sampled. The C. amazonicus lineage possessed unique codon deletions in the genes PTPN12 and RBMX while the remaining Coleodactylus species had unique codon deletions in RAG1. Topology tests could not reject a monophyletic Coleodactylus, but we show that short internal branch lengths decreased the accuracy of topology tests because there were not enough data along short branches to support one phylogenetic hypothesis over another. Morphological data corroborated results of the molecular phylogeny, with Coleodactylus exhibiting substantial morphological heterogeneity. We identified a suite of unique craniofacial features that differentiate C. amazonicus not only from other Coleodactylus species, but also from all other geckos. We describe this novel sphaerodactyl lineage as a new genus, Chatogekko gen. nov. We present a detailed osteology of Chatogekko, characterizing osteological correlates of miniaturization that provide a framework for future studies in sphaerodactyl systematics and biology.

A molecular phylogeny of Equatorial African Lacertidae, with the description of a new genus and species from eastern Democratic Republic of the Congo

Currently, four species of the lacertid lizard genus Adolfus are known from Central and East Africa. We sequenced up to 2,825 bp of two mitochondrial (16S and cyt b) and two nuclear (c-mos and RAG1) genes from 41 samples of Adolfus (representing every species), two species each of Gastropholis and Holaspis, and in separate analyses combined this data with GenBank sequences of all other Eremiadini genera and four Lacertini outgroups. Data from DNA sequences were analyzed with maximum parsimony (PAUP), maximum-likelihood (RAxML) and Bayesian inference (MrBayes) criteria. Results demonstrated that Adolfus is not monophyletic: A. africanus (type species), A. alleni and A. jacksoni are sister taxa, whereas A. vauereselli and a new species from the Itombwe Plateau of Democratic Republic of the Congo are in a separate lineage. Holaspis and Gastropholis were recovered in separate clades. Based on this molecular data, relatively substantial sequence divergence and multiple morphological differences, we describe a new genus of lacertid for the lineage including A. vauereselli and the new Itombwe species. The recognition of this new, endemic genus underscores the conservation importance of the Albertine Rift, especially the Itombwe Plateau, a unique region that is severely threatened by unchecked deforestation, mining and poaching.

Evolutionary relationships of the old world fruit bats (Chiroptera, Pteropodidae): another star phylogeny?

BACKGROUND

The family Pteropodidae comprises bats commonly known as megabats or Old World fruit bats. Molecular phylogenetic studies of pteropodids have provided considerable insight into intrafamilial relationships, but these studies have included only a fraction of the extant diversity (a maximum of 26 out of the 46 currently recognized genera) and have failed to resolve deep relationships among internal clades. Here we readdress the systematics of pteropodids by applying a strategy to try to resolve ancient relationships within Pteropodidae, while providing further insight into subgroup membership, by 1) increasing the taxonomic sample to 42 genera; 2) increasing the number of characters (to >8,000 bp) and nuclear genomic representation; 3) minimizing missing data; 4) controlling for sequence bias; and 5) using appropriate data partitioning and models of sequence evolution.

RESULTS

Our analyses recovered six principal clades and one additional independent lineage (consisting of a single genus) within Pteropodidae. Reciprocal monophyly of these groups was highly supported and generally congruent among the different methods and datasets used. Likewise, most relationships within these principal clades were well resolved and statistically supported. Relationships among the 7 principal groups, however, were poorly supported in all analyses. This result could not be explained by any detectable systematic bias in the data or incongruence among loci. The SOWH test confirmed that basal branches' lengths were not different from zero, which points to closely-spaced cladogenesis as the most likely explanation for the poor resolution of the deep pteropodid relationships. Simulations suggest that an increase in the amount of sequence data is likely to solve this problem.

CONCLUSIONS

The phylogenetic hypothesis generated here provides a robust framework for a revised cladistic classification of Pteropodidae into subfamilies and tribes and will greatly contribute to the understanding of character evolution and biogeography of pteropodids. The inability of our data to resolve the deepest relationships of the major pteropodid lineages suggests an explosive diversification soon after origin of the crown pteropodids. Several characteristics of pteropodids are consistent with this conclusion, including high species diversity, great morphological diversity, and presence of key innovations in relation to their sister group.

Mesozoic retroposons reveal parrots as the closest living relatives of passerine birds

The relationships of passerines (such as the well-studied zebra finch) with non-passerine birds is one of the great enigmas of avian phylogenetic research, because decades of extensive morphological and molecular studies yielded highly inconsistent results between and within data sets. Here we show the first application of the virtually homoplasy-free retroposon insertions to this controversy. Our study examined ~200,000 retroposon-containing loci from various avian genomes and retrieved 51 markers resolving early bird phylogeny. Among these, we obtained statistically significant evidence that parrots are the closest and falcons the second-closest relatives of passerines, together constituting the Psittacopasserae and the Eufalconimorphae, respectively. Our new and robust phylogenetic framework has substantial implications for the interpretation of various conclusions drawn from passerines as model organisms. This includes insights of relevance to human neuroscience, as vocal learning (that is, birdsong) probably evolved in the psittacopasseran ancestor, >30 million years earlier than previously assumed.

Homoplastic microinversions and the avian tree of life

BACKGROUND

Microinversions are cytologically undetectable inversions of DNA sequences that accumulate slowly in genomes. Like many other rare genomic changes (RGCs), microinversions are thought to be virtually homoplasy-free evolutionary characters, suggesting that they may be very useful for difficult phylogenetic problems such as the avian tree of life. However, few detailed surveys of these genomic rearrangements have been conducted, making it difficult to assess this hypothesis or understand the impact of microinversions upon genome evolution.

RESULTS

We surveyed non-coding sequence data from a recent avian phylogenetic study and found substantially more microinversions than expected based upon prior information about vertebrate inversion rates, although this is likely due to underestimation of these rates in previous studies. Most microinversions were lineage-specific or united well-accepted groups. However, some homoplastic microinversions were evident among the informative characters. Hemiplasy, which reflects differences between gene trees and the species tree, did not explain the observed homoplasy. Two specific loci were microinversion hotspots, with high numbers of inversions that included both the homoplastic as well as some overlapping microinversions. Neither stem-loop structures nor detectable sequence motifs were associated with microinversions in the hotspots.

CONCLUSIONS

Microinversions can provide valuable phylogenetic information, although power analysis indicates that large amounts of sequence data will be necessary to identify enough inversions (and similar RGCs) to resolve short branches in the tree of life. Moreover, microinversions are not perfect characters and should be interpreted with caution, just as with any other character type. Independent of their use for phylogenetic analyses, microinversions are important because they have the potential to complicate alignment of non-coding sequences. Despite their low rate of accumulation, they have clearly contributed to genome evolution, suggesting that active identification of microinversions will prove useful in future phylogenomic studies.

Evolution of endogenous retroviruses in the Suidae: evidence for different viral subpopulations in African and Eurasian host species

Background

Porcine endogenous retroviruses (PERVs) represent remnants of an exogenous form that have become integrated in the domestic pig (Sus scrofa) genome. Although they are usually inactive, the capacity of γ1 ERVs to infect human cells in vitro has raised concerns about xenotransplantation because the viruses could cross the species barrier to humans. Here we have analyzed the evolution of γ1 ERVs in ten species of Suidae (suids, pigs and hogs) from Eurasia and Africa using DNA sequences for their coding domains (gag, pro/pol and env genes). For comparison with γ1 PERVs, we have also analysed γ2 ERVs which in domestic pigs are known to be inactive and do not pose a risk to xenotransplantation.

Results

Phylogenetic analysis using Bayesian inference showed that γ1 and γ2 ERVs have distinctive evolutionary histories. Firstly, two different viral lineages of γ1 ERVs were found and a coevolutionary analysis demonstrated that they correspond broadly to their host phylogeny, one of Eurasian and another of African species, and show no evidence of horizontal transmission. γ2 ERVs, however, show a bush-like evolution, suggesting a rapid viral radiation from a single common ancestor with no correspondence between host and viral evolutionary trees. Furthermore, though γ1 ERV env genes do not possess frequent stop codons, γ2 env genes do. To understand whether γ1 suid ERVs may be still replicating, we have also evaluated their likely mechanism of proliferation by statistically testing internal to terminal branches using nonsynonymous versus synonymous substitution ratios. Our results suggest that γ1 ERVs are increasing in copy number by reinfection, which requires the translocation of the virus from one cell to another.

Conclusions

Evidence of at least two viral subpopulations was observed in γ1 ERVs from Eurasian and African host species. These results should be taken into account in xenotransplantation since γ1 ERVs appear to be codiverging with their host and maintaining ongoing capacity to infect somatic and germ cells.

Coming to America: multiple origins of New World geckos

Geckos in the Western Hemisphere provide an excellent model to study faunal assembly at a continental scale. We generated a time-calibrated phylogeny, including exemplars of all New World gecko genera, to produce a biogeographical scenario for the New World geckos. Patterns of New World gecko origins are consistent with almost every biogeographical scenario utilized by a terrestrial vertebrate with different New World lineages showing evidence of vicariance, dispersal via temporary land bridge, overseas dispersal or anthropogenic introductions. We also recovered a strong relationship between clade age and species diversity, with older New World lineages having more species than more recently arrived lineages. Our data provide the first phylogenetic hypothesis for all New World geckos and highlight the intricate origins and ongoing organization of continental faunas. The phylogenetic and biogeographical hypotheses presented here provide an historical framework to further pursue research on the diversification and assembly of the New World herpetofauna.

Identifying the rooted species tree from the distribution of unrooted gene trees under the coalescent

Gene trees are evolutionary trees representing the ancestry of genes sampled from multiple populations. Species trees represent populations of individuals-each with many genes-splitting into new populations or species. The coalescent process, which models ancestry of gene copies within populations, is often used to model the probability distribution of gene trees given a fixed species tree. This multispecies coalescent model provides a framework for phylogeneticists to infer species trees from gene trees using maximum likelihood or Bayesian approaches. Because the coalescent models a branching process over time, all trees are typically assumed to be rooted in this setting. Often, however, gene trees inferred by traditional phylogenetic methods are unrooted. We investigate probabilities of unrooted gene trees under the multispecies coalescent model. We show that when there are four species with one gene sampled per species, the distribution of unrooted gene tree topologies identifies the unrooted species tree topology and some, but not all, information in the species tree edges (branch lengths). The location of the root on the species tree is not identifiable in this situation. However, for 5 or more species with one gene sampled per species, we show that the distribution of unrooted gene tree topologies identifies the rooted species tree topology and all its internal branch lengths. The length of any pendant branch leading to a leaf of the species tree is also identifiable for any species from which more than one gene is sampled.

A multigene phylogeny of the fly superfamily Asiloidea (Insecta): Taxon sampling and additional genes reveal the sister-group to all higher flies (Cyclorrhapha)

Asiloidea are a group of 9 lower brachyceran fly families, considered to be the closest relative to the large Metazoan radiation Eremoneura (Cyclorrhapha+Empidoidea). The evidence for asiloid monophyly is limited, and few characters define the relationships between the families of Asiloidea and Eremoneura. Additionally, enigmatic genera, Hilarimorpha and Apystomyia, retain morphological characters of both asiloids and higher flies. We use the nuclear protein-coding gene CAD and 28S rDNA to test the monophyly of Asiloidea and to resolve its relationship to Eremoneura. We explore the effects of taxon sampling on support values and topological stability, the resolving power of additional genes, and hypothesis testing using four-cluster likelihood mapping. We find that: (1) the 'asiloid' genus Apystomyia is sister to Cyclorrhapha, (2) the remaining asiloids are monophyletic at the exclusion of the family Bombyliidae, and (3) our best estimate of relationships places the asiloid flies excluding Bombyliidae as the sister-group to Eremoneura, though high support is lacking.

A multi-locus phylogeny of Nectogalini shrews and influences of the paleoclimate on speciation and evolution

Nectogaline shrews are a major component of the small mammalian fauna of Europe and Asia, and are notable for their diverse ecology, including utilization of aquatic habitats. So far, molecular phylogenetic analyses including nectogaline species have been unable to infer a well-resolved, well-supported phylogeny, thus limiting the power of comparative evolutionary and ecological analyses of the group. Here, we employ Bayesian phylogenetic analyses of eight mitochondrial and three nuclear genes to infer the phylogenetic relationships of nectogaline shrews. We subsequently use this phylogeny to assess the genetic diversity within the genus Episoriculus, and determine whether adaptation to aquatic habitats evolved independently multiple times. Moreover, we both analyze the fossil record and employ Bayesian relaxed clock divergence dating analyses of DNA to assess the impact of historical global climate change on the biogeography of Nectogalini. We infer strong support for the polyphyly of the genus Episoriculus. We also find strong evidence that the ability to heavily utilize aquatic habitats evolved independently in both Neomys and Chimarrogale+Nectogale lineages. Our Bayesian molecular divergence analysis suggests that the early history of Nectogalini is characterized by a rapid radiation at the Miocene/Pliocene boundary, thus potentially explaining the lack of resolution at the base of the tree. Finally, we find evidence that nectogalines once inhabited northern latitudes, but the global cooling and desiccating events at the Miocene/Pliocene and Pliocene/Pleistocene boundaries and Pleistocene glaciation resulted in the migration of most Nectogalini lineages to their present day southern distribution.

Explosive morphological diversification of spiny-finned teleost fishes in the aftermath of the end-Cretaceous extinction

The spiny-finned teleost fishes (Acanthomorpha) include nearly one-third of all living vertebrate species and assume a bewildering array of bodyplans, but the macroevolutionary assembly of modern acanthomorph biodiversity remains largely unexplored. Here, I reconstruct the trajectory of morphological diversification in this major radiation from its first appearance in the Late Cretaceous to the Miocene using a geometric morphometric database comprising more than 600 extinct species known from complete body fossils. The anatomical diversity (disparity) of acanthomorphs is low throughout the Cretaceous, increases sharply and significantly in the wake of the Cretaceous-Palaeogene (K-P) extinction, and shows little change throughout subsequent Cenozoic intervals. This pattern of morphological diversification appears robust to two potential biasing factors: the 'Lagerstätten effect', and the non-random segregation of rare and common taxa along phenotypic axes. Dissecting the trajectory of acanthomorph radiation along phylogenetic lines reveals that the abrupt post-extinction increase in disparity is driven largely by the proliferation of trophically diverse modern groups within Percomorpha, a spiny-fin subclade containing more than 15 000 living species and identified as showing a substantially elevated diversification rate relative to background vertebrate levels. A major component of the Palaeogene acanthomorph radiation reflects colonization of morphospace previously occupied by non-acanthomorph victims of the K-P. However, other aspects of morphological diversification cannot be explained by this simple ecological release model, suggesting that multiple factors contributed to the prolific anatomical radiation of acanthomorphs.

Phylogenetic relationships among families of Gadiformes (Teleostei, Paracanthopterygii) based on nuclear and mitochondrial data

Phylogenetic hypotheses among Gadiformes fishes at the suborder, family, and subfamily levels are controversial. To address this problem, we analyze nuclear and mitochondrial DNA (mtDNA) sequences for the most extensive taxonomic sampling compiled to date, representing all of the recognized families and subfamilies in the order (except the monotypic family Lyconidae). Our study sampled 117 species from 46 genera, comprising around 20% of the species described for the order (more than 60% of all genera in the order) and produced 2740 bp of DNA sequence data for each species. Our analysis was successful in confirming the monophyly of Gadiformes and most of the proposed families for the order, but alternative hypotheses of sister-group relationships among families were poorly resolved. Our results are consistent with dividing Gadiformes into 12 families in three suborders, Muraenolepidoidei, Macrouroidei, and Gadoidei. Muraenolepidoidei contains the single family Muraenolepididae. The suborder Macrouroidei includes at least three families: Macrouridae, Macruronidae and Steindachneriidae. Macrouridae is deeply divided into two well-supported subfamilies: Macrourinae and Bathygadinae, suggesting that Bathygadinae may be ranked at the family level. The suborder Gadoidei includes the families: Merlucciidae, Melanonidae, Euclichthyidae, Gadidae, Ranicipitidae, and Bregmacerotidae. Additionally, Trachyrincinae could be ranked at family level including two subfamilies: Trachyrincinae and Macrouroidinae within Gadoidei. Further taxonomic sampling and sequencing efforts are needed in order to corroborate these relationships.

Seven new dolphin mitochondrial genomes and a time-calibrated phylogeny of whales

BACKGROUND

The phylogeny of Cetacea (whales) is not fully resolved with substantial support. The ambiguous and conflicting results of multiple phylogenetic studies may be the result of the use of too little data, phylogenetic methods that do not adequately capture the complex nature of DNA evolution, or both. In addition, there is also evidence that the generic taxonomy of Delphinidae (dolphins) underestimates its diversity. To remedy these problems, we sequenced the complete mitochondrial genomes of seven dolphins and analyzed these data with partitioned Bayesian analyses. Moreover, we incorporate a newly-developed "relaxed" molecular clock to model heterogenous rates of evolution among cetacean lineages.

RESULTS

The "deep" phylogenetic relationships are well supported including the monophyly of Cetacea and Odontoceti. However, there is ambiguity in the phylogenetic affinities of two of the river dolphin clades Platanistidae (Indian River dolphins) and Lipotidae (Yangtze River dolphins). The phylogenetic analyses support a sister relationship between Delphinidae and Monodontidae + Phocoenidae. Additionally, there is statistically significant support for the paraphyly of Tursiops (bottlenose dolphins) and Stenella (spotted dolphins).

CONCLUSION

Our phylogenetic analysis of complete mitochondrial genomes using recently developed models of rate autocorrelation resolved the phylogenetic relationships of the major Cetacean lineages with a high degree of confidence. Our results indicate that a rapid radiation of lineages explains the lack of support the placement of Platanistidae and Lipotidae. Moreover, our estimation of molecular divergence dates indicates that these radiations occurred in the Middle to Late Oligocene and Middle Miocene, respectively. Furthermore, by collecting and analyzing seven new mitochondrial genomes, we provide strong evidence that the delphinid genera Tursiops and Stenella are not monophyletic, and the current taxonomy masks potentially interesting patterns of morphological, physiological, behavioral, and ecological evolution.

Reconciliation with non-binary species trees

Reconciliation extracts information from the topological incongruence between gene and species trees to infer duplications and losses in the history of a gene family. The inferred duplication-loss histories provide valuable information for a broad range of biological applications, including ortholog identification, estimating gene duplication times, and rooting and correcting gene trees. While reconciliation for binary trees is a tractable and well studied problem, there are no algorithms for reconciliation with non-binary species trees. Yet a striking proportion of species trees are non-binary. For example, 64% of branch points in the NCBI taxonomy have three or more children. When applied to non-binary species trees, current algorithms overestimate the number of duplications because they cannot distinguish between duplication and incomplete lineage sorting. We present the first algorithms for reconciling binary gene trees with non-binary species trees under a duplication-loss parsimony model. Our algorithms utilize an efficient mapping from gene to species trees to infer the minimum number of duplications in O(|V(G) | x (k(S) + h(S))) time, where |V(G)| is the number of nodes in the gene tree, h(S) is the height of the species tree and k(S) is the size of its largest polytomy. We present a dynamic programming algorithm which also minimizes the total number of losses. Although this algorithm is exponential in the size of the largest polytomy, it performs well in practice for polytomies with outdegree of 12 or less. We also present a heuristic which estimates the minimal number of losses in polynomial time. In empirical tests, this algorithm finds an optimal loss history 99% of the time. Our algorithms have been implemented in NOTUNG, a robust, production quality, tree-fitting program, which provides a graphical user interface for exploratory analysis and also supports automated, high-throughput analysis of large data sets.

Phylogenetic distributions of British birds of conservation concern

Recent studies suggest that species' life histories and ecology can be used to forecast future extinction risk. Threatened species often share similar traits such that if a trait predisposing a species to decline or extinction is evolutionarily conserved, then close relatives of threatened species are themselves likely to be at risk. The phylogenetic distribution of current threat has been argued to provide insight into the species that could be threatened in the future when trait data are not available. Conservation criteria are typically based on multiple indices that capture different symptoms of threat including population trends and range contraction. However, there is no reason to assume consistent phylogenetic distributions of different symptoms. I construct a molecular phylogeny of 249 species of British birds (more than 93% of the breeding and wintering species) and use this to show that the species that are threatened due to population declines are phylogenetically more closely related than expected by chance alone. However, species that are listed for other reasons, including range contraction, are distributed randomly with respect to phylogeny. I suggest that while phylogeny can be informative with respect to identifying clades that are susceptible to some measures of extinction risk, such patterns are likely to be idiosyncratic with respect to symptom and taxa.

The anatomical relationships between the avian eye, orbit and sclerotic ring: implications for inferring activity patterns in extinct birds

Activity pattern, or the time of day when an animal is awake and active, is highly associated with that animal's ecology. There are two principal activity patterns: diurnal, or awake during the day in a photopic, or high light level, environment; and nocturnal, awake at night in scotopic, or low light level, conditions. Nocturnal and diurnal birds exhibit characteristic eye shapes associated with their activity pattern, with nocturnal bird eyes optimized for visual sensitivity with large corneal diameters relative to their eye axial lengths, and diurnal birds optimized for visual acuity, with larger axial lengths of the eye relative to their corneal diameters. The current study had three aims: (1) to quantify the nature of the relationship between the avian eye and its associated bony anatomy, the orbit and the sclerotic ring; (2) to investigate how activity pattern is reflected in that bony anatomy; and (3) to identify how much bony anatomy is required to interpret activity pattern reliably for a bird that does not have the soft tissue available for study, specifically, for a fossil. Knowledge of extinct avian activity patterns would be useful in making palaeoecological interpretations. Here eye, orbit and sclerotic ring morphologies of 140 nocturnal and diurnal bird species are analysed in a phylogenetic context. Although there is a close relationship between the avian eye and orbit, activity pattern can only be reliably interpreted for bony-only specimens, such as a fossil, that include both measurements of the sclerotic ring and orbit depth. Any missing data render the fossil analysis inaccurate, including fossil specimens that are flat and therefore do not have an orbit depth available. For example, activity pattern cannot be determined with confidence for Archaeopteryx lithographica, which has a complete sclerotic ring but no orbit depth measurement. Many of the bird fossils currently available that retain a good sclerotic ring tend to be flat specimens, while three-dimensionally preserved bird fossils tend not to have a well-preserved sclerotic ring or a well-defined optic foramen, necessary for delimiting the orbit depth.