Pinniped phylogeny and a new hypothesis for their origin and dispersal

The relationships and the zoogeography of the three extant pinniped families, Otariidae (sea lions and fur seals), Odobenidae (one extant species, the walrus), and Phocidae (true seals), have been contentious. Here, we address these topics in a molecular study that includes all extant species of true seals and sea lions, four fur seals and the walrus. Contrary to prevailing morphological views the analyses conclusively showed monophyletic Pinnipedia with a basal split between Otarioidea (Otariidae+Odobenidae) and Phocidae. The northern fur seal was the sister to all remaining otariids and neither sea lions nor arctocephaline fur seals were recognized as monophyletic entities. The basal Phocidae split between Monachinae (monk seals and southern true seals) and Phocinae (northern true seals) was strongly supported. The phylogeny of the Phocinae suggests that the ancestors of Cystophora (hooded seal) and the Phocini (e.g. harp seal, ringed seal) adapted to Arctic conditions and ice-breeding before 12 MYA (million years ago) as supported by the white natal coat of these lineages. The origin of the endemic Caspian and Baikal seals was dated well before the onset of major Pleistocene glaciations. The current findings, together with recent advances in pinniped paleontology, allow the proposal of a new hypothesis for pinniped origin and early dispersal. The hypothesis posits that pinnipeds originated on the North American continent with early otarioid and otariid divergences taking place in the northeast Pacific and those of the phocids in coastal areas of southeast N America for later dispersal to colder environments in the N Atlantic and the Arctic Basin, and in Antarctic waters.     

Historical biogeography and phylogeny of monachine seals (Pinnipedia: Phocidae) based on mitochondrial and nuclear DNA data

Aim To determine the origin and diversification of monachine seals using a phylogenetic framework. Methods Molecular sequence data from three mitochondrial genes (cyt b, ND1 and 12S), and one nuclear marker (an intron from the α‐lactalbumin gene) were examined from all extant species of monachine seals. Maximum likelihood and partitioned Bayesian inference were used to analyse separate and combined (mitochondrial + nuclear) data sets. Divergence times were estimated from the resultant phylogeny using nonparametric rate smoothing as implemented by the program r8s. Results Mirounga, Monachus and the Lobodontini form three well‐supported clades within a monophyletic Monachinae. Lobodontini + Mirounga form a clade sister to Monachus. Molecular divergence dates indicate that the first split within the Monachinae (Lobodontini + Mirounga clade and Monachus) occurred between 11.8 and 13.8 Ma and Mirounga, Monachus and the Lobodontini originated 2.7–3.4, 9.1–10.8 and 10.0–11.6 Ma, respectively. Main conclusions Two main clades exist within Monachinae, Monachus and Lobodontini + Mirounga. Monachus, a warm water clade, originated in the North Atlantic and maintained the temperate water affinities of their ancestors as they diversified in the subtropic regions of the Northern Hemisphere. The cold‐water clade, Lobodontini + Mirounga, dispersed southward to the cooler climates of the Southern Hemisphere. The Lobodontini continued south until reaching the Antarctic region where they diversified into the present‐day fauna. Mirounga shows an anti‐tropical distribution either reflective of a once cosmopolitan range that was separated by warming waters in the tropics or of transequatorial dispersal.     

A phylogeny of the extant Phocidae inferred from complete mitochondrial DNA coding regions

Despite extensive interest in the systematics of Pinnipedia, questions remain concerning phylogenetic relationships within the Phocidae or "true" seals. Relationships within the phocids and their placement relative to the remaining pinnipeds and major lineages of arctoid carnivores were examined using a large molecular data set consisting of 12 mitochondrial protein coding genes. Phylogenetic analysis including 15 extant species of the Phocidae, and representatives of the Otariidae, Odobenidae, Ursidae, Mustelidae, Canidae, and Felidae confirmed the monophyletic origins of the Pinnipedia within the Arctoidea. Slightly more support was found for an ursid affinity of the pinnipeds, however, this relationship remains contentious. The Phocidae were placed as the sister group to a common odobenid-otariid clade. Within the family Phocidae, strong support for the traditionally accepted subfamilies Phocinae (northern seals), and Monachinae (southern seals plus monk seals) was found. In contrast to recent suggestions, a monophyletic Monachus was strongly supported and was placed in a deep branching position within the Monachinae. Evidence from sequence divergence under a maximum likelihood model illustrated that the rarely used tribal distinction within the Monachinae are comparable, in terms of evolutionary distance, to accepted tribal distinctions within the Phocinae. In addition, results suggest that Pagophilus should be accepted as a genus within the Phocini. Sequence divergence between Phoca, Pusa, and Halichoerus is minimal, supporting a taxonomic reclassification of the three genera into an emended genus Phoca, without subgeneric distinctions.     

