Phylogenomic analysis of Stylops reveals the evolutionary history of a Holarctic Strepsiptera radiation parasitizing wild bees.

Holarctic Stylops is the largest genus of the enigmatic insect order Strepsiptera, twisted winged parasites. Members of Stylops are obligate endoparasites of Andrena mining bees and exhibit extreme sexual dimorphism typical of Strepsiptera. So far, molecular studies on Stylops have focused on questions on species delimitation. Here, we utilize the power of whole genome sequencing to infer the phylogeny of this morphologically challenging genus from thousands of loci. We use a species tree method, concatenated maximum likelihood analysis and Bayesian analysis with a relaxed clock model to reconstruct the phylogeny of 46 Stylops species, estimate divergence times, evaluate topological consistency across methods and infer the root position. Furthermore, the biogeographical history and coevolutionary patterns with host species are assessed. All methods recovered a well resolved topology with close to all nodes maximally supported and only a handful of minor topological variations. Based on the result, we find that included species can be divided into 12 species groups, seven of them including only Palaearctic species, three Nearctic and two were geographically mixed. We find a strongly supported root position between a clade formed by the spreta, thwaitesi and gwynanae species groups and the remaining species and that the sister group of Stylops is Eurystylops or Eurystylops + Kinzelbachus. Our results indicate that Stylops originated in the Western Palaearctic or Western Palaearctic and Nearctic in the early Neogene or late Paleogene, with four independent dispersal events to the Nearctic. Cophylogenetic analyses indicate that the diversification of Stylops has been shaped by both significant coevolution with the mining bee hosts and host-shifting. The well resolved and strongly supported phylogeny will provide a valuable phylogenetic basis for further studies into the fascinating world of Strepsipterans.

The formation of the Indo-Pacific montane avifauna.

The processes generating the earth's montane biodiversity remain a matter of debate. Two contrasting hypotheses have been advanced to explain how montane populations form: via direct colonization from other mountains, or, alternatively, via upslope range shifts from adjacent lowland areas. We seek to reconcile these apparently conflicting hypotheses by asking whether a species' ancestral geographic origin determines its mode of mountain colonization. Island-dwelling passerine birds at the faunal crossroads between Eurasia and Australo-Papua provide an ideal study system. We recover the phylogenetic relationships of the region's montane species and reconstruct their ancestral geographic ranges, elevational ranges, and migratory behavior. We also perform genomic population studies of three super-dispersive montane species/clades with broad island distributions. Eurasian-origin species populated archipelagos via direct colonization between mountains. This mode of colonization appears related to ancestral adaptations to cold and seasonal climates, specifically short-distance migration. Australo-Papuan-origin mountain populations, by contrast, evolved from lowland ancestors, and highland distribution mostly precludes their further colonization of island mountains. Our study explains much of the distributional variation within a complex biological system, and provides a synthesis of two seemingly discordant hypotheses for montane community formation.

High-Throughput Sequencing and Unsupervised Analysis of Formyltetrahydrofolate Synthetase (FTHFS) Gene Amplicons to Estimate Acetogenic Community Structure.

The formyltetrahydrofolate synthetase (FTHFS) gene is a molecular marker of choice to study the diversity of acetogenic communities. However, current analyses are limited due to lack of a high-throughput sequencing approach for FTHFS gene amplicons and a dedicated bioinformatics pipeline for data analysis, including taxonomic annotation and visualization of the sequence data. In the present study, we combined the barcode approach for multiplexed sequencing with unsupervised data analysis to visualize acetogenic community structure. We used samples from a biogas digester to develop proof-of-principle for our combined approach. We successfully generated high-throughput sequence data for the partial FTHFS gene and performed unsupervised data analysis using the novel bioinformatics pipeline "AcetoScan" presented in this study, which resulted in taxonomically annotated OTUs, phylogenetic tree, abundance plots and diversity indices. The results demonstrated that high-throughput sequencing can be used to sequence the FTHFS amplicons from a pool of samples, while the analysis pipeline AcetoScan can be reliably used to process the raw sequence data and visualize acetogenic community structure. The method and analysis pipeline described in this paper can assist in the identification and quantification of known or potentially new acetogens. The AcetoScan pipeline is freely available at https://github.com/abhijeetsingh1704/AcetoScan.

