BUSCADOR RECOLECTA

All images correspond to projections of confocal sections. (A-D) At early stages, overexpressed GFP-Kkv accumulates apically (white arrow in B) and in intracellular punctae (blue arrow in B), as endogenous Kkv (white and blue arrows in A), but also in the whole cell. At later stages, GFP-Kkv shows a pattern in stripes that corresponds to the taenidial folds (inset in D), in a comparable pattern to the endogenous Kkv (inset in C). Endogenous Kkv at late stages localises mainly apical and almost no intracellular punctae are detected (C). GFP-Kkv also localises mainly apical, but in addition, Kkv intracellular punctae are also detected (blue arrow in D). (E-H) In trachea, the simultaneous overexpression of reb and GFP-kkv anticipates chitin deposition (compare E and F). At later stages, this results in different morphogenetic defects like short and straight tubes and defects in branch fusion (compare H and G). (I) In salivary glands, the coexpression of reb and GFP-kkv promotes accumulation of chitin in the lumen. Scale bars: 10 μm., Peer reviewed

Under a Creative Commons license BY NC ND 4.0., Figure S1: Calibration mix chromatogram obtained for ASTM D2887 method. Figure S2: Calibration curve obtained for ASTM D2887 method. Figure S3: Chromatograms obtained for the a) TPO, the first distillation of TPO b) light fraction (LF), c) heavy fraction (HF) and the second distillation of the light fraction of TPO d) LF and e) HF. Table S1: Percentage of relative area obtained with the quantification ion (m/z) for the TPO by GC-MS according to the NIST2020 library (BTEX=benzene, toluene, ethylbenzene, xylenes; CC= cyclic compounds, AAA= alkanes, alkenes, alkynes, SB= Substituted benzenes, no BTEX, HC= heterocyclic compounds, I= indenes, PAH= polycyclic aromatic hydrocarbons, O= others). Table S2: Percentage of relative area obtained with the quantification ion (m/z) for the first distillation of the TPO, LF-1, by GC-MS according to the NIST2020 library. Table S3: Percentage of relative area obtained with the quantification ion (m/z) for the first distillation of the TPO, HF-1, by GC-MS according to the NIST2020 library. Table S4: Percentage of relative area obtained with the quantification ion (m/z) for the second distillation of the TPO, LF-2, by GC-MS according to the NIST2020 library. Table S5: Percentage of relative area obtained with the quantification ion (m/z) for the second distillation of the TPO, HF-2, by GC-MS according to the NIST2020 library., This work is part of the BLACKCYCLE project (For the circular economy of tyre domain: recycling end of life tyres into secondary raw materials or tyres and other product applications) which has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement N° 869625. The authors would also like to thank the Regional Government of Aragon (DGA) for the support provided under the research groups support programme and CSIC for the interdisciplinary thematic platform SUSPLAST., Peer reviewed

Under creative commons license CC-BY 4.0., Combined data on the migration of the studied standards (m.d. or RF) and UV spectra on silica gel. Migration of standards under different conditions and their recorded on-silica gel UV spectra. UV spectra of the separate peaks were consistent with those of the model compounds studied.-- Fig. 1S Migration distance (m.d. in mm, RF = migration distance of analyte / migration distance of the elution front) of studied standards in classic SARA. Fig. 2S Migration distance (m.d. in mm, RF = migration distance of analyte / migration distance of the elution front) of studied standards in 23-step THF‒DCM‒C7 AMD gradient. Fig. 3S On silica UV spectra of some studied standards: A) 1-octadecene, B) Tetralin, C) 9,10-Dihydrophenanthrene, D) Ethyl-naphthalene, E) Pyrene, F) Benzo (ghi) perylene., Authors thank TOTAL RM (France) for financial support for this project., Peer reviewed

