Uncovering Broad Macroecological Patterns by Comparing the Shape of Species' Distributions along Environmental Gradients.

AbstractSpecies' distributions can take many different forms. For example, fat-tailed or skewed distributions are very common in nature, as these can naturally emerge as a result of individual variability and asymmetric environmental tolerances, respectively. Studying the basic shape of distributions can teach us a lot about the ways climatic processes and historical contingencies shape ecological communities. Yet we still lack a general understanding of how their shapes and properties compare to each other along gradients. Here, we use Bayesian nonlinear models to quantify range shape properties in empirical plant distributions. With this approach, we are able to distil the shape of plant distributions and compare them along gradients and across species. Studying the relationship between distribution properties, we revealed the existence of broad macroecological patterns along environmental gradients-such as those expected from Rapoport's rule and the abiotic stress limitation hypothesis. We also find that some aspects of the shape of observed ranges-such as kurtosis and skewness of the distributions-could be intrinsic properties of species or the result of their historical contexts. Overall, our modeling approach and results untangle the general shape of plant distributions and provide a mapping of how this changes along environmental gradients.

Which contributes more to the relict flora distribution pattern in East Asia, geographical processes or climate change? New evidence from the phylogeography of Rehderodendron kwangtungense.

BACKGROUND: Relict species are important for enhancing the understanding of modern biogeographic distribution patterns. Although both geological and climatic changes since the Cenozoic have affected the relict flora in East Asia, the contributions of geographical processes remain unclear. In this study, we employed restriction-site associated DNA sequencing (RAD-seq) and shallow genome sequencing data, in conjunction with ecological niche modeling (ENM), to investigate the spatial genetic patterns and population differentiation history of the relict species Rehderodendron kwangtungense Chun. RESULTS: A total of 138 individuals from 16 populations were collected, largely covering the natural distribution of R. kwangtungense. The genetic diversity within the R. kwangtungense populations was extremely low (HO = 0.048 ± 0.019; HE = 0.033 ± 0.011). Mantel tests revealed isolation-by-distance pattern (R2 = 0.38, P < 0.001), and AMOVA analysis showed that the genetic variation of R. kwangtungense occurs mainly between populations (86.88%, K = 7). Between 23 and 21 Ma, R. kwangtungense underwent a period of rapid differentiation that coincided with the rise of the Himalayas and the establishment of the East Asian monsoon. According to ENM and population demographic history, the suitable area and effective population size of R. kwangtungense decreased sharply during the glacial period and expanded after the last glacial maximum (LGM). CONCLUSION: Our study shows that the distribution pattern of southern China mountain relict flora may have developed during the panplain stage between the middle Oligocene and the early Miocene. Then, the flora later fragmented under the force of orogenesis, including intermittent uplift during the Cenozoic Himalayan orogeny and the formation of abundant rainfall associated with the East Asian monsoon. The findings emphasized the predominant role of geographical processes in shaping relict plant distribution patterns.

The potential range of west Asian apple species Malus orientalis Uglitzk. under climate change.

The wild relatives of cultivated apples would be an ideal source of diversity for breeding new varieties, which could potentially grow in diverse habitats shaped by climate change. However, there is still a lack of knowledge about the potential distribution of these species. The aim of the presented work was the understand the impacts of climate change on the potential distribution and habitat fragmentation of Caucasian crab apple (Malus orientalis Uglitzk.) and the designation of areas of high interest according to climatic conditions. We used the MaxEnt models and Morphological-Spatial Analysis (MSPA) to evaluate the potential distribution, suitability changes, habitat fragmentation, and connectivity throughout the species range in Turkey, Armenia, Georgia, Russia, and Iran. The results revealed that the potentially suitable range of M. orientalis encompasses 858,877 km², 635,279 km² and 456,795 km² under the present, RCP4.5 and RCP8.5 scenario, respectively. The range fragmentation analysis demonstrated a notable shift in the edge/core ratio, which increased from 50.95% in the current scenario to even 67.70% in the future. The northern part of the range (Armenia, northern Georgia, southern Russia), as well as the central and western parts of Hyrcania will be a core of the species range with suitable habitats and a high connectivity between M. orientalis populations and could work as major refugia for the studied species. However, in the Zagros and central Turkey, the potential range will shrink due to the lack of suitable climatic conditions, and the edge/core ratio will grow. In the southern part of the range, a decline of M. orientalis habitats is expected due to changing climatic conditions. The future outlook suggests that the Hyrcanian forest and the Caucasus region could serve as important refuges for M. orientalis. This study helps to understand spatial changes in species' range in response to climate change and can help develop conservation strategies. This is all the more important given the species' potential use in future breeding programs aimed at enriching the gene pool of cultivated apple varieties.

