Tropilaelaps mercedesae parasitism changes behavior and gene expression in honey bee workers.

Tropilaelaps mercedesae is one of the most problematic honey bee parasites and has become more threatening to the beekeeping industry. Tropilaelaps can easily parasitize immature honey bees (larvae and pupae) and have both lethal and sublethal effects on the individual worker bees. Our study for the first time experimentally assessed the effects of T. mercedesae on olfactory learning, flight ability, homing ability as well as transcriptional changes in parasitized adult honey bees. T. mercedesae infestation had negative impacts on olfactory associated function, flight ability, and homing rate. The volume of the mushroom body significantly increased in infested honey bees, which may be correlated to the lower sucrose responsiveness as well as lower learning ability in the infested bees. The gene expression involved in immune systems and carbohydrate transport and metabolism were significantly different between infested bees and non-infested bees. Moreover, genes function in cell adhesion play an essential role in olfactory sensory in honey bees. Our findings provide a comprehensive understanding of European honey bees in response to T. mercedesae infestation, and could be used to further investigate the complex molecular mechanisms in honey bees under parasitic stress.

Acute and chronic effects of sublethal neonicotinoid thiacloprid to Asian honey bee (Apis cerana cerana).

Pesticide pollution is one of the most important factors for global bee declines. Despite many studies have revealed that the most important Chinese indigenous species，Apis cerana, is presenting a high risk on exposure to neonicotinoids, the toxicology information on Apis cerana remain limited. This study was aimed to determine the acute and chronic toxic effects of thiacloprid (IUPAC name: {(2Z)-3-[(6-Chloro-3-pyridinyl)methyl]-1,3-thiazolidin-2-ylidene}cyanamide) on behavioral and physiological performance as well as genome-wide transcriptome in A. cerana. We found the 1/5 LC50 of thiacloprid significantly impaired learning and memory abilities after both acute and chronic exposure, nevertheless, has no effects on the sucrose responsiveness and phototaxis climbing ability of A. cerana. Moreover, activities of detoxification enzyme P450 monooxygenases and CarE were increased by short-term exposure to thiacloprid, while prolonged exposure caused suppression of CarE activity. Neither acute nor chronic exposure to thiacloprid altered honey bee AChE activities. To further study the potential defense molecular mechanisms in Asian honey bee under pesticide stress, we analyzed the transcriptomes of honeybees in response to thiacloprid stress. The transcriptomic profiles revealed consistent upregulation of immune- and stress-related genes by both acute or chronic treatments. Our results suggest that the chronic exposure to thiacloprid produced greater toxic effects than a single administration to A. cerana. Altogether, our study deepens the understanding of the toxicological characteristic of A. cerana against thiacloprid, and could be used to further investigate the complex molecular mechanisms in Asian honey bee under pesticide stress.

Microplastic Polystyrene Ingestion Promotes the Susceptibility of Honeybee to Viral Infection.

Microplastics (MPs) are an emerging threat to ecological conservation and biodiversity; however, little is known of the types and possible impacts of MPs in pollinators. To examine whether MPs were present in honeybees, we analyzed the honeybee samples collected in fields from six provinces in China. Four types MPs were identified in honeybee including polystyrene (PS) by Raman spectroscopic analysis, and these plastic polymers were detected in 66.7% bee samples. Then, we assessed the physical and biological impacts of PS of three sizes (0.5, 5, and 50 µm) on bees for 21 days. Next, we measured how the presence of PS affected the Israeli acute paralysis virus proliferation, a small RNA virus associated with bee colony decline. Experimental evidence showed that a large mass of PS was ingested and accumulated within the midgut and enhanced the susceptibility of bees to viral infection. Not only histological analysis showed that PS, especially 0.5 µm PS, damaged the midgut tissue and was subsequently transferred to the hemolymph, trachea, and Malpighian tubules, but also qPCR and transcriptomic results indicated that genes correlated with membrane lipid metabolism, immune response, detoxification, and the respiratory system were significantly regulated after PS ingestion. Our results highlight neglected MP contamination to the bees, a pollination ecosystem stressed by the anthropogenic pollution, and have implications for human health via ingestion of bee products.

Flupyradifurone reduces nectar consumption and foraging but does not alter honey bee recruitment dancing.

