Accumulation over evolutionary time as a major cause of biodiversity hotspots in conifers

Chromosomal evolution, environmental heterogeneity, and migration drive spatial patterns of species richness in Calochortus (Liliaceae)

We used nuclear genomic data and statistical models to evaluate the ecological and evolutionary processes shaping spatial variation in species richness in Calochortus (Liliaceae, 74 spp.). Calochortus occupies diverse habitats in the western United States and Mexico and has a center of diversity in the California Floristic Province, marked by multiple orogenies, winter rainfall, and highly divergent climates and substrates (including serpentine). We used sequences of 294 low-copy nuclear loci to produce a time-calibrated phylogeny, estimate historical biogeography, and test hypotheses regarding drivers of present-day spatial patterns in species number. Speciation and species coexistence require reproductive isolation and ecological divergence, so we examined the roles of chromosome number, environmental heterogeneity, and migration in shaping local species richness. Six major clades-inhabiting different geographic/climatic areas, and often marked by different base chromosome numbers (n = 6 to 10)-began diverging from each other ~10.3 Mya. As predicted, local species number increased significantly with local heterogeneity in chromosome number, elevation, soil characteristics, and serpentine presence. Species richness is greatest in the Transverse/Peninsular Ranges where clades with different chromosome numbers overlap, topographic complexity provides diverse conditions over short distances, and several physiographic provinces meet allowing immigration by several clades. Recently diverged sister-species pairs generally have peri-patric distributions, and maximum geographic overlap between species increases over the first million years since divergence, suggesting that chromosomal evolution, genetic divergence leading to gametic isolation or hybrid inviability/sterility, and/or ecological divergence over small spatial scales may permit species co-occurrence.

A revised, annotated checklist of Mexican non-biting midges (Diptera, Chironomidae)

An updated checklist of Mexican non-biting midges (Chironomidae) is presented. A total of 110 species of Chironomidae are known for Mexico: 52 species in 25 genera belong to the subfamily Chironominae, 30 species in 13 genera to Orthocladiinae, 21 species in nine genera to Tanypodinae, five species in two genera to Telmatogetoninae, and two species in one genus to Diamesinae. In addition, 41 genera without identified species are listed. The highest number of species (29) is recorded from the state of Campeche, while 19 species have been found in Veracruz and 15 in Nuevo León. Few or no records exist for states in Central and Northern Mexico, or those on the Pacific coast. The type localities for 34 species are in Mexico; of these, 27 species (25% of the total number of species recorded in the country) are endemic. Twenty-nine species recorded in Mexico have a Neotropical distribution, 15 a Nearctic distribution, and 39 species are distributed in both the Neotropical and Nearctic regions or more widely. It has been suggested that as many as 1000 species might occur in Mexico; so only a little more than 10% of the expected diversity has so far been recorded.

Phylogenetic risk assessment is robust for forecasting the impact of European insects on North American conifers

Some introduced species cause severe damage, although the majority have little impact. Robust predictions of which species are most likely to cause substantial impacts could focus efforts to mitigate those impacts or prevent certain invasions entirely. Introduced herbivorous insects can reduce crop yield, fundamentally alter natural and managed forest ecosystems, and are unique among invasive species in that they require certain host plants to succeed. Recent studies have demonstrated that understanding the evolutionary history of introduced herbivores and their host plants can provide robust predictions of impact. Specifically, divergence times between hosts in the native and introduced ranges of a nonnative insect can be used to predict the potential impact of the insect should it establish in a novel ecosystem. However, divergence time estimates vary among published phylogenetic datasets, making it crucial to understand if and how the choice of phylogeny affects prediction of impact. Here, we tested the robustness of impact prediction to variation in host phylogeny by using insects that feed on conifers and predicting the likelihood of high impact using four different published phylogenies. Our analyses ranked 62 insects that are not established in North America and 47 North American conifer species according to overall risk and vulnerability, respectively. We found that results were robust to the choice of phylogeny. Although published vascular plant phylogenies continue to be refined, our analysis indicates that those differences are not substantial enough to alter the predictions of invader impact. Our results can assist in focusing biosecurity programs for conifer pests and can be more generally applied to nonnative insects and their potential hosts by prioritizing surveillance for those insects most likely to be damaging invaders.

At least it is a dry cold: the global distribution of freeze-thaw and drought stress and the traits that may impart poly-tolerance in conifers

