BUSCADOR RECOLECTA

Supplementary Material 1: Supplemental Figure S1. Climate chamber in which the E. timida slugs were kept in artificial sea water in plastic cups as aquariums. The green tubes provided the air supply. Supplemental Table S1.Databases and tools which were used while operating InterProScan version 5.64-96.0 [101]. Supplemental Table S2.Table of PKS and fatty acid synthase (FAS) sequences from Torres et al. (2020) [34] including the animal species they were received from and the accession number. Supplemental Figure S2. Genome size estimation of E. timida using flow cytometry. The histogram shows the relative propidium iodide fluorescence intensity obtained after simultaneous analysis of E. timida 2C (in green) and the house cricket A. domesticus 2C as an internal standard reference (in red). The PI fluorescent dyes were excited with a solid-state laser emitting at 488 nm. The y-axis gives the counts of propidium iodide (PI) stained nuclei. The x-axis displays the relative red PI fluorescence signal. To obtain the mean relative red PI fluorescence signals, the peaks were enclosed by line segments. The percentages in brackets are the portions of all events in the histogram enclosed by the respective line segments. Supplemental Table S3. Genome size estimates from two individuals of E. timida. The measured individual is given in brackets. Chopping buffer was prepared as described by Galbraith et al. (1983) [152]. Propidium iodide was used as a fluorescent dye. We used the house cricket A. domesticus as standard reference (genome size: 2000 Mb). Supplemental Figure S3. K-mer profile and estimates based on HiFi reads. Supplemental Table S4. Sacoglossan heterozygosity values. The heterozygosity values from all species except for E. timida, were inferred by Theisen & Jensen (1991) [128]. Supplemental Tables S5. PacBio ultra-low library preparation based on PCR amplification with KOD Xtreme™ Hot Start DNA Polymerase (Merck). Supplemental Table S6. Sequencing output and subread mean length of the PacBio low- and ultra-input libraries. Supplemental Figure S4. HiFi read length distribution and statistics. Standard PacBio ultra-low input libraries are listed as SMRT1 and SMRT2. PacBio ultra-low libraries amplified with KOD polymerase are shown as SMRT3 (Sequel IIe) and SMRT4 (Revio). N50 values are presented in bp. Supplemental Table S7. FCS-GX contamination summary. Supplemental Table S8. FCS-GX action summary. Supplemental Figure S5. Blobplot of the assembly after polishing and purging. At this stage of the assembly process, contamination filtering with FCS was already conducted. Supplemental Table S9. Blobtools taxonomic assignment. The table shows all contigs classified as Chlorophyta by “bestsumorder”, which were filtered out among others. Sequences marked with asterisk were categorized as “HICOV” by purge_dups. Supplemental Figure S6. Blobplot of the final genome assembly. Supplemental Figure S7: Maximum likelihood phylogenetic tree of FAS, PKS1 and PKS2 transcripts from E. timida, E. chlorotica, E. diomedea and P. ocellatus. For the alignment the transcriptomic data from the sequences listed in Table 5 and Supplemental Table S2 were used. The branches are labelled with their length and scaled according to the number of substitutions per site. The percentage of trees in which the associated data clustered together is shown next to the branches. The transcript of EtPKS1 was manually constructed based on sequence homology to EcPKS1, EdPKS1 and PoPKS1 as described previously. The transcripts from E. timida are labelled with an asterisk. Supplemental Table S10. Number of blast hits with taxid of Acetabularia acetabulum or Ulva compressa against contigs of the polished E. timida genome assembly. Supplemental Table S11. Number of blast hits for targets with a taxid of Acetabularia acetabulum or Ulva compressa. All target sequences originate from a chloroplast. Supplemental Table S12. Section of the agp file from scaffolding including the 5 Chlorophyta sequences showing all resulting scaffolds containing these Chlorophyta sequences. Except for splitting one of the sequences, none was linked to other nuclear sequences of the E. timida genome assembly. Supplemental Figure S8. The genes encoding EtFAS, EtPKS1 and EtPKS2 are annotated in the genome of E. timida. The exons are labelled in black on the excerpt of the genomic sequence. The arrows present the transcript of the a) EtFAS, b) EtPKS1 and c) EtPKS2. The domains of the enzymes are presented in bubbles below the arrow. The gene encoding EtPKS1 was annotated manually based on sequence homology with the EcPKS1. The label with an x0 indicates an inactive domain. Supplemental Figure S9. Isotopic patterns of the putative polypropionates produced by E. timida and identified by HPLC-ESI-HRMS analysis. Supplemental Figure S10. HPLC-MS data of E. timida extracts. Base peak chromatogram (BPC) and extracted ion chromatograms (EICs) of polypropionates shown in fig. 4. The BPC shows all detected ions present in the crude extract and the EICs show the peaks corresponding to putative polypropionates. Supplemental Table S13. Result of the polypropionate blast search in the annotations of E. timida,E. chlorotica, E. diomedea and P. ocellatus. Both unfiltered (F-) as well as filtered (F+) blast hits are shown. Supplemental Figure S11. Blobplot of the assembly after removing sequences identified as contamination and before HiC scaffolding., Senckenberg Gesellschaft für Naturforschung (SGN) (3507), Peer reviewed