Foraging is essential for honey bee colony fitness and is enhanced by the waggle dance, a recruitment behavior in which bees can communicate food location and quality. We tested if the consumption of nectar (sucrose solution) with a field-realistic concentration of 4 ppm flupyradifurone (FPF) could alter foraging behavior and recruitment dancing in Apis mellifera. Foragers were repelled by FPF. They visited the FPF feeder less often and spent less time imbibing sucrose solution (2.5 M, 65% w/w) with FPF. As a result, bees feeding on the FPF treatment consumed 16% less nectar. However, FPF did not affect dancing: there were no effects on unloading wait time, the number of dance bouts per nest visit, or the number of dance circuits performed per dance bout. FPF could therefore deter bees from foraging on contaminated nectar. However, the willingness of bees to recruit nestmates for nectar with FPF is concerning. Recruitment can rapidly amplify the number of foragers and could overcome the decrease in consumption of FPF-contaminated nectar, resulting in a net inflow of pesticide to the colony. FPF also significantly altered the expression of 116 genes, some of which may be relevant for the olfactory learning deficits induced by FPF and the toxicity of FPF.

IAPV-Induced Paralytic Symptoms Associated with Tachypnea via Impaired Tracheal System Function.

Although it had been reported that Israeli acute paralysis virus (IAPV) can cause systemic infection in honey bees, little is known about how it establishes this infection and results in the typical symptoms, paralysis and trembling. Here, we used our previously constructed IAPV infectious clone to investigate viral loads in different tissues of honey bees and further identify the relation between tissue tropism and paralytic symptoms. Our results showed that tracheae showed a greater concentration of viral abundance than other tissues. The abundance of viral protein in the tracheae was positively associated with viral titers, and was further confirmed by immunological and ultrastructural evidence. Furthermore, higher viral loads in tracheae induced remarkable down-regulation of succinate dehydrogenase and cytochrome c oxidase genes, and progressed to causing respiratory failure of honey bees, resulting in the appearance of typical symptoms, paralysis and body trembling. Our results showed that paralysis symptoms or trembling was actually to mitigate tachypnea induced by IAPV infection due to the impairment of honey bee tracheae, and revealed a direct causal link between paralysis symptoms and tissue tropism. These findings provide new insights into the understanding of the underlying mechanism of paralysis symptoms of honey bees after viral infection and have implications for viral disease prevention and specific therapeutics in practice.

Systematic identification of circular RNAs and corresponding regulatory networks unveil their potential roles in the midguts of eastern honeybee workers.

Currently, knowledge of circular RNAs (circRNAs) in insects including honeybee is extremely limited. Here, differential expression profiles and regulatory networks of circRNAs in the midguts of Apis cerana cerana workers were comprehensively investigated using transcriptome sequencing and bioinformatics. In total, 9589 circRNAs (201-800 nt in length) were identified from 8-day-old and 11-day-old workers' midguts (Ac1 and Ac2); among them, 5916 (61.70%) A. cerana cerana circRNAs showed conservation with our previously indentified circRNAs in Apis mellifera ligucstica workers' midguts (Xiong et al., Acta Entomologica Sinica 61:1363-1375, 2018). Five circRNAs were confirmed by RT-PCR and Sanger sequencing. Interestingly, novel_circ_003723, novel_circ_002714, novel_circ_002451, and novel_circ_001980 were highly expressed in both Ac1 and Ac2. In addition, the source genes of circRNAs were involved in 34 GO terms including organelle and cellular process and 141 pathways such as endocytosis and Wnt signaling pathway. Moreover, 55 DEcircRNAs including 34 upregulated and 21 downregulated circRNAs were identified in Ac2 compared with Ac1. circRNA-miRNA regulatory networks indicated that 1060 circRNAs can target 74 miRNAs; additionally, the DEcircRNA-miRNA-mRNA networks suggested that 13 downregulated circRNAs can bind to eight miRNAs and 29 miRNA-targeted mRNAs, while 16 upregulated circRNAs can link to 9 miRNAs and 29 miRNA-targeted mRNAs. These results indicated that DEcircRNAs as ceRNAs may play a comprehensive role in the growth, development, and metabolism of the worker's midgut via regulating source genes and interacting with miRNAs. Notably, eight DEcircRNAs targeting miR-6001-y were likely to be key participants in the midgut development. Our findings not only offer a valuable resource for further studies on A. cerana cerana circRNA and novel insights into understanding the molecular mechanisms underlying the midgut development of eastern honeybee but also provide putative circRNA candidates for functional research in the near future and novel biomarkers for identification of eastern honeybee species including A. cerana cerana and honeybee diseases such as chalkbrood and microsporidiosis.

