Small RNAs couple embryonic developmental programs to gut microbes

Embryogenesis has long been known for its robustness to environmental factors. Although developmental tuning of embryogenesis to the environment experienced by the parent may be beneficial, little is understood on whether and how developmental patterns proactively change. Here, we show that Caenorhabditis elegans undergoes alternative embryogenesis in response to maternal gut microbes. Harmful microbes result in altered endodermal cell divisions; morphological changes, including left-right asymmetric development; double association between intestinal and primordial germ cells; and partial rescue of fecundity. The miR-35 microRNA family, which is controlled by systemic endogenous RNA interference and targets the β-transducin repeat-containing protein/cell division cycle 25 (CDC25) pathway, transmits intergenerational information to regulate cell divisions and reproduction. Our findings challenge the widespread assumption that C. elegans has an invariant cell lineage that consists of a fixed cell number and provide insights into how organisms optimize embryogenesis to adapt to environmental changes through epigenetic control.


INTRODUCTION
The process of embryogenesis has long been known to be extraordinarily invariable, or "canalized" (1), against extrinsic perturbations in many model systems (2)(3)(4). In viviparous species, such as mammals, embryogenesis proceeds while directly receiving substances from the mother, and it can often respond to the maternal environment, as seen in embryonic diapause upon starvation (5). However, because embryogenesis in most organisms occurs in isolation from the external environment, it is considered difficult for the embryo to directly perceive environmental factors in many cases, with few exceptions such as temperature (6). Therefore, little has been investigated about the extent to which the embryo alters its developmental patterns to adapt to various environmental changes.
The development of the nematode Caenorhabditis elegans is particularly renowned for its robustness; it is widely assumed to have exactly 959 somatic cells (1033 in males), of which 558 (560 in males) are produced during embryogenesis (4). Such highly invariable development allows for the assessment of developmental changes induced by environmental changes, genetic variations, and stochastic noises, making it useful for the study of developmental robustness and plasticity. C. elegans has been reported to alter its development in the presence of environmental stressors. When the environment is harsh for its growth during the first larval stage, it forms a dauer larva, an alternative third larval stage specialized for survival (7). Stresses such as high temperature and starvation also cause cell fate changes in vulval cells (8) and seam cells (9), although these changes are presently considered to be developmental "errors" (8,9) rather than adaptive developmental tuning. The ability of C. elegans to change its developmental patterns during embryogenesis has not been explored.
Many studies have reported that parental experiences of various environmental factors affect the progeny in terms of metabolism, behavior, and stress resistance, and environmental cues converge onto chromatin modifications and small RNAs, such as short interfering RNAs (siRNAs) and PIWI-interacting RNAs (piRNAs), in the germ line to generate and transmit epigenetic information (10,11). siRNAs are amplified via an RNA-directed RNA polymerase chain reaction, and in some animals, piRNAs are also amplified via a "ping-pong" cycle, allowing them to function as heritable epigenetic memory to transmit information across generations (10,11). However, without a known amplification mechanism, it is unclear whether microRNAs (miRNAs) can serve as heritable epigenetic memories without being diluted. Therefore, much attention has been directed toward examining how miRNAs may trigger long-lasting chromatin modifications [e.g., (12)].
Gut microbes have attracted a great deal of attention in the past decade owing to their widespread impact on health and disease risk (13). Recent studies in mice have shown that depletion of the maternal gut microbiota by antibiotic treatment or germ-free rearing leads to deficiencies in various maternal metabolites, which can cause abnormal neurodevelopment (14) and metabolic disorders (15) in progeny. Although it is unclear whether these microbe-free conditions can reflect the natural environment in which the animal lives, these studies raise the possibility that changes in symbiotic gut microbes can influence the embryonic development of progeny.

