Divergence Date Estimates Research Articles

Papio Cranium from the Hominin-Bearing Site of Malapa: Implications for the Evolution of Modern Baboon Cranial Morphology and South African Plio-Pleistocene Biochronology.

A new partial cranium (UW 88-886) of the Plio-Pleistocene baboon Papio angusticeps from Malapa is identified, described and discussed. UW 88-886 represents the only non-hominin primate yet recovered from Malapa and is important both in the context of baboon evolution as well as South African hominin site biochronology. The new specimen may represent the first appearance of modern baboon anatomy and coincides almost perfectly with molecular divergence date estimates for the origin of the modern P. hamadryas radiation. The fact that the Malapa specimen is dated between ~2.026–2.36 million years ago (Ma) also has implications for the biochronology of other South African Plio-Pleistocene sites where P. angusticeps is found.

Molecular Phylogeny of Grassland Caterpillars (Lepidoptera: Lymantriinae: Gynaephora) Endemic to the Qinghai-Tibetan Plateau.

Gynaephora (Lepidoptera Erebidae: Lymantriinae) is a small genus, consisting of 15 nominated species, of which eight species are endemic to the Qinghai-Tibetan Plateau (QTP). In this study, we employed both mitochondrial and nuclear loci to infer a molecular phylogeny for the eight QTP Gynaephora spp. We used the phylogeny to estimate divergence dates in a molecular dating analysis and to delimit species. This information allowed us to investigate associations between the diversification history of the eight QTP species and geological and climatic events. Phylogenetic analyses indicated that the eight QTP species formed a monophyletic group with strong supports in both Bayesian and maximum likelihood analyses. The low K2P genetic distances between the eight QTP species suggested that diversification occurred relatively quickly and recently. Out of the eight species, five species were highly supported as monophyletic, which were also recovered by species delimitation analyses. Samples of the remaining three species (G. aureata, G. rouergensis, and G. minora) mixed together, suggesting that further studies using extensive population sampling and comprehensive morphological approaches are necessary to clarify their species status. Divergence time estimation results demonstrated that the diversification and speciation of Gynaephora on the QTP began during the late Miocene/early Pliocene and was potentially affected by the QTP uplift and associated climate changes during this time.

Comparative genomics and phylogenetic discordance of cultivated tomato and close wild relatives.

Background. Studies of ancestry are difficult in the tomato because it crosses with many wild relatives and species in the tomato clade that have diverged very recently. As a result, the phylogeny in relation to its closest relatives remains uncertain. By using the coding sequence from Solanum lycopersicum, S. galapagense, S. pimpinellifolium, S. corneliomuelleri, and S. tuberosum and the genomic sequence from S. lycopersicum ‘Heinz’, an heirloom line, S. lycopersicum ‘Yellow Pear’, and two of cultivated tomato’s closest relatives, S. galapagense and S. pimpinellifolium, we have aimed to resolve the phylogenies of these closely related species as well as identify phylogenetic discordance in the reference cultivated tomato.Results. Divergence date estimates suggest that the divergence of S. lycopersicum, S. galapagense, and S. pimpinellifolium happened less than 0.5 MYA. Phylogenies based on 8,857 coding sequences support grouping of S. lycopersicum and S. galapagense, although two secondary trees are also highly represented. A total of 25 genes in our analysis had sites with evidence of positive selection along the S. lycopersicum lineage. Whole genome phylogenies showed that while incongruence is prevalent in genomic comparisons between these genotypes, likely as a result of introgression and incomplete lineage sorting, a primary phylogenetic history was strongly supported.Conclusions. Based on analysis of these genotypes, S. galapagense appears to be closely related to S. lycopersicum, suggesting they had a common ancestor prior to the arrival of an S. galapagense ancestor to the Galápagos Islands, but after divergence of the sequenced S. pimpinellifolium. Genes showing selection along the S. lycopersicum lineage may be important in domestication or selection occurring post-domestication. Further analysis of intraspecific data in these species will help to establish the evolutionary history of cultivated tomato. The use of an heirloom line is helpful in deducing true phylogenetic information of S. lycopersicum and identifying regions of introgression from wild species.