A Reverse Genetics System for the Israeli Acute Paralysis Virus and Chronic Bee Paralysis Virus.

Honey bee viruses are associated with honey bee colony decline. Israeli acute paralysis virus (IAPV) is considered to have a strong impact on honey bee survival. Phylogenetic analysis of the viral genomes from several regions of the world showed that various IAPV lineages had substantial differences in virulence. Chronic bee paralysis virus (CBPV), another important honey bee virus, can induce two significantly different symptoms. However, the infection characteristics and pathogenesis of IAPV and CBPV have not been completely elucidated. Here, we constructed infectious clones of IAPV and CBPV using a universal vector to provide a basis for studying their replication and pathogenesis. Infectious IAPV and CBPV were rescued from molecular clones of IAPV and CBPV genomes, respectively, that induced typical paralysis symptoms. The replication levels and expression proteins of IAPV and CBPV in progeny virus production were confirmed by qPCR and Western blot. Our results will allow further dissection of the role of each gene in the context of viral infection while helping to study viral pathogenesis and develop antiviral drugs using reverse genetics systems.

Transcriptomic investigation of immune responses of the Apis cerana cerana larval gut infected by Ascosphaera apis.

Chalkbrood is the most common fungal disease in honeybees. The objective of this study was to reveal immune responses in the Apis cerana cerana larval gut following Ascosphaera apis invasion. Combining a previously assembled transcriptome of A. c. cerana larval gut and the high-throughput sequencing data obtained in this study, 6152 differentially expressed genes (DEGs) were clustered into eight profiles. Trend analysis showed three significant up-regulated profiles (p ≤ 0.05) and three down-regulated profiles. Gene Ontology (GO) term analysis suggested that DEGs within significant up-regulated and down-regulated clusters were enriched in 46 and 38 functional groups, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis indicated a majority of DEGs were involved in ribosome structure or function, carbon metabolism, biosynthesis of amino acids, and oxidative phosphorylation. In addition, 142 and 14 DEGs were annotated in the cellular immune- and humoral immune-related pathways, respectively. Further investigation indicated that DEGs up-regulated in cellular immune and humoral immune pathways outnumbered those that were down-regulated. Moreover, immune responses of A. c. cerana and Apis mellifera ligustica larvae were compared and studied to decipher resistance of eastern honeybee larvae to A. apis. These results demonstrated that a large number of genes involved in immunity-related pathways were activated by A. apis. Our findings provided valuable information for elucidating the molecular mechanisms underlying immune responses of A. c. cerana larvae to A. apis infection and pathogen-host interactions during chalkbrood infection.

The effect of Bt Cry9Ee toxin on honey bee brood and adults reared in vitro, Apis mellifera (Hymenoptera: Apidae).

The effects of Bt Cry9Ee toxin on honey bee, Apis mellifera L., survival, developmental rate, larval weight, pollen consumption, and midgut bacterial diversity were tested in the laboratory. Honey bee larvae and adults were reared in vitro and fed a diet that contained Cry9Ee toxin at 0.01, 0.1, 1, and 10 mg/L. Cry9Ee toxin 0.01, 0.1, and 1 mg/L in diet used in this study may represent a value closer to field relevance and the highest concentration is unlikely to be encountered in the field and thus represent a worst case scenario. The dependent variables were compared for groups of honey bees feeding on treated diet and those feeding on negative control (no addition of a test substance), solvent control (0.01 mM Na2CO3), and positive control diet (dimethoate 45 mg/L). Bt Cry9Ee toxin did not affect survival or larval weight, and the result was great confidence in accepting the null hypothesis by power analysis. The effect on development rates and pollen consumption were the inconclusive results because the post-hoc power was less than 0.8. Furthermore, the midgut bacterial structure and compositions were determined using high-throughput sequencing targeting the V3-V4 regions of the 16S rDNA. All core honey bee intestinal bacterial class such as γ-Proteobacteria, Actinobacteria, α-Proteobacteria, Bacilli, ß-Proteobacteria, and Bacteroidia were detected, and no significant changes were found in the species diversity and richness between Cry9Ee treatments and laboratory control.