S2. Redox charge density & thickness: Figure S2. Cyclic voltamogramms of the different P3HT film series with different film thicknesses (indicated by different colors) taken at a scan rate of 20 mV/s. A) P3HT-CHCl3 B) P3HT-THF C) P3HT-CHCl3(NP) D) P3HT-THF(NP). Respective figures E), F) ,G) and H) represent the linear relations found between resulting redox charge densities and measured thicknesses (colors of each film thickness according to the corresponding cyclic voltammogram). Redox charge density values in units of mCcm-2 are calculated by the integration of the anodic and cathodic CV curves for each thickness followed by dividing the obtained values by the scan rate and the probed surface area.-- Figure S3. Transmittance values at 520 nm for the P3HT film series taken in their A) oxidized transparent states and B) neutral colored states as a function of the calculated redox charge density. Symbols represent experimental data points. Lines represent the fitted curves according to exponential decay functions. Regression coefficients r2 are indicated for all the film series. The transmittance curve for P3HT-CHCl3 in the neutral state not only shows overall higher transmittance values but also its exponential fitting curve with lowest regression coefficient strongly deviates from the behavior of all the other films. The different and rather poor transmittance behavior reflects the non-continuous island-like coating of the ITO substrate obtained for this sample by the employed spray-coating process, as demonstrated by the SEM and profilometry results in the main manuscript (Figure 5D and 6D, respectively). The poor transmittance behavior of the spray-coated P3HT-CHCl3 sample in the neutral state thus accounts for the largely different contrast behavior compared to the other P3HT samples, with an apparent shift of the optimum contrast range towards rather high, i.e. out-of-range redox charge density values, as a consequence of unsatisfying fitting results.--, Figure S4. Transmittance spectra of a P3HT film in its neutral and oxidized states.-- Figure S5. A) Transmittance vs. time plot a P3HT-THF film. B) Contrast vs. pulse length extracted from the previous data and corresponding fitting function from which τ values can be obtained.-- Figure S6. A) Transmittance vs. time plot for a P3HT-THF film. B) Contrast vs. number of cycles extracted from the previous data and corresponding fitting function from which N80 values can be obtained.-- Figure S7. Characterization of spin-coated P3HT films deposited from chloroform. A) SEM image obtained at 30 kX magnification (scale bar 200 nm). B) Profilometry of a representative film with average thickness of 77 nm. C) Transmittance in transparent and colored states, together with resulting contrast, versus redox charge density. D) Switching speeds, represented by t90 values versus redox charge density. E) Cycling stability, represented by the number of cycles corresponding to a 20 % loss of the initial contrast value, i.e. N80 value. F) Images of delaminated films after cycling stability tests. Table S1. Electrochromic performance parameters for spin-coated P3HT-CHCl3 series.--, Under a Creative Commons license BY-NC-ND 4.0., S1. UV-vis absorption spectra. S2. Redox charge density & thickness. S3. Stationary transmittance at 520 nm. S4. Transmittance spectra. S5. Switching speed. S6. Stability test. S7. Spin-coated P3HT-CHCl3 film. References., S1. UV-vis absorption spectra: The UV-vis spectra (Figure S1) of P3HT-THF and CHCl3 solutions show a featureless broad π-π* transition absorption band with its maximum at 445 nm, typical for amorphous P3HT. This band is red-shifted to 510 nm for the nanoparticle polymer P3HT (NP) dispersions. The spectra for these dispersions also show the appearance of peaks at 520 nm, 560 nm, and 620 nm, which indicate the existence of vibronic transitions caused by the internal aggregation of the P3HT chains inside the nanoparticles.[1,2] The acquired aggregate structure with its electronic transitions of the P3HT (NPs) in dispersion is maintained when deposited in the form of films onto substrates.-- S2. Redox charge density & thickness: Figure S2 shows the cyclic voltammograms of the different P3HT film probed for different film thicknesses at a scan rate of 20 mV/s in the potential window from -0.3 to 1.1 V vs. Ag/AgCl reference electrode (RE), calibrated at 0.45 V vs. ferrocene. The surface area exposed to the electrolyte is about 1 cm2.-- S3. Stationary Transmittance: Figure S4 show the stationary transmittance curves of the P3HT film series taken at 520 nm in the oxidized and neutral state as a function of the calculated redox charge densities. Experimental data points are fitted by exponential decay functions. The P3HT film series show similar transmittance curves in the oxidized state (Figure S3A), while those in the neutral state (Figure S3B) exhibit larger deviations. The difference between transmittance values in the oxidized and neutral state then provides the contrast curve as a function of the redox charge density as shown in Figure 2 of the main manuscript.--, S4. Transmittance spectra: Figure S4 shows the transmittance spectra of P3HT-THF film acquired in its neutral and oxidized states, reflecting its magenta and transparent pale blue colors, respectively. The transmittance minimum is obtained at 520 nm for the neutral state and provides the reference value at which maximum contrast, i.e. transmittance differences between the oxidized and neutral state is calculated for the different P3HT film series.-- S5. Switching speed: The switching speed of the P3HT films has been determined following the experimental procedure described in the experimental section of the main manuscript. Here the films are submitted to potential steps of variable pulse lengths of 15, 10, 5, 2, 1, 0.5 and 0.25 s between -0.3 and 1.1 V. A representative case study for a P3HT-THF film is depicted in Figure S5.-- S6. Cycling stability: The cycling stability of the P3HT films has been determined following the experimental procedure described in the experimental section of the main manuscript. Here the films are submitted to a number of potential steps between -0.3 and 1.1 V: 300 cycles of 10 s for each step were applied. A representative case study for a P3HT-THF film is depicted in Figure S6.-- S7. Spin-coated P3HT-CHCl3 film: Spin-coating of non-nanostructured P3HT-CHCl3 dispersions, provides a continuous film coverage of the ITO substrate, as can be seen by SEM (Figure S7A) and the profilometry curve of a representative film with average film thickness of 77 nm (Figure S7B). Therefore, the electrochromic transmittance and contrast behavior at 520 nm (Figure S7C) now shows more consistent results, comparable to those of the spray-coated films of the other P3HT series. This especially refers to the optimum redox charge density and maximum contrast being achieved. Equally, t90 switching speed (Figure S7D), as well as the cycling stability (Figure S7E) in the optimum redox window reveal a behavior close to the ones of the other spray-coated film series. However, the spin-coated P3HT-CHCl3 film is prone to delamination issues (Figure S7F), compromising the mechanical integrity of the film, and thus as well the contact to the underlying substrate. This results in the non-linear enhancement of the t90 switching and the decrease of the cycle stability, beyond the established optimum redox charge densities window, as can be seen in Figure S7D and S7E, respectively. The overall electrochromic performance parameters for the spin-coated P3HT-CHCl3 series are summarized in Table S1., This work was funded by Ministerio de Ciencia e Innovación-Agencia Estatal de Investigación (MCIN-AEI, Spain) under Grant numbers PID2019-104272RB-C55/AEI/10.13039/501100011033, PID2019-104272RB-C51/AEI/10.13039/501100011033 and TED2021-129609B-I00/MCIN/AEI/10.13039/501100011033 (co-funded by European Union NextGenerationEU/PRTR). A.C-S acknowledges financial support from UPCT-Banco Santander through a research grant (“Iniciación en investigacion” Program 2021). W.M. and A.M.B acknowledge financial support from Gobierno de Aragon (DGA) under project “Grupos de Investigación Reconocidos” T03_23R. E.C. is grateful for his PhD grant from MINECO (FPI BES2017-080020) and associated European Social Funds (ESF)., Peer reviewed