Differential introgression and effective size of marker type influence phylogenetic inference of a recently divergent avian group (Phasianidae: Tympanuchus)

Life history strategies can influence the effective population size (Ne) of loci differently based on their mode of inheritance. Recognizing how this may affect the rate of lineage sorting among marker types is important for studies focused on resolving phylogenetic relationships among recently divergent taxa. In this study, we use gene tree, coalescent-based species tree, and isolation-with-migration analyses to explore the differences between marker types (autosomal, Z-linked, and mitochondrial) in resolving phylogenetic relationships among North American prairie grouse (Tympanuchus). We found that Z-linked loci were more likely to identify monophyletic relationships among prairie grouse species compared to autosomal and mtDNA loci in both species and gene tree analyses, with species tree analyses outperforming gene trees. These results were further supported with isolation-with-migration analyses, where Z-linked loci largely followed a strict isolation model while autosomal loci were more likely to fit a model with gene flow between species following population divergence. While accounting for differences in inheritance pattern (or Ne) for marker type, results suggest that additional factors, such as strong sexual selection and sex-biased introgression (i.e., male-biased postzygotic hybrid behavioral isolation or “unsexy son”), may further explain the decreased diversity levels and increased rate of lineage sorting observed with the Z-linked loci relative to autosomal and mtDNA loci. In fact, to our knowledge no hybrid male prairie grouse have been observed breeding in the wild, yet hybrid females along with backcross females are known to produce viable offspring. Overall, this study highlights that more work is needed to determine how complex models of gene flow (i.e., sex biased introgression) and differences in the effective size among marker types based on differing life history strategies influence divergence date estimation and species delimitation.

A dated molecular phylogeny of manta and devil rays (Mobulidae) based on mitogenome and nuclear sequences

Manta and devil rays are an iconic group of globally distributed pelagic filter feeders, yet their evolutionary history remains enigmatic. We employed next generation sequencing of mitogenomes for nine of the 11 recognized species and two outgroups; as well as additional Sanger sequencing of two mitochondrial and two nuclear genes in an extended taxon sampling set. Analysis of the mitogenome coding regions in a Maximum Likelihood and Bayesian framework provided a well-resolved phylogeny. The deepest divergences distinguished three clades with high support, one containing Manta birostris, Manta alfredi, Mobula tarapacana, Mobula japanica and Mobula mobular; one containing Mobula kuhlii, Mobula eregoodootenkee and Mobula thurstoni; and one containing Mobula munkiana, Mobula hypostoma and Mobula rochebrunei. Mobula remains paraphyletic with the inclusion of Manta, a result that is in agreement with previous studies based on molecular and morphological data. A fossil-calibrated Bayesian random local clock analysis suggests that mobulids diverged from Rhinoptera around 30 Mya. Subsequent divergences are characterized by long internodes followed by short bursts of speciation extending from an initial episode of divergence in the Early and Middle Miocene (19–17 Mya) to a second episode during the Pliocene and Pleistocene (3.6 Mya – recent). Estimates of divergence dates overlap significantly with periods of global warming, during which upwelling intensity – and related high primary productivity in upwelling regions – decreased markedly. These periods are hypothesized to have led to fragmentation and isolation of feeding regions leading to possible regional extinctions, as well as the promotion of allopatric speciation. The closely shared evolutionary history of mobulids in combination with ongoing threats from fisheries and climate change effects on upwelling and food supply, reinforces the case for greater protection of this charismatic family of pelagic filter feeders.

Revisiting the phylogenetic relationships, biogeography, and taxonomy of spider monkeys (genus Ateles) in light of new molecular data

