Unraveling the effects of climate change on the Patagonian genus Onuris (Brassicaceae: Eudemeae) at different phylogenetic scales
Creators: name: Salariato, Diego
name: Zuloaga, Fernando O.
Publisher: name: figshare
Year: 2024
Subjects: classificationCode: 310405
subject: Evolutionary ecology
classificationCode: 310410
subject: Phylogeny and comparative analysis
classificationCode: 310406
subject: Evolutionary impacts of climate change

description: <b>Appendix A</b> Taxa and GenBank accession numbers for the ITS, <i>trnL</i>-<i>F</i>, <i>trnH</i>-<i>psbA</i>, <i>rps16</i> intron, and <i>trnG</i> intron sequences used in the phylogenetic analyses. For ingroup sequences (<i>Onuris</i>) geographic coordinates and the acronym used in the phylogenetic analysis is also supplied. New sequences generated for this study have underlined accession numbers.<b>Fig. S1.</b> Study area used in climate niche comparisons, species distribution modeling, and climate change impact assessments. Study area is defined by a minimum convex polygon enclosing all specimen records of <i>Onuris</i> with a 50-km buffer zone. Dots represent specimen records for the genus.<b>Fig. S2.</b> Axis values for <i>Onuris</i> specimens obtained from linear discriminant analyses (LDA) when climatic niche differentiation was explored using values of the five climatic components for the specimens as predictor, and the different phylogenetic groups as response (i.e., species, ESUs, and lineages). (a) classification by lineages. (b) classification by ESUs. (c) classification by species. (d) performance values of the discriminant models for each phylogenetic category, evaluated with the Accuracy and Cohen's kappa coefficients and leave-one-out cross-validation.<b>Fig. S3.</b> Results from the species distribution modeling for <i>Onuris</i> species. Predicted suitable climatic conditions (logistic output) for current projections obtained using ensemble of small models (ESM) technique and MaxEnt. OA: <i>O. alismatifolia</i>; OG: <i>O. graminifolia</i>; OH: <i>O. hatcheriana</i>; OP: <i>O. papillosa</i>; OS: <i>O. spegazziniana</i>.<b>Fig. S4.</b> Suitability change maps (future suitability –current suitability) for each species and SSP. Suitability change values for each species and SPP correspond to the mean obtained by averaging the individual values obtained with the five GCMs assessed (GFDL-ESM4, IPSL-CM6A-LR, MPI-ESM1-2-HR, MRI-ESM2-0, and UKESM1-0-LL). OA: <i>O. alismatifolia</i>; OG: <i>O. graminifolia</i>; OH: <i>O. hatcheriana</i>; OP: <i>O. papillosa</i>; OS: <i>O. spegazziniana</i>.<b>Fig. S5.</b> Barplots showing for each species the percentage of net change in cell occupancy for the future (% cell gain – % cell loss). Values are mapped on the <i>Onuris</i> MCC species tree. Bars for each species correspond to values recovered for each of the three SSPs evaluated (1-2.6, 3-7.0, 5-8.5) averaged by the five GCMs and the two binarization thresholds used. NL: northern lineage; OA: <i>O. alismatifolia</i>; OG: <i>O. graminifolia</i>; OH: <i>O. hatcheriana</i>; OP: <i>O. papillosa</i>; OS: <i>O. spegazziniana</i>; SL: southern lineage.<b>Fig. S6.</b> Histograms showing the impact of climate change at the genus level in 5 arc-min cells, and using the MPA90 threshold for binarization. (a) Cell frequency for predicted change in richness (number of species). (b) Cell frequency for predicted change in phylogenetic diversity, estimated with the mean pairwise distance metric (mpd) and the MCC species tree. **: Future <sub>richness/MPD</sub> - Current <sub>richness/MPD</sub> &lt; 0 at α=0.01 with the paired Wilcoxon Signed Rank Test.<b>Fig. S7.</b> Correlation plots for shared pattern of local extirpation (cell loss) among <i>Onuris</i> species, estimated in 5 arc-min cells, using the MPA90 threshold for binarization, and captured by unconstrained ordination using a generalized linear latent variable model. Positive correlation is depicted in blue, negative correlation in red, and for each SPP evaluated the UPGMA tree corresponding to this correlation matrix is displayed. OA: <i>O. alismatifolia</i>; OG: <i>O. graminifolia</i>; OH: <i>O. hatcheriana</i>; OP: <i>O. papillosa</i>; OS: <i>O. spegazziniana</i>.<b>Table S1</b> Specimen occurrence data for <i>Onuris</i> compiled for this study.<b>Table S2</b> Loadings on the first five components obtained from the principal component analysis using bioclimatic variables (PCA-env) extracted from the study area.<b>Table S3</b> Features of the DNA alignments used in the phylogenetic analyses.<b>Table S4</b> Performance of distribution models (mean ± sd) obtained for <i>Onuris</i> species using the ensemble of small models (ESM) technique and ten cross-validations.<b>Table S5.</b> Multivariate environmental similarity surface (MESS) analysis to assess the effect of model extrapolation on predictor values that fall outside the calibration range.<b>Table S6</b> Variable contribution for distribution models of <i>Onuris</i> species.<b>Table S7</b> Loadings on the first two components obtained from the principal component analysis using suitability change obtained for each <i>Onuris</i> species and averaged for the five GCMs evaluated.<b>ESM1.zip</b> DNA alignments used in phylogenetic analyses.<b>ESM2.txt</b> Occurrence data used in niche and SDM analyses.