Microbiota Dynamics in Lionfish (Pterois): Insights into Invasion and Establishment in the Mediterranean Sea

Microbiota Dynamics in Lionfish (Pterois): Insights into Invasion and Establishment in the Mediterranean Sea | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Microbiota Dynamics in Lionfish (Pterois): Insights into Invasion and Establishment in the Mediterranean Sea Dalit Meron , Maya Lalzar , Shevy Bat-Sheva Rothman , Yael Kroin , Elizabeth Kaufman , Kimani Kitson-Walters , Tal Zvi-Kedem , Eli Shemesh , Rami Tsadok , Hagai Nativ , Shai Einbinder , Dan Tchernov doi: https://doi.org/10.1101/2025.02.03.635474 Dalit Meron 1 Morris Kahn Marine Research Station, University of Haifa , Sdot Yam, Israel 2 Department of Marine Biology, Leon H. Charney School of Marine Sciences, University of Haifa , 3498838, Israel Find this author on Google Scholar Find this author on PubMed Search for this author on this site For correspondence: dmeron{at}univ.haifa.ac.il Maya Lalzar 3 Bioinformatics Services Unit, University of Haifa , Israel Find this author on Google Scholar Find this author on PubMed Search for this author on this site Shevy Bat-Sheva Rothman 4 The Steinhardt, Museum of Natural History, Tel Aviv University , Tel Aviv, Israel Find this author on Google Scholar Find this author on PubMed Search for this author on this site Yael Kroin 1 Morris Kahn Marine Research Station, University of Haifa , Sdot Yam, Israel Find this author on Google Scholar Find this author on PubMed Search for this author on this site Elizabeth Kaufman 5 Israel Aquarium, Tisch Family Zoological Garden in Jerusalem , Jerusalem, Israel Find this author on Google Scholar Find this author on PubMed Search for this author on this site Kimani Kitson-Walters 6 Caribbean Netherlands Science Institute (CNSI) , St. Eustatius, Caribbean Netherlands 7 NIOZ Royal Netherlands Institute for Sea Research, Utrecht University , Den Burg, Texel, Netherlands Find this author on Google Scholar Find this author on PubMed Search for this author on this site Tal Zvi-Kedem 2 Department of Marine Biology, Leon H. Charney School of Marine Sciences, University of Haifa , 3498838, Israel Find this author on Google Scholar Find this author on PubMed Search for this author on this site Eli Shemesh 2 Department of Marine Biology, Leon H. Charney School of Marine Sciences, University of Haifa , 3498838, Israel Find this author on Google Scholar Find this author on PubMed Search for this author on this site Rami Tsadok 1 Morris Kahn Marine Research Station, University of Haifa , Sdot Yam, Israel 2 Department of Marine Biology, Leon H. Charney School of Marine Sciences, University of Haifa , 3498838, Israel Find this author on Google Scholar Find this author on PubMed Search for this author on this site Hagai Nativ 1 Morris Kahn Marine Research Station, University of Haifa , Sdot Yam, Israel 2 Department of Marine Biology, Leon H. Charney School of Marine Sciences, University of Haifa , 3498838, Israel Find this author on Google Scholar Find this author on PubMed Search for this author on this site Shai Einbinder 1 Morris Kahn Marine Research Station, University of Haifa , Sdot Yam, Israel 2 Department of Marine Biology, Leon H. Charney School of Marine Sciences, University of Haifa , 3498838, Israel Find this author on Google Scholar Find this author on PubMed Search for this author on this site Dan Tchernov 1 Morris Kahn Marine Research Station, University of Haifa , Sdot Yam, Israel 2 Department of Marine Biology, Leon H. Charney School of Marine Sciences, University of Haifa , 3498838, Israel Find this author on Google Scholar Find this author on PubMed Search for this author on this site Abstract Full Text Info/History Metrics Supplementary material Preview PDF Abstract Lionfishes ( Pterois spp.), originally native to the Indo-Pacific and Red Sea, have become one of the most invasive marine species globally, including the recent establishment in the Mediterranean Sea. This study investigates the microbiota of lionfish to explore its potential role in their invasion success and establishment. Using high-throughput sequencing and microbiota analyses, we characterized the species-specific core microbiome and identified habitat-specific markers across different regions (Red Sea, Mediterranean Sea, Caribbean, and aquarium populations) and organs. Focusing on the Mediterranean invasion, we tracked lionfish distribution and population dynamics along the Israeli coastline from 2017 to 2023, monitoring size, seasonal trends, and depth preferences. Our findings reveal that lionfish initially established themselves in deeper waters before expanding to shallower habitats, with a gradual increase in population size and body length over time. From a microbial aspect, we compared the microbiota of lionfish organs and identified a similar pattern ( Photobacterium) , to Earlier Lessepsian migrants fish species. This study provides novel insights into the interactions between microbiota and host ecology, shedding light on the mechanisms that may support the successful invasion. This study contributes to the understanding of lionfish invasion dynamics in the Mediterranean. It highlights the microbiota as an integral component for studying the ecological and biological mechanisms underpinning invasive species’ success and establishment of lionfish. Introduction Lionfish (LF) species ( Pterois volitans and P. miles ), are reef fish native to the Indo-Pacific and Red Sea, are considered among marine ecosystems’ most notorious invasive species ( Côté & Smith 2018 ). They have become significant and successful invaders with global distribution, including the northwestern Atlantic, the eastern coast of the USA, Bermuda, the Caribbean region and the Gulf of Mexico ( Schofield, 2009 ; Côté & Smith, 2018 ), Brazil ( Ferreira et al., 2015 ) and more recently, the Mediterranean Sea (Kleitou et al., 2016). In the Mediterranean Sea, P. miles were first recorded near Israel in 1991 ( Golani & Sonin, 1992 ), and the invasion notably accelerated from the mid-2010s and onwards. The Suez Canal has served as a key pathway for this introduction, also termed “Lessepsian migration” facilitating the species’ establishment in new habitats ( Bariche et al., 2017 ; Stern et al., 2019 ). The unprecedented speed and scope of their spread highlight their adaptability and the significant impact they impose on invaded ecosystems. Several biological and ecological factors contribute to the LF’s invasion success ( Côté & Smith 2018 ; Bottacin et al., 2024). Lionfish are generalist predators whose diet includes a wide variety of fish and invertebrates. A stomach content analysis has revealed that LF in the Mexican Caribbean preyed on nearly 50 different local fish species ( Arredondo-Chávez et al., 2016 ), underscoring their lack of dietary selectivity and facilitating their adaptation to new environments. Their early maturation (within the first year) combined with a robust strategy (asynchronous egg release) and fast growth, contribute to the rapid increase in LF density and their successful establishment in invaded areas ( Morris et al., 2011 ; Côté et al., 2013 ; Albin 2013). Additionally, the absence of natural predators in the invaded regions allows LF populations to grow unhindered, thus significantly disrupting the local marine ecosystems. This invasion leads to a decline in local species and changes biodiversity and community structure, affecting the reef ecosystems’ resilience ( Albins & Hixon, 2008 ; Green et al., 2012 ; Côté et al., 2018). Invasive species are often characterized by their remarkable adaptability, which enables them to thrive in novel environments with different abiotic and biotic conditions compared to their native habitats. The microbiota may influence various host functions, including immune responses, nutrient acquisition, and environmental resilience, providing an invasive species with significant advantages in unfamiliar territories ( Aires et al., 2016 ; Ronanki, 2024 ; Coats et al., 2014 ). Another aspect concerns the symbionts of invasive species in the new environments. Numerous studies have shown that various stressors, such as temperature changes, pollution, and fluctuations in resource availability, can potentially impact the dynamics of bacterial invasion in aquatic ecosystems ( Liu et al., 2012 ; Bagra et al., 2022 ). These stressors may not only affect the composition of bacteria and disrupt established symbiotic relationships but also create new niches that invading microorganisms can exploit ( Amalfitano et al., 2015 ; Cavicchioli et al., 2019 ; Paerl et al.,2016 ). Although research has explored the effects of environmental changes on microbiota, the specific role of microbiota in supporting the survival and establishment of invasive species remains inconclusive. The role of symbiotic bacteria in invasive species has been investigated mainly in terrestrial ecosystems ( Lu et al ., 2016 ; Bang et al ., 2018 ). Research focused on insects for example, has demonstrated that altering the microbiota can influence various traits, including: suppressing native competitors ( Vilcinskas et al., 2013 ), enhancing nutrient acquisition ( Adams et al. 2011 ), promoting pest resilience ( Brown et al., 2014 ) and increasing fecundity and survival rates ( Himler et al., 2011 ). These traits contribute to establishing and spreading invasive species within new environments. However, only a few studies have demonstrated the contribution of the microbiota to invasion within marine ecosystems. Among them, a study on the invading alga Caulerpa racemosa showed that an increase in bacteria involved in sediment biogeochemical processes contributed to the invasion’s success by displacing native seagrasses ( Aires et al., 2013 ). Another study demonstrated that symbiotic relationships between chemoautotrophic bacteria and marine invertebrates can influence metabolic capabilities and expand the adaptation to a variety of ecological niches for the host ( Lee et al., 1999 ). Thus far, only a few studies have examined the LF microbiota aspect and have mainly focused on the Indo-Pacific (native) and western Atlantic (invaded) regions ( Stevens & Olson, 2013 ; Stevens & Olson, 2015 ; Stevens et al., 2016 ). However, these were conducted about a decade ago, using methods that were less advanced and comprehensive than those available today, such as high-throughput amplicon sequencing and metagenomics Invasive species may serve as ideal models for investigating the host-microbe-environment interactions. By comparing the microbiota in native and invaded environments, we can investigate three aspects: 1) how the new environmental conditions shape and change the composition and diversity of the bacterial communities; 2) what is the stable core microbiota (that persists despite the changing conditions), and 3) to examine whether the microbiota impacts the success of invasion. This study compares the microbiome aspect of LF from different origins: The Mediterranean and the Caribbean (invaders), the Red Sea (native), and local aquarium. In addition, we chose to focus on the LF in the Mediterranean Sea, a relatively recent and successful invader, making it an ideal case study for examining its invasion, establishment, and interactions with microbial symbionts. To this end, fish surveys, which included fish counts and visually estimated lengths, were conducted over a five-year period. In parallel, for comparison, microbial profiling and trophic position (TP) of the two earlier Lessepsian invasive species, Sargocentron rubrum and Siganus rivulatus , were also described. Methods Fish sampling Between 2019 and 2022, a total of 90 fish specimens were sampled, comprising five individuals of S. rubrum , nine of S. rivulatus , and 79 of Pterois spp. (lionfish) (Table S1). Lionfish were collected from the Gulf of Eilat in the Red Sea (RS), the Mediterranean coast of Israel (MS), along the Israeli coast, Sint Eustatius in the Caribbean Sea (CA) and from captive populations in the Israel Aquarium in Jerusalem (AQ_J) ( https://www.israel-aquarium.org.il ). LF from the latter were originally captured off the coast of Kenya and kept in the aquarium since 2018. Fish specimens from the wild (from MS, RS and CA) were caught using spearguns and transported on ice to the laboratory. In contrast, captive fish specimens were collected with (hand) nets, sampled immediately, and returned to the aquarium. Size (total length) and sex (when applicable) indices were recorded for each individual. Sterile swabs were used to collect samples from four organs of each specimen: cloaca, mouth, gills, and skin. Surrounding water was also sampled, 1.5 liters filtered using a single-use “Nalgene RapidFlow Filters” 0.2μm (Thermo Scientific, cat no. 566–0020, Israel). The swabs and filters were stored at −20°C for further analysis. In total, 378 (365 fish and 13 water) samples were collected and analyzed (Table S2). Lionfish surveys in the Mediterranean Sea As part of the Morris Kahn Marine Research Station (MKMRS) monitoring program, fish surveys were conducted by SCUBA in central (Sdot-Yam) and northern (Achziv) Israel during the Spring (April to May) and Fall (October to November) at depths of 10, 25, and 45 meters. At each depth, an area of 250 square meters was surveyed. Surveys included fish counts and visually estimated lengths. The data on LF from these surveys were analyzed for the period of 2017 - 2023. The fish survey method and data are available at https://med-lter.haifa.ac.il/database-hub/ . DNA extraction, PCR amplification and amplicon sequencing DNA was extracted from all samples using the DNeasy PowerSoil Pro Kit (Qiagen, www.qiagen.com ) following the manufacturer’s instructions. Partial sequences of the 16S rRNA gene at the V4 hypervariable region were amplified using the primers 518F (ACCAGCAGCCGCGGTAATACG) and 806R (GGACTACNVGGGTWTCTAAT) (based on Apprill et al .., 2015 ; Parada et al., 2016 with few modifications) that contained 5’ common sequence tags (CS1 and CS2) ( Moonsamy et al., 2013 ). Amplicons were generated using a two-stage PCR amplification protocol described by Naqib et al. (2018) . Cycling conditions for the first stage PCR were 94°C for 15s, 50°C for 20s, and 72°C for 20s for 15 cycles, with an additional 15 cycles with 62°C for 15s (annealing stage). Subsequent steps were carried out at the Genome Research Core (GRC) within the Research Resources Center (RRC) at the University of Illinois at Chicago (UIC). A second amplification was performed for each sample, with a separate primer pair with a unique 10-base barcode obtained from the Access Array Barcode Library for Illumina (Fluidigm, South San Francisco, CA; Item# 100-4876). Cycling conditions were 95 °C for 5 minutes, followed by eight cycles of 95 °C for 30 sec, 60 °C for 30 sec and 72 °C for 30 sec. The amplified barcoded PCR products were pooled and purified using an AMPure XP cleanup protocol (0.6X, vol/vol; Agencourt, Beckmann-Coulter). With a 15% phiX spike-in, the pooled libraries were loaded onto an Illumina MiniSeq mid-output flow cell and sequence (2×153 bases paired-end reads). Barcode sequences were used for sequence read de-multiplexing of raw data, which was then recovered as FASTQ-formatted files. Raw sequence data is available in the NCBI SRA database under PRJNA1213691. Sequence processing The Dada2 pipeline ( https://www.nature.com/articles/nmeth.3869 ) (dada2 package version 1.20.0) was used for sequence data processing. The processing was conducted for each set of samples and each run separately for the following steps: Sequences were filtered and trimmed for quality using the ‘filterAndTrim’ command with the parameters maxN set to zero, maxEE set to 2, trimLeft set to 21 and 20 bases for the forward and reverse reads, respectively. The sequence error estimation model was calculated using the ‘learnErrors’ option, and then, the dada2 algorithm for error correction was applied using default parameters. Forward and reverse reads were then merged with minimum overlap set at eight bp. Suspected chimaera was detected and removed using the command ‘removeBimeraDenovo’, and a count table was produced. To obtain a taxonomic assignment for each ASV, ASV sequences were aligned to the ARB-Silva small subunit rRNA database (version Silva_nr_138.1) using the command ‘assignTaxonomy’ with default parameters. However, minBoot set at 80% (Table S2). Data analysis Microbiota analysis procedures were conducted in R (version 4.3.1). For analysis of microbiota composition, counts data had been normalized by the cumulative sums squares (CSS) method using the R package ‘metagenomSeq’ (version 1.42.0). Normalized counts were then summed sample-wise, based on taxonomic assignment, up to the genus taxonomy level. Core microbiota composition was determined based on the prevalence of the different genera detected. A genus was considered a core member if its prevalence at each site tested was above 25%. To test the effect of the collection site (Sea) or organ on the composition of the microbiota, a permutational analysis of variance (PERMANOVA) test was performed using the ‘adonis2’ function in the R package ‘vegan’ (version 2.6.4). Non-metric multidimensional scaling (NMDS) ordination was calculated using the R package ‘vegan’ based on Bray-Curtis dissimilarities. To extract microbiota markers, Linear Discriminant Analysis (LDA) Effect Size (LEfSe) was calculated using the function ‘run_lefse’ from package ‘microbiomeMarker’ (version 1.6.0). The association was significant for Bejnamini-Hochberg adjusted P value <0.05 and LDA effect size of 2. For the calculation of the effect of sex, PERMANOVA was performed separately, using 112 samples corresponding to 34 fish specimens for which sex could be determined. Additionally, the contribution of sex to variance was tested by canonical correspondence analysis (CCA) in ‘vegan’ considering the model microbiota∼sex+site+organ. CCA model and items were considered significant at P <0.05. For alpha-diversity parameters calculations, the count data was subsampled to 4,000 reads per sample using the function ‘rrarefy’ from the R package ‘vegan’ package (version 2.6.4). Using the rarified data, Shannon H’s index of diversity, and the Simpson’s index for evenness and richness, calculated as observe number of genera were calculated with ‘vegan’. To test the effect of site, organ and their interaction on alpha diversity parameters, the aligned rank-transformed (ART) ANOVA test was applied, with the R package ‘ARTool’ (version 0.11.1). Post-hoc pair-wise tests were also conducted with package ‘ARTool’ using the ‘art.con’ function and the Benjamini-Hochberg method for P value adjustment against false discovery. Results Lionfish microbiota composition Lionfish microbiota composition was compared among four sites (MS, RS, CA and AQ_J) and among the four organs: mouth, gills, cloaca, and skin (Total 307 samples; Table S1). To define the LF core microbiota, we examined the prevalence of the different genera identified and considered those with prevalence > 25% at all of the examined sites. Based on these criteria, 18 bacterial genera were recognized as the LF core microbiota ( Fig. 1A , Table S2). The most prevalent genera were Photobacterium , Vibrio , Pseudomonas and Shewanella (prevalence above 80% of the samples), all belonging to the class Gammaproteobacteria. While a robust core microbiota was thus suggested, an inspection of the dominant genera in each organ at each site (> 4% mean relative abundance) indicated major differences in composition between sampling sites ( Fig. 1B ). The contribution of the site to variation in microbiota composition, as assessed by PERMANOVA was highly significant ( P < 0.001, R 2 =0.133). Pair-wise comparisons further indicated unique microbiota composition for each site (Table S3). This was corroborated by NMDS ordination ( Fig. 2A ). Linear decomposition model (LDA) effect size analysis was used to identify key genera as markers of each site (Table S4_LDA). The most noted genera identified (LDA score >3) were Psyhromonas (AQ_J), Pseudomonas (CA), Photobacterium and Vibrio (MS), Cetobacterium , including Clostridiaceae and Caedibacter (RS) ( Fig. 2B ). The contribution of the organ to variation was examined for each site separately, using PERMANOVA. Organ effect was significant for LF sampled at site RS ( P <0.001, R 2 =0.087) and CA ( P <0.001, R 2 =0.066), as was indicated from inspection of composition ( Fig. 1B ). Post-hoc pairwise tests for organs in RS concluded a significant difference between each pair, with the greatest effect between the gills and the cloaca microbiota (gills-cloaca: P < 0.001, R 2 =0.089; gills-mouth: P =0.005, R 2 =0.041; gills-skin: P < 0.001, R 2 =0.072; mouth-cloaca: P <0.001, R 2 =0.066; mouth-skin: P =0.006, R 2 =0.046; skin-cloaca: P= 0.022, R 2 =0.039) while in site CA significant difference was observed only between skin and cloaca ( P <0.029, R 2 =0.039). Download figure Open in new tab Figure 1: Composition of lionfish microbiota. Lionfish microbiota was determined by amplicon sequencing of 16S rRNA gene fragments from total DNA extracted from 307 samples, representing the gills, cloaca, skin and mouth of specimens collected at the Mediterranean Sea (MS), Red Sea (RS), Caribbean Sea (CA) and the Jerusalem aquarium (AQ_J). ( A) Core microbiota composition. Includes genera for which prevalence was >25% at each of the sampled sites. ( B) Composition of dominant genera (relative abundance >5% at least at one site). Download figure Open in new tab Figure 2: Variation in lionfish microbiota composition among sites. ( A) Non-metric multidimensional scaling (NMDS) ordination, based on Bray-Curtis distances (Stress K=3 =0.16). ( B) Microbial markers for sites, based on linear decomposition model effect size (LEfSe) analysis. Presented are significant markers (Benjamini-Horchberg adjusted P value 2.5. The effects of site and organ were also reflected in microbiota diversity, examined at the genus level of taxonomy. ART-ANOVA test confirmed the effect of the site, without significant interaction for Shannon H’ and Simpson indices ( Table 1 ). View this table: View inline View popup Download powerpoint Table 1: Effect of site and organ on lionfish microbial community structure. Alpha diversity indices Shannon H’, Simpson and richness (the number of species observed (Sobs)) were calculated and compared by aligned rank transformed ANOVA. Richness, calculated as the observed number of genera, was significant for site and organ as well as their interaction Based on pairwise comparisons among sites the CA site had the lower diversity compared to all other sites, and MS diversity was lower than that of AQ_J. Pairwise comparison among organs calculated significantly lower diversity for the cloaca compared to the mouth and the skin (Table S5). Another aspect examined for effect on microbiota composition was sex. The sex could be identified for a subset of 34 fish specimens (112 samples). PERMANOVA test indicated no significant effect for sex. However, canonical correspondence analysis (CCA) did significantly link sex to microbiota composition (χ 2 =0.31, P =0.005) and explains ∼1.1% of inertia ( Fig. 3 ). Download figure Open in new tab Figure 3: Contribution of sex to variation in lionfish microbiota composition. Canonical correspondence analysis (CCA) analysis determined sex as a significant factor ( P <0.05) contributing 2.8% of variance in microbiota composition. Presented is CCA ordination plot. Lessepsian invasive fish (RS - MS) Interestingly, when analyzing the bacterial overlap at the genus level across paired organs, it was found that over 67% of the bacteria identified in MS were shared between the two regions. Conversely, a notably higher percentage of bacteria was unique to RS, with 47-41%, compared to only 22-23% of unique bacteria in MS ( Fig. 4 ). Noting the major difference in the relative abundance of the Photobacterium in the MS vs. RS samples ( Figs. 1B and 2B ), we extended our microbiota comparison to additional Lessepsian invasive fish species. We chose Sargocentron rubrum and Siganus rivulatus , known as Earlier Lessepsian migrants ( S. rubrum reported by Haas & Steinitz in 1947 and S. rivulatus by W. Steinitz in 1927 ), and characterized their microbiota. Overall relative abundance of Photobacterium ranged differently levels among species (0.6% - 59.2% for S. rubrum ; 0% - 24.2% for S. rivulatus ). Nevertheless, regardless of the organ, Photobacterium was much more dominant in MS specimens compared to RS ones ( Fig. 5 ). Interestingly, the Photobacterium relative abundance was negligible (less than 1%) in the MS water samples (Table S2). Download figure Open in new tab Figure 4: Venn diagrams presenting shared and unique microbial genera among sites for different organs of Lionfish. Numbers in brackets represent the percentage from total number identified genera at the specific site and organ. Download figure Open in new tab Figure 5: Dominance of Photobacterium across species, sites and organs. Relative abundances of Photobacterium are presented as pie charts. Lionfish MS - invasion and establishment Based on the fish surveys, we analyzed the LF data and compared the number of fish observation and their size by site, year, season, and depth ( Fig. 6 ). LF were first recorded in fish surveys in 2017 in the northern site (Achziv). The first sightings of LF in the Sdot-Yam area occurred only two years later. At both sites, initial observations of LF were recorded from the deepest transect (45 m), with sightings at shallower depths (25 m and 10 m) appearing only two to three years later. A consistent upward trend in LF population size was observed at Achziv. In contrast, no such trend was noted at Sdot-yam ( Fig. 6A ). Nevertheless, at both sites, the majority of observations were in the fall ( Fig. 6B ). Download figure Open in new tab Figure 6: Lionfish establishment in the Mediterranean Sea (2017-2022). Numbers of observations of lionfish at different seafloor depths along the Achziv and Sdot-Yam transects (A) or at different seasons (B) . Taken from MKMRS database https://med-lter.haifa.ac.il/database-hub/ Fish surveys also assessed the approximate size of the fish (in terms of length) and indicated a substantial increase in size observation between 2017 and 2023 ( Fig. 7A ). Length measurement of collected fish specimens from MS (2019-2021), RS (2020- 2022) and CA (2021) confirmed that during the first years of invasion, LF at MS were smaller in size compared to RS or CA ( Fig 7B ). Considering that smaller sizes may be accounted for by the quality of nutrition, we examined the TP value. In LF, the value of TP was the lowest in the CA. However, no differences were observed between RS and MS ( Fig. 7C ). Interestingly, the opposite trend was found for the two other Lassepsian invaders examined, S. rubrum and S. rivulatus, TP values were lower in MS compared to RS. Download figure Open in new tab Figure 7: Lionfish establishment at the Mediterranean Sea (2017-2023). (A) Size estimate distribution from visual fish surveys. (B) Size distribution of captured lionfish specimens at different sites. (C) Trophic position of captured lionfish specimens, calculated based on Trophic isotopic discrimination of amino acids (glutamic acid and phenylalanine) ( Takizawa & Chikaraishi, 2024 ). Discussion Lionfish are widely regarded as one of the most invasive and successfully established species in the marine environment. Originally native to the Indo-Pacific Ocean and the Red Sea, they have spread across the entire tropical western Atlantic, including the Brazilian coast and the Mediterranean Sea ( Côté & Smith 2018 ; Soares et al., 2022, 2023; Kletou et al., 2016 ). As previously described, various parameters such as the ability to survive, reproduction rate, the absence of predators, and its adaptation to a wide range of conditions, contributed to the success of its invasion and establishment ( Bottacini et al., 2024 ). In this study, we focused on the microbiome aspect of LF, aiming to characterize the core microbiome shared by the species and to examine differences and variations across different habitats. We also followed the process of the Lessepsian migration and the establishment of LF populations in the Mediterranean Sea. Lionfish microbiota The microbiota of invasive species is crucial to their success and interactions in new ecosystems. Microbial diversity and functional stability can influence how likely an invasive species is to establish itself ( Amalfitano, 2015 ; Huang 2019 ). We characterized the core microbiome of LF based on samples collected from four sites, and various organs. The Gammaproteobacteria class exhibited the highest relative abundance within the core microbiome, with dominant genera such as Photobacterium , Vibrio , Pseudomonas , and Shewanella ( Fig. 1A and 2C ) . The information on the microbiota of LF is very limited; however, Stevens and Olson (2015) and Stevens et al. (2016) reported similar findings regarding the bacterial groups in LF from the Indo-Pacific and Western Atlantic. Notably, Pseudoalteromonas and Photobacterium , identified as key markers in the invaded regions (CA and MS, respectively, Fig. 1B ) are known to exhibit broad antibacterial activity. Compared to that of local fish, this robust antibacterial capability may provide LF with a competitive advantage, enhancing their resilience and success in new habitats ( Stevens et al. 2016 ). Various environmental parameters, such as temperature, pH, and salinity, have been shown to significantly influence the diversity and abundance of microbiota ( Thompson et al., 2017 ; Coutinho et al., 2015 ; Hicks et al., 2018 ). In this study, we observed notable differences in microbiota composition across the analyzed sites, with each site characterized by unique bacterial markers. Identifying site-specific markers may facilitate tracing the origin of LF habitats, providing a valuable tool for monitoring their distribution. This knowledge is particularly relevant in cases of species invasion, where understanding habitat origins can inform management strategies and control measures. In the RS, which represents a native habitat for LF, we identified the highest number of significant microbial markers (with scores exceeding 2.5). This high number of markers may also reflect the influence of the organ factor, which contributed to the microbiota variation in this region, adaptation and fostering niche-specific microbiota associations within different organs ( Kuang et al., 2020 ; Meron et al., 2020 ). In contrast, fish introduced to new habitats, such as in invasive contexts, may only develop distinct organ-specific microbiota at later stages of establishment. The unique and diverse bacterial communities in each organ of RS lionfish suggest the evolution of stable and well-adapted microbiota, shaped over generations by local environmental pressures. In the CA region, a partial trend was observed, potentially reflecting an intermediate stage in establishing lionfish populations. First documented in the Caribbean in the early 2000s, lionfish were officially recognized as a successful invasive species around 2004– 2005 ( Schofield, 2009 ). This indicates that LF in the CA region are beginning to adapt to the local environment. By contrast, the relatively recent invasion of LF in the MS region suggests that adaptation and microbiome stabilization are still in their early stages. In the case of LF from AQ_J, although their origin traces back to the coast of Kenya (a native habitat for LF), the aquarium environment—characterized by artificial water in a closed system—differs significantly from natural conditions. This divergence is also reflected in the composition of the bacterial community. This may explain the observed similarity in microbiota composition across organs and the absence of organ-specific effects in the AQ_J site, a pattern similar to those observed in invaded regions. It is also noteworthy that lionfish invasions in the Atlantic Ocean and Caribbean Sea are largely attributed to releases from aquariums. Such conditions may have shaped the microbiota of released lionfish, resulting in differences from their wild counterparts. While data for direct comparisons are unavailable, this highlights the potential role of aquarium environments in influencing microbiota composition in ways that could facilitate invasion success. Moreover, the environment serves as a microbial reservoir to which organisms are continuously exposed, facilitating the acquisition and exchange of bacteria between the host and its surroundings. Stevens and Olson (2013) demonstrated that eggs extracted from the ovaries of pregnant LF were free of bacteria, indicating the absence of vertical transmission and highlighting the environment’s role as a source of bacterial acquisition for the microbiota composition. The Sex factor was also found to significantly influence the microbiota of LF as reported previously in other marine fish. Most studies have shown sex-based differences in the microbiota of the intestinal organs ( Montalvo & Vázquez, 2018 ; He & Zhang, 2019 ; Hernández & Arrieta, 2020 ). However, microbiota differences have also been documented in other organs, such as the gills ( Wang & Liu, 2021 ), skin ( González-Fernández & Navas, 2020 ) and liver ( Zhao & Xu, 2022 ). These findings, similarly observed in LF, highlight the role of physiological and hormonal differences between sexes in shaping microbial communities, potentially influencing immune responses and adaptation to environmental conditions ( Liu & Yu, 2020 ; Givens & Ransom, 2018 ). Lessepsian invasive fish (MS - RS) Early stages of establishment of the lionfish ( Pterois miles ) in the Mediterranean was reported around 2012 from Lebanon ( Bariche et al., 2013 ), although a single observation from Israel dates back to 1991 ( Golani & Sonin, 1992 ). The second sighting off the Israeli coast occurred in 2013 within the northern Achziv marine reserve, and during 2015 the lionfish has progressed southward to central Israel (Stern & Rothman, 2018). By 2024, P. miles had spread westward to Rhodes and Sicily, successfully establishing populations throughout the eastern Mediterranean ( Bottacini et al., 2024 ). Genetic studies have confirmed that this invasion originated from the Red Sea via the Suez Canal ( Bariche et al., 2017 ; Dimitriou et al., 2019 ). The Suez Canal serves as one of the world’s busiest marine trade routes and provides easy access for non-indigenous species to migrate from the RS to MS through ballast water or as fouling communities on ships ( Williams et al., 1988 ; Galil et al., 2015 ; Lavoie et al., 1999 ; Costello et al., 2022 ). The establishment of lionfish in the Mediterranean follows a well-documented pattern of biological invasion, where populations stabilize in nearby regions before expanding to more distant areas. This pattern is also characteristic of other invasive species, such as Diadema setosum, which have similarly followed the Lessepsian migration route through the Suez Canal ( Zirler et al., 2023 ). During the establishment process, invasive populations grow, adapt to local environmental conditions, and increase their potential for further expansion (Kolar & Lodge, 2000; Zirler et al., 2023 ). We monitored the presence and distribution of LF ( P. miles ) along the Israeli coast beginning in 2017, focusing on two sites, two seasons and three depth zones. Over the five years, we observed, in Israeli coast, a spreading process, north to south and from deeper to shallower waters ( Fig 6A, B ), similar to the findings of Stern & Rothman (2018). The initial invasion of LF along the northern coast of Israel (Achziv), can be attributed by the cyclonic circulation pattern characterizing the Levant Basin in the eastern Mediterranean ( Pinardi et al., 2006 ); however, it is more likely explained by the continuous complex rocky habitats, which may serve as a suitable recruitment area and are characteristic of the initial habitats where LF populations first established in Cyprus, Turkey, and Greece (Stern & Rothman 2018). Notable differences in LF density and distribution were observed between seasons, with higher numbers recorded in the fall. These differences may be driven by adaptations to environmental changes, such as temperature, or by biological factors like predation, food supply, and reproductive cycles. A similar pattern was previously found for other invasive fish which are more abundant in rocky reefs during the fall, when water temperature is higher ( Lazarus et al., 2022 ). These dynamics can lead to shifts in habitat preferences, like increased depth or more sheltered areas (Corey, 2016). Seasonal health status may also play a role in fish distribution and behavior; for example, Minnow fish species exhibited significantly better health status (lower HAI score) and lower parasite abundance and diversity in autumn compared to spring ( Almeida et al., 2024 ). The seasonal trends observed in LF densities and distribution indicate a need for further research to better understand how environmental and biological factors affect their seasonal behavior and habitat preferences. The initial observations of the LF survey were recorded at a depth of around 45 meters before they began appearing in shallower waters (25 and 10m), similar to observations previously reported by Stern & Rothman (2018) . This distribution pattern may suggest that LF initially established themselves in deeper which likely offered a more protected environment with reduced competition, abundant food resources, and lower predation pressure ( jones et al., 2014 ; Tomáš et al.,2018, Chaikin et al., 2022 ). Airey et al., (2023) demonstrated the connection between deep and shallow habitats by examining the post-settlement dispersal patterns of invasive LF. The study found that approximately 34.5% of the LF population in shallow waters consisted of individuals who had migrated from deeper habitats. In comparison, only about 4% had moved in the opposite direction, from shallow to deep waters. This pattern suggests that the primary direction of LF dispersal after settlement is from deeper regions to shallower habitats, emphasizing the role of deep-water areas as potential initial settlement zones before LF spreads into more coastal regions. In addition, their ability to exploit a range of depths likely contributes to their success as an invasive species ( Kletou et al., 2016 ). An increasing trend in the average body size of LF has been observed over the years in surveys conducted in the Mediterranean Sea. The average length (TL) of LF sampled in the RS region (where LF are native) and in CA (an earlier invasion area compared to MS) was approximately 30 cm. In contrast, the average length of LF sampled in the MS region was significantly smaller and only after five years (in 2022) did the LF population in this region reach this size. This observation may indicate the species’ establishment process and thriving in newly invaded regions. Furthermore, studies show that larger body sizes offer several advantages, including improved resources competition, enhanced survival in various environments, and increased reproductive rates. As a result, species with larger body sizes are generally more successful as invaders. ( Roy et al.,2002 ; Schröder et al., 2009 ). Dietary flexibility plays a crucial role in the success of invasive species as it enables them to exploit new resources in diverse environments. Previously studies that compared Lessepsian carnivorous invaders and their counterparts in the RS showed significant differences in isotopic patterns. These differences were more pronounced in S. rubrum , which is an earlier invader, compared to Pterois miles ( Tsadok et al., 2023 ). A similar trend was seen in the calculated TP, where significant differences appeared only in the earlier-established invaders ( Fig. 7c ). In contrast, no difference was seen in TP of LF, which may suggest that the LF are still in the process of establishment and adaptation into the new environment. However, the TP of Caribbean LF was significantly lower from both MS and RS, which again emphasizes the environmental factor. In both S. rubrum (carnivores) and s. rivulatus , (herbivorous) the TP was higher in the RS than in the MS. This reflects how the Mediterranean environment influences food web dynamics and necessitates dietary adjustments by invasive fish during their acclimatization and suggests that Lessepsian species capable of dietary adaptation may be more likely to become successful colonizers ( Golani, 1993 ). It is important to note that fish size was not considered in the TP analyses in this study and it may have affected our results, since fish diet, and particularly LF’s diet shifts as they grow ( Dahl & Patterson 2014 ). Continued study of the LF in the Mediterranean Sea will be able to track the shifts that appear over time in isotope patterns and TP. Photobacterium Photobacterium is widely found in marine environments, where it plays a significant role in ecosystem dynamics, particularly through its interactions with marine fish. This genus, part of the Vibrionaceae family, is notable for its symbiotic and pathogenic interactions with various marine organisms (see review: Labella et al., 2017 ). Symbiotically, certain Photobacterium species, such as P. phosphoreum and P. leiognathi , produce bioluminescence in light organs of fish, contributing to communication, predation, camouflage, and predator protection. Additionally, these species contribute to nutrient cycling by releasing essential fatty acids, which provide nutritional resources within marine food webs and support broader ecological interactions ( Moi 2017 ; Urbanczyk, 2011 , Zarubin, 2012 ). However, much of the literature on Photobacterium species focuses on their potential pathogenicity. This genus is frequently linked to diseases in various marine species, especially fish. For instance, P. damselae , particularly its subspecies P. damselae and P. piscicida , poses significant threats to marine aquaculture and wild fish populations. These bacteria are responsible for causing photobacteriosis, a disease characterized by skin lesions, septicemia, and high mortality rates among affected fish, such as groupers and sea bream. ( Rivas et al., 2013 ; Labella et al., 2017 , Austin & Austin, 2016 ). In our study, Photobacterium was found in all samples and constituted a dominant component of the LF microbiome core. However, compared to the RS region, the relative abundance of Photobacterium was significantly higher in the MS, where it exceeded 50% across all examined organs. This high abundance may explain the lower number of taxa (Genus level) observed in the MS compared to the RS region, in all examined organs ( Fig 4 ). Interestingly, further comparative analysis of other Lessepsian invasive fish species, S. rubrum S. rivulatus, revealed a similar trend, with a significantly higher relative abundance of Photobacterium in the fish organs from the MS compared to the RS region. Invasive species may acquire beneficial microbes from co-occurring native species, which can accelerate their adaptation to local conditions and reduce invasion lag times ( Martignoni & Kolodny, 2023 ). Two dominant sequences were identified among the Photobacterium found in the MS and RS samples, corresponding to Photobacterium damselae subsp . damselae , according to BLAST results from NCBI. Interestingly, despite identifying the pathogenic species, no external or internal signs of disease were observed in the fish analyzed. This aligns with similar observations in additional studies on wild Mediterranean fish and sharks, where Photobacterium species presence did not correlate with visible disease symptoms ( Meron et al., 2020 ; Bregman et al, 2023)). The high relative abundance of Photobacterium in Mediterranean lionfish (LF) raises further questions about its potential role in the fish’s resistance and adaptation. For instance, are lionfish resistant to these pathogens, or do these bacteria contribute to their acclimatization to new environmental conditions, thereby enhancing their success as invasive species? These questions underscore the need for further research to explore the dynamics and ecological mechanisms underlying the interactions between lionfish and Photobacterium. Summary This study focuses on the microbiota of the lionfish (LF), characterizing its species-specific core microbiome and identifying habitat-specific markers while discussing its potential role in invasion and establishment. The findings suggest that lionfish (LF) are still establishing themselves in the Mediterranean environment but are exhibiting a growing presence. The relatively recent invasion of LF from the Red Sea to the Mediterranean offers a unique opportunity to monitor microbial and ecological establishment processes, including population distribution, body size, and dietary patterns. This understanding can offer valuable insights into ecological impacts and support the development of effective management strategies for invasive species. Funding The research leading to these results received partial funding from the European Union’s Horizon 2020 research and innovation program under grant agreement No 10730, ASSEMBLE Plus project. Acknowledgements We would like to thank the staff at Caribbean Netherlands Science Institute (CNSI), The Gottesman Family Israel Aquarium, Jerusalem and Morris Kahn Marine Research Station (MKMRS) for their assistance and support. References 1. ↵ Adams , A. S. , Jordan , M. S. , Adams , S. M. , Suen , G. , Goodwin , L. A. , Davenport , K. W. , … & Raffa , K. F . ( 2011 ). Cellulose-degrading bacteria associated with the invasive woodwasp Sirex noctilio . The ISME journal , 5 ( 8 ), 1323 – 1331 . OpenUrl CrossRef PubMed 2. ↵ Aires , T. , Serrão , E. A. , & Engelen , A. H . ( 2016 ). Host and environmental specificity in bacterial communities associated to two highly invasive marine species (genus Asparagopsis ) . Frontiers in microbiology , 7 , 559 . OpenUrl PubMed 3. ↵ Aires , T. , Serrão , E. A. , Kendrick , G. , Duarte , C. M. , & Arnaud-Haond , S . ( 2013 ). Invasion is a community affair: clandestine followers in the bacterial community associated to green algae, Caulerpa racemosa , track the invasion source . PLoS One , 8 ( 7 ), e68429 . OpenUrl CrossRef PubMed 4. ↵ Airey , M. E. , Fogg , A. Q. , & Drew , J. A . ( 2023 ). Invasive lionfish dispersal between shallow-and deep-water habitats within coastal Floridian waters . Biological Invasions , 25 ( 12 ), 3983 – 3991 . OpenUrl CrossRef 5. ↵ Albins , M. A. , & Hixon , M. A . ( 2008 ). Invasive Indo-Pacific lionfish Pterois volitans reduce recruitment of Atlantic coral-reef fishes . Marine Ecology Progress Series , 367 , 233 – 238 . OpenUrl CrossRef 6. Albins , M. A. , & Hixon , M. A . ( 2013 ). Worst case scenario: potential long-term effects of invasive predatory lionfish ( Pterois volitans ) on Atlantic and Caribbean coral-reef communities . Environmental Biology of Fishes , 96 , 1151 – 1157 . OpenUrl CrossRef 7. ↵ Almeida , D. , Alcaraz-Hernández , J. D. , Cruz , A. , Lantero , E. , Fletcher , D. H. , & García-Berthou , E . ( 2024 ). Seasonal effects on health status and parasitological traits of an invasive minnow in Iberian waters . Animals , 14 ( 10 ), 1502 . OpenUrl CrossRef PubMed 8. ↵ Amalfitano , S. , Coci , M. , Corno , G. , & Luna , G. M . ( 2015 ). A microbial perspective on biological invasions in aquatic ecosystems . Hydrobiologia , 746 , 13 – 22 . OpenUrl CrossRef 9. ↵ Apprill , A. , McNally , S. , Parsons , R. , & Weber , L . ( 2015 ). Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton . Aquatic Microbial Ecology , 75 ( 2 ), 129 – 137 . OpenUrl CrossRef PubMed 10. ↵ Arredondo-Chavez , A. T. , Sanchez-Jimenez , J. A. , Avila-Morales , O. G. , Torres-Chavez , P. , Herrerias-Diego , Y. , Medina-Nava , M. , … & Caballero-Vazquez , J. A . ( 2016 ). Spatio-temporal variation in the diet composition of red lionfish, Pterois volitans (Actinopterygii: Scorpaeniformes: Scorpaenidae ), in the Mexican Caribbean: insights into the ecological effect of the alien invasion . Acta Ichthyologica et Piscatoria , 46 ( 3 ). 11. ↵ Austin , B. , & Austin , D. A . ( 2016 ). Bacterial fish pathogens: disease of farmed and wild fish . Springer . 12. Azzurro , E. , Stancanelli , B. , Di Martino , V. , & Bariche , M. ( 2017 ). Range expansion of the common lionfish Pterois miles (Bennett, 1828) in the Mediterranean Sea: an unwanted new guest for Italian waters . BioInvasions Records , 6 ( 2 ), 95 – 98 . OpenUrl CrossRef 13. ↵ Bagra , K. , Bellanger , X. , Merlin , C. , Singh , G. , Berendonk , T. U. , & Klümper , U . ( 2022 ). Exposure to environmental stress decreases the resistance of river microbial communities towards invasion by antimicrobial resistant bacteria . BioRxiv , 2022 – 11 . 14. ↵ Bang , C. , Dagan , T. , Deines , P. , Dubilier , N. , Duschl , W. J. , Fraune , S. , … & Bosch , T. C . ( 2018 ). Metaorganisms in extreme environments: do microbes play a role in organismal adaptation? . Zoology , 127 , 1 – 19 . OpenUrl CrossRef PubMed 15. ↵ Bariche , M. , Kleitou , P. , Kalogirou , S. , & Bernardi , G . ( 2017 ). Genetics reveal the identity and origin of the lionfish invasion in the Mediterranean Sea . Scientific reports , 7 ( 1 ), 6782 . OpenUrl CrossRef PubMed 16. ↵ Bariche , M. , Torres , M. , & Azzurro , E. ( 2013 ). The presence of the invasive Lionfish Pterois miles in the Mediterranean Sea . Mediterranean Marine Science , 292 – 294 . 17. Ben-Tuvia , A . ( 1985 ). The impact of the Lessepsian (Suez Canal) fish migration on the eastern Mediterranean ecosystem . In Mediterranean marine ecosystems (pp. 367 – 375 ). Boston, MA : Springer US . 18. ↵ Bottacini , D. , Pollux , B. J. , Nijland , R. , Jansen , P. A. , Naguib , M. , & Kotrschal , A . ( 2024 ). Lionfish ( Pterois miles ) in the Mediterranean Sea: a review of the available knowledge with an update on the invasion front . NeoBiota , 92 , 233 – 257 . OpenUrl CrossRef 19. ↵ Brown , A. M. , Huynh , L. Y. , Bolender , C. M. , Nelson , K. G. , & McCutcheon , J. P . ( 2014 ). Population genomics of a symbiont in the early stages of a pest invasion . Molecular Ecology , 23 ( 6 ), 1516 – 1530 . OpenUrl CrossRef Web of Science 20. ↵ Cavicchioli , R. , Ripple , W. J. , Timmis , K. N. , Azam , F. , Bakken , L. R. , Baylis , M. , … & Webster , N. S . ( 2019 ). Scientists’ warning to humanity: microorganisms and climate change . Nature Reviews Microbiology , 17 ( 9 ), 569 – 586 . OpenUrl CrossRef PubMed 21. ↵ Chaikin , S. , Dubiner , S. , & Belmaker , J . ( 2022 ). Cold-water species deepen to escape warm water temperatures . Global Ecology and Biogeography , 31 ( 1 ), 75 – 88 . OpenUrl CrossRef 22. Cheng , D. , Chen , S. , Huang , Y. , Pierce , N. E. , Riegler , M. , Yang , F. , … & Xu , Y. ( 2019 ). Symbiotic microbiota may reflect host adaptation by resident to invasive ant species . PLoS pathogens , 15 ( 7 ), e1007942 . OpenUrl CrossRef PubMed 23. ↵ Coats , V. C. , & Rumpho , M. E . ( 2014 ). The rhizosphere microbiota of plant invaders: an overview of recent advances in the microbiomics of invasive plants . Frontiers in microbiology , 5 , 368 . OpenUrl PubMed 24. ↵ Costello , K. E. , Lynch , S. A. , McAllen , R. , O’Riordan , R. M. , & Culloty , S. C . ( 2022 ). Assessing the potential for invasive species introductions and secondary spread using vessel movements in maritime ports . Marine Pollution Bulletin , 177 , 113496 . OpenUrl CrossRef PubMed 25. Côté , I. M. , & Maljković , A . ( 2010 ). Predation rates of Indo-Pacific lionfish on Bahamian coral reefs . Marine Ecology Progress Series , 404 , 219 – 225 . OpenUrl CrossRef 26. ↵ Côté , I. M. , & Smith , N. S . ( 2018 ). The lionfish Pterois sp. invasion: Has the worst-case scenario come to pass? . Journal of Fish Biology , 92 ( 3 ), 660 – 689 . OpenUrl CrossRef PubMed 27. ↵ Côté , I. M. , Green , S. J. , & Hixon , M. A . ( 2013 ). Predatory fish invaders: Insights from Indo-Pacific lionfish in the western Atlantic and Caribbean . Biological Conservation , 164 , 50 – 61 . OpenUrl CrossRef 28. ↵ Coutinho , F. H. , Meirelles , P. M. , Moreira , A. P. B. , Paranhos , R. P. , Dutilh , B. E. , & Thompson , F. L . ( 2015 ). Niche distribution and influence of environmental parameters in marine microbial communities: a systematic review . PeerJ , 3 , e1008 . OpenUrl CrossRef PubMed 29. ↵ Dahl , K. A. , & Patterson III , W. F . ( 2014 ). Habitat-specific density and diet of rapidly expanding invasive red lionfish, Pterois volitans, populations in the northern Gulf of Mexico . PloS one , 9 ( 8 ), e105852 . OpenUrl CrossRef PubMed 30. ↵ Dimitriou , A. C. , Chartosia , N. , Hall-Spencer , J. M. , Kleitou , P. , Jimenez , C. , Antoniou , C. , … & Sfenthourakis , S. ( 2019 ). Genetic data suggest multiple introductions of the lionfish ( Pterois miles ) into the Mediterranean Sea . Diversity , 11 ( 9 ), 149 . OpenUrl CrossRef 31. Eddy , C. , Pitt , J. , Morris Jr , J. A. , Smith , S. , Goodbody-Gringley , G. , & Bernal , D . ( 2016 ). Diet of invasive lionfish ( Pterois volitans and P. miles ) in Bermuda . Marine Ecology Progress Series , 558 , 193 – 206 . OpenUrl CrossRef 32. ↵ Ferreira , C. E. , Luiz , O. J. , Floeter , S. R. , Lucena , M. B. , Barbosa , M. C. , Rocha , C. R. , & Rocha , L. A . ( 2015 ). First record of invasive lionfish ( Pterois volitans ) for the Brazilian coast . PloS one , 10 ( 4 ), e0123002 . OpenUrl CrossRef PubMed 33. ↵ Galil , B. S. , Boero , F. , Campbell , M. L. , Carlton , J. T. , Cook , E. , Fraschetti , S. , … & Ruiz , G. M . ( 2015 ). ‘Double trouble’: the expansion of the Suez Canal and marine bioinvasions in the Mediterranean Sea . Biological Invasions , 17 , 973 – 976 . OpenUrl CrossRef 34. ↵ Givens , C. E. , & Ransom , B . ( 2018 ). Microbiome differences in the reproductive stages of female and male fish . Fish Physiology and Biochemistry , 44 ( 2 ), 409 – 422 . OpenUrl 35. ↵ Golani , D . ( 1993 ). Trophic adaptation of Red Sea fishes to the eastern Mediterranean environment - a review and new data . Israel Journal of Zoology , 39 ( 4 ), 391 – 402 . OpenUrl Web of Science 36. ↵ Golani , D. , & Sonin , O . ( 1992 ). New records of the Red Sea fishes, Pterois miles (Scorpaenidae) and Pteragogus pelycus (Labridae) from the eastern Mediterranean Sea . Japanese Journal of Ichthyology , 39 ( 2 ), 167 – 169 . OpenUrl 37. ↵ González-Fernández , C. , & Navas , J. M . ( 2020 ). Skin microbiome of fish: A review of its role in health and disease . Journal of Fish Diseases , 43 ( 4 ), 431 – 441 . OpenUrl CrossRef PubMed 38. ↵ Green , S. J. , Akins , J. L. , Maljković , A. , & Côté , I. M . ( 2012 ). Invasive lionfish drive Atlantic coral reef fish declines . PloS one , 7 ( 3 ), e32596 . OpenUrl CrossRef PubMed 39. ↵ Haas G , Steinitz H ( 1947 ) Erythrean fishes on the Mediterranean coast of Palestine . Nature 160 ( 4053 ): 28 OpenUrl 40. Havel , J. E. , Kovalenko , K. E. , Thomaz , S. M. , Amalfitano , S. , & Kats , L. B . ( 2015 ). Aquatic invasive species: challenges for the future . Hydrobiologia , 750 , 147 – 170 . OpenUrl CrossRef PubMed 41. ↵ He , Y. , & Zhang , Q . ( 2019 ). Differences in gut microbiota between male and female Japanese sea bass ( Lateolabrax japonicus ) . Aquaculture Research , 50 ( 7 ), 1958 – 1968 . OpenUrl 42. ↵ Hernández , P. E. , & Arrieta , C . ( 2020 ). Microbiome comparison between male and female Seriola lalandi (California yellowtail) in a hatchery environment . Aquaculture Research , 51 ( 6 ), 2214 – 2223 . OpenUrl 43. ↵ Hicks , N. , Liu , X. , Gregory , R. , Kenny , J. , Lucaci , A. , Lenzi , L. , … & Duncan , K. R . ( 2018 ). Temperature driven changes in benthic bacterial diversity influences biogeochemical cycling in coastal sediments . Frontiers in microbiology , 9 , 1730 . OpenUrl CrossRef PubMed 44. ↵ Himler , A. G. , Adachi-Hagimori , T. , Bergen , J. E. , Kozuch , A. , Kelly , S. E. , Tabashnik , B. E. , … & Hunter , M. S . ( 2011 ). Rapid spread of a bacterial symbiont in an invasive whitefly is driven by fitness benefits and female bias . science , 332 ( 6026 ), 254 – 256 . OpenUrl Abstract / FREE Full Text 45. Seckbach , J. , Chapman , D Hoffman , R. , Dubinsky , Z . ( 2010 ). Invasive and Alien Rhodophyta in the Mediterranean and along the Israeli Shores . In: Seckbach , J. , Chapman , D . (eds) Red Algae in the Genomic Age. Cellular Origin, Life in Extreme Habitats and Astrobiology, vol 13 . Springer , Dordrecht . 46. ↵ Huang , J. , Zhang , Y. , & Yan , S . ( 2019 ). Microbial plasticity and the invasion success of non-native species . Microbiome , 7 ( 1 ), 34 . OpenUrl CrossRef PubMed 47. ↵ Jones , L. A. , & Ricciardi , A . ( 2014 ). The influence of pre-settlement and early post-settlement processes on the adult distribution and relative dominance of two invasive mussel species . Freshwater Biology , 59 ( 5 ), 1086 – 1100 . OpenUrl CrossRef 48. Jůza , T. , Blabolil , P. , Baran , R. , Draštík , V. , Holubová , M. , Kočvara , L. , … & Ketelaars , H. A . ( 2018 ). Comparison of two passive methods for sampling invasive round goby (Neogobius melanostomus) populations at different depths in artificial lakes . Fisheries Research , 207 , 175 – 181 . OpenUrl CrossRef 49. Katsanevakis , S. , Coll , M. , Piroddi , C. , Steenbeek , J. , Ben Rais Lasram , F. , Zenetos , A. , & Cardoso , A. C. ( 2014 ). Invading the Mediterranean Sea: biodiversity patterns shaped by human activities . Frontiers in Marine Science , 1 , 32 . OpenUrl 50. ↵ Kletou , D. , Hall-Spencer , J. M. , & Kleitou , P . ( 2016 ). A lionfish (Pterois miles) invasion has begun in the Mediterranean Sea . Marine Biodiversity Records , 9 , 1 – 7 . OpenUrl CrossRef 51. Kolar , C. S. , & Lodge , D. M . ( 2001 ). Progress in invasion biology: predicting invaders . Trends in ecology & evolution , 16 ( 4 ), 199 – 204 . OpenUrl CrossRef PubMed 52. ↵ Kuang , T. , He , A. , Lin , Y. , Huang , X. , Liu , L. , & Zhou , L . ( 2020 ). Comparative analysis of microbial communities associated with the gill, gut, and habitat of two filter-feeding fish . Aquaculture Reports , 18 , 100501 OpenUrl CrossRef 53. ↵ Labella , A. M. , Arahal , D. R. , Castro , D. , Lemos , M. L. , & Borrego , J. J . ( 2017 ). Revisiting the genus Photobacterium: taxonomy, ecology and pathogenesis . Int Microbiol , 20 ( 1 ), 1 – 10 . OpenUrl CrossRef PubMed 54. ↵ Lavoie , D. M. , Smith , L. D. , & Ruiz , G. M . ( 1999 ). The potential for intracoastal transfer of non-indigenous species in the ballast water of ships . Estuarine, Coastal and Shelf Science , 48 ( 5 ), 551 – 564 . OpenUrl CrossRef Web of Science 55. ↵ Lazarus M , Frid-Landau O , Yahel R ( 2022 ) Marine BioBlitz-biological survey of the Israelis MPAs in the Mediterranean Sea-Summary report for 2015, 2017, 2019, 2021 - Chapter B- Composition of fish community . Jerusalem, Israel Nature and Parks Authority (in Hebrew) , p. 90 56. ↵ Lee , R. W. , Robinson , J. J. , & Cavanaugh , C. M . ( 1999 ). Pathways of inorganic nitrogen assimilation in chemoautotrophic bacteria–marine invertebrate symbioses: expression of host and symbiont glutamine synthetase . Journal of Experimental Biology , 202 ( 3 ), 289 – 300 . OpenUrl Abstract 57. ↵ Liu , M. , Bjørnlund , L. , Rønn , R. , Christensen , S. , & Ekelund , F . ( 2012 ). Disturbance promotes non-indigenous bacterial invasion in soil microcosms: analysis of the roles of resource availability and community structure . PLOS ONE , 7 ( 10 ), e45306 . OpenUrl CrossRef PubMed 58. ↵ Liu , W. , & Yu , J . ( 2020 ). Differences in gut microbiota between male and female zebrafish ( Danio rerio ) and their association with stress responses . Fish & Shellfish Immunology , 97 , 245 – 251 . OpenUrl 59. ↵ Lu , M. , Hulcr , J. , & Sun , J . ( 2016 ). The role of symbiotic microbes in insect invasions . Annual Review of Ecology, Evolution, and Systematics , 47 ( 1 ), 487 – 505 . OpenUrl CrossRef 60. ↵ Martignoni , M. M. , & Kolodny , O . ( 2023 ). Microbiome transfer from native to invasive species may increase invasion risk and shorten invasion lag . bioRxiv , 2023 – 08 . 61. ↵ Meron , D. , Davidovich , N. , Ofek-Lalzar , M. , Berzak , R. , Scheinin , A. , Regev , Y. , … & Morick , D. ( 2020 ). Specific pathogens and microbial abundance within liver and kidney tissues of wild marine fish from the Eastern Mediterranean Sea . Microbial Biotechnology , 13 ( 3 ), 770 – 780 . OpenUrl CrossRef PubMed 62. ↵ Moi , I. M. , Roslan , N. N. , Leow , A. T. C. , Ali , M. S. M. , Rahman , R. N. Z. R. A. , Rahimpour , A. , & Sabri , S . ( 2017 ). The biology and the importance of Photobacterium species . Applied microbiology and biotechnology , 101 ( 11 ), 4371 – 4385 . OpenUrl CrossRef 63. ↵ Montalvo , T. , & Vázquez , M . ( 2018 ). Sex-specific differences in the gut microbiota of the marine fish Mugil cephalus (Linnaeus) . Microbial Ecology , 75 ( 4 ), 1015 – 1028 . OpenUrl 64. ↵ Moonsamy , P. V. , Williams , T. , Bonella , P. , Holcomb , C. L. , Höglund , B. N. , Hillman , G. , … & Erlich , H. A . ( 2013 ). High throughput HLA genotyping using 454 sequencing and the Fluidigm Access Array™ system for simplified amplicon library preparation . Tissue antigens , 81 ( 3 ), 141 – 149 . OpenUrl CrossRef PubMed Web of Science 65. ↵ Morris Jr , J. A. , Thomas , A. , Rhyne , A. L. , Breen , N. , Akins , L. , & Nash , B. ( 2011 ). Nutritional properties of the invasive lionfish: a delicious and nutritious approach for controlling the invasion . 66. ↵ N. Raghavachari and N. Garcia-Reyero Naqib , A. , Poggi , S. , Wang , W. , Hyde , M. , Kunstman , K. , and Green , S. J . ( 2018 ). Making and sequencing heavily multiplexed, high-throughput 16S ribosomal RNA gene amplicon libraries using a flexible, two-stage PCR protocol, in Gene Expression Analysis: Methods and Protocols . eds. N. Raghavachari and N. Garcia-Reyero ( New York, NY : Springer New York ), 149 – 169 . 67. ↵ Paerl , H. W. , Gardner , W. S. , Havens , K. E. , Joyner , A. R. , McCarthy , M. J. , Newell , S. E. , … & Scott , J. T . ( 2016 ). Mitigating cyanobacterial harmful algal blooms in aquatic ecosystems impacted by climate change and anthropogenic nutrients . Harmful Algae , 54 , 213 – 222 . OpenUrl CrossRef PubMed 68. ↵ Parada , A. E. , Needham , D. M. , & Fuhrman , J. A . ( 2016 ). Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples . Environmental microbiology , 18 ( 5 ), 1403 – 1414 . OpenUrl CrossRef PubMed 69. ↵ Pinardi , N. , et al. ( 2006 ). The Mediterranean Sea circulation and the impact of climate variability . Deep Sea Research Part II: Topical Studies in Oceanography , 53 ( 3-4 ), 417 – 432 . OpenUrl 70. ↵ Rivas , A. J. , Lemos , M. L. , & Osorio , C. R . ( 2013 ). Photobacterium damselae subsp. damselae , a bacterium pathogenic for marine animals and humans . Frontiers in microbiology , 4 , 283 . OpenUrl PubMed 71. ↵ Ronanki , J . ( 2024 ). Metabolic Adaptations in the Gut Microbiome of Invasive Marine Fish Species . NATURALISTA CAMPANO , 28 ( 1 ), 3108 – 3118 . OpenUrl 72. ↵ Roy , K. , Jablonski , D. , & Valentine , J. W . ( 2002 ). Body size and invasion success in marine bivalves . Ecology Letters , 5 ( 2 ), 163 – 167 . OpenUrl CrossRef Web of Science 73. Sampson , T. R. , & Mazmanian , S. K . ( 2015 ). Control of brain development, function, and behavior by the microbiome . Cell host & microbe , 17 ( 5 ), 565 – 576 . OpenUrl CrossRef PubMed 74. ↵ Schofield , P. J . ( 2009 ). Geographic extent and chronology of the invasion of non-native lionfish ( Pterois volitans and P. miles ) in the Western North Atlantic and Caribbean Sea . Aquatic Invasions , 4 ( 3 ), 473 – 479 . OpenUrl CrossRef 75. ↵ Schröder , A. , Nilsson , K. A. , Persson , L. , Van Kooten , T. , & Reichstein , B. ( 2009 ). Invasion success depends on invader body size in a size-structured mixed predation– competition community . Journal of Animal Ecology , 78 ( 6 ), 1152 – 1162 . OpenUrl CrossRef PubMed 76. ↵ Steinitz , W ., Beiträge zur Kenntnis der Küstenfauna Palästinas, Erster Teil , Pubbl. Stn. Zool. Napoli , 1927 , vol. 8 , pp. 311 – 353 . OpenUrl 77. ↵ Stern , N. , Jimenez , C. , Huseyinoglu , M. F. , Andreou , V. , Hadjioannou , L. , Petrou , A. , … & Rothman , S. B . ( 2019 ). Constructing the genetic population demography of the invasive lionfish Pterois miles in the Levant Basin , Eastern Mediterranean. Mitochondrial Dna Part A , 30 ( 2 ), 249 – 255 . OpenUrl CrossRef 78. ↵ Stern , N. , Rothman , S. B. , Hüseyinoglu , M. F. , & Öztürk , B . ( 2018 ). Iron Lion Zion: the successful, albeit lingered, invasion of the lionfish in the Israeli Mediterranean Sea. Lionfish Invasion and Its Management in the Mediterranean Sea. Turkish Marine Research Foundation (TUDAV) Publication Istanbul , Turkey , 49 , 51 – 56 . OpenUrl 79. ↵ Stevens , J. L. , & Olson , J. B . ( 2013 ). Invasive lionfish harbor a different external bacterial community than native Bahamian fishes . Coral Reefs , 32 , 1113 – 1121 . OpenUrl CrossRef 80. ↵ Stevens , J. L. , & Olson , J. B . ( 2015 ). Bacterial communities associated with lionfish in their native and invaded ranges . Marine Ecology Progress Series , 531 , 253 – 262 . OpenUrl CrossRef 81. ↵ Stevens , J. L. , Jackson , R. L. , & Olson , J. B . ( 2016 ). Bacteria associated with lionfish ( Pterois volitans/miles complex) exhibit antibacterial activity against known fish pathogens . Marine Ecology Progress Series , 558 , 167 – 180 . OpenUrl CrossRef 82. ↵ Takizawa , Y. , & Chikaraishi , Y . ( 2024 ). Trophic isotopic discrimination of amino acids (glutamic acid and phenylalanine) in the leopard gecko Eublepharis macularius : a report based on a controlled feeding experiment . Researches in Organic Geochemistry , 40 ( 1 ), 1 – 7 . OpenUrl 83. ↵ Thompson , L. R. , Williams , G. J. , Haroon , M. F. , Shibl , A. , Larsen , P. , Shorenstein , J. , … & Stingl , U. ( 2017 ). Metagenomic covariation along densely sampled environmental gradients in the Red Sea . The ISME journal , 11 ( 1 ), 138 – 151 . OpenUrl CrossRef PubMed 84. ↵ Tsadok , R. , Zemah-Shamir , Z. , Shemesh , E. , Martinez , S. , Ramon , D. , Kolski , I. , … & Tchernov , D. ( 2023 ). Dietary habits change of Lessepsian migrants’ fish from the Red Sea to the Eastern Mediterranean Sea . Aquatic Invasions , 18 ( 4 ), 521 – 531 . OpenUrl CrossRef 85. ↵ Urbanczyk , H. , Ast , J. C. , & Dunlap , P. V . ( 2011 ). Phylogeny, genomics, and symbiosis of Photobacterium . FEMS microbiology reviews , 35 ( 2 ), 324 – 342 . OpenUrl CrossRef PubMed 86. ↵ Vilcinskas , A. , Stoecker , K. , Schmidtberg , H. , Röhrich , C. R. , & Vogel , H . ( 2013 ). Invasive harlequin ladybird carries biological weapons against native competitors . Science , 340 ( 6134 ), 862 – 863 . OpenUrl Abstract / FREE Full Text 87. ↵ Wang , Y. , & Liu , Y . ( 2021 ). Gill microbiota diversity in marine fish is affected by sex and season . Fish & Shellfish Immunology , 113 , 174 – 183 . OpenUrl 88. ↵ Williams , R. J. , Griffiths , F. B. , Van der Wal , E. J. , & Kelly , J. ( 1988 ). Cargo vessel ballast water as a vector for the transport of non-indigenous marine species . Estuarine, Coastal and Shelf Science , 26 ( 4 ), 409 – 420 . OpenUrl CrossRef Web of Science 89. ↵ Zarubin , M. , Belkin , S. , Ionescu , M. , & Genin , A . ( 2012 ). Bacterial bioluminescence as a lure for marine zooplankton and fish . Proceedings of the National Academy of Sciences , 109 ( 3 ), 853 – 857 . OpenUrl Abstract / FREE Full Text 90. ↵ Zhao , Y. , & Xu , Y . ( 2022 ). Differences in liver microbiome between male and female Japanese sea bass ( Lateolabrax japonicus ) under stress . Aquaculture Research , 53 ( 2 ), 543 – 556 . OpenUrl 91. ↵ Zirler , R. , Leck , L. A. , Farkash , T. F. , Holzknecht , M. , Kroh , A. , Gerovasileiou , V. , … & Bronstein , O. ( 2023 ). Gaining a (tube) foothold–trends and status following two decades of the long-spined echinoid Diadema setosum (Leske, 1778) invasion to the Mediterranean Sea . Frontiers in Marine Science , 10 , 1152584 . OpenUrl CrossRef View the discussion thread. Back to top Previous Next Posted February 07, 2025. Download PDF Supplementary Material Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Microbiota Dynamics in Lionfish (Pterois): Insights into Invasion and Establishment in the Mediterranean Sea Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Microbiota Dynamics in Lionfish (Pterois): Insights into Invasion and Establishment in the Mediterranean Sea Dalit Meron , Maya Lalzar , Shevy Bat-Sheva Rothman , Yael Kroin , Elizabeth Kaufman , Kimani Kitson-Walters , Tal Zvi-Kedem , Eli Shemesh , Rami Tsadok , Hagai Nativ , Shai Einbinder , Dan Tchernov bioRxiv 2025.02.03.635474; doi: https://doi.org/10.1101/2025.02.03.635474 Share This Article: Copy Citation Tools Microbiota Dynamics in Lionfish (Pterois): Insights into Invasion and Establishment in the Mediterranean Sea Dalit Meron , Maya Lalzar , Shevy Bat-Sheva Rothman , Yael Kroin , Elizabeth Kaufman , Kimani Kitson-Walters , Tal Zvi-Kedem , Eli Shemesh , Rami Tsadok , Hagai Nativ , Shai Einbinder , Dan Tchernov bioRxiv 2025.02.03.635474; doi: https://doi.org/10.1101/2025.02.03.635474 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Ecology Subject Areas All Articles Animal Behavior and Cognition (7622) Biochemistry (17648) Bioengineering (13870) Bioinformatics (41880) Biophysics (21423) Cancer Biology (18553) Cell Biology (25458) Clinical Trials (138) Developmental Biology (13364) Ecology (19866) Epidemiology (2067) Evolutionary Biology (24290) Genetics (15589) Genomics (22475) Immunology (17711) Microbiology (40327) Molecular Biology (17145) Neuroscience (88472) Paleontology (666) Pathology (2826) Pharmacology and Toxicology (4815) Physiology (7635) Plant Biology (15114) Scientific Communication and Education (2044) Synthetic Biology (4286) Systems Biology (9815) Zoology (2268)