BUSCADOR RECOLECTA

Feral goat eating Chamaerops humilis inflorescences recorded by camera-trapping (camera-trap LTL ACORN 5310A, detection range = 18 m) on 04/27/2020 in Mallorca (Balearic Islands, Spain). Note that it is a composition of video fragments at x1.5 speed., Peer reviewed

2 files. -- File 1 includes: Figure S1. Palm study populations in Mallorca (Balearic Islands, Spain): EB, Ermita de Betlem; PS, Puig de Santuïri; PF, Platja de Formentor; CV, Campament de la Victòria. -- Figure S2. The fate of female inflorescences was sampled to estimate the effect of (a) moth herbivory and (b) goat herbivory on fruit initiation and fruit development for two sampling years (2019 and 2020). -- Figure S3. Isolated effects of goat herbivory on palm reproduction. -- Figure S4. Interaction effects of moth herbivory and goat herbivory on palm reproduction. -- Table S1. Other inflorescence visitors observed in the dwarf palms during the visual censuses carried out in Mallorca during the springs of 2019 and 2020. -- File 2 includes: Figure S1. Figure 3 from the main document including raw data points. -- Figure S2. Figure 4a from the main document including raw data points for (a) 2019 and (b) 2020. -- Figure S3. Figure 4b from the main document including raw data points., Appendix_S1_Non-additive.pdf, Appendix_S2_Non-additive.pdf, Peer reviewed

17 pages. -- Figure S1: Autocorrelation of original signals and residual components of ENSO NOAA oceanic indexes (a)-(c)-(e)-(g), of TMI study cases shown in Fig. 4 (b)-(d)-(f)-(h) and of ENSO indexes derived by Takahashi et al. (2011) (i)-(j), and Sullivan et al. (2016) (k)-(l)-(m). -- Figure S2: (top panel) Hovmöller diagram depicting the zonal propagation of Kelvin waves along the equatorial Pacific Ocean between years 1993 and 2017 (both included). (bottom panel) Time series of the zonally-averaged Kelvin sea level anomalies over the equatorial Pacific Ocean (between 145°E and the Pacific coast of Colombia between -4°S and 4°N of latitude). -- Figure S3: (top panel) Hovmöller diagram illustrating the zonal propagation of wind-forcing associated to Kelvin waves along equatorial Pacific Ocean between years 1993 and 2017 (both included). (bottom panel) Time series of zonally-averaged Kelvin sea level anomalies over the equatorial Pacific Ocean (between 145°E and the Pacific coast of Colombia between -4°S and 4°N of latitude). -- Figure S4: (a)-(c) Maps of cross-correlation between the TMI and the time series of CHIRPS rain associated with every grid cell between years 1983 and 2017. -- Figure S5: (a)-(c) Maps of cross-correlation between the ENSO 1+2 and the time series of CHIRPS rain associated with every grid cell between years 1983 and 2017. -- Figure S6: (a)-(c) Maps of cross-correlation between the ENSO 3 and the time series of CHIRPS rain associated with every grid cell between years 1983 and 2017. -- Figure S7: (a)-(c) Maps of cross-correlation between the ENSO 3.4 and the time series of CHIRPS rain associated with every grid cell between years 1983 and 2017. -- Figure S8: (a)-(c) Maps of cross-correlation between the ENSO 4 and the time series of CHIRPS rain associated with every grid cell between years 1983 and 2017. -- Figure S9: Map showing the mean SST field (shading), and mean surface wind (black arrows) during years 1988–2017. Data is described in section 3. For the sake of clarity, only 1 of every 7 arrows are shown. -- Figure S10: Mean (panels a-c-e-g) and standard deviation (panels b-d-f-h) composites of rainfall anomalies for different ENSO types. -- Figure S11: (a) Time series of 20◦C isotherm depth (z20◦C ) at 95◦W from Tropical Atmosphere Ocean project (TAO) data; (b) Time series of z20◦C at 110◦W from TAO; (c) Time series of the difference between z20◦C at 110◦W and at 95◦W; (d) Scatter diagram between z20◦C at 95◦W and at 110◦W, where the red line represents a linear fit between both depths; (e) original time series of z20◦C at 110◦W (black points) and the reconstructed time series OF z20◦C at 110◦W (sky blue solid line) after filling gaps with 95◦W data and performing a quadratic interpolation. --Figure S12: Time series of TMI4∗ (a, 20◦C isotherm depth instead of sea level), and TMI5 (b, 20◦C isotherm depth is added to the original TMI4. -- Figure S13: (a)-(b) Maps of cross-correlation between TMI4∗ (which includes the 20◦C isotherm depth instead of sea level) and the time series of CHIRPS rain associated with every grid cell between years 1983 and 2017. -- Figure S14: (a)-(b) Maps of cross-correlation between TMI5 and the time series of CHIRPS rain associated with every grid cell between years 1983 and 2017. -- Table S1: Annual mean (x) and standard deviation (σ) of rainfall from CHIRPS data during the period 1983–2017 (both included). -- Table S2: Number of CHIRPS grid cells (Ng) with time lags within the interval [-6, 0) (time lag = 0 not included) for TMI, TMI4 and ENSO oceanic indices: ENSO 1+2, ENSO 3, ENSO 3.4, ENSO 4. -- Table S3: Number of CHIRPS grid cells (Ng) with time lags within the interval [-6, 0) (time lag = 0 not included) for TMI4 and Takahashi et al. (2011) (C, P) and Sullivan et al. (2016). -- Table S4: Cross-correlation of original TMI4 and the new TMI4 (renamed as TMI4 ∗ ) and tMI5, against selected ENSO indices: ENSO 1+2, ENSO 3, ENSO 3.4 and ENSO 4, and other indices more adapted to regionally identify Central Pacific and Eastern Pacific ENSO events such as those of Takahashi et al. (2011) (C, E) and Sullivan et al. (2016) (CP, EP, Mixed), which are described in Section 3.4. -- Table S5: Same as Table S4 but for the residual components (denoted by the subscript r). -- Table S6: Number of CHIRPS grid cells (Ng) with time lags within the interval [-6, 0) (time lag = 0 not included) for the original TMI4 (without z20◦C ), TMI4∗ (which includes z20◦C instead of sea level) and TMI5., Peer reviewed

Figure S1. Weekly average of maximum and minimum daily temperatures (ºC) and weekly-accumulated rainfall and ETo (mm) during the two agricultural seasons of this study (2016/17-2017/18) at the experimental site. Figure S2. A) Scheme of the experimental design, showing the distribution of tillage system (CTR and ZTR) and irrigation strategy (FI and RDI) treatments in each block. B) Arrangement of the three consecutive PVC collars in an elementary plot. Each group of three collars [+F, L, -F] in consecutive positions constitutes a measuring site, or unit of measure, for the two main factors (tillage system and irrigation strategy). Table S1. Number of irrigations and fortnight irrigation amount for FI and RDI treatments during the 2016 and 2017 irrigation seasons. Table S2. Daily mean, and minimum and maximum diurnal soil C-CO2 effluxes (mg m-2 h-1) measured during the three diurnal cycles in the CTR and ZTR tillage systems, and relative difference in relation to the 24-h mean for the minimum and maximum diurnal values (%; negative value = underestimation; positive values = overestimation). CV = coefficient of variation (%)., Peer reviewed