Additional file 1., Peer reviewed

Additional file 2., Peer reviewed

Additional file 3., Peer reviewed

The PDF file includes: Supplementary Figures 1-9. -- Table S1. -- Supplementary Figure 1. Prediction of the capacity to display flagellar motility in bacterial taxa based on presence/absence information of 21 genes involved in flagellar assembly. -- Supplementary Figure 2. Taxonomic distribution of genomes with empirically determined capacity for flagellar motility that were used as training data for a boosted regression machine learning model to predict the capacity for flagellar motility based on the presence/absence of 21 flagellar genes. -- Supplementary Figure 3. Predictive accuracy across phyla of a boosted regression machine learning model for the prediction of the capacity for flagellar motility in any given bacterial genome based on the presence/absence of 21 flagellar genes. -- Supplementary Figure 4. Prevalence of flagellar motility across bacterial families containing more than 100 high-quality genomes in the Genome Taxonomy Database (GTDB r207; [27]). -- Supplementary Figure 5. Variation in flagellar motility status across taxonomic ranks. -- Supplementary Figure 6. Distribution of the total number of 16S rRNA gene copies per genome and genome size across the 6 phyla with even proportions of taxa predicted to be flagellated and non-flagellated. -- Supplementary Figure 7. Estimated prevalence of flagellar motility in bacterial communities from 9 soil depth profiles collected across the USA (Surface, 0-20cm; Subsurface, 20-90cm, N = 66; [64]). -- Supplementary Figure 8. Difference in the prevalence of flagellar motility in rhizosphere and bulk soil bacterial communities collected from citrus species across the globe (N = 10; [67]). -- Supplementary Figure 9. Taxonomic composition of the Amplicon Sequence Variants (ASVs) that responded to glucose amendment in soil. -- Table S1. Summary of the shotgun sequencing reads generated from a soil incubation experiment where soils were amended or not with glucose over a 117-day period [70]., Peer reviewed

Supplementary Table 1. Details of the study sites: habitat type, number of transects, total length of transects at each study site and habitat type (m), water depth (m), and latitude and longitude coordinates. -- Supplementary Table 2. List of fish species recorded within the five study sites in the island of Menorca during May 2022. Tropic category (macrophagous, MA; mesophagous, ME; microphagous, MI; omnivorous, O; and herbivorous) is indicated. a and b parameters refer to the weight-length relationships used to estimate fish biomass. -- Supplementary Table 3. PERMANOVA main test results in the effects of habitat type (seagrass meadows, seagrass escarpments and rocky substrate) on fish abundance (number of fish m-2 ), total number of species (number of species m-2 ), total biomass (g m-2 ), Shannon´s diversity Index (H’), structural complexity, trophic category (macrophagous, mesophagous, microphagous, omnivorous and herbivorous), and fish size categories at 5 cm intervals. df = degrees of freedom, SS = sum of squares, MS = mean squares, P(perm) = p value using permutations. Significant P(perm) values < 0.05 are indicated in bold. -- Supplementary Table 4. Post-hoc pair-wise comparisons showing the effect of habitat type (seagrass meadows, seagrass escarpments and rocky substrate) on fish abundance (number of fish m2 ), total number of fish species, trophic category (macrophagous, mesophagous, microphagous, omnivorous and herbivorous), fish size categories, total number of species (number of species m-2 ), total biomass (g m-2 ), Shannon´s diversity Index (H’), and structural complexity index. P(perm) = p value using permutations. Significant P(perm) values., Peer reviewed