Ecological and metabolic implications of the nurse effect of Maihueniopsis camachoi in the Atacama Desert.

Plant-plant positive interactions are key drivers of community structure. Yet, the underlying molecular mechanisms of facilitation processes remain unexplored. We investigated the 'nursing' effect of Maihueniopsis camachoi, a cactus that thrives in the Atacama Desert between c. 2800 and 3800 m above sea level. We hypothesised that an important protective factor is thermal amelioration of less cold-tolerant species with a corresponding impact on molecular phenotypes. To test this hypothesis, we compared plant cover and temperatures within the cactus foliage with open areas and modelled the effect of temperatures on plant distribution. We combined eco-metabolomics and machine learning to test the molecular consequences of this association. Multiple species benefited from the interaction with M. camachoi. A conspicuous example was the extended distribution of Atriplex imbricata to colder elevations in association with M. camachoi (400 m higher as compared to plants in open areas). Metabolomics identified 93 biochemical markers predicting the interaction status of A. imbricata with 79% accuracy, independently of year. These findings place M. camachoi as a key species in Atacama plant communities, driving local biodiversity with an impact on molecular phenotypes of nursed species. Our results support the stress-gradient hypothesis and provide pioneer insights into the metabolic consequences of facilitation.

Soil properties constrain predicted poleward migration of plants under climate change.

Many plant species are predicted to migrate poleward in response to climate change. Species distribution models (SDMs) have been widely used to quantify future suitable habitats, but they often neglect soil properties, despite the importance of soil for plant fitness. As soil properties often change along latitudinal gradients, higher-latitude soils might be more or less suitable than average conditions within the current ranges of species, thereby accelerating or slowing potential poleward migration. In this study, we built three SDMs - one with only climate predictors, one with only soil predictors, and one with both - for each of 1870 plant species in Eastern North America, in order to investigate the relative importance of soil properties in determining plant distributions and poleward shifts under climate change. While climate variables were the most important predictors, soil properties also had a substantial influence on continental-scale plant distributions. Under future climate scenarios, models including soil predicted much smaller northward shifts in distributions than climate-only models (c. 40% reduction). Our findings strongly suggest that high-latitude soils are likely to impede ongoing plant migration, and they highlight the necessity of incorporating soil properties into models and predictions for plant distributions and migration under environmental change.

Effective dispersal patterns in prairie plant species across human-modified landscapes.

Effective dispersal among plant populations is dependent on vector behaviour, landscape features and availability of adequate habitats. To capture landscape feature effects on dispersal, studies must be conducted at scales reflecting single-generation dispersal events (mesoscale). Many studies are conducted at large scales where genetic differentiation is due to dispersal occurring over multiple generations, making it difficult to interpret the effects of specific landscape features on vector behaviour. Genetic structure at the mesoscale may be determined by ecological and evolutionary processes, such as the consequences of vector behaviour on patterns of gene flow. We used chloroplast haplotypes and nuclear genome SNP surveys to identify landscape features influencing seed and pollen dispersal at a mesoscale within the Rogue River Valley in southern Oregon. We evaluated biotic and abiotic vector behaviour by contrasting two annual species with differing dispersal mechanisms; Achyrachaena mollis (Asteraceae) is a self-pollinating and anemochoric species, and Plectritis congesta (Caprifoliaceae) is biotically pollinated with barochoric seeds. Using landscape genetics methods, we identified features of the study region that conduct or restrict dispersal. We found chloroplast haplotypes were indicative of historic patterns of gene flow prior to human modification of landscapes. Seed dispersal of A. mollis was best supported by models of isolation by distance, while seed-driven gene flow of P. congesta was determined by the distribution of preserved natural spaces and quality habitat. Nuclear genetic structure was driven by both pollen and seed dispersal, and both species responded to contemporary landscape changes, such as urban and agricultural conversion, and habitat availability.