Under a Creative Commons license by-nc-nd 4.0, Fig. S1. STEM (left) and TEM (right) images of Cu-N-C-Ac2 acquired using a Tecnay F30 microscope (300 kV). Fig. S2. XRD pattern of Cu-N-C-Ac2-3h obtained using the same experimental details than for Cu-N-CAc2 but with a lixiviation treatment of 3 hours. Fig. S3. N1s XPS spectra of the Cu-N-C. The material Cu-N-C-Ac1-BPT before the two-step post-treatment (acid leaching and 2nd carbonization) is also included for comparison purposes. Fig. S4. Effect of the two-step post-treatment (acid leaching and 2nd carbonization) on the nitrogen adsorption-desorption isotherms for Cu-N-C-Ac1. Fig. S5. Correlation between the Cu oxidation state and the energy position of the XANES spectrum, determined as the first maximum of the first derivative spectrum of Cu-N-C-Ac0.5 (red circle) and different copper reference compounds (black circles).Fig. S6. Variation with time of pH of the electrolyte (0.1 M KHCO3 aqueous solution) as a function of CO2 partial pressure. Balance with an inert gas for the equilibrium CO2 (aq) + H2O ↔ HCO3– + H+, pKa = 6.4, room temperature, 1 atm total pressure. Fig. S7. Oxidation of formic acid and methanol at the Pt ring of an RRDE, disc turned off, 1600 rpm. Fig. S8. Comparison of the CV in 0.1 M KHCO3 saturated with N2 (blue curve), methane oxidation voltammogram after saturation of the electrolyte with CH4 at the open circuit potential (red curve) and CH4 stripping voltammogram (Ead = 0.4 V vs. RHE, green curve) at the Pt ring of an RRDE, disc turned off, 1600 rpm. CV conditions from 0.05 to 1.6 V vs. RHE, 50 mV s-1. Room temperature. Fig. S9. CVs of the disc at 50 mV s-1 and 1600 rpm in CO2-saturated 0.1 M KHCO3 of Cu-N-C catalysts. Fig. S10. Chronoamperometric experiments at the disc for the Cu-N-C catalysts. Electrolyte: Fig. S11. CVs at the platinum ring for disc potentials reported in the legend. Table S1. Bulk chemical composition (wt. %) from EA and ICP of Cu-N-C electrocatalysts before and after the two-step post-treatment (acid leaching and 2nd carbonization). Table S2. N content (at. %) from XPS of the Cu-N-C electrocatalysts. Relative areas (at. %) of the deconvoluted peaks of N1s. Table S3. Textural parameters of Cu-N-C electrocatalysts. Table S4. CO, CH4 and C2H4 formation rate (mmol g-1 h-1) at selected applied potentials (V vs. RHE) and molar H2/CO ratio. Product analysis in a H-type electrochemical cell.  , Authors acknowledge Grant PID2020-115848RB-C21 “STORELEC” project, and TED2021-129694B-C22 “DEFY-CO2” project funded by MCIN/AEI/10.13039/501100011033. They also acknowledge LMP253_21 project funded by Gobierno de Aragón. Sara Pérez-Rodríguez thanks Grant IJC2019-041874-I funded by the MCIN/AEI/10.13039/501100011033. Ana Cristina Giménez thanks CSIC for her JAE Intro ICU 2021-ICB-04 grant. David Ríos-Ruiz acknowledges the Y2020/EMT-6419 “CEOTRES” project funded by the Comunidad Autonoma de Madrid., Peer reviewed

