Andrew F. Hugall

Calibration Choice, Rate Smoothing, and the Pattern of Tetrapod Diversification According to the Long Nuclear Gene RAG-1

A phylogeny of tetrapods is inferred from nearly complete sequences of the nuclear RAG-1 gene sampled across 88 taxa encompassing all major clades, analyzed via parsimony and Bayesian methods. The phylogeny provides support for Lissamphibia, Theria, Lepidosauria, a turtle-archosaur clade, as well as most traditionally accepted groupings. This tree allows simultaneous molecular clock dating for all tetrapod groups using a set of well-corroborated calibrations. Relaxed clock (PLRS) methods, using the amniote = 315 Mya (million years ago) calibration or a set of consistent calibrations, recovers reasonable divergence dates for most groups. However, the analysis systematically underestimates divergence dates within archosaurs. The bird-crocodile split, robustly documented in the fossil record as being around approximately 245 Mya, is estimated at only approximately 190 Mya, and dates for other divergences within archosaurs are similarly underestimated. Archosaurs, and particulary turtles have slow apparent rates possibly confounding rate modeling, and inclusion of calibrations within archosaurs (despite their high deviances) not only improves divergence estimates within archosaurs, but also across other groups. Notably, the monotreme-therian split ( approximately 210 Mya) matches the fossil record; the squamate radiation ( approximately 190 Mya) is younger than suggested by some recent molecular studies and inconsistent with identification of approximately 220 and approximately 165 Myo (million-year-old) fossils as acrodont iguanians and approximately 95 Myo fossils colubroid snakes; the bird-lizard (reptile) split is considerably older than fossil estimates (< or = 285 Mya); and Sphenodon is a remarkable phylogenetic relic, being the sole survivor of a lineage more than a quarter of a billion years old. Comparison with other molecular clock studies of tetrapod divergences suggests that the common practice of enforcing most calibrations as minima, with a single liberal maximal constraint, will systematically overestimate divergence dates. Similarly, saturation of mitochondrial DNA sequences, and the resultant greater compression of basal branches means that using only external deep calibrations will also lead to inflated age estimates within the focal ingroup.

Molecular Support for Vicariance as a Source of Diversity in Rainforest

The origin of high diversity in tropical rainforests is attributed to divergence amongst relatively mesic, late Pleistocene refuges. To test this hypothesis we analysed sequence variation within the mitochondrial DNA of populations of six rainforest-dwelling vertebrate species (one lizard and five birds) in the Wet Tropical rainforests of northeastern Australia. Vicariance among historical refuges was indicated by geographically congruent variation across a historical climatic barrier previously predicted by modelling. Sequence divergence across the barrier varied widely, being highest in species with lowest vagility and greatest restriction to rainforest. A high altitude, rainforest-restricted species was exceptional in lacking mtDNA variation. These data support the suggested role of vicariance in promoting evolutionary diversity in rainforests but also indicate variation in the timing of vicariance events, with most speciation or divergence events occurring well before the late Pleistocene. The relation between the amount of molecular divergence and current ecology suggests a sequence of isolation events, perhaps of increasing severity, occurring between the late Tertiary and the late Pleistocene.

Accelerated speciation in colour-polymorphic birds.

Colour polymorphism exemplifies extreme morphological diversity within populations. It is taxonomically widespread but generally rare. Theory suggests that where colour polymorphism does occur, processes generating and maintaining it can promote speciation but the generality of this claim is unclear. Here we confirm, using species-level molecular phylogenies for five families of non-passerine birds, that colour polymorphism is associated with accelerated speciation rates in the three groups in which polymorphism is most prevalent. In all five groups, colour polymorphism is lost at a significantly greater rate than it is gained. Thus, the general rarity and phylogenetic dispersion of colour polymorphism is accounted for by a combination of higher speciation rate and higher transition rate from polymorphism to monomorphism, consistent with theoretical models where speciation is driven by fixation of one or more morphs. This is corroborated by evidence from a species-level molecular phylogeny of passerines, incorporating 4,128 (66.5%) extant species, that polymorphic species tend to be younger than monomorphic species. Our results provide empirical support for the general proposition, dating from classical evolutionary theory, that colour polymorphism can increase speciation rates.

PHYLOGENETIC ANALYSIS OF THE MONOGENEA AND THEIR RELATIONSHIPS WITH DIGENEA AND EUCESTODA INFERRED FROM 28S RDNA SEQUENCES