Spider monkeys (Ateles) are one of the most endangered groups of primates in the Neotropics. The genus is widely distributed from Mexico to the north of Bolivia and includes many morphologically distinct forms in terms of pelage color and patterning. The taxonomy, phylogenetic relationships, and biogeographic history of the genus have been subject to much debate, making scientific communication difficult and creating challenges for conservation actions. We extracted DNA from samples of all currently recognized species of spider monkeys collected from across the geographic range of the genus, sequenced ∼3.5kilobases of coding sequence from the mitochondrial genome, and used this large dataset to (a) infer the phylogenetic relationships among the different forms of spider monkeys, (b) evaluate whether currently recognized species of spider monkeys form reciprocally monophyletic groups that are concordant with contemporary classifications, and (c) estimate divergence dates among the different lineages of Ateles. We found that all proposed species of spider monkeys for which we have samples from multiple localities indeed appear to form monophyletic groups. However, in contrast to previous studies, several of our analyses robustly inferred Ateles marginatus from northeast Brazil as the sister taxon to all other spider monkeys. A Bayesian dating analysis suggests that the most recent common ancestor of extant Ateles dates to ∼6.7Ma, in the late Miocene, and most species-level splits within the genus took place in the late Pliocene, suggesting that the modern diversity in spider monkeys cannot be explained principally by isolation and divergence of populations in forest refugia during the Pleistocene. Based on our new phylogenetic inference and dating analysis, we propose a revised biogeographic scenario for the evolution of this genus.

Ancient dates or accelerated rates? Morphological clocks and the antiquity of placental mammals.

Analyses of a comprehensive morphological character matrix of mammals using 'relaxed' clock models (which simultaneously estimate topology, divergence dates and evolutionary rates), either alone or in combination with an 8.5 kb nuclear sequence dataset, retrieve implausibly ancient, Late Jurassic-Early Cretaceous estimates for the initial diversification of Placentalia (crown-group Eutheria). These dates are much older than all recent molecular and palaeontological estimates. They are recovered using two very different clock models, and regardless of whether the tree topology is freely estimated or constrained using scaffolds to match the current consensus placental phylogeny. This raises the possibility that divergence dates have been overestimated in previous analyses that have applied such clock models to morphological and total evidence datasets. Enforcing additional age constraints on selected internal divergences results in only a slight reduction of the age of Placentalia. Constraining Placentalia to less than 93.8 Ma, congruent with recent molecular estimates, does not require major changes in morphological or molecular evolutionary rates. Even constraining Placentalia to less than 66 Ma to match the 'explosive' palaeontological model results in only a 10- to 20-fold increase in maximum evolutionary rate for morphology, and fivefold for molecules. The large discrepancies between clock- and fossil-based estimates for divergence dates might therefore be attributable to relatively small changes in evolutionary rates through time, although other explanations (such as overly simplistic models of morphological evolution) need to be investigated. Conversely, dates inferred using relaxed clock models (especially with discrete morphological data and MrBayes) should be treated cautiously, as relatively minor deviations in rate patterns can generate large effects on estimated divergence dates.

The origins and radiation of Australian Coptotermes termites: From rainforest to desert dwellers

The termite genus Coptotermes (Rhinotermitidae) is found in Asia, Africa, Central/South America and Australia, with greatest diversity in Asia. Some Coptotermes species are amongst the world’s most damaging invasive termites, but the genus is also significant for containing the most sophisticated mound-building termites outside the family Termitidae. These mound-building Coptotermes occur only in Australia. Despite its economic and evolutionary significance, the biogeographic history of the genus has not been well investigated, nor has the evolution of the Australian mound-building species. We present here the first phylogeny of the Australian Coptotermes to include representatives from all described species. We combined our new data with previously generated data to estimate the first phylogeny to include representatives from all continents where the genus is found. We also present the first estimation of divergence dates during the evolution of the genus. We found the Australian Coptotermes to be monophyletic and most closely related to the Asian Coptotermes, with considerable genetic diversity in some Australian taxa possibly representing undescribed species. The Australian mound-building species did not form a monophyletic clade. Our ancestral state reconstruction analysis indicated that the ancestral Australian Coptotermes was likely to have been a tree nester, and that mound-building behaviour has arisen multiple times. The Australian Coptotermes were found to have diversified ∼13million years ago, which plausibly matches with the narrowing of the Arafura Sea allowing Asian taxa to cross into Australia. The first diverging Coptotermes group was found to be African, casting doubt on the previously raised hypothesis that the genus has an Asian origin.

