Heteroplasmy is rare in plant mitochondria compared to plastids despite similar mutation rates.

Plant cells harbor two membrane-bound organelles containing their own genetic material -plastids and mitochondria. Although the two organelles co-exist and co-evolve within the same plant cells, they differ in genome copy number, intracellular organization, and mode of segregation. How these attributes affect the time to fixation, or conversely, loss of neutral alleles is currently unresolved. Here we show that mitochondria and plastids share the same mutation rate yet plastid alleles remain in a heteroplasmic state significantly longer compared to mitochondrial alleles. By analyzing genetic variants across populations of the marine flowering plant Zostera marina and simulating organelle allele dynamics, we examine the determinants of allele segregation and allele fixation. Our results suggest that bottlenecks on the cell population, e.g., during branching or seeding, and stratification of the meristematic tissue, are important determinants of mitochondrial allele dynamics. Furthermore, we suggest that the prolonged plastid allele dynamics are due to a yet unknown active plastid partition mechanism. The dissimilarity between plastid and mitochondrial novel allele fixation at different levels of organization may manifest in differences in adaptation processes. Our study uncovers fundamental principles of organelle population genetics that are essential for further investigations of long-term evolution and molecular dating of divergence events.

A somatic genetic clock for clonal species.

Age and longevity are key parameters for demography and life-history evolution of organisms. In clonal species, a widespread life history among animals, plants, macroalgae and fungi, the sexually produced offspring (genet) grows indeterminately by producing iterative modules, or ramets, and so obscure their age. Here we present a novel molecular clock based on the accumulation of fixed somatic genetic variation that segregates among ramets. Using a stochastic model, we demonstrate that the accumulation of fixed somatic genetic variation will approach linearity after a lag phase, and is determined by the mitotic mutation rate, without direct dependence on asexual generation time. The lag phase decreased with lower stem cell population size, number of founder cells for the formation of new modules, and the ratio of symmetric versus asymmetric cell divisions. We calibrated the somatic genetic clock on cultivated eelgrass Zostera marina genets (4 and 17 years respectively). In a global data set of 20 eelgrass populations, genet ages were up to 1,403 years. The somatic genetic clock is applicable to any multicellular clonal species where the number of founder cells is small, opening novel research avenues to study longevity and, hence, demography and population dynamics of clonal species.

The archaeome in metaorganism research, with a focus on marine models and their bacteria-archaea interactions.

Metaorganism research contributes substantially to our understanding of the interaction between microbes and their hosts, as well as their co-evolution. Most research is currently focused on the bacterial community, while archaea often remain at the sidelines of metaorganism-related research. Here, we describe the archaeome of a total of eleven classical and emerging multicellular model organisms across the phylogenetic tree of life. To determine the microbial community composition of each host, we utilized a combination of archaea and bacteria-specific 16S rRNA gene amplicons. Members of the two prokaryotic domains were described regarding their community composition, diversity, and richness in each multicellular host. Moreover, association with specific hosts and possible interaction partners between the bacterial and archaeal communities were determined for the marine models. Our data show that the archaeome in marine hosts predominantly consists of Nitrosopumilaceae and Nanoarchaeota, which represent keystone taxa among the porifera. The presence of an archaeome in the terrestrial hosts varies substantially. With respect to abundant archaeal taxa, they harbor a higher proportion of methanoarchaea over the aquatic environment. We find that the archaeal community is much less diverse than its bacterial counterpart. Archaeal amplicon sequence variants are usually host-specific, suggesting adaptation through co-evolution with the host. While bacterial richness was higher in the aquatic than the terrestrial hosts, a significant difference in diversity and richness between these groups could not be observed in the archaeal dataset. Our data show a large proportion of unclassifiable archaeal taxa, highlighting the need for improved cultivation efforts and expanded databases.

Seagrass genomes reveal ancient polyploidy and adaptations to the marine environment.

We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the 'savannahs of the sea' are of major concern in times of climate change and loss of biodiversity.