Parallel Evolution of Bower-Building Behavior in Two Groups of Bowerbirds Suggested by Phylogenomics.

The bowerbirds in New Guinea and Australia include species that build the largest and perhaps most elaborately decorated constructions outside of humans. The males use these courtship bowers, along with their displays, to attract females. In these species, the mating system is polygynous and the females alone incubate and feed the nestlings. The bowerbirds also include 10 species of the socially monogamous catbirds in which the male participates in most aspects of raising the young. How the bower-building behavior evolved has remained poorly understood, as no comprehensive phylogeny exists for the family. It has been assumed that the monogamous catbird clade is sister to all polygynous species. We here test this hypothesis using a newly developed pipeline for obtaining homologous alignments of thousands of exonic and intronic regions from genomic data to build a phylogeny. Our well-supported species tree shows that the polygynous, bower-building species are not monophyletic. The result suggests either that bower-building behavior is an ancestral condition in the family that was secondarily lost in the catbirds, or that it has arisen in parallel in two lineages of bowerbirds. We favor the latter hypothesis based on an ancestral character reconstruction showing that polygyny but not bower-building is ancestral in bowerbirds, and on the observation that Scenopoeetes dentirostris, the sister species to one of the bower-building clades, does not build a proper bower but constructs a court for male display. This species is also sexually monomorphic in plumage despite having a polygynous mating system. We argue that the relatively stable tropical and subtropical forest environment in combination with low predator pressure and rich food access (mostly fruit) facilitated the evolution of these unique life-history traits. [Adaptive radiation; bowerbirds; mating system, sexual selection; whole genome sequencing.].

Identifying Bird Remains Using Ancient DNA Barcoding.

Bird remains that are difficult to identify taxonomically using morphological methods, are common in the palaeontological record. Other types of challenging avian material include artefacts and food items from endangered taxa, as well as remains from aircraft strikes. We here present a DNA-based method that enables taxonomic identification of bird remains, even from material where the DNA is heavily degraded. The method is based on the amplification and sequencing of two short variable parts of the 16S region in the mitochondrial genome. To demonstrate the applicability of this approach, we evaluated the method on a set of Holocene and Late Pleistocene postcranial bird bones from several palaeontological and archaeological sites in Europe with good success.

Metagenomic Analysis of the Indian Ocean Picocyanobacterial Community: Structure, Potential Function and Evolution.

Unicellular cyanobacteria are ubiquitous photoautotrophic microbes that contribute substantially to global primary production. Picocyanobacteria such as Synechococcus and Prochlorococcus depend on chlorophyll a-binding protein complexes to capture light energy. In addition, Synechococcus has accessory pigments organized into phycobilisomes, and Prochlorococcus contains chlorophyll b. Across a surface water transect spanning the sparsely studied tropical Indian Ocean, we examined Synechococcus and Prochlorococcus occurrence, taxonomy and habitat preference in an evolutionary context. Shotgun sequencing of size fractionated microbial communities from 0.1 µm to 20 µm and subsequent phylogenetic analysis indicated that cyanobacteria account for up to 15% of annotated reads, with the genera Prochlorococcus and Synechococcus comprising 90% of the cyanobacterial reads, even in the largest size fraction (3.0-20 mm). Phylogenetic analyses of cyanobacterial light-harvesting genes (chl-binding pcb/isiA, allophycocyanin (apcAB), phycocyanin (cpcAB) and phycoerythin (cpeAB)) mostly identified picocyanobacteria clades comprised of overlapping sequences obtained from Indian Ocean, Atlantic and/or Pacific Oceans samples. Habitat reconstructions coupled with phylogenetic analysis of the Indian Ocean samples suggested that large Synechococcus-like ancestors in coastal waters expanded their ecological niche towards open oligotrophic waters in the Indian Ocean through lineage diversification and associated streamlining of genomes (e.g. loss of phycobilisomes and acquisition of Chl b); resulting in contemporary small celled Prochlorococcus. Comparative metagenomic analysis with picocyanobacteria populations in other oceans suggests that this evolutionary scenario may be globally important.

