THE INFLUENCE OF SAMPLING DESIGN ON SPECIES TREE INFERENCE: A NEW RELATIONSHIP FOR THE NEW WORLD CHICKADEES (AVES: POECILE)

In this study, we explore the long‐standing issue of how many loci are needed to infer accurate phylogenetic relationships, and whether loci with particular attributes (e.g., parsimony informativeness, variability, gene tree resolution) outperform others. To do so, we use an empirical data set consisting of the seven species of chickadees (Aves: Paridae), an analytically tractable, recently diverged group, and well‐studied ecologically but lacking a nuclear phylogeny. We estimate relationships using 40 nuclear loci and mitochondrial DNA using four coalescent‐based species tree inference methods (BEST, *BEAST, STEM, STELLS). Collectively, our analyses contrast with previous studies and support a sister relationship between the Black‐capped and Carolina Chickadee, two superficially similar species that hybridize along a long zone of contact. Gene flow is a potential source of conflict between nuclear and mitochondrial gene trees, yet we find a significant, albeit low, signal of gene flow. Our results suggest that relatively few loci with high information content may be sufficient for estimating an accurate species tree, but that substantially more loci are necessary for accurate parameter estimation. We provide an empirical reference point for researchers designing sampling protocols with the purpose of inferring phylogenies and population parameters of closely related taxa.     

Comparison of Y chromosome and mtDNA phylogenies leads to unique inferences of macaque evolutionary history

We report here the results of one of the first analyses to use male-specific nuclear markers in elucidating primate phylogenetic relationships at the intrageneric level. Two closely linked Y chromosome markers, TSPY and SRY, were sequenced for a total of 3100 bases. Forty-four macaques, representing 18 of the 19 recognized species, were sequenced for the full 3.1 kb, as was 1 individual from each of the following outgroup genera: Papio, Theropithecus, Mandrillus, Allenopithecus,Cercopithecus, Trachypithecus, Presbytis, and Homo. In contrast to recent mtDNA phylogenies, Y chromosome loci support four monophyletic species groups, including a sinica group containing M. arctoides-a classification largely congruent with those of Fooden and Delson. Comparison of mtDNA and Y chromosome phylogenies highlight (1) a potential hybrid origin of Macaca arctoides from M. fascicularis and proto-M. assamensis/thibetana and (2) cases of mitochondrial paraphyly in macaque species whose Y chromosome lineages are monophyletic-a probable evolutionary consequence of philopatric females vs dispersing males. These results raise the question of whether a phylogenetic tree should be a topology of species origins or a depiction of more current species relationships, including subsequent episodes of introgression.     

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.     

Gene tree parsimony of multilocus snake venom protein families reveals species tree conflict as a result of multiple parallel gene loss

The proliferation of gene data from multiple loci of large multigene families has been greatly facilitated by considerable recent advances in sequence generation. The evolution of such gene families, which often undergo complex histories and different rates of change, combined with increases in sequence data, pose complex problems for traditional phylogenetic analyses, and in particular, those that aim to successfully recover species relationships from gene trees. Here, we implement gene tree parsimony analyses on multicopy gene family data sets of snake venom proteins for two separate groups of taxa, incorporating Bayesian posterior distributions as a rigorous strategy to account for the uncertainty present in gene trees. Gene tree parsimony largely failed to infer species trees congruent with each other or with species phylogenies derived from mitochondrial and single-copy nuclear sequences. Analysis of four toxin gene families from a large expressed sequence tag data set from the viper genus Echis failed to produce a consistent topology, and reanalysis of a previously published gene tree parsimony data set, from the family Elapidae, suggested that species tree topologies were predominantly unsupported. We suggest that gene tree parsimony failure in the family Elapidae is likely the result of unequal and/or incomplete sampling of paralogous genes and demonstrate that multiple parallel gene losses are likely responsible for the significant species tree conflict observed in the genus Echis. These results highlight the potential for gene tree parsimony analyses to be undermined by rapidly evolving multilocus gene families under strong natural selection.     

Species tree inference in a recent radiation of orioles (Genus Icterus): multiple markers and methods reveal cytonuclear discordance in the northern oriole group