Effects of Tropilaelaps mercedesae on midgut bacterial diversity of Apis mellifera.

Tropilaelaps mercedesae is an ectoparasite of Apis mellifera in Asia and is considered a major threat to honey bee health. Herein, we used the Illumina MiSeq platform 16S rDNA Amplicon Sequencing targeting the V3-V4 regions and analysed the effects on the midgut bacterial communities of honey bees infested with T. mercedesae. The overall bacterial community in honey bees infested with T. mercedesae were observed at different developmental stages. Honey bee core intestinal bacterial genera such as Gilliamella, Lactobacillus and Frischella were detected. Tropilaelapsmercedesae infestation changed the bacterial communities in the midgut of A. mellifera. Tropilaelapsmercedesae-infested pupae had greatly increased relative abundances of Micrococcus and Sphingomonas, whereas T. mercedesae-infested 15-day-old workers had significantly reduced relative abundance of non-core microbes: Corynebacterium, Sphingomonas, Acinetobacter and Enhydrobacter compared to T. mercedesae-infested newly emerged bees. The bacterial community was significantly changed at the various T. mercedesae-infested developmental stages of A. mellifera. Tropilaelapsmercedesae infestation also changed the non-core bacterial community from larvae to newly emerged honey bees. Bacterial communities were significantly different between T. mercedesa-infested and non-mite-infested 15-day-old workers. Lactobacillus was dominant in T. mercedesae-infested 15-day-old workers compared to non-mite-infested 15-day-old workers.

Identification of long non-coding RNAs in the chalkbrood disease pathogen Ascospheara apis.

Ascospheara apis is a widespread fungal pathogen that exclusively invades honeybee larvae. Thus far, non-coding RNA in A. apis has not yet been documented. In this study, we sequenced A. apis using strand specific cDNA library construction and Illumina RNA sequencing methods, and identified 379 lncRNAs, including antisense lncRNAs, lincRNAs, intronic lncRNAs and sense lncRNAs. Additionally, these lncRNAs were found to be shorter in length and have fewer exons and transcript isoforms than protein-coding genes, similar to those identified in mammals and plants. Furthermore, the existence of 15 predicted lncRNAs of A. apis was confirmed using RT-PCR and expression levels of 11 were lower than those of adjacent protein-coding genes. Our findings not only enlarge the lncRNA database for fungi, but also lay a foundation for further investigation of potential lncRNA-mediated regulation of genes in A. apis.

The discovery of Halictivirus resolves the Sinaivirus phylogeny.

By providing pollination services, bees are among the most important insects, both in ecological and economical terms. Combined next-generation and classical sequencing approaches were applied to discover and study new insect viruses potentially harmful to bees. A bioinformatics virus discovery pipeline was used on individual Illumina transcriptomes of 13 wild bees from three species from the genus Halictus and 30 ants from six species of the genera Messor and Aphaenogaster. This allowed the discovery and description of three sequences of a new virus termed Halictus scabiosae Adlikon virus (HsAV). Phylogenetic analyses of ORF1, RNA-dependent RNA-polymerase (RdRp) and capsid genes showed that HsAV is closely related to (+)ssRNA viruses of the unassigned Sinaivirus genus but distant enough to belong to a different new genus we called Halictivirus. In addition, our study of ant transcriptomes revealed the first four sinaivirus sequences from ants (Messor barbarus, M. capitatus and M. concolor). Maximum likelihood phylogenetic analyses were performed on a 594 nt fragment of the ORF1/RdRp region from 84 sinaivirus sequences, including 31 new Lake Sinai viruses (LSVs) from honey bees collected in five countries across the globe and the four ant viral sequences. The phylogeny revealed four main clades potentially representing different viral species infecting honey bees. Moreover, the ant viruses belonged to the LSV4 clade, suggesting a possible cross-species transmission between bees and ants. Lastly, wide honey bee screening showed that all four LSV clades have worldwide distributions with no obvious geographical segregation.

Complete sequences of the RNA 1 and RNA 2 segments of chronic bee paralysis virus strain CBPV-BJ detected in China.

