BUSCADOR RECOLECTA

Figures S1–S2: 1H-NMR and 13C-NMR of compound 1a; Figure S3: 13C-NMR of compound 2; Figures S4–S5: 13C-NMR of compounds 4–5; Figures S6–S7: 1H-NMR an 13C-NMR of compound 7a; Figures S8–S10: 13C-NMR of compounds 8–10; Figures S11–S20: 1H-NMR and 13C-NMR of compounds 11–14 and 15–16; Figure S21: 1H-NMR of compound 17; Figures S22–S27: 1H-NMR and 13C-NMR of compounds 18,18a and 19., Peer reviewed

This PDF file includes: Figs. S1 to S13., Peer reviewed

Table S1: Test statistic and p-values of Kruskall–Wallis tests comparing differences in oviposition on four different oviposition substrates (cabbage, rape, pea, and aluminium foil) for the strains DBM-W, DBM-NOQA, DBM-C, and DBM-P of P. xylostella. The significance values have been adjusted by the Bonferroni correction for multiple tests. Significant p-values (p ≤ 0.05) are shown in bold type; Table S2: Comparison of the preference between abaxial and adaxial leaf surfaces among the three P. xylostella strains reared on cabbage (DBM-C), artificial diet (DBM-G88), and pea (DBM-P). Oviposition preference data were analyzed using a one-tailed, two-sample test of proportions (p ≤ 0.05) comparing the percentages of the total number of eggs laid on the abaxial side of leaves (n = 3–96, except in the case of C. papaya for DBM-C, where n = 2).; and Table S3: Total glucosinolate content (TOT) and content of aliphatic glucosinolates with sulfur-containing side chains (AS), other aliphatic glucosinolates (AO), benzenic glucosinolates (BEN), and indolic glucosinolates (IN) for each of the plant types tested (A). Glucosinolate richness (S), Shannon’s diversity index for the four glucosinolate classes (HA), Shannon’s diversity index for the relative concentrations of all individual glucosinolates (HB), and chemical complexity index for glucosinolates (CCI) for each of the plant types tested (B). Values based on means across replicates taken from Badenes-Pérez et al. 2020 [10]., Peer reviewed

Table S1: Pairwise comparisons in OPI between plant species after conducting Kruskal–Wallis tests; Table S2: Comparison between P. rapae and P. xylostella for the percentage of eggs laid on plant species tested compared to A. thaliana (n = 3). Significant differences are shown in bold type; Table S3: Significance of correlations between oviposition preference index in two-choice tests (OPI), total oviposition in no-choice tests (TO), and larval survival (LS) for P. rapae and P. xylostella in the plants tested; Table S4: Mean ± SE glucosinolate content (µmol g-1 plant dry weight) in the plants tested [37]; Table S5: Total glucosinolate content (TOT) and content of aliphatic glucosinolates with sulfur-containing side chains (AS), other aliphatic glucosinolates (AO), benzenic glucosinolates (BEN), and indolic glucosinolates (IN) for each of the plant types tested (A). Glucosinolate richness (S), Shannon’s diversity index for the four glucosinolate classes (HA), Shannon’s diversity index for the relative concentrations of all individual glucosinolates (HB), and chemical complexity index for glucosinolates (CCI) for each of the plant types tested (B) [37]; Table S6: Correlations between oviposition preference index (OPI), total oviposition (TO), and larval survival (LS) and glucosinolate richness (S), Shannon’s diversity index for the four glucosinolate classes (HA), Shannon’s diversity index for the relative concentrations of all individual glucosinolates (HB), glucosinolate complexity index (GCI), total glucosinolate content (TOT), aliphatic glucosinolates with sulfur-containing side chains (AS), other aliphatic glucosinolates (AO), benzenic (BEN), and indolic glucosinolates (IN) as shown by CATPCA analysis; Table S7: Pairwise comparisons in total oviposition in no-choice tests (TO) between plant species after conducting Kruskal–Wallis tests; Table S8: Comparison between P. rapae and P. xylostella for the percentage of eggs laid on the abaxial side of the leaves in the plant species tested (n = 3–96, except in the case of A. argenteum, C. bursa-pastoris, and T. majus for P. rapae and in the case of C. papaya and M. oleifera for P. xylostella, in which n = 2); Table S9: Pairwise comparisons in total oviposition in no-choice tests (TO) between plant species after conducting Kruskal–Wallis tests; Table S8: Pairwise comparisons in larval survival (LS) between plant species after conducting Kruskal–Wallis tests; Table S10: Origin of the seeds of the plant species tested; Figure S1: CATPCA plots showing the relationship between oviposition preference index (OPI) (A), total oviposition (TO) (B), and larval survival (LS) (C) and total glucosinolate content (TOT) and content of aliphatic glucosinolates with sulfur-containing side chains (AS), other aliphatic glucosinolates (AO), benzenic glucosinolates (BEN), and indolic glucosinolates (IN) in the plant species tested. Reference [37] is cited in the Supplementary Materials., Peer reviewed