The role of Pleistocene dispersal in shaping species richness of sky island wintergreens from the Himalaya-Hengduan Mountains.

In addition to topography and climate, biogeographic dispersal has been considered to influence plant diversity in the Himalaya-Hengduan Mountains (HHM), yet, the mode and tempo of sky island dispersal and its influence on species richness has been little explored. Through phylogenetic analysis of Gaultheria ser. Trichophyllae, a sky island alpine clade within the HHM, we test the hypothesis that dispersal has affected current local species richness. We inferred the dynamics of biogeographic dispersal with correlation tests on direction, distance, occurrence time, and regional species richness. We found that G. ser. Trichophyllae originated at the end of the Miocene and mostly dispersed toward higher longitudes (eastward). In particular, shorter intra-regional eastward dispersals and longer inter-regional westward dispersals were most frequently observed. We detected a prevalence of eastward intra-region dispersals in both glacial periods and interglacials. These dispersals may have been facilitated by the reorganization of paleo-drainages and monsoon intensification through time. We suggest that the timing of dispersal corresponding to glacial periods and the prevalence of intra-region dispersal, rather than dispersal frequency, most influenced the pattern of species richness of G. ser. Trichophyllae. This study facilitates a more comprehensive understanding of biodiversity in the sky islands within the HHM.

Evaluating the definition and distribution of spring ephemeral wildflowers in eastern North America.

PREMISE: The herbaceous layer accounts for the majority of plant biodiversity in eastern North American forests, encompassing substantial variation in life history strategy and function. One group of early-season herbaceous understory species, colloquially referred to as spring ephemeral wildflowers, are ecologically and culturally important, but little is known about the prevalence and biogeographic patterns of the spring ephemeral strategy. METHODS: We used observations collected by the Global Biodiversity Information Facility (GBIF) to quantify the ephemerality of 559 understory forb species across eastern North America and classify them according to a continuous ephemerality index (ranging from 0 = never ephemeral to 1 = always ephemeral). We then used this information to model where ephemeral forbs were most common across the landscape with the goal of identifying geographic and environmental drivers important to their distributions and ranges. RESULTS: Only 3.4% of all understory wildflower species were spring ephemerals in all parts of their range, and 18.4% (103 species) were ephemeral in at least part of their range. Spring ephemerals peaked in absolute species richness and relative proportion at mid latitudes. CONCLUSIONS: Spring ephemeral phenology is an important shade-avoidance strategy for a large segment of the total understory species in temperate deciduous forests. In North America, the strategy is relatively most important for forest understories at mid latitudes. The definitions of spring ephemerality we provide here serve as an important ecological context for conservation priorities and to evaluate responses of this biodiverse group to future environmental change.

Ancient plant DNA reveals High Arctic greening during the Last Interglacial.

Summer warming is driving a greening trend across the Arctic, with the potential for large-scale amplification of climate change due to vegetation-related feedbacks [Pearson et al., Nat. Clim. Chang. (3), 673-677 (2013)]. Because observational records are sparse and temporally limited, past episodes of Arctic warming can help elucidate the magnitude of vegetation response to temperature change. The Last Interglacial ([LIG], 129,000 to 116,000 y ago) was the most recent episode of Arctic warming on par with predicted 21st century temperature change [Otto-Bliesner et al., Philos. Trans. A Math. Phys. Eng. Sci. (371), 20130097 (2013) and Post et al., SciAdv (5), eaaw9883 (2019)]. However, high-latitude terrestrial records from this period are rare, so LIG vegetation distributions are incompletely known. Pollen-based vegetation reconstructions can be biased by long-distance pollen transport, further obscuring the paleoenvironmental record. Here, we present a LIG vegetation record based on ancient DNA in lake sediment and compare it with fossil pollen. Comprehensive plant community reconstructions through the last and current interglacial (the Holocene) on Baffin Island, Arctic Canada, reveal coherent climate-driven community shifts across both interglacials. Peak LIG warmth featured a ∼400-km northward range shift of dwarf birch, a key woody shrub that is again expanding northward. Greening of the High Arctic-documented here by multiple proxies-likely represented a strong positive feedback on high-latitude LIG warming. Authenticated ancient DNA from this lake sediment also extends the useful preservation window for the technique and highlights the utility of combining traditional and molecular approaches for gleaning paleoenvironmental insights to better anticipate a warmer future.