Metagenomic Analysis from the Interior of a Speleothem in Tjuv-Ante's Cave, Northern Sweden.

Speleothems are secondary mineral deposits normally formed by water supersaturated with calcium carbonate percolating into underground caves, and are often associated with low-nutrient and mostly non-phototrophic conditions. Tjuv-Ante's cave is a shallow-depth cave formed by the action of waves, with granite and dolerite as major components, and opal-A and calcite as part of the speleothems, making it a rare kind of cave. We generated two DNA shotgun sequencing metagenomic datasets from the interior of a speleothem from Tjuv-Ante's cave representing areas of old and relatively recent speleothem formation. We used these datasets to perform i) an evaluation of the use of these speleothems as past biodiversity archives, ii) functional and taxonomic profiling of the speleothem's different formation periods, and iii) taxonomic comparison of the metagenomic results to previous microscopic analyses from a nearby speleothem of the same cave. Our analyses confirm the abundance of Actinobacteria and fungi as previously reported by microscopic analyses on this cave, however we also discovered a larger biodiversity. Interestingly, we identified photosynthetic genes, as well as genes related to iron and sulphur metabolism, suggesting the presence of chemoautotrophs. Furthermore, we identified taxa and functions related to biomineralization. However, we could not confidently establish the use of this type of speleothems as biological paleoarchives due to the potential leaching from the outside of the cave and the DNA damage that we propose has been caused by the fungal chemical etching.

Phylogeny, evolution and classification of gall wasps: the plot thickens.

Gall wasps (Cynipidae) represent the most spectacular radiation of gall-inducing insects. In addition to true gall formers, gall wasps also include phytophagous inquilines, which live inside the galls induced by gall wasps or other insects. Here we present the first comprehensive molecular and total-evidence analyses of higher-level gall wasp relationships. We studied more than 100 taxa representing a rich selection of outgroups and the majority of described cynipid genera outside the diverse oak gall wasps (Cynipini), which were more sparsely sampled. About 5 kb of nucleotide data from one mitochondrial (COI) and four nuclear (28S, LWRh, EF1alpha F1, and EF1alpha F2) markers were analyzed separately and in combination with morphological and life-history data. According to previous morphology-based studies, gall wasps evolved in the Northern Hemisphere and were initially herb gallers. Inquilines originated once from gall inducers that lost the ability to initiate galls. Our results, albeit not conclusive, suggest a different scenario. The first gall wasps were more likely associated with woody host plants, and there must have been multiple origins of gall inducers, inquilines or both. One possibility is that gall inducers arose independently from inquilines in several lineages. Except for these surprising results, our analyses are largely consistent with previous studies. They confirm that gall wasps are conservative in their host-plant preferences, and that herb-galling lineages have radiated repeatedly onto the same set of unrelated host plants. We propose a revised classification of the family into twelve tribes, which are strongly supported as monophyletic across independent datasets. Four are new: Aulacideini, Phanacidini, Diastrophini and Ceroptresini. We present a key to the tribes and discuss their morphological and biological diversity. Until the relationships among the tribes are resolved, the origin and early evolution of gall wasps will remain elusive.

Local hopping mobile DNA implicated in pseudogene formation and reductive evolution in an obligate cyanobacteria-plant symbiosis.

BACKGROUND: Insertion sequences (ISs) are approximately 1 kbp long "jumping" genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I ("cut-and-paste") Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium 'Nostoc azollae' 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of 'Nostoc azollae' 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria. RESULTS: Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location ("local hopping") and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of 'Nostoc azollae' 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style. CONCLUSIONS: The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than genome wide this likely explains the proximity between IS elements and pseudogenes. These findings suggest that ISs facilitate the reductive evolution in for instance in the symbiotic cyanobacterium 'Nostoc azollae' 0708 and in other obligate prokaryotic symbionts.

Functional tradeoffs underpin salinity-driven divergence in microbial community composition.