Table S1. Soil samples collected per plant species, forest (Jaén JA; Segura SE), dominant mycorrhizal type (Myc) of plants (arbuscular mycorrhizal AM; ectomycorrhizal ECM) and season (Autumn; Spring). Enzymatic activities (pmol mg-1 h-1 SOM-1) related to carbon (C), nitrogen (N) and phosphorous (P) cycles in the soils collected under the different plant species are shown. The enzymes corresponding to each cycle are: Labile C = β-glucosidase + cellobiohydrolase; Non-labile C = β-xylosidase + β-glucuronidase + laccase; N = chitinase + leucine-aminopeptidase; P = acid + alkaline phosphatases. Table S2. Linear Mixed-Effect Models testing the effect of season, forest and plant species on soil properties and C and nutrient cycling. Season and forest were fixed factors, while the site nested in forest and the plant species nested in site and forest were random factors. Significant interactions among factors were not detected. Labile C = β-glucosidase + cellobiohydrolase; Non-labile C = β-glucuronidase + β-xylosidase + laccase; N cycling = chitinase + leucine-aminopeptidase; P cycling = acid + alkaline phosphatase. Soil properties: pH, SOM = soil organic matter, soil moisture. χ12 values for fixed and random factors and Bonferroni corrected p-values ‘***’ p < 0.001, ‘**’ p < 0.01, ‘*’ p < 0.05, ‘.’ p < 0.10, are given; significant effects are noted in bold. The coefficient of determination (pseudo-R2, i.e., variance explained) is shown for both pools of fixed and random factors. Table S3. Plant phylogenetic composition and morpho-functional traits predictors of carbon and nutrient cycling in soils, analysed by Redundancy Analysis. A Redundancy Analysis was carried out per forest, a) Jaén and b) Segura, to determine the non-redundant phylogenetic PCoA axes and morpho-functional traits that explained responses of carbon and nutrient cycling. ANOVA results (F statistic and p-value) were extracted from each Redundancy Analysis and those significant (p < 0.05) compiled in this table (see also Figure 4). Fig. S1. Relationship between N cycling-related enzymatic activities with C cycling-related ones. Pearson correlations were carried out to relate N-related enzymes (chitinase + leucine amino-peptidase) with hydrolytic (β-glucosidase + cellobiohydrolase + β-glucuronidase + β-xylosidase) and oxidative (laccase) C-related enzymes. Each plot shows the correlation coefficient (r) with its p-value. Fig. S2. Local Moran’s Index (LMI) calculated for carbon and nutrient cycling and cycling ratios and mapped onto the plant phylogeny per forest, a) Jaén and b) Segura. Positive/negative LMI values indicate that closely relatives tend to show similar/dissimilar values for a given variable. Values marked in red indicate significant and marginal autocorrelations, i.e., p < 0.05 and 0.05 > p < 0.10, respectively. Hotspots of positive (i.e., phylogenetic conservatism) and negative (i.e., phylogenetic divergence) autocorrelation are differentiated with triangles and all are represented by LMI with p < 0.05. See correspondence between acronyms and plant names in Table S1. Labile C = β-glucosidase and cellobiohydrolase; Non-labile C = β-xylosidase, β-glucuronidase and laccase; N = chitinase and leucine-aminopeptidase; P = acid and alkaline phosphatases. Figure S3. Phylogenetic composition of plant communities. Principal Coordinates Analysis (PCoA) was run over the phylogenetic distance matrix of plant species studied per forest, Jaén and Segura. Eigenvector values (PCoA-n) that significantly explained responses of carbon and nutrient cycling were used in RDA analyses. Scores of PCoA axis were plotted by plant species at each forest, a) Jaén and b) Segura, and may locate taxa across the phylogeny on the PCoA axes. The total percentage of variance explained by significant PCoA axes was 3.4 % for Jaén and 0.2 % for Segura., Peer reviewed