Tethyan vicariance, relictualism and speciation: evidence from a global molecular phylogeny of the opisthobranch genus Bulla

Aim Our aims were: (1) to reconstruct a molecular phylogeny of the cephalaspidean opisthobranch genus Bulla, an inhabitant of shallow sedimentary environments; (2) to test if divergence times are consistent with Miocene and later vicariance among the four tropical marine biogeographical provinces; (3) to examine the phylogenetic status of possible Tethyan relict species; and (4) to infer the timing and causes of speciation events. Location Tropical and warm‐temperate regions of the Atlantic, Indo‐West Pacific, Australasia and eastern Pacific. Methods Ten of the 12 nominal species of Bulla were sampled, in a total sample of 65 individuals, together with cephalaspidean outgroups. Phylogenetic relationships were inferred by Bayesian analysis of partial sequences of the mitochondrial cytochrome c oxidase I (COI) and 16S rRNA and nuclear 28S rRNA genes. Divergence times and rates of evolution were estimated using uncorrelated relaxed‐clock Bayesian methods with fossil calibrations (based on literature review and examination of fossil specimens), implemented in beast . The geographical pattern of speciation was assessed by estimating the degree of overlap between sister lineages. Results Four clades were supported: Indo‐West Pacific (four species), Australasia (one species), Atlantic plus eastern Pacific (three species) and Atlantic (two species), with estimated mean ages of 35–46 Ma. Nominal species were monophyletic, but deep divergences were found within one Indo‐West Pacific and one West Atlantic species. Species‐level divergences occurred in the Miocene or earlier. The age of a sister relationship across the Isthmus of Panama was estimated at 7.9–32.1 Ma, and the divergence of a pair of sister species on either side of the Atlantic Ocean occurred 20.4–27.2 Ma. Main conclusions Fossils suggest that Bulla originated in the Tethys realm during the Middle Eocene. Average ages of the four main clades fall in the Eocene, and far pre‐date the 18–19 Ma closure of the Tethys Seaway. This discrepancy could indicate earlier vicariant events, selective extinction or errors of calibration. Similarly, the transisthmian divergence estimate far pre‐dates the uplift of the Panamanian Isthmus at about 3 Ma. Speciation events occurred in the Miocene, consistent with tectonic events in the central Indo‐West Pacific, isolation of the Arabian Sea by upwelling and westward trans‐Atlantic dispersal. Differences in habitat between sister species suggest that ecological speciation may also have played a role. The basal position of the Australasian species supports its interpretation as a Tethyan relict.     

Gondwanan radiation of the Southern Hemisphere crayfishes (Decapoda: Parastacidae): evidence from fossils and molecules

Aim The sequential break‐up of Gondwana is thought to be a dominant process in the establishment of shared biota across landmasses of the Southern Hemisphere. Yet similar distributions are shared by taxa whose radiations clearly post‐date the Gondwanan break‐up. Thus, determining the contribution of vicariance versus dispersal to seemingly Gondwanan biota is complex. The southern freshwater crayfishes (family Parastacidae) are distributed on Australia and New Guinea, South America, Madagascar and New Zealand and are unlikely to have dispersed via oceans, owing to strict freshwater limitations. We test the hypotheses that the break‐up of Gondwana has led to (1) a predominately east–west (((Australia, New Zealand: 80 Ma) Madagascar: 160–121 Ma) South America: 165–140 Ma), or (2) a southern (((Australia, South America: 52–35 Ma) New Zealand: 80 Ma) Madagascar: 160–121 Ma) pattern for parastacid crayfish. Further, we examine the evidence for a complete drowning of New Zealand and subsequent colonization by freshwater crayfish. Location Southern Hemisphere. Methods The evolutionary relationships among the 15 genera of Parastacidae were reconstructed using mitochondrial [16S, cytochrome c oxidase subunit I (COI)] and nuclear (18S, 28S) sequence data and maximum likelihood and Bayesian methods of phylogenetic reconstruction. A Bayesian (multidivtime ) molecular dating method using six fossil calibrations and phylogenetic inference was used to estimate divergence time among crayfish clades on Gondwanan landmasses. Results The South American crayfish are monophyletic and a sister group to all other southern crayfish. Australian crayfish are not monophyletic, with two Tasmanian genera, Spinastacoides and Ombrastacoides, forming a clade with New Zealand and Malagasy crayfish (both monophyletic). Divergence of crayfish among southern landmasses is estimated to have occurred around the Late Jurassic to Early Cretaceous (109–178 Ma). Main conclusions The estimated phylogenetic relationships and time of divergence among the Southern Hemisphere crayfishes were consistent with an east–west pattern of Gondwanan divergence. The divergence between Australia and New Zealand (109–160 Ma) pre‐dated the rifting at around 80 Ma, suggesting that these lineages were established prior to the break‐up. Owing to the age of the New Zealand crayfish, we reject the hypothesis that there was a complete drowning of New Zealand crayfish habitat.     

