Understanding mammalian evolution using Bayesian phylogenetic inference

Indel evolution of mammalian introns and the utility of non-coding nuclear markers in eutherian phylogenetics

Nuclear DNA intron sequences are increasingly used to investigate evolutionary relationships among closely related organisms. The phylogenetic usefulness of intron sequences at higher taxonomic levels has, however, not been firmly established and very few studies have used these markers to address evolutionary questions above the family level. In addition, the mechanisms driving intron evolution are not well understood. We compared DNA sequence data derived from three presumably independently segregating introns (THY, PRKC I and MGF) across 158 mammalian species. All currently recognized extant eutherian mammalian orders were included with the exception of Cingulata, Dermoptera and Scandentia. The total aligned length of the data was 6366 base pairs (bp); after the exclusion of autapomorphic insertions, 1511 bp were analyzed. In many instances the Bayesian and parsimony analyses were complementary and gave significant posterior probability and bootstrap support (>80) for the monophyly of Afrotheria, Euarchontoglires, Laurasiatheria and Boreoeutheria. Apart from finding congruent support when using these methods, the intron data also provided several indels longer than 3 bp that support, among others, the monophyly of Afrotheria, Paenungulata, Ferae and Boreoeutheria. A quantitative analysis of insertions and deletions suggested that there was a 75% bias towards deletions. The average insertion size in the mammalian data set was 16.49 bp +/- 57.70 while the average deletion was much smaller (4.47 bp +/- 14.17). The tendency towards large insertions and small deletions is highlighted by the observation that out of a total of 17 indels larger than 100 bp, 15 were insertions. The majority of indels (>60% of all events) were 1 or 2 bp changes. Although the average overall indel substitution rate of 0.00559 per site is comparable to that previously reported for rodents and primates, individual analyses among different evolutionary lineages provide evidence for differences in the formation rate of indels among the different mammalian groups.

In silico identification and Bayesian phylogenetic analysis of multiple new mammalian kallikrein gene families

Kallikrein gene families have been identified previously in genomes of the human, the mouse, and the rat, and individual kallikrein-like genes have been found in many more species. This study presents the in silico identification of kallikrein gene families in the recently sequenced genomes of four additional mammalian species, the chimpanzee, the dog, the pig, and the opossum. Phylogenies were constructed with gene sequences from all seven mammalian families, using Bayesian analysis, which clarified the evolutionary relationships between these genes. Individual gene sequences, as well as concatenated constructs of multiple sequences, were used. Fifteen kallikrein genes were located in the chimpanzee (Pan troglodytes) genome, while only 14 were identified in the canine (Canis familiaris) genome as no orthologue to human KLK3 was found. Thirteen genes were identified from the pig (Sus scrofa) genome, which lacked homologues to KLK2 and KLK3, and 11 genes, orthologous to human KLK5 through KLK15, were found in the opossum (Monodelphis domestica) genome. No kallikrein genes were identified from the available genome sequences of the chicken (Gallus gallus) or African clawed frog (Xenopus tropicalis). Within the family of kallikreins several subfamilies were suggested by phylogenetic analysis. One consisted of KLK4, KLK5, and KLK14; another of KLK9, KLK11, and KLK15; a third of KLK10 and KLK12; a fourth of KLK6 and KLK13; and finally one of KLK8 and the classical kallikreins (KLK1, KLK2, and KLK3).

The Phylogeny of Early Eureptiles: Comparing Parsimony and Bayesian Approaches in the Investigation of a Basal Fossil Clade

For the first time the phylogenetic relationships of early eureptiles, consisting of captorhinids, diapsids, and protorothyridids, are investigated in a modern phylogenetic context using both parsimony and Bayesian approaches. Ninety parsimony-informative characters and 25 taxa were included in the analyses. The Bayesian analysis was run with and without a gamma-shape parameter allowing for variable rates across characters. In addition, we ran two more Bayesian analyses that included 42 autapomorphies and thus parsimony-uninformative characters in order to test the effect of variable branch lengths. The different analyses largely converged to the same topology, suggesting that the "protorothyridid" Coelostegus is the sister taxon of all other eureptiles and that the remaining "protorothyridids" are paraphyletic. Also, there is a close relationship between diapsids and Anthracodromeus, Cephalerpeton, and Protorothyris, a grouping of Thuringothyris with captorhinids, and a variable position of the "protorothyridids" Brouffia, Hylonomus, and Paleothyris. The lack of resolution in some parts of the tree might be due to "hard polytomies" and short divergence times between the respective taxa. The tree topology is consistent with the hypothesis that the temporal fenestrations of diapsid reptiles appear to be the consequence of a more lightly built skeleton, indicating a significant ecological shift in the early stages of diapsid evolution. Bayesian analysis is a very useful additional approach in studies of fossil taxa in which more traditional statistical support like the bootstrap is often weak. However, the exclusive use of the Mk model appears suitable only if autapomorphic characters are included, whereas the Mk+gamma model performed well with or without autapomorphies.