Population genomic and historical analysis suggests a global invasion by bridgehead processes in Mimulus guttatus

Lightsheet microscopy integrates single-cell optical visco-elastography and fluorescence cytometry of 3D live tissues

Most common cytometry methods, including flow cytometry, observe suspended or fixed cells and cannot evaluate their structural roles in 3D tissues. However, cellular physical interactions are critical in physiological, developmental, and pathological processes. Here, we present a novel optical visco-elastography that characterizes single-cellular physical interactions by applying in-situ micro-mechanical perturbation to live microtissues under 3D lightsheet microscopy. The 4D digital image correlation (DIC) analysis of ~20,000 nodes tracked the compressive deformation of 3D tissues containing ~500 cells. The computational 3D image segmentation allowed cell-by-cell qualitative observation and statistical analysis, directly correlating multi-channel fluorescence and viscoelasticity. To represent epithelia-stroma interactions, we used a 3D organoid model of maternal-fetal interface and visualized solid-like, well-aligned displacement and liquid-like random motion between individual cells. The statistical analysis through our unique cytometry confirmed that endometrial stromal fibroblasts stiffen in response to decidualization. Moreover, we demonstrated in the 3D model that interaction with placental extravillous trophoblasts partially reverses the attained stiffness, which was supported by the gene expression analysis. Placentation shares critical cellular and molecular significance with various fundamental biological events such as cancer metastasis, wound healing, and gastrulation. Our analysis confirmed existing beliefs and discovered new insights, proving the broad applicability of our method.

Temporal matches between monarch butterfly and milkweed population changes over the past 25,000 years

In intimate ecological interactions, the interdependency of species may result in correlated demographic histories. For species of conservation concern, understanding the long-term dynamics of such interactions may shed light on the drivers of population decline. Here, we address the demographic history of the monarch butterfly, Danaus plexippus, and its dominant host plant, the common milkweed Asclepias syriaca (A. syriaca), using broad-scale sampling and genomic inference. Because genetic resources for milkweed have lagged behind those for monarchs, we first release a chromosome-level genome assembly and annotation for common milkweed. Next, we show that despite its enormous geographic range across eastern North America, A. syriaca is best characterized as a single, roughly panmictic population. Using approximate Bayesian computation with random forests (ABC-RF), a machine learning method for reconstructing demographic histories, we show that both monarchs and milkweed experienced population expansion during the most recent recession of North American glaciers 10,000-20,000 years ago. Our data also identify concurrent population expansions in both species during the large-scale clearing of eastern forests (∼200 years ago). Finally, we find no evidence that either species experienced a reduction in effective population size over the past 75 years. Thus, the well-documented decline of monarch abundance over the past 40 years is not visible in our genomic dataset, reflecting a possible mismatch of the overwintering census population to effective population size in this species.

Demographic consequences of an extreme heat wave are mitigated by spatial heterogeneity in an annual monkeyflower

Heat waves are becoming more frequent and intense with climate change, but the demographic and evolutionary consequences of heat waves are rarely investigated in herbaceous plant species. We examine the consequences of a short but extreme heat wave in Oregon populations of the common yellow monkeyflower (Mimulus guttatus) by leveraging a common garden experiment planted with range-wide populations and observational studies of 11 local populations. In the common garden, 89% of seedlings died during the heat wave including >96% of seedlings from geographically local populations. Some populations from hotter and drier environments had higher fitness, however, others from comparable environments performed poorly. Observational studies of local natural populations drastically differed in the consequences of the heat wave-one population was completely extirpated and nearly half had a >50% decrease in fitness. However, a few populations had greater fitness during the heat wave year. Differences in mortality corresponded to the impact of the heat wave on soil moisture-retention of soil moisture throughout the heat wave led to greater survivorship. Our results suggest that not all populations experience the same intensity or degree of mortality during extreme events and such heterogeneity could be important for genetic rescue or to facilitate the distribution of adaptive variants throughout the region.

Trait evolution during a rapid global weed invasion despite little genetic differentiation

