What causes latitudinal gradients in species diversity? Evolutionary processes and ecological constraints on swallowtail biodiversity

Tropical Origin, Global Diversification, and Dispersal in the Pond Damselflies (Coenagrionoidea) Revealed by a New Molecular Phylogeny

The processes responsible for the formation of Earth's most conspicuous diversity pattern, the latitudinal diversity gradient (LDG), remain unexplored for many clades in the Tree of Life. Here, we present a densely sampled and dated molecular phylogeny for the most speciose clade of damselflies worldwide (Odonata: Coenagrionoidea) and investigate the role of time, macroevolutionary processes, and biome-shift dynamics in shaping the LDG in this ancient insect superfamily. We used process-based biogeographic models to jointly infer ancestral ranges and speciation times and to characterize within-biome dispersal and biome-shift dynamics across the cosmopolitan distribution of Coenagrionoidea. We also investigated temporal and biome-dependent variation in diversification rates. Our results uncover a tropical origin of pond damselflies and featherlegs ~105 Ma, while highlighting the uncertainty of ancestral ranges within the tropics in deep time. Even though diversification rates have declined since the origin of this clade, global climate change and biome-shifts have slowly increased diversity in warm- and cold-temperate areas, where lineage turnover rates have been relatively higher. This study underscores the importance of biogeographic origin and time to diversify as important drivers of the LDG in pond damselflies and their relatives, while diversification dynamics have instead resulted in the formation of ephemeral species in temperate regions. Biome-shifts, although limited by tropical niche conservatism, have been the main factor reducing the steepness of the LDG in the last 30 Myr. With ongoing climate change and increasing northward range expansions of many damselfly taxa, the LDG may become less pronounced. Our results support recent calls to unify biogeographic and macroevolutionary approaches to improve our understanding of how latitudinal diversity gradients are formed and why they vary across time and among taxa.

Responses of marine trophic levels to the combined effects of ocean acidification and warming

Marine organisms are simultaneously exposed to anthropogenic stressors associated with ocean acidification and ocean warming, with expected interactive effects. Species from different trophic levels with dissimilar characteristics and evolutionary histories are likely to respond differently. Here, we perform a meta-analysis of controlled experiments including both ocean acidification and ocean warming factors to investigate single and interactive effects of these stressors on marine species. Contrary to expectations, we find that synergistic interactions are less common (16%) than additive (40%) and antagonistic (44%) interactions overall and their proportion decreases with increasing trophic level. Predators are the most tolerant trophic level to both individual and combined effects. For interactive effects, calcifying and non-calcifying species show similar patterns. We also identify climate region-specific patterns, with interactive effects ranging from synergistic in temperate regions to compensatory in subtropical regions, to positive in tropical regions. Our findings improve understanding of how ocean warming, and acidification affect marine trophic levels and highlight the need for deeper consideration of multiple stressors in conservation efforts.

Comprehensive phylogeny of Pieridae butterflies reveals strong correlation between diversification and temperature

Temperature is thought to be a key factor influencing global species richness patterns. We investigate the link between temperature and diversification in the butterfly family Pieridae by combining next generation DNA sequences and published molecular data with fine-grained distribution data. We sampled nearly 600 pierid butterfly species to infer the most comprehensive molecular phylogeny of the family and curated a distribution dataset of more than 800,000 occurrences. We found strong evidence that species in environments with more stable daily temperatures or cooler maximum temperatures in the warm seasons have higher speciation rates. Furthermore, speciation and extinction rates decreased in tandem with global temperatures through geological time, resulting in a constant net diversification.

Effects of climate and environmental heterogeneity on the phylogenetic structure of regional angiosperm floras worldwide

The tendency of species to retain ancestral ecological distributions (phylogenetic niche conservatism) is thought to influence which species from a species pool can persist in a particular environment. Thus, investigating the relationships between measures of phylogenetic structure and environmental variables at a global scale can help understand the variation in species richness and phylogenetic structure in biological assemblages across the world. Here, we analyze a comprehensive data set including 341,846 species in 391 angiosperm floras worldwide to explore the relationships between measures of phylogenetic structure and environmental variables for angiosperms in regional floras across the world and for each of individual continental (biogeographic) regions. We find that the global phylogenetic structure of angiosperms shows clear and meaningful relationships with environmental factors. Current climatic variables have the highest predictive power, especially on phylogenetic metrics reflecting recent evolutionary relationships that are also related to current environmental heterogeneity, presumably because this favors plant speciation in various ways. We also find evidence that past climatic conditions, and particularly refugial conditions, play an important role in determining the phylogenetic structure of regional floras. The relationships between environmental conditions and phylogenetic metrics differ between continents, reflecting the different evolutionary histories of their floras.

Taxonomic Structure and Wing Pattern Evolution in the Parnassius mnemosyne Species Complex (Lepidoptera, Papilionidae)

In our study, using the analysis of DNA barcodes and morphology (wing color, male genitalia, and female sphragis shape), we show that the group of species close to P. mnemosyne comprises the western and eastern phylogenetic lineages. The eastern lineage includes P. stubbendorfii, P. glacialis, and P. hoenei. The western lineage includes three morphologically similar species: P. mnemosyne (Western Eurasia), P. turatii (southwestern Europe), and P. nubilosusstat. nov. (Turkmenistan and NE Iran), as well as the morphologically differentiated P. ariadne (Altai). The latter species differs from the rest of the group in the presence of red spots on the wings. Parnassius mnemosyne s.s. is represented by four differentiated mitochondrial clusters that show clear association with specific geographic regions. We propose to interpret them as subspecies: P. mnemosyne mnemosyne (Central and Eastern Europe, N Caucasus, N Turkey), P. mnemosyne adolphi (the Middle East), P. mnemosyne falsa (Tian Shan), and P. mnemosyne gigantea (Gissar-Alai in Central Asia). We demonstrate that in P. ariadne, the red spots on the wing evolved as a reversion to the ancestral wing pattern. This reversion is observed in Altai, where the distribution areas of the western lineage, represented by P. ariadne, and the eastern lineage, represented by P. stubbendorfii, overlap. These two species hybridize in Altai, and we hypothesize that the color change in P. ariadne is the result of reinforcement of prezygotic isolation in the contact zone. The lectotype of Parnassius mnemosyne var. nubilosus Christoph, 1873, is designated.

Genomics of the relict species Baronia brevicornis sheds light on its demographic history and genome size evolution across swallowtail butterflies

Relict species, like coelacanth, gingko, tuatara, are the remnants of formerly more ecologically and taxonomically diverse lineages. It raises the questions of why they are currently species-poor, have restrained ecology, and are often vulnerable to extinction. Estimating heterozygosity level and demographic history can guide our understanding of the evolutionary history and conservation status of relict species. However, few studies have focused on relict invertebrates compared to vertebrates. We sequenced the genome of Baronia brevicornis (Lepidoptera: Papilionidae), which is an endangered species, the sister species of all swallowtail butterflies, and is the oldest lineage of all extant butterflies. From a dried specimen, we were able to generate both long-read and short-read data and assembled a genome of 406 Mb for Baronia. We found a fairly high level of heterozygosity (0.58%) compared to other swallowtail butterflies, which contrasts with its endangered and relict status. Taking into account the high ratio of recombination over mutation, demographic analyses indicated a sharp decline of the effective population size initiated in the last million years. Moreover, the Baronia genome was used to study genome size variation in Papilionidae. Genome sizes are mostly explained by transposable elements activities, suggesting that large genomes appear to be a derived feature in swallowtail butterflies as transposable elements activity is recent and involves different transposable elements classes among species. This first Baronia genome provides a resource for assisting conservation in a flagship and relict insect species as well as for understanding swallowtail genome evolution.

