Kristin R. Duffield1*
John Hunt2
Ben M. Sadd3
Scott K. Sakaluk3Brenda Oppert4Karyna Rosario5
Robert W. Behle1José L. Ramirez1- 1Crop BioProtection Research Unit, Agricultural Research Service, United States Department of Agriculture, National Center for Agricultural Utilization Research, Peoria, IL, United States
- 2School of Science, Western Sydney University, Richmond, NSW, Australia
- 3School of Biological Sciences, Illinois State University, Normal, IL, United States
- 4Stored Product Insect and Engineering Research Unit, Agricultural Research Service, United States Department of Agriculture, Center for Grain and Animal Health Research, Manhattan, KS, United States
- 5Marine Genomics Laboratory, University of South Florida, St. Petersburg, FL, United States
Interest in developing food, feed, and other useful products from farmed insects has gained remarkable momentum in the past decade. Crickets are an especially popular group of farmed insects due to their nutritional quality, ease of rearing, and utility. However, production of crickets as an emerging commodity has been severely impacted by entomopathogenic infections, about which we know little. Here, we identified and characterized an unknown entomopathogen causing mass mortality in a lab-reared population of Gryllodes sigillatus crickets, a species used as an alternative to the popular Acheta domesticus due to its claimed tolerance to prevalent entomopathogenic viruses. Microdissection of sick and healthy crickets coupled with metagenomics-based identification and real-time qPCR viral quantification indicated high levels of cricket iridovirus (CrIV) in a symptomatic population, and evidence of covert CrIV infections in a healthy population. Our study also identified covert infections of Acheta domesticus densovirus (AdDNV) in both populations of G. sigillatus. These results add to the foundational research needed to better understand the pathology of mass-reared insects and ultimately develop the prevention, mitigation, and intervention strategies needed for economical production of insects as a commodity.
Introduction
Insect production is a rapidly growing industry globally. While the practice of farming insects has been around for millennia (e.g., silkworm farming and apiculture) (Defoliart, 1995), applications for mass-produced insects continue to expand beyond traditional uses (Castro-López et al., 2020; van Huis, 2020b), to include chitin production (Hahn et al., 2020), waste management and valorization (Surendra et al., 2016, 2020; Gasco et al., 2020), and use as feed for both pets, including cats and dogs (Bosch et al., 2014), and agricultural animals (Makkar et al., 2014; Henry et al., 2015; Tomberlin et al., 2015). Moreover, edible insects reared for direct human consumption (e.g., mealworms and crickets) are becoming an increasingly attractive solution to address the world’s critical need for novel and environmentally sustainable protein sources (van Huis et al., 2013; Hawkey et al., 2021). Worldwide, farmed insects have the potential to become a critically important commodity, buffering against food insecurities, providing additional revenue streams for rural and urban farmers, and offering sustainably produced resources in a growing number of applications (van Huis et al., 2013). Crickets (family: Gryllidae) are an especially popular group of insects due to their ease of rearing and nutritional profile (Wang et al., 2004; Zielińska et al., 2015; Stull et al., 2018) and are already used in many foods, including nutritional and functional additives (Hall et al., 2017; Osimani et al., 2018; Udomsil et al., 2019), and feed applications (van Huis, 2020a; Magara et al., 2021).
The demand for mass-produced insects is steadily increasing due to their utility with relatively low associated costs (Wilkie, 2018); however, there are several obstacles that hinder farmed insects from becoming an extensively utilized resource, including a dearth of rigorous empirical data (van Huis, 2017; Stull and Patz, 2020). Critical among these gaps is a lack of research on the entomopathogenic microbes that negatively impact reared insect colony health and production despite infectious disease outbreaks plaguing modern insect farms for decades (Eilenberg et al., 2018; Maciel-Vergara et al., 2021). As with traditional animal livestock, farmed insects are susceptible to parasites and pathogens, including viruses, bacteria, fungi, microsporidia, and nematodes (Kaya and Vega, 2012; Eilenberg et al., 2015). Pathogenic viruses can devastate reared insect populations (Maciel-Vergara and Ros, 2017) and especially so within cricket colonies, which are often reared in environments ideal for virus transmission (i.e., crowded, humid, and warm). For example, the Acheta domesticus densovirus (AdDNV), a small parvovirus, is responsible for causing severe epizootics throughout European, North American, and Asian cricket farms, resulting in massive mortality and product losses of the commonly utilized house cricket, A. domesticus (Styer and Hamm, 1991; Liu et al., 2011; Szelei et al., 2011; Weissman et al., 2012; Pham et al., 2013b). As a direct response to these outbreaks, many producers switched to farming alternative species, including Gryllodes sigillatus in North America due to reports that they were less susceptible to AdDNV (Weissman et al., 2012).
Beyond a single study (Weissman et al., 2012), we know little about viral disease and susceptibility in G. sigillatus, yet the list of entomopathogenic viruses infecting other farmed cricket species continues to grow to include iridoviruses, nudiviruses, and other dicistroviruses. Moreover, because there have not been wide-scale systematic surveillance efforts (de Miranda et al., 2021b), it is very likely that there are other pathogenic viruses among reared crickets still to be discovered. Here, we provide one of the first empirical reports of the identification and characterization of a viral pathogen in a diseased colony of lab-reared G. sigillatus crickets. We screened diseased colonies for known cricket viruses by measuring viral prevalence across sexes using both real-time PCR (qPCR) as well as Sanger sequencing and contrasted these results with a related population with no apparent signs of infection. Furthermore, we screened for novel pathogens by performing a non-targeted shotgun metagenomic analysis on guts and hemolymph from adult crickets. Metagenomic data corroborate our qPCR results and provide a complete genome of the suspected disease-causing agent, cricket iridovirus (CrIV) with significant similarities to the previously reported lizard–cricket iridovirus (Liz–CrIV) (Papp and Marschang, 2019). Our results add to the critical, but presently scarce, research on insect pathology of farmed crickets and highlight the importance of understanding viral infection and transmission dynamics in reared insect colonies.
Materials and Methods
Cricket Colonies
Experimental G. sigillatus crickets came from either of two populations (“Diseased”: an apparently diseased population, or “Healthy”: an apparently disease-free population) of lab-reared colonies (20 individuals of each sex within each population; 80 crickets total). These populations were descendants from the same ancestral wild-caught crickets collected from Las Cruces, New Mexico (United States) and have been cultured in a lab setting since 2001. Populations were split and maintained in separate labs since 2007. Symptoms present in the Diseased colony were high, intermittent mortality among late-instar nymphs and adults, a strong putrid odor within rearing containers, milky white hemolymph which appeared iridescent under illuminated magnification, increased cuticle and tissue frailty, and underdeveloped or absent ovaries in some adult females (Figure 1). Rearing methods followed standard cricket rearing protocol within a research laboratory setting (Duffield et al., 2019). Briefly, about 500 crickets were housed in 55 L plastic storage bins with ventilated lids packed with egg carton to increase rearing surface area. They were provisioned with a standard diet (roughly equal parts Mazuri® Rat and Mouse Diets and Purina® Cat Chow Complete pellets) and water (glass vials plugged with moist cotton) ad libitum. All individuals were housed in an environmental chamber at 32°C on a 16 h:8 h light:dark cycle. Experimental individuals were at least 1-week-old post-eclosion when they were killed by freezing at −80°C.
Figure 1. Micro- and macroscopic images of pathology discovered in adult Gryllodes sigillatus. (A) Ventral view of crickets from the Healthy population with no apparent signs of infection (from left to right: female, male, female). (B) Ventral view of crickets from the Diseased population with clear signs of infection (from left to right: female, male, female). (C) Ventral view of hemocoel of male crickets from the Healthy (left) and Diseased (right) populations. (D) Dissected guts from a female from the Diseased population. Note the presence of nodules (black arrow) as well as the characteristic iridescent sheen (white arrows).
DNA, RNA Extraction, and cDNA Synthesis
To screen for DNA and RNA viruses, we extracted both DNA and RNA from whole body homogenates. Previously frozen (−80°C) crickets were placed individually in tubes with 1 mL sterile 1x PBS (pH 7.2) and two 3.2 mm diameter sterile stainless-steel beads and macerated using a TissueLyser II (Qiagen, Germany). The resulting liquid homogenate was removed (about 0.9 mL) and placed in a new sterile tube for DNA and RNA extraction.
DNA was extracted from 200 μL of cricket homogenate using the DNeasy Blood and Tissue kit (Qiagen) following the “Purification of total DNA from insects” protocol. We adapted these methods to include an addition of 200 μL of ATL buffer to the homogenate with 20 μL proteinase K and kept this mixture at 56°C overnight. RNA was extracted from 100 μL of cricket homogenate using the RNeasy Mini prep kit (Qiagen) following the “Purification of Total RNA from Animal Tissues” protocol. DNA and RNA were estimated via a NanoDrop OneC Microvolume UV-Vis Spectrophotometer (Thermo Fisher Scientific). All samples were stored at −20°C until further use.
Real-Time PCR (qPCR) Detection and Quantification
For quantification purposes, we designed a primer targeting a tubulin-like reference gene (“tubu3,” GS-tubu-F3 5′-TGCGAGATCGTATTCCGTGG-3′ and GS-tubu-R3 5′-ACCTCGGGAGAGTCAATCCA-3′, amplicon size = 137 bp) using Primer-BLAST (NCBI) and used this as our reference gene target throughout (all primers from IDT, Inc., United States). Prior to conducting qPCR assays, we normalized DNA samples to 100 ng/μL, based on NanoDrop estimates. RNA was normalized to 1 μg, treated with DNA Wipeout, and then converted to cDNA using the QuantiTect Reverse Transcription Kit (Qiagen).
Viral screening targeted 8 known viruses identified in reared cricket populations, based on primary literature searches (Table 1): cricket iridovirus (CrIV) (Jakob et al., 2002; Papp et al., 2014), A. domesticus densovirus (AdDNV) (Szelei et al., 2011), A. domesticus mini ambidensovirus (AdMADV) (Pham et al., 2013c), Gryllus bimaculatus nudivirus (GbNV) (Huger, 1985; Wang et al., 2007), A. domesticus volvovirus (AdVVV) (Pham et al., 2013a), cricket paralysis virus (CrPV) (Wang et al., 2019), A. domesticus iflavirus (AdIV) (de Miranda et al., 2021b), and A. domesticus virus (AdV) (Valles and Chen, 2006). While most of these viruses are associated with high mortality, the pathology associated with AdMADV, AdVVV, and AdIV is currently unknown (Fernandez-Cassi et al., 2019; de Miranda et al., 2021b).
Table 1. Real-time PCR (qPCR) and Sanger sequencing (in gray) primers used to target known viruses infecting rearing crickets (family: Gryllidae) in this study.
Absolute quantification was performed using known standards via double-stranded DNA (dsDNA) fragments (IDT, Inc., gBlocks Gene Fragments) designed to target the genes of interest (Table 1). The absolute abundance was expressed as the ratio of viral copies to cricket genome (tubu3) copies. For single stranded viruses, calculated copies were divided by 2. Real-time qPCR reactions were run on a Quant-Studio 6 Real-Time PCR instrument (Thermo Fisher Scientific, United States), and included a melt-curve stage to confirm product specificity. The identity of the viruses detected in our assays was confirmed via Basic Local Alignment Search Tool (BLAST, NCBI) following traditional PCR and Sanger sequencing (Table 1). One microliter of DNA or cDNA product was used in a 10 μL qPCR reaction using gene specific primers (Table 1) and PowerUp SYBR green Master mix kit (Qiagen). qPCR cycling conditions consisted of holding at 50°C for 2 min and 95°C for 2 min and 40 cycles of 1 s at 95°C and 30 min at 60°C. Standard curve efficiencies were 99.03% (R2 = 0.9989) for AdDNV and 91.31% (R2 = 0.9985) for CrIV. The limit of quantification for each qPCR assay was 1 copy/μL for both AdDNV and CrIV.
Metagenomic Analysis
To screen for novel viruses, we performed non-targeted shotgun metagenomic sequencing of four adult Diseased crickets. DNA extracts from the hemolymph of each specimen and two dissected guts were used for library preparation and multiplexing using the Illumina Nextera DNA Flex kit. Briefly, libraries were prepared using the Nextera DNA flex library standard protocol according to the manufacturer’s instructions (Illumina, Nextera DNA Flex Library Prep Reference guide). Each sample was diluted to 240 ng of gDNA in 30 μL of nuclease-free water for the Tagmentation reaction, using 10 μL of Bead-Linked Transposase (BLT) and 10 μL of TB1 solution with an incubation period of 15 min at 55°C and then held at 10°C. Following Tagmentation, the DNA-BLT complex was washed three times using Tagmentation wash buffer and the tagmented DNA amplified with 5 cycles of PCR using the Enhanced PCR Mix (EPM) and Nextera DNA CD i5 and i7 index adapters (Illumina). Libraries were cleaned and 30 μL of eluted library was transferred to a new sample plate to measure their concentrations via a Quant-iT™ High-Sensitivity dsDNA Assay Kit (Thermo Fisher Scientific) on a Varioskan Lux (Thermo Fisher Scientific) microplate reader. Each sample library was diluted to 4 nM concentration, and 5 μL of the normalized library was denatured with 5 μL of 0.2 N Sodium acetate. A total of 6 Flex libraries were created (2 from gut samples and 4 from hemolymph samples). Samples were pooled and sequenced on an Illumina MiSeq system (Illumina) using a MiSeq Reagent V3 (2 × 300 bp) sequencing kit (Illumina) at the National Center for Agricultural Utilization Research (Peoria, Illinois, United States). Metagenomic data was submitted to the Sequence Read Archive (SRA) database under BioProject accession number PRJNA764167.
Metagenomic paired reads were processed using different bioinformatic applications available on the University of South Florida high performance computing cluster. Raw sequence reads were quality-filtered using Trimmomatic v 0.36.0 (Bolger et al., 2014) with default parameters. FastQC v 0.11.5 (Andrews, 2010) was used to verify the quality of trimmed sequences and assemblies were performed using metaSPAdes v 3.11.1 with default parameters (Nurk et al., 2017). Assembled contigs were filtered by size on the Galaxy web-based platform (Afgan et al., 2018) to retain contigs larger than 100 bp. Contig sequences were compared against the GenBank non-redundant database using BLASTx as implemented in DIAMOND (Buchfink et al., 2015) to identify viral sequences.
BLAST searches revealed large contigs (∼195 kb) with significant similarities to an iridovirus (Liz-CrIV; accession number MN081869) in each of the six libraries. These contigs were further explored using Geneious Prime v 2021.1.1 given that their size approximated near-complete iridovirus genomes. Since iridovirus genomes can be circularly permutated (Jakob et al., 2001), putative unit length genomes were obtained by identifying terminally redundant sequences for annotation purposes. To evaluate genome coverage, quality-trimmed forward reads from each library were mapped to the identified unit length iridovirus genomes using default parameters within the “Map to Reference” Geneious tool. Genome co-linearity was verified using Mauve with the progressiveMauve algorithm (Darling et al., 2010) as implemented in Geneious Prime. A representative viral genome, named cricket iridovirus isolate Liz-CrIV_USDA_2019 (CrIV_USDA), was submitted to GenBank under accession number OK181107. Genome-wide pairwise identities among genomes were calculated by aligning sequences using EMBOSS Stretcher (Madeira et al., 2019). Alignments were then used to calculate pairwise identities using the formula employed by the Species Demarcation Tool, which has been used to classify viral sequences (Muhire et al., 2014). Specifically, the percent identity values were calculated as [1-(M/A)]*100 where M is the number of mismatching nucleotides and A the total number of aligned positions with no gaps.
Statistical Analysis
To compare viral copy abundance across populations and between sexes, we performed general linear models for the ratio of viral copies normalized to cricket tubu3 copies including population (Diseased or Healthy), sex, and virus (CrIV or AdDNV) and their interactions. Individuals with viral copies lower than the limit of quantification were removed from analysis. All abundance data were log transformed to fit normality assumptions and reported results derive from the best models as determined by corrected Akaike’s information criterion (AICc using the stepAIC function in R; Sugiura, 1978; Hurvich and Tsai, 1989) or before the removal of terms from the final model. All statistical analyses were carried out in R (version 4.0.5, R Core Team, 2021) and graphs were made using GraphPad Prism 9 (version 9.0.0). Viral prevalence and abundance data can be accessed in Supplementary Material Data Sheet 1.
Results
Viral Prevalence and Abundance via qPCR
Across 80 crickets sampled, we detected CrIV and AdDNV in all but one individual for each virus (both females from the Healthy population, 98.75% positivity rate for each target) (Table 2). Based on qPCR, we did not find any evidence that either population was infected with additional DNA (AdMADV, GbNV, AdVVV) or RNA viruses (CrPV, AdIV, AdV) (Table 2).
Table 2. Prevalence of targeted viruses across two populations (“Healthy” and “Diseased”) of male and female adult reared Gryllodes sigillatus using qPCR.
Notably, crickets from the Diseased population of both sexes had a much greater abundance of CrIV copies compared to their Healthy counterparts and significantly more copies of CrIV than AdDNV in both populations (Tables 3, 4 and Figure 2). Males and females had similar viral loads of CrIV across both populations. For AdDNV, males and females from the Diseased population had similar viral loads while males had significant lower viral loads compared to females within the Healthy population (Tables 3, 4 and Figure 2). Both male and female crickets from the Diseased population also had a greater abundance of AdDNV compared to the Healthy population (Tables 3, 4 and Figure 2). Estimated amounts of viral copies per cricket are listed in Supplementary Table 1.
Table 3. Model terms and statistics from generalized linear models for log-transformed absolute abundance (number of viral copies) of cricket iridovirus (CrIV) and Acheta domesticus densovirus (AdDNV) detected across two populations (“Healthy” and “Diseased”) of male and female adult reared Gryllodes sigillatus.
Table 4. Post hoc comparisons of interactions from generalized linear models for log-transformed normalized number of viral copies of cricket iridovirus (CrIV) and Acheta domesticus densovirus (AdDNV) detected across two populations (“Healthy” and “Diseased”) of male and female adult reared Gryllodes sigillatus.
Figure 2. Viral loads of cricket iridovirus (CrIV) (left) and Acheta domesticus densovirus (AdDNV) (right) detected across two populations (“Diseased” and “Healthy”) of male and female adult reared Gryllodes sigillatus (geometric mean ± 95% confidence intervals). The graphs depict the ratio of viral copies normalized to cricket tubu3 copies, shown on a log scale. Post hoc comparisons are indicated on Table 4.
Genome Sequence of CrIV via Shotgun Metagenomic Analysis
Metagenomic analysis confirmed the presence of CrIV across all samples. We identified near-complete CrIV genomes with high coverage in each sample (mean coverage ranged from 53x to 1050x), suggesting high viral concentrations. These genomes shared >99.9% identity (genome size between individuals ranged from 194,811 to 195,316 bp). Genome-wide pairwise identities revealed that CrIV isolated from the Diseased population, referred herein as “CrIV_USDA”, was most closely related to Liz-CrIV (accession no. MN081869) (97%) (Figure 3) followed by Invertebrate iridescent virus 6 (IIV6) (accession no. NC_003038) (71.5%) (Supplementary Figure 1). We were not able to detect AdDNV in our metagenomic analysis which could be due to the fact that we had low concentrations in our samples and AdDNV is a ssDNA virus with a relatively small genome (5,425 bp) (Liu et al., 2011; de Miranda et al., 2021a).
Figure 3. Schematic showing a genome-wide pairwise comparison between CrIV_USDA and lizard–cricket iridescent virus (Liz-CrIV). Each block, including “Identity,” “CrIV_USDA,” and “Liz-Cr-IV,” illustrates the shared identity and organization for each genome, respectively. The numbers at the top each block indicate alignment positions. The “Identity” panel highlights identical sites (green), mismatches (red), and gaps (white) across the alignment. Non-overlapping open reading frames encoding putative proteins larger than 80 amino acids in each genome are highlighted with yellow arrows. Areas of the Liz-CrIV genome annotated as gaps (G) in GenBank are highlighted with gray arrows.
Discussion
Entomopathogenic viruses are known to cause significant losses to the reared insect industry (Maciel-Vergara and Ros, 2017); however, little is known about their diversity, biology, and host association. Using a range of molecular approaches in the present study, we were able to identify an cricket iridovirus (CrIV) as the likely causal agent of a disease affecting a colony of reared G. sigillatus, a species of growing importance for food, feed, and industrial purposes (Weissman et al., 2012). Metagenomic analyses revealed that our CrIV isolate, CrIV_USDA, was most closely related to lizard-cricket iridovirus (Liz-CrIV) (Papp et al., 2014; Papp and Marschang, 2019). Gross pathology included milky hemolymph, decreased fecundity, sluggish behavior, and melanotic lesions, all indicative of an active microbial infection. Our molecular identification and quantification assays indicated that Diseased crickets were supporting large numbers of CrIV. Diseased crickets were estimated to have about 1.5 trillion more copies of CrIV than Healthy crickets, based on the volume of DNA template extracted and used in our qPCR assays (Supplementary Table 1). Although our studies did not include a complete Koch’s postulate to confirm the causal agent infecting this colony of G. sigillatus, our molecular characterization based on shotgun metagenomics of hemolymph and gut samples and qPCR-based absolute quantification strongly suggests that CrIV is the main viral entomopathogen driving the pathology in the Diseased population. Future studies will focus on dissecting the intricacies of host-pathogen interactions, routes of transmission, and identifying the biological factors that might be associated with an active and covert infection. The latter consideration will be essential as we do not yet fully understand why CrIV presents as an active infection in one population but not the other.
Invertebrate iridescent viruses (family Iridoviridae) have non-occluded icosahedral particles, approximately 130 nm diameter, that contain a double-stranded DNA genome of about 140–210 kpb (Williams et al., 2000; İnce et al., 2018). They can infect a broad range of invertebrates, including terrestrial isopods, and have been isolated from several insect taxa (Kleespies et al., 1999; Jakob et al., 2002). Furthermore, certain invertebrate iridoviruses are known to cause disease in reptiles (Just et al., 2001; Weinmann et al., 2007; Marschang, 2011; Papp et al., 2014) and amphibians (Stöhr et al., 2016) via ingestion of infected insects, posing a particular health concern for the pet trade industry. In crickets, CrIV is known to cause disease in several species, including reared Gryllus texensis (Adamo et al., 2014), G. bimaculatus (Just and Essbauer, 2001), G. campestris, A. domesticus (Kleespies et al., 1999), and now G. sigillatus.
Our metagenomic analyses confirm the draft genome of Liz-CrIV from a previous study (Papp and Marschang, 2019) as a distinct virus from Invertebrate Iridescent virus 6 (IIV-6), although formal analyses are needed to determine if Liz-CrIV represents a strain of IIV-6 or a new species of invertebrate iridovirus (Papp and Marschang, 2019). Analyses from the current study also provide the most complete genome of Liz-CrIV. Importantly, our PCR primers were not able to discern between IIV-6 and Liz-CrIV and so future screening efforts should target areas of the genome that distinguish the two (see Supplementary Figure 1). Liz-CrIV was first discovered in Europe from commercially produced crickets in the mid to late 1990’s (Kleespies et al., 1999; Just and Essbauer, 2001). The populations used in the present study, descents of field-caught crickets in New Mexico, have been reared in a research lab setting since 2001, and have had no history of contact with commercially produced crickets. While we currently do not know the origins or the spread of this virus, our study confirms that Liz-CrIV is active across multiple continents.
In addition to active CrIV infections, our study found evidence for covert, or asymptomatic, infections of both CrIV and Acheta domesticus densovirus (AdDNV), the latter of which has been documented in G. sigillatus previously (Weissman et al., 2012). Covert infections of invertebrate iridovirus are reportedly more prevalent than active lethal infections in some insect populations (Williams, 1993, 2008; Tonka and Weiser, 2000). Although our study design did not allow us to examine in detail the effects that covert infections might have on cricket health, covert invertebrate iridovirus infections may have significant fitness consequences for hosts (e.g., increased development times and reduced fecundity) (Marina et al., 1999, 2003), which could have important ramifications for product yield in production facilities. Our study design prevents us from distinguishing between a persistent infection (with low levels of virus replication) or a latent infection and additional experiments (e.g., measuring transcription) could elucidate the nature of these covert infections (Williams et al., 2017).
Interestingly, we found a significant effect of the interaction between sex and population on the abundance of viral copies such that Healthy males had a lower abundance of AdDNV copies compared to Diseased males and females from either population. Although females are the larger sex (Sakaluk et al., 2019), our analysis normalized abundance values such that body size would not account for these differences. Thus, our results could be indicative of a higher tolerance of AdDNV by females, but controlled exposures would be needed to confirm this. We did not see the same sex-effect for CrIV. Few studies, if any, quantify viral loads across sex in reared and farmed insects, which makes generalizing this finding difficult. However, previous studies have demonstrated higher immune activity in female G. sigillatus compared with males (Gershman et al., 2010), which could contribute to differences in viral loads.
A key aspect of this system yet to be determined is the route of virus transmission, although per os is suspected to be the main route for CrIV (Williams et al., 2005). Adamo et al. (2014) found that CrIV can be transmitted horizontally via topical exposure in G. texensis and found no evidence of virus within the testes. Similar results were found by Just and Essbauer (2001) who found no signs of infection in the ovaries or testes in G. bimaculatus. However, both studies assessed transmission among populations experiencing active infections and so it is plausible that CrIV may be transmitted vertically when populations are exhibiting covert infections. Indeed, several classes of pathogenic viruses, including iridoviruses, are known to adopt mixed-mode transmission (both vertical and horizontal transmission) based on the relative fitness gains that are obtained via each transmission strategy (Ebert, 2013). Because horizontal transmission is riskier when hosts are rare but vertical transmission is constrained by host fitness, this mixed-mode transmission balances the likelihood of transmission with the constraints of host fitness. Thus, vertical transmission is favored when host densities are low while horizontal transmission is favored when host density is high. Further studies exploring the potential for mixed-mode transmission of CrIV across populations with active and covert infections will be critical and may explain the discrepancy between populations.
In conclusion, foundational research on the pathology of mass-produced insects, including crickets, will be essential to maintain the health and yield necessary for this growing industry. Here, we report the viral loads across sexes of active and covert infections of two pathogenic viruses (CrIV and AdDNV) in two populations of reared G. sigillatus. Importantly, these findings will inform future work addressing diagnostic, mitigation, and therapeutic interventions in reared insect colonies that could ultimately improve product yield and support this burgeoning sustainable industry.
Data Availability Statement
Full genome of the CrIV isolate (Cricket iridovirus isolate Liz-CrIV_USDA_2019) can be found at GenBank (submission number: OK181107). Metagenomic sequences for hemolymph and guts are available at the NCBI Sequence Read Archive database (study number: PRJNA764167, SRA accession numbers: SAMN21469590, SAMN21469591, SAMN21469592, SAMN21469593, SAMN21469594, and SAMN21469595). The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
KD and JR contributed to conception, design of the study, and collected the data. KR and BO analyzed metagenomic data. KD performed statistical analyses and wrote the first draft of the manuscript. KR and JR wrote sections of the manuscript. BS, SS, JH, and RB provided necessary samples, materials, and additional expertise. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This research was supported in part by the U.S. Department of Agriculture, Agricultural Research Service. SS and BS were supported by a grant from the National Science Foundation (IOS 16–54028). JH was supported by a grant from the Australian Research Council (DP180101708). 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. USDA is an equal opportunity provider and employer.
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.
Acknowledgments
We warmly thank Molly Schumacher and Kylie Hampton for their help with data collection.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.780796/full#supplementary-material
References
Adamo, S. A., Kovalko, I., Easy, R. H., and Stoltz, D. (2014). A viral aphrodisiac in the cricket Gryllus texensis. J. Exp. Biol. 217, 1970–1976. doi: 10.1242/jeb.103408
Afgan, E., Baker, D., Batut, B., Van Den Beek, M., Bouvier, D., Čech, M., et al. (2018). The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res. 46, W537–W544. doi: 10.1093/nar/gky379
Andrews, S. (2010). FastQC: a quality control tool for high throughput sequence data. Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc (accessed January, 2021).
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
Bosch, G., Zhang, S., Oonincx, D. G., and Hendriks, W. H. (2014). Protein quality of insects as potential ingredients for dog and cat foods. J. Nutr. Sci. 3:E29. doi: 10.1017/jns.2014.23
Buchfink, B., Xie, C., and Huson, D. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176
Castro-López, C., Santiago-López, L., Vallejo-Cordoba, B., González-Córdova, A. F., Liceaga, A. M., García, H. S., et al. (2020). An insight to fermented edible insects: a global perspective and prospective. Food Res. Int. 137:109750. doi: 10.1016/j.foodres.2020.109750
Darling, A. E., Mau, B., and Perna, N. T. (2010). progressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One 5:e11147. doi: 10.1371/journal.pone.0011147
de Miranda, J. R., Granberg, F., Onorati, P., Jansson, A., and Berggren, Å (2021b). Virus prospecting in crickets—Discovery and strain divergence of a novel iflavirus in wild and cultivated Acheta domesticus. Viruses 13:364. doi: 10.3390/v13030364
de Miranda, J. R., Granberg, F., Low, M., Onorati, P., Semberg, E., Jansson, A., et al. (2021a). Virus Diversity and Loads in Crickets Reared for Feed: implications for Husbandry. Front. Vet. Sci. 8:642085. doi: 10.3389/fvets.2021.642085
Defoliart, G. R. (1995). Edible insects as minilivestock. Biodivers. Conserv. 4, 306–321. doi: 10.1007/BF00055976
Duffield, K. R., Hampton, K. J., Houslay, T. M., Hunt, J., Sadd, B. M., and Sakaluk, S. K. (2019). Inbreeding alters context-dependent reproductive effort and immunity in male crickets. J. Evol. Biol. 32, 731–741. doi: 10.1111/jeb.13478
Ebert, D. (2013). The epidemiology and evolution of symbionts with mixed-mode transmission. Annu. Rev. Ecol. Evol. Syst. 44, 623–643. doi: 10.1146/annurev-ecolsys-032513-100555
Eilenberg, J., Van Oers, M. M., Jensen, A. B., Lecocq, A., Maciel-Vergara, G., Santacoloma, L. P. A., et al. (2018). Towards a coordination of European activities to diagnose and manage insect diseases in production facilities. J. Insects Food Feed 4, 157–166. doi: 10.3920/JIFF2018.0002
Eilenberg, J., Vlak, J. M., Nielsen-Leroux, C., Cappellozza, S., and Jensen, A. B. (2015). Diseases in insects produced for food and feed. J. Insects Food Feed 1, 87–102. doi: 10.3920/JIFF2014.0022
Fernandez-Cassi, X., Supeanu, A., Vaga, M., Jansson, A., Boqvist, S., and Vagsholm, I. (2019). The house cricket (Acheta domesticus) as a novel food: a risk profile. J. Insects Food Feed 5, 137–157. doi: 10.3920/JIFF2018.0021
Gasco, L., Biancarosa, I., and Liland, N. S. (2020). From waste to feed: a review of recent knowledge on insects as producers of protein and fat for animal feeds. Curr. Opin. Green Sustain. Chem. 23, 67–79. doi: 10.1016/j.cogsc.2020.03.003
Gershman, S. N., Barnett, C. A., Pettinger, A. M., Weddle, C. B., Hunt, J., and Sakaluk, S. K. (2010). Inbred decorated crickets exhibit higher measures of macroparasitic immunity than outbred individuals. Heredity 105, 282–289. doi: 10.1038/hdy.2010.1
Hahn, T., Tafi, E., Paul, A., Salvia, R., Falabella, P., and Zibek, S. (2020). Current state of chitin purification and chitosan production from insects. J. Chem. Technol. Biotechnol. 95, 2775–2795. doi: 10.1002/jctb.6533
Hall, F. G., Jones, O. G., O’haire, M. E., and Liceaga, A. M. (2017). Functional properties of tropical banded cricket (Gryllodes sigillatus) protein hydrolysates. Food Chem. 224, 414–422. doi: 10.1016/j.foodchem.2016.11.138
Hawkey, K. J., Lopez-Viso, C., Brameld, J. M., Parr, T., and Salter, A. M. (2021). Insects: a potential source of protein and other nutrients for feed and food. Annu. Rev. Anim. Biosci. 9, 333–354. doi: 10.1146/annurev-animal-021419-083930
Henry, M., Gasco, L., Piccolo, G., and Fountoulaki, E. (2015). Review on the use of insects in the diet of farmed fish: past and future. Anim. Feed Sci. Technol. 203, 1–22. doi: 10.1016/j.anifeedsci.2015.03.001
Huger, A. M. (1985). A new virus disease of crickets (Orthoptera: Gryllidae) causing macronucleosis of fatbody. J. Invertebr. Pathol. 45, 108–111. doi: 10.1016/0022-2011(85)90055-2
Hurvich, C. M., and Tsai, C. L. (1989). Regression and time series model selection in small samples. Biometrika 76, 297–307.
İnce, İA., Özcan, O., Ilter-Akulke, A. Z., Scully, E. D., and Özgen, A. (2018). Invertebrate iridoviruses: a glance over the last decade. Viruses 10:161. doi: 10.3390/v10040161
Jakob, N., Kleespies, R., Tidona, C., Müller, K., Gelderblom, H., and Darai, G. (2002). Comparative analysis of the genome and host range characteristics of two insect iridoviruses: chilo iridescent virus and a cricket iridovirus isolate. J. Gen. Virol. 83, 463–470. doi: 10.1099/0022-1317-83-2-463
Jakob, N. J., Müller, K., Bahr, U., and Darai, G. (2001). Analysis of the first complete DNA sequence of an invertebrate iridovirus: coding strategy of the genome of Chilo iridescent virus. Virology 286, 182–196. doi: 10.1006/viro.2001.0963
Just, F., Essbauer, S., Ahne, W., and Blahak, S. (2001). Occurrence of an invertebrate iridescent-like virus (Iridoviridae) in reptiles. J. Vet. Med. B. 48, 685–694. doi: 10.1046/j.1439-0450.2001.00495.x
Just, F. T., and Essbauer, S. S. (2001). Characterization of an iridescent virus isolated from Gryllus bimaculatus (Orthoptera: Gryllidae). J. Invertebr. Pathol. 77, 51–61. doi: 10.1006/jipa.2000.4985
Kaya, H. K., and Vega, F. E. (2012). “Scope and basic principles of insect pathology,” in Insect pathology, 2nd Edn, eds F. E. Vega and H. K. Kaya (Amsterdamn: Elsevier), 1–12. doi: 10.1016/B978-0-12-384984-7.00001-4
Kleespies, R., Tidona, C., and Darai, G. (1999). Characterization of a new iridovirus isolated from crickets and investigations on the host range. J. Invertebr. Pathol. 73, 84–90. doi: 10.1006/jipa.1998.4821
Liu, K., Li, Y., Jousset, F.-X., Zadori, Z., Szelei, J., Yu, Q., et al. (2011). The Acheta domesticus densovirus, isolated from the European house cricket, has evolved an expression strategy unique among parvoviruses. J. Virol. 85, 10069–10078. doi: 10.1128/JVI.00625-11
Maciel-Vergara, G., Jensen, A. B., Lecocq, A., and Eilenberg, J. (2021). Diseases in edible insect rearing systems. J. Insects Food Feed 7, 621–638. doi: 10.3920/JIFF2021.0024
Maciel-Vergara, G., and Ros, V. I. D. (2017). Viruses of insects reared for food and feed. J. Invertebr. Pathol. 147, 60–75. doi: 10.1016/j.jip.2017.01.013
Madeira, F., Park, Y. M., Lee, J., Buso, N., Gur, T., Madhusoodanan, N., et al. (2019). The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res. 47, W636–W641. doi: 10.1093/nar/gkz268
Magara, H. J., Niassy, S., Ayieko, M. A., Mukundamago, M., Egonyu, J. P., Tanga, C. M., et al. (2021). Edible crickets (Orthoptera) around the world: distribution, nutritional value, and other benefits—a review. Front. Nutr. 7:537915. doi: 10.3389/fnut.2020.537915
Makkar, H. P. S., Tran, G., Henze, V., and Ankers, P. (2014). State-of-the-art on use of insects as animal feed. Anim. Feed Sci. Technol. 197, 1–33. doi: 10.1016/j.anifeedsci.2014.07.008
Marina, C., Ibarra, J., Arredondo-Jiménez, J., Fernández-Salas, I., Valle, J., and Williams, T. (2003). Sublethal iridovirus disease of the mosquito Aedes aegypti is due to viral replication not cytotoxicity. Med. Vet. Entomol. 17, 187–194. doi: 10.1046/j.1365-2915.2003.00422.x
Marina, C. F., Arredondo-Jiménez, J. I., Castillo, A., and Williams, T. (1999). Sublethal effects of iridovirus disease in a mosquito. Oecologia 119, 383–388. doi: 10.1007/s004420050799
Muhire, B. M., Varsani, A., and Martin, D. P. (2014). SDT: a virus classification tool based on pairwise sequence alignment and identity calculation. PLoS One 9:e108277. doi: 10.1371/journal.pone.0108277
Nurk, S., Meleshko, D., Korobeynikov, A., and Pevzner, P. A. (2017). metaSPAdes: a new versatile metagenomic assembler. Genome Res. 27, 824–834.
Osimani, A., Milanović, V., Cardinali, F., Roncolini, A., Garofalo, C., Clementi, F., et al. (2018). Bread enriched with cricket powder (Acheta domesticus): a technological, microbiological and nutritional evaluation. Innov. Food Sci. Emerg. Technol. 48, 150–163. doi: 10.1016/j.ifset.2018.06.007
Papp, T., and Marschang, R. E. (2019). Detection and characterization of invertebrate iridoviruses found in reptiles and prey insects in Europe over the past two decades. Viruses 11:600. doi: 10.3390/v11070600
Papp, T., Spann, D., Marschang, R. E., Mikrobiologie, F. Ä, and Reptilien, Z. (2014). Development and use of a real-time polymerase chain reaction for the detection of group II invertebrate iridoviruses in pet lizards and prey insects. J. Zoo Wildl. Med. 45, 219–227. doi: 10.1638/2012-0044.1
Pham, H. T., Iwao, H., Szelei, J., Li, Y., Liu, K., Bergoin, M., et al. (2013b). Comparative genomic analysis of Acheta domesticus densovirus isolates from different outbreaks in Europe, North America, and Japan. Genome Announc. 1, e00629–13. doi: 10.1128/genomeA.00629-13
Pham, H. T., Yu, Q., Bergoin, M., and Tijssen, P. (2013c). A novel ambisense densovirus, Acheta domesticus mini ambidensovirus, from crickets. Genome Announc. 1, e00914–13. doi: 10.1128/genomeA.00914-13
Pham, H. T., Bergoin, M., and Tijssen, P. (2013a). Acheta domesticus volvovirus, a novel single-stranded circular DNA virus of the house cricket. Genome Announc. 1, e00079–13. doi: 10.1128/genomeA.00079-13
R Core Team (2021). R: A Language and Environment For Statistical Computing. Vienna: R Foundation for Statistical Computing. Available online at: https://www.R-project.org/
Sakaluk, S. K., Oldzej, J., Poppe, C. J., Harper, J. L., Rines, I. G., Hampton, K. J., et al. (2019). Effects of inbreeding on life-history traits and sexual competency in decorated crickets. Anim. Behav. 155, 241–248. doi: 10.1016/j.anbehav.2019.05.027
Stöhr, A. C., Papp, T., and Marschang, R. E. (2016). Repeated detection of an invertebrate iridovirus in amphibians. J. Herpetol. Med. Surg. 26, 54–58. doi: 10.5818/1529-9651-26.1-2.54
Stull, V., and Patz, J. (2020). Research and policy priorities for edible insects. Sustain. Sci. 15, 633–645. doi: 10.1007/s11625-019-00709-5
Stull, V. J., Finer, E., Bergmans, R. S., Febvre, H. P., Longhurst, C., Manter, D. K., et al. (2018). Impact of edible cricket consumption on gut microbiota in healthy adults, a double-blind, randomized crossover trial. Sci. Rep. 8:10762. doi: 10.1038/s41598-018-29032-2
Styer, E., and Hamm, J. (1991). Report of a densovirus in a commercial cricket operation in the southeastern United States. J. Invertebr. Pathol. 58, 283–285. doi: 10.1016/0022-2011(91)90076-3
Sugiura, N. (1978). Further analysis of the data by Akaike’s information criterion and the finite corrections. Comm. Statist. Theory Methods 7, 13–26. doi: 10.1080/03610927808827599
Surendra, K., Olivier, R., Tomberlin, J. K., Jha, R., and Khanal, S. K. (2016). Bioconversion of organic wastes into biodiesel and animal feed via insect farming. Renew. Energy 98, 197–202. doi: 10.1016/j.renene.2016.03.022
Surendra, K. C., Tomberlin, J. K., van Huis, A., Cammack, J. A., Heckmann, L.-H. L., and Khanal, S. K. (2020). Rethinking organic wastes bioconversion: evaluating the potential of the black soldier fly (Hermetia illucens (L.)) (Diptera: Stratiomyidae) (BSF). Waste Manag. 117, 58–80. doi: 10.1016/j.wasman.2020.07.050
Szelei, J., Woodring, J., Goettel, M., Duke, G., Jousset, F.-X., Liu, K., et al. (2011). Susceptibility of North-American and European crickets to Acheta domesticus densovirus (AdDNV) and associated epizootics. J. Invertebr. Pathol. 106, 394–399. doi: 10.1016/j.jip.2010.12.009
Tomberlin, J. K., van Huis, A., Benbow, M. E., Jordan, H., Astuti, D. A., Azzollini, D., et al. (2015). Protecting the environment through insect farming as a means to produce protein for use as livestock, poultry, and aquaculture feed. J. Insects Food Feed 1, 307–309. doi: 10.3920/JIFF2015.0098
Tonka, T., and Weiser, J. (2000). Iridovirus infection in mayfly larvae. J. Invertebr. Pathol. 76, 229–231. doi: 10.1006/jipa.2000.4971
Udomsil, N., Imsoonthornruksa, S., Gosalawit, C., and Ketudat-Cairns, M. (2019). Nutritional values and functional properties of house cricket (Acheta domesticus) and field cricket (Gryllus bimaculatus). Food Sci. Technol. 25, 597–605. doi: 10.3136/fstr.25.597
Valles, S. M., and Chen, Y. (2006). Serendipitous discovery of an RNA virus from the cricket, Acheta domesticus. Fla. Entomol. 89, 282–283. doi: 10.1653/0015-4040(2006)89[282:SDOARV]2.0.CO;2
van Huis, A. (2017). Edible insects and research needs. J. Insects Food Feed 3, 3–5. doi: 10.3920/JIFF2017.x002
van Huis, A. (2020b). Insects as food and feed, a new emerging agricultural sector: a review. J. Insects Food Feed 6, 27–44. doi: 10.3920/JIFF2019.0017
van Huis, A. (2020a). Edible crickets, but which species? J. Insects Food Feed 6, 91–94. doi: 10.3920/JIFF2020.x001
van Huis, A., Van Itterbeeck, J., Klunder, H., Mertens, E., Halloran, A., Muir, G., et al. (2013). Edible insects: future prospects for food and feed security. Rome: Food and Agriculture Organization of the United Nations.
Wang, D., Bai, Y. Y., Li, J. H., and Zhang, C. X. (2004). Nutritional value of the field cricket (Gryllus testaceus Walker). Insect Sci. 11, 275–283. doi: 10.1111/j.1744-7917.2004.tb00424.x
Wang, L., Cappelle, K., Santos, D., Broeck, J. V., Smagghe, G., and Swevers, L. (2019). Short-term persistence precedes pathogenic infection: infection kinetics of cricket paralysis virus in silkworm-derived Bm5 cells. J. Insect Physiol. 115, 1–11. doi: 10.1016/j.jinsphys.2019.03.004
Wang, Y., Kleespies, R. G., Huger, A. M., and Jehle, J. A. (2007). The genome of Gryllus bimaculatus nudivirus indicates an ancient diversification of baculovirus-related nonoccluded nudiviruses of insects. J. Virol. 81, 5395–5406. doi: 10.1128/JVI.02781-06
Weinmann, N., Papp, T., De Matos, A. P. A., Teifke, J. P., and Marschang, R. E. (2007). Experimental infection of crickets (Gryllus bimaculatus) with an invertebrate iridovirus isolated from a high-casqued chameleon (Chamaeleo hoehnelii). J. Vet. Diagn. Invest. 19, 674–679. doi: 10.1177/104063870701900609
Weissman, D. B., Gray, D. A., Pham, H. T., and Tijssen, P. (2012). Billions and billions sold: pet-feeder crickets (Orthoptera: Gryllidae), commercial cricket farms, an epizootic densovirus, and government regulations make for a potential disaster. Zootaxa 3504, 67–88. doi: 10.11646/zootaxa.3504.1.3
Wilkie, R. M. (2018). ‘Minilivestock’ farming: who is farming edible insects in Europe and North America? J. Sociol. 54, 520–537. doi: 10.1177/1440783318815304
Williams, T. (1993). Covert iridovirus infection of blackfly larvae. Proc. R. Soc. B. 251, 225–230. doi: 10.1098/rspb.1993.0033
Williams, T. (2008). Natural invertebrate hosts of iridoviruses (Iridoviridae). Neotrop. Entomol. 37, 615–632. doi: 10.1590/S1519-566X2008000600001
Williams, T., Barbosa-Solomieu, V., and Chinchar, V. G. (2005). A decade of advances in iridovirus research. Adv. Virus Res. 65, 173–248. doi: 10.1016/S0065-3527(05)65006-3
Williams, T., Chinchar, V. G., Darai, G., Hyatt, A., Kalmakoff, J., and Seligy, V. (2000). “Iridoviridae,” in Virus taxonomy, VII Report of the International Committee on Virus Taxonomy, eds M. H. V. van Regenmortel, C. M. Fauquet, D. H. L. Bishop, E. B. Carstens, and M. K. Estes (New York: Academic Press), 167–182.
Williams, T., Virto, C., Murillo, R., and Caballero, P. (2017). Covert infection of insects by baculoviruses. Front. Microbiol. 8:1337. doi: 10.3389/fmicb.2017.01337
Keywords: reared crickets, insects as food and feed, Gryllodes sigillatus, cricket viruses, entomopathogenic viruses
Citation: Duffield KR, Hunt J, Sadd BM, Sakaluk SK, Oppert B, Rosario K, Behle RW and Ramirez JL (2021) Active and Covert Infections of Cricket Iridovirus and Acheta domesticus Densovirus in Reared Gryllodes sigillatus Crickets. Front. Microbiol. 12:780796. doi: 10.3389/fmicb.2021.780796
Received: 21 September 2021; Accepted: 01 November 2021;
Published: 30 November 2021.
Edited by:
Simon Luke Elliot, Universidade Federal de Viçosa, BrazilReviewed by:
Jiaqin Xie, Chongqing University, ChinaJoachim R. de Miranda, Swedish University of Agricultural Sciences, Sweden
Copyright © 2021 Duffield, Hunt, Sadd, Sakaluk, Oppert, Rosario, Behle and Ramirez. 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: Kristin R. Duffield, Kristin.Duffield@USDA.gov