Invasive species often possess a great capacity to adapt to novel environments in the form of spatial trait variation, as a result of varying selection regimes, genetic drift, or plasticity. We explored the geographic differentiation in several phenotypic traits related to plant growth, reproduction, and defense in the highly invasive Centaurea solstitialis by measuring neutral genetic differentiation (F _ST), and comparing it with phenotypic differentiation (P _ST), in a common garden experiment in individuals originating from regions representing the species distribution across five continents. Native plants were more fecund than non-native plants, but the latter displayed considerably larger seed mass. We found indication of divergent selection for these two reproductive traits but little overall genetic differentiation between native and non-native ranges. The native versus invasive P _ST-F _ST comparisons demonstrated that, in several invasive regions, seed mass had increased proportionally more than the genetic differentiation. Traits displayed different associations with climate variables in different regions. Both capitula numbers and seed mass were associated with winter temperature and precipitation and summer aridity in some regions. Overall, our study suggests that rapid evolution has accompanied invasive success of C. solstitialis and provides new insights into traits and their genetic bases that can contribute to fitness advantages in non-native populations.

Out of East Asia: Early Warning of the Possible Invasion of the Important Bean Pest Stalk-Eyed Seed Bug Chauliops fallax (Heteroptera: Malcidae: Chauliopinae)

The short stay at the beginning of the invasion process is a critical time for invasive species identification and preventing invasive species from developing a wider distribution and significant economic impact. The stalk-eyed seed bug Chauliops fallax is an important agricultural pest of soybean and was first reported to occur outside East Asia. Here, we reported the native evolutionary history, recent invasion history, and potential invasion threats of C. fallax for the first time based on population genetic methods and ecological niche modelling. The results showed that four native East Asian genetic groups (EA, WE, TL, and XZ) were well supported, showing an east-west differentiation pattern consistent with the geographical characteristics of three-step landforms in China. Two main haplotypes existed: Hap1 might have experienced a rapid northwards expansion process after the LGM period, and Hap5 reflected local adaptation to the environment in southeastern China. The Kashmir sample was found to come from the recent invasion of populations in the coastal areas of southern China. Ecological niche modelling results suggested that North America has a high risk of invasion, which might pose a serious threat to local soybean production. In addition, with future global warming, the suitable habitat in Asia will move towards the higher latitude region and gradually deviate from the soybean planting area, which indicates the threat of C. fallax to soybean production in Asia will decrease in the future. The results could provide new insights into the monitoring and management of this agricultural pest in the early invasion stage.

Contrasting continental patterns of adaptive population divergence in the holarctic ectomycorrhizal fungus Boletus edulis

In the hyperdiverse fungi, the process of speciation is virtually unknown, including for the > 20 000 species of ectomycorrhizal mutualists. To understand this process, we investigated patterns of genome-wide differentiation in the ectomycorrhizal porcini mushroom, Boletus edulis, a globally distributed species complex with broad ecological amplitude. By whole-genome sequencing 160 individuals from across the Northern Hemisphere, we genotyped 792 923 single nucleotide polymorphisms to characterize patterns of genome-wide differentiation and to identify the adaptive processes shaping global population structure. We show that B. edulis exhibits contrasting patterns of genomic divergence between continents, with multiple lineages present across North America, while a single lineage dominates Europe. These geographical lineages are inferred to have diverged 1.62-2.66 million years ago, during a period of climatic upheaval and the onset of glaciation in the Pliocene-Pleistocene boundary. High levels of genomic differentiation were observed among lineages despite evidence of substantial and ongoing introgression. Genome scans, demographic inference, and ecological niche models suggest that genomic differentiation is maintained by environmental adaptation, not physical isolation. Our study uncovers striking patterns of genome-wide differentiation on a global scale and emphasizes the importance of local adaptation and ecologically mediated divergence, rather than prezygotic barriers such as allopatry or genomic incompatibility, in fungal population differentiation.

Genetic diversity loss in the Anthropocene