Speciation chronology of rockhopper penguins inferred from molecular, geological and palaeoceanographic data

Aim The Southern Ocean is split into several biogeographical provinces between convergence zones that separate watermasses of different temperatures. Recent molecular phylogenies have uncovered a strong phylogeographic structure among rockhopper penguin populations, Eudyptes chrysocome sensu lato, from different biogeographical provinces. These studies suggested a reclassification as three species in two major clades, corresponding, respectively, to warm, subtropical and cold sub‐Antarctic watermasses rather than to geographic proximity. Such a phylogeographic pattern, also observed in plants, invertebrates and fishes of the Southern Ocean, suggests that past changes in the positions of watermasses may have affected the evolutionary history of penguins. We calculated divergence times among various rockhopper penguin clades and calibrated these data with palaeomagmatic and palaeoceanographic events to generate a speciation chronology in rockhopper penguins. Location Southern Ocean. Methods Divergence times between populations were calculated using five distinct mitochondrial DNA loci, and assuming a molecular clock model as implemented in mdiv . The molecular evolution rate of rockhopper penguins was calibrated using the radiochronological age of St Paul Island and Amsterdam Island in the southern Indian Ocean. Separations within other clades were correlated with palaeoceanographic data using this calibrated rate. Results The split between the Atlantic and Indian populations of rockhopper penguins was dated as 0.25 Ma, using the date of emergence of St Paul and Amsterdam islands, and the divergence between sub‐Antarctic and subtropical rockhopper penguins was dated as c. 0.9 Ma (i.e. during the mid‐Pleistocene transition, a major change in the Earth’s climate cycles). Main conclusions The mid‐Pleistocene transition is known to have caused a major southward shift in watermasses in the Southern Ocean, thus changing the environment around the northernmost rockhopper penguin breeding sites. This ecological isolation of northernmost populations may have caused vicariant speciation, splitting the species into two major clades. After the emergence of St Paul and Amsterdam islands in the subtropical Indian Ocean 0.25 Ma, these islands were colonized by penguins from the subtropical Atlantic, 6000 km away, rather than by penguins from the sub‐Antarctic Indian Ocean, 5000 km closer.     

Molecular lineage diversity and inter‐oceanic biogeographical history in Hiatella (Mollusca,Bivalvia)

Hiatella is one of the most widespread marine bivalve genera, occurring in diverse habitats from the temperate to polar latitudes in both hemispheres, and in fossil strata since almost 150 Myr ago. Despite variation in some biological and morphological traits, characters to resolve the current systematic structure consistently across the range of the genus are not known: all samples are often referred to a single species, Hiatella arctica (L.). Exploring the systematics of Hiatella using partial sequences of three genes (mitochondrial COI, and the nuclear ANT and 28S rRNA), we find high diversity of deep lineages (11–22% p‐distance in COI), and identify at least 13 distinct taxa both by heuristic criteria (congruence of the nuclear and mtDNA data) and by coalescence‐based analyses. At several localities, two or three of these cryptic species were found in sympatry. In the framework of previous fossil evidence and of hypotheses of paleoceanographical connections, scenarios of the phylogeny and biogeographical history of the identified species at a range of different time scales are outlined. The distinction between the main North Pacific and North Atlantic Hiatella clades and systematic diversification within each of them seem to have followed a Miocene trans‐Panamanian invasion. Apart from such earlier intra‐basin diversification, the data suggest that three successive counter‐invasions from the Pacific to the Atlantic via the Arctic Ocean route have later contributed to the current North Atlantic Hiatella diversity. These invasions probably took place in connection with (i) the Great Trans‐Arctic Biotic Interchange in the Pliocene, (ii) the last interglacial period c. 120 kya and (iii) the Holocene, postdating the last glaciation. This sequence of trans‐Arctic invasions is largely analogous to that hypothesized for some other boreal‐arctic bivalves.     