Worldwide Population Genomics Reveal Long-Term Stability of the Mitochondrial Genome Architecture in a Keystone Marine Plant.

Mitochondrial genomes (mitogenomes) of flowering plants are composed of multiple chromosomes. Recombination within and between the mitochondrial chromosomes may generate diverse DNA molecules termed isoforms. The isoform copy number and composition can be dynamic within and among individual plants due to uneven replication and homologous recombination. Nonetheless, despite their functional importance, the level of mitogenome conservation within species remains understudied. Whether the ontogenetic variation translates to evolution of mitogenome composition over generations is currently unknown. Here we show that the mitogenome composition of the seagrass Zostera marina is conserved among worldwide populations that diverged ca. 350,000 years ago. Using long-read sequencing, we characterized the Z. marina mitochondrial genome and inferred the repertoire of recombination-induced configurations. To characterize the mitochondrial genome architecture worldwide and study its evolution, we examined the mitogenome in Z. marina meristematic region sampled in 16 populations from the Pacific and Atlantic oceans. Our results reveal a striking similarity in the isoform relative copy number, indicating a high conservation of the mitogenome composition among distantly related populations and within the plant germline, despite a notable variability during individual ontogenesis. Our study supplies a link between observations of dynamic mitogenomes at the level of plant individuals and long-term mitochondrial evolution.

Population genetic differentiation of the ubiquitous brooding coral Pocillopora acuta along Phuket Island reefs in the Andaman Sea, Thailand.

BACKGROUND: The widespread Indo-Pacific coral species Pocillopora acuta Lamarck, 1816 displays varying levels of asexual versus sexual reproduction, with strong repercussions on genetic diversity, connectivity and genetic structuring within and among populations. For many geographic regions, baseline information on genetic diversity is still lacking, particularly in the Andaman Sea. The region suffered a massive heat-induced bleaching event in 2010 with high coral cover loss of branching coral species such as P. acuta. A subsequent bleaching in 2016, however, revealed a mild bleaching response in pocilloporids compared to other coral taxa in the region, suggesting that rare, heat tolerant genotypes had been selected by the 2010 bleaching event. In order to test whether this potential 'evolutionary rescue' event has led to a low genetic diversity, we conducted a population genetic survey covering a total of nine different P. acuta populations (336 individuals) along a 50 km coastal stretch around Phuket Island, Thailand. We used six microsatellite markers to assess genotypic diversity and to determine the prevalent mode of reproduction (i.e. sexual or asexual recruitment). RESULTS: In contrast to other Indian Ocean P. acuta populations, the majority of corals in this study adopted a sexual reproduction mode (75% across all populations). At the same time, substantial regional gene flow was observed around Phuket Island with strong genetic differentiation as indicated by three genetic clusters that were separated by only a few kilometers. Patterns of isolation by distance over 0.7 - 40 km suggest small-scale genetic barriers, such as changing currents throughout each monsoonal season, potentially contributing to locally restricted dispersal of P. acuta larvae. CONCLUSIONS: The occurrence of distinct genetic clusters within short coastal stretches suggests that the 2010 bleaching event has not led to extreme genetic impoverishment. While more in-depth genomic analyses are necessary to investigate changes in genetic diversity following extreme bleaching events, our results will help guide conservation efforts to maintain genetic diversity of a coral species that likely will be dominant in future, warmer Andaman Sea reefs.

Ocean current patterns drive the worldwide colonization of eelgrass (Zostera marina).

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

Origin matters: Using a local reference genome improves measures in population genomics.