Developmental plasticity in C. elegans embryogenesis
To precisely characterize the effects of environmental factors on embryogenesis, we traced the developmental behaviors of almost all cells in C. elegans embryos using long-term, single-cell level, whole-embryo imaging under various environmental conditions (see Materials and Methods). We found that when C. elegans mothers had experienced mildly harmful microbes, such as Microbacterium nematophilum (16) CBX102 and Enterococcus faecalis (17) OG1RF, the number of intestinal nuclei in their embryos was increased (Fig. 1, A to C). The extra cell nuclei were generated by cell division but not by binucleation of cells because they were surrounded by the cell membrane (40 of 40 nuclei; fig. S1A).
The C. elegans intestine consists of nine intestinal rings (int1 to int9) arranged in an anterior-to-posterior sequence. The first int ring (int1) contains four cells, and the other int rings (int2 to int9) each contain two left and right cells (intXL and intXR) (18). All 20 intestinal cells are generated from the endoderm founder cell termed "E" during embryogenesis (4,18). Our imaging and subsequent cell tracking in embryos whose mothers had experienced CBX102 showed that in most cases, the extra divisions induced by CBX102 occurred at either int7L (Eplpa) or int7R (Eprpa), or both, increasing the number of intestinal cells to 21 (Fig. 1B). These extra cell divisions occurred almost simultaneously with, or slightly after, the increase in the intestinal cell number from 16 to 20 ( Fig. 1B and fig. S1C). In the int7 divisions, there is a left-right asymmetric tendency for int7R to divide along the anterior-posterior axis more often than int7L (fig. S1D). Binucleation, which occurs in all int3 to int7 cells and some int8 and int9 cells during the postembryonic development (18), appeared to occur also in extra cells (18 of 18 cells; fig. S1E).
Changes in cell numbers caused by CBX102 appeared to be specific to the endodermal lineage. We manually traced the entire cell lineage in three embryos with CBX102-induced extra division in the E lineage up to the 500-to 550-cell stage. Other than the extra divisions in the E lineage, all cell divisions and programmed cell deaths up to that stage were normal when compared to Sulston's reference lineage ( fig. S2).
Extra cells were also induced by Sphingobacterium multivorum BIGb0170 and Chryseobacterium scophthalmum JUb44 (Fig. 1C), both of which are part of the natural microbiota of C. elegans that colonize the gut and retard the growth of C. elegans (19), implying that the developmental changes also occur in the natural habitat of C. elegans. Escherichia coli OP50, the standard food for C. elegans in laboratory culture but reported to be weakly toxic (17,20), had a slight but detectable activity to induce extra cells (Fig. 1C). Sulston et al. (4) noticed that an extra cell was "occasionally" produced in the intestine when they described the entire embryonic cell lineage. Our results suggest that their observation is an environmentally induced phenomenon by OP50 rather than intrinsic variability due to less robust control of the lineage. In contrast, no extra cells were observed on the benign bacteria Comamonas aquatica (21) DA1877 and Bacillus subtilis (20) 1A1 (Fig. 1, C and D).
How do microbes induce the developmental changes? Given the diversity of the microbial species, we reasoned that it is more likely to be a physiological state of the worm than common chemicals from the microbes. Highly pathogenic Pseudomonas aeruginosa (17) PA14 and the hard-to-eat bacterium Bacillus megaterium (21) DA1880 did not induce extra cells (Fig. 1C), suggesting that general health, developmental speed, or nutrient deficiency in the mother is not a critical determinant. We noticed that all the microbes that induced extra cells are described as viscous (22,23) and/or cause "constipation," "distention," or "lumen bloating" in the C. elegans intestine ( fig. S3A) (16,(24)(25)(26). It has recently been suggested that bloating of the intestinal lumen is a "danger signal" that signifies microbial colonization (24). We examined defecation-defective mutants that are known to exhibit intestinal bloating, namely, aex-1, aex-2, and aex-5 (27,28), and found that their intestinal cell numbers increased on OP50 (Fig. 1E). Furthermore, restricting the entry of microbes into the gut tube (pharynx, intestine, and rectum) in mothers suppressed extra cells induced by CBX102, and this suppression was unaffected when worms were crossed with males with normal feeding ( fig. S3, B to D), suggesting that the microbe acts in the maternal gut. We next exposed parental worms to B. subtilis (natto). B. subtilis (natto) is classified as the same species as the model organism B. subtilis but is characterized by the production of a hard-to-digest biofilm composed of the biopolymer -poly-dlglutamic acid (PGA), a major constituent of the sticky strings of the Japanese food Natto and an important virulence factor for anthrax (29,30). B. subtilis (natto) strains strongly induced extra cells in C. elegans embryos (Fig. 1D). B. subtilis (natto) mutants defective in PGA synthesis (29) did not cause extra cells (Fig. 1D), indicating that PGA is required for the induction, although we cannot rule out the possibility that the mutations may affect additional aspects of the bacteria. B. subtilis (natto) also slowed the growth of worms in a PGA-dependent manner ( fig. S3E). Together, these results suggest that the maldigestion of maternal gut microbes that are highly proliferative and/or produce a viscous biofilm is an important factor for this developmental plasticity.

Consequences of alternative embryogenesis
To better understand the effect of maternal gut microbes, we further examined the embryonic development and fitness of the progeny that had extra cells. First, extra cells are accompanied by changes in the symmetry of intestinal development. In embryos that had extra cells, the right-side intestinal cells near the mid-intestine moved anteriorly, resulting in a left-right asymmetric cell arrangement (Fig. 2, A and B, fig. S4A, and movies S1 and S2). For example, int5R was paired with int4L ( Fig. 2, B and C) instead of the normal int5L, and int5R was even slightly anterior to int4L in ~60% of embryos ( fig. S4B). This left-right asymmetric cell arrangement seemed to be resolved by the inclusion of int3R in the int2 ring ( fig. S4C). The leftright asymmetry in the directions of the extra int7 division axis (fig. S1D) appears to be consistent with the left-right asymmetric cell arrangement.
Second, embryos with extra cells increase the physical association between the endoderm and primordial germ cells (PGCs), a process that is widely conserved in invertebrate and vertebrate animals (31)(32)(33)(34). From the 88-cell stage to the L1 larval stage, C. elegans has two PGCs (4), namely, Z2 and Z3. In embryos without extra cells whose mothers experienced CBX102, Z2 and Z3 took on an hourglass shape and each inserted a lobe into int5R and int5L, respectively (Fig. 2D), as reported in canonical embryogenesis (4,18,32). In embryos with extra cells, PGCs formed additional lobes, which associate with other mid-intestinal cells ( Third, worms with extra cells show partial rescue of fecundity when CBX102 is present; they produced >50% more progeny than worms without extra cells ( Fig. 2G and fig. S6A). In contrast, on OP50, they did not produce more progeny ( Fig. 2H and fig. S6B). Their developmental speed on CBX102 was comparable to that of worms without the extra cells ( fig. S6C), implying that the reproductive fitness was not due to better general health caused by the increased number of intestinal cells. Although the extent to which altered PGC morphology contributes to fecundity is unknown, the three observations suggest that maternal gut microbes lead to an alternative developmental program with extra cells and reproductive fitness.

Epigenetic control of alternative development
The C. elegans embryo is isolated from the external environment by the eggshell, and thus, information on microbes is highly likely to be transmitted through the maternal germ line. To understand how the information is processed in the mother and the embryo, we investigated the developmental phenotype of mutants of genes implicated in stress responses. We found that endogenous RNA interference (endo-RNAi), RNAi caused by endogenous siRNAs (35,36), is involved in this plasticity. Mutants of rde-4, eri-1, and rrf-3, all of which are components of the enhanced RNAi (ERI) complex required for endo-RNAi biogenesis (37,38), showed increases in intestinal cells on OP50 (Fig. 3A), suggesting that endo-RNAi acts to prevent the extra divisions. This result may be consistent with an increased intestinal cell nuclei number caused by RNAi of eri-5 (39), which is another ERI component (37,38). In contrast, mutants of prg-1, a PIWI Argonaute involved in the piRNA pathway (40,41), did not show an increase in intestinal cells (Fig. 3A). In C. elegans, intercellular transport of double-stranded RNAs (dsRNAs) leads to the propagation of RNAi from cell to cell (systemic RNAi), including transmission from soma to germ line (42). Systemic RNAi requires the systemic RNA interference deficiency-1 (SID-1) dsRNA importer (43)(44)(45)(46), the worm ortholog of mammalian SID-1 transmembrane family member 1 (SIDT1) and SIDT2, and SID-1-dependent intercellular propagation has been suggested to function in endo-RNAi (47,48). Mutants of sid-1 also showed increases in intestinal cells (Fig. 3A).
Conditional knockout of rrf-3 using somatic CRISPR-Cas9 (49) in the maternal intestine, but not in maternal neurons or body wall muscles, induced extra cells on OP50 (Fig. 3B), making the intestine a candidate as the source of endo-RNAi. Furthermore, the sid-1 mutant phenotype was rescued by germline-specific expression of sid-1 driven from a single-copy transgene (Si[pie-1 prom ::sid-1]) (Fig. 3C), implying a role for soma-to-germline transmission of dsRNAs. The rescue effect was observed in transgene-free embryos (−/−) born to heterozygous mothers (Si/−) (Fig. 3C), and sid-1 mutants crossed with wild-type males also showed an increase in extra cells ( fig.  S7A), indicating that the maternal action of sid-1 is important. Last, mutations in rde-4, rrf-3, or sid-1 did not enhance the extra cells induced by CBX102 ( fig. S7B), which is consistent with the notion that systemic endo-RNAi is the major response pathway to the microbe.
Next, we focused on miRNAs of the miR-35 family (miR-35-42, hereafter referred to as miR-35 fam ) as potential epigenetic regulators, because their expression is highly enriched in the germ line (50) and they positively regulate intestinal cell numbers (51). In mir- [35][36][37][38][39][40][41] and mir-35-41(gk262) deletion mutants (51,52), both of which lack a gene cluster consisting of seven of the eight family members, CBX102-induced extra cells were completely suppressed (Fig. 4A). Furthermore, introduction of nEx1187, a multicopy mir-35 expression array (52,53), induced extra cells on OP50 (Fig. 4B). This induction of extra cells was observed in array-negative embryos born to array-positive mothers (Fig. 4B), indicating that increased maternal mir-35 expression is sufficient for the regulation. Notably, nEx1187 also led to an improved fecundity on harmful microbes. We cultured mothers with the nEx1187 array on OP50 and isolated array-negative embryos with or without extra cells. When they were raised on CBX102, worms with extra cells produced more progeny ( Fig. 4C and fig. S8A). Reverse transcription quantitative polymerase chain reaction (RT-qPCR) showed that mature miR-35-3p and miR-40-3p, which are miRNAs expressed from the mir-35-41 cluster, were increased in the germ line when worms had experienced CBX102 or B. subtilis (natto) BVb750 (Fig. 4D). The miRNAs were not detected in the mir-35-41(nDf50) deletion mutants, indicating the specificity of our qPCR (Fig. 4D). The introduction of nEx1187 also increased the miR-35-3p expression in the maternal germ line ( fig. S8B), which may be due to leaky expression from the array that evades germline silencing. Mutants of the endo-RNAi pathway and the nEx1187 strain showed normal defecation ( fig. S8, C and D), showing that these genes act in parallel or downstream of intestinal stress. Together, these results show that the expression level of miR-35 fam in the maternal germ line increases in response to the microbes and that miR-35 fam up-regulation activates the alternative developmental program including extra cells and fitness.
Next, we examined zygotic expression of miR-35 fam using wwIs8[mir-35-41 prom ::GFP], a transcriptional reporter driven by the promoter of the mir-35-41 cluster (54), which is silenced in the maternal germ line but active in embryos. The expression of the reporter was elevated in embryos whose mothers had experienced CBX102 (Fig. 4E and fig. S9A). Furthermore, the reporter expression was reduced in the mir-35-41(nDf50) mutant background (Fig. 4F). Thus, similar to the positive feedback-based epigenetic inheritance of chromatin structure (55), the long-lasting response transmitted from mother to mid-stage embryo may be enabled by an autoregulatory feedback loop of miR-35 fam , where miR-35 fam positively controls its own expression. The reporter expression was increased in mutants of rrf-3 and sid-1 (Fig. 4G) S9D). While the lin-23(zb11) mutation did not lead to any gross phenotype, it, as well as a hypomorphic mutation in cdc-25.2, suppressed the nEx1187-induced extra cells (Fig. 4B), suggesting that the LIN-23/CDC-25 pathway acts downstream of miR-35 fam (Fig. 4H), although the incomplete suppression of extra

DISCUSSION
We propose here that maternal exposure to harmful gut microbes can lead to an alternative embryogenesis program with adaptive advantage and that the alternative program is epigenetically controlled by the endo-RNAi pathway and the miR-35 miRNA family (Fig. 4H). These findings challenge the widespread assumption that C. elegans has an invariant cell lineage. Our results suggest that the maldigestion of maternal gut microbes, which causes intestinal distention and poses a great survival risk to worms (16,(24)(25)(26), triggers the developmental plasticity. However, because of our current technological Our results suggest that intergenerationally inherited miR-35 fam responses induce specific endodermal cell divisions and adaptive advantages by controlling cell cycle regulators. Assuming that intergenerationally induced miRNA responses can be an entity of epigenetic memory by a simple transcriptional feedback loop, they may greatly expand the current repertoire of intergenerational epigenetic controls that allow sequence-specific gene regulation. While the cell cycle regulators LIN-23 and CDC-25 are ubiquitously expressed in the embryo to exert pleiotropic functions, small changes in their activities result in organ-specific abnormalities in endodermal cell divisions (57,58). Our results are consistent with these known phenotypes and the notion that miRNAs tune the expression of their targets within a range (59,60). The cells that divide excessively in response to maternal gut microbes (fig. S1B) are sisters of the cells that divide when the number of endodermal cells increases from 16 to 20 (Fig. 1B). Such a lineage relationship may be key to creating the specificity of the developmental changes.
The interactions between animals and microbes in their environments are complex. Some are beneficial, such as maternal gut microbes in mice providing essential metabolites for normal embryonic development (14,15). In contrast, many studies, including a report that germ-free mice can live 1.5 times longer than normal (61), have suggested that animals are adversely affected by chronic exposure to some indigenous microbes. Our results show that such stress can induce developmental tuning in animals with mosaic (determinate) development. Given the obvious plasticity of regulative development and the ability of the microbiome to affect systemic signaling (13), our study raises the question of whether such exposure could induce adaptive tuning of embryogenesis in humans and other complex animals.
In most animals, including humans, the development of the egg (and the embryo in viviparous animals) occurs near the intestine. This proximity makes the mother-to-embryo nutrition transport efficient, but it also has the disadvantage that the effects of harmful gut microbes can easily reach the next generation. For example, in humans, there is a great risk of Listeria monocytogenes vertically infecting the fetus in utero from the intestine, causing miscarriage and neonatal death (62). Future research may further reveal the mechanisms that protect the next generation of cells from harmful gut microbes through developmental changes.

Strains and culture
Standard methods were used to construct and culture C. elegans strains (63)(64)(65)(66). Worms were cultured on nematode growth medium (NGM) plates containing nystatin (50 U/ml; Sigma-Aldrich) at 20° to 21°C. Bristol N2 was used as the wild type. The strains used in this study are listed in table S1.

Embryo mounting and fluorescence microscopy
Preparation and mounting of embryos were performed as previously described (67)(68)(69) with modifications. Briefly, ~20 gravid adult worms were picked to ~20 l of Boyd buffer [5 mM Hepes (pH 7.2), NaCl (3.5 g/liter), KCl (2.4 g/liter), Na 2 HPO 4 (0.4 g/liter), CaCl 2 (0.2 g/liter), MgCl 2 (0.2 g/liter), and 0.2% glucose] and washed three times with Boyd buffer. Worms were then cut with a 25-gauge needle (BD PrecisionGlide), and early embryos at the 2-to 30-cell stage were collected with an eyelash. Embryos were transferred with a glass needle and mounted in ~1.5 l of Boyd buffer containing ~50 polystyrene-based microspheres (20 m in diameter, Polysciences) between a 24 mm by 50 mm glass slide and an 18 mm by 18 mm glass coverslip, which was sealed at the edges with melted Vaseline.
Images of live embryos were acquired with an Olympus UPlanSAPO 60× silicone oil objective [numerical aperture (NA), 1.30] or an Olympus UPlanSAPO 40× silicone oil objective (NA, 1.25). The microscopes used were (i) a spinning-disk confocal microscope (Quorum Technologies Inc.) consisting of a Zeiss Axio Observer Z1, a Yokogawa CSU-X1 spinning-disk unit, and two Hamamatsu C9100-13 EM-CCD cameras and (ii) an instant structured illumination microscope (Visitech iSIM) consisting of an Olympus IX73 and a Hamamatsu Flash 4.0v2 scientific complementary metaloxide-semiconductor (sCMOS) camera. Total fluorescence of the mir-35-41::GFP (green fluorescent protein) reporter (Fig. 4, E to G) was measured in ~30-cell stage embryos using Fiji software (National Institutes of Health).

Cell lineage analyses
For cell lineage analyses, embryos were imaged with 1-m z-steps across 29 m and with 75-s intervals. Positions of cell nuclei and timing of cell divisions were determined with StarryNite and AceTree software (68)(69)(70). Directions of cell divisions and distance between cells were calculated from the central coordinates of the cell nuclei. Three-dimensional reconstruction of cell nuclei positions in embryos was performed by WormGUIDES software (69,71).

Exposure of worms to a microbial lawn
For preculture of microbial strains, a freshly streaked colony was inoculated and grown at 30°C in LB broth with the exceptions of OG1RF, which was cultured in brain heart infusion broth (Difco), and JUb44 and BIGb0170, which were cultured at 25°C. After 24 hours, 100 l of preculture was spread onto a 60-mm NGM plate and then incubated at 20° to 21°C for 24 hours.
NGM plates for culturing parental worms were seeded with OP50 and incubated at 20° to 21°C for 2 days. The plates were stored at 4°C for 2 to 20 days until the day of use. Parental worms were cultivated on the OP50-seeded plates until L4 larval to young adult stages. The parental worms were collected and washed five times with M9 buffer [KH 2 PO 4 (3 g/liter), Na 2 HPO 4 (6 g/liter), NaCl (5 g/liter), 1 mM MgSO 4 , and 0.03% gelatin], and then ~50 worms were transferred to a plate seeded with each indicated microbial strain. After 24 hours of incubation at 20° to 21°C, their embryos were collected and mounted as described above. Ivermectin was first dissolved in dimethyl sulfoxide at a concentration of 10 g/ml and added to melted NGM agar to give a final concentration of 1 ng/ml before pouring.

Quantification of intestinal cell numbers
Numbers of intestinal cells were counted in 1.5-to 3-fold stage embryos using the fluorescence of stIs10453[elt-2 prom ::elt-2::gfp] or the combination of rrIs1[elt-2 prom ::nls::gfp] and ujIs113[pie-1 prom ::mCherry::H2B, nhr-2 prom ::mCherry::his -24]. Because the expression of rrIs1 was variable and occasionally undetectable in some cells, the additional fluorescence of ujIs113 was required to accurately count intestinal cells. The distribution of the number of intestinal cells in each trial is shown in data S1.
Intestinal cell numbers in cross progeny (figs. S3C and S7A) were counted as follows. him-8(−) males containing the stIs10453[elt-2 prom :: elt-2::gfp] marker and him-8(+) hermaphrodites lacking any fluorescent marker were raised until mid-L4 stage on OP50-seeded plates. The males and hermaphrodites were separately collected, washed five times with M9 buffer, and transferred to NGM plates on which CBX102 or OP50 was grown for 18 hours. After 6 hours of separated preexposure, 50 males and 10 hermaphrodites were picked to a new NGM plate on which CBX102 or OP50 was grown for 24 hours. After 24 hours of incubation at 20° to 21°C, embryos were collected from the hermaphrodites and mounted as described above. Intestinal cells were observed only in the cross progeny, which carried the stIs10453 marker.
Brood size, growth, and defecation assays L4 and young adult worms (P0) grown on OP50 were collected and washed five times with M9 buffer, and then ~50 worms were transferred to a plate seeded with CBX102 (Fig. 2, G and H, and fig. S6) or with OP50 ( Fig. 4C and fig. S8A). After 30 hours, 1.5-fold embryos (F 1 ) were directly picked into ~50 l of M9 buffer dropped on a 24 mm by 50 mm glass slide, and their ELT-2::GFP fluorescence was observed with the spinning-disk confocal microscope described above. Embryos without extra intestinal cell (E20), embryos with one extra intestinal cell produced by either int7L or int7R divisions (E21), and embryos with two extra intestinal cells produced by both int7L and int7R divisions (E22) were isolated. The F 1 embryos were individually washed with 1% sodium hypochlorite to avoid microbial contamination and then plated onto NGM plates seeded with OP50 or onto NGM plates seeded with a mixture of stationary-phase cultures of CBX102 and OP50 mixed in a volumetric ratio of 1:9 (10% CBX102 plates).
For brood size assays, the isolated F 1 embryos were grown until young adults, and all the F 2 hatchlings were killed and counted every 24 hours until the F 1 worms no longer produced progeny. For growth assays on CBX102 ( fig. S6C), the isolated F 1 embryos were allowed to develop at 20°C on 10% CBX102 plates, and their developmental stages were determined after 72 and 80 hours.
For growth assays on the E. coli and B. subtilis strains ( fig. S3E), early embryos were prepared by hypochlorite treatment of gravid hermaphrodites raised on OP50 and were transferred onto a new plate seeded with an indicated E. coli or B. subtilis strain. The embryos were allowed to develop at 20°C for 65 hours, and their developmental stages were determined. Defecation cycle length, the duration between the first muscular contraction of one defecation and the first muscular contraction of the next defecation, was measured in 1-day-old young adult worms at 20°C under a dissecting microscope.

Isolation of gonadal miRNAs and RT-qPCR
For isolation of gonadal miRNAs, young adults were transferred into egg buffer [25 mM Hepes (pH 7.2), 118 mM NaCl, 48 mM KCl, 2 mM CaCl 2 , and 2 mM MgCl 2 ] and washed five times with egg buffer immediately before dissection. The gonads were isolated from 100 worms by cutting the worms behind the pharyngeal bulb and in front of the spermatheca with a 25-gauge needle. The isolated gonads were placed into 300 l of TRIzol on ice. Sixty microliters of chloroform was added to it, and after centrifugation, the aqueous phase was transferred into a new tube. After adding 1.5 volumes of ethanol, the samples were loaded into a miRNeasy column (Qiagen) and purified according to the supplier's instructions. The purified samples were reverse-transcribed using the miScript reverse transcriptase and HiSpec buffer (Qiagen).
For isolation of mRNAs, approximately 20 1-day-old gravid adults were placed into 300 l of TRIzol and frozen. After thawing, total RNA was purified as described above with an RNeasy column (Qiagen). The purified samples were reverse-transcribed using Quantitect reverse transcription (Qiagen).

Statistical analysis
Statistical analyses were performed with a statistic package (Prism v.9.0.2, GraphPad software). Error bars indicate SEM. Statistical comparisons were performed with two-tailed one-way analysis of variance (ANOVA) with Dunnett's posttest or Tukey's posttest, two-tailed Kruskal-Wallis with Dunn's posttest, two-tailed t test, or two-tailed Mann-Whitney test. All experiments were repeated three to seven times with similar results. Raw data that support the findings of this study are available on Figshare with the DOI 10.6084/m9. figshare.18133622.v2.