Recent computational advances provide novel opportunities to infer species trees based on multiple independent loci. Thus, single gene trees no longer need suffice as proxies for species phylogenies. Several methods have been developed to deal with the challenges posed by incomplete and stochastic lineage sorting. In this study, we employed four Bayesian methods to infer the phylogeny of a clade of 11 recently diverged oriole species within the genus Icterus. We obtained well-resolved and mostly congruent phylogenies using a set of seven unlinked nuclear intron loci and sampling multiple individuals per species. Most notably, Bayesian concordance analysis generally agreed well with concatenation; the two methods agreed fully on eight of nine nodes. The coalescent-based method ^∗BEAST further supported six of these eight nodes. The fourth method used, BEST, failed to converge despite exhaustive efforts to optimize the tree search. Overall, the results obtained by new species tree methods and concatenation generally corroborate our findings from previous analyses and data sets. However, we found striking disagreement between mitochondrial and nuclear DNA involving relationships within the northern oriole group. Our results highlight the danger of reliance on mtDNA alone for phylogenetic inference. We demonstrate that in spite of low variability and incomplete lineage sorting, multiple nuclear loci can produce largely congruent phylogenies based on multiple species tree methods, even for very closely-related species.     

A comparative analysis of Y chromosome and mtDNA phylogenies of the Hylobates gibbons

ABSTRACT: BACKGROUND: The evolutionary relationships of closely related species have long been of interest to biologists since these species experienced different evolutionary processes in a relatively short period of time. Comparison of phylogenies inferred from DNA sequences with differing inheritance patterns, such as mitochondrial, autosomal, and X and Y chromosomal loci, can provide more comprehensive inferences of the evolutionary histories of species. Gibbons, especially the genus Hylobates, are particularly intriguing as they consist of multiple closely related species which emerged rapidly and live in close geographic proximity. Our current understanding of relationships among Hylobates species is largely based on data from the maternally-inherited mitochondrial DNAs (mtDNAs). RESULTS: To infer the paternal histories of gibbon taxa, we sequenced multiple Y chromosomal loci from 26 gibbons representing 10 species. As expected, we find levels of sequence variation some five times lower than observed for the mitochondrial genome (mtgenome). Although our Y chromosome phylogenetic tree shows relatively low resolution compared to the mtgenome tree, our results are consistent with the monophyly of gibbon genera suggested by the mtgenome tree. In a comparison of the molecular dating of divergences and on the branching patterns of phylogeny trees between mtgenome and Y chromosome data, we found: 1) the inferred divergence estimates were more recent for the Y chromosome than for the mtgenome, 2) the species H. lar and H. pileatus are reciprocally monophyletic in the mtgenome phylogeny but a H. pileatus individual falls into the H. lar Y chromosome clade. CONCLUSIONS: Based on the ~6.4 kb of Y chromosomal DNA sequence data generated for each of the 26 individuals in this study, we provide molecular inferences on gibbon and particularly on Hylobates evolution complementary to those from mtDNA data. Overall, our results illustrate the utility of comparative studies of loci with different inheritance patterns for investigating potential sex specific processes on the evolutionary histories of closely related taxa, and emphasize the need for further sampling of gibbons of known provenance.     

Discordant mtDNA and cpDNA phylogenies indicate geographic speciation and reticulation as driving factors for the diversification of the genus Picea

The phylogeny of the genus Picea was investigated by sequencing three loci from the paternally inherited chloroplast genome (trnK, rbcL and trnTLF) and the intron 2 of the maternally transmitted mitochondrial gene nad1 for 35 species. Significant topological differences were found between the trnK tree and the rbcL and trnTLF phylogenetic trees, and between cpDNA and mtDNA phylogenies. None of the phylogenies matched morphological classifications. The mtDNA phylogeny was geographically more structured than cpDNA phylogenies, reflecting the different inheritance of the two cytoplasmic genomes in the Pinaceae and their differential dispersion by seed only and seed and pollen, respectively. Most North American taxa formed a monophyletic group on the mtDNA tree, with topological patterns suggesting geographic speciation by range fragmentation or by dispersal and isolation. Similar patterns were also found among Asian taxa. Such a trend towards geographic speciation is anticipated in other Pinaceae genera with similar life history, autecology and reproductive system. Incongruences between organelle phylogenies suggested the occurrence of mtDNA capture by invading cpDNA. Incongruences between cpDNA partitions further suggested heterologous recombination presumably also linked to ancient reticulate evolution. Whilst cpDNA appears potentially valuable for molecular taxonomy and systematics purposes, these results emphasize the reduced value of cpDNA to infer vertical descent and the speciation history for plants with paternal transmission and high dispersal of their chloroplast genome.     