A molecular view of pinniped relationships with particular emphasis on the true seals

Phylogenetic analysis of conservative nucleotide substitutions in 18 complete sequences of the mitochondrial cytochrome gene of Phocidae (true seals), Odobenidae (walruses), and Otariidae (sea lions and fur seals), plus three ursid and three felid sequences, identified the pinnipeds as monophyletic with Otariidae and Odobenidae on a common evolutionary branch. Analysis of total nucleotide differences separated the evolutionary lineages of northern and southern phocids. Both lineages are distinct from the most ancestral phocid genus, Monachus (monk seals), represented by the Hawaiian monk seal. The inclusion of the Hawaiian monk seal in the subfamily Monachinae makes the subfamily paraphyletic. Among the northern phocids, the hooded seal (genus Cystophora, chromosome number 2n = 34) is sister taxon to the Phoca complex. The Phoca complex, which is characterized by the chromosome number 2n = 32, includes genus Phoca and the monotypic genus Halichoerus (grey seal). The comparison does not support a generic distinction of Halichoerus within the Phoca complex. The present data suggest that Cystophora and Phoca separated 6 million years ago. Among the southern phocids the close molecular relationship of the Weddell and leopard seals relative to their morphological distinction exemplifies rapid adaptation to different ecological niches. This result stands in contrast to the limited morphological differentiation relative to the pronounced molecular distinctions that may occur within the Phoca complex.Correspondence to: Ú. Árnason     

Bayesian relaxed clock estimation of divergence times in foraminifera

Accurate and precise estimation of divergence times during the Neo-Proterozoic is necessary to understand the speciation dynamic of early Eukaryotes. However such deep divergences are difficult to date, as the molecular clock is seriously violated. Recent improvements in Bayesian molecular dating techniques allow the relaxation of the molecular clock hypothesis as well as incorporation of multiple and flexible fossil calibrations. Divergence times can then be estimated even when the evolutionary rate varies among lineages and even when the fossil calibrations involve substantial uncertainties. In this paper, we used a Bayesian method to estimate divergence times in Foraminifera, a group of unicellular eukaryotes, known for their excellent fossil record but also for the high evolutionary rates of their genomes. Based on multigene data we reconstructed the phylogeny of Foraminifera and dated their origin and the major radiation events. Our estimates suggest that Foraminifera emerged during the Cryogenian (650-920 Ma, Neo-Proterozoic), with a mean time around 770 Ma, about 220 Myr before the first appearance of reliable foraminiferal fossils in sediments (545 Ma). Most dates are in agreement with the fossil record, but in general our results suggest earlier origins of foraminiferal orders. We found that the posterior time estimates were robust to specifications of the prior. Our results highlight inter-species variations of evolutionary rates in Foraminifera. Their effect was partially overcome by using the partitioned Bayesian analysis to accommodate rate heterogeneity among data partitions and using the relaxed molecular clock to account for changing evolutionary rates. However, more coding genes appear necessary to obtain more precise estimates of divergence times and to resolve the conflicts between fossil and molecular date estimates.     

Dinucleotide microsatellite markers from the Antarctic seals and their use in other Pinnipeds

Twenty‐four microsatellite loci were isolated from three species of Antarctic seals (Subfamily Monachinae, Tribe Lobodontini). Eleven loci were cloned from Weddell seal, Leptonychotes weddellii, seven from leopard seal, Hydrurga leptonyx, and six from crabeater seal, Lobodon carcinophagus. Variability was assessed in Weddell seals collected in McMurdo Sound, leopard seals from Bird Island, South Georgia, and crabeater seals sampled in the eastern Ross Sea. All loci were variable in the three species used for cloning and 22 of these loci amplified variable products in the Ross seal, Ommatophoca rossii. Cross‐species amplification was largely successful, with an average of 19 loci amplifying products in other phocids.     