This in situ dataset collected during the ARCTIC-MON field campaign will be linked in the Barcelona Polar Lab (BPL) webpage: https://polar.icm.csic.es, The ARCTIC-MON campaign conducted at the Canadian High Arctic Research Station (CHARS) collected brightness temperatures from L-band radiometric measurements of sea ice, using the ARIEL portable L-band radiometer. This data is key to improve the emission modelling that is required for the development of the sea ice thickness products from remote sensing observations, This project is funded with the ARCTIC-MON project (PID2021-125324OB-I00) and also Programación Conjunta Internacional project called ICE-MOD, with reference PCI2019-111844-2, from the Agencia Estatal de Investigación in Spain. With the institutional support of the ‘Severo Ochoa Centre of Excellence’ accreditation (CEX2019-000928-S), Brightness temperature measured over sea ice at each stop done during each days transects. Each stop is labelled within each day’s measurements, Peer reviewed

Soil characteristics, geographical coordinates, date of collection and abundance (read number) OTU table of arbuscular mycorrhizal fungi present in the island Peberholm (Denmark) from 2010 to 2016., This project has received funding from the European Commission's Horizon 2020 research and innovation program (European Union) under the Marie S. Curie grant agreement No 708530 – DISPMIC. ALG was also funded by the grant ref. RYC2022-038499-I, funded by MCIN/AEI/10.13039/501100011033 and FSE+. JT would like to thank financial supports from Villum Young Investigator (No. VIL53048) and Danish National Research Foundation within the Center for Volatile Interactions (VOLT, DNRF168)., Peer reviewed

Under a Creative Commons license BY-NC-ND 4.0, Figure S1 XRD analysis to the solid obtained from the hydrothermal reaction (prior to calcination at 1000ºC); Figure S2 POR-Hg results obtained for the catalyst support (MCaO); Figure S3 POR-Hg results obtained for the 5 %Ni-MCaO catalyst; Figure S4 N2 adsorption isotherms obtained for the catalyst support (MCaO, red lines) and the 5 %Ni-MCaO catalyst (5%Ni-MCaO, black lines)Figure S5 XRD diffractogram of the 5%Ni-MCaO catalyst after being used under Test 1 conditions along 6 reaction-regeneration cycles., This work has received funding from project TED2021–131768B-I00 funded by MCIN/AEI/10.13039/501100011033 and by the “European Union NextGenerationEU/PRTR”. Carlos Ordóñez thanks the grant PRE2019–090024 funded by MCIN/AEI/ 10.13039/501100011033 and EI ‘‘ESF Investing in your future”. Also, authors acknowledge the partial financial support of the Regional Government of Aragón through the Research Groups Funding Program., Peer reviewed

Additional file 1: Table S1. Chemical analysis of the water used for the bottled water condition. Table S2. Parameters, description, notation and formula of the lifetables. Table S3. Pairwise comparisons of each group in Kaplan Meier survival analysis of the immature survival (from L1 to adult emergence). Table S4. Summary table of the Cox analysis with increasing salinity concentrations and distilled condition as control. Table S5. Summary table of the Cox analysis describing the effects of distilled vs. bottled water in the survival., ’la Caixa’ Foundation MoSquITo project Consejo Superior de Investigaciones Cientificas (CSIC), Peer reviewed