An improved phylogeny of the Andean tit-tyrants (Aves, Tyrannidae): more characters trump sophisticated analyses

The phylogeny of the flycatcher genus Anairetes was previously inferred using short fragments of mitochondrial DNA and parsimony and distance-based methods. The resulting topology spurred taxonomic revision and influenced understanding of Andean biogeography. More than a decade later, we revisit the phylogeny of Anairetes tit-tyrants using more mtDNA characters, seven unlinked loci (three mitochondrial genes, six nuclear loci), more closely related outgroup taxa, partitioned Bayesian analyses, and two coalescent species-tree approaches (Bayesian estimation of species trees, BEST; Bayesian evolutionary analysis by sampling trees, (*)BEAST). Of these improvements in data and analyses, the fourfold increase in mtDNA characters was both necessary and sufficient to incur a major shift in the topology and near-complete resolution. The species-tree analyses, while theoretically preferable to concatenation or single gene approaches, yielded topologies that were compatible with mtDNA but with weaker statistical resolution at nodes. The previous results that had led to taxonomic and biogeographic reappraisal were refuted, and the current results support the resurrection of the genus Uromyias as the sister clade to Anairetes. The sister relationship between these two genera corresponds to an ecological dichotomy between a depauperate humid cloud forest clade and a diverse dry-tolerant clade that has diversified along the latitudinal axis of the Andes. The species-tree results and the concatenation results each reaffirm the primacy of mtDNA to provide phylogenetic signal for avian phylogenies at the species and subspecies level. This is due in part to the abundance of informative characters in mtDNA, and in part to its lower effective population size that causes it to more faithfully track the species tree.     

Evolutionary history of the odd-nosed monkeys and the phylogenetic position of the newly described Myanmar snub-nosed monkey Rhinopithecus strykeri

Odd-nosed monkeys represent one of the two major groups of Asian colobines. Our knowledge about this primate group is still limited as it is highlighted by the recent discovery of a new species in Northern Myanmar. Although a common origin of the group is now widely accepted, the phylogenetic relationships among its genera and species, and the biogeographic processes leading to their current distribution are largely unknown. To address these issues, we have analyzed complete mitochondrial genomes and 12 nuclear loci, including one X chromosomal, six Y chromosomal and five autosomal loci, from all ten odd-nosed monkey species. The gene tree topologies and divergence age estimates derived from different markers were highly similar, but differed in placing various species or haplogroups within the genera Rhinopithecus and Pygathrix. Based on our data, Rhinopithecus represent the most basal lineage, and Nasalis and Simias form closely related sister taxa, suggesting a Northern origin of odd-nosed monkeys and a later invasion into Indochina and Sundaland. According to our divergence age estimates, the lineages leading to the genera Rhinopithecus, Pygathrix and Nasalis+Simias originated in the late Miocene, while differentiation events within these genera and also the split between Nasalis and Simias occurred in the Pleistocene. Observed gene tree discordances between mitochondrial and nuclear datasets, and paraphylies in the mitochondrial dataset for some species of the genera Rhinopithecus and Pygathrix suggest secondary gene flow after the taxa initially diverged. Most likely such events were triggered by dramatic changes in geology and climate within the region. Overall, our study provides the most comprehensive view on odd-nosed monkey evolution and emphasizes that data from differentially inherited markers are crucial to better understand evolutionary relationships and to trace secondary gene flow.     

Thermal biology of African lovebirds and Australian grass parakeets