Forest elephant mitochondrial genomes reveal that elephantid diversification in Africa tracked climate transitions

Among elephants, the phylogeographic patterns of mitochondrial (mt) and nuclear markers are often incongruent. One hypothesis attributes this to sex differences in dispersal and in the variance of reproductive success. We tested this hypothesis by examining the coalescent dates of genetic markers within elephantid lineages, predicting that lower dispersal and lower variance in reproductive success among females would have increased mtDNA relative to nuclear coalescent dates. We sequenced the mitochondrial genomes of two forest elephants, aligning them to mitogenomes of African savanna and Asian elephants, and of woolly mammoths, including the most divergent mitogenomes within each lineage. Using fossil calibrations, the divergence between African elephant F and S clade mitochondrial genomes (originating in forest and savanna elephant lineages, respectively) was estimated as 5.5 Ma. We estimated that the (African) ancestor of the mammoth and Asian elephant lineages diverged 6.0 Ma, indicating that four elephantid lineages had differentiated in Africa by the Miocene–Pliocene transition, concurrent with drier climates. The coalescent date for forest elephant mtDNAs was c. 2.4 Ma, suggesting that the decrease in tropical forest cover during the Pleistocene isolated distinct African forest elephant lineages. For all elephantid lineages, the ratio of mtDNA to nuclear coalescent dates was much greater than 0.25. This is consistent with the expectation that sex differences in dispersal and in variance of reproductive success would have increased the effective population size of mtDNA relative to nuclear markers in elephantids, contributing to the persistence of incongruent mtDNA phylogeographic patterns.     

Phylogeographic analysis of Pimoidae (Arachnida: Araneae) inferred from mitochondrial cytochrome c oxidase subunit I and nuclear 28S rRNA gene regions

Using mitochondrial DNA cytochrome c oxidase subunit I and nuclear DNA 28S rRNA data, we explored the phylogenetic relationships of the family Pimoidae (Arachnida: Araneae) and tested the North America to Asia dispersal hypothesis. Sequence data were analysed using maximum parsimony and Bayesian inference. A phylogenetic analysis suggested that vicariance, instead of dispersal, better explained the present distribution pattern of Pimoidae. Times of divergence events were estimated using penalized likelihood method. The dating analysis suggested that the emergence time of Pimoidae was approximately 140 million years ago (Ma). The divergence time of the North American and Asian species of Pimoa was approximately 110 Ma. Our phylogenetic hypothesis supports the current morphology‐based taxonomy and suggests that the cave dwelling might have played an important role in the speciation of pimoids in arid areas.     

Origins of the Greenland shark (Somniosus microcephalus): Impacts of ice‐olation and introgression

Herein, we use genetic data from 277 sleeper sharks to perform coalescent‐based modeling to test the hypothesis of early Quaternary emergence of the Greenland shark (Somniosus microcephalus) from ancestral sleeper sharks in the Canadian Arctic‐Subarctic region. Our results show that morphologically cryptic somniosids S. microcephalus and Somniosus pacificus can be genetically distinguished using combined mitochondrial and nuclear DNA markers. Our data confirm the presence of genetically admixed individuals in the Canadian Arctic and sub‐Arctic, and temperate Eastern Atlantic regions, suggesting introgressive hybridization upon secondary contact following the initial species divergence. Conservative substitution rates fitted to an Isolation with Migration (IM) model indicate a likely species divergence time of 2.34 Ma, using the mitochondrial sequence DNA, which in conjunction with the geographic distribution of admixtures and Pacific signatures likely indicates speciation associated with processes other than the closing of the Isthmus of Panama. This time span coincides with further planetary cooling in the early Quaternary period followed by the onset of oscillating glacial‐interglacial cycles. We propose that the initial S. microcephalus–S. pacificus split, and subsequent hybridization events, were likely associated with the onset of Pleistocene glacial oscillations, whereby fluctuating sea levels constrained connectivity among Arctic oceanic basins, Arctic marginal seas, and the North Atlantic Ocean. Our data demonstrates support for the evolutionary consequences of oscillatory vicariance via transient oceanic isolation with subsequent secondary contact associated with fluctuating sea levels throughout the Quaternary period—which may serve as a model for the origins of Arctic marine fauna on a broad taxonomic scale.     