In this document, we provide supplementary materials for the work ‘Climate change impacts on winter chill in Mediterranean temperate fruit orchards’ by Eduardo Fernandez and co-authors. The study is published in the journal Regional Environmental Change under the doi: 10.1007/s10113-022-02006-x. We conducted this work in collaboration with researchers from northern and southern Spain, Tunisia, Morocco and Germany under the umbrella of an international project (AdaMedOr) funded by the Partnership for Research and Innovation in the Mediterranean Area (PRIMA). Compared to previous similar studies, we provide now an analysis that combines the spatial interpolation of winter chill accumulation in the Mediterranean region under future scenarios with expert knowledge regarding the impacts of climate change on temperate orchards as well as future concerns of farmers cultivating temperate species. Our approach allowed us to frame and contextualize the results of our chill estimations, potentially contributing to the development of management strategies to adapt Mediterranean orchards to future climate conditions. We offer figures that were not included in the main manuscript, as well as additional information about the weather stations used for the analysis., We conducted this work in collaboration with researchers from northern and southern Spain, Tunisia, Morocco and Germany under the umbrella of an international project (AdaMedOr) funded by the Partnership for Research and Innovation in the Mediterranean Area (PRIMA)., Weather stations used in the analysis For this study, we used 387 weather stations as primary sources of minimum and maximum temperature records between 1974 and 2020. In the following table (Table S1), we provide the name, location (coordinates) and percentage of data complete for each weather station. Climate models used in the projections In Table S2, we show the 15 climate models used in the analysis to obtain future temperature data from the ClimateWizard data base. As described in the main manuscript, we later grouped these models into “pessimistic”, “intermediate” and “optimistic” classes according to Safe Winter Chill distributions. Correction model As described in the main manuscript, we implemented a spatial interpolation and used a 3D model to correct for large errors that originated from the initial Kriging procedure. This 3D correction model (Fig. S1) consisted of the relationship between the monthly minimum and maximum temperatures in January (x- and y-axis, respectively) and the observed chill in each weather station (color surface). This allowed us to identify the combination of temperatures that was associated with a given amount of chill. We later used this combination to estimate chill values from the co-variables (mean daily minimum and maximum temperatures) from both data sources (weather stations and WorldClim) and obtain a chill correction map. Additional figures In the following figures, we show the expected change in Safe Winter Chill compared to the baseline period (median SWC across the historic simulated scenarios) for the “pessimistic” and “optimistic” climate model classes for the RCP4.5 and RCP8.5 scenarios by 2050 and 2085. As expected, major chill losses will occur under the “pessimistic” version of the RCP8.5 scenario by 2085, whereas minor changes may be expected by the near future under the RCP4.5 scenario., Peer reviewed