Platyhelminth phylogeny is controversial. Phylogenetic analyses of the partial domain C1 and the full domains D1 and C2 (358 nucleotides) from the 28S ribosomal RNA gene for 21 species from the Monogenea, Digenea, Cestoda, and, as the outgroup, Tricladida reveal major departures from prevailing theory. The Digenea and not the Monogenea (Monopisthocotylea and Polyopisthocotylea) form the sister group of the cestodes; the Monopisthocotylea and Polyopisthocotylea are each monophyletic, but the Monogenea do not form a monophylum; the sister group of the Digenea + Cestoda is the Polyopisthocotylea; and Monopisthocotylea are the sister group of all other parasitic flatworms.

Phylogenomic Resolution of the Class Ophiuroidea Unlocks a Global Microfossil Record

Our understanding of the origin, evolution, and biogeography of seafloor fauna is limited because we have insufficient spatial and temporal data to resolve underlying processes. The abundance and wide distribution of modern and disarticulated fossil Ophiuroidea, including brittle stars and basket stars, make them an ideal model system for global marine biogeography if we have the phylogenetic framework necessary to link extant and fossil morphology in an evolutionary context. Here we construct a phylogeny from a highly complete 425-gene, 61-taxa transcriptome-based data set covering 15 of the 18 ophiuroid families and representatives of all extant echinoderm classes. We calibrate our phylogeny with a series of novel fossil discoveries from the early Mesozoic. We confirm the traditional paleontological view that ophiuroids are sister to the asteroids and date the crown group Ophiuroidea to the mid-Permian (270 ± 30 mega-annum). We refute all historical classification schemes of the Ophiuroidea based on gross structural characters but find strong congruence with schemes based on lateral arm plate microstructure and the temporal appearance of various plate morphologies in the fossil record. The verification that these microfossils contain phylogenetically informative characters unlocks their potential to advance our understanding of marine biogeographical processes.

Phylogeny of snakes (Serpentes): combining morphological and molecular data in likelihood, Bayesian and parsimony analyses

Abstract The phylogeny of living and fossil snakes is assessed using likelihood and parsimony approaches and a dataset combining 263 morphological characters with mitochondrial (2693 bp) and nuclear (1092 bp) gene sequences. The ‘no common mechanism’ (NCMr) and ‘Markovian’ (Mkv) models were employed for the morphological partition in likelihood analyses; likelihood scores in the NCMr model were more closely correlated with parsimony tree lengths. Both models accorded relatively less weight to the molecular data than did parsimony, with the effect being milder in the NCMr model. Partitioned branch and likelihood support values indicate that the mtDNA and nuclear gene partitions agree more closely with each other than with morphology. Despite differences between data partitions in phylogenetic signal, analytic models, and relative weighting, the parsimony and likelihood analyses all retrieved the following widely accepted groups: scolecophidians, alethinophidians, cylindrophiines, macrostomatans (sensu lato) and caenophidians. Anilius alone emerged as the most basal alethinophidian; the combined analyses resulted in a novel and stable position of uropeltines and cylindrophiines as the second‐most basal clade of alethinophidians. The limbed marine pachyophiids, along with Dinilysia and Wonambi, were always basal to all living snakes. Other results stable in all combined analyses include: Xenopeltis and Loxocemus were sister taxa (fide morphology) but clustered with pythonines (fide molecules), and Ungaliophis clustered with a boine‐erycine clade (fide molecules). Tropidophis remains enigmatic; it emerges as a basal alethinophidian in the parsimony analyses (fide molecules) but a derived form in the likelihood analyses (fide morphology), largely due to the different relative weighting accorded to data partitions.

THE LIKELIHOOD NODE DENSITY EFFECT AND CONSEQUENCES FOR EVOLUTIONARY STUDIES OF MOLECULAR RATES

