Ira Loren Wiggins

Abstract Barriers to reproduction are often how progress in speciation is measured. Nonetheless, an unresolved question concerns the extent to which reproductive barriers diminish gene flow between incipient species. The Sierra Nevada foothill endemic Mimulus glaucescens and the widespread M. guttatus are considered distinct species based on striking differences in vegetative morphology, but barriers to reproduction have not been previously identified, nor has gene flow between species been characterized. Here, we examined 15 potential reproductive barriers within a Northern California area of broad sympatry. Most barriers, with the exception of ecogeographic isolation, were weak or absent, and total isolation for each species was incomplete. Population genomic analyses of range-wide and broadly sympatric accessions revealed extensive gene flow between these taxa, particularly in sympatry. Despite widespread introgression, Mimulus glaucescens, emerged as monophyletic and largely comprised a single ancestry that was found at intermediate frequency within M. guttatus. This result, along with observed ecological and phenotypic differentiation, suggests that natural selection may contribute to the maintenance of distinct phenotypic forms in the earliest stages of speciation. Integrating estimates of barrier strength with direct estimates of gene flow can strengthen a more nuanced interpretation of the process of speciation in natural communities.

Species climate niche models (CNMs) have been widely used for assessing climate change impact, developing conservation strategies and guiding assisted migration for adaptation to future climates. However, the CNMs built based on species occurrence data only reflect the species’ realized niche, which can overestimate the potential loss of suitable habitat of existing forests and underestimate the potential of assisted migration to mitigate climate change. In this study, we explored building a fundamental climate niche model using widely available species occurrence data with two important forest tree species, lodgepole pine (Pinus contorta Dougl. ex Loud.) and Douglas-fir (Pseudotsuga menziesii Franco.), which were introduced to many countries worldwide. We first compared and optimized three individual modeling techniques and their ensemble by adjusting the ratio of presence to absence (p/a) observations using an innovative approach to predict the realized climate niche of the two species. We then extended the realized climate niches to their fundamental niches by determining a new cut-off threshold based on species occurrence data beyond the native distributions. We found that the ensemble model comprising Random Forest and Maxent had the best performance and identified a common cut-off threshold of 0.3 for predicting the fundamental climate niches of the two species, which is likely applicable to other species. We then predicted the fundamental climate niches of the two species under current and future climate conditions. Our study demonstrated a novel approach for predicting species’ fundamental climate niche with high accuracy using only species occurrence data, including both presence and absence data points. It provided a new tool for assessing climate change impact on the future loss of existing forests and implementing assisted migration for better adapting to future climates.

Paleoclimate reconstructions can provide a window into the environmental conditions in Earth history when atmospheric carbon dioxide concentrations were higher than today. In the eastern USA, paleoclimate reconstructions are sparse, because terrestrial sedimentary deposits are rare. Despite this, the eastern USA has the largest population and population density in North America, and understanding the effects of current and future climate change is of vital importance. Here, we provide terrestrial paleoclimate reconstructions of the eastern USA from Miocene fossil floras. Additionally, we compare proxy paleoclimate reconstructions from the warmest period in the Miocene, the Miocene Climatic Optimum (MCO), to those of an MCO Earth System Model. Reconstructed Miocene temperatures and precipitation north of 35°N are higher than modern. In contrast, south of 35°N, temperatures and precipitation are similar to today, suggesting a poleward amplification effect in eastern North America. Reconstructed Miocene rainfall seasonality was predominantly higher than modern, regardless of latitude, indicating greater variability in intra-annual moisture transport. Reconstructed climates are almost uniformly in the temperate seasonal forest biome, but heterogeneity of specific forest types is evident. Reconstructed Miocene terrestrial temperatures from the eastern USA are lower than modeled temperatures and coeval Atlantic sea surface temperatures. However, reconstructed rainfall is consistent with modeled rainfall. Our results show that during the Miocene, climate was most different from modern in the northeastern states, and may suggest a drastic reduction in the meridional temperature gradient along the North American east coast compared to today.