Supporting Information Table S5. RT-qPCR primers for the thirteen selected genes and the two reference genes used in this study., The impact of the versatile biocontrol and plant-growth-promoting rhizobacteria Pseudomonas simiae PICF7 on the banana holobiont under controlled conditions was investigated. We examine the fate of this biological control agent (BCA) upon introduction in the soil, the effect on the banana root microbiota, and the influence on specific host genetic defense responses. While the presence of strain PICF7 significantly altered neither the composition nor the structure of the root microbiota, a significant shift in microbial community interactions through co-occurrence network analysis was observed. Despite the fact that PICF7 did not constitute a keystone, the topology of this network was significantly modified—the BCA being identified as a constituent of one of the main network modules in bacterized plants. Gene expression analysis showed the early suppression of several systemic acquired resistance and induced systemic resistance (ISR) markers. This outcome occurred at the time in which the highest relative abundance of PICF7 was detected. The absence of major and permanent changes on the banana holobiont upon PICF7 introduction poses advantages regarding the use of this beneficial rhizobacteria under field conditions. Indeed a BCA able to control the target pathogen while altering as little as possible the natural host-associated microbiome should be a requisite when developing effective bio-inoculants., Peer reviewed

Supplementary Table S4. Pseudomonas simiae strain PICF7 (GenBank: CP0052975.1) 16S rRNA gene sequence comparison with Amplicon Sequence Variants belonging to Pseudomonas species found in Grand Naine Control and Grand Naine PICF7 samples., The impact of the versatile biocontrol and plant-growth-promoting rhizobacteria Pseudomonas simiae PICF7 on the banana holobiont under controlled conditions was investigated. We examine the fate of this biological control agent (BCA) upon introduction in the soil, the effect on the banana root microbiota, and the influence on specific host genetic defense responses. While the presence of strain PICF7 significantly altered neither the composition nor the structure of the root microbiota, a significant shift in microbial community interactions through co-occurrence network analysis was observed. Despite the fact that PICF7 did not constitute a keystone, the topology of this network was significantly modified—the BCA being identified as a constituent of one of the main network modules in bacterized plants. Gene expression analysis showed the early suppression of several systemic acquired resistance and induced systemic resistance (ISR) markers. This outcome occurred at the time in which the highest relative abundance of PICF7 was detected. The absence of major and permanent changes on the banana holobiont upon PICF7 introduction poses advantages regarding the use of this beneficial rhizobacteria under field conditions. Indeed a BCA able to control the target pathogen while altering as little as possible the natural host-associated microbiome should be a requisite when developing effective bio-inoculants., Peer reviewed