Methods for the quantitative comparison of molecular estimates of clade age and the fossil record.

Approaches quantifying the relative congruence, or incongruence, of molecular divergence estimates and the fossil record have been limited. Previously proposed methods are largely node specific, assessing incongruence at particular nodes for which both fossil data and molecular divergence estimates are available. These existing metrics, and other methods that quantify incongruence across topologies including entirely extinct clades, have so far not taken into account uncertainty surrounding both the divergence estimates and the ages of fossils. They have also treated molecular divergence estimates younger than previously assessed fossil minimum estimates of clade age as if they were the same as cases in which they were older. However, these cases are not the same. Recovered divergence dates younger than compared oldest known occurrences require prior hypotheses regarding the phylogenetic position of the compared fossil record and standard assumptions about the relative timing of morphological and molecular change to be incorrect. Older molecular dates, by contrast, are consistent with an incomplete fossil record and do not require prior assessments of the fossil record to be unreliable in some way. Here, we compare previous approaches and introduce two new descriptive metrics. Both metrics explicitly incorporate information on uncertainty by utilizing the 95% confidence intervals on estimated divergence dates and data on stratigraphic uncertainty concerning the age of the compared fossils. Metric scores are maximized when these ranges are overlapping. MDI (minimum divergence incongruence) discriminates between situations where molecular estimates are younger or older than known fossils reporting both absolute fit values and a number score for incompatible nodes. DIG range (divergence implied gap range) allows quantification of the minimum increase in implied missing fossil record induced by enforcing a given set of molecular-based estimates. These metrics are used together to describe the relationship between time trees and a set of fossil data, which we recommend be phylogenetically vetted and referred on the basis of apomorphy. Differences from previously proposed metrics and the utility of MDI and DIG range are illustrated in three empirical case studies from angiosperms, ostracods, and birds. These case studies also illustrate the ways in which MDI and DIG range may be used to assess time trees resultant from analyses varying in calibration regime, divergence dating approach or molecular sequence data analyzed.

The fossilized birth-death process for coherent calibration of divergence-time estimates.

Time-calibrated species phylogenies are critical for addressing a wide range of questions in evolutionary biology, such as those that elucidate historical biogeography or uncover patterns of coevolution and diversification. Because molecular sequence data are not informative on absolute time, external data--most commonly, fossil age estimates--are required to calibrate estimates of species divergence dates. For Bayesian divergence time methods, the common practice for calibration using fossil information involves placing arbitrarily chosen parametric distributions on internal nodes, often disregarding most of the information in the fossil record. We introduce the "fossilized birth-death" (FBD) process--a model for calibrating divergence time estimates in a Bayesian framework, explicitly acknowledging that extant species and fossils are part of the same macroevolutionary process. Under this model, absolute node age estimates are calibrated by a single diversification model and arbitrary calibration densities are not necessary. Moreover, the FBD model allows for inclusion of all available fossils. We performed analyses of simulated data and show that node age estimation under the FBD model results in robust and accurate estimates of species divergence times with realistic measures of statistical uncertainty, overcoming major limitations of standard divergence time estimation methods. We used this model to estimate the speciation times for a dataset composed of all living bears, indicating that the genus Ursus diversified in the Late Miocene to Middle Pliocene.

Horizontal transfer of an adaptive chimeric photoreceptor from bryophytes to ferns.

Ferns are well known for their shade-dwelling habits. Their ability to thrive under low-light conditions has been linked to the evolution of a novel chimeric photoreceptor--neochrome--that fuses red-sensing phytochrome and blue-sensing phototropin modules into a single gene, thereby optimizing phototropic responses. Despite being implicated in facilitating the diversification of modern ferns, the origin of neochrome has remained a mystery. We present evidence for neochrome in hornworts (a bryophyte lineage) and demonstrate that ferns acquired neochrome from hornworts via horizontal gene transfer (HGT). Fern neochromes are nested within hornwort neochromes in our large-scale phylogenetic reconstructions of phototropin and phytochrome gene families. Divergence date estimates further support the HGT hypothesis, with fern and hornwort neochromes diverging 179 Mya, long after the split between the two plant lineages (at least 400 Mya). By analyzing the draft genome of the hornwort Anthoceros punctatus, we also discovered a previously unidentified phototropin gene that likely represents the ancestral lineage of the neochrome phototropin module. Thus, a neochrome originating in hornworts was transferred horizontally to ferns, where it may have played a significant role in the diversification of modern ferns.