Genome sequencing enables answering fundamental questions about the genetic basis of adaptation, population structure and epigenetic mechanisms. Yet, we usually need a suitable reference genome for mapping population-level resequencing data. In some model systems, multiple reference genomes are available, giving the challenging task of determining which reference genome best suits the data. Here, we compared the use of two different reference genomes for the three-spined stickleback (Gasterosteus aculeatus), one novel genome derived from a European gynogenetic individual and the published reference genome of a North American individual. Specifically, we investigated the impact of using a local reference versus one generated from a distinct lineage on several common population genomics analyses. Through mapping genome resequencing data of 60 sticklebacks from across Europe and North America, we demonstrate that genetic distance among samples and the reference genomes impacts downstream analyses. Using a local reference genome increased mapping efficiency and genotyping accuracy, effectively retaining more and better data. Despite comparable distributions of the metrics generated across the genome using SNP data (i.e. π, Tajima's D and FST ), window-based statistics using different references resulted in different outlier genes and enriched gene functions. A marker-based analysis of DNA methylation distributions had a comparably high overlap in outlier genes and functions, yet with distinct differences depending on the reference genome. Overall, our results highlight how using a local reference genome decreases reference bias to increase confidence in downstream analyses of the data. Such results have significant implications in all reference-genome-based population genomic analyses.

A trait-based framework for seagrass ecology: Trends and prospects.

In the last three decades, quantitative approaches that rely on organism traits instead of taxonomy have advanced different fields of ecological research through establishing the mechanistic links between environmental drivers, functional traits, and ecosystem functions. A research subfield where trait-based approaches have been frequently used but poorly synthesized is the ecology of seagrasses; marine angiosperms that colonized the ocean 100M YA and today make up productive yet threatened coastal ecosystems globally. Here, we compiled a comprehensive trait-based response-effect framework (TBF) which builds on previous concepts and ideas, including the use of traits for the study of community assembly processes, from dispersal and response to abiotic and biotic factors, to ecosystem function and service provision. We then apply this framework to the global seagrass literature, using a systematic review to identify the strengths, gaps, and opportunities of the field. Seagrass trait research has mostly focused on the effect of environmental drivers on traits, i.e., "environmental filtering" (72%), whereas links between traits and functions are less common (26.9%). Despite the richness of trait-based data available, concepts related to TBFs are rare in the seagrass literature (15% of studies), including the relative importance of neutral and niche assembly processes, or the influence of trait dominance or complementarity in ecosystem function provision. These knowledge gaps indicate ample potential for further research, highlighting the need to understand the links between the unique traits of seagrasses and the ecosystem services they provide.

Detecting clonemate pairs in multicellular diploid clonal species based on a shared heterozygosity index.

Clonal reproduction, the formation of nearly identical individuals via mitosis in the absence of genetic recombination, is a very common reproductive mode across plants, fungi and animals. To detect clonal genetic structure, genetic similarity indices based on shared alleles are widely used, such as the Jaccard index, or identity by state. Here we propose a new pairwise genetic similarity index, the SH index, based on segregating genetic marker loci (typically single nucleotide polymorphisms) that are identically heterozygous for pairs of samples (NSH ). To test our method, we analyse two old seagrass clones (Posidonia australis, estimated to be around 8500 years old; Zostera marina, >750 years old) along with two young Z. marina clones of known age (17 years old). We show that focusing on shared heterozygosity amplifies the power to distinguish sample pairs belonging to different clones compared to methods focusing on all shared alleles. Our proposed workflow can successfully detect clonemates at a location dominated by a single clone. When the collected samples involve two or more clones, the SH index shows a clear gap between clonemate pairs and interclone sample pairs. Ideally NSH should be on the order of approximately ≥3000, a number easily achievable via restriction-site associated DNA (RAD) sequencing or whole-genome resequencing. Another potential application of the SH index is to detect possible parent-descendant pairs under selfing. Our proposed workflow takes advantage of the availability of the larger number of genetic markers in the genomic era, and improves the ability to distinguish clonemates from nonclonemates in multicellular diploid clonal species.

Effects of changing environmental conditions on plastic ingestion and feeding ecology of a benthopelagic fish (Gadus morhua) in the Southwest Baltic Sea.

This study documents how the abundance of microplastics (<5 mm) in the Atlantic cod, Gadus morhua, relates to the changes of the fish diet during years with contrasting levels of anoxia for example following years of low or high major Baltic inflows (MBI). A MultiNet Maxi trawl and CTD were deployed annually to collect microplastic samples alongside oxygen, temperature, and salinity conditions. Microplastics were homogenously distributed both within the water column and across years. Gadus morhua diet shifted from dominantly benthic invertebrates (61 %) under oxygenated conditions to dominantly Sprattus sprattus (81 %) under anoxic conditions. The proportion of G. morhua with microplastics in their digestive tract increased when they fed on pelagic fish (38 %) versus on benthic invertebrates (15 %). The proportion of S. sprattus which ingested microplastics (~18 %) did not vary. As anoxia at depth is expected to increase due to climate change, microplastic ingestion by G. morhua will potentially increase.

Invasion genomics uncover contrasting scenarios of genetic diversity in a widespread marine invader.

Invasion rates have increased in the past 100 y irrespective of international conventions. What characterizes a successful invasion event? And how does genetic diversity translate into invasion success? Employing a whole-genome perspective using one of the most successful marine invasive species world-wide as a model, we resolve temporal invasion dynamics during independent invasion events in Eurasia. We reveal complex regionally independent invasion histories including cases of recurrent translocations, time-limited translocations, and stepping-stone range expansions with severe bottlenecks within the same species. Irrespective of these different invasion dynamics, which lead to contrasting patterns of genetic diversity, all nonindigenous populations are similarly successful. This illustrates that genetic diversity, per se, is not necessarily the driving force behind invasion success. Other factors such as propagule pressure and repeated introductions are an important contribution to facilitate successful invasions. This calls into question the dominant paradigm of the genetic paradox of invasions, i.e., the successful establishment of nonindigenous populations with low levels of genetic diversity.

Improved chromosome-level genome assembly and annotation of the seagrass, Zostera marina (eelgrass).

Background: Seagrasses (Alismatales) are the only fully marine angiosperms.  Zostera marina (eelgrass) plays a crucial role in the functioning of coastal marine ecosystems and global carbon sequestration. It is the most widely studied seagrass and has become a marine model system for exploring adaptation under rapid climate change. The original draft genome (v.1.0) of the seagrass  Z. marina (L.) was based on a combination of Illumina mate-pair libraries and fosmid-ends. A total of 25.55 Gb of Illumina and 0.14 Gb of Sanger sequence was obtained representing 47.7× genomic coverage. The assembly resulted in ~2000 unordered scaffolds (L50 of 486 Kb), a final genome assembly size of 203MB, 20,450 protein coding genes and 63% TE content. Here, we present an upgraded chromosome-scale genome assembly and compare v.1.0 and the new v.3.1, reconfirming previous results from Olsen et al. (2016), as well as pointing out new findings.   Methods: The same high molecular weight DNA used in the original sequencing of the Finnish clone was used. A high-quality reference genome was assembled with the MECAT assembly pipeline combining PacBio long-read sequencing and Hi-C scaffolding.  Results: In total, 75.97 Gb PacBio data was produced. The final assembly comprises six pseudo-chromosomes and 304 unanchored scaffolds with a total length of 260.5Mb and an N50 of 34.6 MB, showing high contiguity and few gaps (~0.5%). 21,483 protein-encoding genes are annotated in this assembly, of which 20,665 (96.2%) obtained at least one functional assignment based on similarity to known proteins.  Conclusions: As an important marine angiosperm, the improved  Z. marina genome assembly will further assist evolutionary, ecological, and comparative genomics at the chromosome level. The new genome assembly will further our understanding into the structural and physiological adaptations from land to marine life.

Evolution via somatic genetic variation in modular species.

Somatic genetic variation (SoGV) may play a consequential yet underappreciated role in long-lived, modular species among plants, animals, and fungi. Recent genomic data identified two levels of genetic heterogeneity, between cell lines and between modules, that are subject to multilevel selection. Because SoGV can transfer into gametes when germlines are sequestered late in ontogeny (plants, algae, and fungi and some basal animals), sexual and asexual processes provide interdependent routes of mutational input and impact the accumulation of genetic load and molecular evolution rates of the integrated asexual/sexual life cycle. Avenues for future research include possible fitness effects of SoGV, the identification and implications of multilevel selection, and modeling of asexual selective sweeps using approaches from tumor evolution.

Phenotypic plasticity under rapid global changes: The intrinsic force for future seagrasses survival.

Coastal oceans are particularly affected by rapid and extreme environmental changes with dramatic consequences for the entire ecosystem. Seagrasses are key ecosystem engineering or foundation species supporting diverse and productive ecosystems along the coastline that are particularly susceptible to fast environmental changes. In this context, the analysis of phenotypic plasticity could reveal important insights into seagrasses persistence, as it represents an individual property that allows species' phenotypes to accommodate and react to fast environmental changes and stress. Many studies have provided different definitions of plasticity and related processes (acclimation and adaptation) resulting in a variety of associated terminology. Here, we review different ways to define phenotypic plasticity with particular reference to seagrass responses to single and multiple stressors. We relate plasticity to the shape of reaction norms, resulting from genotype by environment interactions, and examine its role in the presence of environmental shifts. The potential role of genetic and epigenetic changes in underlying seagrasses plasticity in face of environmental changes is also discussed. Different approaches aimed to assess local acclimation and adaptation in seagrasses are explored, explaining strengths and weaknesses based on the main results obtained from the most recent literature. We conclude that the implemented experimental approaches, whether performed with controlled or field experiments, provide new insights to explore the basis of plasticity in seagrasses. However, an improvement of molecular analysis and the application of multi-factorial experiments are required to better explore genetic and epigenetic adjustments to rapid environmental shifts. These considerations revealed the potential for selecting the best phenotypes to promote assisted evolution with fundamental implications on restoration and preservation efforts.

A novel metabarcoding primer pair for environmental DNA analysis of Cephalopoda (Mollusca) targeting the nuclear 18S rRNA region.

Cephalopods are pivotal components of marine food webs, but biodiversity studies are hampered by challenges to sample these agile marine molluscs. Metabarcoding of environmental DNA (eDNA) is a potentially powerful technique to study oceanic cephalopod biodiversity and distribution but has not been applied thus far. We present a novel universal primer pair for metabarcoding cephalopods from eDNA, Ceph18S (Forward: 5'-CGC GGC GCT ACA TAT TAG AC-3', Reverse: 5'-GCA CTT AAC CGA CCG TCG AC-3'). The primer pair targets the hypervariable region V2 of the nuclear 18S rRNA gene and amplifies a relatively short target sequence of approximately 200 bp in order to allow the amplification of degraded DNA. In silico tests on a reference database and empirical tests on DNA extracts from cephalopod tissue estimate that 44-66% of cephalopod species, corresponding to about 310-460 species, can be amplified and identified with this primer pair. A multi-marker approach with the novel Ceph18S and two previously published cephalopod mitochondrial 16S rRNA primer sets targeting the same region (Jarman et al. 2006 Mol. Ecol. Notes. 6, 268-271; Peters et al. 2015 Mar. Ecol. 36, 1428-1439) is estimated to amplify and identify 89% of all cephalopod species, of which an estimated 19% can only be identified by Ceph18S. All sequences obtained with Ceph18S were submitted to GenBank, resulting in new 18S rRNA sequences for 13 cephalopod taxa.

From ecosystems to socio-economic benefits: A systematic review of coastal ecosystem services in the Baltic Sea.

Seagrass meadows, algal forests and mussel beds are widely regarded as foundation species that support communities providing valuable ecosystem services in many coastal regions; however, quantitative evidence of the relationship is scarce. Using the Baltic Sea as a case study, a region of significant socio-economic importance in the northern hemisphere, we systematically synthesized the primary literature and summarized the current knowledge on ecosystem services derived from seagrass, macroalgae, and mussels (see animated video summary of the manuscript: Video abstract). We found 1740 individual ecosystem service records (ESR), 61% of which were related to macroalgae, 26% to mussel beds and 13% to seagrass meadows. The most frequently reported ecosystem services were raw material (533 ESR), habitat provision (262 ESR) and regulation of pollutants (215 ESR). Toxins (356 ESR) and nutrients (302 ESR) were the most well-documented pressures to services provided by coastal ecosystems. Next, we assessed the current state of knowledge as well as knowledge transfer of ecosystem services to policies through natural, social, human and economic dimensions, using a systematic scoring tool, the Eco-GAME matrix. We found good quantitative information about how ecosystems generated the service but almost no knowledge of how they translate into socio-economic benefits (8 out of 657 papers, 1.2%). While we are aware that research on Baltic Sea socio-economic benefits does exist, the link with ecosystems providing the service is mostly missing. To close this knowledge gap, we need a better analytical framework that is capable of directly linking existing quantitative information about ecosystem service generation with human benefit.

Modeling eelgrass spatial response to nutrient abatement measures in a changing climate.

For many coastal areas including the Baltic Sea, ambitious nutrient abatement goals have been set to curb eutrophication, but benefits of such measures were normally not studied in light of anticipated climate change. To project the likely responses of nutrient abatement on eelgrass (Zostera marina), we coupled a species distribution model with a biogeochemical model, obtaining future water turbidity, and a wave model for predicting the future hydrodynamics in the coastal area. Using this, eelgrass distribution was modeled for different combinations of nutrient scenarios and future wind fields. We are the first to demonstrate that while under a business as usual scenario overall eelgrass area will not recover, nutrient reductions that fulfill the Helsinki Commission's Baltic Sea Action Plan (BSAP) are likely to lead to a substantial areal expansion of eelgrass coverage, primarily at the current distribution's lower depth limits, thereby overcompensating losses in shallow areas caused by a stormier climate.

The Native Microbiome is Crucial for Offspring Generation and Fitness of Aurelia aurita.

All multicellular organisms are associated with microbial communities, ultimately forming a metaorganism. Several studies conducted on well-established model organisms point to immunological, metabolic, and behavioral benefits of the associated microbiota for the host. Consequently, a microbiome can influence the physiology of a host; moreover, microbial community shifts can affect host health and fitness. The present study aimed to evaluate the significance and functional role of the native microbiota for life cycle transitions and fitness of the cnidarian moon jellyfish Aurelia aurita A comprehensive host fitness experiment was conducted studying the polyp life stage and integrating 12 combinations of treatments with microbiota modification (sterile conditions, foreign food bacteria, and potential pathogens). Asexual reproduction, e.g., generation of daughter polyps, and the formation and release of ephyrae were highly affected in the absence of the native microbiota, ultimately resulting in a halt of strobilation and ephyra release. Assessment of further fitness traits showed that health, growth, and feeding rate were decreased in the absence and upon community changes of the native microbiota, e.g., when challenged with selected bacteria. Moreover, changes in microbial community patterns were detected by 16S rRNA amplicon sequencing during the course of the experiment. This demonstrated that six operational taxonomic units (OTUs) significantly correlated and explained up to 97% of fitness data variability, strongly supporting the association of impaired fitness with the absence/presence of specific bacteria. Conclusively, our study provides new insights into the importance and function of the microbiome for asexual reproduction, health, and fitness of the basal metazoan A. auritaIMPORTANCE All multicellular organisms are associated with a diverse and specific community of microorganisms; consequently, the microbiome is of fundamental importance for health and fitness of the multicellular host. However, studies on microbiome contribution to host fitness are in their infancy, in particular, for less well-established hosts such as the moon jellyfish Aurelia aurita Here, we studied the impact of the native microbiome on the asexual reproduction and on further fitness traits (health, growth, and feeding) of the basal metazoan due to induced changes in its microbiome. We observed significant impact on all fitness traits analyzed, in particular, in the absence of the protective microbial shield and when challenged with marine potentially pathogenic bacterial isolates. Notable is the identified crucial importance of the native microbiome for the generation of offspring, consequently affecting life cycle decisions. Thus, we conclude that the microbiome is essential for the maintenance of a healthy metaorganism.