Clifford W. Cunningham

Combining data in phylogenetic analysis.

Systematists have access to multiple sources of character information in phylogenetic analysis. For example, it is not unusual to have nucleotide sequences from several different genes, or to have molecular and morphological data. How should diverse data be analyzed in phylogenetic analysis? Several methods have been proposed for the treatment of partitioned data: the total evidence, separate analysis, and conditional combination approaches. Here, we review some of the advantages and disadvantages of the different approaches, with special concentration on which methods help us to discern the evolutionary process and provide the most accurate estimates of phylogeny.

Arthropod relationships revealed by phylogenomic analysis of nuclear protein-coding sequences

The remarkable antiquity, diversity and ecological significance of arthropods have inspired numerous attempts to resolve their deep phylogenetic history, but the results of two decades of intensive molecular phylogenetics have been mixed. The discovery that terrestrial insects (Hexapoda) are more closely related to aquatic Crustacea than to the terrestrial centipedes and millipedes (Myriapoda) was an early, if exceptional, success. More typically, analyses based on limited samples of taxa and genes have generated results that are inconsistent, weakly supported and highly sensitive to analytical conditions. Here we present strongly supported results from likelihood, Bayesian and parsimony analyses of over 41 kilobases of aligned DNA sequence from 62 single-copy nuclear protein-coding genes from 75 arthropod species. These species represent every major arthropod lineage, plus five species of tardigrades and onychophorans as outgroups. Our results strongly support Pancrustacea (Hexapoda plus Crustacea) but also strongly favour the traditional morphology-based Mandibulata (Myriapoda plus Pancrustacea) over the molecule-based Paradoxopoda (Myriapoda plus Chelicerata). In addition to Hexapoda, Pancrustacea includes three major extant lineages of ‘crustaceans’, each spanning a significant range of morphological disparity. These are Oligostraca (ostracods, mystacocarids, branchiurans and pentastomids), Vericrustacea (malacostracans, thecostracans, copepods and branchiopods) and Xenocarida (cephalocarids and remipedes). Finally, within Pancrustacea we identify Xenocarida as the long-sought sister group to the Hexapoda, a result confirming that ‘crustaceans’ are not monophyletic. These results provide a statistically well-supported phylogenetic framework for the largest animal phylum and represent a step towards ending the often-heated, century-long debate on arthropod relationships.

Reconstructing ancestral character states: a critical reappraisal

Using parsimony to reconstruct ancestral character states on a phylogenetic tree has become a popular method for testing ecological and evolutionary hypotheses. Despite its popularity, the assumptions and uncertainties of reconstructing the ancestral states of a single character have received less attention than the much less challenging endeavor of reconstructing phylogenetic trees from many characters. Recent research suggests that parsimony reconstructions are often sensitive to violations of the almost universal assumption of equal probabilities of gains and losses. In addition, maximum likelihood has been developed as an alternative to parsimony reconstruction, and has also revealed a surprising amount of uncertainty in ancestral reconstructions.

PHYLOGEOGRAPHY AND HISTORICAL ECOLOGY OF THE NORTH ATLANTIC INTERTIDAL

Abstract Recent glaciation covered the full extent of rocky intertidal habitat along the coasts of New England and the Canadian Maritimes. To test whether this glaciation in fact caused wholesale extinction of obligate rocky intertidal invertebrates, and thus required a recolonization from Europe, we compared American and European populations using allelic diversity and techniques adapted from coalescent theory. Mitochondrial DNA sequences were collected from amphi-Atlantic populations of three cold-temperate obligate rocky intertidal species (a barnacle, Semibalanus balanoides, and two gastropods, Nucella lapillus and Littorina obtusata) and three cold-temperate habitat generalist species (a seastar, Asterias rubens; a mussel, Mytilus edulis, and an isopod, Idotea balthica). For many of these species we were able to estimate the lineage-specific mutation rate based on trans-Arctic divergences between Pacific and Atlantic taxa. These data indicate that some obligate rocky intertidal taxa have colonized New England from European populations. However, the patterns of persistence in North America indicate that other life-history traits, including mechanisms of dispersal, may be more important for surviving dramatic environmental and climatic change. Corresponding Editor: R. Burton