A multimillion-year-old record of Greenland vegetation and glacial history preserved in sediment beneath 1.4 km of ice at Camp Century.

Understanding the history of the Greenland Ice Sheet (GrIS) is critical for determining its sensitivity to warming and contribution to sea level; however, that history is poorly known before the last interglacial. Most knowledge comes from interpretation of marine sediment, an indirect record of past ice-sheet extent and behavior. Subglacial sediment and rock, retrieved at the base of ice cores, provide terrestrial evidence for GrIS behavior during the Pleistocene. Here, we use multiple methods to determine GrIS history from subglacial sediment at the base of the Camp Century ice core collected in 1966. This material contains a stratigraphic record of glaciation and vegetation in northwestern Greenland spanning the Pleistocene. Enriched stable isotopes of pore-ice suggest precipitation at lower elevations implying ice-sheet absence. Plant macrofossils and biomarkers in the sediment indicate that paleo-ecosystems from previous interglacial periods are preserved beneath the GrIS. Cosmogenic 26Al/10Be and luminescence data bracket the burial of the lower-most sediment between <3.2 ± 0.4 Ma and >0.7 to 1.4 Ma. In the upper-most sediment, cosmogenic 26Al/10Be data require exposure within the last 1.0 ± 0.1 My. The unique subglacial sedimentary record from Camp Century documents at least two episodes of ice-free, vegetated conditions, each followed by glaciation. The lower sediment derives from an Early Pleistocene GrIS advance. 26Al/10Be ratios in the upper-most sediment match those in subglacial bedrock from central Greenland, suggesting similar ice-cover histories across the GrIS. We conclude that the GrIS persisted through much of the Pleistocene but melted and reformed at least once since 1.1 Ma.

Modeling the impact of sea level rise on endangered deer habitat.

Numerous unique flora and fauna inhabit the Lower Florida Keys, including the endangered Florida Key deer, found nowhere else. In this vulnerable habitat of flat islands with low elevation, accelerated sea level rise poses a threat. Predicting the impact of sea level rise on vegetation and wildlife is crucial. This study used 5 Intergovernmental Panel on Climate Change (IPCC) sea level rise scenarios to assess their effects on No Name Key, Florida. The goal was to estimate changes in the Florida Key deer population relative to sea level rise using a lidar-derived elevation data and a vegetation map. The method used 2 cases to model the sea level rise impact. In Case 1, total non-submerged area at current sea level was determined. Using 5 IPCC scenarios, a new total non-submerged land area was estimated, and deer numbers were predicted for each scenario. In Case 2, upward migration of coastal vegetation combined with the coastal squeeze process was modeled. A distinct elevation range for each vegetation type at the current sea level was determined. Vegetation ranges were redistributed based on respective elevation ranges in the sea level rise scenarios. Areas for each vegetation type were recalculated, and Key deer numbers were estimated for each sea level rise scenario. Results under the worst emission scenario showed the following: (1) for case 1, the land area was reduced to 30 % of the current land area, corresponding to having about 27 deer, and (2) for case 2, the land area was reduced to 70 % of the current land area, having about 54 deer on No Name Key. The results indicated reduced non-submerged land area and less upland vegetation, particularly hardwoods/hammocks, by the year 2100. As less land area is available, a decline in Key deer population is expected as sea levels rise. Since Key deer favor upland vegetation, habitat affected by sea level rise will likely support a smaller deer population. The findings emphasize the need for precise, timely predictions of sea level rise impacts and long-term conservation strategies. Specifically designed measures are required to protect and maintain endangered wildlife, such as the Florida Key deer, residing on these vulnerable islands.