Global determinants of insect mitochondrial genetic diversity

Understanding global patterns of genetic diversity is essential for describing, monitoring, and preserving life on Earth. To date, efforts to map macrogenetic patterns have been restricted to vertebrates, which comprise only a small fraction of Earth's biodiversity. Here, we construct a global map of predicted insect mitochondrial genetic diversity from cytochrome c oxidase subunit 1 sequences, derived from open data. We calculate the mitochondrial genetic diversity mean and genetic diversity evenness of insect assemblages across the globe, identify their environmental correlates, and make predictions of mitochondrial genetic diversity levels in unsampled areas based on environmental data. Using a large single-locus genetic dataset of over 2 million globally distributed and georeferenced mtDNA sequences, we find that mitochondrial genetic diversity evenness follows a quadratic latitudinal gradient peaking in the subtropics. Both mitochondrial genetic diversity mean and evenness positively correlate with seasonally hot temperatures, as well as climate stability since the last glacial maximum. Our models explain 27.9% and 24.0% of the observed variation in mitochondrial genetic diversity mean and evenness in insects, respectively, making an important step towards understanding global biodiversity patterns in the most diverse animal taxon.

Up and down from North to South: Latitudinal Distribution of Flea Beetle Genera in Continental Africa (Coleoptera, Chrysomelidae, Galerucinae, Alticini)

The distribution of global biodiversity can be investigated based on comprehensive datasets and many methods to process them. The taxonomic diversity of phytophagous insects is typically linked to plant diversity, which increases from temperate to tropical latitudes. In this paper, we explored the latitudinal distribution of the flea beetle genera (Coleoptera, Chrysomelidae, Galerucinae, Alticini) on the African continent. We divided the area into latitudinal belts and looked for possible correlations with the number and types of vegetational divisions, the area of each belt, and the bioclimatic variables. The number of flea beetle genera is related to the number and types of vegetation divisions rather than the area of each belt. Some bioclimatic variables are highly related to the number of genera, which is higher within those belts where climate factors limit the oscillation of temperature over the year and favor high precipitations, especially in the warmest months. These biotic and abiotic factors lead to a two-peak trend in the taxonomic richness of flea beetle genera from north to south. Genera endemic to restricted areas are linked to the presence of high mountain systems and increase the taxonomic richness of the belt they belong to.

A comprehensive phylogeny and revised taxonomy illuminate the origin and diversification of the global radiation of Papilio (Lepidoptera: Papilionidae)

The swallowtail genus Papilio (Lepidoptera: Papilionidae) is species rich, distributed worldwide, and has broad morphological habits and ecological niches. Because of its elevated species richness, it has been historically difficult to reconstruct a densely sampled phylogeny for this clade. Here we provide a taxonomic working list for the genus, resulting in 235 Papilio species, and assemble a molecular dataset of seven gene fragments representing ca. 80% of the currently described diversity. Phylogenetic analyses reconstructed a robust tree with highly supported relationships within subgenera, although a few nodes in the early history of the Old World Papilio remain unresolved. Contrasting with previous results, we found that Papilio alexanor is sister to all Old World Papilio and that the subgenus Eleppone is no longer monotypic. The latter includes the recently described Fijian Papilio natewa with the Australian Papilio anactus and is sister to subgenus Araminta (formerly included in subgenus Menelaides) occurring in Southeast Asia. Our phylogeny also includes rarely studied (P. antimachus, P. benguetana) or endangered species (P. buddha, P. chikae). Taxonomic changes resulting from this study are elucidated. Molecular dating and biogeographic analyses indicate that Papilio originated ca. 30 million years ago (Oligocene), in a northern region centered on Beringia. A rapid early Miocene radiation in the Paleotropics is revealed within Old World Papilio, potentially explaining their low early branch support. Most subgenera originated in the early to middle Miocene followed by synchronous southward biogeographic dispersals and repeated local extirpations in northern latitudes. This study provides a comprehensive phylogenetic framework for Papilio with clarification of subgeneric systematics and species taxonomic changes enumerated, which will facilitate further studies to address questions on their ecology and evolutionary biology using this model clade.

Species Richness of Papilionidae Butterflies (Lepidoptera: Papilionoidea) in the Hengduan Mountains and Its Future Shifts under Climate Change

The family of Papilionidae (Lepidoptera: Papilionoidea) is a group of butterflies with high ecological and conservation value. The Hengduan Mountains (HMDs) in Southwest China is an important diversity centre for these butterflies. However, the spatial distribution pattern and the climate vulnerability of Papilionidae butterflies in the HDMs remain unknown to date. The lack of such knowledge has already become an obstacle in formulating effective butterfly conservation strategies. The present research compiled a 59-species dataset with 1938 occurrence points. The Maxent model was applied to analyse the spatial pattern of species richness in subfamilies Parnassiinae and Papilioninae, as well as to predict the response under the influence of climate change. The spatial pattern of both subfamilies in the HDMs has obvious elevation prevalence, with Parnassiinae concentrated in the subalpine to alpine areas (2500-5500 m) in western Sichuan, northwestern Yunnan and eastern Tibet, while Papilioninae is concentrated in the low- to medium-elevation areas (1500-3500 m) in the river valleys of western Yunnan and western Sichuan. Under the influence of climate change, both subfamilies would exhibit northward and upward range shifts. The majority of Parnassiinae species would experience drastic habitat contraction, resulting in lower species richness across the HDMs. In contrast, most Papilioninae species would experience habitat expansion, and the species richness would also increase significantly. The findings of this research should provide new insights and a clue for butterfly diversity and climatic vulnerability in southwestern China. Future conservation efforts should be focused on species with habitat contraction, narrow-ranged distribution and endemicity with both in situ and ex situ measures, especially in protected areas. Commercialised collecting targeting these species must also be regulated by future legislation.

Phylogenetic diversity of eastern Asia-eastern North America disjunct plants is mainly associated with divergence time

The underlying causes of biodiversity disparities among geographic regions have long been a fundamental theme in ecology and evolution. However, the patterns of phylogenetic diversity (PD) and phylogenetic beta diversity (PBD) of congeners that are disjunctly distributed between eastern Asia-eastern North America (EA-ENA disjuncts) and their associated factors remain unknown. Here we investigated the standardized effect size of PD (SES-PD), PBD, and potentially associated factors in 11 natural mixed forest sites (five in EA and six in ENA) where abundant EA-ENA disjuncts occur. We found that the disjuncts in ENA possessed higher SES-PD than those in EA at the continental scale (1.96 vs -1.12), even though the number of disjunct species in ENA is much lower than in EA (128 vs 263). SES-PD of the EA-ENA disjuncts tended to decrease with increasing latitude in 11 sites. The latitudinal diversity gradient of SES-PD was stronger in EA sites than in ENA sites. Based on the unweighted unique fraction metric (UniFrac) distance and the phylogenetic community dissimilarity, PBD showed that the two northern sites in EA were more similar to the six-site ENA group than to the remaining southern EA sites. Based on the standardized effect size of mean pairwise distances (SES-MPD), nine of eleven studied sites showed a neutral community structure (-1.96 ≤ SES-MPD ≤ 1.96). Both Pearson's r and structural equation modeling suggested that SES-PD of the EA-ENA disjuncts was mostly associated with mean divergence time. Moreover, SES-PD of the EA-ENA disjuncts was positively correlated with temperature-related climatic factors, although negatively correlated with mean diversification rate and community structure. By applying approaches from phylogenetics and community ecology, our work sheds light on historical patterns of the EA-ENA disjunction and paves the way for further research.