We compared aspects of the thermal biology of two groups of small parrots, of similar body mass, each derived from a range of habitat types, varying in aridity, but indigenous to either southern Africa or Australia. By accounting for phylogenetic differences, we were able to question whether arid zone species have lower metabolic rates and greater thermal tolerances than mesic species in relation to the “pre-adapted” and “post-arrival adaptation” hypotheses. Four species of African lovebird (Agapornis) and four species of Australian grass parakeet (one Neopsephotus and three Neophema species) were investigated. The Rosy-faced Lovebird (Agapornis roseicollis), Bourke's Parakeet (Neopsephotus bourkii) and the Scarlet-chested Parakeet (Neophema splendida) were categorised as arid zone species, Fischer's Lovebird (Agapornis fischeri), the Black-masked Lovebird (Agapornis personatus) and the Elegant Parakeet (Neophema elegans) as semi-arid zone species, and the Black-cheeked Lovebird (Agapornis nigrigenis) and the Turquoise Parakeet (Neophema pulchella) as mesic zone species. Conventional and phylogenetically independent statistical methods yielded no significant differences in the basal metabolic rates of birds from different habitats or between the species assemblages from Africa and Australia.     

A Total Evidence Phylogeny of the Arctoidea (Carnivora: Mammalia): Relationships Among Basal Taxa

A total evidence phylogenetic analysis was performed for 14 extant and 18 fossil caniform genera using a data matrix of 5.6 kbp of concatenated sequence data from six independent loci and 80 morphological characters from the cranium and dentition. Maximum parsimony analysis recovered a single most parsimonious cladogram (MPC). The topology of the extant taxa in the MPC agreed with previous molecular phylogenies. Phylogenetic positions for fossil taxa indicate that several taxa previously described as early members of extant families (e.g., Bathygale and Plesictis) are likely stem taxa at the base of the Arctoidea. Taxa in the “Paleomustelidae” were found to be paraphyletic, but a monophyletic Oligobuninae was recovered within this set of taxa. This clade was closely related to the extant genera Gulo and Martes, therefore, nested within the extant radiation of the family Mustelidae. This analysis provides a resolution to several discrepancies between phylogenies considering either fossil taxa or extant taxa separately, and provides a framework for incorporating fossil and extant taxa into comprehensive combined evidence analyses.     

Multilocus analysis of a taxonomically densely sampled dataset reveal extensive non-monophyly in the avian family Locustellidae

The phylogeny of most of the species in the avian passerine family Locustellidae is inferred using a Bayesian species tree approach (Bayesian Estimation of Species Trees, BEST), as well as a traditional Bayesian gene tree method (MrBayes), based on a dataset comprising one mitochondrial and four nuclear loci. The trees inferred by the different methods agree fairly well in topology, although in a few cases there are marked differences. Some of these discrepancies might be due to convergence problems for BEST (despite up to 1×10(9) iterations). The phylogeny strongly disagrees with the current taxonomy at the generic level, and we propose a revised classification that recognizes four instead of seven genera. These results emphasize the well known but still often neglected problem of basing classifications on non-cladistic evaluations of morphological characters. An analysis of an extended mitochondrial dataset with multiple individuals from most species, including many subspecies, suggest that several taxa presently treated as subspecies or as monotypic species as well as a few taxa recognized as separate species are in need of further taxonomic work.     

Total evidence: molecules, morphology, and the phylogenetics of cichlid fishes

We present a most comprehensive phylogenetic analysis of the family Cichlidae. New data analyzed include mitochondrial 16S rRNA sequences and two nuclear loci (Tmo-M27 and Tmo-4C4) for a large taxonomic sampling with emphasis on South American species. We also incorporate a published morphological data set for a total evidence analysis. Character congruence among mitochondrial (74 taxa) and nuclear data (50 taxa) was high. However, partition-homogeneity tests suggest significant heterogeneity among molecular and morphological data. In agreement with results obtained from molecular data alone, total evidence analysis (1,460 characters for 34 taxa) supports a robust phylogenetic hypothesis for the family Cichlidae that is congruent with drift-vicariance events associated with the fragmentation of Gondwana. Our analyses confirm the placement of Malagasy/Indian cichlids as the most basal lineages, with a sister-group relationship to the monophyletic African and Neotropical clades. Total evidence suggests that the controversial African genus Heterochromis is at the base of the African radiation. Among more than 50 Neotropical genera analyzed, Retroculus is identified as the basal taxon, with successive branching of Cichla, Astronotus, geophagines (including crenicichlines) + chaetobranchines, and cichlasomines + heroines. Relative rate tests applied to mitochondrial DNA suggest significantly higher rates of genetic variation in Neotropical than in African taxa, and both mitochondrial and nuclear sequences show that rate heterogeneity among Neotropical lineages is confined to the geophagine cichlids.     