The genome and preliminary single-nuclei transcriptome of Lemna minuta reveals mechanisms of invasiveness.

The ability to trace every cell in some model organisms has led to the fundamental understanding of development and cellular function. However, in plants the complexity of cell number, organ size, and developmental time makes this a challenge even in the diminutive model plant Arabidopsis (Arabidopsis thaliana). Duckweed, basal nongrass aquatic monocots, provide an opportunity to follow every cell of an entire plant due to their small size, reduced body plan, and fast clonal growth habit. Here we present a chromosome-resolved genome for the highly invasive Lesser Duckweed (Lemna minuta) and generate a preliminary cell atlas leveraging low cell coverage single-nuclei sequencing. We resolved the 360 megabase genome into 21 chromosomes, revealing a core nonredundant gene set with only the ancient tau whole-genome duplication shared with all monocots, and paralog expansion as a result of tandem duplications related to phytoremediation. Leveraging SMARTseq2 single-nuclei sequencing, which provided higher gene coverage yet lower cell count, we profiled 269 nuclei covering 36.9% (8,457) of the L. minuta transcriptome. Since molecular validation was not possible in this nonmodel plant, we leveraged gene orthology with model organism single-cell expression datasets, gene ontology, and cell trajectory analysis to define putative cell types. We found that the tissue that we computationally defined as mesophyll expressed high levels of elemental transport genes consistent with this tissue playing a role in L. minuta wastewater detoxification. The L. minuta genome and preliminary cell map provide a paradigm to decipher developmental genes and pathways for an entire plant.

A climatic dipole drives short- and long-term patterns of postfire forest recovery in the western United States.

Researchers are increasingly examining patterns and drivers of postfire forest recovery amid growing concern that climate change and intensifying fires will trigger ecosystem transformations. Diminished seed availability and postfire drought have emerged as key constraints on conifer recruitment. However, the spatial and temporal extent to which recurring modes of climatic variability shape patterns of postfire recovery remain largely unexplored. Here, we identify a north-south dipole in annual climatic moisture deficit anomalies across the Interior West of the US and characterize its influence on forest recovery from fire. We use annually resolved establishment models from dendrochronological records to correlate this climatic dipole with short-term postfire juvenile recruitment. We also examine longer-term recovery trajectories using Forest Inventory and Analysis data from 989 burned plots. We show that annual postfire ponderosa pine recruitment probabilities in the northern Rocky Mountains (NR) and the southwestern US (SW) track the strength of the dipole, while declining overall due to increasing aridity. This indicates that divergent recovery trajectories may be triggered concurrently across large spatial scales: favorable conditions in the SW can correspond to drought in the NR that inhibits ponderosa pine establishment, and vice versa. The imprint of this climatic dipole is manifest for years postfire, as evidenced by dampened long-term likelihoods of juvenile ponderosa pine presence in areas that experienced postfire drought. These findings underscore the importance of climatic variability at multiple spatiotemporal scales in driving cross-regional patterns of forest recovery and have implications for understanding ecosystem transformations and species range dynamics under global change.

The spatiotemporal spread of human migrations during the European Holocene.

The European continent was subject to two major migrations of peoples during the Holocene: the northwestward movement of Anatolian farmer populations during the Neolithic and the westward movement of Yamnaya steppe peoples during the Bronze Age. These movements changed the genetic composition of the continent's inhabitants. The Holocene was also characterized by major changes in vegetation composition, which altered the environment occupied by the original hunter-gatherer populations. We aim to test to what extent vegetation change through time is associated with changes in population composition as a consequence of these migrations, or with changes in climate. Using ancient DNA in combination with geostatistical techniques, we produce detailed maps of ancient population movements, which allow us to visualize how these migrations unfolded through time and space. We find that the spread of Neolithic farmer ancestry had a two-pronged wavefront, in agreement with similar findings on the cultural spread of farming from radiocarbon-dated archaeological sites. This movement, however, did not have a strong association with changes in the vegetational landscape. In contrast, the Yamnaya migration speed was at least twice as fast and coincided with a reduction in the amount of broad-leaf forest and an increase in the amount of pasture and natural grasslands in the continent. We demonstrate the utility of integrating ancient genomes with archaeometric datasets in a spatiotemporal statistical framework, which we foresee will enable future studies of ancient populations' movements, and their putative effects on local fauna and flora.