Global Patterns of Earwig Species Richness

One of the most investigated patterns in species diversity is the so-called latitudinal gradient, that is, a decrease in species richness from the equator to the poles. However, few studies investigated this pattern in insects at a global scale because of insufficient taxonomic and biogeographical information. Using estimates of earwig species richness at country level, their latitudinal diversity gradient was modelled globally and for the two hemispheres separately after correcting for differences in country areas. Separate analyses were also conducted for mainland and island countries. All analyses clearly indicated the existence of latitudinal gradients. The most plausible explanation for the observed pattern is the so-called tropical conservatism hypothesis, which postulates (1) a tropical origin of many extant clades, (2) a longer time for cladogenesis in tropical environments thanks to their environmental stability, and (3) a limited ability of historically tropical lineages to adapt to temperate climates. Earwigs probably evolved on Gondwana and secondarily colonized the Northern Hemisphere. This colonization was hampered by both geographical and climatic factors. The Himalayan orogenesis obstructed earwig dispersal into the Palearctic region. Additionally, earwig preferences for warm/hot and humid climates hampered the colonization of temperate regions. Pleistocene glaciation further contributed to reducing diversity at northern latitudes.

A phylogenetic study to assess the link between biome specialization and diversification in swallowtail butterflies

The resource-use hypothesis, proposed by E.S. Vrba, states that habitat fragmentation caused by climatic oscillations would affect particularly biome specialists (species inhabiting only one biome), which might show higher speciation and extinction rates than biome generalists. If true, lineages would accumulate biome-specialist species. This effect would be particularly exacerbated for biomes located at the periphery of the global climatic conditions, namely, biomes that have high/low precipitation and high/low temperature such as rainforest (warm-humid), desert (warm-dry), steppe (cold-dry) and tundra (cold-humid). Here, we test these hypotheses in swallowtail butterflies, a clade with more than 570 species, covering all the continents but Antarctica, and all climatic conditions. Swallowtail butterflies are among the most studied insects, and they are a model group for evolutionary biology and ecology studies. Continental macroecological rules are normally tested using vertebrates, this means that there are fewer examples exploring terrestrial invertebrate patterns at global scale. Here, we compiled a large Geographic Information System database on swallowtail butterflies' distribution maps and used the most complete time-calibrated phylogeny to quantify diversification rates (DRs). In this paper, we aim to answer the following questions: (1) Are there more biome-specialist swallowtail butterflies than biome generalists? (2) Is DR related to biome specialization? (3) If so, do swallowtail butterflies inhabiting extreme biomes show higher DRs? (4) What is the effect of species distribution area? Our results showed that swallowtail family presents a great number of biome specialists which showed substantially higher DRs compared to generalists. We also found that biome specialists are unevenly distributed across biomes. Overall, our results are consistent with the resource-use hypothesis, species climatic niche and biome fragmentation as key factors promoting isolation.

Genome assembly of Luehdorfia taibai, an endangered butterfly endemic to Qinling Moutains in China with extremely small populations

Conservation genomic resources over the past decade has drastically improved, since genomes can be used to predict diverse parameters vital to conservation management. Luehdorfia taibai is an endemic butterfly only found in restricted aeras in middle-west China and is critically endangered. It was classfied as a vunerlable (VN) species in the “China species red list.” Here we generated 34.38 Gb of raw DNA sequencing reads and obtained a high-qualified draft genome assembly of L. taibai. The final genome is ~683.3 Mb, with contig N50 size of 10.19 Mb. Further, 98.6% of single-copy orthologous genes have been recovered by BUSCO. An estimated 42.34% of the genome of L. taibai consists of repetitive elements. Combined with gene prediction and transcriptome sequencing, genome annotation produced 15,968 protein-coding genes. Additionally, a nearly 1:1 orthology ratio of syntenic blocks between L. taibai and its closest genome Luehdorfia chinensis suggested that the genome structures have not changed much after speciation. The genome of L. taibai have not undergone a whole genome duplication event. Population dynamics analyses indicates that L. taibai has an extremely low heterozygosity of 0.057%, and its population size has declined dramatically over the past 10 thousand years. Our study describes a draft genome assembly of the L. taibai, the first implication of this species. We consider the globally overexploited of the host plants is not the main reason to threaten L. taibai. The genome will provide advice for the conservation to the economically important Luehdorfia lineage and this specific species.

From monocots to dicots: host shifts in Afrotropical derelomine weevils shed light on the evolution of non-obligatory brood pollination mutualism

Weevils from the tribe Derelomini (Curculionidae: Curculioninae) are specialized brood pollinators engaged in mutualistic relationships with several angiosperm lineages. In brood pollination systems, reproductive plant tissues are used for the development of insect larval stages, whereas adult insects pollinate their plant hosts as a reward. The evolutionary history of derelomines in relationship to their hosts is poorly understood and potentially contrasts with other brood pollination systems, wherein a pollinator lineage is usually associated with a single host plant family. In the case of Afrotropical Derelomini, host records indicate a diverse host repertoire consisting of several families of monocot and dicot plants. In this study, we investigate their phylogenetic relationships, timing of diversification and evolution of host use. Our results suggest that derelomine lineages started their diversification ~40 Mya. Reconstructions of host use evolution support an ancestral association with the monocotyledonous palm family (Arecaceae), followed by several shifts towards other plant families in Afrotropical lineages, especially to dicotyledonous plants from the family Ebenaceae (on the genus Euclea L.). Some level of phylogenetic conservatism of host use is recovered for the lineages associated with either palms or Euclea. Multiple instances of sympatric weevil assemblages on the same plant are also unravelled, corresponding to either single or independent colonization events. Overall, the diversity of hosts colonized and the frequency of sympatric assemblages highlighted in non-obligatory plant–derelomine brood pollination systems contrast with what is generally expected from plant–insect brood pollination systems.

Conserved ancestral tropical niche but different continental histories explain the latitudinal diversity gradient in brush-footed butterflies

The global increase in species richness toward the tropics across continents and taxonomic groups, referred to as the latitudinal diversity gradient, stimulated the formulation of many hypotheses to explain the underlying mechanisms of this pattern. We evaluate several of these hypotheses to explain spatial diversity patterns in a butterfly family, the Nymphalidae, by assessing the contributions of speciation, extinction, and dispersal, and also the extent to which these processes differ among regions at the same latitude. We generate a time-calibrated phylogeny containing 2,866 nymphalid species (~45% of extant diversity). Neither speciation nor extinction rate variations consistently explain the latitudinal diversity gradient among regions because temporal diversification dynamics differ greatly across longitude. The Neotropical diversity results from low extinction rates, not high speciation rates, and biotic interchanges with other regions are rare. Southeast Asia is also characterized by a low speciation rate but, unlike the Neotropics, is the main source of dispersal events through time. Our results suggest that global climate change throughout the Cenozoic, combined with tropical niche conservatism, played a major role in generating the modern latitudinal diversity gradient of nymphalid butterflies.