Premise Researchers often use ecological niche models to predict where species might establish and persist under future or novel climate conditions. However, these predictive methods assume species have stable niches across time and space. Furthermore, ignoring the time of occurrence data can obscure important information about species reproduction and ultimately fitness. Here, we assess compare ecological niche models generated from full-year averages to seasonal models Methods In this study, we generate full-year and monthly ecological niche models for Capsella bursa-pastoris in Europe and North America to see if we can detect changes in the seasonal niche of the species after long-distance dispersal. Key Results We find full-year ecological niche models have low transferability across continents and there are continental differences in the climate conditions that influence the distribution of C. bursa-pastoris. Monthly models have greater predictive accuracy than full-year models in cooler seasons, but no monthly models are able to predict North American summer occurrences very well. Conclusions The relative predictive ability of European monthly models compared to North American monthly models suggests a change in the seasonal timing between the native range to the non-native range. These results highlight the utility of ecological niche models at finer temporal scales in predicting species distributions and unmasking subtle patterns of evolution.

Ophiuroidea is one of the most suitable marine groups for exploring diversity partitioning in the ocean due to its wide distribution and particular lifestyles. Nevertheless, diversity and its variation have yet to be investigated, and even basic information for large areas such as the eastern tropical Pacific (ETP) is still lacking. The present contribution explores α, β, and γ-diversity patterns of Ophiuroidea from the ETP at four spatial scales (Operational Geographic Units, Ecoregions, Provinces, and Realms). Based on literature records, databases, and scientific collections, an occurrence matrix was constructed for 69 shallow water (0–200 m) Ophiuroidea of the ETP (Mexico–Peru). Diversity evaluation based on rarefaction curves indicated that the observed richness tends to reach the asymptote. At the province and the ecoregion levels, β-diversity was the most important component explaining γ-diversity. The components that mainly contributed to the differentiation between provinces and ecoregions were the intersection of nestedness and β-diversity. PERMANOVA and SIMPER results showed that species composition presented significant differences at all spatial levels. The PCO ordination indicated that the first component (PCO1) explained the variation in species composition in a longitudinal gradient between coastal and oceanic ecoregions, while PCO2 showed a latitudinal gradient. The shade plot yielded three clusters (northern, southern, and widely distributed species). In general, α-diversity was explained by differences in sampling effort and methods; in contrast, β-diversity and its components were mainly explained by patterns and processes occurring at different spatial scales (provinces and ecoregions) such as oceanographic conditions, geographic extension, dispersal, and environmental heterogeneity. This work represents the first attempt to analyze the distribution patterns of shallow-water Ophiuroidea from the ETP.

Aim Museum and herbarium specimen records are frequently used to assess the conservation status of species and their responses to climate change. Typically, occurrences with imprecise geolocality information are discarded because they cannot be matched confidently to environmental conditions and are thus expected to increase uncertainty in downstream analyses. However, using only precisely georeferenced records risks undersampling of the environmental and geographical distributions of species. We present two related methods to allow the use of imprecisely georeferenced occurrences in biogeographical analysis. Innovation Our two procedures assign imprecise records to the (1) locations or (2) climates that are closest to the geographical or environmental centroid of the precise records of a species. For virtual species, including imprecise records alongside precise records improved the accuracy of ecological niche models projected to the present and the future, especially for species with c. 20 or fewer precise occurrences. Using only precise records underestimated loss of suitable habitat and overestimated the amount of suitable habitat in both the present and the future. Including imprecise records also improves estimates of niche breadth and extent of occurrence. An analysis of 44 species of North American Asclepias (Apocynaceae) yielded similar results. Main conclusions Existing studies examining the effects of spatial imprecision typically compare outcomes based on precise records against the same records with spatial error added to them. However, in real-world cases, analysts possess a mix of precise and imprecise records and must decide whether to retain or discard the latter. Discarding imprecise records can undersample the geographical and environmental distributions of species and lead to mis-estimation of responses to past and future climate change. Our method, for which we provide a software implementation in the enmSdmX package for R, is simple to use and can help leverage the large number of specimen records that are typically deemed “unusable” because of spatial imprecision in their geolocation.