Conifers inhabit some of the most challenging landscapes where multiple abiotic stressors (e.g., aridity, freezing temperatures) often co-occur. Physiological tolerance to multiple stressors ('poly-tolerance') is thought to be rare because exposure to one stress generally limits responses to another through functional trade-offs. However, the capacity to exhibit poly-tolerance may be greater when combined abiotic stressors have similar physiological impacts, such as the disruption of hydraulic function imposed by drought or freezing. Here, we reviewed empirical data in light of theoretical expectations for conifer adaptations to drought and freeze-thaw cycles with particular attention to hydraulic traits of the stem and leaf. Additionally, we examined the commonality and spatial distribution of poly-stress along indices of these combined stressors. We found that locations with the highest values of our poly-stress index (PSi) are characterized by moderate drought and moderate freeze-thaw, and most of the global conifer distribution occupies areas of moderate poly-stress. Among traits examined, we found diverse responses to the stressors. Turgor loss point did not correlate with freeze-thaw or drought stress individually, but did with the PSi, albeit inverse to what was hypothesized. Leaf mass per area was more strongly linked with drought stress than the poly-stress and not at all with freeze-thaw stress. In stems, the water potential causing 50% loss of hydraulic conductivity became more negative with increasing drought stress and poly-stress but did not correlate with freeze-thaw stress. For these traits, we identified a striking lack of coverage for substantial portions of species ranges, particularly at the upper boundaries of their respective PSis, demonstrating a critical gap in our understanding of trait prevalence and plasticity along these stress gradients. Future research should investigate traits that confer tolerance to both freeze-thaw and drought stress in a wide range of species across broad geographic scales.

Neontological and paleontological congruence in the evolution of Podocarpaceae (coniferales) reproductive morphology

Introduction

Podocarpaceae are a diverse, primarily tropical conifer family that commonly produce large leaves and highly reduced, fleshy seed cones bearing large seeds. These features may result from relatively recent adaptation to closed-canopy angiosperm forests and bird-mediated seed dispersal, although determining precisely when shifts in leaf and seed cone morphology occurred is difficult due to a sparse fossil record and relatively few surviving deep lineages.

Methods

We compare the fossil record of Podocarpaceae with results from ancestral state reconstruction methods and correlated character models using neontological data and a previously published molecular time-tree.

Results

Ancestral state reconstructions suggest that small leaves, small seeds, and multi-seeded cones are ancestral in crown Podocarpaceae, with reduced cones bearing few seeds appearing in the Early Cretaceous and the correlated evolution of large leaves and large seeds occurring from the Late Cretaceous onwards. The exact timing of these shifts based on neontological data alone are poorly constrained, however, and estimates of leaf and seed size are imprecise.

Discussion

The fossil record is largely congruent with results based on the molecular time-tree, but provide important constraints on the range of leaf and seed sizes that were present in Cretaceous Podocarpaceae and the time by which changes in cone morphology and seed size likely occurred. We suggest in particular that reduced seed cones appeared in the Early Cretaceous and are linked to the contemporaneous diversification of small bodied avialans (birds), with shifts to larger seed sizes occurring after the Cretaceous in association with the spread of closed-canopy angiosperm forests.

Species richness patterns and the determinants of larch forests in China

Larch forests are important for species diversity, as well as soil and water conservation in mountain regions. In this study, we determined large-scale patterns of species richness in larch forests and identified the factors that drive these patterns. We found that larch forest species richness was high in southern China and low in northern China, and that patterns of species richness along an elevational gradient depend on larch forest type. In addition, we found that patterns of species richness in larch forests are best explained by contemporary climatic factors. Specifically, mean annual temperature and annual potential evapotranspiration were the most important factors for species richness of tree and shrub layers, while mean temperature of the coldest quarter and anomaly of annual precipitation from the Last Glacial Maximum to the present were the most important for that of herb layer and the whole community. Community structural factors, especially stand density, are also associated with the species richness of larch forests. Our findings that species richness in China's larch forests is mainly affected by energy availability and cold conditions support the ambient energy hypothesis and the freezing tolerance hypothesis.

Incipient speciation, high genetic diversity, and ecological divergence in the alligator bark juniper suggest complex demographic changes during the Pleistocene

