<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Insect Sci.</journal-id>
<journal-title>Frontiers in Insect Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Insect Sci.</abbrev-journal-title>
<issn pub-type="epub">2673-8600</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/finsc.2023.1059046</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Insect Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Sequencing to identify pathogens in <italic>Tenebrio molitor</italic>: Implications in insects farmed for food and feed</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Leierer</surname>
<given-names>Dewey</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2076590"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Olmstead</surname>
<given-names>Morgan</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Oppert</surname>
<given-names>Brenda</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/205598"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Department of Biochemistry and Molecular Biophysics, Kansas State University</institution>, <addr-line>Manhattan, KS</addr-line>, <country>United States</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>United States Department of Agriculture Agricultural (USDA) Research Service, Center for Grain and Animal Health Research</institution>, <addr-line>Manhattan, KS</addr-line>, <country>United States</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Sudeshna Mazumdar-Leighton, University of Delhi, India</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Hongyu Zhang, Huazhong Agricultural University, China; Thomas Lefebvre, Ynsect, France</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Brenda Oppert, <email xlink:href="mailto:Brenda.oppert@usda.gov">Brenda.oppert@usda.gov</email>
</p>
</fn>
<fn fn-type="other" id="fn002">
<p>This article was submitted to Insect Molecular Genetics, a section of the journal Frontiers in Insect Science</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>30</day>
<month>01</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>3</volume>
<elocation-id>1059046</elocation-id>
<history>
<date date-type="received">
<day>30</day>
<month>09</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>04</day>
<month>01</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Leierer, Olmstead and Oppert</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Leierer, Olmstead and Oppert</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>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.</p>
</license>
</permissions>
<abstract>
<sec>
<title>Introduction</title>
<p>The farmed insect industry is increasing in number and size to meet the demand for sustainably-produced protein. Larger insect farms are prone to losses due to pathogens, and more information is needed regarding the health of insects reared for food and feed.</p>
</sec>
<sec>
<title>Methods</title>
<p>In this study, high throughput sequencing was used to identify potential pathogens in a colony of <italic>Tenebrio molitor</italic> (yellow mealworm, Coleoptera: Tenebrionidae) that exhibited increased mortality in immature stages with eventual colony collapse. Sequences also were obtained from a healthy new colony of <italic>T. molitor</italic>, as well as a recovered individual from the collapsed colony.</p>
</sec>
<sec>
<title>Results</title>
<p>Screening of sequences obtained from the colonies and their rearing diet indicated that the collapsed colony had low diversity in microbial taxa, with predominantly sequences from the families Staphylococcaeceae and Streptococcaceae constituting from 53 to 88% of the total microbial reads. Conversely, in the new colony and their rearing diet, microbial sequences were from more than 15 different taxa, with Lactobacilleceae the most prevalent but representing only 21% of the total microbial reads. Evidence indicates that <italic>Bacillus thuringiensis</italic> may have been involved in the collapse of the colony, leading to sepsis and microbial dysbiosis, although the source of the bacteria was not identified. Sequences from the recovered individual reflected a microbial flora profile that was intermediate between those of the diseased collapsed and new colonies.</p>
</sec>
<sec>
<title>Discussion</title>
<p>These findings have implications for insects reared in confined environments and provide a rapid method to screen insect colonies by sequencing healthy and potentially diseased individuals.</p>
</sec>
</abstract>
<kwd-group>
<kwd>
<italic>Tenebrio molitor</italic>
</kwd>
<kwd>yellow mealworm</kwd>
<kwd>insects for food</kwd>
<kwd>insect pathogens</kwd>
<kwd>shotgun metagenomics</kwd>
<kwd>farmed insects</kwd>
</kwd-group>
<counts>
<fig-count count="3"/>
<table-count count="3"/>
<equation-count count="0"/>
<ref-count count="28"/>
<page-count count="7"/>
<word-count count="3491"/>
</counts>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>The global population is growing exponentially, with current forecasts exceeding 9 billion people by 2050 (<xref ref-type="bibr" rid="B1">1</xref>). To feed a population of that size, annual food production must increase by 70% and annual meat production must increase by 200 million tons. Obstacles to increased food production include climate change, depleting water sources, and restrictions on expanding farmable land. In addition, animal feed production costs from protein-rich ingredients like fishmeal and soybean meal are increasing (<xref ref-type="bibr" rid="B2">2</xref>). Therefore, sustainable alternative protein sources are needed to supplement food and feed, and insects are promising candidates. Insects are high in nutrients like insoluble fiber, proteins, vitamins, and minerals (<xref ref-type="bibr" rid="B3">3</xref>). Insects also consume less water per gram of protein produced and require lower input feeds than most traditional livestock (<xref ref-type="bibr" rid="B4">4</xref>). The yellow mealworm, <italic>Tenebrio molitor</italic>, is one of the most popular and well-researched insects for supplementing animal feed. Companies already mass-rear <italic>T. molitor</italic> on an industrial scale for feeding domestic birds, fish, and reptiles (<xref ref-type="bibr" rid="B5">5</xref>).</p>
<p>There is considerable research on diseases affecting <italic>T. molitor</italic> because of its notoriety as a pest of poultry houses and its use as a lab model for coleopteran biochemical research. <italic>Pseudomonas</italic> and <italic>Serratia</italic> species are opportunistic bacterial pathogens that infect <italic>T. molitor</italic> colonies under stress (<xref ref-type="bibr" rid="B6">6</xref>, <xref ref-type="bibr" rid="B7">7</xref>). Gregarines are known apicomplexan parasites of <italic>T. molitor</italic>, but infection typically only results in a decrease in the adult lifespan, not affecting the fitness of larval stages (<xref ref-type="bibr" rid="B8">8</xref>, <xref ref-type="bibr" rid="B9">9</xref>). Entomopathogenic fungi, including those from the genera <italic>Metarhizium</italic> and <italic>Beauveria</italic> (<xref ref-type="bibr" rid="B10">10</xref>), as well as nematodes (<xref ref-type="bibr" rid="B11">11</xref>) also have been documented as potential biological control agents for <italic>T. molitor</italic>. Our research has used <italic>T. molitor</italic> as a model for studies of the mode of action of <italic>Bacillus thuringiensis</italic> Cry3Aa toxins in coleopterans (<xref ref-type="bibr" rid="B12">12</xref>&#x2013;<xref ref-type="bibr" rid="B16">16</xref>).</p>
<p>Often insects in rearing facilities are in population-dense environments with a high potential for pathogen spread throughout the colony (<xref ref-type="bibr" rid="B17">17</xref>). Insects do not have an adaptive immune system, so disease prevention methods that rely on learned immunity, like vaccines, cannot be used to prevent illness in insects as they do with traditional livestock. Although some entomopathogens may be prevented or treated with antibiotics, overuse of antibiotics in other farming industries has contributed to the development of antibiotic-resistant bacteria (<xref ref-type="bibr" rid="B18">18</xref>). Farmed insects represent considerable investment, and as such there is increased interest in identifying diseases that may affect insects produced for feed, and development of prophylactics and treatments for insect diseases in reared systems (<xref ref-type="bibr" rid="B19">19</xref>, <xref ref-type="bibr" rid="B20">20</xref>). Current prevention and treatment recommendations are suitable for smaller rearing facilities or for those who rear insects at home but are not practical for larger facilities (<xref ref-type="bibr" rid="B5">5</xref>). Our long-term goal is to assist larger facilities in avoiding the practice of removing infected colonies and replacing with new stock, especially with the costs of genetically-selected or modified strains. More research is needed on the detection, treatment, and prevention of pathogens that may affect mass-rearing efficiency and colony health of mealworms.</p>
<p>This case study uses DNA sequencing to identify pathogens that may be responsible for the collapse of a <italic>T. molitor</italic> colony that had decreased survival of egg to adult and subsequent mite infestation. The data identifies presumably beneficial bacteria in a healthy mealworm microbiome and pathogenic bacteria that were associated with the <italic>T. molitor</italic> colony collapse.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>Methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Colony rearing</title>
<p>
<italic>T. molitor</italic> colonies were reared on a diet consisting of 90% organic Golden Buffalo&#x2122; whole wheat flour, 5% brewer&#x2019;s yeast, and 5% stabilized wheat germ for more than 20 years. Slices of sweet potatoes were added to colonies to provide moisture and additional nutrients. Colonies were kept at 27 &#xb0;C, 65 R.H in 2 L glass jars. When increased larval mortality was observed over time, the colony was moved to a separate chamber to be reared in quarantine (25 &#xb0;C, 65 R.H), and a new colony was established in October, 2015 from commercial sources (reptilefood.com, vita-mealie.weebly.com, and rainbowmealworms.net) in the original chamber (27 &#xb0;C, 65 R.H). The new colony was reared in ventilated plastic boxes, ranging in size: 19 x 12.5 x 8 cm (small larvae, 0-1 mo), 35 x 20 x 12 cm (medium larvae, 1-3 mo, and 39 x 34 x 8 cm (&gt;3 mo, large larvae, pupae, adults). After the collapse of the original colony, the few surviving larvae of the quarantined colony (approximately 15) were placed in a 100 mL plastic cup with fresh diet. These larvae were kept in the lab at room temperature at ambient lab humidity.</p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>DNA extraction and sequencing</title>
<p>DNA was extracted from samples including: samples from the collapsed colony - <italic>T. molitor</italic> adult (C-Adult), frass (C-Frass), and Diet (C-Diet); a recovered larva from a colony established from the remaining members of the collapsed colony (R-Larva); a pupa (N-Pupa) and adult (N-Adult) from a new colony and a diet sample (N-Diet) from the rearing container of the new colony containing both flour and frass (<xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>). At the time of sequencing, the original colony already was collapsing with few survivors. One dead adult was sequenced, the diet it had been reared on, and the frass from the colony. The remaining larvae (about 15) were placed on fresh diet to reduce mite infestation and to see if recovery was possible &#x2013; one of the larvae (R-Larva) was sequenced. An extraction of DNA from a sample of pre-rearing diet (i.e., &#x201c;fresh diet&#x201d;, F-Diet) also was sequenced. Samples were stored at -20&#xb0;C until extraction.</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Summary of sample type, description, and collection date.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="left">Abbreviation</th>
<th valign="middle" align="center">Sample</th>
<th valign="middle" align="center">Description</th>
<th valign="middle" align="center">Collection Date</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">C-Adult</td>
<td valign="middle" align="center">Collapsed Colony Adult</td>
<td valign="middle" align="center">Adult collected from colony at time of collapse (dead).</td>
<td valign="middle" align="left">May 2021</td>
</tr>
<tr>
<td valign="middle" align="left">C-Frass</td>
<td valign="middle" align="center">Collapsed Colony Frass</td>
<td valign="middle" align="center">Frass collected from colony at time of collapse.</td>
<td valign="middle" align="left">May 2021</td>
</tr>
<tr>
<td valign="middle" align="left">C-Diet</td>
<td valign="middle" align="center">Collapsed Colony Diet</td>
<td valign="middle" align="center">Flour substrate collected from collapsed colony at time of collapse.</td>
<td valign="middle" align="left">May 2021</td>
</tr>
<tr>
<td valign="middle" align="left">R-Larva</td>
<td valign="middle" align="center">Recovered Larva</td>
<td valign="middle" align="center">Larva collected from the recovered colony (living).</td>
<td valign="middle" align="left">June 2021</td>
</tr>
<tr>
<td valign="middle" align="left">N-Adult</td>
<td valign="middle" align="center">New Colony Adult</td>
<td valign="middle" align="center">Adult collected from the newly established colony (living).</td>
<td valign="middle" align="left">January 2022</td>
</tr>
<tr>
<td valign="middle" align="left">N-Diet</td>
<td valign="middle" align="center">New Colony Diet</td>
<td valign="middle" align="center">Rearing diet from the new colony</td>
<td valign="middle" align="left">January 2022</td>
</tr>
<tr>
<td valign="middle" align="left">N-Pupa</td>
<td valign="middle" align="center">New Colony Pupa</td>
<td valign="middle" align="center">Pupa collected from newly established colony (living).</td>
<td valign="middle" align="left">January 2022</td>
</tr>
<tr>
<td valign="middle" align="left">F-Diet</td>
<td valign="middle" align="center">Fresh Diet</td>
<td valign="middle" align="center">Sample of pre-rearing diet.</td>
<td valign="middle" align="left">August 2022</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Samples were ground in kit buffer (E.Z.N.A. Insect DNA kit, Omega Biotek, Norcross, GA, USA) with a micropestle in a 1.5 mL microcentrifuge tube and DNA was extracted using 4 &#x3bc;l RNase (instead of 2), an optional column equilibrium step, and an additional wash step prior to matrix drying, all variations from the recommended protocol. The extraction kit was optimized for extracting from samples with high polysaccharide content. Sequencing libraries were generated using an NGS Library Preparation Enzymatic Fragmentation kit (Twist Biosciences, CA, USA, version 2.0) following manufacturer&#x2019;s protocol. Average library size was 460 bp for the Collapsed, Recovered, and New colony samples and 363 bp for the F-Diet sample. Sequencing was on a MiSeq (Illumina, San Diego, CA) with v3 600 cycle and v2 300 cycle nano kits, 151X151 bp paired end. Sequencing metrics are provided in <xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>.</p>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Sequencing metrics and classification of reads by Kraken2 (<xref ref-type="bibr" rid="B21">21</xref>).</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="left">Sample</th>
<th valign="middle" align="center">Total reads</th>
<th valign="middle" align="center">Total classified reads</th>
<th valign="middle" align="center">Percent Classified reads</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">C-Adult</td>
<td valign="middle" align="center">7,322,142</td>
<td valign="middle" align="center">671,310</td>
<td valign="middle" align="center">9.17%</td>
</tr>
<tr>
<td valign="middle" align="left">C-Frass</td>
<td valign="middle" align="center">7,584,608</td>
<td valign="middle" align="center">2,163,158</td>
<td valign="middle" align="center">28.52%</td>
</tr>
<tr>
<td valign="middle" align="left">C-Diet</td>
<td valign="middle" align="center">8,990,803</td>
<td valign="middle" align="center">388,244</td>
<td valign="middle" align="center">4.32%</td>
</tr>
<tr>
<td valign="middle" align="left">R-Larva</td>
<td valign="middle" align="center">4,514,734</td>
<td valign="middle" align="center">93,404</td>
<td valign="middle" align="center">2.07%</td>
</tr>
<tr>
<td valign="middle" align="left">N-Adult</td>
<td valign="middle" align="center">3,774,632</td>
<td valign="middle" align="center">202,973</td>
<td valign="middle" align="center">5.38%</td>
</tr>
<tr>
<td valign="middle" align="left">N-Diet</td>
<td valign="middle" align="center">9,253,710</td>
<td valign="middle" align="center">348,601</td>
<td valign="middle" align="center">5.75%</td>
</tr>
<tr>
<td valign="middle" align="left">N-Pupa</td>
<td valign="middle" align="center">4,459,590</td>
<td valign="middle" align="center">531,933</td>
<td valign="middle" align="center">7.82%</td>
</tr>
<tr>
<td valign="middle" align="left">F-Diet</td>
<td valign="middle" align="center">23,137,004</td>
<td valign="middle" align="center">663,383</td>
<td valign="middle" align="center">2.87%</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Data analysis</title>
<p>Fastq files were converted to fasta files and were submitted to Kraken2 analysis (<xref ref-type="bibr" rid="B21">21</xref>) in Omicsbox (Biobam, Valencia, Spain, version 2.1.4 confidence interval filter of 0.4, RefSeq 2021.04). Classified reads from the superkingdom Bacteria were extracted into a separate file and resubmitted to Kraken2 for further analysis.</p>
<p>A bacterial species count table was exported from Omicsbox to compare the number of reads associated with each species across the collapsed, recovered, and new colony samples (<xref ref-type="supplementary-material" rid="SF1">
<bold>Supplementary File 1</bold>
</xref>). Reads from each sample were combined into categories (C-Adult, C-Frass, C-Diet in the collapsed colony and N-Adult, N-Diet, and N-Pupa in the new colony). Only species with counts over 100 reads in at least one colony group were used in the calculations. Classified microbial reads that were identified to the species level as <italic>B. thuringiensis</italic> were selected and analyzed by blastn at NCBI (accessed July 28, 2022, default parameters).</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<p>This study investigated a <italic>T. molitor</italic> colony that demonstrated reduced viability and eventual collapse. In May 2021, the colony suffered a flour mite infestation that likely contributed to the collapse by competing for resources (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1A</bold>
</xref>), with only a few <italic>T. molitor</italic> larvae surviving. Larvae from the collapsed colony had reduced movement with darkening and melanization of the gut, evidence of bacterial infections in <italic>T. molitor</italic> (20, <xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;1B</bold>
</xref>). When problems were identified with the original colony, <italic>T. molitor</italic> were obtained from three commercial sources and eventually combined into a colony labeled &#x201c;new&#x201d;. The approximately 15 surviving larvae from the collapsed colony were transferred to a new container with fresh diet and remained quarantined from the newly established colony. With careful monitoring, the survivors of the collapsed colony, labeled &#x201c;recovered&#x201d;, were able to pupate, reproduce, and began to reestablish. As time passed, there were fewer larvae exhibiting signs of infection, and the mite population was brought under control with careful hygiene.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Observations of the collapsed colony of <italic>T. molitor.</italic> <bold>(A)</bold> Mites (unidentified) isolated from collapsed colony. <bold>(B)</bold> Larva from the newly established colony (left) and larvae recovered from the collapsed colony. Melanized gut is encircled.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="finsc-03-1059046-g001.tif"/>
</fig>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Family taxa found in different samples found in <italic>T. molitor</italic> DNA sequences. Abbreviations for samples are in <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="finsc-03-1059046-g002.tif"/>
</fig>
<p>We sequenced DNA from an adult from the collapsed colony (C-Adult), as well as their frass (C-Frass) and diet (C-Diet), and DNA from a recovered larva (R-Larva) (<xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>). We also sequenced DNA from a new colony adult (N-Adult) and their frass and diet combined (N-Diet). Sequencing provided from 3.8 to 9.3 million reads for downstream analysis (<xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>). Additionally, we sequenced over 23 million reads from fresh diet (F-Diet). Reads were analyzed for similarity to sequences in microbial databases. The majority of reads in all samples were unclassified, as they were likely host (<italic>T. molitor</italic>, <xref ref-type="supplementary-material" rid="SF1">
<bold>Additional File S1A</bold>
</xref>). At the superkingdom level, most reads were bacterial (<xref ref-type="supplementary-material" rid="SF1">
<bold>Additional File S1B</bold>
</xref>), the focus of this study. Diets had a higher percentage (16.4-29.3%) of eukaryotic reads than insect samples (0.99-9.93%) due mostly to contribution from reads similar to <italic>Saccharomycetes</italic> and likely due to the yeast in the diet (<xref ref-type="supplementary-material" rid="SF1">
<bold>Additional File S1A</bold>
</xref>). Some reads were assigned to <italic>Mammalia</italic>, but likely those were insect reads aligning to orthologs in the database. Fewer reads were associated with archea and viruses.</p>
<p>The lowest percentage of classified reads was from R-Larva DNA (2%), whereas 28% of reads from the C-Frass sample were classified (<xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>). We compared reads from the major bacterial classifications among the samples, focusing on of those were assigned to families at 1% or more of the total reads (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>). Overall, the taxonomic distribution of the bacterial reads from the new colony samples were more similar, and those of the collapsed colony were similar. The taxonomic distribution of the bacterial reads from the R-Larva sample was intermediate to the distributions of those from the new colony and the collapsed colony samples.</p>
<p>Reads derived from the family Staphylococcaeceae were predominant in the microbial reads from the collapsed colony samples, constituting 53, 76, and 88% of the total classified reads in C-Adult, C-Frass, and C-Diet, respectively (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>). In contrast, reads associated with Staphylococcaeceae in the new colony samples were 12, 13, and 30% for N-Diet, N-Adult, and N-Pupa, respectively. Reads from the bacterial family Streptococcaceae also were prevalent in the collapsed colony, constituting 4-23% of classified reads compared to 0-6% in the new colony. The newly established colony had higher percentages of bacteria from the families Lactobacillaceae (21-47%), Paenibacillacea<italic>e</italic> (4-9%), and Petrotogacea<italic>e</italic> (4-19%) than the collapsed colony (2-10, 0, 0-2%, respectively).</p>
<p>At the genus level, large increases in reads from the collapsed colony samples were observed from the genera <italic>Mammaliicoccus</italic>, <italic>Lactococcus</italic>, <italic>Staphylococcus</italic> compared to those from the new colony samples (<xref ref-type="supplementary-material" rid="SF2">
<bold>Additional File S2</bold>
</xref>). However, the collapsed colony samples had fewer reads from the genera <italic>Lactobacillus</italic> and <italic>Lacticaseibacillus</italic> than those from the new colony samples, in agreement with observations at the family level.</p>
<p>We were particularly interested in the family Bacillaceae since the phenotype of some of the collapsed insects suggested an active <italic>B. thuringiensis</italic> infection (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>). Therefore, we examined sequences in the genus <italic>Bacillus</italic> through a targeted blast search to the <italic>B. thuringiensis</italic> database at NCBI. There were a larger number of reads from the Bacillaceae family in the samples from the new colony, but most of these were from different <italic>Bacillus</italic> species and other Bacillaceae genera (<xref ref-type="supplementary-material" rid="SF3">
<bold>Additional File S3</bold>
</xref>). More reads with high scoring blastn matches to <italic>B. thuringiensis</italic> were found in the collapsed colony samples, especially in the C-Frass sample (<xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>). The number of reads in the new colony samples were too low (<xref ref-type="bibr" rid="B1">1</xref>, <xref ref-type="bibr" rid="B2">2</xref>) to confirm the presence of <italic>B. thuringiensis</italic>. However, 240 reads from C-Frass aligned to <italic>Bacillus</italic> spp, and 74 aligned to <italic>B. thuringiensis</italic> serovar morrisoni str. 4AA1 (taxid ID 1461024) with evalues &lt; 10<sup>-40</sup> (<xref ref-type="supplementary-material" rid="SF3">
<bold>Additional File S3</bold>
</xref>).</p>
<table-wrap id="T3" position="float">
<label>Table&#xa0;3</label>
<caption>
<p>Microbial reads from <italic>T. molitor</italic> samples that were similar to <italic>B. thuringiensis</italic> serovars, with more reads aligning to serovar <italic>Morrisonii</italic> in bold.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="left">Sample</th>
<th valign="middle" align="center">Number of reads aligned to <italic>Bacillus thuringiensis</italic>
</th>
<th valign="middle" align="center">Serovar (number of reads)</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">N-Pupa</td>
<td valign="top" align="center">1</td>
<td valign="top" align="left">None</td>
</tr>
<tr>
<td valign="top" align="left">N-Adult</td>
<td valign="top" align="center">3</td>
<td valign="top" align="left">
<italic>Alesti</italic> (2)</td>
</tr>
<tr>
<td valign="top" align="left">N-Diet</td>
<td valign="top" align="center">8</td>
<td valign="top" align="left">
<italic>Monterrey</italic> (1)</td>
</tr>
<tr>
<td valign="top" align="left">R-Larva</td>
<td valign="top" align="center">1</td>
<td valign="top" align="left">None</td>
</tr>
<tr>
<td valign="top" align="left">C-Adult</td>
<td valign="top" align="center">8</td>
<td valign="top" align="left">
<italic>Morrisoni</italic> (4)</td>
</tr>
<tr>
<td valign="top" align="left">C-Diet</td>
<td valign="top" align="center">6</td>
<td valign="top" align="left">
<italic>Andoulensis</italic> (1)</td>
</tr>
<tr>
<td valign="top" align="left">C-Frass</td>
<td valign="top" align="center">149</td>
<td valign="top" align="left">
<italic>Alesti</italic> (6)<break/>
<italic>Andalousiensis</italic> (8)<break/>
<italic>Huazhongensis</italic> (1)<break/>
<italic>Kurstaki</italic> (2)<break/>
<bold>
<italic>Morrisonii</italic> (74)</bold>
<break/>
<italic>Pakistani</italic> (24)<break/>
<italic>Pulsiensis</italic> (2)</td>
</tr>
<tr>
<td valign="top" align="left">F-Diet</td>
<td valign="top" align="center">0</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Fresh samples of rearing diet (F-Diet) were sequenced to confirm that was not the original source of <italic>B. thuringiensis</italic> in our collapsed <italic>T. molitor</italic> colony. The number of F-Diet sample reads were more than 4-fold greater than those in the other samples to ensure that <italic>B. thuringiensis</italic> would be identified if present (<xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>). However, none of the reads in F-Diet were identified as <italic>B. thuringiensis</italic>, indicating that the rearing diet likely was not the source of <italic>B. thuringiensis</italic> found in the other samples (<xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>).</p>
<p>A comparison of microbial sequences from diets demonstrated the progression of loss of microbial diversity in pre-rearing diet to diet collected from the collapsed colony (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>). The pre-rearing diet (F-Diet) only had 0.4% of microbial reads but from a more diverse bacteria flora, whereas the bacterial profile in diet from the new colony (N-Diet), also 0.4% of total reads, had shifted to primarily Lactobacilleceae. In contrast, microbial sequences in the diet from the collapsed colony (C-Diet), with 2.5% of the classified sequences as bacteria, had shifted to Staphylococceae with fewer other bacterial sequences found.</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Comparison of bacterial taxa in sequences from pre-rearing diet (F-Diet) compared to diets from the new colony (N-Diet) and collapsed colony (C-Diet).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="finsc-03-1059046-g003.tif"/>
</fig>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<p>High throughput sequencing can be an effective tool to understand the microbial profiles in insects, and this approach can be even more effective in diagnostics of diseased colonies or routine monitoring of farmed insects. Successful results are dependent on sequencing depth and the accuracy of microbial databases. We were limited by sample size in the collapsed colony, as only a few individuals survived. While we were able to identify meaningful and actionable data from the small number of samples in this study, improvements could be made with additional samples from different developmental stages and increased sequencing depth. As more insect and microbial sequences are acquired by data depositories, the accuracy of predicting microbial profiles will be improved.</p>
<p>In this study, we found that reads from the families Staphylococcaeceae and Streptococcaceae collectively constituted 76% of microbial families in the adult sequenced from the collapsed colony, with the bacterial communities of the frass and diet from that colony also containing large numbers of reads (93% and 86%, respectively, of the total bacterial reads in each sample). This narrow microbial diversity represented a shift from the more diverse community found in the adult from the new colony, with only 21% of bacteria from Lactobacilleceae, and the remainder of microbial reads spread over more than 15 different families. Similar profiles were found in the diet from the new colony. There also was a difference in the overall relative abundance of bacterial reads found in the two colonies, as the total number of microbial reads in the new colony was 62-68% for the sequenced larva and adult, whereas more microbial reads (85%) were found in the collapsed adult.</p>
<p>The phenotype of the diseased adult in the collapsed colony suggested there was an active <italic>B. thuringiensis</italic> infection, and our sequencing data support this observation. There were 74 reads from the collapsed colony diet that aligned to <italic>B. thuringiensis</italic> serovar <italic>morrisoni</italic>, and this serovar typically expresses Cry3Aa toxins that are insecticidal to Coleoptera (<xref ref-type="bibr" rid="B22">22</xref>, <xref ref-type="bibr" rid="B23">23</xref>). However, it is unclear how the infection occurred, although finding a small number of reads aligning to <italic>B. thuringiensis</italic> in the new colony adult and diet suggests that these bacteria were present, possibly <italic>via</italic> the sweet potatoes that are introduced to larvae and adults during rearing.</p>
<p>These data do not inform us as to whether the bacteria identified in the collapsed adult contributed to the colony collapse, but more likely they were a result of the sepsis that occurs in the intestinal tract during infection with the insecticidal toxins from <italic>B. thuringiensis. Staphylococcus aureus</italic> was found in the intestinal tract of <italic>Musca domesticus</italic> and is speculated to be contained in that insect <italic>via</italic> excretion, digestion, and antimicrobial peptides, indicating that insects may have an active system to contain harmful bacteria (<xref ref-type="bibr" rid="B24">24</xref>). Streptococci, specifically enterococci, have been identified in more than half of insects surveyed, and it was speculated that insects may serve as an overwintering host for these bacteria (<xref ref-type="bibr" rid="B25">25</xref>). Given that the adult sequenced from the adult colony was dead, we cannot know if the microbiota populations were a result of postmortem changes in the insect. On the other hand, our data suggests that Lactobacilli may be beneficial to mealworms as they were found in 21% of the microbial reads from the healthy adult in the new colony. In Hymenoptera, evidence supports that hypothesis, in that some bees have host-specific Lactobacilli, and may be important for bee health (<xref ref-type="bibr" rid="B26">26</xref>). Lactobacilli can check the growth of other bacteria and are often components of probiotics, and these bacteria were among those suggested as beneficial for farmed insects (<xref ref-type="bibr" rid="B27">27</xref>). Our data supports the hypotheses that all of these bacteria are part of the normal flora of the mealworm gut tract and only become pathogenic when other factors disrupt their balance in the microbiome.</p>
<p>A healthy microbiome has a large diversity of species over highly diverged phylogenetic taxa (<xref ref-type="bibr" rid="B28">28</xref>). Divergent taxa fill unique roles within a system, but when closely related taxa increase in a biological system, these roles break down and a dysbiosis event occurs. Our data indicates that during the bottleneck event of the collapsed colony, a dysbiosis occurred where there was a major shift in the microbiome, either leading to the colony collapse or a factor of the collapse.</p>
<p>While the high mite infestation was competing for resources and emphasizes the need for frequent culturing of reared insects, we never found the mites directly feeding on <italic>T. molitor</italic> developmental stages. There are no references that indicate that mites can pass on pathogens to insects, but that possibility remains.</p>
<p>More data are needed on the natural and diseased profiles of insects reared for downstream food applications. Our efforts are in developing insect-specific pathogen databases that can rapidly and accurately identify insect pathogens harmful to farmed insects (disease diagnosis) as well as those that can identify potentially harmful pathogens in insect food and feed for downstream applications (disease prevention).</p>
</sec>
<sec id="s5" sec-type="data-availability">
<title>Data availability statement</title>
<p>The datasets presented in this study can be found in online repositories. The name of the repository and accession number can be found below: NCBI; PRJNA883751.</p>
</sec>
<sec id="s6" sec-type="author-contributions">
<title>Author contributions</title>
<p>DL, MO, BO were involved in all aspects of the project and manuscript preparation. All authors contributed to the article and approved the submitted version.</p>
</sec>
</body>
<back>
<sec id="s7" sec-type="funding-information">
<title>Funding</title>
<p>Funding was provided by ARS CRIS 3020-43000-034-000D.</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>We wish to thank Abby Perkins and Molly Edeburn for technical help. This research was conducted under Kansas State University Institutional Biosafety Permit #1456. 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.</p>
</ack>
<sec id="s8" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>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.</p>
</sec>
<sec id="s9" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>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.</p>
</sec>
<sec id="s10" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/finsc.2023.1059046/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/finsc.2023.1059046/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="DataSheet_1.zip" id="SF1" mimetype="application/zip">
<label>Additional File S1</label>
<caption>
<p>Details from the analysis of sequence reads in Kraken2 (<xref ref-type="bibr" rid="B21">21</xref>). <bold>(A)</bold> Number of of reads from the collapsed colony (C-Adult, C-Diet, C-Frass), the new colony (N-Adult, N-Pupa, N-Diet), the recovered colony (R-Larva), and fresh diet (F-Diet) that were unclassified or classified by taxa rank. Classification at the phylum level were color coded: bacteria, yellow; archea, grey; eukaryote, green; and virus, orange. <bold>(B)</bold> Percentage of reads classified at the superkingdom rank.</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="DataSheet_1.zip" id="SF2" mimetype="application/zip">
<label>Additional File S2</label>
<caption>
<p>Species identified from reads in this study, with all reads associated with the collapsed colony (C-Adult, C-Diet, C-Frass, &#x201c;Collapsed Colony Reads&#x201d;) combined, and those associated with the new colony (N-Adult, N-Pupa, N-Diet, &#x201c;New Colony Reads&#x201d;) also combined, with a comparison of the number of reads from the collapsed and new colony reads (&#x201c;Collapsed/New&#x201d;).</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="DataSheet_1.zip" id="SF3" mimetype="application/zip">
<label>Additional File S3</label>
<caption>
<p>Reads with similarity to the genus <italic>Bacillus</italic>, with <italic>B. thuringiensis</italic> serovar <italic>morrisoni</italic> 4AA1 reads in yellow. Blast e-value, Blast top hit taxonomy name, and Blast hit accession are provided.</p>
</caption>
</supplementary-material>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1</label>
<citation citation-type="book">
<person-group person-group-type="author">
<collab>Food and Agricultural Organization (FAO)</collab>
</person-group>. <source>How to feed the world in 2050</source>. <publisher-loc>Rome, Italy</publisher-loc>: <publisher-name>Food and Agriculture Organization of the United Nations</publisher-name> (<year>2009</year>).</citation>
</ref>
<ref id="B2">
<label>2</label>
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Van Huis</surname> <given-names>A</given-names>
</name>
<name>
<surname>Van Itterbeeck</surname> <given-names>J</given-names>
</name>
<name>
<surname>Klunder</surname> <given-names>H</given-names>
</name>
<name>
<surname>Mertens</surname> <given-names>E</given-names>
</name>
<name>
<surname>Halloran</surname> <given-names>A</given-names>
</name>
<name>
<surname>Muir</surname> <given-names>G</given-names>
</name>
<etal/>
</person-group>. <source>Edible insects: future prospects for food and feed security</source>. <publisher-loc>Rome, Italy</publisher-loc>: <publisher-name>Food and Agriculture Organization of the United Nations</publisher-name> (<year>2013</year>).</citation>
</ref>
<ref id="B3">
<label>3</label>
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Akhtar</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Isman</surname> <given-names>MB</given-names>
</name>
</person-group>. <article-title>Insects as an alternative protein source</article-title>. In: <person-group person-group-type="editor">
<name>
<surname>Yada</surname> <given-names>RY</given-names>
</name>
</person-group>, editor. <source>Proteins in food processing</source>. <publisher-loc>Duxford</publisher-loc>: <publisher-name>Woodhead Publishing</publisher-name> (<year>2018</year>).</citation>
</ref>
<ref id="B4">
<label>4</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Miglietta</surname> <given-names>PP</given-names>
</name>
<name>
<surname>De Leo</surname> <given-names>F</given-names>
</name>
<name>
<surname>Ruberti</surname> <given-names>M</given-names>
</name>
<name>
<surname>Massari</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>Mealworms for food: A water footprint perspective</article-title>. <source>Water</source> (<year>2015</year>) <volume>7</volume>(<issue>11</issue>):<page-range>6190&#x2013;203</page-range>. doi: <pub-id pub-id-type="doi">10.3390/w7116190</pub-id>
</citation>
</ref>
<ref id="B5">
<label>5</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Grau</surname> <given-names>T</given-names>
</name>
<name>
<surname>Vilcinskas</surname> <given-names>A</given-names>
</name>
<name>
<surname>Joop</surname> <given-names>G</given-names>
</name>
</person-group>. <article-title>Sustainable farming of the mealworm tenebrio molitor for the production of food and feed</article-title>. <source>Z F&#xfc;r Naturforschung C</source> (<year>2017</year>) <volume>72</volume>(<issue>9-10</issue>):<page-range>337&#x2013;49</page-range>. doi: <pub-id pub-id-type="doi">10.1515/znc-2017-0033</pub-id>
</citation>
</ref>
<ref id="B6">
<label>6</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Maciel-Vergara</surname> <given-names>G</given-names>
</name>
<name>
<surname>Jensen</surname> <given-names>AB</given-names>
</name>
<name>
<surname>Eilenberg</surname> <given-names>J</given-names>
</name>
</person-group>. <article-title>Cannibalism as a possible entry route for opportunistic pathogenic bacteria to insect hosts, exemplified by pseudomonas aeruginosa, a pathogen of the giant mealworm zophobas morio</article-title>. <source>Insects</source> (<year>2018</year>) <volume>9</volume>(<issue>3</issue>):<fpage>88</fpage>. doi: <pub-id pub-id-type="doi">10.3390/insects9030088</pub-id>
</citation>
</ref>
<ref id="B7">
<label>7</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dupriez</surname> <given-names>F</given-names>
</name>
<name>
<surname>Rejasse</surname> <given-names>A</given-names>
</name>
<name>
<surname>Rios</surname> <given-names>A</given-names>
</name>
<name>
<surname>Lefebvre</surname> <given-names>T</given-names>
</name>
<name>
<surname>Nielsen-LeRoux</surname> <given-names>C</given-names>
</name>
</person-group>. <article-title>Impact and persistence of serratia marcescens in tenebrio molitor larvae and feed under optimal and stressed mass rearing conditions</article-title>. <source>Insects</source> (<year>2022</year>) <volume>13</volume>(<issue>5</issue>):<fpage>458</fpage>. doi: <pub-id pub-id-type="doi">10.3390/insects13050458</pub-id>
</citation>
</ref>
<ref id="B8">
<label>8</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rodriguez</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Omoto</surname> <given-names>CK</given-names>
</name>
<name>
<surname>Gomulkiewicz</surname> <given-names>R</given-names>
</name>
</person-group>. <article-title>Individual and population effects of eugregarine, gregarina niphandrodes (Eugregarinida: Gregarinidae), on tenebrio molitor (Coleoptera: Tenebrionidae)</article-title>. <source>Environ Entomol</source> (<year>2007</year>) <volume>36</volume>(<issue>4</issue>):<page-range>689&#x2013;93</page-range>. doi: <pub-id pub-id-type="doi">10.1603/0046-225X(2007)36[689:IAPEOE]2.0.CO;2</pub-id>
</citation>
</ref>
<ref id="B9">
<label>9</label>
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Clopton</surname> <given-names>RE</given-names>
</name>
<name>
<surname>Janovy</surname> <given-names>JJ</given-names> <suffix>Jr</suffix>
</name>
</person-group>. <article-title>Developmental niche structure in the gregarine assemblage parasitizing <italic>Tenebrio molitor</italic>
</article-title>. <source>J Parasitology</source> <volume>79</volume>:<page-range>701&#x2013;9</page-range>.</citation>
</ref>
<ref id="B10">
<label>10</label>
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Lestari</surname> <given-names>AS</given-names>
</name>
<name>
<surname>Rao</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>Laboratory bioassays of metarhizium spp. and beauveria spp. against tenebrio molitor larvae</article-title>. In: <source>The 6th international symposium for sustainable humanosphere, humanosphere science school</source> <publisher-loc>Bangor, Indonesia</publisher-loc>: <publisher-name>Humanosphere Science School</publisher-name> (<year>2016</year>). p. <page-range>15&#x2013;6</page-range>.</citation>
</ref>
<ref id="B11">
<label>11</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shapiro-Ilan</surname> <given-names>D</given-names>
</name>
<name>
<surname>Rojas</surname> <given-names>MG</given-names>
</name>
<name>
<surname>Morales-Ramos</surname> <given-names>JA</given-names>
</name>
<name>
<surname>Lewis</surname> <given-names>EE</given-names>
</name>
<name>
<surname>Tedders</surname> <given-names>WL</given-names>
</name>
</person-group>. <article-title>Effects of host nutrition on virulence and fitness of entomopathogenic nematodes: Lipid-and protein-based supplements in tenebrio molitor diets</article-title>. <source>J Nematol</source> (<year>2008</year>) <volume>40</volume>(<issue>1</issue>):<page-range>13&#x2013;9</page-range>.</citation>
</ref>
<ref id="B12">
<label>12</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fabrick</surname> <given-names>J</given-names>
</name>
<name>
<surname>Oppert</surname> <given-names>C</given-names>
</name>
<name>
<surname>Lorenzen</surname> <given-names>MD</given-names>
</name>
<name>
<surname>Oppert</surname> <given-names>B</given-names>
</name>
<name>
<surname>Jurat-Fuentes</surname> <given-names>JL</given-names>
</name>
</person-group>. <article-title>A novel tenebrio molitor cadherin is a functional receptor for <italic>Bacillus thuringiensis</italic> toxin Cry3Aa</article-title>. <source>J Biol Chem</source> (<year>2009</year>) <volume>284</volume>:<page-range>18401&#x2013;10</page-range>. doi: <pub-id pub-id-type="doi">10.1074/jbc.M109.001651</pub-id>
</citation>
</ref>
<ref id="B13">
<label>13</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oppert</surname> <given-names>B</given-names>
</name>
</person-group>. <article-title>Rapid bioassay to screen potential biopesticides in <italic>Tenebrio molitor</italic> larvae</article-title>. <source>Biopestic Int</source> (<year>2010</year>) <volume>6</volume>:<fpage>67</fpage>&#x2013;<lpage>73</lpage>.</citation>
</ref>
<ref id="B14">
<label>14</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oppert</surname> <given-names>B</given-names>
</name>
<name>
<surname>Morgan</surname> <given-names>T</given-names>
</name>
<name>
<surname>Kramer</surname> <given-names>KJ</given-names>
</name>
</person-group>. <article-title>Efficacy of <italic>Bacillus thuringiensis</italic> Cry3Aa protoxin and protease inhibitors toward coleopteran storage pests</article-title>. <source>Pest Manage Sci</source> (<year>2011</year>) <volume>67</volume>:<page-range>568&#x2013;73</page-range>. doi: <pub-id pub-id-type="doi">10.1002/ps.2099</pub-id>
</citation>
</ref>
<ref id="B15">
<label>15</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oppert</surname> <given-names>B</given-names>
</name>
<name>
<surname>Dowd</surname> <given-names>SE</given-names>
</name>
<name>
<surname>Bouffard</surname> <given-names>P</given-names>
</name>
<name>
<surname>Li</surname> <given-names>L</given-names>
</name>
<name>
<surname>Conesa</surname> <given-names>A</given-names>
</name>
<name>
<surname>Lorenzen</surname> <given-names>MD</given-names>
</name>
<etal/>
</person-group>. <article-title>Transcriptome profiling of the intoxication response of tenebrio molitor larvae to <italic>Bacillus thuringiensis</italic> Cry3Aa protoxin</article-title>. <source>PloS One</source> (<year>2012</year>):<elocation-id>e34624</elocation-id>. doi: <pub-id pub-id-type="doi">10.1371/journal.pone.0034624</pub-id>
</citation>
</ref>
<ref id="B16">
<label>16</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oppert</surname> <given-names>B</given-names>
</name>
<name>
<surname>Martynov</surname> <given-names>AG</given-names>
</name>
<name>
<surname>Elpidina</surname> <given-names>EN</given-names>
</name>
</person-group>. <article-title>
<italic>Bacillus thuringiensis</italic> Cry3Aa intoxication of <italic>Tenebrio molitor</italic> induces widespread changes in the expression of serine peptidase transcripts</article-title>. <source>Comp Biochem Physiol</source> (<year>2012</year>) <volume>7D</volume>:<page-range>233&#x2013;42</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.cbd.2012.03.005</pub-id>
</citation>
</ref>
<ref id="B17">
<label>17</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vigneron</surname> <given-names>A</given-names>
</name>
<name>
<surname>Jehan</surname> <given-names>C</given-names>
</name>
<name>
<surname>Rigaud</surname> <given-names>T</given-names>
</name>
<name>
<surname>Moret</surname> <given-names>Y</given-names>
</name>
</person-group>. <article-title>Immune defenses of a beneficial pest: the mealworm beetle, tenebrio molitor</article-title>. <source>Front Physiol</source> (<year>2019</year>) <volume>10</volume>:<elocation-id>138</elocation-id>. doi: <pub-id pub-id-type="doi">10.3389/fphys.2019.00138</pub-id>
</citation>
</ref>
<ref id="B18">
<label>18</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mathew</surname> <given-names>AG</given-names>
</name>
<name>
<surname>Cissell</surname> <given-names>R</given-names>
</name>
<name>
<surname>Liamthong</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>Antibiotic resistance in bacteria associated with food animals: a united states perspective of livestock production</article-title>. <source>Foodborne Pathog Dis</source> (<year>2007</year>) <volume>4</volume>(<issue>2</issue>):<page-range>115&#x2013;33</page-range>. doi: <pub-id pub-id-type="doi">10.1089/fpd.2006.0066</pub-id>
</citation>
</ref>
<ref id="B19">
<label>19</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Maciel-Vergara</surname> <given-names>G</given-names>
</name>
<name>
<surname>Jensen</surname> <given-names>AB</given-names>
</name>
<name>
<surname>Lecocq</surname> <given-names>A</given-names>
</name>
<name>
<surname>Eilenberg</surname> <given-names>J</given-names>
</name>
</person-group>. <article-title>Diseases in edible insect rearing systems</article-title>. <source>J Insects Food Feed</source> (<year>2021</year>) <volume>7</volume>:<page-range>621&#x2013;38</page-range>. doi: <pub-id pub-id-type="doi">10.3920/JIFF2021.0024</pub-id>
</citation>
</ref>
<ref id="B20">
<label>20</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Eilenberg</surname> <given-names>J</given-names>
</name>
<name>
<surname>Vlak</surname> <given-names>JM</given-names>
</name>
<name>
<surname>Nielsen-LeRoux</surname> <given-names>C</given-names>
</name>
<name>
<surname>Cappellozza</surname> <given-names>S</given-names>
</name>
<name>
<surname>Jensen</surname> <given-names>AB</given-names>
</name>
</person-group>. <article-title>Diseases in insects produced for food and feed</article-title>. <source>J Insects as Food Feed</source> (<year>2015</year>) <volume>1</volume>(<issue>2</issue>):<fpage>87</fpage>&#x2013;<lpage>102</lpage>. doi: <pub-id pub-id-type="doi">10.3920/JIFF2014.0022</pub-id>
</citation>
</ref>
<ref id="B21">
<label>21</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wood</surname> <given-names>DE</given-names>
</name>
<name>
<surname>Lu</surname> <given-names>J</given-names>
</name>
<name>
<surname>Langmead</surname> <given-names>B</given-names>
</name>
</person-group>. <article-title>Improved metagenomic analysis with kraken 2</article-title>. <source>Genome Biol</source> (<year>2019</year>) <volume>20</volume>(<issue>1</issue>):<fpage>1</fpage>&#x2013;<lpage>3</lpage>. doi: <pub-id pub-id-type="doi">10.1186/s13059-019-1891-0</pub-id>
</citation>
</ref>
<ref id="B22">
<label>22</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Donovan</surname> <given-names>WP</given-names>
</name>
<name>
<surname>Gonzalez</surname> <given-names>JM</given-names> <suffix>Jr</suffix>
</name>
<name>
<surname>Gilbert</surname> <given-names>MP</given-names>
</name>
<name>
<surname>Dankocsik</surname> <given-names>C</given-names>
</name>
</person-group>. <article-title>Isolation and characterization of EG2158, a new strain of <italic>Bacillus thuringiensis</italic> toxic to coleopteran larvae, and nucleotide sequence of the toxin gene</article-title>. <source>Mol Gen Genet</source> (<year>1988</year>) <volume>214</volume>:<page-range>365&#x2013;72</page-range>. doi: <pub-id pub-id-type="doi">10.1007/BF00330468</pub-id>
</citation>
</ref>
<ref id="B23">
<label>23</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Crickmore</surname> <given-names>N</given-names>
</name>
<name>
<surname>Berry</surname> <given-names>C</given-names>
</name>
<name>
<surname>Panneerselvam</surname> <given-names>S</given-names>
</name>
<name>
<surname>Mishra</surname> <given-names>R</given-names>
</name>
<name>
<surname>Connor</surname> <given-names>TR</given-names>
</name>
<name>
<surname>Bonning</surname> <given-names>BC</given-names>
</name>
</person-group>. <article-title>A structure-based nomenclature for bacillus thuringiensis and other bacteria-derived pesticidal proteins</article-title>. <source>J Invertebr Pathol</source> (<year>2021</year>) <volume>186</volume>:<elocation-id>107438</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jip.2020.107438</pub-id>
</citation>
</ref>
<ref id="B24">
<label>24</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nayduch</surname> <given-names>D</given-names>
</name>
<name>
<surname>Cho</surname> <given-names>H</given-names>
</name>
<name>
<surname>Joyner</surname> <given-names>C</given-names>
</name>
</person-group>. <article-title>
<italic>Staphylococcus aureus</italic> in the house fly: temporospatial fate of bacteria and expression of the antimicrobial peptide defensin</article-title>. <source>J Med Entomol</source> (<year>2013</year>) <volume>50</volume>:<page-range>171&#x2013;8</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1603/me12189</pub-id>
</citation>
</ref>
<ref id="B25">
<label>25</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mundt</surname> <given-names>JO</given-names>
</name>
</person-group>. <article-title>The ecology of the streptococci</article-title>. <source>Microb Ecol</source> (<year>1982</year>) <volume>8</volume>:<page-range>355&#x2013;69</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/BF02010675</pub-id>
</citation>
</ref>
<ref id="B26">
<label>26</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McFrederick</surname> <given-names>QS</given-names>
</name>
<name>
<surname>Cannone</surname> <given-names>JJ</given-names>
</name>
<name>
<surname>Gutell</surname> <given-names>RR</given-names>
</name>
<name>
<surname>Kellner</surname> <given-names>K</given-names>
</name>
<name>
<surname>Plowes</surname> <given-names>RM</given-names>
</name>
<name>
<surname>Mueller</surname> <given-names>UG</given-names>
</name>
</person-group>. <article-title>Specificity between lactobacilli and hymenopteran hosts is the exception rather than the rule</article-title>. <source>Appl Environ Microbiol</source> (<year>2013</year>) <volume>79</volume>:<page-range>1803&#x2013;12</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/AEM.03681-12</pub-id>
</citation>
</ref>
<ref id="B27">
<label>27</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Savio</surname> <given-names>C</given-names>
</name>
<name>
<surname>Mugo-Kamiri</surname> <given-names>L</given-names>
</name>
<name>
<surname>Upfold</surname> <given-names>JK</given-names>
</name>
</person-group>. <article-title>Bugs in bugs: The role of probiotics and prebiotics in maintenance of health in mass-reared insects</article-title>. <source>Insects</source> (<year>2022</year>) <volume>13</volume>:<elocation-id>376</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/insects13040376</pub-id>
</citation>
</ref>
<ref id="B28">
<label>28</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Robinson</surname> <given-names>CJ</given-names>
</name>
<name>
<surname>Bohannan</surname> <given-names>BJ</given-names>
</name>
<name>
<surname>Young</surname> <given-names>VB</given-names>
</name>
</person-group>. <article-title>From structure to function: the ecology of host-associated microbial communities</article-title>. <source>Microbiol Mol Biol Rev</source> (<year>2010</year>) <volume>74</volume>(<issue>3</issue>):<page-range>453&#x2013;76</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/MMBR.00014-10</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>