Bacterial community composition and functional potential change subtly across gradients in the surface ocean. In contrast, while there are significant phylogenetic divergences between communities from freshwater and marine habitats, the underlying mechanisms to this phylogenetic structuring yet remain unknown. We hypothesized that the functional potential of natural bacterial communities is linked to this striking divide between microbiomes. To test this hypothesis, metagenomic sequencing of microbial communities along a 1,800 km transect in the Baltic Sea area, encompassing a continuous natural salinity gradient from limnic to fully marine conditions, was explored. Multivariate statistical analyses showed that salinity is the main determinant of dramatic changes in microbial community composition, but also of large scale changes in core metabolic functions of bacteria. Strikingly, genetically and metabolically different pathways for key metabolic processes, such as respiration, biosynthesis of quinones and isoprenoids, glycolysis and osmolyte transport, were differentially abundant at high and low salinities. These shifts in functional capacities were observed at multiple taxonomic levels and within dominant bacterial phyla, while bacteria, such as SAR11, were able to adapt to the entire salinity gradient. We propose that the large differences in central metabolism required at high and low salinities dictate the striking divide between freshwater and marine microbiomes, and that the ability to inhabit different salinity regimes evolved early during bacterial phylogenetic differentiation. These findings significantly advance our understanding of microbial distributions and stress the need to incorporate salinity in future climate change models that predict increased levels of precipitation and a reduction in salinity.

Dating the diversification of the major lineages of Passeriformes (Aves).

BACKGROUND: The avian Order Passeriformes is an enormously species-rich group, which comprises almost 60% of all living bird species. This diverse order is believed to have originated before the break-up of Gondwana in the late Cretaceous. However, previous molecular dating studies have relied heavily on the geological split between New Zealand and Antarctica, assumed to have occurred 85-82 Mya, for calibrating the molecular clock and might thus be circular in their argument. RESULTS: This study provides a time-scale for the evolution of the major clades of passerines using seven nuclear markers, five taxonomically well-determined passerine fossils, and an updated interpretation of the New Zealand split from Antarctica 85-52 Mya in a Bayesian relaxed-clock approach. We also assess how different interpretations of the New Zealand-Antarctica vicariance event influence our age estimates. Our results suggest that the diversification of Passeriformes began in the late Cretaceous or early Cenozoic. Removing the root calibration for the New Zealand-Antarctica vicariance event (85-52 Mya) dramatically increases the 95% credibility intervals and leads to unrealistically old age estimates. We assess the individual characteristics of the seven nuclear genes analyzed in our study. Our analyses provide estimates of divergence times for the major groups of passerines, which can be used as secondary calibration points in future molecular studies. CONCLUSIONS: Our analysis takes recent paleontological and geological findings into account and provides the best estimate of the passerine evolutionary time-scale currently available. This time-scale provides a temporal framework for further biogeographical, ecological, and co-evolutionary studies of the largest bird radiation, and adds to the growing support for a Cretaceous origin of Passeriformes.

Genome erosion in a nitrogen-fixing vertically transmitted endosymbiotic multicellular cyanobacterium.

BACKGROUND: An ancient cyanobacterial incorporation into a eukaryotic organism led to the evolution of plastids (chloroplasts) and subsequently to the origin of the plant kingdom. The underlying mechanism and the identities of the partners in this monophyletic event remain elusive. METHODOLOGY/PRINCIPAL FINDINGS: To shed light on this evolutionary process, we sequenced the genome of a cyanobacterium residing extracellularly in an endosymbiosis with a plant, the water-fern Azolla filiculoides Lam. This symbiosis was selected as it has characters which make it unique among extant cyanobacterial plant symbioses: the cyanobacterium lacks autonomous growth and is vertically transmitted between plant generations. Our results reveal features of evolutionary significance. The genome is in an eroding state, evidenced by a large proportion of pseudogenes (31.2%) and a high frequency of transposable elements (approximately 600) scattered throughout the genome. Pseudogenization is found in genes such as the replication initiator dnaA and DNA repair genes, considered essential to free-living cyanobacteria. For some functional categories of genes pseudogenes are more prevalent than functional genes. Loss of function is apparent even within the 'core' gene categories of bacteria, such as genes involved in glycolysis and nutrient uptake. In contrast, serving as a critical source of nitrogen for the host, genes related to metabolic processes such as cell differentiation and nitrogen-fixation are well preserved. CONCLUSIONS/SIGNIFICANCE: This is the first finding of genome degradation in a plant symbiont and phenotypically complex cyanobacterium and one of only a few extracellular endosymbionts described showing signs of reductive genome evolution. Our findings suggest an ongoing selective streamlining of this cyanobacterial genome which has resulted in an organism devoted to nitrogen fixation and devoid of autonomous growth. The cyanobacterial symbiont of Azolla can thus be considered at the initial phase of a transition from free-living organism to a nitrogen-fixing plant entity, a transition process which may mimic what drove the evolution of chloroplasts from a cyanobacterial ancestor.