Geographic Pattern of Bryophyte Species Richness in China: The Influence of Environment and Evolutionary History

How contemporary environment interacts with macroevolutionary processes to generate the geographic pattern of bryophyte species is still unresolved. China is very rich in bryophytes, with more than 3,000 bryophytes covering 70% of the families in the world. In this study, we assessed the effects of the contemporary environment (average temperature of the coldest season TCQ, precipitation of the warmest season PWQ, and elevational range) and the recent diversification rates (estimated as mean species number per genus, MSG) on the geographical pattern of species richness for bryophytes and two groups (i.e., liverworts and mosses) in China. We compiled the provincial level distribution of bryophyte species and estimated the geographic pattern of the recent diversification rate by MSG for species in China. Univariate, multivariate regressions and path model analyses were used to assess the relationships between species richness, MSG, and their potential environmental drivers. Species richness of all bryophytes and liverworts significantly increased with the increase of MSG, either in regressions or path analyses, indicating that provinces with high bryophyte richness were mainly inhabited by species (especially liverworts) from lineages with particularly high MSG. In contrast, the species richness of mosses was insignificantly decreased with MSG in univariate regression or insignificantly increased with MSG in path analysis. Both species richness and MSG of all bryophytes and liverworts increased with the increase in energy and water availability. In contrast, for mosses, the species richness significantly increased with the increase of energy and water availability, while MSG decreased with the increase of energy and water availability. The MSG of liverworts increase with the increase of elevational range but the MSG of mosses decrease with the increase of elevational range. Our study suggests that the humid tropical and subtropical mountains in China are not only diversity hotspots for bryophytes, but also cradles for high recent diversification of liverworts, and refuges for mosses to hold many monotypic and oligotypic genera.

Out of the extratropics: the evolution of the latitudinal diversity gradient of Cenozoic marine plankton

Many ecological and evolutionary hypotheses have been proposed to explain the latitudinal diversity gradient, i.e. the increase in species richness from the poles to the tropics. Among the evolutionary hypotheses, the 'out of the tropics' (OTT) hypothesis has received considerable attention. The OTT posits that the tropics are both a cradle and source of biodiversity for extratropical regions. To test the generality of the OTT hypothesis, we explored the spatial biodiversity dynamics of unicellular marine plankton over the Cenozoic era (the last 66 Myr). We find large-scale climatic changes during the Cenozoic shaped the diversification and dispersal of marine plankton. Origination was generally more likely in the extratropics and net dispersal was towards the tropics rather than in the opposite direction, especially during the warmer climates of the early Cenozoic. Although migration proportions varied among major plankton groups and climate phases, we provide evidence that the extratropics were a source of tropical microplankton biodiversity over the last 66 Myr.

Phylogenetic and Functional Traits Verify the Combined Effect of Deterministic and Stochastic Processes in the Community Assembly of Temperate Forests along an Elevational Gradient

Explaining community assembly mechanisms along elevational gradients dominated by deterministic processes or stochastic processes is a pressing challenge. Many studies suggest that phylogenetic and functional diversity are significant indicators of the process. In this study, we analyzed the structure and beta diversity of phylogenetic and functional traits along an elevational gradient and discussed the effects of environmental and spatial factors. We found that the phylogenetic and functional traits showed inconsistent changes, and their variations were closely related to the abiotic environment. The results suggested that the community assembly of woody plants was obviously affected by the combined effect of deterministic processes and the stochastic hypothesis (primarily by the latter). Phylogenetic and functional traits had a certain relationship but changed according to different rules. These results enhance our understanding of the assembly mechanism of forest communities by considering both phylogenetic and functional traits.

Backbone phylogeny and evolution of Apioideae (Apiaceae): New insights from phylogenomic analyses of plastome data

Traditional phylogenies inferred from chloroplast DNA fragments have not obtained a well-resolved evolutionary history for the backbone of Apioideae, the largest subfamily of Apiaceae. In this study, we applied the genome skimming approach of next-generation sequencing to address whether the lack of resolution at the tip of the Apioideae phylogenetic tree is due to limited information loci or the footprint of ancient radiation. A total of 90 complete chloroplast genomes (including 23 newly sequenced genomes and covering 20 major clades of Apioideae) were analyzed (RAxML and MrBayes) to provide a phylogenomic reconstruction of Apioideae. Dating analysis was also implemented using BEAST to estimate the origin and divergence time of the major clades. As a result, the early divergences of Apioideae have been clarified but the relationship among its distally branching clades (Group A) was only partially resolved, with short internal branches pointing to an ancient radiation scenario. Four major clades, Tordyliinae I, Pimpinelleae I, Apieae and Coriandreae, were hypothesized to have originated from chloroplast capture events induced by early hybridization according to the incongruence between chloroplast-based and nrDNA-based phylogenetic trees. Furthermore, the variable and nested distribution of junction positions of LSC (Large single copy region) and IR_B (inverted repeat region B) in Group A may reflect incomplete lineage sorting within this group, which possibly contributed to the unclear phylogenetic relationships among these clades inferred from plastome data. Molecular clock analysis revealed the chloroplast capture events mainly occurred during the middle to late Miocene, providing a geological and climate context for the evolution of Apioideae.

Global warming is causing a more pronounced dip in marine species richness around the equator

The latitudinal gradient in species richness, with more species in the tropics and richness declining with latitude, is widely known and has been assumed to be stable over recent centuries. We analyzed data on 48,661 marine animal species since 1955, accounting for sampling variation, to assess whether the global latitudinal gradient in species richness is being impacted by climate change. We confirm recent studies that show a slight dip in species richness at the equator. Moreover, richness across latitudinal bands was sensitive to temperature, reaching a plateau or declining above a mean annual sea surface temperature of 20 °C for most taxa. In response, since the 1970s, species richness has declined at the equator relative to an increase at midlatitudes and has shifted north in the northern hemisphere, particularly among pelagic species. This pattern is consistent with the hypothesis that climate change is impacting the latitudinal gradient in marine biodiversity at a global scale. The intensification of the dip in species richness at the equator, especially for pelagic species, suggests that it is already too warm there for some species to survive.

Phylogenetic relatedness of woody angiosperm assemblages and its environmental determinants along a subtropical elevational gradient in China

The species composition of plant communities is determined by a number of factors, including current environmental conditions as well as biogeographical and evolutionary history. Despite evidence that plant diversity decreases and species relatedness increases along latitudinal and environmental gradients (e.g., low temperatures), it remains unclear whether these same patterns occur along elevational gradients, especially in the subtropical mountainous areas harboring rich biodiversity. In this study, we explored the pattern of phylogenetic relatedness of woody angiosperm assemblages and examined the effects of temperature variables on the phylogenetic relatedness among angiosperm woody plants using generalized linear model in subtropical forest communities along a broad elevational gradient in the Dulong Valley of Yunnan Province, China. Our results showed that woody angiosperm species in local forest plots tend to be more phylogenetically related at higher elevations and in areas with lower temperatures. Additionally, winter average temperature, rather than mean annual temperature, is a major predictor of the pattern of increasing phylogenetic relatedness with increasing elevation. This finding is consistent with the prediction of 'Tropical Niche Conservatism' hypothesis, which highlights the role of niche constraints in driving phylogenetic community assembly along an elevational gradient.