Chronic bee paralysis virus (CBPV) is an important viral pathogen that affects adult bees. Although several CBPV strains have been reported, little information has been obtained from China. In this study, two major segments of the CBPV Chinese isolate CBPV-BJ, RNA 1 and RNA 2, were determined to be 3657 and 2267 nucleotides (nt) in length, respectively. RNA 1 and RNA 2 contained three and four open reading frames (ORFs), respectively, which agreed with known reference strains (EU122229 and EU122230). The RNA 1 had 98% nucleotide sequence identity to a known Chinese strain (KU950353), and RNA 2 had 97% nucleotide sequence identity to another Chinese strain (KU950354). Although the lengths of the RNA 1 and RNA 2 sequences were 17 nt and 38 nt shorter than those of the CBPV reference strains EU122229 and EU122230, respectively, the complete CBPV-BJ RNA 1 and RNA 2 sequences shared 91% and 92% identity with them. Phylogenetic analysis based on the sequences of the RNA-dependent RNA polymerase (RdRp) and putative structural proteins (pSPs) showed that CBPV-BJ was most closely related to the other two Chinese isolate (KU950353 and KU950354) and clustered with most Asian strains. These data provide new information that will lead to a better understanding of the diversity of the CBPV genome.

Current knowledge of detoxification mechanisms of xenobiotic in honey bees.

The western honey bee Apis mellifera is the most important managed pollinator species in the world. Multiple factors have been implicated as potential causes or factors contributing to colony collapse disorder, including honey bee pathogens and nutritional deficiencies as well as exposure to pesticides. Honey bees' genome is characterized by a paucity of genes associated with detoxification, which makes them vulnerable to specific pesticides, especially to combinations of pesticides in real field environments. Many studies have investigated the mechanisms involved in detoxification of xenobiotics/pesticides in honey bees, from primal enzyme assays or toxicity bioassays to characterization of transcript gene expression and protein expression in response to xenobiotics/insecticides by using a global transcriptomic or proteomic approach, and even to functional characterizations. The global transcriptomic and proteomic approach allowed us to learn that detoxification mechanisms in honey bees involve multiple genes and pathways along with changes in energy metabolism and cellular stress response. P450 genes, is highly implicated in the direct detoxification of xenobiotics/insecticides in honey bees and their expression can be regulated by honey/pollen constitutes, resulting in the tolerance of honey bees to other xenobiotics or insecticides. P450s is also a key detoxification enzyme that mediate synergism interaction between acaricides/insecticides and fungicides through inhibition P450 activity by fungicides or competition for detoxification enzymes between acaricides. With the wide use of insecticides in agriculture, understanding the detoxification mechanism of insecticides in honey bees and how honeybees fight with the xenobiotis or insecticides to survive in the changing environment will finally benefit honeybees' management.

Exploring poisonous mechanism of honeybee, Apis mellifera ligustica Spinola, caused by pyrethroids.

As the important intracellular secondary messengers, calcium channel is the target of many neurotoxic pesticides as calcium homeostasis in the neuroplasm play important role in neuronal functions and behavior in insects. This study investigated the effect of deltamethrin (DM) on calcium channel in the brain nerve cells of adult workers of Apis mellifera ligustica Spinola that were cultured in vitro. The results showed that the intracellular calcium concentration was significantly elevated even with a very low concentration of the DM (3.125×10-2mg/L). Further testing revealed that T-type voltage-gated calcium channels (VGCCs), except for sodium channels, was one of the target of DM on toxicity of Apis mellifera, while DM has no significant effect on the L-type VGCCs, N-methyl-d-aspartate receptor-gated calcium channels and calcium store. These results suggesting that the DM may act on T-type VGCCs in brain cells of honeybees and result in behavioral abnormalities including swarming, feeding, learning, and acquisition.

First detection of Apis mellifera filamentous virus in Apis cerana cerana in China.

Although many honey bee RNA viruses have been correlated with colony declines, little is known regarding the potential role of DNA viruses. Here, we examined seemingly healthy and crawling bee samples from China using PCR to identify whether Apis mellifera filamentous virus (AmFV) was present in A. cerana cerana. The highest AmFV infection percentage among Chinese provinces occurred in crawling bees from Gansu province (85.48%), and the lowest was in bees from Beijing (31.58%). A phylogenetic analysis showed that the Chinese isolate of AmFV exhibited a high genetic similarity with isolates from Belgium, Switzerland and USA. This is the first report of AmFV infections in Chinese A. cerana cerana populations.