Primate phylogenetic relationships and divergence dates inferred from complete mitochondrial genomes

The origins and the divergence times of the most basal lineages within primates have been difficult to resolve mainly due to the incomplete sampling of early fossil taxa. The main source of contention is related to the discordance between molecular and fossil estimates: while there are no crown primate fossils older than 56Ma, most molecule-based estimates extend the origins of crown primates into the Cretaceous. Here we present a comprehensive mitogenomic study of primates. We assembled 87 mammalian mitochondrial genomes, including 62 primate species representing all the families of the order. We newly sequenced eleven mitochondrial genomes, including eight Old World monkeys and three strepsirrhines. Phylogenetic analyses support a strong topology, confirming the monophyly for all the major primate clades. In contrast to previous mitogenomic studies, the positions of tarsiers and colugos relative to strepsirrhines and anthropoids are well resolved. In order to improve our understanding of how fossil calibrations affect age estimates within primates, we explore the effect of seventeen fossil calibrations across primates and other mammalian groups and we select a subset of calibrations to date our mitogenomic tree. The divergence date estimates of the Strepsirrhine/Haplorhine split support an origin of crown primates in the Late Cretaceous, at around 74Ma. This result supports a short-fuse model of primate origins, whereby relatively little time passed between the origin of the order and the diversification of its major clades. It also suggests that the early primate fossil record is likely poorly sampled.

Revisiting the Great American Biotic Interchange through analyses of amphitropical bees

The Great American Biotic Interchange (GABI) is zoogeographic event characterized by the exchange of taxa between North and South America, typically associated with the rise of the Isthmus of Panama in the late Pliocene. Recent geologic evidence suggests the connections between North and South America may be much older, and that the interchange of organisms between the two continents could have therefore happened much earlier than 3 Ma. Most of the research investigating the GABI has come from tropical vertebrate taxa; little work has been done on invertebrates or on non‐tropical species. To investigate how the GABI shaped the distribution of arid‐adapted species, particularly those with amphitropical distributions (i.e. taxa found in South and North American xeric regions yet absent from the tropics), we examine the historical biogeography of the bee genus Diadasia using a hypothesis of Diadasia phylogenetic relationships. Nuclear and mitochondrial genetic loci are used to reconstruct a phylogeny of Diadasia, which is then used to estimate divergence dates and reconstruct ancestral area relationships. Our analyses suggest the divergence between North and South American Diadasia species occurred between 20.5 and 15 Ma, long before the formation of the Isthmus of Panama. This study is the first to show a Miocene connection for an amphitropically‐distributed insect group. It suggests that the biotic connection between continents is more complicated than previously thought and may have initiated long before the late Pliocene.

Assignment of Calibration Information to Deeper Phylogenetic Nodes is More Effective in Obtaining Precise and Accurate Divergence Time Estimates.

Divergence time estimation has become an essential tool for understanding macroevolutionary events. Molecular dating aims to obtain reliable inferences, which, within a statistical framework, means jointly increasing the accuracy and precision of estimates. Bayesian dating methods exhibit the propriety of a linear relationship between uncertainty and estimated divergence dates. This relationship occurs even if the number of sites approaches infinity and places a limit on the maximum precision of node ages. However, how the placement of calibration information may affect the precision of divergence time estimates remains an open question. In this study, relying on simulated and empirical data, we investigated how the location of calibration within a phylogeny affects the accuracy and precision of time estimates. We found that calibration priors set at median and deep phylogenetic nodes were associated with higher precision values compared to analyses involving calibration at the shallowest node. The results were independent of the tree symmetry. An empirical mammalian dataset produced results that were consistent with those generated by the simulated sequences. Assigning time information to the deeper nodes of a tree is crucial to guarantee the accuracy and precision of divergence times. This finding highlights the importance of the appropriate choice of outgroups in molecular dating.