Anthropogenic habitat loss and climate change are reducing species' geographic ranges, increasing extinction risk and losses of species' genetic diversity. Although preserving genetic diversity is key to maintaining species' adaptability, we lack predictive tools and global estimates of genetic diversity loss across ecosystems. We introduce a mathematical framework that bridges biodiversity theory and population genetics to understand the loss of naturally occurring DNA mutations with decreasing habitat. By analyzing genomic variation of 10,095 georeferenced individuals from 20 plant and animal species, we show that genome-wide diversity follows a mutations-area relationship power law with geographic area, which can predict genetic diversity loss from local population extinctions. We estimate that more than 10% of genetic diversity may already be lost for many threatened and nonthreatened species, surpassing the United Nations' post-2020 targets for genetic preservation.

Digitized collections elucidate invasion history and patterns of awn polymorphism in Microstegium vimineum

PREMISE

Digitized collections can help illuminate the mechanisms behind the establishment and spread of invasive plants. These databases provide a record of traits in space and time that allows for investigation of abiotic and biotic factors that influence invasive species.

METHODS

Over 1100 digitized herbarium records were examined to investigate the invasion history and trait variation of Microstegium vimineum. Presence-absence of awns was investigated to quantify geographic patterns of this polymorphic trait, which serves several functions in grasses, including diaspore burial and dispersal to germination sites. Floret traits were further quantified, and genomic analyses of contemporary samples were conducted to investigate the history of M. vimineum's introduction and spread into North America.

RESULTS

Herbarium records revealed similar patterns of awn polymorphism in native and invaded ranges of M. vimineum, with awned forms predominating at higher latitudes and awnless forms at lower latitudes. Herbarium records and genomic data suggested initial introduction and spread of the awnless form in the southeastern United States, followed by a putative secondary invasion and spread of the awned form from eastern Pennsylvania. Awned forms have longer florets, and floret size varies significantly with latitude. There is evidence of a transition zone with short-awned specimens at mid-latitudes. Genomic analyses revealed two distinct clusters corresponding to awnless and awned forms, with evidence of admixture.

CONCLUSIONS

Our results demonstrate the power of herbarium data to elucidate the invasion history of a problematic weed in North America and, together with genomic data, reveal a possible key trait in introduction success: presence or absence of an awn.

Traces of Genetic but Not Epigenetic Adaptation in the Invasive Goldenrod Solidago canadensis Despite the Absence of Population Structure

Biological invasions may result from multiple introductions, which might compensate for reduced gene pools caused by bottleneck events, but could also dilute adaptive processes. A previous common-garden experiment showed heritable latitudinal clines in fitness-related traits in the invasive goldenrod Solidago canadensis in Central Europe. These latitudinal clines remained stable even in plants chemically treated with zebularine to reduce epigenetic variation. However, despite the heritability of traits investigated, genetic isolation-by-distance was non-significant. Utilizing the same specimens, we applied a molecular analysis of (epi)genetic differentiation with standard and methylation-sensitive (MSAP) AFLPs. We tested whether this variation was spatially structured among populations and whether zebularine had altered epigenetic variation. Additionally, we used genome scans to mine for putative outlier loci susceptible to selection processes in the invaded range. Despite the absence of isolation-by-distance, we found spatial genetic neighborhoods among populations and two AFLP clusters differentiating northern and southern Solidago populations. Genetic and epigenetic diversity were significantly correlated, but not linked to phenotypic variation. Hence, no spatial epigenetic patterns were detected along the latitudinal gradient sampled. Applying genome-scan approaches (BAYESCAN, BAYESCENV, RDA, and LFMM), we found 51 genetic and epigenetic loci putatively responding to selection. One of these genetic loci was significantly more frequent in populations at the northern range. Also, one epigenetic locus was more frequent in populations in the southern range, but this pattern was lost under zebularine treatment. Our results point to some genetic, but not epigenetic adaptation processes along a large-scale latitudinal gradient of S. canadensis in its invasive range.

Population genomic consequences of life-history and mating system adaptation to a geothermal soil mosaic in yellow monkeyflowers