Simultaneous determination of trace migration of phthalate esters in honey and royal jelly by GC-MS.

A simple, rapid, and reliable liquid-liquid extraction coupled to GC-MS method was developed and validated for the quantification of 22 phthalate esters (PAEs) in honey and royal jelly. Instrument parameters for GC-MS were tested to obtain the satisfactory separation between 22 PAEs with high sensitivity. The extraction procedure was optimized in order to achieve the best recovery. The following criteria were used to validate the developed method: linearity, LOD, lower LOQ, precision, accuracy, matrix effect and carry-over. Correlation coefficients were >0.999 by applying the linear regression model based on the least-squares method with a weighting factor (1/x). The intra- and interday precision were within 12.7% in terms of RSD, and the accuracy was within -11.8% in terms of relative error. The mean extraction recoveries ranged between 80.1 and 110.9% for honey and royal jelly. No significant matrix effect and carry-over for PAEs were observed for the analysis of honey and royal jelly samples. A total of 20 real samples were analyzed for a mini-survey using the developed method. Seven PAEs in honey samples and five PAEs in royal jelly samples were found, indicating potential contamination with several PAEs.

Thiacloprid impairs honeybee worker learning and memory with inducing neuronal apoptosis and downregulating memory-related genes.

Neonicotinoids are among the most widely used insecticides in the world and are recognized as a potential cause of pollinator decline. Previous studies have demonstrated that the neonicotinoid thiacloprid has adverse effects on foraging and memory behaviors. However, there is no direct evidence linking thiacloprid-induced neuronal cell damage in the brains of honeybees to learning and memory dysfunction. Adult honeybee (Apis mellifera L.) workers were chronically exposed to sub-lethal concentrations of thiacloprid. We discovered that thiacloprid negatively affected their survival, food consumption, and body weight. In addition, sucrose sensitivity and memory performance were impaired. We evaluated the apoptosis of honeybee brain cells using TUNEL (Terminal deoxynucleotidyl transferase-mediated digoxigenin-dUTP-biotin nick-end labeling) and Caspase-3 assays, which revealed that thiacloprid increases the dose-dependent apoptosis of neurons in the mushroom bodies (MB) and antennal lobes (AL). We also determined the abnormal transcripts of multiple genes, including vitellogenin (Vg), immune system genes (apidaecin and catalase), and memory-associated genes (pka, creb, Nmdar1, Dop2, Oa1, Oa-2R, and Oa-3R). These results indicate that exposure to sublethal concentrations of thiacloprid cause abnormal expression of memory-related genes and apoptosis of brain cells in the AL and MB, which may contribute to the memory disorder induced by thiacloprid exposure.

Exposure of chlorothalonil and acetamiprid reduce the survival and cause multiple internal disturbances in Apis mellifera larvae reared in vitro.

Background: Chlorothalonil and acetamiprid are chemical pesticides commonly used in agricultural production and have been shown to have negative effects on bee's fitness. Despite many studies have revealed that honey bee (Apis mellifera L.) larvae are posting a high risk on exposure to pesticides, but the toxicology information of chlorothalonil and acetamiprid on bee larvae remain limited. Results: The no observed adverse effect concentration (NOAEC) of chlorothalonil and acetamiprid for honey bee larvae were 4 µg/mL and 2 µg/mL, respectively. Except for CarE, the enzymic activities of GST and P450 were not influenced by chlorothalonil at NOAEC, while chronic exposure to acetamiprid slightly increased the activities of the three tested enzymes at NOAEC. Further, the exposed larvae showed significantly higher expression of genes involved in a series of different toxicologically relevant process following, including caste development (Tor (GB44905), InR-2 (GB55425), Hr4 (GB47037), Ac3 (GB11637) and ILP-2 (GB10174)), immune system response (abaecin (GB18323), defensin-1 (GB19392), toll-X4 (GB50418)), and oxidative stress response (P450, GSH, GST, CarE). Conclusion: Our results suggest that the exposure to chlorothalonil and acetamiprid, even at concentrations below the NOAEC, showed potentially effects on bee larvae's fitness, and more important synergistic and behavioral effects that can affect larvae fitness should be explored in the further.