Biogeography of Iberian Ants (Hymenoptera: Formicidae)

Ants are highly diverse in the Iberian Peninsula (IP), both in species richness (299 cited species) and in number of endemic species (72). The Iberian ant fauna is one of the richest in the broader Mediterranean region, it is similar to the Balkan Peninsula but lower than Greece or Israel, when species richness is controlled by the surface area. In this first general study on the biogeography of Iberian ants, we propose seven chorological categories for grouping thems. Moreover, we also propose eight biogeographic refugium areas, based on the criteria of “refugia-within-refugium” in the IP. We analysed species richness, occurrence and endemism in all these refugium areas, which we found to be significantly different as far as ant similarity was concerned. Finally, we collected published evidence of biological traits, molecular phylogenies, fossil deposits and geological processes to be able to infer the most probable centre of origin and dispersal routes followed for the most noteworthy ants in the IP. As a result, we have divided the Iberian myrmecofauna into four biogeographical groups: relict, Asian-IP disjunct, Baetic-Rifan and Alpine. To sum up, our results support biogeography as being a significant factor for determining the current structure of ant communities, especially in the very complex and heterogenous IP. Moreover, the taxonomic diversity and distribution patterns we describe in this study highlight the utility of Iberian ants for understanding the complex evolutionary history and biogeography of the Iberian Peninsula.

Genome-wide macroevolutionary signatures of key innovations in butterflies colonizing new host plants

The mega-diversity of herbivorous insects is attributed to their co-evolutionary associations with plants. Despite abundant studies on insect-plant interactions, we do not know whether host-plant shifts have impacted both genomic adaptation and species diversification over geological times. We show that the antagonistic insect-plant interaction between swallowtail butterflies and the highly toxic birthworts began 55 million years ago in Beringia, followed by several major ancient host-plant shifts. This evolutionary framework provides a valuable opportunity for repeated tests of genomic signatures of macroevolutionary changes and estimation of diversification rates across their phylogeny. We find that host-plant shifts in butterflies are associated with both genome-wide adaptive molecular evolution (more genes under positive selection) and repeated bursts of speciation rates, contributing to an increase in global diversification through time. Our study links ecological changes, genome-wide adaptations and macroevolutionary consequences, lending support to the importance of ecological interactions as evolutionary drivers over long time periods.

The Phylogeny and Evolution of the Flashiest of the Armored Harvestmen (Arachnida: Opiliones)

Gonyleptoidea, largely restricted to the Neotropics, constitutes the most diverse superfamily of Opiliones and includes the largest and flashiest representatives of this arachnid order. However, the relationships among its main lineages (families and subfamilies) and the timing of their origin are not sufficiently understood to explain how this tropical clade has been able to colonize the temperate zone. Here, we used transcriptomics and divergence time dating to investigate the phylogeny of Gonyleptoidea. Our results support the monophyly of Gonyleptoidea and all of its families with more than one species represented. Resolution within Gonyleptidae s.s. is achieved for many clades, but some subfamilies are not monophyletic (Gonyleptinae, Mitobatinae, and Pachylinae), requiring taxonomic revision. Our data show evidence for one colonization of today's temperate zone early in the history of Gonyleptidae, during the Paleogene, at a time when the Neotropical area extended poleward into regions now considered temperate. This provides a possible mechanism for the colonization of the extratropics by a tropical group following the Paleocene-Eocene Thermal Maximum, explaining how latitudinal diversity gradients can be established. Taxonomic acts: Ampycidae Kury 2003 is newly ranked as family; Neosadocus Mello-Leitão is transferred to Progonyleptoidellinae (new subfamilial assignment). [Arachnids; biogeography; phylogenomics; transcriptomics.].

Sampling effort and the drivers of plant species richness in the Brazilian coastal regions

The causes of the gradients in species richness remain contentious because of multiple competing hypotheses, significant knowledge gaps, and regional effects of environmental and historical factors on species pools. Coastal zones are subject to particular sets of environmental constraints, thus identifying the drivers of species richness therein should shed light on the regional gradients of species diversity. Here, we investigate the geographic patterns and drivers of plant diversity across coastal regions while allowing for pervasive sampling deficiencies. Based on 142708 records of flowering plant occurrences, we mapped species richness and estimated the level of knowledge across the coastal zone of Brazil. Based on inventory completeness, we used linear regression models to test the predictive power of environmental variables that represent different environmental hypotheses. Few cells (25%) were well-surveyed, reflecting little knowledge about the distribution and diversity of flowering plants on the highly-populated Brazilian coast. Still, we found support for the habitat heterogeneity hypothesis as the best explanation of the variation in species richness of flowering plants in this region. Soil properties and water constraints are also important factors. Although our work emphasises the paucity of information on plant diversity in tropical and human-dominated areas, we show that knowledge limitations should not curb our capability of addressing hypotheses about species diversity.

Ancient tropical extinctions at high latitudes contributed to the latitudinal diversity gradient

Global biodiversity currently peaks at the equator and decreases toward the poles. Growing fossil evidence suggest this hump-shaped latitudinal diversity gradient (LDG) has not been persistent through time, with similar diversity across latitudes flattening out the LDG during past greenhouse periods. However, when and how diversity declined at high latitudes to generate the modern LDG remains an open question. Although diversity-loss scenarios have been proposed, they remain mostly undemonstrated. We outline the "asymmetric gradient of extinction and dispersal" framework that contextualizes previous ideas behind the LDG under a time-variable scenario. Using phylogenies and fossils of Testudines, Crocodilia, and Lepidosauria, we find that the hump-shaped LDG could be explained by (1) disproportionate extinctions of high-latitude tropical-adapted clades when climate transitioned from greenhouse to icehouse, and (2) equator-ward biotic dispersals tracking their climatic preferences when tropical biomes became restricted to the equator. Conversely, equivalent diversification rates across latitudes can account for the formation of an ancient flat LDG. The inclusion of fossils in macroevolutionary studies allows revealing time-dependent extinction rates hardly detectable from phylogenies only. This study underscores that the prevailing evolutionary processes generating the LDG during greenhouses differed from those operating during icehouses.

Phylogenomic test of mitochondrial clues to archaic ancestors in a group of hybridizing swallowtail butterflies

Genomics has revolutionized our understanding of hybridization and introgression, but most of the early evidence for these processes came from studies of mitochondrial introgression. To expand these evolutionary insights from mitochondrial patterns, we evaluate phylogenetic discordance across the nuclear genomes of a hybridizing system, the Papilio machaon group of swallowtail butterflies. This species group contains three hybrid lineages (P. brevicauda, P. joanae, and P. m. kahli) that are geographically disjunct across North America and have complete fixation of a mitochondrial lineage that is otherwise primarily found in P. m. hudsonianus, a boreal subspecies of the Holarctic P. machaon. Genome-wide nuclear markers place the three hybrid lineages as a monophyletic group that is sister to P. polyxenes/P. zelicaon rather than P. machaon, although ancient hybridization between a subspecies of P. machaon and the ancestor of these three lineages is also shown by their greater nuclear affinity to P. m. hudsonianus than to other subspecies of P. machaon. Individuals from contemporary hybrid swarms in Alberta, where mitochondrial DNA fixation has not occurred, were more intermediate between their respective parent species, demonstrating diversity in mito-nuclear discordance following hybrid interactions. Our new phylogenetic findings for the P. machaon species group also include: subspecific paraphyly within P. machaon itself across its Holarctic distribution; paraphyly of P. zelicaon relative to P. polyxenes; and more divergent placement of a Mediterranean species, P. hospiton. These results provide the first comprehensive genomic evaluation of relationships within this species group and provide insight into the evolutionary dynamics of hybridization and mitochondrial introgression.