The most recent glacial cycles of the Pleistocene affected the distribution, population sizes, and levels of genetic structure of temperate-forest species in the main Mexican mountain systems. Our objective was to investigate the effects these cycles had on the genetic structure and distribution of a dominant species of the "mexical" vegetation across North and Central America. We studied the genetic diversity of Juniperus deppeana, a conifer distributed from the Southwestern United States to the highlands of Central America. We combined information of one plastid marker and two nuclear markers to infer phylogeographic structure, genetic diversity and demographic changes. We also characterized the climatic niche for each variety to infer the plausible area of suitability during past climatic conditions and to evaluate climatic niche discontinuities along with the species distribution. We found a marked phylogeographic structure separating the populations North and South of the Isthmus of Tehuantepec, with populations to the South of this barrier forming a distinct genetic cluster corresponding to Juniperus deppeana var. gamboana. We also found signals of population expansion in the Northern genetic cluster. Ecological niche modeling results confirmed climatic niche differences and discontinuities among J. deppeana varieties and heterogeneous responses to climatic oscillations. Overall, J. deppeana's genetic diversity has been marked by distribution shifts, population growth and secondary contact the North, and in situ permanence in the South since the last interglacial to the present. High genetic variation suggests a wide and climatically diverse distribution during climatic oscillations. We detected the existence of two main genetic clusters, supporting previous proposals that Juniperus deppeana and Juniperus gamboana may be considered two separate species.

An integrated high-resolution mapping shows congruent biodiversity patterns of Fagales and Pinales

The documentation of biodiversity distribution through species range identification is crucial for macroecology, biogeography, conservation, and restoration. However, for plants, species range maps remain scarce and often inaccurate. We present a novel approach to map species ranges at a global scale, integrating polygon mapping and species distribution modelling (SDM). We develop a polygon mapping algorithm by considering distances and nestedness of occurrences. We further apply an SDM approach considering multiple modelling algorithms, complexity levels, and pseudo-absence selections to map the species at a high spatial resolution and intersect it with the generated polygons. We use this approach to construct range maps for all 1957 species of Fagales and Pinales with data compilated from multiple sources. We construct high-resolution global species richness maps of these important plant clades, and document diversity hotspots for both clades in southern and south-western China, Central America, and Borneo. We validate the approach with two representative genera, Quercus and Pinus, using previously published coarser range maps, and find good agreement. By efficiently producing high-resolution range maps, our mapping approach offers a new tool in the field of macroecology for studying global species distribution patterns and supporting ongoing conservation efforts.

Climate Stability Index maps, a global high resolution cartography of climate stability from Pliocene to 2100

Climate changes are top biodiversity shapers, both during the past and future. Mapping the most climatic stable and unstable zones on Earth could improve our understanding of biodiversity distribution and evolution. Here, we present a set of maps based on a global scale, high resolution (ca. 5 km) new Climate Stability Index (CSI). The CSI considers bioclimatic variables for two different time ranges: (1) from Pliocene (3.3 Ma) to the present (CSI-past map set), using 12 time periods of PaleoClim representing warm and cold cycles; and (2) from present to the year 2100 (CSI-future), using nine general circulation models of climate change of four periods available from WorldClim. We calculated standard deviation of the variables and selected an uncorrelated set for summing, normalizing and obtaining the CSI maps. Our approach is useful for fields such as biogeography, earth sciences, agriculture, or sociology. However, CSI is an index that can be re-calculated according to particular criteria and objectives (e.g. temperature variables); maps are, therefore, customizable to every user.

Phylogenetic diversity and conservation of crop wild relatives in Colombia

Crop wild relatives (CWR) are an important agricultural resource as they contain genetic traits not found in cultivated species due to localized adaptation to unique environmental and climatic conditions. Phylogenetic diversity (PD) measures the evolutionary relationship of species using the tree of life. Our knowledge of CWR PD in neotropical regions is in its infancy. We analysed the distribution of CWR PD across Colombia and assessed its conservation status. The areas with the largest concentration of PD were identified as being in the northern part of the central and western Andean mountain ranges and the Pacific region. These centres of high PD were comprised of predominantly short and closely related branches, mostly of species of wild tomatoes and black peppers. In contrast, the CWR PD in the lowland ecosystems of the Amazon and Orinoquia regions had deeply diverging clades predominantly represented by long and distantly related branches (i.e. tuberous roots, grains and cacao). We categorized 50 (52.6%) of the CWR species as 'high priority', 36 as 'medium priority' and nine as 'low priority' for further ex-situ and in situ conservation actions. New areas of high PD and richness with large ex-situ gap collections were identified mainly in the northern part of the Andes of Colombia. We found that 56% of the grid cells with the highest PD values were unprotected. These baseline data could be used to create a comprehensive national strategy of CWR conservation in Colombia.

Phylogenomic and ecological analyses reveal the spatiotemporal evolution of global pines