Molecular phylogeny of the Notostraca

We used a combined analysis of one nuclear (28S rDNA) and three mitochondrial markers (COI, 12S rDNA, 16S rDNA) to infer the molecular phylogeny of the Notostraca, represented by samples from the six continents that are inhabited by this group of branchiopod crustaceans. Our results confirm the monophyly of both extant notostracan genera Triops and Lepidurus with good support in model based and maximum parsimony analyses. We used branchiopod fossils as a calibration to infer divergence times among notostracan lineages and accounted for rate heterogeneity among lineages by applying relaxed-clock models. Our divergence date estimates indicate an initial diversification into the genera Triops and Lepidurus in the Mesozoic, most likely at a minimum age of 152.3–233.5Ma, i.e., in the Triassic or Jurassic. Implications for the interpretation of fossils and the evolution of notostracan morphology are discussed. We further use the divergence date estimates to formulate a biogeographic hypothesis that explains distributions of extant lineages predominantly by overland dispersal routes.We identified an additional hitherto unrecognised highly diverged lineage within Lepidurus apus lubbocki and three additional previously unknown major lineages within Triops. Within T. granarius we found deep differentiation, with representatives distributed among three major phylogenetic lineages. One of these major lineages comprises T. cancriformis, the T. mauritanicus species group and two hitherto unrecognised T. granarius lineages. Samples that were morphologically identified as T. granarius diverged from the most basal nodes within this major lineage, and divergence dates suggested an approximate age of 23.7–49.6Ma for T. cancriformis, indicating the need for a taxonomic revision of Triassic and Permian fossils that are currently attributed to the extant T. cancriformis. We thus elevate T. cancriformis minor to full species status as Triops minorTrusheim, 1938 and include in this species the additional Upper Triassic samples that were attributed to T. cancriformis. We further elevate T. cancriformis permiensis to full species status as Triops permiensisGand et al., 1997.

Did southern Western Ghats of peninsular India serve as refugia for its endemic biota during the Cretaceous volcanism?

The Western Ghats (WG) of south India, a global biodiversity hotspot, has experienced complex geological history being part of Gondwana landmass and encountered extensive volcanic activity at the end of Cretaceous epoch. It also has a climatically and topographically heterogeneous landscape. Thus, the WG offer a unique setting to explore the influence of ecological and geological processes on the current diversity and distribution of its biota. To this end, three explicit biogeographical scenarios were hypothesized to evaluate the distribution and diversification of wet evergreen species of the WG – (1) southern WG was a refuge for the wet evergreen species during the Cretaceous volcanism, (2) phylogenetic breaks in the species phylogeny would correspond to geographic breaks (i.e., the Palghat gap) in the WG, and (3) species from each of the biogeographic subdivisions within the WG would form distinct clades. These hypotheses were tested on the centipede genus Digitipes from the WG which is known to be an ancient, endemic, and monophyletic group. The Digitipes molecular phylogeny was subjected to divergence date estimation using Bayesian approach, and ancestral areas were reconstructed using parsimony approach for each node in the phylogeny. Ancestral-area reconstruction suggested 13 independent dispersal events to explain the current distribution of the Digitipes species in the WG. Among these 13 dispersals, two dispersal events were at higher level in the Digitipes phylogeny and were from the southern WG to the central and northern WG independently in the Early Paleocene, after the Cretaceous Volcanism. The remaining 11 dispersal events explained the species’ range expansions of which nine dispersals were from the southern WG to other biogeographic subdivisions in the Eocene-Miocene in the post-volcanic periods where species-level diversifications occurred. Taken together, these results suggest that southern WG might have served as a refuge for Digitipes species during Cretaceous volcanism.

Phylogeographic analyses strongly suggest cryptic speciation in the giant spiny frog (Dicroglossidae: Paa spinosa) and interspecies hybridization in Paa.