Pleistocene genetic connectivity in a widespread,open‐habitat‐adapted mosquito in the Indo‐Oriental region

Aim The environmental effect of Pleistocene climatic change in the Indo‐Oriental region has resulted in allopatric fragmentation and the generation of diversity in forest‐associated species. The aim of this study was to determine the extent to which Pleistocene climatic change has resulted in the fragmentation and speciation of an open‐habitat‐adapted mosquito, Anopheles vagus s.l., across its range. Location Anopheles vagus s.l. was sampled across the Indo‐Oriental region. Methods We generated 116 mitochondrial cytochrome c oxidase subunit I (COI) and 121 nuclear internal transcribed spacer 2 (ITS2) DNA sequences from 18 populations. Relationships between mitochondrial haplotypes were reconstructed using minimum spanning networks, and population structure was examined using analyses of molecular variance. The population history, including lineage divergence times, population expansion and gene flow, was inferred using beast and the isolation with migration (IM) model. Results There was no evidence to support the presence of the endemic Philippines species, A. limosus; instead, Philippine populations were closely related to, and derived from, A. vagus on the eastern Southeast Asian mainland. The most distinct populations were those from Java and East Timor, which differed from all other populations by all individuals having a 4‐bp insertion in the ITS2 sequence. The corresponding mitochondrial haplotypes had an estimated divergence time of 2.6 Ma [95% confidence interval (CI) 1.9–3.6 Ma]. Haplotype networks and analysis of molecular variance for COI supported western (Sri Lanka, India and Myanmar) and eastern (Thailand, Singapore, Cambodia, Vietnam and the Philippines) population groupings. This grouping structure results from the divergence of an eastern and a western mitochondrial lineage, estimated to have occurred 0.37 Ma (95% CI 0.26–0.55 Ma). Subsequent migration from the east to the west (0.16 Ma) is inferred to have created an admixture zone in Myanmar and Thailand. Main conclusions With the possible exception of populations from Java and East Timor, A. vagus appears to be one widespread genetically diverse taxon across its extensive range. The abundance of grassland during long interglacial periods may have facilitated population connectivity and range expansion across the Oriental and western Australasian regions.     

The Mediterranean monk seal Monachus monachus: status,biology, threats,and conservation priorities

相似文献   

Fossil calibration of molecular clocks and the divergence times of geminate species pairs separated by the isthmus of panama

Calibration of nucleotide sequence divergence rates provides an important method by which to test many hypotheses of evolution. In the absence of an adequate fossil record, geological events, rather than the first appearances of sister taxa in the geological record, are often used to calibrate molecular clocks. The formation of the Isthmus of Panama, which isolated the tropical western Atlantic and eastern Pacific oceans, is one such event that is frequently used to infer rates of nucleotide sequence divergence. Isthmian calibrations assume that morphologically similar "geminate" species living now on either side of the isthmus were isolated geographically by the latest stages of seaway closure 3.1-3.5 MYA. Here, I have applied calibration dates from the fossil record to cytochrome c oxidase-1 (CO1) and nuclear histone-3 (H3) divergences among six pairs of geminates in the Arcidae to test this hypothesis. Analysis of CO1 first and third positions yield geminate divergences that predate final seaway closure, and on the basis of CO1 first positions, times for all six geminates are significantly greater than 3.5 Myr. H3 sequences produce much more recent geminate divergences, some that are younger than 3.1 Myr. But H3-derived estimates for all arcid geminates are not significantly different from both 0 and 15 Myr. According to CO1, one of the two most divergent pairs, Arca mutabilis and A. imbricata, split more than 30 MYA. This date is compatible with the fossil record, which indicates that these species were morphologically distinct at least 16-21 MYA. Across all CO1 nucleotide sites, divergence rates for arcids are slower than the rates reported for other taxa on the basis of isthmian calibrations, with the exception of rates determined from the least divergent species pair in larger surveys of multiple transisthmian pairs. Rate differences between arcids and some taxa may be real, but these data suggest that divergence rates can be greatly overestimated when dates corresponding to final closure of the Central American Seaway are used to calibrate the molecular clocks of marine organisms.     

Variation in blubber fatty acid composition among marine mammals in the Canadian Arctic