Under a Creative Commons license by-nc-nd 4.0., Table S1. Gas composition of ethanol at 250 and 400 ºC. Figure S2. TGA of a) corncob and b) Agave angustifolia bagasse in N2 atmosphere. Table S2. Carbon balances of the different products after HDO.Figure S2. Van Krevelen diagram comparing HDO of C-250 bio-oil with and without Mo2C/CNF based-catalyst at 350 ºC. Figure S3. Van Krevelen diagram zoom of the general influence of temperature of HDO in a Mo2C/CNF based-catalyst., This study was supported by MCIN with funding from European Union NextGenerationEU (PRTR-C17. I1) within the Green Hydrogen and Energy Program- CSIC, as part of the CSIC Interdisciplinary Thematic Platform (PTI+) Transición Energética Sostenible+(PTI-TRANSENER+), and the financial support of the I+D+i project PID2020–115053RB-I00, funded by MCIN/ AEI/10.13039/501100011033. Authors also acknowledge the financial support of DGAPA-PAPIIT UNAM through grant IN107923 “Licuefacción hidrotérmica de biomasa residual” and Fondo Sectorial CONACYT-SENER-Sustentabilidad Energética through Grant 207450 and Centro Mexicano de Innovación en Energía Solar (CeMIE-Sol), within strategic project No. 120. DT is grateful for the Juan de la Cierva Incorporación (JdC-I) fellowship (Grant Number: IJC2020–045553-I) funded by MCIN/AEI/ 10.13039/501100011033 and by “European Union NextGenerationEU/PRTR”., Peer reviewed