Are the Yellow and Red Marked Club-Tail Losaria coon the Same Species?

Losaria coon (Fabricius, 1793) is currently comprised of ten subspecies, which were originally described under two names, Papilio coon and P. doubledayi before 1909, when they were combined as one species. The main difference between them is the colour of abdomen and hindwing subterminal spots-yellow in coon and red in doubledayi. Wing morphology, male and female genitalia, and molecular evidence (DNA barcodes) were analysed for multiple subspecies of L. coon and three other Losaria species-rhodifer, neptunus, and palu. Our molecular data support the separation of L. coon and L. doubledayi stat. rev. as two distinct species, with L. rhodifer positioned between them in phylogenetic analyses. Wing morphology and genitalic structures also confirm the molecular conclusions. Our findings divide L. coon into two species occupying different geographic ranges: with L. coon restricted to southern Sumatra, Java, and Bawean Island, while L. doubledayi occurs widely in regions from North India to northern Sumatra, including Hainan and Nicobar Islands. Hence, future conservation efforts must reassess the status and threat factors of the two species to form updated strategies.

The latitudinal diversity gradient of tetrapods across the Permo-Triassic mass extinction and recovery interval

The decline in species richness from the equator to the poles is referred to as the latitudinal diversity gradient (LDG). Higher equatorial diversity has been recognized for over 200 years, but the consistency of this pattern in deep time remains uncertain. Examination of spatial biodiversity patterns in the past across different global climate regimes and continental configurations can reveal how LDGs have varied over Earth history and potentially differentiate between suggested causal mechanisms. The Late Permian-Middle Triassic represents an ideal time interval for study, because it is characterized by large-scale volcanic episodes, extreme greenhouse temperatures and the most severe mass extinction event in Earth history. We examined terrestrial and marine tetrapod spatial biodiversity patterns using a database of global tetrapod occurrences. Terrestrial tetrapods exhibit a bimodal richness distribution throughout the Late Permian-Middle Triassic, with peaks in the northern low latitudes and southern mid-latitudes around 20-40° N and 60° S, respectively. Marine reptile fossils are known almost exclusively from the Northern Hemisphere in the Early and Middle Triassic, with highest diversity around 20° N. Reconstructed terrestrial LDGs contrast strongly with the generally unimodal gradients of today, potentially reflecting high global temperatures and prevailing Pangaean super-monsoonal climate system during the Permo-Triassic.

Phylogenetic niche conservatism and plant diversification in South American subtropical grasslands along multiple climatic dimensions

Phylogenetic niche conservatism can be investigated at multiple scales on an explicit geographical context. Haplotype-based comparative analyses of lineages occupying the same region, and thus subjected to similar environmental factors, allow decoupling shared evolutionary and ecological patterns, as well as multiple dimensions of adaptive diversification. Here we aimed to assess the role of environmental drivers on diversification of subtropical grassland, based on haplotypic diversity of two plant genera. We sampled two closely related and co-distributed grassland plant genera, Petunia and Calibrachoa, across their entire distribution area. Eigenvectors extracted from pairwise distances based on chloroplast DNA haplotypes were used to fit Phylogenetic Signal-Representation (PSR) curves to estimate evolutionary patterns in 19 bioclimatic variables and altitude. The PSR curves showed that altitude, precipitation, and temperature variables changed at different rates with haplotype differentiation. Altitude and temperature traits evolved under conditions closer to a neutral dynamics, whereas precipitation traits differentiated following more complex models. Our results indicated that the diversification in the two genera was more limited by precipitation conditions. Based on these novel findings, we suggest that future studies should test the possible impact of precipitation variables on the process of ecological differentiation in these genera.

Regional assemblages shaped by historical and contemporary factors: Evidence from a species-rich insect group

Understanding diversity patterns requires accounting for the roles of both historical and contemporary factors in the assembly of communities. Here, we compared diversity patterns of two moth assemblages sampled from Taihang and Yanshan mountains in Northern China and performed ancestral range reconstructions using the Multi-State Speciation and Extinction model, to track the origins of these patterns. Further, we estimated diversification rates of the two moth assemblages and explored the effects of contemporary ecological factors. From 7,788 specimens we identified 835 species belonging to 23 families, using both DNA barcode analysis and morphology. Moths in Yanshan mountains showed higher species diversity than in Taihang mountains. Ancestral range analysis indicated Yanshan as the origin, with significant historical dispersals from Yanshan to Taihang. Asymmetrical diversification, population expansion, along with frequent and considerable gene flow were detected between communities. Moreover, dispersal limitation or the joint effect of environment filtering and dispersal limitation were inferred as main driving forces shaping current diversity patterns. In summary, we demonstrate that a multiscale (community, population and species level) analysis incorporating both historical and contemporary factors can be useful in delineating factors contributing to community assembly and patterning in diversity.

Angiosperm speciation cools down in the tropics

Recent evidence has questioned whether the Latitudinal Diversity Gradient (LDG), whereby species richness increases towards the Equator, results in higher rates of speciation in the tropics. Allowing for time heterogeneity in speciation rate estimates for over 60,000 angiosperm species, we found that the LDG does not arise from variation in speciation rates because lineages do not speciate faster in the tropics. These results were consistently retrieved using two other methods to test the association between occupancy of tropical habitats and speciation rates. Our speciation rate estimates were robust to the effects of both undescribed species and missing taxa. Overall, our results show that speciation rates follow an opposite pattern to global variation in species richness. Greater ecological opportunity in the temperate zones, stemming from less saturated communities, higher species turnover or greater environmental change, may ultimately explain these results.

Priors and Posteriors in Bayesian Timing of Divergence Analyses: The Age of Butterflies Revisited

The need for robust estimates of times of divergence is essential for downstream analyses, yet assessing this robustness is still rare. We generated a time-calibrated genus-level phylogeny of butterflies (Papilionoidea), including 994 taxa, up to 10 gene fragments and an unprecedented set of 12 fossils and 10 host-plant node calibration points. We compared marginal priors and posterior distributions to assess the relative importance of the former on the latter. This approach revealed a strong influence of the set of priors on the root age but for most calibrated nodes posterior distributions shifted from the marginal prior, indicating significant information in the molecular data set. Using a very conservative approach we estimated an origin of butterflies at 107.6 Ma, approximately equivalent to the latest Early Cretaceous, with a credibility interval ranging from 89.5 Ma (mid Late Cretaceous) to 129.5 Ma (mid Early Cretaceous). In addition, we tested the effects of changing fossil calibration priors, tree prior, different sets of calibrations and different sampling fractions but our estimate remained robust to these alternative assumptions. With 994 genera, this tree provides a comprehensive source of secondary calibrations for studies on butterflies.

A shared genetic basis of mimicry across swallowtail butterflies points to ancestral co-option of doublesex