Diversity and distribution of the genus Tetmemorus (Desmidiaceae, Zygnematophyceae) in Brazil.

In the present study, a taxonomic review was conducted on representatives of the genus Tetmemorus (Desmidiaceae, Zygnematophyceae) documented within Brazilian territory. This review involved compiling data from the literature and analyzing samples collected throughout the Bahia State, updating our knowledge about this genus in Brazil. For each identified taxon, we provided information such as description, distribution across biomes and states, watersheds, ecological aspects (including habitat and community types), a list of examined (and excluded) materials, and taxonomic comments. Additionally, a taxonomic key for all species reported in Brazil was provided. Through this comprehensive review, we identified a total of eight Tetmemorus taxa occurring in Brazilian territory, comprising five species (T. brebissonii, T. furcatus, T. granulatus, T. laevis, T. planctonicus) and three non-typical varieties (T. brebissonii var. minor, T. laevis var. borgei, T. laevis var. minutus).

Pollen-based reconstruction reveals the impact of the onset of agriculture on plant functional trait composition.

The onset of agriculture improved the capacity of ecosystems to produce food, but inadvertently altered other vital ecosystem functions. Plant traits play a central role in determining ecosystem properties, therefore we investigated how the onset of agriculture in Europe changed plant trait composition using 78 pollen records. Using a novel Bayesian approach for reconstructing plant trait composition from pollen records, we provide a robust method that can account for trait variability within pollen types. We estimate an overall four-fold decrease in plant size through agriculture and associated decreases in leaf and seed size. We show an increase in niche space towards the resource-acquisitive end of the leaf economic spectrum. Decreases in leaf phosphorus might have been caused by nutrient depletion through grazing and burning. Our results show that agriculture, from its start, has likely been gradually impacting biogeochemical cycles through altered vegetation composition.

Maintenance of Adaptive Dynamics and No Detectable Load in a Range-Edge Outcrossing Plant Population.

During range expansion, edge populations are expected to face increased genetic drift, which in turn can alter and potentially compromise adaptive dynamics, preventing the removal of deleterious mutations and slowing down adaptation. Here, we contrast populations of the European subspecies Arabidopsis lyrata ssp. petraea, which expanded its Northern range after the last glaciation. We document a sharp decline in effective population size in the range-edge population and observe that nonsynonymous variants segregate at higher frequencies. We detect a 4.9% excess of derived nonsynonymous variants per individual in the range-edge population, suggesting an increase of the genomic burden of deleterious mutations. Inference of the fitness effects of mutations and modeling of allele frequencies under the explicit demographic history of each population predicts a depletion of rare deleterious variants in the range-edge population, but an enrichment for fixed ones, consistent with the bottleneck effect. However, the demographic history of the range-edge population predicts a small net decrease in per-individual fitness. Consistent with this prediction, the range-edge population is not impaired in its growth and survival measured in a common garden experiment. We further observe that the allelic diversity at the self-incompatibility locus, which ensures strict outcrossing and evolves under negative frequency-dependent selection, has remained unchanged. Genomic footprints indicative of selective sweeps are broader in the Northern population but not less frequent. We conclude that the outcrossing species A. lyrata ssp. petraea shows a strong resilience to the effect of range expansion.

Mechanisms underlying interaction frequencies and robustness in a novel seed dispersal network: lessons for restoration.

As human-caused extinctions and invasions accumulate across the planet, understanding the processes governing ecological functions mediated by species interactions, and anticipating the effect of species loss on such functions become increasingly urgent. In seed dispersal networks, the mechanisms that influence interaction frequencies may also influence the capacity of a species to switch to alternative partners (rewiring), influencing network robustness. Studying seed dispersal interactions in novel ecosystems on O'ahu island, Hawai'i, we test whether the same mechanisms defining interaction frequencies can regulate rewiring and increase network robustness to simulated species extinctions. We found that spatial and temporal overlaps were the primary mechanisms underlying interaction frequencies, and the loss of the more connected species affected networks to a greater extent. Further, rewiring increased network robustness, and morphological matching and spatial and temporal overlaps between partners were more influential on network robustness than species abundances. We argue that to achieve self-sustaining ecosystems, restoration initiatives can consider optimal morphological matching and spatial and temporal overlaps between consumers and resources to maximize chances of native plant dispersal. Specifically, restoration initiatives may benefit from replacing invasive species with native species possessing characteristics that promote frequent interactions and increase the probability of rewiring (such as long fruiting periods, small seeds and broad distributions).