The composition of predator adipose stores can provide important insights into foraging patterns and the ecological relationships among species. We determined the fatty acid (FA) composition of 843 blubber samples from 80 bearded seals (Erignathus barbatus), 33 harbor seals (Phoca vitulina), 239 harp seals (Pagophilus groenlandicus), 32 hooded seals (Cystophora cristata), 281 ringed seals (Phoca hispida), 53 walruses (Odobenus rosmarus rosmarus), 105 beluga whales (Delphinapterus leucas), and 20 narwhals (Monodon monoceros) across the Canadian Arctic to examine patterns of variability among and within species. FA signatures accurately distinguished phocid seals, walruses, and whales. Belugas and narwhals had the most similar FA signatures of any two species, suggesting substantial overlap in their diets, especially in the narwhal‐wintering area off eastern Baffin Island. Among phocid seals, harp and hooded seals had the most similar FA signatures. Bearded seals were most similar to walruses, which was consistent with the benthic feeding habits of both species. Within species, geographic differences in FA signatures were found over both large (>4,000 km) and small (<100 km) spatial scales. Overall, within‐species differences were smaller than among‐species differences. In general, FA signature patterns were consistent with previous studies of the ecology and diets of arctic marine mammals.     

Using molecules and morphology to infer the phylogenetic relationships and evolutionary history of the Dirini (Nymphalidae: Satyrinae), a tribe of butterflies endemic to Southern Africa

The first empirically supported phylogenetic hypothesis of relationships for the southern African endemic butterfly tribe Dirini is presented. Data derived from the morphology and ecology of the adults and immature stages (33 characters), and portions of the mitochondrial gene cytochrome oxidase I (COI) and the nuclear genes elongation factor 1α (EF1α) and wingless (WG) (totalling 1734 bp) were used to infer the relationships of the in‐group genera. An expanded molecular dataset using four genera from the Nymphalini and Satyrini to root the tree, and three genera from the Melanitini to test the monophyly of the tribe, was analysed using parsimony and Bayesian methods. Estimates of divergence times were calculated using two fossil calibrations under a relaxed molecular clock model. The monophyly of the tribe and each in‐group genus were strongly supported. Key findings are the sister‐taxon relationship of Aeropetes and Tarsocera, the apparent simultaneous or nearly simultaneous radiation of four lineages, the polyphyly of the species within Torynesis, and the apparent trans‐Atlantic dispersal of the ancestors of Manataria about 40 Ma. Estimates of divergence times indicate that the tribe has undergone two major radiations since its origin: the first when they left forest habitats in the mid–late Oligocene, shortly after the radiation of the grasses (Poaceae), and the second in the early‐middle Pliocene, coinciding with the aridification of southern Africa and the spread of conditions that favoured C₄ grasses over the C₃ grasses that dirine larvae prefer to eat. The high species diversity within the tribe appears to be partly a taxonomic artefact that may have resulted from the misinterpretation of climate‐related phenotypic variation within extant species. Relocation and breeding experiments should test this hypothesis.     

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.     

Perpendicular axes of differentiation generated by mitochondrial introgression

Differential introgression of mitochondrial vs. nuclear DNA generates discordant patterns of geographic variation and can promote population divergence and speciation. We examined a potential case of mitochondrial introgression leading to two perpendicular axes of differentiation. The Eastern Yellow Robin Eopsaltria australis, a widespread Australian bird, shows a deep mitochondrial split that is perpendicular to north–south nuclear DNA and plumage colour differentiation. We propose a scenario to explain this pattern: (i) first, both nuclear and mitochondrial genomes differentiated in concert during north–south population divergence; (ii) later, their histories disconnected after two mitochondrial introgression events resulting in a deep mitochondrial split perpendicular to the nuclear DNA structure. We explored this scenario by coalescent modelling of ten mitochondrial genes and 400 nuclear DNA loci. Initial mitochondrial and nuclear genome divergences were estimated to have occurred in the early Pleistocene, consistent with the proposed scenario. Subsequent climatic transitions may have driven later mitochondrial introgression. We consider neutral introgression unlikely and instead propose that the evidence is more consistent with adaptive mitochondrial introgression and selection against incompatible mitochondrial‐nuclear combinations. This likely generated an axis of coastal‐inland mitochondrial differentiation in the face of nuclear gene flow, perpendicular to the initial north–south axis of differentiation (reflected in genomewide nuclear DNA and colour variation).