- 1Department of Biology, Department of Biology, Utah State University, Logan, UT, United States
- 2United States Department of Agriculture, Agricultural Research Service, Pollinating Insects Research Unit, Logan, UT, United States
Maternal effects are an important source of phenotypic variance, whereby females influence offspring developmental trajectory beyond direct genetic contributions, often in response to changing environmental conditions. However, relatively little is known about the mechanisms by which maternal experience is translated into molecular signals that shape offspring development. One such signal may be maternal RNA transcripts (mRNAs and miRNAs) deposited into maturing oocytes. These regulate the earliest stages of development of all animals, but are understudied in most insects. Here we investigated the effects of female internal (body condition) and external (time of season) environmental conditions on maternal RNA in the maturing oocytes and 24-h-old eggs (24-h eggs) of alfalfa leafcutting bees. Using gene expression and WGCNA analysis, we found that females adjust the quantity of mRNAs related to protein phosphorylation, transcriptional regulation, and nuclease activity deposited into maturing oocytes in response to both poor body condition and shorter day lengths that accompany the late season. However, the magnitude of these changes was higher for time of season. Females also adjusted miRNA deposition in response to seasonal changes, but not body condition. We did not observe significant changes in maternal RNAs in response to either body condition or time of season in 24-h eggs, which were past the maternal-to-zygotic transition. Our results suggest that females adjust the RNA transcripts they provide for offspring to regulate development in response to both internal and external environmental cues. Variation in maternal RNAs may, therefore, be important for regulating offspring phenotype in response to environmental change.
1 Introduction
Females can influence offspring development in ways that are independent of direct genetic inheritance, typically in response to changing environmental cues (Bernardo, 1996a; Wolf and Wade, 2009). The mechanisms underpinning such maternal effects include both pre- and post-zygotic functions. For example, females can influence the developmental rate, size, or sex of their offspring through post-zygotic mechanisms such as choice of nest site, incubation frequency, or provisioning rate (Klostermeyer et al., 1973; Torchio and Tepedino, 1980; Bernardo, 1996a; Meaney, 2001). Females can also influence the earliest stages of offspring development by adjusting the transcriptional and endocrine profiles of maturing oocytes (Wolf and Wade, 2009; Groothuis et al., 2019; Vastenhouw et al., 2019). Although these pre-zygotic mechanisms are likely to have early and ongoing effects on offspring phenotype, they have been relatively understudied compared to post-zygotic mechanisms, especially in insects. This is particularly true with regard to how females interpret environmental cues and translate them into molecular signals that influence offspring development (Huestis and Marshall, 2006; Lee and Duvall, 2022). Furthermore, although it is likely that maternal effects are influenced by multiple cues (Marshall and Uller, 2007), these cues are typically studied in isolation (though see Potticary and Duckworth (2020)). Studies assessing the effects of multiple cues—a more realistic view of environmental influences—will yield a better understanding of non-genetic drivers of phenotypic variance, which provide the raw material for adaptive evolution and are an important source of ecological diversity (Mousseau and Dingle, 1991; Bernardo, 1996a; Mousseau and Fox, 1998; Räsänen and Kruuk, 2007).
Cues influencing pre-zygotic maternal effects on offspring phenotype can come from the external or internal environment. External cues experienced by most females are both biotic and abiotic, including interactions with natural enemies (Agrawal et al., 1999; Rolff, 1999; Mitchell and Read, 2005; Sharda et al., 2021) and exposure to changing weather conditions (Bernardo, 1996a; Burgess and Marshall, 2011). Photoperiod is a common environmental cue, because it can serve as a reliable indicator of seasonal change (Bradshaw and Holzapfel, 2007), and can thus serve as a coordinating mechanism for critical life history events such as reproduction, mating, migration, and diapause (Mousseau and Dingle, 1991; Bradshaw and Holzapfel, 2007). For example, maternal photoperiod influences egg size, diapause, and survival outcomes of offspring in mosquitoes (Lacour et al., 2014; Lee and Duvall, 2022). Similarly, maternal photoperiod is the primary determinant of whether larvae enter diapause in multiple fly species (Saunders et al., 1986; McWatters and Saunders, 1997). There is also evidence for transgenerational effects of photoperiod on offspring diapause, development time, and several morphological traits in a parasitoid wasp (Tougeron et al., 2020). Together, the results of these studies suggest female response to external cues such as photoperiod is an important driver of maternal effects in insects.
Environmental cues that shape maternal effects may also stem from changes in one or more internal conditions. These may include hormones (for reviews see Edwards et al., 2021; Groothuis and Schwabl 2008; Groothuis et al., 2019; Meylan et al., 2012), body size (Steiger, 2013), and body condition (de Zwaan et al., 2019). Interestingly, maternal nutritional status can also impact offspring development. In non-biting midges, females reared in food stressed conditions yielded offspring that developed faster and had decreased fecundity relative to those given resources in excess (Colombo et al., 2014). Similar effects of parental nutritional status have been observed in mosquitoes, where offspring of nutritionally-stressed parents were more likely to transmit the dengue virus (Zirbel et al., 2018). In neriid flies, offspring of females in poor body condition developed faster than those from females in good body condition (Bonduriansky and Head, 2007). Finally, parasitoid wasps experiencing high levels of competition, which may decrease body condition, were more likely to produce diapausing offspring (Tougeron et al., 2018). Based on these studies, it is clear that females can influence the same set of offspring traits (e.g., development rate, fecundity, and even diapause) in response to different sets of internal and external environmental cues. How females integrate these cues to influence offspring development via maternal effects is unknown. Resolving this relationship will require perspective on how cues from the environment are translated into signals that regulate development in offspring.
Understanding how maternal experience (i.e., internal and external cues) is translated into variation in offspring phenotype requires insight into the mechanisms by which females can influence the developmental trajectory of their offspring. One of the earliest stages at which females can influence offspring phenotype is during oogenesis. During this stage, small shifts in maternal input of hormones (Navara et al., 2006), yolk (or nutrition) (Bernardo, 1996b), and RNA (Vastenhouw et al., 2019) can have large, organizational effects on the developing zygote. In insects, maternal RNAs are transcribed in the supporting nurse cell and transferred into the growing oocyte via cytoskeletal machinery (Spradling, 1993). Maternal RNAs can include both protein-coding (mRNA) and small, non-coding regulatory RNAs (e.g., microRNA [miRNA]), and are necessary for the final stages of oocyte maturation and activation and the earliest stages of embryonic development following fertilization (Winata and Korzh, 2018; Vastenhouw et al., 2019). These maternally-derived molecules are responsible for regulating critical processes in early embryogenesis such as cell structure and division, biosynthesis, blastula formation, and gastrulation (Weeks and Melton, 1987; Tang et al., 2007; Baroux et al., 2008; Tadros and Lipshitz, 2009; Pauli et al., 2011; Harvey et al., 2013; Paranjpe et al., 2013; Winata and Korzh, 2018; Torres-Paz et al., 2019). Maternal RNAs are programmatically degraded and cleared in a dynamic process leading up to the maternal-to-zygotic transition (MZT), which occurs during the blastoderm stage in insects (Sung et al., 2013; Pires et al., 2016). After the MZT, developmental processes, including RNA transcription, are controlled by the embryo.
Given their influential role in the earliest stages of development, variation in maternal RNAs are a potential way that females translate environmental cues into maternal effects on offspring phenotype. Research with fish has revealed extensive variation in how females deposit RNA transcripts in their eggs. For example, there are significant differences in the relative abundance of maternally-deposited transcripts among zebrafish females, but almost no variation within a female’s clutch (Rauwerda et al., 2016). This could suggest that female zebrafish make consistent adjustments in how they deposit transcripts during oogenesis based on their internal condition. In support of this, bacterial supplementation treatments in zebrafish modified maternal condition and influenced the composition of maternal mRNAs (Miccoli et al., 2017, 2015). Additionally, in the round goby, females adjust their RNA contributions to embryos in response to water temperature (Adrian-Kalchhauser et al., 2018). In the annual killifish, females respond to environmental cues via maternal programming of mRNAs and miRNAs deposited in oocytes that determine whether young embryos develop directly or enter diapause (Romney and Podrabsky, 2017). Lastly, in both cichlids and cavefish, differences in maternal RNA provisioning play a causal role in generating the phenotypic novelty that promotes local adaptation and species diversification (Ahi et al., 2018; Torres-Paz et al., 2019). Overall, the potential for fish to alter offspring development and impact phenotype through environmentally-sensitive maternal RNA provisioning is evident.
There is accumulating evidence that insects also transmit signals about the environment to developing offspring via maternally-derived RNAs. Egg diapause in the Asian tiger mosquito (Aedes albopictus) is maternally regulated based on exposure of adults to short day lengths (Wang, 1966; Mori et al., 1981). Experimental studies have demonstrated that females adjust the composition of mRNAs (Poelchau et al., 2011), but not miRNAs (Batz et al., 2017), deposited into mature oocytes in response to photoperiod regime, and these mRNA adjustments coincide with an increased likelihood of diapause among offspring (Poelchau et al., 2011). Further research suggests these differences in maternally-deposited mRNAs may lead to increasingly divergent expression profiles between diapause-destined and non-diapaused destined embryos, even past the MZT (Poelchau et al., 2013). Similarly, mRNA differences have been found in the ovaries of female locusts (Locusta migratoria) exposed to short- and long-days (Hao et al., 2019), which is known to trigger maternally-mediated diapause in locust eggs (Tanaka, 1994). This suggests that maternally provisioned RNAs may be a common factor in the initiation of diapause in the egg stage. However, insects can diapause during any stage of the life cycle and there is a diverse array of molecular signals that trigger diapause across species (Denlinger, 2002, 2022). A complete understanding of maternal effects on insect diapause requires investigating species that diapause at different life stages.
We investigated the mechanistic underpinnings of maternal effects in alfalfa leafcutting bees (Megachile rotundata). Various aspects of alfalfa leafcutting bee biology make them an ideal species for investigating how maternal perception of the environment stemming from internal and external cues can shape maternal effects on offspring development. First, M. rotundata naturally exhibit facultative diapause as late-stage larvae (“prepupae”) (Krunic, 1972; Hobbs and Richards, 1976; Tepedino and Parker, 1988, 1986; Pitts-Singer and Cane, 2011; Pitts-Singer, 2020); hence, a readily-observable dichotomy in developmental trajectory already exists among offspring. Second, maternal influence has long been recognized as a factor driving this facultative diapause (Parker and Tepedino, 1982; Tepedino and Parker, 1986; Johnson, 2022). Therefore, their life cycle offers a unique opportunity for evaluating how maternal environmental cues are transmitted to the offspring. Third, two maternal cues, photoperiod and lipid stores (i.e., external and internal, respectively) are known to influence the probability of offspring diapause (Pitts-Singer, 2020; Wilson et al., 2021; Johnson, 2022). Adult M. rotundata females experience seasonal changes (early to late summer) in photoperiod, and manipulations of maternal day length experience has been shown to alter patterns of diapause induction in the offspring (Pitts-Singer, 2020; Wilson et al., 2021). Likewise, maternal condition is also a significant determinant of diapause fate in progeny (Johnson, 2022). Females regulate offspring size via the amount of food provided in the brood cell (Klostermeyer et al., 1973), which impacts the probability of survival during diapause (Fischman et al., 2017). Because brood cell construction and provisioning is energetically demanding (Klostermeyer and Gerber, 1969; Klostermeyer et al., 1973), maternal body condition is likely to influence offspring diapause destiny. Jointly, these features make M. rotundata useful for disentangling how females integrate external and internal cues of the environment to direct offspring development via maternal RNA.
To test how external signals and internal condition influence the mechanisms of maternal effects, we explored the impacts of time of season (i.e., photoperiod) and maternal body condition (i.e., depletion of lipid stores) on maternal RNA in M. rotundata. Here, time of season includes the combined effects of all changes that occur during a season, including photoperiod, but also pests, temperature, floral resources, etc. To determine whether and how these cues influenced maternal RNAs, we quantified mRNA and miRNA from maturing oocytes and eggs that were 24 h post-oviposition, and, therefore, past the maternal-to-zygotic transition (hereafter referred to as “24-h eggs”). In our previous research, we found that females with experimentally reduced lipid stores had fewer diapausing offspring in both the early season and the late season. However, the probability of offspring diapause was significantly higher in the late season than in the early season. We thus hypothesized that M. rotundata females manipulate their RNA provisions as a response to both external and internal cues of the environment in a way that is consistent with known maternal effects on offspring diapause outcomes. Our results reveal a potential mechanism for observable maternal effects in bees.
2 Methods
2.1 Study organism and field collections
Megachile rotundata are solitary bees that are intensively managed for alfalfa seed production (Pitts-Singer and Cane, 2011; Calderone, 2012; Pitts-Singer, 2020). They readily nest in above-ground, artificial tunnels made of wood or polystyrene when they are used as commercial pollinators (Fairey and Lieverse, 1986; Pitts-Singer and Cane, 2011). In the western U.S., most offspring overwinter as cocooned, non-feeding fifth instar larvae (prepupae), but almost half the progeny avert this process, emerging as adults in the same summer as the parent generation (Krunic, 1972; Hobbs and Richards, 1976; Tepedino and Parker, 1988, 1986; Pitts-Singer and Cane, 2011; Pitts-Singer, 2020).
For this study, we received bee cells containing diapausing prepupae from farms in the alfalfa seed-growing areas of Box Elder County, Utah, USA during October 2019 and kept them in cold storage at 4°C until diapause was broken by incubating the cocoons at 29°C. To manipulate female body condition, we interrupted development during incubation by temporarily moving bees into cool storage (18°C) for either 1 day (control) or 14 days (poor condition) before returning them to 29°C until adult emergence. This type of interruption is a commonly used management technique for aligning the timing of bee emergence and field release with alfalfa bloom (Richards, 1984). Previous research has shown that bees can survive a prolonged setback in temperature for 2 weeks at 18°C (Yocum et al., 2010). Our recent research demonstrated that this developmental pause significantly reduces lipid quantity at emergence, and females who have experienced this interruption produce significantly fewer diapausing offspring and provide smaller provision masses for offspring than females that experience only a 1-day interruption (Johnson, 2022). Additionally, late season females produce significantly more diapausing offspring than early season females, most of whose offspring are non-diapausers (Johnson, 2022). To assess seasonal effects at the molecular level, we released bees from each body condition treatment in both the early and late summer. Average day length was 15.22 min longer during our early season experiment (
Our experimental design allowed the release of females from both the poor body condition and control treatments in both early and late season. On the day of emergence, adult females were given a unique thoracic paint mark (Testors, Rockford, IL, USA) and then released into two 6.1 × 6.1 × 1.8 m screened cages erected over flowering alfalfa in Logan, Utah, United States. For both early and late season releases, all bees were placed into their respective cages within a 1–2 days period of each other. Male M. rotundata were released simultaneously in a ratio of 2:1 male to female (Rossi et al., 2010).
Each cage was provided with a small section of prefabricated, polystyrene nesting block mounted 1.1 m above the ground in the center of the cage. This allowed us to easily observe nesting activity. Cage assignments were random, but with equal numbers of control and poor body condition treatments in each cage. To compare how maternal RNA abundances respond to variation in body condition and time of season, we collected maturing oocytes and 24-h eggs that were past the MZT. Nests were monitored three times daily to verify the identity of each nest owner. Additionally, we conducted hourly nest checks between 10:00 and 19:00 to see if any new eggs had been laid or were about to be laid. These hourly checks included watching for signs of nest initiation and pollen provisioning that are done just before egg laying. If an egg was laid, but the cell holding the egg was not fully closed, we waited to collect the female until after she capped the cell. Once a female was known to have laid an egg, she was captured as soon as possible and flash frozen in liquid nitrogen within minutes. We later dissected the most mature oocyte from her ovaries (see below). Eggs still in the nest were left undisturbed for 24 h, then frozen in liquid nitrogen, and later removed from the nest on dry ice prior to RNA isolation. We chose 24 h for egg collection as this is post-MZT transition for insect eggs, including honey bees (Tadros and Lipshitz, 2009; Ninova et al., 2016; Pires et al., 2016; Vastenhouw et al., 2019). These experimental methods provided us with maturing stage 4 oocytes and 24-h eggs from females who experienced experimental lipid reduction and control females, in both the early and late season.
2.2 Sample processing
All samples were maintained at -80°C until dissection. Prior to dissections, abdomens were detached from the thorax while on dry ice and incubated in pre-chilled RNAlater-ICE Frozen Tissue Transition Solution (Thermo Fisher Scientific, Waltham, MA, USA) at -20°C. After a 16–20 h incubation period, abdomens were dissected in RNAlater-ICE using a Leica M80 stereomicroscope with an IC80HD camera attached (Leica Microsystems, Buffalo Grove, IL, USA). While dissecting, we identified the longest terminal oocyte, removed it from the ovary, and transferred the oocyte to dry ice. Dissected oocytes were stored at -80°C until RNA isolation. Oocyte, trophocyte, and total egg chamber length was measured with the Leica Application Suite (v.4.5) software. We calculated the proportion of the egg chamber occupied by the oocyte and used these percentages, along with various egg chamber characteristics, to assign oocyte maturation stage as described by Kapheim and Johnson (2017). We selected oocytes from the same stage of maturation for each experimental group. This yielded 22 stage 4 oocytes occupying 79–95% of the egg chamber: ‘early control’, 85.8–92.8%, n = 6; ‘early poor condition’, 78.9–89.7%, n = 4; ‘late control’, 84.0–94.6%, n = 6; ‘late poor condition’, 84.1–94.9%, n = 6.
We isolated total RNA from a single oocyte (n = 22) and egg (n = 22) from each female in our study to yield 44 total samples. We used the mirVana miRNA Isolation Kit with phenol (Ambion, Austin, TX, United States) according to the manufacturer’s protocol. The total RNA yield of isolates ranged from 0.52 to 1.54 μg (
All dissections, ovary measurements, and RNA isolations were conducted without knowledge of sample experimental treatment group.
2.3 Sequencing
2.3.1 mRNA
RNAseq libraries were prepared with the Illumina TruSeq Stranded mRNAseq Sample Prep kit (Illumina). The libraries were pooled, quantitated by qPCR, and sequenced on two SP lanes for 151 cycles from both ends of the fragments on a NovaSeq 6000. This generated 1,832,084,298 total 100 nt reads with a mean of 41,638,279.50 (±579,462.15 standard error) reads per sample. FASTQ files were generated and demultiplexed with the bcl2fastq v2.20 Conversion software (Illumina).
2.3.2 miRNA
Small RNA libraries were prepared using the Qiagen Small RNA Sample Prep kit. Libraries were pooled, quantitated by qPCR and sequenced on one SP lane for 51 cycles from one end of the fragments on a NovaSeq 6000. This generated 519,046,055 total 50 nt reads with a mean of 11,796,501.25 (±248,977.88 standard error) reads per sample. FASTQ files were generated and demultiplexed with the bcl2fastq v2.20 Conversion Software (Illumina). The 3’ adapters were trimmed.
2.4 Alignment and quantification
2.4.1 mRNA
After quality checks with MultiQC v1.5 (Ewels et al., 2016), we trimmed reads with Trimmomatic v.0.39 (Bolger et al., 2014) in Paired-End mode. We used the Illuminaclip function to remove adapters and the sliding window function to trim reads with an average quality score below 20, based on 4 nt windows. We used the ‘keep both reads’ option and filtered reads shorter than 51 nts. Trimmed reads were then aligned to the M. rotundata genome (Kapheim et al., 2015) using STAR v2.7.8a (Dobin et al., 2013) in paired end multi mode. We did not clip any bases from the read ends during alignment. Finally, we quantified the number of reads aligned to each CDS feature in the M. rotundata annotation v.1.1 (Kapheim et al., 2015) using the Python v3.6.3 script HTSeq-count v0.9.1 (Anders et al., 2015) in union mode with the stranded option set to reverse. The selection of these pipeline tools was made based on a comparative analysis of performance (Corchete et al., 2020).
2.4.2 miRNA
After quality checks with MultiQC v1.5 (Ewels et al., 2016), we implemented miRDeep2 for the quantification, alignment and identification of miRNA’s in M. rotundata. First, the M. rotundata genome (Kapheim et al., 2015) was indexed using bowtie-build. Then we ran each sample’s FASTQ miRNA through the mapper.pl script to align them to the M. rotundata genome. Reads shorter than 18 bases were discarded, non-nucleotide letters were removed (anything other than the following characters: ATCGUNatcgun), and reads were then collapsed. Next, we identified known miRNAs and discovered novel miRNAs using known precursor sequences from M. rotundata, and the known mature sequences from M. rotundata, Bombus impatiens, B. terrestris, Apis mellifera, Megalopta genalis and Nomia melanderi (Kapheim et al., 2020) as inputs alongside our sample sequences to the miRDeep2.pl script of miRDeep2 (Friedländer et al., 2011). Next, we ran the quantifier.pl script of miRDeep2 on found known and novel miRNAs for each M. rotundata individual sample to obtain read quantities for each pair of precursor and mature miRNA. We then filtered precursor-mature pairs to keep only those that were not rRNA/tRNAs, had a minimum of five reads each on mature and star strands of the hairpin sequence, and a randfold p < 0.05, following Kapheim et al. (2020). We then re-ran miRDeep.pl and quantifier.pl on the original sample set of mapped miRNAs. For the second run, we first compiled a list of the unique novel precursor miRNA sequences from across all samples which passed the aforementioned filtering criteria for the output from the first quantifier.pl run. This list of novel precursor miRNAs was added to the list of known precursor miRNAs of M. rotundata from Kapheim et al. (2020), and the composite was used as the species-specific precursor miRNA input for the second run of miRDeep2.pl, without changing any other inputs. We then applied the same filters to compile our final list of miRNAs.
For downstream analysis of differential miRNA expression, we treated miRNAs with identical mature sequences, but non-overlapping precursor locations, as separate miRNAs. Overlapping precursors with identical mature sequences were merged and considered the same miRNA. (Mapping to multiple precursor sequences did not occur within individuals.) miRNA homologs were identified with miRBase v22.1 (Griffiths-Jones et al., 2006). We used the search tool, powered by RNA Central, to identify the top hit for each mature miRNA sequence. We kept only those hits with query and target sequence match
2.5 Target prediction of miRNAs
We ran the target prediction software, miRanda v3.3 (Enright et al., 2003), on the final set of unique mature miRNAs from all samples (minimum energy threshold -20, minimum score 140, strict alignment to seed region -en -20 -sc 140 -strict). In addition, we analyzed this same set of mature miRNAs for potential target sequences using RNAhybrid v2.12 (Krüger and Rehmsmeier, 2006) (minimum free energy threshold -20). We kept only miRNA-target gene pairs that were predicted by both programs with p < 0.01.
2.6 Statistical analysis
2.6.1 mRNA
We combined all feature counts from the HTSeq-count output in R v4.0.2 (R Core Team, 2020) and used these to make a DGEList with the EdgeR package v3.30.3 (Robinson et al., 2010). Visual inspection of samples on a PCA (Supplementary Figure S1) revealed clear separation between oocyte and egg samples and the presence of one outlier (egg sample 202198). We therefore proceeded to analyze these separately after removing the outlier. For each dataset (oocyte and egg), we used the filterByExpr function to filter genes based on counts per million (CPM) and number of samples in each group. This kept 8,276 (64.8%) and 9,200 (72.0%) out of 12,770 genes for oocytes and eggs, respectively. This concurs with the proportion of protein-coding genes known to be maternally transcribed and deposited into maturing oocytes in Drosophila (Laver et al., 2015). Libraries were then normalized with the trimmed mean of M values (TMM) method. We evaluated gene expression differences as a function of seasonal period and maternal body condition in a non-intercept model (∼0 + group). The model design was implemented with the lmFit and eBayes functions applied to contrasts of interest after controlling the variance with voom (R package limma v3.44.3) (Ritchie et al., 2015). These analysis steps were chosen following current best practices (Law et al., 2016, 2020).
2.6.2 Shared significant DEGs with previous diapause studies
We compared genes responding to time of season or body condition in our study and those associated with diapause in other species. We compared our genes of interest with those differentially expressed in diapausing queens in B. terrestris (Supplementary Table S3D in Amsalem et al. (2015)), those related to diapause in larvae of the tropical oil-collecting bee T. diversipes (Supplementary Tables S1, S2 in Santos et al. (2018)), as well as those related to early and late diapause induction in pupae of M. rotundata in November (supplement S2 in Yocum et al. (2018)). We identified orthologs as reciprocal best BLAST hits (evalue cutoff of 10 × 10−5), using amino acid sequences from Mrot v1.1 and amino acid sequences converted from nucleotide sequences with TransDecoder v5.5.0 (https://github.com/TransDecoder/TransDecoder) or nucleotide sequences B. terrestris. We calculated representation factors and assessed significance of gene overlap using hypergeometric tests with the phyper function in R v. 4.2.1 (Hankin, 2016).
2.6.3 Weighted gene coexpression network analysis
We conducted a weighted gene correlation network analysis (WGCNA) on our expression data to explore relationships between clusters of highly coexpressed genes (i.e., gene modules) and our experimental variables (Langfelder and Horvath, 2008). To do so, we constructed scale-free gene coexpression networks using the R package WGCNA v1.70–3 (Langfelder and Horvath, 2008) for mRNA expression in egg (n = 22) and oocyte (n = 22) samples independently. Using the default filtering parameters of the goodSamplesGenes function, we removed 852 genes from the expression data set due to excessive missing data or zero variance across samples (n = 11,918 genes were retained for subsequent analyses). After filtering, we subset the data by sample type for tissue-specific analyses. We used hierarchical clustering (stats v3.6.1) to assess outliers. No outliers were identified among the oocyte samples, but two egg samples (202196 and 202198) were identified as outliers and removed (remaining: oocytes, n = 22; eggs, n = 20). Gene networks for oocytes and eggs were constructed at a soft power of 4, which appeared best suited for balancing scale independence (R2 = 0.624 and R2 = 0.529, respectively) and mean connectivity (303 and 458, respectively). In both cases, we computed a topological overlap matrix from the adjacency matrix using WGCNA default parameters. Then we began module assignment. WGCNA gene modules are clusters of densely interconnected genes (Langfelder and Horvath, 2008) defined here using an unsupervised hierarchical approach. Module membership for each gene is assigned using an eigengene-based connectivity approach; each gene is assigned to a single module by correlating the gene’s expression profile with the eigengene value of each module. Values close to one or -1 suggest high connectivity to a particular module (Langfelder and Horvath, 2008). We initially assigned modules with a minimum module threshold of 30 genes and deepSplit = 2. Modules with similar coexpression, i.e. those with a correlation
We assessed module overlap between oocytes and eggs for those modules that were significantly associated with either season or body condition (henceforth termed ‘modules of interest’; n = 4). One large egg module, eggmod4, had significant (p < 0.05) overlap with three of those oocyte modules of interest. We used bootstrapping methods to estimate the likelihood that this pattern of overlap could have occurred by chance. To do so, we randomly sampled (without replacement) gene IDs from the list of all genes placed into module assignments equal to the number of genes in the respective modules of interest over 10,000 iterations. For each sample run, the proportion of genes sampled from eggmod4 out of the total number of genes sampled was calculated. We then calculated the probability that the randomly sampled gene sets equal to the size of oocyte modules had a greater proportion of eggmod4 genes than the observed proportion for each of modules of interest.
2.6.4 miRNA
miRNA expression was analyzed similar to the mRNA expression (R v4.0.2 R Core Team 2020). We used the matrix of counts to create a DGEList object with the EdgeR package v3.30.3 (Robinson et al., 2010). Visual inspection of unsupervised clustering plots (Supplementary Figure S1) revealed clear separation between oocyte and egg samples. We therefore proceeded to analyze these separately. For each dataset (oocyte and egg), we used the filterByExpr function to filter genes based on counts per million (CPM) and number of samples in each group. This kept 46 (61.3%) and 50 (66.7%) out of 75 miRNAs for oocytes and eggs, respectively. Libraries were then normalized with the TMM method. We evaluated gene expression differences as a function of time of season and maternal body condition in a non-intercept model (∼0 + group). The model design was implemented with the lmFit and eBayes functions applied to contrasts of interest after controlling the variance with voom (R package limma v3.44.3) (Ritchie et al., 2015). These analysis steps were chosen following current best practices (Law et al., 2016, 2020).
2.6.5 Gene Ontology enrichment analyses: mRNA and miRNA
Prior to enrichment analyses, we re-annotated the M. rotundata gene set v.1.1 (Kapheim et al., 2015) using InterProScan v5.56-89.0 (Jones et al., 2014) with Pfam (Mistry et al., 2020), SUPERFAMILY (Gough et al., 2001), and PANTHER (Thomas et al., 2022) analyses selected (-appl Pfam, SUPERFAMILY, PANTHER). Of the 12,770 genes, GO terms were assigned to 6,354 genes. Using those annotations, we performed Gene Ontology enrichment analyses for differentially expressed (Benjamini–Hochberg (BH)-adjusted p < 0.1) mRNA, predicted targets of differentially expressed miRNA (BH-adjusted p < 0.1), and the genes belonging to WGCNA modules of interest using TopGO v2.46.0 (Alexa and Rahnenfuhrer, 2020) with default method (‘weight01’) and node size (nodesize = 1) parameters.
3 Results
3.1 Differentially expressed mRNAs
We detected 8,276 genes in oocytes and 9,200 genes in eggs. Only 28.5% of genes detected in oocytes and 26.5% of genes detected in eggs have orthologs that were identified in a previous study of honey bee maternal RNA (Supplementary File S1) (Pires et al., 2016). Most of these had been previously classified as zygotic in origin (class II) in honey bees. Differences in maternal experience were evident in genes and miRNAs expressed in oocytes, but not in eggs (Table 1; Figure 1). These effects were more apparent for differences in time of season than they were for maternal body condition (Figure 1, Supplementary Figure S2). In oocytes, we detected 19 differentially expressed genes in response to poor body condition (8 up-regulated, 11 down-regulated; BH-adjusted p < 0.05). These included genes involved in embryogenesis, including orthologs of branchless, knot, and juvenile hormone-inducible protein 26 (Supplementary File S1). To investigate additional genes that were less consistently differentially expressed, i.e., those missing the arbitrary p < 0.05 threshold, we relaxed α to BH-adjusted p < 0.1. Subsequently, the list grew to 35 differentially expressed genes (15 up-regulated, 20 down-regulated) (Table 1). These genes were enriched (p < 0.01) for enzyme regulator and proteasome-activating activity, as well as a transcription factor core complex (Supplementary Table S3; Supplementary File S2). Seasonal effects in the oocyte were evident with four genes significantly up-regulated in the late season (BH-adjusted p < 0.05). These included genes involved in circadian rhythms and metabolic processes, including orthologs of arylalkylamine N-acetyltransferase 1, knockdown, and chloride channel-b. When we expanded α to BH-adjusted p < 0.1, we identified 413 genes up-regulated and 543 down-regulated in the late season (956 total). These genes were enriched (p < 0.01) for protein phosphorylation, dephosphorylation, and processing, regulation of transcription, and protein tyrosine phosphatase activity, as well as signal transduction and proton transmembrane transport (Supplementary Table S4; Supplementary File S2). We did not detect any genes responding differently to body condition as a function of time of season, as would be indicated by the interaction term in our model. Only a single gene (Mrot07161) was identified as differentially expressed in eggs (Table 1). This was up-regulated in response to poor body condition. This result did not change when we adjusted the α to BH-adjusted p < 0.1. The D. melanogaster ortholog of Mrot07161 (omega) is involved in proteolysis, among other functions.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1064332/fgene-13-1064332-HTML-r1/image_m/fgene-13-1064332-t001.jpg)
TABLE 1. Number of differentially expressed mRNAs and miRNAs at a threshold of a BH-adjusted p < 0.05 and (in parentheses) p < 0.1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1064332/fgene-13-1064332-HTML-r1/image_m/fgene-13-1064332-g001.jpg)
FIGURE 1. Principal component analysis showing mRNA and miRNA expression profiles for (A,B) oocyte (blue color scheme) and (C,D) egg (green color scheme) samples. Each data point represents an individual sample. Manipulations of maternal body condition included both control (circle) and poor condition (triangle) treatment groups. Time of season includes early (light shade) and late (dark shade). Percent of variance captured by each principal component (PC) is included in the axes labels for PC1 and PC2.
Some of the maternal mRNAs that were differentially expressed in oocytes in response to maternal experience have been previously associated with the initiation or maintenance of diapause in insects. We compared genes differentially expressed between treatments in our oocyte tissue to those found differentially expressed in diapausing B. terrestris queens (Amsalem et al., 2015), diapausing T. diversipes larvae (Santos et al., 2018), and between fall-sampled M. rotundata prepupae that had entered diapause either in early and late summer (Yocum et al., 2018). The only significant overlap that we found was between our study and those related to diapausing B. terrestris queens (representation factor: 1.39; hypergeometric test: p = 0.07; Supplementary File S3).
3.2 Differentially expressed miRNAs
More than half (39) of the 72 unique mature miRNA sequences identified in our oocyte and egg samples of M. rotundata were previously found in miRNA expression data from brains of six bee species (Kapheim et al., 2020). Many (34) of the 72 unique mature miRNA sequences identified in our oocyte and egg samples were previously found in miRNA expression data of M. rotundata brains (Supplementary Figure S3) (Kapheim et al., 2020). While only 72 unique mature miRNA sequences were identified, we characterized 75 distinct miRNAs for differential expression analysis. miRNAs with precursor sequences which did not overlap, but were paired with a mature miRNA that met all aforementioned precursor-mature pair filtering criteria were considered separate miRNAs for differential expression analysis, even if the mature sequence of one pair was identical to another. Of the 75 miRNAs expressed in our dataset, all except two had homologs in miRBase, and 80% (60) were detected in previous studies of maternal RNA or the MZT in insects (Supplementary File S4) (Marco, 2015; Pires et al., 2016).
Differential expression analysis of miRNAs revealed a pattern similar to mRNAs (Figure 1). Poor body condition did not yield any significant changes in miRNA expression in oocytes or eggs. We identified 13 miRNAs that responded significantly (BH-adjusted p < 0.05) to time of season in oocytes, with two up-regulated and 11 down-regulated in the late season, as compared to the early season (Table 1). While most of these 13 were involved in gene silencing, other functions of these miRNAs included regulation of the circadian rhythm, embryogenesis, and the cellular response to hypoxia (Supplementary File S4). Of these 13 miRNAs, 85% (11) were homologous to maternal miRNAs expressed in mature honey bee oocytes and degraded in early embryogenesis (Class I or III) (Pires et al., 2016), and 62% (8) were homologous to maternal miRNAs expressed in D. melanogaster oocytes (Marco, 2015) (Supplementary File S4). We did not detect any miRNAs responding differently to body condition or time of season in the eggs (Table 1). Nor did we detect any miRNAs responding differently to body condition as a function of time of season in the oocytes or eggs, as would be indicated by the interaction term in our model.
We performed Gene Ontology enrichment of the set of predicted targets for the miRNAs differentially expressed across the season in the oocytes. This included a set of 164 unique genes, but there was no significant enrichment for any GO term (Supplementary File S5). There was no overlap between this set of predicted targets and the mRNA genes differentially expressed in oocytes in early and late season (BH-adjusted p < 0.1). However, more than half (56.3%) of the 375 unique predicted targets of miRNAs expressed in the eggs were genes expressed in oocytes, as would be expected if these miRNAs expressed in the early stages of embryogenesis function in clearing maternal transcripts (Marco, 2015).
3.3 Gene coexpression networks
The gene coexpression network for oocytes included 30 modules (Supplementary Table S1, Supplementary File S6), none of which were significantly correlated with cage or maternal body condition (Supplementary Figure S4). However, we identified four modules (modules of interest) that were significantly associated with time of season (oocyteMod4, r = -0.57, p = 0.005; oocyteMod29, r = -0.51, p = 0.02; oocyteMod16, r = -0.46, p = 0.03; oocyteMod20, r = -0.45, p = 0.04; Supplementary Figure S4). Genes belonging to these modules were enriched for various aspects of proton transmembrane transport, cell differentiation and cell fate determination, semaphorin receptor binding, vesicle-mediated transport, and tRNA activities, various types of substrate binding, and regulation of transcription (Supplementary Tables S5-S8, Supplementary File S7). The gene coexpression network for eggs included 27 modules, none of which were significantly correlated with cage, maternal body condition, or time of season (Supplementary Table S2; Supplementary Figure S5; Supplementary File S8).
To investigate how the function of maternal RNAs shifts across the MZT, we evaluated the degree of overlap between gene coexpression modules in the oocytes and the eggs. We found that the structure of oocyte modules were not maintained in 24-h eggs (Figure 2). However, one large egg module was more likely to include genes from the four oocyte modules of interest than expected by chance. Three oocyte modules (oocyteMod16, oocyteMod4, and oocyteMod29) had significant overlap with eggMod4, with 51.9% (p = 2.08x10−20), 45.1% (p = 0.00114), and 65.2% (p = 6.84x10−10) of their genes being assigned to this module in egg tissue, respectively (Supplementary Table S9, Supplementary Figure S6, Supplementary File S9). Permutation tests suggest this degree of overlap was unlikely due to the large size of eggMod4 alone (p < 0.0001 for oocyteMod16, p < 0.0016 for oocyteMod4, and p < 0.0001 for oocyteMod29).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1064332/fgene-13-1064332-HTML-r1/image_m/fgene-13-1064332-g002.jpg)
FIGURE 2. Dendrograms depicting gene module relationships for oocyte and egg tissues. Lines (moving left to right) indicate where genes from oocyte modules significantly correlated with time of season (oocyteMod4, oocyteMod16, oocyteMod20, oocyteMod29) are distributed throughout egg modules. The line color represents oocyte module membership, including gold, dark blue, turquoise, and magenta as oocyteMod4, oocyteMod16, oocyteMod20, oocyteMod29, respectively. Egg modules that shared significant overlap across these four oocyte modules are indicated with a gray box. The shaded box highlights (eggMod4), the egg module that shared a significant number of genes with three of these oocyte modules (oocyteMod16, oocyteMod4, and oocyteMod29). Line thickness scales to the percentage of genes from the oocyte modules that belong to each respective egg module. Line scale (from smallest to largest): <1%, 1–10%, 10–20%, 20–30%, 30–40%, 40–50%, 50–60%, and 60–70%.
4 Discussion
Maternal effects are an important source of phenotypic diversity, but the mechanisms by which females interpret cues from the environment and translate them into molecular signals that influence offspring development are understudied. We found that alfalfa leafcutting bee females respond to cues from both the internal and external environment by adjusting the composition of RNA they deposit into maturing oocytes. These adjustments included both mRNAs and miRNAs previously known to be maternally-derived and playing a role in the earliest stages of development in insects. These effects were no longer detectable within 24 h after egg-laying, which is consistent with current understanding of the maternal-to-zygotic transition. Our results shed light on the dynamic nature of how maternal experience can influence the developmental trajectory of her offspring.
Maternal mRNA provisioning is one of the earliest ways that a female can influence the phenotype of her offspring. These molecular signals are required for zygote genome activation in all animals (Berry, 1982; Lee et al., 2014) and have been shown to affect survival after egg fertilization in mice (Christians et al., 2000), cell differentiation in nematodes (Mello et al., 1992), polarity in both fruit flies and nematodes (Evans et al., 1994), and cell fate determination in frogs (Xanthos et al., 2001). Given these large, organizational effects on offspring development, we hypothesized that adjustments to maternal RNA composition or abundance is an important mechanism underlying maternal effects. Indeed, research with mosquitoes and locusts suggests that female insects can influence offspring diapause outcome by adjusting mRNA provisioning of maturing oocytes in response to environmental cues (Poelchau et al., 2011; Hao et al., 2019). We found evidence that female alfalfa leafcutting bees also make dynamic adjustments to maternal RNAs, suggesting this could be a widespread mechanism by which females alter offspring development using cues from the environment. This is compelling, because M. rotundata diapause at a later stage of development (prepupae) than mosquitoes and locusts, which both diapause as pharate larvae inside the egg chorion. Our study reveals that development of later-diapausing species may also be influenced by maternal effects on the oocytes.
Given the role of maternal RNAs in processes critical to development (i.e., polarity and cell differentiation), we anticipated relatively few of these genes to be responsive to environmental cues. Instead, we anticipated most would have a narrow range of expression constrained by their role in canonical functions. This is consistent with our results, whereby fewer than 20 mRNAs and miRNAs showed highly significant differences in expression in response to female body condition or time of season. However, far more genes had less consistent or lower magnitude changes in expression, as indicated by the larger number of genes with significant differences in expression at a higher α (BH-adjusted p < 0.1; Table 1). These genes were enriched for functions that indicate high-level organizational effects on cellular differentiation, including processes related to transcription and translation. Thus, changes in the expression of even a small number of these genes could have potentially large effects on developmental trajectory. Notably, our assay only detected mRNAs active at the time of sample collection. Maternal RNA expression is dynamically regulated leading up to the MZT by the shortening and lengthening of polyA tails, a process that silences and activates translation, respectively (Winata and Korzh, 2018). In this study, we used bead selection of polyA tails to generate RNA libraries for transcriptome sequencing. Hence, our results only capture the active portion of the maternal RNAs being passed to offspring. Additionally, we were unable to detect any adjustments that females made to the relative abundance of transcripts that were silenced until translation at later stages of the MZT. Thus, our results present a conservative picture of the dynamic changes possible for maternal RNAs.
The highly canonical role of maternal RNAs in development also suggests that the genes involved would be highly conserved across species. We found this to be the case for the miRNAs we detected in our samples, but less so for the mRNAs. Less than a third of the genes we detected in alfalfa leafcutting bee oocytes and eggs have orthologs that were identified in a previous study of honey bee maternal RNA (Pires et al., 2016). In contrast, the vast majority of the miRNAs we detected have homologs that were previously characterized as maternal in origin in honey bees or fruit flies (Marco, 2015; Pires et al., 2016). One potential explanation for this is that there are fewer conserved protein-coding mRNAs than miRNAs in the M. rotundata genome, making it more likely to detect conserved miRNAs due to chance. However, the M. rotundata genome has roughly equal proportions of lineage-specific miRNAs and protein-coding genes (Kapheim et al., 2015, 2020). An alternative explanation is that the early developmental processes regulated by maternally-deposited miRNAs are more conserved than those regulated by maternally-deposited mRNAs. Across organisms, maternal miRNAs act, in part, to clear maternal mRNAs during the MZT by targeting maternal mRNA transcripts that were active during oocyte maturation, but are no longer needed during embryogenesis (Bashirullah et al., 1999; Winata and Korzh, 2018). For example, in D. melanogaster the zygotically-expressed miR-309 cluster targets approximately 400 maternal mRNAs for degradation (Bushati et al., 2008). Consistent with this function, we found that more than half of the predicted targets of the miRNAs expressed in our 24-h eggs were genes expressed in the oocytes. Interestingly, this suggests that the function of maternal miRNAs are conserved between alfalfa leafcutting bees and other insects (Marco, 2015).
Conversely, we did not find evidence that maternal miRNAs and mRNAs share a coordinated response to environmental cues. There was no overlap between oocyte genes that were differentially expressed according to time of season and the predicted targets of oocyte miRNAs that were differentially expressed according to time of season. One potential explanation for this is a temporal lag in when these differences in expression would be detected. For instance, because it takes time for new target mRNAs to accumulate in the cytoplasm, a decrease in a mature miRNA might not correlate with an increase in target mRNAs. Thus, sampling mRNA profiles at additional time points would likely be necessary to observe the effects of environmentally-induced changes in miRNA expression on its target genes. However, our results are also consistent with previous work in Drosophila showing that only a small proportion of maternal transcripts are targeted by miRNAs during oocyte maturation (Nakahara et al., 2005). Further, there is some evidence that maternal transcripts have been selected to avoid targeting by maternal miRNAs, likely to preserve their function in oocyte maturation (Marco, 2015).
Although our results offer compelling evidence that females adjust the mRNA and miRNA transcripts that they deposit into maturing oocytes in response to environmental cues, how this translates into offspring phenotype remains an open question. Answering this question is necessary to determine whether variation in maternal RNA composition is a mechanism underlying maternal effects. There is some evidence in various species of fish that variation in maternal RNAs alters the developmental trajectory of the embryos, resulting in different offspring phenotypes (Romney and Podrabsky, 2017; Ahi et al., 2018; Torres-Paz et al., 2019). In mosquitoes, maternal adjustments to oocyte mRNA are concomitant with increased diapause incidence in offspring (Poelchau et al., 2011), and in locusts, knockdowns of maternal genes in the FOXO pathway have been causally linked to changes in offspring diapause (Hao et al., 2019). While direct evaluation of offspring phenotypes was outside the scope of our study, our results provide indirect evidence that variation in maternal RNAs in response to environmental cues could shape the developmental trajectory of alfalfa leafcutting bees. This evidence comes from what is already known about maternal effects in alfalfa leafcutting bees.
Whether leafcutting bee larvae enter diapause or not has long been recognized as being regulated by maternal effects (Johansen and Eves, 1973; Bitner, 1976; Parker and Tepedino, 1982; Rank and Goerzen, 1982; Tepedino and Parker, 1986; Kemp and Bosch, 2001), potentially triggered by changes in photoperiod (Pitts-Singer, 2020; Wilson et al., 2021). In a recent study using an experimental design similar to that used in this study, we found that the proportion of offspring entering diapause was significantly higher among females nesting in the late season and significantly lower among females in poor body condition following experimental reduction of lipid stores (Johnson, 2022). However, the mechanism regulating these effects appeared to vary. Females in poor body condition made smaller food provisions for their larvae, and size of food provision was a strong predictor of the diapause trajectory for individual larvae (Fischman et al., 2017; Johnson, 2022). This could indicate that post-zygotic mechanisms of maternal effects (i.e., nutrient provisioning) act in concert with pre-zygotic effects (i.e., maternal RNA provisioning) to influence offspring phenotype.
These differences in the size of larval food provisions were not evident between early and late season females, suggesting that the increased propensity for late season females to produce more diapausing offspring was more heavily influenced by pre-zygotic effects (i.e., maternal RNA provisioning) (Johnson, 2022). Leafcutting bee larvae develop inside leaf-wrapped cell cups inside dark cavities, and are therefore unable to sense photoperiod directly. Moreover, changes in temperature have been ruled out as potential cues that regulate diapause (Pitts-Singer, 2020; Wilson et al., 2021). While we cannot rule out other factors that correlate with time of season, such as changes in nutrient availability, these patterns are highly consistent with our findings—time of season leads to substantial significant differences in how females deposit mRNA and miRNA transcripts into maturing oocytes, but fewer significant differences associated with female body condition. Moreover, the genes differentially expressed in response to time of season had significant overlap with genes involved in diapause in another bee species and included genes and miRNAs with a known role in circadian rhythms and metabolism. Together, these results suggest that photoperiod-induced differences in maternally-provisioned RNA influence offspring propensity for diapause. Similar to this, environmentally-responsive maternal programming via mRNAs and miRNAs deposited in annual killifish oocytes determined whether young embryos develop directly or enter diapause (Romney and Podrabsky, 2017). Overall, this suggests that diapause, as a developmental process, may be particularly sensitive to variation in maternal RNAs across species.
In contrast to oocytes, there was virtually no signal of variation in maternal environment observed in 24-h eggs. This seems to suggest that whatever transcriptional adjustments females make to maturing oocytes do not persist beyond the MZT. This calls into question whether maternal RNAs are likely to have lasting effects on offspring phenotype. One possible explanation for this pattern is that the biochemical signatures of early embryogenesis dominated the transcriptomic profiles of 24-h eggs, making any signal of maternal effects either difficult to detect or non-existent. For example, previous research in butterflies suggests that hormone pulses demarcating developmental transitions elicit a greater influence on gene expression profiles than environmental influences, such as temperature (Tian and Monteiro, 2022). Our WGCNA analyses also showed signs of embryogenesis dominating transcriptional signatures in 24-h eggs. While four modules of coexpressed genes correlated significantly with time of season in the oocytes, there were none in the eggs that correlated with either treatment. Further, when comparing patterns of co-regulated genes, major reorganization occurs between oocytes and eggs, such that no module of co-regulated genes remained entirely intact when progressing from oocyte to egg (Figure 2). However, the fact that a significant proportion of genes belonging to the oocyte modules that responded to seasonal effects were coexpressed in the same egg module (eggMod4) could suggest that the maternal signal in oocytes is at least partially conserved in eggs. eggMod4 was exceptionally large, with the genes from those four oocyte modules only comprising about one-tenth of those belonging to it. Yet, this is still more than was expected by chance (see permutation tests). Together, these results suggest that, while molecular signatures of maternal experience did not appear to be maintained past the MZT (i.e., in the 24-h eggs), certain patterns of co-regulated genes may have persisted undetected against the overwhelming signal of transcriptional processes associated with embryogenesis.
However, it is also possible that the organizational effects of differential expression in maturing oocytes are no longer evident in transcriptomic data, but would be most readily observed in downstream processes. As such, effects of differentially expressed maternal RNAs in maturing oocytes might be more evident by examining metabolomic, proteomic, or endocrine profiles of eggs if expression changes have altered protein or metabolic production post-MZT. Future studies that incorporate additional types of biochemical data will be informative for understanding the effects of differences in maternal RNA expression on embryogenesis.
Understanding how maternal experiences are translated into developmental changes in offspring requires investigating the molecular mechanisms underlying maternal effects. We found that both internal and external environmental cues elicit changes in mRNA and miRNA deposition by female bees into their maturing oocytes. This is consistent with our prior work showing that maternal photoperiod and body condition impacts offspring diapause induction in this species. Our current findings reveal how those environmental signals may facilitate facultative diapause in an economically important pollinator. Dynamic adjustment of maternal RNA in response to environmental cues has previously been documented in egg-diapausing insects, but it has heretofore been unclear whether similar mechanisms operate in species that diapause at a later life stage. Thus, our results also provide new insights about insect developmental biology, as well as mechanisms contributing to phenotypic variation in a changing world.
Data availability statement
The datasets presented in this study can be found in the article/Supplementary Material. Sequences have been deposited at NCBI (PRJNA887543). Snakefiles (Sequence Alignment and Quantification and Reciprocal Best Hit Analysis) and R markdown (Statistical Analysis) files are available at Github: www.github.com/kapheimlab/mrotundata_maternaleffects_RNA.
Author contributions
Conceptualization: KK and TP-S; Data curation: MH and FH; Formal analysis: MH, FH, TD, and KK; Funding acquisition: KK and TP-S; Investigation: MH, FH, and MJ; Methodology: KK and TP-S; Project administration: KK; Resources: KK and TP-S; Supervision: KK; Visualization: MH and KK; Writing—original draft: MH, FH, TD, and KK; Writing—Review & Editing: MJ and TP-S.
Funding
This work was supported by the USDA NIFA Foundational Programs award #2018-67014-27542 (KK and TS), the NSF Graduate Research Fellowship Program (MJ), and Utah State University. This research was supported by the Utah Agricultural Experiment Station, Utah State University (Project #1297 (KK)) and approved as journal paper #9614. This research was supported by the U.S. Department of Agriculture, Agricultural Research Service.
Acknowledgments
We are grateful to the farmers for cultivating and allowing us access to their alfalfa fields. We thank E. Klomps and N. K. Boyle for valuable discussion and assistance in the field and A. Hernandez for RNA sequencing. We thank R. A. Weinstock for their feedback on early manuscript drafts. We thank K. Hartfelder, S. M. Strobel, M. C. Womack, L. L. Surber, J. R. Phillips, R. A. Weinstock for helpful discussion of our results. We thank P. K. Santos for translated protein FASTA files of T. diversipes. This research was supported by the U.S. Department of Agriculture, Agricultural Research Service. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.1064332/full#supplementary-material
References
Adrian-Kalchhauser, I., Walser, J.-C., Schwaiger, M., and Burkhardt-Holm, P. (2018). Rna sequencing of early round goby embryos reveals that maternal experiences can shape the maternal rna contribution in a wild vertebrate. BMC Evol. Biol. 18, 34–14. doi:10.1186/s12862-018-1132-2
Agrawal, A. A., Laforsch, C., and Tollrian, R. (1999). Transgenerational induction of defences in animals and plants. Nature 401, 60–63. doi:10.1038/43425
Ahi, E. P., Singh, P., Lecaudey, L. A., Gessl, W., and Sturmbauer, C. (2018). Maternal mrna input of growth and stress-response-related genes in cichlids in relation to egg size and trophic specialization. EvoDevo 9, 23–17. doi:10.1186/s13227-018-0112-3
Alexa, A., and Rahnenfuhrer, J. (2020). topGO: Enrichment Analysis for Gene Ontology. R package version 2.40.0.
Ambros, V., Bartel, B., Bartel, D. P., Burge, C. B., Carrington, J. C., Chen, X., et al. (2003). A uniform system for microrna annotation. Rna 9, 277–279. doi:10.1261/rna.2183803
Amsalem, E., Galbraith, D. A., Cnaani, J., Teal, P. E., and Grozinger, C. M. (2015). Conservation and modification of genetic and physiological toolkits underpinning diapause in bumble bee queens. Mol. Ecol. 24, 5596–5615. doi:10.1111/mec.13410
Anders, S., Pyl, P. T., and Huber, W. (2015). Htseq—A python framework to work with high-throughput sequencing data. bioinformatics 31, 166–169. doi:10.1093/bioinformatics/btu638
Baroux, C., Autran, D., Gillmor, C., Grimanelli, D., and Grossniklaus, U. (2008). “The maternal to zygotic transition in animals and plants,” in Cold Spring Harbor symposia on quantitative biology (Cold Spring Harbor Laboratory Press), 73, 89.
Bashirullah, A., Halsell, S. R., Cooperstock, R. L., Kloc, M., Karaiskakis, A., Fisher, W. W., et al. (1999). Joint action of two rna degradation pathways controls the timing of maternal transcript elimination at the midblastula transition in Drosophila melanogaster. EMBO J. 18, 2610–2620. doi:10.1093/emboj/18.9.2610
Batz, Z. A., Goff, A. C., and Armbruster, P. A. (2017). Micrornas are differentially abundant during aedes albopictus diapause maintenance but not diapause induction. Insect Mol. Biol. 26, 721–733. doi:10.1111/imb.12332
Bernardo, J. (1996a). Maternal effects in animal ecology. Am. Zool. 36, 83–105. doi:10.1093/icb/36.2.83
Bernardo, J. (1996b). The particular maternal effect of propagule size, especially egg size: Patterns, models, quality of evidence and interpretations. Am. Zool. 36, 216–236. doi:10.1093/icb/36.2.216
Berry, S. J. (1982). Maternal direction of oogenesis and early embryogenesis in insects. Annu. Rev. Entomol. 27, 205–227. doi:10.1146/annurev.en.27.010182.001225
Bitner, R. M. (1976). Ecological management of the alfalfa leafcutter bee, Megachile pacifica (Panzer), with emphasis on diapause induction. Logan, UT: Dissertation, Utah State University.
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: A flexible trimmer for illumina sequence data. Bioinformatics 30, 2114–2120. doi:10.1093/bioinformatics/btu170
Bonduriansky, R., and Head, M. (2007). Maternal and paternal condition effects on offspring phenotype in Telostylinus angusticollis (diptera: Neriidae). J. Evol. Biol. 20, 2379–2388. doi:10.1111/j.1420-9101.2007.01419.x
Bradshaw, W. E., and Holzapfel, C. M. (2007). Evolution of animal photoperiodism. Annu. Rev. Ecol. Evol. Syst. 38, 1–25. doi:10.1146/annurev.ecolsys.37.091305.110115
Burgess, S. C., and Marshall, D. J. (2011). Temperature-induced maternal effects and environmental predictability. J. Exp. Biol. 214, 2329–2336. doi:10.1242/jeb.054718
Bushati, N., Stark, A., Brennecke, J., and Cohen, S. M. (2008). Temporal reciprocity of mirnas and their targets during the maternal-to-zygotic transition in drosophila. Curr. Biol. 18, 501–506. doi:10.1016/j.cub.2008.02.081
Calderone, N. W. (2012). Insect pollinated crops, insect pollinators and us agriculture: Trend analysis of aggregate data for the period 1992–2009. PloS one 7, e37235. doi:10.1371/journal.pone.0037235
Christians, E., Davis, A., Thomas, S., and Benjamin, I. (2000). Maternal effect of hsf1 on reproductive success. Nature 407, 693–694. doi:10.1038/35037669
Colombo, V., Pettigrove, V. J., Golding, L. A., and Hoffmann, A. A. (2014). Transgenerational effects of parental nutritional status on offspring development time, survival, fecundity, and sensitivity to zinc in Chironomus tepperi midges. Ecotoxicol. Environ. Saf. 110, 1–7. doi:10.1016/j.ecoenv.2014.07.037
Corchete, L. A., Rojas, E. A., Alonso-López, D., De Las Rivas, J., Gutiérrez, N. C., and Burguillo, F. J. (2020). Systematic comparison and assessment of rna-seq procedures for gene expression quantitative analysis. Sci. Rep. 10, 19737–19815. doi:10.1038/s41598-020-76881-x
de Zwaan, D. R., Camfield, A. F., MacDonald, E. C., and Martin, K. (2019). Variation in offspring development is driven more by weather and maternal condition than predation risk. Funct. Ecol. 33, 447–456. doi:10.1111/1365-2435.13273
Denlinger, D. L. (2002). Regulation of diapause. Annu. Rev. Entomol. 47, 93–122. doi:10.1146/annurev.ento.47.091201.145137
Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., et al. (2013). Star: Ultrafast universal rna-seq aligner. Bioinformatics 29, 15–21. doi:10.1093/bioinformatics/bts635
Edwards, P. D., Lavergne, S. G., McCaw, L. K., Wijenayake, S., Boonstra, R., McGowan, P. O., et al. (2021). Maternal effects in mammals: Broadening our understanding of offspring programming. Front. Neuroendocrinol. 62, 100924. doi:10.1016/j.yfrne.2021.100924
Enright, A., John, B., Gaul, U., Tuschl, T., Sander, C., and Marks, D. (2003). Microrna targets in Drosophila. Genome Biol. 4, 1–27. doi:10.1186/gb-2003-5-1-r1
Evans, T. C., Crittenden, S. L., Kodoyianni, V., and Kimble, J. (1994). Translational control of maternal glp-1 mrna establishes an asymmetry in the C. elegans embryo. Cell. 77, 183–194. doi:10.1016/0092-8674(94)90311-5
Ewels, P., Magnusson, M., Lundin, S., and Käller, M. (2016). Multiqc: Summarize analysis results for multiple tools and samples in a single report. Bioinformatics 32, 3047–3048. doi:10.1093/bioinformatics/btw354
Fairey, D., and Lieverse, J. (1986). Cell production by the alfalfa leafcutting bee (megachile rotundata f.) in new and used wood and polystyrene nesting materials. J. Appl. Entomology 102, 148–153. doi:10.1111/j.1439-0418.1986.tb00905.x
Fischman, B. J., Pitts-Singer, T. L., and Robinson, G. E. (2017). Nutritional regulation of phenotypic plasticity in a solitary bee (hymenoptera: Megachilidae). Environ. Entomol. 46, 1070–1079. doi:10.1093/ee/nvx119
Friedländer, M. R., Mackowiak, S. D., Li, N., Chen, W., and Rajewsky, N. (2011). miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 40, 37–52. doi:10.1093/nar/gkr688
Gough, J., Karplus, K., Hughey, R., and Chothia, C. (2001). Assignment of homology to genome sequences using a library of hidden markov models that represent all proteins of known structure. J. Mol. Biol. 313, 903–919. doi:10.1006/jmbi.2001.5080
Griffiths-Jones, S., Grocock, R. J., Van Dongen, S., Bateman, A., and Enright, A. J. (2006). mirbase: microrna sequences, targets and gene nomenclature. Nucleic Acids Res. 34, D140–D144. doi:10.1093/nar/gkj112
Groothuis, T. G., Hsu, B.-Y., Kumar, N., and Tschirren, B. (2019). Revisiting mechanisms and functions of prenatal hormone-mediated maternal effects using avian species as a model. Philos. Trans. R. Soc. Lond. B Biol. Sci. 374, 20180115. doi:10.1098/rstb.2018.0115
Groothuis, T. G., and Schwabl, H. (2008). Hormone-mediated maternal effects in birds: Mechanisms matter but what do we know of them? Philos. Trans. R. Soc. Lond. B Biol. Sci. 363, 1647–1661. doi:10.1098/rstb.2007.0007
Hao, K., Jarwar, A. R., Ullah, H., Tu, X., Nong, X., and Zhang, Z. (2019). Transcriptome sequencing reveals potential mechanisms of the maternal effect on egg diapause induction of locusta migratoria. Int. J. Mol. Sci. 20, 1974. doi:10.3390/ijms20081974
Harvey, S. A., Sealy, I., Kettleborough, R., Fenyes, F., White, R., Stemple, D., et al. (2013). Identification of the zebrafish maternal and paternal transcriptomes. Development 140, 2703–2710. doi:10.1242/dev.095091
Hobbs, G., and Richards, K. (1976). Selection for a univoltine strain of Megachile (eutricharaea) pacifica (hymenoptera: Megachilidae). Can. Entomol. 108, 165–167. doi:10.4039/ent108165-2
Huestis, D. L., and Marshall, J. L. (2006). Interaction between maternal effects and temperature affects diapause occurrence in the cricket Allonemobius socius. Oecologia 146, 513–520. doi:10.1007/s00442-005-0232-z
Johansen, C. A., and Eves, J. D. (1973). Effects of chilling, humidity and seasonal conditions on emergence of the alfalfa leafcutting bee. Environ. Entomol. 2, 23–26. doi:10.1093/ee/2.1.23
Johnson, M. M. (2022). Maternal effects and management of alfalfa leafcutting bees. Logan, UT: Utah State University. Masters thesis.
Jones, P., Binns, D., Chang, H.-Y., Fraser, M., Li, W., McAnulla, C., et al. (2014). InterProScan 5: Genome-scale protein function classification. Bioinformatics 30, 1236–1240. doi:10.1093/bioinformatics/btu031
Kapheim, K. M., and Johnson, M. M. (2017). Juvenile hormone, but not nutrition or social cues, affects reproductive maturation in solitary alkali bees (nomia melanderi). J. Exp. Biol. 220, 3794–3801. doi:10.1242/jeb.162255
Kapheim, K. M., Jones, B. M., Søvik, E., Stolle, E., Waterhouse, R. M., Bloch, G., et al. (2020). Brain micrornas among social and solitary bees. R. Soc. Open Sci. 7, 200517. doi:10.1098/rsos.200517
Kapheim, K. M., Pan, H., Li, C., Salzberg, S. L., Puiu, D., Magoc, T., et al. (2015). Social evolution. Genomic signatures of evolutionary transitions from solitary to group living. Science 348, 1139–1143. doi:10.1126/science.aaa4788
Kemp, W., and Bosch, J. (2001). Postcocooning temperatures and diapause in the alfalfa pollinator Megachile rotundata (hymenoptera: Megachilidae). an. 94, 244–250. doi:10.1603/0013-8746(2001)094[0244:ptadit]2.0.co;2
Klostermeyer, E., and Gerber, H. S. (1969). Nesting behavior of Megachile rotundata (hymenoptera: Megachilidae) monitored with an event recorder. Ann. Entomological Soc. Am. 62, 1321–1325. doi:10.1093/aesa/62.6.1321
Klostermeyer, E., Mech, S. J., and Rasmussen, W. B. (1973). Sex and weight of Megachile rotundata (hymenoptera: Megachilidae) progeny associated with provision weights. J. Kans. Entomological Soc., 536–548.
Krüger, J., and Rehmsmeier, M. (2006). RNAhybrid: microRNA target prediction easy, fast and flexible. Nucleic Acids Res. 34, W451–W454. doi:10.1093/nar/gkl243
Krunic, M. (1972). Voltinism in Megachile rotundata (megachilidae: Hymenoptera) in southern alberts. Can. Entomol. 104, 185–188. doi:10.4039/ent104185-2
Lacour, G., Vernichon, F., Cadilhac, N., Boyer, S., Lagneau, C., and Hance, T. (2014). When mothers anticipate: Effects of the prediapause stage on embryo development time and of maternal photoperiod on eggs of a temperate and a tropical strains of Aedes albopictus (diptera: Culicidae). J. Insect Physiol. 71, 87–96. doi:10.1016/j.jinsphys.2014.10.008
Langfelder, P., and Horvath, S. (2008). Wgcna: An r package for weighted correlation network analysis. BMC Bioinforma. 9, 559–613. doi:10.1186/1471-2105-9-559
Laver, J. D., Marsolais, A. J., Smibert, C. A., and Lipshitz, H. D. (2015). Regulation and function of maternal gene products during the maternal-to-zygotic transition in Drosophila. Curr. Top. Dev. Biol. 113, 43–84. doi:10.1016/bs.ctdb.2015.06.007
Law, C. W., Alhamdoosh, M., Su, S., Dong, X., Tian, L., Smyth, G. K., et al. (2016). Rna-seq analysis is easy as 1-2-3 with limma, glimma and edger. F1000Res. 5, 1408. doi:10.12688/f1000research.9005.3
Law, C. W., Zeglinski, K., Dong, X., Alhamdoosh, M., Smyth, G. K., and Ritchie, M. E. (2020). A guide to creating design matrices for gene expression experiments. F1000Res. 9, 1444. doi:10.12688/f1000research.27893.1
Lee, I. H., and Duvall, L. B. (2022). Maternally instigated diapause in aedes albopictus: Coordinating experience and internal state for survival in variable environments. Front. Behav. Neurosci. 121, 8264. doi:10.3389/fnbeh.2022.778264
Lee, M. T., Bonneau, A. R., and Giraldez, A. J. (2014). Zygotic genome activation during the maternal-to-zygotic transition. Annu. Rev. Cell. Dev. Biol. 30, 581–613. doi:10.1146/annurev-cellbio-100913-013027
Marco, A. (2015). Selection against maternal microrna target sites in maternal transcripts. G3 5, 2199–2207. doi:10.1534/g3.115.019497
Marshall, D. J., and Uller, T. (2007). When is a maternal effect adaptive? Oikos 116, 1957–1963. –1963. doi:10.1111/j.2007.0030-1299.16203.x
McWatters, H., and Saunders, D. (1997). Inheritance of the photoperiodic response controlling larval diapause in the blow fly, Calliphora vicina. J. Insect Physiol. 43, 709–717. doi:10.1016/s0022-1910(97)00051-6
Meaney, M. J. (2001). Maternal care, gene expression, and the transmission of individual differences in stress reactivity across generations. Annu. Rev. Neurosci. 24, 1161–1192. doi:10.1146/annurev.neuro.24.1.1161
Mello, C. C., Draper, B. W., Krause, M., Weintraub, H., and Priess, J. R. (1992). The pie-1 and mex-1 genes and maternal control of blastomere identity in early C. elegans embryos. Cell. 70, 163–176. doi:10.1016/0092-8674(92)90542-k
Meylan, S., Miles, D. B., and Clobert, J. (2012). Hormonally mediated maternal effects, individual strategy and global change. Philos. Trans. R. Soc. Lond. B Biol. Sci. 367, 1647–1664. doi:10.1098/rstb.2012.0020
Miccoli, A., Dalla Valle, L., and Carnevali, O. (2017). The maternal control in the embryonic development of zebrafish. Gen. Comp. Endocrinol. 245, 55–68. doi:10.1016/j.ygcen.2016.03.028
Miccoli, A., Gioacchini, G., Maradonna, F., Benato, F., Skobo, T., and Carnevali, O. (2015). Beneficial bacteria affect Danio rerio development by the modulation of maternal factors involved in autophagic, apoptotic and dorsalizing processes. Cell. Physiol. biochem. 35, 1706–1718. doi:10.1159/000373983
Mistry, J., Chuguransky, S., Williams, L., Qureshi, M., Salazar, G., Sonnhammer, E. L. L., et al. (2020). Pfam: The protein families database in 2021. Nucleic Acids Res. 49, D412–D419. doi:10.1093/nar/gkaa913
Mitchell, S. E., and Read, A. F. (2005). Poor maternal environment enhances offspring disease resistance in an invertebrate. Proc. Biol. Sci. 272, 2601–2607. doi:10.1098/rspb.2005.3253
Mori, A., Oda, T., Wada, Y., et al. (1981). Studies on the egg diapause and overwintering of aedes albopictus in nagasaki. Trop. Med. 23, 79.
Mousseau, T. A., and Dingle, H. (1991). Maternal effects in insect life histories. Annu. Rev. Entomol. 36, 511–534. doi:10.1146/annurev.en.36.010191.002455
Mousseau, T. A., and Fox, C. W. (1998). Maternal effects as adaptations. New York, NY: Oxford University Press.
Nakahara, K., Kim, K., Sciulli, C., Dowd, S. R., Minden, J. S., and Carthew, R. W. (2005). Targets of microrna regulation in the Drosophila oocyte proteome. Proc. Natl. Acad. Sci. U. S. A. 102, 12023–12028. doi:10.1073/pnas.0500053102
Navara, K., Siefferman, L., Hill, G., and Mendonca, M. (2006). Yolk androgens vary inversely to maternal androgens in eastern bluebirds: An experimental study. Funct. Ecol. 20, 449–456. doi:10.1111/j.1365-2435.2006.01114.x
Ninova, M., Ronshaugen, M., and Griffiths-Jones, S. (2016). Microrna evolution, expression, and function during short germband development in Tribolium castaneum. Genome Res. 26, 85–96. doi:10.1101/gr.193367.115
Paranjpe, S. S., Jacobi, U. G., van Heeringen, S. J., and Veenstra, G. J. C. (2013). A genome-wide survey of maternal and embryonic transcripts during Xenopus tropicalis development. BMC genomics 14, 762–817. doi:10.1186/1471-2164-14-762
Parker, F. D., and Tepedino, V. J. (1982). Maternal influence on diapause in the alfalfa leafcutting Bee1 (hymenoptera: Megachilidae)2. Ann. Entomological Soc. Am. 75, 407–410. doi:10.1093/aesa/75.4.407
Pauli, A., Rinn, J. L., and Schier, A. F. (2011). Non-coding rnas as regulators of embryogenesis. Nat. Rev. Genet. 12, 136–149. doi:10.1038/nrg2904
Pires, C. V., Freitas, F. C. d. P., Cristino, A. S., Dearden, P. K., and Simões, Z. L. P. (2016). Transcriptome analysis of honeybee (apis mellifera) haploid and diploid embryos reveals early zygotic transcription during cleavage. PLoS One 11, e0146447. doi:10.1371/journal.pone.0146447
Pitts-Singer, T. L. (2020). Photoperiod effect on megachile rotundata (hymenoptera: Megachilidae) female regarding diapause status of progeny: The importance of data scrutiny. Environ. Entomol. 49, 516–527. doi:10.1093/ee/nvaa004
Pitts-Singer, T. L., and Cane, J. H. (2011). The alfalfa leafcutting bee, Megachile rotundata: The world’s most intensively managed solitary bee. Annu. Rev. Entomol. 56, 221–237. doi:10.1146/ANNUREV-ENTO-120709-144836
Poelchau, M. F., Reynolds, J. A., Denlinger, D. L., Elsik, C. G., and Armbruster, P. A. (2011). A de novo transcriptome of the Asian tiger mosquito, aedes albopictus, to identify candidate transcripts for diapause preparation. BMC genomics 12, 619–19. doi:10.1186/1471-2164-12-619
Poelchau, M. F., Reynolds, J. A., Elsik, C. G., Denlinger, D. L., and Armbruster, P. A. (2013). Deep sequencing reveals complex mechanisms of diapause preparation in the invasive mosquito, aedes albopictus. Proc. Biol. Sci. 280, 20130143. doi:10.1098/rspb.2013.0143
Potticary, A. L., and Duckworth, R. A. (2020). Multiple environmental stressors induce an adaptive maternal effect. Am. Nat. 196, 487–500. doi:10.1086/710210
R Core Team (2020). R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.
Rank, G., and Goerzen, D. (1982). Effect of incubation temperatures on emergence of Megachile rotundata (hymenoptera: Megachilidae). J. Econ. Entomology 75, 467–471. doi:10.1093/jee/75.3.467
Räsänen, K., and Kruuk, L. (2007). Maternal effects and evolution at ecological time-scales. Funct. Ecol. 21, 408–421. doi:10.1111/j.1365-2435.2007.01246.x
Rauwerda, H., Wackers, P., Pagano, J. F., de Jong, M., Ensink, W., Dekker, R., et al. (2016). Mother-specific signature in the maternal transcriptome composition of mature, unfertilized zebrafish eggs. PLoS One 11, e0147151. doi:10.1371/journal.pone.0147151
Richards, K. W. (1984). Alfalfa leafcutter bee management in western Canada. Ottawa: Agriculture Canada Publication 149EMinistry of Supply and Services.
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma powers differential expression analyses for rna-sequencing and microarray studies. Nucleic Acids Res. 43, e47. doi:10.1093/nar/gkv007
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edger: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi:10.1093/bioinformatics/btp616
Rolff, J. (1999). Parasitism increases offspring size in a damselfly: Experimental evidence for parasite-mediated maternal effects. Anim. Behav. 58, 1105–1108. doi:10.1006/anbe.1999.1240
Romney, A. L., and Podrabsky, J. E. (2017). Transcriptomic analysis of maternally provisioned cues for phenotypic plasticity in the annual killifish, Austrofundulus limnaeus. Evodevo 8, 6–20. doi:10.1186/s13227-017-0069-7
Rossi, B. H., Nonacs, P., and Pitts-Singer, T. L. (2010). Sexual harassment by males reduces female fecundity in the alfalfa leafcutting bee, megachile rotundata. Anim. Behav. 79, 165–171. doi:10.1016/j.anbehav.2009.10.023
Santos, P. K. F., de Souza Araujo, N., Françoso, E., Zuntini, A. R., and Arias, M. C. (2018). Diapause in a tropical oil-collecting bee: Molecular basis unveiled by rna-seq. BMC genomics 19, 305–311. doi:10.1186/s12864-018-4694-x
Saunders, D., Macpherson, J., and Cairncross, K. (1986). Maternal and larval effects of photoperiod on the induction of larval diapause in two species of fly, Calliphora vicina and Lucilia sericata. Exp. Biol. 46, 51
Sharda, S., Zuest, T., Erb, M., and Taborsky, B. (2021). Predator-induced maternal effects determine adaptive antipredator behaviors via egg composition. Proc. Natl. Acad. Sci. U. S. A. 118, e2017063118. doi:10.1073/pnas.2017063118
Spradling, A. C. (1993). Germline cysts: Communes that work. Cell. 72, 649–651. doi:10.1016/0092-8674(93)90393-5
Steiger, S. (2013). Bigger mothers are better mothers: Disentangling size-related prenatal and postnatal maternal effects. Proc. Biol. Sci. 280, 20131225. doi:10.1098/rspb.2013.1225
Sung, H.-w., Spangenberg, S., Vogt, N., and Großhans, J. (2013). Number of nuclear divisions in the Drosophila blastoderm controlled by onset of zygotic transcription. Curr. Biol. 23, 133–138. doi:10.1016/j.cub.2012.12.013
Tadros, W., and Lipshitz, H. D. (2009). The maternal-to-zygotic transition: A play in two acts. Development 136, 3033–3042. doi:10.1242/dev.033183
Tanaka, H. (1994). Embryonic diapause and life cycle in the migratory locust, locusta migratoria l.(orthoptera: Acrididae), in kyoto. Appl. Entomol. Zool. 29, 179–191. doi:10.1303/aez.29.179
Tang, F., Kaneda, M., O’Carroll, D., Hajkova, P., Barton, S. C., Sun, Y. A., et al. (2007). Maternal micrornas are essential for mouse zygotic development. Genes. Dev. 21, 644–648. doi:10.1101/gad.418707
Tepedino, V. J., and Parker, F. D. (1988). Alternation of sex ratio in a partially bivoltine bee, Megachile rotundata (hymenoptera: Megachilidae). Ann. Entomological Soc. Am. 81, 467–476. doi:10.1093/aesa/81.3.467
Tepedino, V. J., and Parker, F. D. (1986). Effect of rearing temperature on mortality, second-generation emergence, and size of adult in megachile rotundata (hymenoptera: Megachilidae). J. Econ. Entomology 79, 974–977. doi:10.1093/jee/79.4.974
Thomas, P. D., Ebert, D., Muruganujan, A., Mushayahama, T., Albou, L.-P., and Mi, H. (2022). Panther: Making genome-scale phylogenetics accessible to all. Protein Sci. 31, 8–22. doi:10.1002/pro.4218
Tian, S., and Monteiro, A. (2022). A transcriptomic atlas underlying developmental plasticity of seasonal forms of Bicyclus anynana butterflies. Mol. Biol. Evol. 39, msac126. doi:10.1093/molbev/msac126
Torchio, P. F., and Tepedino, V. J. (1980). Sex ratio, body size and seasonality in a solitary bee, Osmia lignaria propinqua cresson (hymenoptera: Megachilidae). Evolution 34, 993–1003. doi:10.2307/2408004
Torres-Paz, J., Leclercq, J., and Rétaux, S. (2019). Maternally regulated gastrulation as a source of variation contributing to cavefish forebrain evolution. Elife 8, e50160. doi:10.7554/eLife.50160
Tougeron, K., Devogel, M., van Baaren, J., Le Lann, C., and Hance, T. (2020). Trans-generational effects on diapause and life-history-traits of an aphid parasitoid. J. Insect Physiol. 121, 104001. doi:10.1016/j.jinsphys.2019.104001
Tougeron, K., Hraoui, G., Le Lann, C., van Baaren, J., and Brodeur, J. (2018). Intraspecific maternal competition induces summer diapause in insect parasitoids. Insect Sci. 25, 1080–1088. doi:10.1111/1744-7917.12491
Vastenhouw, N. L., Cao, W. X., and Lipshitz, H. D. (2019). The maternal-to-zygotic transition revisited. Development 146, dev161471. doi:10.1242/dev.161471
Wang, R. (1966). Observations on the influence of photoperiod on egg diapause in aedes albopictus skuse. Acta Entomol. Sin. 15, 75–77.
Weeks, D., and Melton, D. (1987). A maternal mrna localized to the vegetal hemisphere in Xenopus eggs codes for a growth factor related to tgf-β. Cell. 51, 861–867. doi:10.1016/0092-8674(87)90109-7
Wilson, E. S., Murphy, C. E., Wong, C., Rinehart, J. P., Yocum, G. D., and Bowsher, J. H. (2021). Environmental impacts on diapause and survival of the alfalfa leafcutting bee, Megachile rotundata. Plos one 16, e0254651. doi:10.1371/journal.pone.0254651
Winata, C. L., and Korzh, V. (2018). The translational regulation of maternal mrna s in time and space. FEBS Lett. 592, 3007–3023. doi:10.1002/1873-3468.13183
Wolf, J. B., and Wade, M. J. (2009). What are maternal effects (and what are they not)? Philos. Trans. R. Soc. Lond. B Biol. Sci. 364, 1107–1115. doi:10.1098/rstb.2008.0238
Xanthos, J. B., Kofron, M., Wylie, C., and Heasman, J. (2001). Maternal vegt is the initiator of a molecular network specifying endoderm in Xenopus laevis. Development 128, 167–180. doi:10.1242/dev.128.2.167
Yocum, G. D., Childers, A. K., Rinehart, J. P., Rajamohan, A., Pitts-Singer, T. L., Greenlee, K. J., et al. (2018). Environmental history impacts gene expression during diapause development in the alfalfa leafcutting bee, Megachile rotundata. J. Exp. Biol. 221, jeb173443. doi:10.1242/jeb.173443
Yocum, G. D., Rinehart, J. P., West, M., and Kemp, W. P. (2010). Interrupted incubation and short-term storage of the alfalfa pollinator Megachile rotundata (hymenoptera: Megachilidae): A potential tool for synchronizing bees with bloom. J. Econ. Entomol. 103, 234–241. doi:10.1603/ec09351
Keywords: maternal effects, development, maternal RNA, microRNA, oogenesis, diapause, photoperiod, maternal-to-zygotic transition
Citation: Hagadorn MA, Hunter FK, DeLory T, Johnson MM, Pitts-Singer TL and Kapheim KM (2023) Maternal body condition and season influence RNA deposition in the oocytes of alfalfa leafcutting bees (Megachile rotundata). Front. Genet. 13:1064332. doi: 10.3389/fgene.2022.1064332
Received: 08 October 2022; Accepted: 28 November 2022;
Published: 04 January 2023.
Edited by:
Hongmei Li-Byarlay, Central State University, United StatesReviewed by:
Chuan Ma, Institute of Apiculture Research (CAAS), ChinaJulie A. Reynolds, The Ohio State University, United States
Copyright © 2023 Hagadorn, Hunter, DeLory, Johnson, Pitts-Singer and Kapheim. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Karen M. Kapheim , karen.kapheim@usu.edu
†These authors contributed equally to this work and share first authorship