Abstract Many molecular phylogenies show longer root-to-tip path lengths in species-rich groups, encouraging hypotheses linking cladogenesis with accelerated molecular evolution. However, the pattern can also be caused by an artifact called the node density effect (NDE): this effect occurs when the method used to reconstruct a tree underestimates multiple hits that would have been revealed by extra nodes, leading to longer root-to-tip path lengths in clades with more terminal taxa. Here we use a twofold approach to demonstrate that maximum likelihood and Bayesian methods also suffer from the NDE known to affect parsimony. First, simulations deliberately mismatching the simulation and reconstruction models show that the greater the model disparity, the greater the gap between actual and reconstructed tree lengths, and the greater the NDE. Second, taxon sampling manipulation with empirical data shows that NDE can still be present when using optimized models: across 12 datasets, 70 out of 109 sister path comparisons showed significant evidence of NDE. Unless the model fairly accurately reconstructs the real tree length—and given the complexity of real sequence evolution this may be uncommon—it will consistently produce a node density artifact. At commonly encountered divergence levels, a 10% underestimation of tree length results in ≥ 80% of simulated phylogenies showing a positive NDE. Bayesian trees have a slight but consistently stronger effect. This pervasive methodological artifact increases apparent rate heterogeneity, and can compromise investigations of factors influencing molecular evolutionary rate that use path lengths in topologically asymmetric trees.

SEXUAL SELECTION AND THE EVOLUTION OF COMPLEX COLOR PATTERNS IN DRAGON LIZARDS

Many species have elaborate and complex coloration and patterning, which often differ between the sexes. Sexual selection may increase the size or intensity of color patches (elaboration) in one sex or drive the evolution of novel signal elements (innovation). The latter potentially increases color pattern complexity. Color pattern complexity may also be influenced by ecological factors related to predation and environment; however, very few studies have investigated the effects of both sexual and natural selection on color pattern complexity across species. We used a phylogenetic comparative approach to examine these effects in 85 species and subspecies of Australian dragon lizards (family Agamidae). We quantified color pattern complexity by adapting the Shannon–Wiener diversity index. There were clear sex differences in color pattern complexity, which were positively correlated with both sexual dichromatism and sexual size dimorphism, consistent with the idea that sexual selection plays a significant role in the evolution of color pattern complexity. By contrast, we found little evidence of a link between environmental factors and color pattern complexity on body regions exposed to predators. Our results suggest that sexual selection rather than natural selection has led to increased color pattern complexity in males.

An Exon-Capture System for the Entire Class Ophiuroidea

Exon-capture studies have typically been restricted to relatively shallow phylogenetic scales due primarily to hybridization constraints. Here, we present an exon-capture system for an entire class of marine invertebrates, the Ophiuroidea, built upon a phylogenetically diverse transcriptome foundation. The system captures approximately 90% of the 1,552 exon target, across all major lineages of the quarter-billion-year-old extant crown group. Key features of our system are 1) basing the target on an alignment of orthologous genes determined from 52 transcriptomes spanning the phylogenetic diversity and trimmed to remove anything difficult to capture, map, or align; 2) use of multiple artificial representatives based on ancestral state reconstructions rather than exemplars to improve capture and mapping of the target; 3) mapping reads to a multi-reference alignment; and 4) using patterns of site polymorphism to distinguish among paralogy, polyploidy, allelic differences, and sample contamination. The resulting data give a well-resolved tree (currently standing at 417 samples, 275,352 sites, 91% data-complete) that will transform our understanding of ophiuroid evolution and biogeography.

A molecular phylogeny of rainbow skinks (Scincidae: Carlia): taxonomic and biogeographic implications

The phylogenetic relationships amongst 29 species of Carlia and Lygisaurus were estimated using a 726-base-pair segment of the protein-coding mitochondrial ND4 gene. Results do not support the recent resurrection of the genus Lygisaurus. Although most Lygisaurus species formed a single clade, this clade is nested within Carlia and includes Carlia parrhasius. Due to this new molecular evidence, and the paucity of diagnostic morphological characters separating the genera, Lygisaurus de Vis 1884 is re-synonymised with Carlia Gray 1845. Our analysis is also inconsistent with a previous suggestion that Lygisaurus timlowi should be removed to Menetia, a genus that is distantly related relative to outgroups used here. Intraspecific variation in Carlia is, in several instances, greater than interspecific distance. The most strikingly divergent lineages are found within C. rubrigularis, which appears to be paraphyletic, with southern populations more closely related to C. rhomboidalis than to northern populations of C. rubrigularis. The two C. rubrigularis-C. rhomboidalis lineages form part of a major polytomy at an intermediate level of divergence. Lack of resolution at this level, however, does not appear to be due to saturation or loss of phylogenetic signal. Rather, the polytomy probably reflects a period of relatively rapid diversification that occurred sometime during the Miocene.