Discordant mitochondrial and nuclear gene phylogenies in emydid turtles: implications for speciation and conservation

Do phylogenies and branch lengths based on mitochondrial DNA (mtDNA) provide a reasonable approximation to those based on multiple nuclear loci? In the present study, we show widespread discordance between phylogenies based on mtDNA (two genes) and nuclear DNA (nucDNA; six loci) in a phylogenetic analysis of the turtle family Emydidae. We also find an unusual type of discordance involving the unexpected homogeneity of mtDNA sequences across species within genera. Of the 36 clades in the combined nucDNA phylogeny, 24 are contradicted by the mtDNA phylogeny, and six are strongly contested by each data set. Two genera (Graptemys, Pseudemys) show remarkably low mtDNA divergence among species, whereas the combined nuclear data show deep divergences and (for Pseudemys) strongly supported clades. These latter results suggest that the mitochondrial data alone are highly misleading about the rate of speciation in these genera and also about the species status of endangered Graptemys and Pseudemys species. In addition, despite a strongly supported phylogeny from the combined nuclear genes, we find extensive discordance between this tree and individual nuclear gene trees. Overall, the results obtained illustrate the potential dangers of making inferences about phylogeny, speciation, divergence times, and conservation from mtDNA data alone (or even from single nuclear genes), and suggest the benefits of using large numbers of unlinked nuclear loci. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 445–461.     

A multilocus molecular phylogeny of the parrots (Psittaciformes): support for a Gondwanan origin during the cretaceous

The question of when modern birds (Neornithes) first diversified has generated much debate among avian systematists. Fossil evidence generally supports a Tertiary diversification, whereas estimates based on molecular dating favor an earlier diversification in the Cretaceous period. In this study, we used an alternate approach, the inference of historical biogeographic patterns, to test the hypothesis that the initial radiation of the Order Psittaciformes (the parrots and cockatoos) originated on the Gondwana supercontinent during the Cretaceous. We utilized broad taxonomic sampling (representatives of 69 of the 82 extant genera and 8 outgroup taxa) and multilocus molecular character sampling (3,941 bp from mitochondrial DNA (mtDNA) genes cytochrome oxidase I and NADH dehydrogenase 2 and nuclear introns of rhodopsin intron 1, tropomyosin alpha-subunit intron 5, and transforming growth factor ss-2) to generate phylogenetic hypotheses for the Psittaciformes. Analyses of the combined character partitions using maximum parsimony, maximum likelihood, and Bayesian criteria produced well-resolved and topologically similar trees in which the New Zealand taxa Strigops and Nestor (Psittacidae) were sister to all other psittaciforms and the cockatoo clade (Cacatuidae) was sister to a clade containing all remaining parrots (Psittacidae). Within this large clade of Psittacidae, some traditionally recognized tribes and subfamilies were monophyletic (e.g., Arini, Psittacini, and Loriinae), whereas several others were polyphyletic (e.g., Cyclopsittacini, Platycercini, Psittaculini, and Psittacinae). Ancestral area reconstructions using our Bayesian phylogenetic hypothesis and current distributions of genera supported the hypothesis of an Australasian origin for the Psittaciformes. Separate analyses of the timing of parrot diversification constructed with both Bayesian relaxed-clock and penalized likelihood approaches showed better agreement between geologic and diversification events in the chronograms based on a Cretaceous dating of the basal split within parrots than the chronograms based on a Tertiary dating of this split, although these data are more equivocal. Taken together, our results support a Cretaceous origin of Psittaciformes in Gondwana after the separation of Africa and the India/Madagascar block with subsequent diversification through both vicariance and dispersal. These well-resolved molecular phylogenies will be of value for comparative studies of behavior, ecology, and life history in parrots.     