Table S1. Location, altitude and dominant vegetation at each study site. Table S2. Samples distribution by species, season, region and site. Figure S3. Sequencing rarefaction curves and sampling completeness. Methods S4. Root processing, DNA extraction, sequencing and bioinformatic analyses. Table S5. A recruitment network represented by its recruitment matrix. Table S6. Description of main network metrics and indices used throughout the work. Methods S7. General metrics describing the plant-AMF networks studied. Methods S8. Null models for modularity. Table S9. Maximum likelihood tree of the sequences defining the AMF virtual taxa sampled in this study. Table S10. PERMANOVA exploring the effect of region, season, site and plant species on AMF community composition. Table S11. Bayesian glmm testing the influence of the degree of canopy plant species on the dissimilarity in AMF communities between canopy and recruit species., Peer reviewed

Derivation of the dispersion relation in high-q approximation considering the TM modes only; description of the method of the calculation of the isofrequency curves by minimizing the losses; derivation of the expressions for the directions of the group velocity and the imaginary part of the wavevector; comparison of the IFC obtained with the different approaches; analytical expression for the field distribution; conditions for the surface or volume type of the mode; and derivation of the dispersion relation of the TE modes in the ultrathin-slab limit., Peer reviewed

[EN] WETYDAS v.2.0.0. (Weather Types Data set of Spain) contains 143 TXT archives localized by their coordinates in 20-century reanalyses grid (from 36º-44 lat North and 350º-5º longitude) covering 1836 to 2015; information include Year, Month, Day, Direction (in deg), WT.dir (8 classes, only dependent on direction), Hyb (separate between directional type, Pure A or C types, or Hybrid type), WT (final classification with 26 classes), WT.num (same as WT column but the names of the 26 WT classes are converted to numbers from 1 to 26; classes are ordered from NE = 1), WT10 (idem as WT but reclassified to 10 classes only: 8 directional types + C and A), WT10.num (idem as WT10 but converted in numbers from 1 to 10), SF (geostrophic southerly flow index), WF (geostrophic westerly flow index), F (total flow index: square root of SF^2+WF2^2), ZS (southerly shear vorticity), ZW (easterly shear vorticity) and Z (shear vorticity: ZS + ZW). For details of the computation see Trigo and DaCamara (2000) Int. Jr. Clim. https://doi.org/10.1002/1097-0088(20001115)20:13<1559:AID-JOC555>3.0.CO;2-5; [ES] WETYDAS v.2.0.0. (Weather Types Data set of Spain), contiene 143 archivos TXT con la clasificación de tipos de tiempo en los puntos de la malla del reanálisis 20-century (36º-44 lat North y 350º-5º longitude) period 1836- 2015; la información se incluye en columnas con el año, mes y día, la dirección (en grados), WT.dir (8 clasificación en ocho clases solamente acorde la direción del flujo), Hyb (identificación entre tipos direccionales, A Puro, C Puro o tipo Híbrido), WT (clasificación final en 26 tipos), WT.num (igual que WT pero numeradas las clases con NE = 1), WT10 (idem a WT con clases reclasificadas a 10 tipos, 8 direcionlaes + C y A), WT10.num (idem a WT10 convirtiendo a número de 1 a 10), SF (índice del flujo geostrófico sur), WF (índice del flujo geostrófico oeste), F (iínidce flujo total: raíz cuadrada de SF^2+WF2^2), ZS (shear vorticity meridional), ZW (shear vorticity del este) y Z (shear vorticity: ZS + ZW). Detalles de cálculo según Trigo and DaCamara (2000) Int. Jr. Clim. https://doi.org/10.1002/1097-0088(20001115)20:13<1559:AID-JOC555>3.0.CO;2-5. Validación del método aplicado a la malla del reanálisis en Fernández-Granja, J.A., Brands, S., Bedia, J. et al. Exploring the limits of the Jenkinson–Collison weather types classification scheme: a global assessment based on various reanalyses. Clim Dyn 61, 1829–1845 (2023). https://doi.org/10.1007/s00382-022-06658-7, [ES] Hemos calculado la nueva base de datos de alta resolución espacial empleando la clasificación de Jenkinson & Collison clasificación de tipos de tiempo a la Península Ibérica en la malla de presiones diaria del reanálisis del siglo (Enero-1836 a Diciembre 2015). Dada su resolución, 1ºx1º lat/long, WETYDAS v.2.0.0. contiene 143 archivos acorde las coordenadas de la ventana analizada (36º-44 lat Norte/ 350º-5º longitud). Los tipos de tiempo incluyen 8 direccionales puros, el tipo Anticiclónico y Ciclónico puros más 8 Ciclónicos híbridos y 8 Anticiclónicos Híbridos (total 26 tipos). Se recalculan reclasificaciones para evitar los tipos híbridos. Los casos no determinados fueron diseminados entre la clasificación., [EN] It has been applied the Jenkinson & Collison classification of Weather Types to Iberian Peninsula and Balearic Island by using the daily NOAA/CIRES/DOE 20th Century Reanalysis (V3) dataset (January-1836/December-2015). WETYDAS v.2.0.0. grid resolution/1.0ºx1.0 lat/long produces 143 series. Weather Types classification includes 8 directional pure, Anticyclonic and Cyclonic pure types, and combination of previous ones in the hybrid types. Non determines cases were spread homogeneously., No