DIVERSITY IN THE WEAPONS OF SEXUAL SELECTION: HORN EVOLUTION IN THE BEETLE GENUS ONTHOPHAGUS (COLEOPTERA: SCARABAEIDAE)

Abstract Both ornaments and weapons of sexual selection frequently exhibit prolific interspecific diversity of form. Yet, most studies of this diversity have focused on ornaments involved with female mate choice, rather than on the weapons of male competition. With few exceptions, the mechanisms of divergence in weapon morphology remain largely unexplored. Here, we characterize the evolutionary radiation of one type of weapon: beetle horns. We use partial sequences from four nuclear and three mitochondrial genes to develop a phylogenetic hypothesis for a worldwide sample of 48 species from the dung beetle genus Onthophagus (Coleoptera: Scarabaeidae). We then use these data to test for multiple evolutionary origins of horns and to characterize the evolutionary radiation of horns. Although our limited sampling of one of the worlds most species-rich genera almost certainly underestimates the number of evolutionary events, our phylogeny reveals prolific evolutionary lability of these exaggerated sexually selected weapons (more than 25 separate gains and losses of five different horn types). We discuss these results in the context of the natural history of these beetles and explore ways that sexual selection and ecology may have interacted to generate this extraordinary diversity of weapon morphology.

A comparative study of asymmetric migration events across a marine biogeographic boundary.

In many nonclonal, benthic marine species, geographic distribution is mediated by the dispersal of their larvae. The dispersal and recruitment of marine larvae may be limited by temperature gradients that can affect mortality or by ocean currents that can directly affect the movements of pelagic larvae. We focus on Point Conception, a well‐known biogeographic boundary between the Californian and Oregonian biogeographic provinces, to investigate whether ocean currents affect patterns of gene flow in intertidal marine invertebrates. The predominance of pelagically dispersing species with northern range limits at Point Conception suggests that ocean currents can affect species distributions by erecting barriers to the dispersal of planktonic larvae. In this paper, we investigate whether the predominantly southward currents have left a recognizable genetic signature in species with pelagically dispersing larvae whose ranges span Point Conception. We use patterns of genetic diversity and a new method for inferring cladistic migration events to test the hypothesis that southward currents increase southward gene flow for species with pelagically dispersing larvae. We collected mitochondrial DNA (mtDNA) sequence data for the barnacles Balanus glandula and Chthamalus fissus and also reanalyzed a previously published mtDNA dataset (Strongylocentrotus purpuratus, Edmands et al. 1996). For all three species, our cladistic approach identified an excess of southward migration events across Point Conception. In data from a fourth species with nondispersing larvae (Nucella emarginata, Marko 1998), our method suggests that ocean currents have not played a role in generating genetic structure.

Local adaptation and species segregation in two mussel ( Mytilus edulis × × × Mytilus trossulus ) hybrid zones

Few marine hybrid zones have been studied extensively, the major exception being the hybrid zone between the mussels Mytilus edulis and Mytilus galloprovincialis in southwestern Europe. Here, we focus on two less studied hybrid zones that also involve Mytilus spp.; Mytilus edulis and Mytilus trossulus are sympatric and hybridize on both western and eastern coasts of the Atlantic Ocean. We review the dynamics of hybridization in these two hybrid zones and evaluate the role of local adaptation for maintaining species boundaries. In Scandinavia, hybridization and gene introgression is so extensive that no individuals with pure M. trossulus genotypes have been found. However, M. trossulus alleles are maintained at high frequencies in the extremely low salinity Baltic Sea for some allozyme genes. A synthesis of reciprocal transplantation experiments between different salinity regimes shows that unlinked Gpi and Pgm alleles change frequency following transplantation, such that post‐transplantation allelic composition resembles native populations found in the same salinity. These experiments provide strong evidence for salinity adaptation at Gpi and Pgm (or genes linked to them). In the Canadian Maritimes, pure M. edulis and M. trossulus individuals are abundant, and limited data suggest that M. edulis predominates in low salinity and sheltered conditions, whereas M. trossulus are more abundant on the wave‐exposed open coasts. We suggest that these conflicting patterns of species segregation are, in part, caused by local adaptation of Scandinavian M. trossulus to the extremely low salinity Baltic Sea environment.