Differentiating between hypotheses of lineage sorting and introgression in New Zealand alpine cicadas (Maoricicada Dugdale)

Lineage sorting and introgression can lead to incongruence among gene phylogenies, complicating the inference of species trees for large groups of taxa that have recently and rapidly radiated. In addition, it can be difficult to determine which of these processes is responsible for this incongruence. We explore these issues with the radiation of New Zealand alpine cicadas of the genus Maoricicada Dugdale. Gene trees were estimated from four putative independent loci: mitochondrial DNA (2274 nucleotides), elongation factor 1-alpha (1275 nucleotides), period (1709 nucleotides), and calmodulin (678 nucleotides). We reconstructed phylogenies using maximum likelihood and Bayesian methods from 44 individuals representing the 19 species and subspecies of Maoricicada and two outgroups. Species-level relationships were reconstructed using a novel extension of gene tree parsimony, whereby gene trees were weighted by their Bayesian posterior probabilities. The inferred gene trees show marked incongruence in the placement of some taxa, especially the enigmatic forest and scrub dwelling species, M. iolanthe. Using the species tree estimated by gene tree parsimony, we simulated coalescent gene trees in order to test the null hypothesis that the nonrandom placement of M. iolanthe among gene trees has arisen by chance. Under the assumptions of constant population size, known generation time, and panmixia, we were able to reject this null hypothesis. Furthermore, because the two alternative placements of M. iolanthe are in each case with species that share a similar song structure, we conclude that it is more likely that an ancient introgression event rather than lineage sorting has caused this incongruence.     

Multigene phylogeny of the Old World mice, Murinae, reveals distinct geographic lineages and the declining utility of mitochondrial genes compared to nuclear genes

Despite its great diversity and biomedical importance, the rodent subfamily Murinae is poorly resolved phylogenetically. We present the first cladistic analysis sampling multiple representatives of most major groups based on DNA sequence for three nuclear (GHR, RAG1, and AP5) and one mitochondrial (COII and parts of COI and ATPase 8) fragments. Analyzed separately, the four partitions agree broadly with each other and the combined analysis. The basal split is between a clade of Philippine Old Endemics and all remaining murines. Within the latter, rapid radiation led to at least seven geographically distinct lineages, including a Southeast Asian Rattus clade; a diverse Australo-Papuan and Philippine clade; an African arvicanthine group including the otomyines; an African Praomys group; and three independent genera from Africa and Asia, Mus, Apodemus, and Malacomys. The murines appear to have originated in Southeast Asia and then rapidly expanded across all of the Old World. Both nuclear exons provide robust support at all levels. In contrast, the bootstrap proportions from mitochondrial data decline rapidly with increasing depth in the tree, together suggesting that nuclear genes may be more useful even for relatively recent divergences (< 10MYA).     

A hitchhikers guide to the Galápagos: co-phylogeography of Galápagos mockingbirds and their parasites

Background

Parasites are evolutionary hitchhikers whose phylogenies often track the evolutionary history of their hosts. Incongruence in the evolutionary history of closely associated lineages can be explained through a variety of possible events including host switching and host independent speciation. However, in recently diverged lineages stochastic population processes, such as retention of ancestral polymorphism or secondary contact, can also explain discordant genealogies, even in fully co-speciating taxa. The relatively simple biogeographic arrangement of the Galápagos archipelago, compared with mainland biomes, provides a framework to identify stochastic and evolutionary informative components of genealogic data in these recently diverged organisms.

Results

Mitochondrial DNA sequences were obtained for four species of Galápagos mockingbirds and three sympatric species of ectoparasites - two louse and one mite species. These data were complemented with nuclear EF1α sequences in selected samples of parasites and with information from microsatellite loci in the mockingbirds. Mitochondrial sequence data revealed differences in population genetic diversity between all taxa and varying degrees of topological congruence between host and parasite lineages. A very low level of genetic variability and lack of congruence was found in one of the louse parasites, which was excluded from subsequent joint analysis of mitochondrial data. The reconciled multi-species tree obtained from the analysis is congruent with both the nuclear data and the geological history of the islands.