Uncovering whether convergent adaptations share a genetic basis is consequential for understanding the evolution of phenotypic diversity. This information can help us understand the extent to which shared ancestry or independent evolution shape adaptive phenotypes. In this study, we first ask whether the same genes underlie polymorphic mimicry in Papilio swallowtail butterflies. By comparing signatures of genetic variation between polymorphic and monomorphic species, we then investigate how ancestral variation, hybridization, and independent evolution contributed to wing pattern diversity in this group. We report that a single gene, doublesex (dsx), controls mimicry across multiple taxa, but with species-specific patterns of genetic differentiation and linkage disequilibrium. In contrast to widespread examples of phenotypic evolution driven by introgression, our analyses reveal distinct mimicry alleles. We conclude that mimicry evolution in this group was likely facilitated by ancestral polymorphism resulting from early co-option of dsx as a mimicry locus, and that evolutionary turnover of dsx alleles may underlie the wing pattern diversity of extant polymorphic and monomorphic lineages.

Molecular phylogenetics disentangles the limits of a group of related Neotropical doryctine parasitoid wasp genera (Hymenoptera : Braconidae)

The Doryctinae is one of the most specious subfamilies of the mainly parasitoid wasp family Braconidae. Members of this subfamily are characterised by having a very heterogeneous external morphology, which has made establishing the limits of several of its genera extremely difficult. The Neotropical Callihormius Ashmead, Platydoryctes Barbalho & Penteado-Dias, and Aphelopsia Marsh are prime examples of this, since they show extensive interspecific morphological variation but lack consistent diagnostic features among them. We carried out a comprehensive phylogenetic study among several representative species belonging to the above genera as well as other morphologically similar doryctine taxa based on three nuclear and two mitochondrial (mt)DNA sequence markers. Based on the relationships recovered, we assessed the generic limits among the examined taxa and investigated the times of origin and diversification in this group. Our best estimate of phylogeny did not recover Callihormius as monophyletic, and thus we propose a taxonomic arrangement where we elevate the status of Ca. (Dmitriohormius) Belokobylskij, Zaldívar-Riverón & Coronado-Blanco to genus level, stat. nov., and also describe four new genera: Macrometasoma, gen. nov., Caputlenis, gen. nov., Platyhormius, gen. nov. and Caputrugosus, gen. nov. Platydoryctes is delimited morphologically, and the body flatness is proposed to have independently evolved more than once within the group of genera examined. Aphelopsia was not recovered as monophyletic, though the relationships involved were weakly supported and therefore we maintain its species composition. Panama Marsh, stat. rev., which was considered a junior synonym of Leluthia Cameron, is resurrected. Basal relationships in the phylogeny had very short branch lengths, which could be due to ancient rapid radiation events that occurred during the early evolution of the group. The estimated times of divergence showed that the clade containing Callihormius and its related genera originated between the late Oligocene and middle Miocene, whereas its early diversification events probably occurred between the middle to late Miocene.

Phylogenetic dispersion and diversity in regional assemblages of seed plants in China

Species assemble into communities through ecological and evolutionary processes. Phylogenetic niche conservatism-the tendency of species to retain ancestral ecological distributions-is thought to influence which species from a regional species pool can persist in a particular environment. We analyzed data for seed plants in China to test hypotheses about the distribution of species within regional floras. Of 16 environmental variables, actual evapotranspiration, minimum temperature of the coldest month, and annual precipitation most strongly influenced regional species richness, phylogenetic dispersion, and phylogenetic diversity for both gymnosperms (cone-bearing plants) and angiosperms (flowering plants). For most evolutionary clades at, and above, the family level, the relationships between metrics of phylogenetic dispersion (i.e., average phylogenetic distance among species), or phylogenetic diversity, and the 3 environmental variables were consistent with the tropical niche conservatism hypothesis, which predicts closer phylogenetic relatedness and reduced phylogenetic diversity with increasing environmental stress. The slopes of the relationships between phylogenetic relatedness and the 3 environmental drivers identified in this analysis were steeper for primarily tropical clades, implying greater niche conservatism, than for primarily temperate clades. These observations suggest that the distributions of seed plants across large-scale environmental gradients in China are constrained by conserved adaptations to the physical environment, i.e., phylogenetic niche conservatism.

Host and geography together drive early adaptive radiation of Hawaiian planthoppers

The interactions between insects and their plant host have been implicated in driving diversification of both players. Early arguments highlighted the role of ecological opportunity, with the idea that insects "escape and radiate" on new hosts, with subsequent hypotheses focusing on the interplay between host shifting and host tracking, coupled with isolation and fusion, in generating diversity. Because it is rarely possible to capture the initial stages of diversification, it is particularly difficult to ascertain the relative roles of geographic isolation versus host shifts in initiating the process. The current study examines genetic diversity between populations and hosts within a single species of endemic Hawaiian planthopper, Nesosydne umbratica (Hemiptera, Delphacidae). Given that the species was known as a host generalist occupying unrelated hosts, Clermontia (Campanulaceae) and Pipturus (Urticaceae), we set out to determine the relative importance of geography and host in structuring populations in the early stages of differentiation on the youngest islands of the Hawaiian chain. Results from extensive exon capture data showed that N. umbratica is highly structured, both by geography, with discrete populations on each volcano, and by host plant, with parallel radiations on Clermontia and Pipturus leading to extensive co-occurrence. The marked genetic structure suggests that populations can readily become established on novel hosts provided opportunity; subsequent adaptation allows monopolization of the new host. The results support the role of geographic isolation in structuring populations and with host shifts occurring as discrete events that facilitate subsequent parallel geographic range expansion.

Whole Genome Shotgun Phylogenomics Resolves the Pattern and Timing of Swallowtail Butterfly Evolution

Evolutionary relationships have remained unresolved in many well-studied groups, even though advances in next-generation sequencing and analysis, using approaches such as transcriptomics, anchored hybrid enrichment, or ultraconserved elements, have brought systematics to the brink of whole genome phylogenomics. Recently, it has become possible to sequence the entire genomes of numerous nonbiological models in parallel at reasonable cost, particularly with shotgun sequencing. Here, we identify orthologous coding sequences from whole-genome shotgun sequences, which we then use to investigate the relevance and power of phylogenomic relationship inference and time-calibrated tree estimation. We study an iconic group of butterflies—swallowtails of the family Papilionidae—that has remained phylogenetically unresolved, with continued debate about the timing of their diversification. Low-coverage whole genomes were obtained using Illumina shotgun sequencing for all genera. Genome assembly coupled to BLAST-based orthology searches allowed extraction of 6621 orthologous protein-coding genes for 45 Papilionidae species and 16 outgroup species (with 32% missing data after cleaning phases). Supermatrix phylogenomic analyses were performed with both maximum-likelihood (IQ-TREE) and Bayesian mixture models (PhyloBayes) for amino acid sequences, which produced a fully resolved phylogeny providing new insights into controversial relationships. Species tree reconstruction from gene trees was performed with ASTRAL and SuperTriplets and recovered the same phylogeny. We estimated gene site concordant factors to complement traditional node-support measures, which strengthens the robustness of inferred phylogenies. Bayesian estimates of divergence times based on a reduced data set (760 orthologs and 12% missing data) indicate a mid-Cretaceous origin of Papilionoidea around 99.2 Ma (95% credibility interval: 68.6–142.7 Ma) and Papilionidae around 71.4 Ma (49.8–103.6 Ma), with subsequent diversification of modern lineages well after the Cretaceous-Paleogene event. These results show that shotgun sequencing of whole genomes, even when highly fragmented, represents a powerful approach to phylogenomics and molecular dating in a group that has previously been refractory to resolution.