Plastid and nuclear DNA markers reveal intricate relationships at subfamilial and tribal levels in the soapberry family (Sapindaceae).

The economically important soapberry family (Sapindaceae) comprises about 1900 species mainly found in the tropical regions of the world, with only a few genera being restricted to temperate areas. The infrafamilial classification of the Sapindaceae and its relationships to the closely related Aceraceae and Hippocastanaceae - which have now been included in an expanded definition of Sapindaceae (i.e., subfamily Hippocastanoideae) - have been debated for decades. Here we present a phylogenetic analysis of Sapindaceae based on eight DNA sequence regions from the plastid and nuclear genomes and including 85 of the 141 genera defined within the family. Our study comprises 997 new sequences of Sapindaceae from 152 specimens. Despite presenting 18.6% of missing data our complete data set produced a topology fully congruent with the one obtained from a subset without missing data, but including fewer markers. The use of additional information therefore led to a consistent result in the relative position of clades and allowed the definition of a new phylogenetic hypothesis. Our results confirm a high level of paraphyly and polyphyly at the subfamilial and tribal levels and even contest the monophyletic status of several genera. Our study confirms that the Chinese monotypic genus Xanthoceras is sister to the rest of the family, in which subfamily Hippocastanoideae is sister to a clade comprising subfamilies Dodonaeoideae and Sapindoideae. On the basis of the strong support demonstrated in Sapindoideae, Dodonaeoideae and Hippocastanoideae as well as in 14 subclades, we propose and discuss informal groupings as basis for a new classification of Sapindaceae.

Tracing the impact of the Andean uplift on Neotropical plant evolution.

Recent phylogenetic studies have revealed the major role played by the uplift of the Andes in the extraordinary diversification of the Neotropical flora. These studies, however, have typically considered the Andean uplift as a single, time-limited event fostering the evolution of highland elements. This contrasts with geological reconstructions indicating that the uplift occurred in discrete periods from west to east and that it affected different regions at different times. We introduce an approach for integrating Andean tectonics with biogeographic reconstructions of Neotropical plants, using the coffee family (Rubiaceae) as a model group. The distribution of this family spans highland and montane habitats as well as tropical lowlands of Central and South America, thus offering a unique opportunity to study the influence of the Andean uplift on the entire Neotropical flora. Our results suggest that the Rubiaceae originated in the Paleotropics and used the boreotropical connection to reach South America. The biogeographic patterns found corroborate the existence of a long-lasting dispersal barrier between the Northern and Central Andes, the "Western Andean Portal." The uplift of the Eastern Cordillera ended this barrier, allowing dispersal of boreotropical lineages to the South, but gave rise to a huge wetland system ("Lake Pebas") in western Amazonia that prevented in situ speciation and floristic dispersal between the Andes and Amazonia for at least 6 million years. Here, we provide evidence of these events in plants.

Allopolyploidy in Fragariinae (Rosaceae): comparing four DNA sequence regions, with comments on classification.