Conclusions

The gene genealogies of Galápagos mockingbirds and two of their ectoparasites show strong phylogeographic correlations, with instances of incongruence mostly explained by ancestral genetic polymorphism. A third parasite genealogy shows low levels of genetic diversity and little evidence of co-phylogeny with their hosts. These differences can mostly be explained by variation in life-history characteristics, primarily host specificity and dispersal capabilities. We show that pooling genetic data from organisms living in close ecological association reveals a more accurate phylogeographic history for these taxa. Our results have implications for the conservation and taxonomy of Galápagos mockingbirds and their parasites.     

Implementing a cumulative supermatrix approach for a comprehensive phylogenetic study of the Teloschistales (Pezizomycotina, Ascomycota)

The resolution of the phylogenetic relationships within the order Teloschistales (Ascomycota, lichen-forming-fungi), with nearly 2000 known species and outstanding phenotypic diversity, has been hindered by the limitation in the resolving power that single-locus or two-locus phylogenetic studies have provided to date. In this context, an extensive taxon sampling within the Teloschistales with more loci (especially nuclear protein-coding genes) was needed to confront the current taxonomic delimitations and to understand evolutionary trends within this order. Comprehensive maximum likelihood and bayesian analyses were performed based on seven loci using a cumulative supermatrix approach, including protein-coding genes RPB1 and RPB2 in addition to nuclear and mitochondrial ribosomal RNA-coding genes. We included 167 taxa representing 12 of the 15 genera recognized within the currently accepted Teloschistineae, 22 of the 43 genera within the Physciineae, 49 genera of the closely related orders Lecanorales, Lecideales, and Peltigerales, and the dubiously placed family Brigantiaeaceae and genus Sipmaniella. Although the progressive addition of taxa (cumulative supermatrix approach) with increasing amounts of missing data did not dramatically affect the loss of support and resolution, the monophyly of the Teloschistales in the current sense was inconsistent, depending on the loci-taxa combination analyzed. Therefore, we propose a new, but provisional, classification for the re-circumscribed orders Caliciales and Teloschistales (previously referred to as Physciineae and Teloschistineae, respectively). We report here that the family Brigantiaeaceae, previously regarded as incertae sedis within the subclass Lecanoromycetidae, and Sipmaniella, are members of the Teloschistales in a strict sense. Within this order, one lineage led to the diversification of the mostly epiphytic crustose Brigantiaeaceae and Letrouitiaceae, with a circumpacific center of diversity and found mostly in the tropics. The other main lineage led to another epiphytic crustose family, mostly tropical, and with an Australasian center of diversity--the Megalosporaceae--which is sister to the mainly rock-inhabiting, cosmopolitan, and species rich Teloschistaceae, with a diversity of growth habits ranging from crustose to fruticose. Our results confirm the use of a cumulative supermatrix approach as a viable method to generate comprehensive phylogenies summarizing relationships of taxa with multi-locus to single locus data.     

Impact of deep coalescence on the reliability of species tree inference from different types of DNA markers in mammals

An important challenge for phylogenetic studies of closely related species is the existence of deep coalescence and gene tree heterogeneity. However, their effects can vary between species and they are often neglected in phylogenetic analyses. In addition, a practical problem in the reconstruction of shallow phylogenies is to determine the most efficient set of DNA markers for a reliable estimation. To address these questions, we conducted a multilocus simulation study using empirical values of nucleotide diversity and substitution rates obtained from a wide range of mammals and evaluated the performance of both gene tree and species tree approaches to recover the known speciation times and topological relationships. We first show that deep coalescence can be a serious problem, more than usually assumed, for the estimation of speciation times in mammals using traditional gene trees. Furthermore, we tested the performance of different sets of DNA markers in the determination of species trees using a coalescent approach. Although the best estimates of speciation times were obtained, as expected, with the use of an increasing number of nuclear loci, our results show that similar estimations can be obtained with a much lower number of genes and the incorporation of a mitochondrial marker, with its high information content. Thus, the use of the combined information of both nuclear and mitochondrial markers in a species tree framework is the most efficient option to estimate recent speciation times and, consequently, the underlying species tree.