Supplementary Tables: Table S1. Taxonomical classification of the bacterial species isolated from solar panel samples according to 16S rDNA sequence similarity. Table S2. Summary of sequencing statistics from the 16S/18S profile analysis. Table S3. Summary of sequencing statistics from the shotgun metagenomic sequencing of solar panels 1 and 3 in 2014. Table S4. Differentially expressed proteins between day- and night-collected samples., Supplementary Figures: Figure S1. Microbial colonies growing on LB incubated at room temperature for two weeks (A). Case example of growth restoration in nearby-grown isolates under conditions of extreme pH. A local buffering of the pH of the plate is observed (B). Figure S2. Stress tests results data matrix. The growth of the 53 strains isolated from the solar panels under particular stress conditions was compared to that of a control strain (XL1-Blue E. coli strain). Green colour (from light to dark) indicates better growth than the control (from slightly to strongly better growth). Symbols ‘+’ indicate cases of growth restoration by another isolate. Two independent experiments were performed for each test. Strains are numbered according to Table S1. Figure S3. Taxonomic diversity of one of the panels (panel 1) sampled in the summer solstice of 2014 as deduced from shotgun metagenomic sequencing., Peer reviewed

Figure S1: predicted secondary structure of circRAJ31. Figure S2: analysis of the secondary structure of the 5 UTR and the first 30 bp of the coding sequence of gfp or cat. Figure S3: electrophoretic assay and sequencing results of RT-PCR after in vitro transcription. Figure S4: putative sequence of the RT-PCR misprocessing reaction. Figure S5: response of the system to varying concentrations of aTc (which controls expression of the riboregulator). Figure S6: characterization of mutant systems. Figure S7: growth of cells cotransformed with plasmids expressing circRAJ31 and CamR on plates without chloramphenicol (Cam) and on plates with chloramphenicol and different inducers. Figure S8: growth curves with a cis-repressed gene coding for chloramphenicol acetyltransferase (CamR). Figure S9: additional growth curves for different concentrations of Cam. Table S2: sequences of the primers used in this work., Peer reviewed

Under a Creative Commons license by-nc-nd 4.0, Figure S1: Thermostated electrochemical reactor with a typical three-electrode configuration. Figure S2: Deconvoluted spectra of C1s peak for a) T-GFD4.6 (100 cycles in 2 M H2SO4), b) T-GFD4.6 (immersion 2 M H2SO4 0.4 M VOSO4) and c) T-GFD4.6 (100 cycles 2 M H2SO4 0.4 M VOSO4). Figure S3: Deconvoluted spectra of O1s peak for a) T-GFD4.6 (100 cycles in 2 M H2SO4), b) T-GFD4.6 (immersion 2 M H2SO4 0.4 M VOSO4) and c) T-GFD4.6 (100 cycles 2 M H2SO4 0.4 M VOSO4). Figure S4: Deconvoluted spectra of N1s peak for a) T-GFD4.6 electrode after being immersed 5 days in 2 M H2SO4 + 0.4 M VOSO4 electrolyte b) the T-GFD4.6 activated after 100 CV cycles in 2 M H2SO4. Figure S5: Tafel plots of the a) GFD4.6-EA, GFA6 and Avcarb6 and b) T-GFD4.6-EA, T-GFA6 and T-Avcarb6 obtained from the Linear Sweep Voltammetry curves. Figure S6: Linear relationship between the logarithm of the redox peak current and the logarithm of the scan rate for the thermal pre-treated electrodes; a)T-GFD4.6, b) T-AvCarb6 and c) T-GFA6. Table S1. Table S2. Table S3., This research is part of the CSIC program for the Spanish Recovery, Transformation and Resilience Plan funded by the Recovery and Resilience Facility of the European Union, established by the Regulation (EU) 2020/2094. Support of the CSIC Interdisciplinary Thematic Platform (PTI+) Transición Energética Sostenible+ (PTI-TRANSENER+). Funding provided by the Spanish Ministry of Science and Innovation under project MAFALDA (PID2021–126001OB-C32), as well as the support of the Regional Government of Aragon to the Fuel Conversion Research Group (T06_23R)., Peer reviewed