Crab shell-crushing predation and gastropod architectural defense

Abstract The shell-breaking behavior of the crabs Ozius verreauxii Saussure 1853 and Eriphia squamata , Stimpson 1859 from the Bay of Panama is described. The master claws of both these crabs are well designed for breaking shells. Small shells, relative to the size of a crab predator, are crushed by progressively breaking off larger segments of a shells apex, while larger shells are peeled by inserting a large dactyl molar into the aperture of a shell and progressively chipping away the lip of the shell. Heavy gastropod shells are shown to be less vulnerable to crab predators than lighter shells, and narrow shell apertures and axial shell sculpture are demonstrated to be architectural features that deter crab predation. The incidence of architectural features which deter crab predation appears to be higher for smaller gastropod species than for larger gastropods which are too large for most crab predators. Large fish predators prey upon both gastropods and shell-crushing crabs. To avoid fish predators, both these prey groups seek refuge under rocks when covered by the tide. Fish predation thus appears to enforce a close sympatry between smaller gastropods and their crab predators.

BEST-FIT MAXIMUM-LIKELIHOOD MODELS FOR PHYLOGENETIC INFERENCE : EMPIRICAL TESTS WITH KNOWN PHYLOGENIES

Despite the proliferation of increasingly sophisticated models of DNA sequence evolution, choosing among models remains a major problem in phylogenetic reconstruction. The choice of appropriate models is thought to be especially important when there is large variation among branch lengths. We evaluated the ability of nested models to reconstruct experimentally generated, known phylogenies of bacteriophage T7 as we varied the terminal branch lengths. Then, for each phylogeny we determined the best‐fit model by progressively adding parameters to simpler models. We found that in several cases the choice of best‐fit model was affected by the parameter addition sequence. In terms of phylogenetic performance, there was little difference between models when the ratio of short: long terminal branches was 1:3 or less. However, under conditions of extreme terminal branch‐length variation, there were not only dramatic differences among models, but best‐fit models were always among the best at overcoming long‐branch attraction. The performance of minimum‐evolution‐distance methods was generally lower than that of discrete maximum‐likelihood methods, even if maximum‐likelihood methods were used to generate distance matrices. Correcting for among‐site rate variation was especially important for overcoming long‐branch attraction. The generality of our conclusions is supported by earlier simulation studies and by a preliminary analysis of mitochondrial and nuclear sequences from a well‐supported four‐taxon amniote phylogeny.

INDEPENDENT CONTRASTS SUCCEED WHERE ANCESTOR RECONSTRUCTION FAILS IN A KNOWN BACTERIOPHAGE PHYLOGENY

Abstract.— Methods of ancestor reconstruction are important tools for evolutionary inference that are difficult to test empirically because ancestral states are rarely known with certainty. We evaluated reconstruction methods for continuous phenotypic characters using taxa from an experimentally generated bacteriophage phylogeny. Except for one slowly evolving character, the estimated ancestral states of continuous phenotypic characters were highly inaccurate and biased, even when including a known ancestor at the root. This error was caused by a directional trend in character evolution and by rapid rates of character evolution. Computer simulations confirmed that such factors affect reconstruction of continuous characters in general. We also used phenotypic viral characters to evaluate two methods that attempt to estimate the correlation between characters during evolution. Whereas a nonphylogenetic regression was relatively inaccurate and biased, independent contrasts accurately estimated the correlation between characters with little bias.