Potential events of allopolyploidy may be indicated by incongruences between separate phylogenies based on plastid and nuclear gene sequences. We sequenced two plastid regions and two nuclear ribosomal regions for 34 ingroup taxa in Fragariinae (Rosaceae), and six outgroup taxa. We found five well supported incongruences that might indicate allopolyploidy events. The incongruences involved Aphanes arvensis, Potentilla miyabei, Potentilla cuneata, Fragaria vesca/moschata, and the Drymocallis clade. We evaluated the strength of conflict and conclude that allopolyploidy may be hypothesised in the four first cases. Phylogenies were estimated using Bayesian inference and analyses were evaluated using convergence diagnostics. Taxonomic implications are discussed for genera such as Alchemilla, Sibbaldianthe, Chamaerhodos, Drymocallis and Fragaria, and for the monospecific Sibbaldiopsis and Potaninia that are nested inside other genera. Two orphan Potentilla species, P. miyabei and P. cuneata are placed in Fragariinae. However, due to unresolved topological incongruences they are not reclassified in any genus.

Accounting for phylogenetic uncertainty in biogeography: a Bayesian approach to dispersal-vicariance analysis of the thrushes (Aves: Turdus).

The phylogeny of the thrushes (Aves: Turdus) has been difficult to reconstruct due to short internal branches and lack of node support for certain parts of the tree. Reconstructing the biogeographic history of this group is further complicated by the fact that current implementations of biogeographic methods, such as dispersal-vicariance analysis (DIVA; Ronquist, 1997), require a fully resolved tree. Here, we apply a Bayesian approach to dispersal-vicariance analysis that accounts for phylogenetic uncertainty and allows a more accurate analysis of the biogeographic history of lineages. Specifically, ancestral area reconstructions can be presented as marginal distributions, thus displaying the underlying topological uncertainty. Moreover, if there are multiple optimal solutions for a single node on a certain tree, integrating over the posterior distribution of trees often reveals a preference for a narrower set of solutions. We find that despite the uncertainty in tree topology, ancestral area reconstructions indicate that the Turdus clade originated in the eastern Palearctic during the Late Miocene. This was followed by an early dispersal to Africa from where a worldwide radiation took place. The uncertainty in tree topology and short branch lengths seems to indicate that this radiation took place within a limited time span during the Late Pliocene. The results support the role of Africa as a probable source area for intercontinental dispersals as suggested for other passerine groups, including basal diversification within the songbird tree.

AWTY (are we there yet?): a system for graphical exploration of MCMC convergence in Bayesian phylogenetics.

UNLABELLED: A key element to a successful Markov chain Monte Carlo (MCMC) inference is the programming and run performance of the Markov chain. However, the explicit use of quality assessments of the MCMC simulations-convergence diagnostics-in phylogenetics is still uncommon. Here, we present a simple tool that uses the output from MCMC simulations and visualizes a number of properties of primary interest in a Bayesian phylogenetic analysis, such as convergence rates of posterior split probabilities and branch lengths. Graphical exploration of the output from phylogenetic MCMC simulations gives intuitive and often crucial information on the success and reliability of the analysis. The tool presented here complements convergence diagnostics already available in other software packages primarily designed for other applications of MCMC. Importantly, the common practice of using trace-plots of a single parameter or summary statistic, such as the likelihood score of sampled trees, can be misleading for assessing the success of a phylogenetic MCMC simulation. AVAILABILITY: The program is available as source under the GNU General Public License and as a web application at http://ceb.scs.fsu.edu/awty.

Phylogenetic relationships of typical antbirds (Thamnophilidae) and test of incongruence based on Bayes factors.