Morphological Innovations and Vast Extensions of Mountain Habitats Triggered Rapid Diversification Within the Species-Rich Irano-Turanian Genus Acantholimon (Plumbaginaceae)

The Irano-Turanian floristic region spans a topographically complex and climatically continental territory, which has served as a source of xerophytic taxa for neighboring regions and is represented by a high percent of endemics. Yet, a comprehensive picture of the abiotic and biotic factors that have driven diversification within this biota remains to be established due to the scarcity of phylogenetic studies. Acantholimon is an important component of the subalpine steppe flora of the Irano-Turanian region, containing c. 200 cushion-forming sub-shrubby pungent-leaved species. Our recent molecular phylogenetic study has led to enlarging the circumscription of this genus to include eight mono- or oligospecific genera lacking the characteristic life-form and leaves. Using the same molecular phylogeny, here we investigate the tempo and mode of diversification as well as the biogeographic patterns in this genus, to test the hypothesis that a combination of key morphological innovations and abiotic factors is behind Acantholimon high species diversity. Molecular dating analysis indicates that Acantholimon s.l. started to diversify between the Late Miocene and the Pliocene and the biogeographic analysis points to an Eastern Iran-Afghanistan origin. Macroevolutionary models support the hypothesis that the high diversity of the genus is explained by accelerated diversification rates in two clades associated with the appearance of morphological key innovations such as a cushion life-form and pungent leaves; this would have favored the colonization of water-stressed, substrate-poor mountainous habitats along the newly uplifted IT mountains during the Mio-Pliocene. Given the apparent similarity of mountain habitats for most species of Acantholimon, we hypothesize that its current high species diversity responds to a scenario of non-adaptive radiation fueled by allopatric speciation rather than evolutionary radiation driven by ecological opportunity. Similar scenarios might underlie the high diversity of other speciose genera in the topographically complex Irano-Turanian landscape, though this remains to be tested with fine-grained distribution and climatic data.

Detecting Environment-Dependent Diversification From Phylogenies: A Simulation Study and Some Empirical Illustrations

Understanding the relative influence of various abiotic and biotic variables on diversification dynamics is a major goal of macroevolutionary studies. Recently, phylogenetic approaches have been developed that make it possible to estimate the role of various environmental variables on diversification using time-calibrated species trees, paleoenvironmental data, and maximum-likelihood techniques. These approaches have been effectively employed to estimate how speciation and extinction rates vary with key abiotic variables, such as temperature and sea level, and we can anticipate that they will be increasingly used in the future. Here we compile a series of biotic and abiotic paleodatasets that can be used as explanatory variables in these models and use simulations to assess the statistical properties of the approach when applied to these paleodatasets. We demonstrate that environment-dependent models perform well in recovering environment-dependent speciation and extinction parameters, as well as in correctly identifying the simulated environmental model when speciation is environment-dependent. We explore how the strength of the environment-dependency, tree size, missing taxa, and characteristics of the paleoenvironmental curves influence the performance of the models. Finally, using these models, we infer environment-dependent diversification in two empirical phylogenies: temperature-dependence in Cetacea and $\delta^{13}C$-dependence in Ruminantia. We illustrate how to evaluate the relative importance of abiotic and biotic variables in these two clades and interpret these results in light of macroevolutionary hypotheses. Given the important role paleoenvironments are presumed to have played in species evolution, our statistical assessment of how environment-dependent models behave is crucial for their utility in macroevolutionary analysis.

Evolution of ontogenic change in color defenses of swallowtail butterflies

Natural selection by visually hunting predators has led to the evolution of color defense strategies such as masquerade, crypsis, and aposematism that reduce the risk of predation in prey species. These color defenses are not mutually exclusive, and switches between strategies with ontogenic development are widespread across taxa. However, the evolutionary dynamics of ontogenic color change are poorly understood. Using comparative phylogenetics, we studied the evolution of color defenses in the complex life cycles of swallowtail butterflies (family Papilionidae). We also tested the relative importance of life history traits, chemical and visual backgrounds, and ancestry on the evolution of protective coloration. We found that vulnerable early- and late-instar caterpillars of species that feed on sparsely vegetated, toxic plants were aposematic, whereas species that feed on densely vegetated, nontoxic plants had masquerading and cryptic caterpillars. Masquerading caterpillars resembled bird droppings at early instars and transitioned to crypsis with an increase in body size at late instars. The immobile pupae-safe from motion-detecting, visually hunting predators-retained the ancestral cryptic coloration in all lineages, irrespective of the toxic nature of the host plant. Thus, color defense strategy (masquerade, crypsis, or aposematism) at a particular lifestage in the life cycle of swallowtail butterflies was determined by the interaction between life history traits such as body size and motion levels, phytochemical and visual backgrounds, and ancestry. We show that ontogenic color change in swallowtail butterflies is an adaptive response to age-dependent vulnerability to predation.

Climate and host-plant associations shaped the evolution of ceutorhynch weevils throughout the Cenozoic

Using molecular phylogenetic data and methods we inferred divergence times and diversification patterns for the weevil subfamily Ceutorhynchinae in the context of host‐plant associations and global climate over evolutionary time. We detected four major diversification shifts that correlate with both host shifts and major climate events. Ceutorhynchinae experienced an increase in diversification rate at ∼53 Ma, during the Early Eocene Climate Optimum, coincident with a host shift to Lamiaceae. A second major diversification phase occurred at the end of the Eocene (∼34 Ma). This contrasts with the overall deterioration in climate equability at the Eocene‐Oligocene boundary, but tracks the diversification of important host plant clades in temperate (higher) latitudes, leading to increased diversification rates in the weevil clades infesting temperate hosts. A third major phase of diversification is correlated with the rising temperatures of the Late Oligocene Warming Event (∼26.5 Ma); diversification rates then declined shortly after the Middle Miocene Climate Transition (∼14.9 Ma). Our results indicate that biotic and abiotic factors together explain the evolution of Ceutorhynchinae better than each of these drivers viewed in isolation.

Macroecology and macroevolution of the latitudinal diversity gradient in ants

The latitudinal diversity gradient—the tendency for more species to occur toward the equator—is the dominant pattern of life on Earth, yet the mechanisms responsible for it remain largely unexplained. Recently, the analysis of global data has led to advances in understanding, but these advances have been mostly limited to vertebrates and trees and have not provided consensus answers. Here we synthesize large-scale geographic, phylogenetic, and fossil data for an exemplar invertebrate group—ants—and investigate whether the latitudinal diversity gradient arose due to higher rates of net diversification in the tropics, or due to a longer time period to accumulate diversity due to Earth’s climatic history. We find that latitudinal affinity is highly conserved, temperate clades are young and clustered within tropical clades, and diversification rate shows no systematic variation with latitude. These results indicate that diversification time—and not rate—is the main driver of the diversity gradient in ants.

Multiple hypotheses have been proposed for the declining biodiversity gradient between the tropics and poles. Here, the authors compile and analyze geographic data for all ant species and large-scale phylogenies, suggesting that diversification time drives the latitudinal diversity gradient in ants.