How coniferous forests evolved in the Northern Hemisphere remains largely unknown. Unlike most groups of organisms that generally follow a latitudinal diversity gradient, most conifer species in the Northern Hemisphere are distributed in mountainous areas at middle latitudes. It is of great interest to know whether the midlatitude region has been an evolutionary cradle or museum for conifers and how evolutionary and ecological factors have driven their spatiotemporal evolution. Here, we investigated the macroevolution of Pinus, the largest conifer genus and characteristic of northern temperate coniferous forests, based on nearly complete species sampling. Using 1,662 genes from transcriptome sequences, we reconstructed a robust species phylogeny and reestimated divergence times of global pines. We found that ∼90% of extant pine species originated in the Miocene in sharp contrast to the ancient origin of Pinus, indicating a Neogene rediversification. Surprisingly, species at middle latitudes are much older than those at other latitudes. This finding, coupled with net diversification rate analysis, indicates that the midlatitude region has provided an evolutionary museum for global pines. Analyses of 31 environmental variables, together with a comparison of evolutionary rates of niche and phenotypic traits with a net diversification rate, found that topography played a primary role in pine diversification, and the aridity index was decisive for the niche rate shift. Moreover, fire has forced diversification and adaptive evolution of Pinus Our study highlights the importance of integrating phylogenomic and ecological approaches to address evolution of biological groups at the global scale.

Dispersal patterns of an introduced wild bee, Megachile sculpturalis Smith, 1853 (Hymenoptera: Megachilidae) in European alpine countries

Biodiversity monitoring programs are the baseline of species abundancy studies, which in case of introduced species are especially critical. Megachile sculpturalis Smith, 1853 native to Eastern-Asia, constitutes the first ever recorded wild bee species accidently introduced in Europe. Since its first discovery in 2008, M. sculpturalis has been spreading across the continent. By initiating a citizen science monitoring program, we aimed to investigate the occurrence pattern of M. sculpturalis. Within only two years after starting the project, 111 new reports from Switzerland, Liechtenstein and Austria were recorded. Comparably to other European countries, the population progressed remarkably fast from year to year expanding its area geographically but also ecologically by increasing its altitudinal range. The distribution pattern indicates human assisted jump-dispersal travelling on the major traffic routes of central Europe.

Diversification in evolutionary arenas-Assessment and synthesis

Understanding how and why rates of evolutionary diversification vary is a key issue in evolutionary biology, ecology, and biogeography. Evolutionary rates are the net result of interacting processes summarized under concepts such as adaptive radiation and evolutionary stasis. Here, we review the central concepts in the evolutionary diversification literature and synthesize these into a simple, general framework for studying rates of diversification and quantifying their underlying dynamics, which can be applied across clades and regions, and across spatial and temporal scales. Our framework describes the diversification rate (d) as a function of the abiotic environment (a), the biotic environment (b), and clade-specific phenotypes or traits (c); thus, d ~ a,b,c. We refer to the four components (a-d) and their interactions collectively as the "Evolutionary Arena." We outline analytical approaches to this framework and present a case study on conifers, for which we parameterize the general model. We also discuss three conceptual examples: the Lupinus radiation in the Andes in the context of emerging ecological opportunity and fluctuating connectivity due to climatic oscillations; oceanic island radiations in the context of island formation and erosion; and biotically driven radiations of the Mediterranean orchid genus Ophrys. The results of the conifer case study are consistent with the long-standing scenario that low competition and high rates of niche evolution promote diversification. The conceptual examples illustrate how using the synthetic Evolutionary Arena framework helps to identify and structure future directions for research on evolutionary radiations. In this way, the Evolutionary Arena framework promotes a more general understanding of variation in evolutionary rates by making quantitative results comparable between case studies, thereby allowing new syntheses of evolutionary and ecological processes to emerge.

Accumulation over evolutionary time as a major cause of biodiversity hotspots in conifers

Biodiversity hotspots are important for understanding how areas of high species richness form, but disentangling the processes that produce them is difficult. We combine geographical ranges, phylogenetic relationships and trait data for 606 conifer species in order to explore the mechanisms underlying richness hotspot formation. We identify eight richness hotspots that overlap known centres of plant endemism and diversity, and find that conifer richness hotspots occur in mountainous areas within broader regions of long-term climate stability. Conifer hotspots are not unique in their species composition, traits or phylogenetic structure; however, a large percentage of their species are not restricted to hotspots and they rarely show either a preponderance of new radiating lineages or old relictual lineages. We suggest that conifer hotspots have primarily formed as a result of lineages accumulating over evolutionary time scales in stable mountainous areas rather than through high origination, preferential retention of relictual lineages or radiation of species with unique traits, although such processes may contribute to nuanced differences among hotspots. Conifers suggest that a simple accumulation of regional diversity can generate high species richness without additional processes and that geography rather than biology may play a primary role in hotspot formation.