BACKGROUND: The typical antbirds (Thamnophilidae) form a monophyletic and diverse family of suboscine passerines that inhabit neotropical forests. However, the phylogenetic relationships within this assemblage are poorly understood. Herein, we present a hypothesis of the generic relationships of this group based on Bayesian inference analyses of two nuclear introns and the mitochondrial cytochrome b gene. The level of phylogenetic congruence between the individual genes has been investigated utilizing Bayes factors. We also explore how changes in the substitution models affected the observed incongruence between partitions of our data set. RESULTS: The phylogenetic analysis supports both novel relationships, as well as traditional groupings. Among the more interesting novel relationship suggested is that the Terenura antwrens, the wing-banded antbird (Myrmornis torquata), the spot-winged antshrike (Pygiptila stellaris) and the russet antshrike (Thamnistes anabatinus) are sisters to all other typical antbirds. The remaining genera fall into two major clades. The first includes antshrikes, antvireos and the Herpsilochmus antwrens, while the second clade consists of most antwren genera, the Myrmeciza antbirds, the "professional" ant-following antbirds, and allied species. Our results also support previously suggested polyphyly of Myrmotherula antwrens and Myrmeciza antbirds. The tests of phylogenetic incongruence, using Bayes factors, clearly suggests that allowing the gene partitions to have separate topology parameters clearly increased the model likelihood. However, changing a component of the nucleotide substitution model had much higher impact on the model likelihood. CONCLUSIONS: The phylogenetic results are in broad agreement with traditional classification of the typical antbirds, but some relationships are unexpected based on external morphology. In these cases their true affinities may have been obscured by convergent evolution and morphological adaptations to new habitats or food sources, and genera like Myrmeciza antbirds and the Myrmotherula antwrens obviously need taxonomic revisions. Although, Bayes factors seem promising for evaluating the relative contribution of components to an evolutionary model, the results suggests that even if strong evidence for a model allowing separate topology parameters is found, this might not mean strong evidence for separate gene phylogenies, as long as vital components of the substitution model are still missing.

Bayesian phylogenetic analysis of combined data.

The recent development of Bayesian phylogenetic inference using Markov chain Monte Carlo (MCMC) techniques has facilitated the exploration of parameter-rich evolutionary models. At the same time, stochastic models have become more realistic (and complex) and have been extended to new types of data, such as morphology. Based on this foundation, we developed a Bayesian MCMC approach to the analysis of combined data sets and explored its utility in inferring relationships among gall wasps based on data from morphology and four genes (nuclear and mitochondrial, ribosomal and protein coding). Examined models range in complexity from those recognizing only a morphological and a molecular partition to those having complex substitution models with independent parameters for each gene. Bayesian MCMC analysis deals efficiently with complex models: convergence occurs faster and more predictably for complex models, mixing is adequate for all parameters even under very complex models, and the parameter update cycle is virtually unaffected by model partitioning across sites. Morphology contributed only 5% of the characters in the data set but nevertheless influenced the combined-data tree, supporting the utility of morphological data in multigene analyses. We used Bayesian criteria (Bayes factors) to show that process heterogeneity across data partitions is a significant model component, although not as important as among-site rate variation. More complex evolutionary models are associated with more topological uncertainty and less conflict between morphology and molecules. Bayes factors sometimes favor simpler models over considerably more parameter-rich models, but the best model overall is also the most complex and Bayes factors do not support exclusion of apparently weak parameters from this model. Thus, Bayes factors appear to be useful for selecting among complex models, but it is still unclear whether their use strikes a reasonable balance between model complexity and error in parameter estimates.

A maximum-likelihood analysis of eight phylogenetic markers in gallwasps (Hymenoptera: Cynipidae): implications for insect phylogenetic studies.

We assessed the utility of eight DNA sequence markers (5.8S rDNA, 18S rDNA, 28S rDNA, ITS regions, long-wavelength opsin, elongation factor 1-alpha, cytochrome b, and cytochrome oxidase I) in reconstructing phylogenetic relationships at various levels of divergence in gallwasps (Hymenoptera: Cynipidae), using a set of eight exemplar taxa. We report sequence divergence values and saturation levels and compare phylogenetic results of these sequences analyzed both separately and combined to a well-corroborated morphological phylogeny. Likelihood ratio tests were used to find the best evolutionary model fitting each of the markers. The likelihood model best explaining the data is, for most loci, parameter rich, with strong A-T bias for mitochondrial loci and strong rate heterogeneity for the majority of loci. Our data suggest that 28S rDNA, elongation factor 1-alpha, and long-wavelength opsin may be potentially useful markers for the resolution of cynipid and other insect within-family-level divergences (circa 50-100 mya old), whereas mitochondrial loci and ITS regions are most useful for lower-level phylogenetics. In contrast, the 18S rDNA marker is likely to be useful for the resolution of above-family-level relationships.