The region of central Veracruz is considered a biodiversity hotspot due to its high species richness and environmental heterogeneity, but only 2% of this region is currently protected. This study aimed to assess the current protected area system’s effectiveness and to identify priority conservation areas for expanding the existing protected area system. We used the distribution models of 1186 species from three kingdoms (Animalia, Plantae, and Fungi) together with ZONATION software, a conservation planning tool, to determine areas that could help expand the current network of protected areas. We applied three different parametrizations (including only species, using the boundary quality penalty, and using corridor connectivity). We found that protecting an additional 15% of the area would increase, between 16.2% and 19.3%, the protection of the distribution area of all species. We propose that the regions with a consensus of the three parametrizations should be declared as new protected areas to expand 374 km2 to the 216 km2 already protected. Doing so would double the protected surface in central Veracruz. The priority areas identified in this study have more species richness, carbon stock values, natural vegetation cover, and less human impact index than the existing protected areas. If our identified priority areas are declared protected, we could expect a future recovery of endangered species populations for Veracruz. The proposed new protected areas are planned and designed as corridors connecting currently isolated protected areas to promote biodiversity protection.

The northern bog lemming Synaptomys borealis (NBL) is a rare small mammal that is undergoing a federal Species Status Assessment (SSA) under the US Endangered Species Act. Despite a wide North American distribution, very little is known about NBL dietary or habitat needs, both of which are germane to the resiliency of this species to climate change. To quantify diet composition of NBL in Alaska, we used DNA metabarcoding from 59 archived specimens to describe the taxonomic richness and relative abundance of foods in recent diets. DNA analyses revealed a broad diet composed of at least 110 families and 92 genera of bryophytes (mosses and liverworts), graminoids, fungi, forbs, and woody shrubs. Nine bryophyte genera and Carex sedges composed the largest portions of NBL diets. To quantify habitat preference, we intersected 467 georeferenced occurrence records of NBL in Alaska with remotely sensed land cover classes and used a compositional analysis framework that accounts for the relative abundance of land cover types. We did not detect significant habitat preferences for specific land cover types, although NBL frequently occurred in evergreen forest, woody wetlands, and adjacent to water. Our research highlights the importance of bryophytes, among a high diversity of dietary components, and describes NBL as boreal habitat generalists. Results will inform the current federal SSA by quantifying the extent to which ecological constraints are likely to affect NBL in a rapidly changing boreal environment.

Climate change has significantly affected global agricultural production, particularly in arid zones of Central Asia. Thus, we analyzed changes in the habitat suitability of cotton in Central Asia under various shared socioeconomic pathway (SSP) scenarios during 2021–2060. The results showed that the average minimum temperature in April, precipitation seasonality, and distance to rivers were the main environmental factors influencing the suitable distribution of cotton. Suitable habitats expanded toward the north and east, reaching a maximum net increase of 10.85 × 104 km2 under the SSP5–8.5 scenario during 2041–2060, while habitats in the southwestern area showed a contracting trend. The maximum decreased and increased habitats were concentrated at approximately 68°E and 87°E, respectively. In addition, their latitudinal distributions were concentrated at approximately 40°N and 44°N. The longitudinal and latitudinal dividing lines of increased and decreased habitats were 69°E and 41°N, respectively. Habitats at the same altitude showed an increasing trend, excluding the elevation range of 125–325 m. Habitat shifts could exacerbate spatial conflicts with forest/grassland and natural reserves. The maximum spatial overlap between them was observed under the SSP5–8.5 scenario during 2041–2060. These findings could provide scientific evidence for rational cotton cultivation planning in global arid zones.

Outside humans, true agriculture was previously thought to be restricted to social insects farming fungus. However, obligate farming of plants by ants was recently discovered in Fiji, prompting a re-examination of plant cultivation by ants. Here, we generate a database of plant cultivation by ants, identify three main types, and show that these interactions evolved primarily for shelter rather than food. We find that plant cultivation evolved at least 65 times independently for crops (~200 plant species), and 15 times in farmer lineages (~37 ant taxa) in the Neotropics and Asia/Australasia. Because of their high evolutionary replication, and variation in partner dependence, these systems are powerful models to unveil the steps in the evolution and ecology of insect agriculture.