21679463.zip contains: Cera et al. PED.R; Supplementary files_25112022.docx; data_Cera et al. Plant Ecology and Diversity.xlsx, Plants growing on extreme soils have mainly been described in relation to their adaptations to edaphic conditions, although herbivores may also be an important factor in these ecosystems. Gypsum soils occur in drylands often where livestock practices occur. However, it is unknown whether plant traits related to gypsum soil constraints are associated with resistance to herbivory. In order to assess whether gypsum specialist species might be favoured at higher grazing levels and to detect the traits involved, we evaluated the responses of gypsum specialists vs. generalists to three intensities of livestock pressure. We analysed the relative cover shifts of species along a livestock gradient, and variation in canopy height, canopy area, leaf carbon (C), nitrogen (N), and sulphur (S), specific leaf area (SLA) and leaf dry matter content (LDMC). We found that gypsum-specialists responded by increasing or maintaining their cover at medium and high grazing pressure, whereas most generalists responded by decreasing it. Gypsum-specialists showed higher leaf S than generalists, regardless of grazing intensity. All species showed similar patterns for traits linked to loss of above-ground biomass when grazing increased. Plant affinity to gypsum soils mediates vulnerability to grazing with foliar S possibly being a defence trait., This work was supported by Gobierno de España [MICINN, CGL2015-71360-P, CGL2016-80783-R and PID2019-111159GB-C31]; by European Union’s Horizon 2020 [H2020-MSCA-RISE-777803]; and by Consejo Superior de Investigaciones Científicas [COOPB20231]. AC and SP were funded by a FPI fellowship [MICINN, BES-2016-076455] and a Ramón y Cajal Fellowship [MICINN, RYC-2013-14164], respectively., Peer reviewed

[ES] MOPREDAS_Annual_Books_1916-1950 se compone de los datos digitalizados de los Libros Resúmenes Anuales editados por los sucesivos Servicios Meteorológicos de España en el periodo 1916-1950; [EN] MOPREDAS_Annual_Books_1916-1950 is a compilation of monthly precipitation data digitalized from Annual Books 1916-1950., [ES] Mopredas_Annual_Books_1916-1950 es el resultado de un proceso de digitalización de los datos publicados en los Libros Resúmenes Anuales editados por los servicios meteorológicos de España, periodo 1916-1950. Combinados con los datos incluidos en el BNDC han sido empleados para el desarrollo de las mallas de alta resolución denominadas MOPREDAS_1916-2015 y subsiguientes versiones, [EN from original) ] Mopredas_Annual_Books_1916-1950 is the result of a long digitalization procedure of Annual Books edited by Spanish Meteorological Services; the compilation includes hundreds of data not at present included in the Banco Nacional de Datos del Clima and increase the spatial and temporal density of precipitation records in pre-1950 period. Combined with BNDC after identification of digitalized data with BNDC observatories, it has been uses to develop the new high spatial density grid MOPREDAS_1916-2015 and versions, CGL2014-52135-C3-3-R Desarrollo de índices de sequía sectoriales: mejora de la monitorización y alerta temprana de las sequías en España DESEMON. Gobierno de España y UE, CGL2011-27574-C02-01 Impactos Hidrológicos del Calentamiento Global en España HIDROCAES Gobierno de España y UE, CGL2008-05112-C02-01/CLI. Proyecto Cambio climático: base de datos de precipitaciones, análisis de tendencias e impactos en los sistemas naturales Gobierno de España y UE, Three .csv files, No, Manual