Local selection can promote phenotypic divergence despite gene flow across habitat mosaics, but adaptation itself may generate substantial barriers to genetic exchange. In plants, life-history, phenology, and mating system divergence have been proposed to promote genetic differentiation in sympatry. In this study, we investigate phenotypic and genetic variation in Mimulus guttatus (yellow monkeyflowers) across a geothermal soil mosaic in Yellowstone National Park (YNP). Plants from thermal annual and nonthermal perennial habitats were heritably differentiated for life-history and mating system traits, consistent with local adaptation to the ephemeral thermal-soil growing season. However, genome-wide genetic variation primarily clustered plants by geographic region, with little variation sorting by habitat. The one exception was an extreme thermal population also isolated by a 200 m geographical gap of no intermediate habitat. Individual inbreeding coefficients (F_IS ) were higher (and predicted by trait variation) in annual plants and annual pairs showed greater isolation by distance at local (<1 km) scales. Finally, YNP adaptation does not reuse a widespread inversion that underlies M. guttatus life-history ecotypes range-wide, suggesting a novel genetic mechanism. Overall, this work suggests that life-history and mating system adaptation strong enough to shape individual mating patterns does not necessarily generate incipient speciation without geographical barriers.

The evolution of thermal performance in native and invasive populations of Mimulus guttatus

The rise of globalization has spread organisms beyond their natural range, allowing further opportunity for species to adapt to novel environments and potentially become invaders. Yet, the role of thermal niche evolution in promoting the success of invasive species remains poorly understood. Here, we use thermal performance curves (TPCs) to test hypotheses about thermal adaptation during the invasion process. First, we tested the hypothesis that if species largely conserve their thermal niche in the introduced range, invasive populations may not evolve distinct TPCs relative to native populations, against the alternative hypothesis that thermal niche and therefore TPC evolution has occurred in the invasive range. Second, we tested the hypothesis that clines of TPC parameters are shallower or absent in the invasive range, against the alternative hypothesis that with sufficient time, standing genetic variation, and temperature‐mediated selection, invasive populations would re‐establish clines found in the native range in response to temperature gradients. To test these hypotheses, we built TPCs for 18 native (United States) and 13 invasive (United Kingdom) populations of the yellow monkeyflower, Mimulus guttatus. We grew clones of multiple genotypes per population at six temperature regimes in growth chambers. We found that invasive populations have not evolved different thermal optima or performance breadths, providing evidence for evolutionary stasis of thermal performance between the native and invasive ranges after over 200 years post introduction. Thermal optimum increased with mean annual temperature in the native range, indicating some adaptive differentiation among native populations that was absent in the invasive range. Further, native and invasive populations did not exhibit adaptive clines in thermal performance breadth with latitude or temperature seasonality. These findings suggest that TPCs remained unaltered post invasion, and that invasion may proceed via broad thermal tolerance and establishment in already climatically suitable areas rather than rapid evolution upon introduction.

Insights into invasive species from whole-genome resequencing

Studies of invasive species can simultaneously inform management strategies and quantify rapid evolution in the wild. The role of genomics in invasion science is increasingly recognised, and the growing availability of reference genomes for invasive species is paving the way for whole-genome resequencing studies in a wide range of systems. Here, we survey the literature to assess the application of whole-genome resequencing data in invasion biology. For some applications, such as the reconstruction of invasion routes in time and space, sequencing the whole genome of many individuals can increase the accuracy of existing methods. In other cases, population genomic approaches such as haplotype analysis can permit entirely new questions to be addressed and new technologies applied. To date whole-genome resequencing has only been used in a handful of invasive systems, but these studies have confirmed the importance of processes such as balancing selection and hybridization in allowing invasive species to reuse existing adaptations and rapidly overcome the challenges of a foreign ecosystem. The use of genomic data does not constitute a paradigm shift per se, but by leveraging new theory, tools, and technologies, population genomics can provide unprecedented insight into basic and applied aspects of invasion science.