Supplementary Figure 1. Principal coordinates analyses of the banana root microbial communities. PCoA of bacterial (A) and fungal (B) communities by treatment on Bray-Curtis dissimilarities. GNC: control plant, GNF7: plant inoculated with Pseudomonas simiae PICF7, DAI: days after inoculation. Supplementary Figure 2. The rarefaction curve. The abscissa is the number of sampling sequences and the ordinate is the number of Amplicon Sequence Variants (ASV) in the sampling sequences. Different colors indicate different groups. Black color indicates control and inoculated plants (GNC and GNF7) and orange color Grand Nain in vitro plants (GNI). The ASV 0 to 100 range has been magnified to show differences. Supplementary Figure 3. Box plot of α-diversity indices: InvSimpson, observed Amplicon Sequence Variants (ASV) and Shannon for the fungal community comparing Grand Nain in vitro (GNI), control (GNC) and Pseudomanas simiae PICF7-inoculated (GNF7) plants. Supplementary Figure 4. Principal coordinates analyses of the banana root microbial communities. PCoA of bacterial by treatment on Weighted Unifrac distances (A) and fungal communities by treatment on Bray-Curtis dissimilarities (B). GNI: Grand Nain in vitro, GNC: control plant, GNF7: plant inoculated with Pseudomonas simiae PICF7. Supplementary Figure 5. Relative expression levels of defence-related genes in roots at five time points upon inoculation with Pseudomonas simiae PICF7. (A) Non-expresser of PR genes 1 (NPR1) and Pathogenesis-related 1 (PR1); (B) Isochorismate synthase 1 (ICS1) and Phenylalanine ammonia-lyase (PAL); (C) Allene oxide synthase (AOS) and Jasmonic acid carboxyl methyltransferase (JMT); (D) Anthranilate synthase alpha subunit 1 (ASA1) and Polyphenol oxidase (PPO); (E) Enhanced disease susceptibility 1 (EDS1) and Phytoalexin deficient 4 (PAD4); (F) Manganese superoxide dismutase 1 (MSD1) and Coronatine insensitive 1 (COI1); (G) Ascorbate peroxidase (APX1). CNRQ: Calibrated Normalized Relative Quantity. Musa acuminata genes EF-1 and L2 were used as internal controls to normalize the expression data. Letters in black and bold type indicate significant differences based on comparison over treatments and time points (LSD; α = 0.05). Letters in grey indicate non-significant differences. Error bars: standard error of the mean (SEM). DAI: days after inoculation. GNC: non-inoculated plants; GNF7: Pseudomonas simiae PICF7-inoculated plants. N control/inoculated=4/4. Supplementary Figure 6. Relative expression levels of defence-related genes in leaves at five time points upon inoculation with Pseudomonas simiae PICF7. (A) Non-expresser of PR genes 1 (NPR1) and Isochorismate synthase 1 (ICS1), (B) Enhanced disease susceptibility 1 (EDS1) and Phytoalexin deficient (PAD4), (C) Coronatine insensitive 1 (COI1) and Anthranilate synthase alpha subunit 1 (ASA1), (D) Phenylalanine ammonia-lyase (PAL) and Polyphenol oxidase (PPO), (E) Ascorbate peroxidase (APX1) and Manganese superoxide dismutase 1 (MSD1). CNRQ: Calibrated Normalized Relative Quantity. Musa acuminata genes EF-1 and L2 were used as internal controls to normalize the expression data. Letters in black and bold type indicate significant differences based on comparison over treatments and time points (LSD; α = 0.05). Letters in grey indicate non-significant differences. Error bars: standard error of the mean (SEM). DAI: days after inoculation. GNC: non-inoculated plants; GNF7: Pseudomonas simiae PICF7-inoculated plants. N control/inoculated=4/4. Supplementary Table S1. Taxonomical profile of Amplicon Sequence Variants (ASV) present in at least 4 of the 10 Grand Nain in vitro samples analysed. Supplementary Table S2. p-values of bacterial and fungal richness (Observed Amplicon Sequence Variants; ASV), Shannon and InvSimpsom α-diversity indices for the different comparisons analysed. Significant p-values and its correspondent significant comparisons are shown in bold type and italics. GNI: Grand Nain in vitro plants, GNC: control plants, GNF7: Pseudomonas simiae-inoculated plants. Supplementary Table S3. p-values of PERMANOVA and BETADISPER analysis of quantitative β-diversity index for the different comparisons analysed. GNI: Grand Nain in vitro plants, GNC: control plants, GNF7: Pseudomonas simiae-inoculated plants. Significant p-values and its correspondent significant comparisons are shown in bold type and italics. Supplementary Table S4. Pseudomonas simiae strain PICF7 (GenBank: CP005975.1) 16S rRNA gene sequence comparison with Amplicon Sequence Variants (ASV) belonging to Pseudomonas species found in Grand Nain Control and Grand Nain PICF7 samples. Supplementary Table S5. List of defence related genes studied by real-time PCR experiments in root and aerial tissues of banana plants (cv. Grand Nain) grown under control conditions upon Pseudomonas simiae PICF7 inoculation. For all transcripts, Relative expression analysis was repeated at least two times in independent (4 replicates each one) real-time qPCR experiments. Gene names, abbreviation, Arabidopsis thaliana orthologs, Musa acuminata paralogs, primers sequences, amplicon size, annealing temperature, PCR efficiencies, correlation coefficients (R2) and linear equations are indicated. Fw: forward, Rv: reverse., The impact of the versatile biocontrol and plant-growth-promoting rhizobacteria Pseudomonas simiae PICF7 on the banana holobiont under controlled conditions was investigated. We examine the fate of this biological control agent (BCA) upon introduction in the soil, the effect on the banana root microbiota, and the influence on specific host genetic defense responses. While the presence of strain PICF7 significantly altered neither the composition nor the structure of the root microbiota, a significant shift in microbial community interactions through co-occurrence network analysis was observed. Despite the fact that PICF7 did not constitute a keystone, the topology of this network was significantly modified—the BCA being identified as a constituent of one of the main network modules in bacterized plants. Gene expression analysis showed the early suppression of several systemic acquired resistance and induced systemic resistance (ISR) markers. This outcome occurred at the time in which the highest relative abundance of PICF7 was detected. The absence of major and permanent changes on the banana holobiont upon PICF7 introduction poses advantages regarding the use of this beneficial rhizobacteria under field conditions. Indeed a BCA able to control the target pathogen while altering as little as possible the natural host-associated microbiome should be a requisite when developing effective bio-inoculants., Peer reviewed

Fig. S1. Best topology design of the MLP based ANN algorithm. Fig. S2. impact of individual parameters on the removal of BPA from polluted water. Table S1. Enzyme kinetics constants of free laccase and Ba-alginate beads containing laccase., Peer reviewed

Document S1. Figures S1–S6 and supplemental experimental procedures. Table S1. Gene expression changes (> 1.6 and < −1.6 fold-change) in cells depleted of PIK3CB, or PIK3CA, or N/C-MYC (72 h) compared with controls cells. Document S2. Article plus supplemental information., Peer reviewed

SEM micrographs of Ag-P, Sr-P, and Zn-P (Figure S1); EDX area mapping showing the homogeneous distribution of Ag, Sr, and Zn in phosphate phases (Figure S2); powder X-ray diffraction (PXRD) patterns of Ag-P, Sr-P, and Zn-P (Figure S3); and normalized average weight and atomic percentages of phosphorus, calcium, silver, strontium, and zinc from the three highly doped phases based on EDX elemental analysis (Table S1)., Peer reviewed