An overlooked dispersal route of Cardueae (Asteraceae) from the Mediterranean to East Asia revealed by phylogenomic and biogeographical analyses of Atractylodes.

BACKGROUND AND AIMS: The East Asian-Tethyan disjunction pattern and its mechanisms of formation have long been of interest to researchers. Here, we studied the biogeographical history of Asteraceae tribe Cardueae, with a particular focus on the temperate East Asian genus Atractylodes DC., to understand the role of tectonic and climatic events in driving the diversification and disjunctions of the genus. METHODS: A total of 76 samples of Atractylodes from 36 locations were collected for RAD-sequencing. Three single nucleotide polymorphism (SNP) datasets based on different filtering strategies were used for phylogenetic analyses. Molecular dating and ancestral distribution reconstruction were performed using both chloroplast DNA sequences (127 Cardueae samples) and SNP (36 Atractylodes samples) datasets. KEY RESULTS: Six species of Atractylodes were well resolved as individually monophyletic, although some introgression was identified among accessions of A. chinensis, A. lancea and A. koreana. Dispersal of the subtribe Carlininae from the Mediterranean to East Asia occurred after divergence between Atractylodes and Carlina L. + Atractylis L. + Thevenotia DC. at ~31.57 Ma, resulting in an East Asian-Tethyan disjunction. Diversification of Atractylodes in East Asia mainly occurred from the Late Miocene to the Early Pleistocene. CONCLUSIONS: Aridification of Asia and the closure of the Turgai Strait in the Late Oligocene promoted the dispersal of Cardueae from the Mediterranean to East China. Subsequent uplift of the Qinghai-Tibet Plateau as well as changes in Asian monsoon systems resulted in an East Asian-Tethyan disjunction between Atractylodes and Carlina + Atractylis + Thevenotia. In addition, Late Miocene to Quaternary climates and sea level fluctuations played major roles in the diversification of Atractylodes. Through this study of different taxonomic levels using genomic data, we have revealed an overlooked dispersal route between the Mediterranean and far East Asia (Japan/Korea) via Central Asia and East China.

Invasive grasses increase fire occurrence and frequency across US ecoregions.

Fire-prone invasive grasses create novel ecosystem threats by increasing fine-fuel loads and continuity, which can alter fire regimes. While the existence of an invasive grass-fire cycle is well known, evidence of altered fire regimes is typically based on local-scale studies or expert knowledge. Here, we quantify the effects of 12 nonnative, invasive grasses on fire occurrence, size, and frequency across 29 US ecoregions encompassing more than one third of the conterminous United States. These 12 grass species promote fire locally and have extensive spatial records of abundant infestations. We combined agency and satellite fire data with records of abundant grass invasion to test for differences in fire regimes between invaded and nearby "uninvaded" habitat. Additionally, we assessed whether invasive grass presence is a significant predictor of altered fire by modeling fire occurrence, size, and frequency as a function of grass invasion, in addition to anthropogenic and ecological covariates relevant to fire. Eight species showed significantly higher fire-occurrence rates, which more than tripled for Schismus barbatus and Pennisetum ciliare. Six species demonstrated significantly higher mean fire frequency, which more than doubled for Neyraudia reynaudiana and Pennisetum ciliare Grass invasion was significant in fire occurrence and frequency models, but not in fire-size models. The significant differences in fire regimes, coupled with the importance of grass invasion in modeling these differences, suggest that invasive grasses alter US fire regimes at regional scales. As concern about US wildfires grows, accounting for fire-promoting invasive grasses will be imperative for effectively managing ecosystems.