Species identification is one of the most basic yet crucial issues in biology with potentially far-reaching implications for fields such as conservation, population ecology, and epidemiology. The widely distributed but threatened frog Paa spinosa has been speculated to represent a complex of multiple species. In this study, 254 individuals representing species of the genus Paa were investigated along the entire range of P. spinosa: 196 specimens of P. spinosa, 8 specimens of P. jiulongensis, 5 specimens of P. boulengeri, 20 specimens of P. exilispinosa, and 25 specimens of P. shini. Approximately 1333 bp of mtDNA sequence data (genes 12S rRNA and 16S rRNA) were used. Phylogenetic analyses were conducted using maximum parsimony, maximum likelihood and Bayesian inference. BEAST was used to estimate divergence dates of major clades. Results suggest that P. spinosa can be divided into three distinct major lineages. Each major lineage totally corresponds to geographical regions, revealing the presence of three candidate cryptic species. Isolation and differentiation among lineages are further supported by the great genetic distances between the lineages. The bifurcating phylogenetic pattern also suggests an east–west dispersal trend during historic cryptic speciation. Dating analysis estimates that P. spinosa from Western China split from the remaining P. spinosa populations in the Miocene and that P. spinosa from Eastern China diverged from Central China in the Pliocene. We also found that P. exilispinosa from Mainland China and Hong Kong might have a complex of multiple species. After identifying cryptic lineages, we then determine the discrepancy between the mtDNA and the morphotypes in several individuals. This discrepancy may have been caused by introgressive hybridization between P. spinosa and P. shini.

Phylogeography of two closely related species of Nolana from the coastal Atacama Desert of Chile: post‐glacial population expansions in response to climate fluctuations

AbstractAimTo investigate the impact of Quaternary climate cycles on the coastal Atacama Desert flora by assessing phylogeographical patterns of the desert shrub Nolana crassulifolia (Solanaceae) and its congener Nolana incana.LocationThe latitudinal aridity gradient from the southern margin of the coastal Atacama Desert to the mediterranean semi‐arid region of Chile (25–33° S).MethodsTwo cpDNA regions were sequenced for 130 individuals in 15 populations of these two closely related species, covering their entire distribution range (1000 km) along the arid to semi‐arid Chilean coast. We explored haplotype relationships in a statistical parsimony network, and assessed population genetic diversity, population differentiation and phylogeographical structure. In addition, we conducted demographic analyses and spatial analysis of genetic variation, identified barriers to gene flow with Monmonier's algorithm, and used Bayesian phylogenetic reconstruction to estimate divergence dates between lineages.ResultsWe found a total of 14 haplotypes – four of them shared by both species – and high levels of genetic differentiation among populations, but no past distribution breaks that could account for major vicariant events. Genetic diversity decreased continuously from north to south, with loss of haplotypes and a greater number of monomorphic populations in the southern range. Landscape analysis revealed greater genetic differentiation in the northernmost populations of both species.Main conclusionsThe documented north–south gradient of declining genetic diversity, the origin and location of ancestral haplotypes in northern sites, and the loss of haplotypes from southern populations all support the hypothesis of post‐glacial range expansion of Nolana southwards during arid/warmer cycles in the Atacama. The higher genetic diversity and greater differentiation of northern populations of both Nolana species support the hypothesis that populations survived in northern arid sites during wetter/colder episodes of the glacial cycles. We suggest that Quaternary arid phases in the Atacama promoted southward expansion of the coastal desert vegetation.

A molecular phylogeny of nephilid spiders: Evolutionary history of a model lineage

The pantropical orb web spider family Nephilidae is known for the most extreme sexual size dimorphism among terrestrial animals. Numerous studies have made Nephilidae, particularly Nephila, a model lineage in evolutionary research. However, a poorly understood phylogeny of this lineage, relying only on morphology, has prevented thorough evolutionary syntheses of nephilid biology. We here use three nuclear and five mitochondrial genes for 28 out of 40 nephilid species to provide a more robust nephilid phylogeny and infer clade ages in a fossil-calibrated Bayesian framework. We complement the molecular analyses with total evidence analysis including morphology. All analyses find strong support for nephilid monophyly and exclusivity and the monophyly of the genera Herennia and Clitaetra. The inferred phylogenetic structure within Nephilidae is novel and conflicts with morphological phylogeny and traditional taxonomy. Nephilengys species fall into two clades, one with Australasian species (true Nephilengys) as sister to Herennia, and another with Afrotropical species (Nephilingis Kuntner new genus) as sister to a clade containing Clitaetra plus most currently described Nephila. Surprisingly, Nephila is also diphyletic, with true Nephila containing N. pilipes+N. constricta, and the second clade with all other species sister to Clitaetra; this “Nephila” clade is further split into an Australasian clade that also contains the South American N. sexpunctata and the Eurasian N. clavata, and an African clade that also contains the Panamerican N. clavipes. An approximately unbiased test constraining the monophyly of Nephilengys, Nephila, and Nephilinae (Nephila, Nephilengys, Herennia), respectively, rejected Nephilengys monophyly, but not that of Nephila and Nephilinae. Further data are therefore necessary to robustly test these two new, but inconclusive findings, and also to further test the precise placement of Nephilidae within the Araneoidea. For divergence date estimation we set the minimum bound for the stems of Nephilidae at 40Ma and of Nephila at 16Ma to accommodate Palaeonephila from Baltic amber and Dominican Nephila species, respectively. We also calibrated and dated the phylogeny under three different interpretations of the enigmatic 165Ma fossil Nephila jurassica, which we suspected based on morphology to be misplaced. We found that by treating N. jurassica as stem Nephila or nephilid the inferred clade ages were vastly older, and the mitochondrial substitution rates much slower than expected from other empirical spider data. This suggests that N. jurassica is not a Nephila nor a nephilid, but possibly a stem orbicularian. The estimated nephilid ancestral age (40–60Ma) rejects a Gondwanan origin of the family as most of the southern continents were already split at that time. The origin of the family is equally likely to be African, Asian, or Australasian, with a global biogeographic history dominated by dispersal events. A reinterpretation of web architecture evolution suggests that a partially arboricolous, asymmetric orb web with a retreat, as exemplified by both groups of “Nephilengys”, is plesiomorphic in Nephilidae, that this architecture was modified into specialized arboricolous webs in Herennia and independently in Clitaetra, and that the web became aerial, gigantic, and golden independently in both “Nephila” groups. The new topology questions previously hypothesized gradual evolution of female size from small to large, and rather suggests a more mosaic evolutionary pattern with independent female size increases from medium to giant in both “Nephila” clades, and two reversals back to medium and small; combined with male size evolution, this pattern will help detect gross evolutionary events leading to extreme sexual size dimorphism, and its morphological and behavioral correlates.

Ancient host shifts followed by host conservatism in a group of ant parasitoids

While ant colonies serve as host to a diverse array of myrmecophiles, few parasitoids are able to exploit this vast resource. A notable exception is the wasp family Eucharitidae, which is the only family of insects known to exclusively parasitize ants. Worldwide, approximately 700 Eucharitidae species attack five subfamilies across the ant phylogeny. Our goal is to uncover the pattern of eucharitid diversification, including timing of key evolutionary events, biogeographic patterns and potential cophylogeny with ant hosts. We present the most comprehensive molecular phylogeny of Eucharitidae to date, including 44 of the 53 genera and fossil-calibrated estimates of divergence dates. Eucharitidae arose approximately 50 Ma after their hosts, during the time when the major ant lineages were already established and diversifying. We incorporate host association data to test for congruence between eucharitid and ant phylogenies and find that their evolutionary histories are more similar than expected at random. After a series of initial host shifts, clades within Eucharitidae maintained their host affinity. Even after multiple dispersal events to the New World and extensive speciation within biogeographic regions, eucharitids remain parasitic on the same ant subfamilies as their Old World relatives, suggesting host conservatism despite access to a diverse novel ant fauna.