Detection of Alphitobius diaperinus by Real-Time Polymerase Chain Reaction With a Single-Copy Gene Target

Use of edible insects as an alternative source of proteins in food and feed is increasing. These last years, numerous companies in Europe have started producing insects for food and feed purposes. In the European Union, the use of edible insects for human consumption falls within Regulation (EU) No. 2015/2283 on novel foods. For feed, Commission Regulation (EU) 2017/893 authorizes seven insect species as processed animal proteins for aquaculture. Methods of authentication are required to check the conformity of the products. In this study, we propose a real-time polymerase chain reaction (PCR) method for the specific detection of the lesser mealworm (Alphitobius diaperinus), one of the species included in the shortlist of authorized insects. The selected target is the cadherin gene with a single-copy (per haploid genome) illustrated by our experimental evidence. The PCR test amplified a 134-bp fragment of the cadherin gene. The qualitative method was assessed toward several performance criteria. Specificity was checked against 54 insect species next to other animal and plant species. The sensitivity, efficiency, robustness, and transferability of the PCR assay were also successfully tested. Finally, the applicability of the test was assessed on real-life processed samples (industrial meals) of A. diaperinus. The study also showed that there seems to be a huge confusion on the correct labeling of the marketed mealworms. We did not succeed to get Alphitobius laevigatus samples. They all appeared to belong to the A. diaperinus taxon.


INTRODUCTION
In recent years, edible insects are becoming an increasing alternative source of proteins. The main reason for it is that the food production rate is expected to be lower than the population growth (1). However, the insect consumption in European countries is limited because of strong impeding barriers associated with texture and appearance, much more than the taste (2,3). Therefore, the incorporation of insects as ingredients in common food items such as sausages, protein bars, pâté, buns, and pastas is a way to overcome this problem (3).
Insects represent a promising strategy for enriching food in some nutrients, thereby achieving a better nutritional balance (4). In fact, the positive nutritional features of edible insects are the presence of high quantities of proteins, essential and non-essential amino acids, lipids, fibers, vitamins, and minerals (4,5). Roncolini et al. (4) conducted a study to evaluate the use of lesser mealworm powder to replace a part of wheat flour as a means to enhance the protein and mineral content of crunchy snacks. The technological, microbiological, nutritional, and sensory characteristics of the fortified rusks were also evaluated in their study. A publication of Lacroix et al. (6) suggested the potential of lesser mealworm protein hydrolysates to serve as functional food ingredients to help improve glycemic regulation. Many insect species have been shown to have a high concentration not only of iron but also of zinc (3,7). The rearing substrate can then modulate the insect nutritional quality (3). Insect growth rates, production efficiencies, and protein quality are influenced by substrates used to feed the insects but also depend on entomological species raised and rearing conditions (8,9).
In animal production, the main protein source comes from either soybean (to feed cattle) or from fishmeal (to feed fish). However, soybean is associated with deforestation, use of genetically modified plant breeds, and massive use of pesticides, whereas industrial fishing to produce fishmeal depletes fish stocks (10). Insects are an ideal feed alternative because they are a good source of nutrients. Compared with conventional livestock, insects require less space for production (11), are expected to use less water (12), and emit less greenhouse gases (13), making them a more sustainable source of animal proteins (8).
Human consumption of edible insects in the European Union falls within Regulation (EU) No. 2015/2283 on novel foods (18). This means that food containing insects and their derived products must be subject to authorization once European Food Safety Authority (EFSA) has performed the risk assessments on their safety (4). In 2015, EFSA published a list of insect species potentially usable as food (19). Lesser mealworm (A. diaperinus) is included in this list. T. molitor larvae and A. domesticus are now assessed as safe novel food by EFSA (20)(21)(22), and dried larvae of the yellow mealworm are already authorized for placing on the market (23) with appropriate labeling of the food in which it is contained.
A. diaperinus is therefore used as feed for rearing animals and might be used for human consumption, once EFSA has given its approval for it as a safe novel food. A second species of the Alphitobius genus, Alphitobius laevigatus F., not listed by EFSA, is also marketed but it is only intended to feed house pets, including reptiles and amphibians.
Methods are required to check the conformity of these novel products. These last years, different methods of insect detection were developed. Mass spectrometry approaches were proposed for the detection of T. molitor and G. assimilis (24), H. illucens (25), A. domesticus, T. molitor, Locusta migratoria L., and A. diaperinus (25,26), as well as the Drosophila genus (Fallén) (27). Mass spectrometry may also help to differentiate insect meals according to their taxonomic groups (28). However, an accurate identification by proteomic methods will only be achievable after more intensive sequencing efforts, given the obvious lack of proteomic data for insect species (24). An adapted sedimentation protocol concentrating the insect particles for their detection in feed by light microscopy was proposed by Veys and Baeten (29).
At present, real-time polymerase chain reaction (PCR) remains the reference technique for DNA detection in food or feed products (30)(31)(32)(33)(34). Therefore, detection methods by real-time PCR were also proposed for insects. Most PCR tests published were developed on mitochondrial DNA. This facilitates detection in processed food and feed in which the DNA may be degraded, as a mitochondrion contains several copies of its genome, and several mitochondria can be present in a single cell (35). This is applied for H. illucens detection in feed (36,37), Oxya chinensis (Thunberg) in food (38), Bombyx mori L. in food (1) and feed (39), and G. sigillatus in food and feed (40). However, this multicopy characteristic is a disadvantage for quantitation purposes as the copy number per cell will be variable, depending on the considered tissue (36). This is why other publications focused on single-copy gene tests such as the real-time PCR test for T. molitor detection in food and feed (41) and B. mori in feed (39). In this study, we chose to focus on a chromosomal target that is a single copy per haploid genome.
This study proposes the first real-time PCR method for specific detection of lesser mealworm (A. diaperinus). The target used is the cadherin gene, which for A. diaperinus was characterized by Hua et al. (42). Its equivalent in T. molitor was already used as a target for the development of a real-time PCR test (41).

Samples
Insects were either collected in natural environment by trained entomologists, purchased from specialized companies, or provided by the Functional and Evolutionary Entomology laboratory of Gembloux Agro-Bio Tech (ULiège, Gembloux, Belgium). Insects were selected in order to cover several taxonomic groups or to have close relatives with a practical interest to the species considered. All insects were dead at the arrival at the laboratory. DNA extraction was mostly performed on a single individual, except for smaller insects for which several individuals were required. Real-life processed samples (four industrial meals from different EU-based companies) of A. diaperinus were obtained through the International Producers of Insects for Food and Feed (IPIFF), but the origin of samples must remain anonymous.
For each A. diaperinus industrial meal, a mix containing 0.1% (in mass fraction) of A. diaperinus in a commercial fish feed (composition: fishmeal, fish oil, wheat gluten, protein concentrate extracted from pea, maize starch, yeast, lecithin, vitamins, and minerals) was prepared.

DNA Extraction
Genomic DNA was extracted and purified from all samples following the CTAB-based method described in Annex A.3.1 of the international standard ISO 21571:2005 (43). Plasmid DNA was isolated from bacterial cultures with the help of the Genopure Plasmid Maxi Kit (Roche Diagnostics GmbH, Mannheim, Germany). The quality and quantity of DNA extracted from samples were estimated spectrophotometrically using a Nanodrop ND-1000 spectrophotometer at 260 nm (A260) and 280 nm (A280) absorbance. DNA purity was determined using the A260/A280 ratio. The amplifiability of the DNA extract was successfully checked by real-time PCR with the 18S target for insects (41), rbcL (44) for plants, and with a generic fish target (45) for fish. Other species were tested with targets developed or evaluated within the European Union Reference Laboratory for Animal Proteins in Feedingstuff 1 (EURL-AP, 2017) (32) or with the 18S target (46,47). Ten nanograms of DNA was used in the PCRs.
The industrial meals and the feed mixes at 0.1% of A. diaperinus were extracted following the method recommended by EURL-AP and based on the adaptation of the protocol of the Wizard Magnetic DNA Purification System for Food kit (Promega, Madison, WI, USA). This method is described in the EURL-AP Standard Operating Procedure 2 (EURL-AP, 2014). The quantities tested for this purpose are also in line with the EURL-AP SOP.

Primers and Probe for the Real-Time PCR
Eurogentec (Seraing, Belgium) synthesized the oligonucleotides. The primers and probe sequences developed for the detection of A. diaperinus are presented in Table 1. The probe for this latter method was labeled with the reporter dye FAM TM (6-carboxyfluorescein) at the 5 ′ end and the quencher dye TAMRA TM (tetramethyl-6-carboxyrhodamine) at the 3 ′ end.

Real-Time PCR Method
Real-time PCR (total reaction volume of 25 µL) was performed on a Lightcycler 480 (Roche Diagnostics Ltd., Rotkreuz, Switzerland) thermocycler using the Universal Mastermix provided by Diagenode (Seraing, Belgium). The reaction mixture included 12.5 µL of Master Mix, 1.7 µL of each primer (5 µM), 1.5 µL of probe (9 µM), 2.6 µL of bidistilled water, and 5 µL of DNA and distributed on 96-well reaction plates (Roche Diagnostics) for specific thermocyclers. Wells were covered with an adhesive film and centrifuged (2 min at 500 revolutions/min) to eliminate any air bubble in the well-bottoms. The thermal program was applied as follows: 2 min at 50 • C; 10 min at 95 • C; 50 cycles of 15 s at 95 • C and 1 min at 60 • C.

Specificity of the PCR Method
The specificity of the method was checked on 54 insect species of different taxonomic orders including 19 Coleoptera taxa other than A. diaperinus, 6 Diptera, 10 Orthoptera, 7 Hemiptera, 4 Hymenoptera, 5 Lepidoptera, 1 Neuroptera, and 1 Blattodea. One sample of larvae labeled buffalo worm and another sample of larvae considered as A. laevigatus were purchased from specialized companies. However, a Sanger sequencing trial (see Results) shows that these two samples belong to the A. diaperinus species and are therefore considered as such for the specificity test. The specificity tests were performed on two arachnids and six crustaceans, which like insects belong to the Arthropoda phylum, 1 mollusk and 33 species of vertebrates (12 terrestrial mammals, 6 sea mammals, 8 birds, 7 fish). The possibility of a cross-reaction with human DNA was also envisaged. Seven plant species frequently used in feed and food were included in the experimental setup ( Table 2). Ten nanograms of DNA was used in the PCRs. Each DNA extract was tested at least in duplicate.
*Sample labeled buffalo worm was purchased from a specialized company but identified as A. diaperinus by Sanger sequencing. **Sample labeled A. laevigatus was purchased from a specialized company but identified as A. diaperinus by Sanger sequencing. ***The Lucanus cervus (protected species) was not collected in the environment but obtained from an old insect box coming from a private collection. For positive samples, mean Cq values (m) and standard deviations (σ ) are given in brackets. Origin of the insect samples is specified with "a" for insects collected by trained entomologists, "b" for insects purchased from specialized companies and "c" for the insects provided by the Functional and Evolutionary Entomology lab of Gembloux Agro-Bio Tech (ULiège, Gembloux, Belgium).

Cloning of the Target, Copy Number Determination of the Plasmid DNA, and Dilutions
The 134-bp target from the A. diaperinus cadherin gene was ligated into the 3.9-kb pCR R 2.1-TOPO plasmid vector (Invitrogen, Merelbeke, Belgium) following the TOPO R TA Cloning R kit instructions (Invitrogen). Plasmid DNA isolated from bacterial cultures was linearized with the HindIII restriction enzyme (Promega) and then purified using phenol-chloroformisoamyl alcohol. The quantity of recovered plasmid DNA was converted into copy numbers as usual (36,48,49), taking into consideration that (i) 1 unit of absorbance at 260 nm corresponds to a concentration of 50 µg/mL of double-stranded DNA, and (ii) the mean molar weight of one base pair is set at 635 Da (50).
The sensitivity, efficiency, and robustness of the PCR test were determined on diluted plasmid DNA. These dilutions were performed in water until an estimated copy number of 10,000 copies/5 µL was reached. Higher dilutions of the target DNA were prepared in a solution containing 50 ng/µL of salmon sperm DNA as background DNA. Low binding tubes were chosen to minimize DNA losses.

Copy Number Determination of A. diaperinus Genomic DNA and Dilutions
The quantity of genomic DNA corresponding to 36,000 target copies was estimated at 10.08 ng for A. diaperinus based on data from the animal genome size database (https://www.genome.size.com/results.php?page=3) at the University of Guelph (Ontario, Canada). The sensitivity (limit of detection [LOD]) and efficiency of the PCR test were determined on diluted genomic DNA. Genomic DNA dilutions were performed in a similar way as described for the cloned target.

Limit of Detection
Target sensitivity was evaluated following the recommendations of the former AFNOR XP V03-020-2 standard (51). This standard no longer exists, but the principles detailed in it are still valid. The absolute LOD was determined for the PCR assay (primers + probe + amplification program) on dilutions of plasmid material and on dilutions of genomic material.
The subsequent dilutions had to contain 50, 20, 10, 5, 2, 1, and 0.1 copies of the target. Six PCRs had to be achieved for each dilution. The method's LOD 6 was the smallest copy number for which the six PCRs were positive, but only if the highest dilution supposed to contain the 0.1 copy per reaction generated a maximum of one positive PCR signal on six replicates. If more than one positive signal was observed for the 0.1 copy, the DNA quantities had to be revised. The copy number corresponding to LOD 6 was then tested 60 times on the same plate (determination of the LOD 95% ). The LOD 95% is validated as equal or below a given copy number if at least 95% of positive signals are recorded out of the 60 replicates. The highest acceptable copy number for LOD 6 and LOD 95% is 20 copies.

Efficiency
The efficiency of the PCR assay was calculated with a dilution series of genomic DNA and plasmid material at target levels of 5,000, 2,500, 1,000, 500, and 100 copies. Each dilution was analyzed in six replicates and on four runs. The efficiency has to be between 90 and 110% (52).

Digital PCR
The number of copies of the nuclear and plasmid DNA dilutions at approximately 500 copies/5 µL (5 µL being the volume of the DNA extracts added in the real-time PCR mix) were checked by digital PCR. Digital PCR was performed on the Biomark

Robustness of the PCR Method
The method robustness was tested by introducing some slight deviations to the standard experimental conditions (54). Parameters considered were as usual (36,52,55): the annealing temperature (60 ± 1 • C), the primer concentrations (standard or reduced by 30%), the probe concentration (standard or reduced by 30%), and the real-time PCR Master Mix volume (standard or ± 1 µL), which involves a final reaction volume of 25 ± 1 µL. Six replicates of the plasmid borne target at 20 copies/5 µL were tested in the conditions described in Table 3. The robustness was performed on two realtime PCR platforms: thermocycler Lightcycler 480 (Roche Diagnostics Ltd.) with Universal Mastermix by Diagenode and thermocycler QuantStudio TM 5 Real-Time PCR system (Applied Biosystems, Thermo Fisher Scientific, Foster City, CA, USA) with ABI TaqMan 2x Universal PCR Master Mix No AmpErase UNG (Applied Biosystems). The acceptance criterion is that all deviations to the standard protocol must give a positive result at a level of 20 copies of the target in the reaction (52).

Applicability of the PCR Method
The applicability of the PCR method was checked in duplicate on four real-life processed samples (industrial meals) of A. diaperinus produced in the EU and on a fish feed containing 0.1% in mass fraction of A. diaperinus industrial meal. Two DNA extracts with two dilutions were tested by PCR in duplicate.

Transferability of the PCR Method
The efficiency and the LOD (LOD 6 and LOD 95% ) of the PCR assay were tested on genomic DNA in the laboratory of Eurofins Biologie Moléculaire France with conditions similar to those of the developer's laboratory (CRA-W, Gembloux, Belgium). The dilution series for the LOD test were carried out by this second laboratory starting from the solution at 500 copies/5 µL checked by digital PCR.
Real-time PCR was performed on thermocycler CFX96 Deep Well Real-time PCR Detection Systems (Bio-rad, Hercules, CA, USA) using the Applied Biosystems TM TaqMan TM Universal PCR Master Mix (Applied Biosystems). Reaction mixtures were distributed on Hard-Shell R 96-Well PCR Plates (Bio-rad) developed for the CFX96 thermocyclers. Wells were covered with adhesive film, and the plates were centrifuged to eliminate any air bubble in the well-bottoms.

Amplicon Preparation for Sanger Sequencing
The PCRs to generate amplicons that had to be checked by Sanger sequencing were done as follows. For the large cadherin target, the PCR contained 1 µL of DNA extract from the sample to be checked, 6 µL of 5× GoTaq R Flexi Buffer (Promega), 3 µL Five microliters of PCR products were run on a 1.2% agarose gel to check amplicon quality. The remaining 25 µL of PCR products was sent to Eurofins Genomics (Constance, Germany) for Sanger sequencing.

RESULTS
The cadherin gene of A. diaperinus was used to select a piece of DNA that is specific to the considered insect species. Appropriate primers and probe were designed to amplify a 134-bp fragment of the cadherin gene ( Table 1). The latter is considered as a singlecopy gene in several insect species such as the Lepidopteran Ostrinia nubilalis (Hübner) (56) or the Coleoptera Diabrotica virgifera virgifera (Le Conte) (57), which is an advantage for quantitation purposes.
The specificity was first investigated in silico using the Blast tool available on the National Center for Biotechnology Information (NCBI) database. The different blasts and the alignments of DNA sequences with other Coleopteran (58) performed indicated that the PCR test should be specific to the target species A. diaperinus. The alignment with the closest relative within the NCBI database is the T. molitor sequence. It shows some similarity but not enough for allowing primers and probes to hybridize (Supplementary File 1). Specificity was also experimentally tested on DNA from A. diaperinus of various origins as well as on 53 other non-target insect species, including 19 Coleoptera. Positive results were obtained only with all samples of A. diaperinus. No signal was obtained with the 53 other insect species, the 42 other animal species (arachnids, crustaceans, mollusk and vertebrates), and the seven plant species tested ( Table 2).
A difficulty that appeared was to know if the test enables to distinguish A. diaperinus from A. laevigatus, a closely related species commercialized as feed product for non-farmed animals. Products labeled as A. laevigatus or as buffalo worm were tested, and a clear positive signal with the PCR test for A. diaperinus was observed. A doubt existed, however, if these samples really belonged to the A. laevigatus species. To check this hypothesis, Sanger sequencing was applied on PCR products of eight collected samples of Alphitobius (Supplementary File 2). Three targets were considered. The first one focused on a larger cadherin target ( Table 1), which completely contains the smaller target of 134 bp. The two other targets focused on COI regions known as more variable and considered as suitable to allow a distinction between A. diaperinus and A. laevigatus according to published data. All Sanger sequencings were successful and showed identical sequences for the eight samples (Supplementary Files 2-5).
The sequences obtained for the cadherin target corresponded to the A. diaperinus sequence published by Hua et al.  (Supplementary File 3). The alignment of sequenced COI fragments showed that all analyzed specimens belong to the A. diaperinus species when compared with A. diaperinus and A. laevigatus reference sequences. The first COI portion corresponded to the A. diaperinus sequence published by Hong et al. (2020-NC_049092.1) 4 (Supplementary File 4); this COI region is, however, not available for A. laevigatus. Reference sequences are available for both Alphitobius species for the second COI region considered. Results for that target showed that all samples of Alphitobius tested belonged to the A. diaperinus species (Supplementary File 5).
The amplification efficiency, LOD, and robustness were evaluated on plasmid DNA. The efficiency and LOD were also determined on genomic DNA.
To check that the number of copies in the dilutions used to assess the performance criteria was correct, the dilutions at approximately 500 copies/5 µL were estimated by digital PCR on a Biomark TM HD system. For the genomic DNA, the average obtained over the 10 measurements by digital PCR was 339 copies/5 µL with a variation coefficient at 9.45% ( Table 4). This mean copy number measured was slightly lower than the expected value (based on the genome size of A. diaperinus and considering the cadherin gene as a single-copy gene). The difference between these values did not exceed a factor of 2, corresponding to the subsequent dilution. These results therefore confirm that the cadherin gene is a single-copy gene per haploid genome as mentioned in other studies (56,57). However, in order to be closer to the expected values to evaluate the efficiency and LOD, new dilutions on the genomic DNA were carried out, taking into account the results obtained in digital PCR. On these new dilutions, the average obtained over the 22 measurements by digital PCR was 498 copies/5 µL with a variation coefficient at 9.21% ( Table 5). The mean copy number measured was estimated at the expected value and it is from this dilution series that the efficiency and sensitivity were evaluated.
The copy number of the plasmid material (linearized) was also checked, and the average copy obtained over the 22 measurements by digital PCR was 446 copies/5 µL with a variation coefficient at 9.26% ( Table 6). The mean copy number measured was estimated at the expected value, and it is from this dilution series that the efficiency, sensitivity, and robustness were evaluated.
The PCR efficiency was evaluated at 102.5% on plasmid DNA and 100.0% on genomic DNA. This was calculated, taking into account the mean Cq (quantification cycle) values obtained at the different copy numbers tested (from 5,000 to 10), and no outliers were encountered (Tables 7, 8). When calculated per plate, the efficiency was always higher than 90% and therefore met the acceptance criterion proposed by Broeders et al. (52). Concerning the sensitivity testing, the LOD 6 was estimated at five copies following the AFNOR XP V03-020-2 standard approach (51) and also at five copies for the LOD 95% with 57/60 positive signals on plasmid DNA and 58/60 positive signals on genomic DNA. Sixty of 60 positive signals were obtained at the level of 10 copies with a mean Cq value of 36.14 cycles on plasmid DNA and 36.78 cycles on genomic DNA. Therefore, the PCR test easily reaches the recommended performance criterion (≤ 20 copies).
The PCR method robustness was also evaluated on plasmid DNA, with success. All tested deviations to the standard protocol delivered positive results at the level of 20 copies in the PCR.
Positive signals were obtained on industrial samples (PAPs of A. diaperinus) showing the applicability of the PCR test on real-life samples ( Table 9).
The applicability was tested on a commercial fish feed adulterated with a low content of processed A. diaperinus. The commercial fish feed was first tested as free of A. diaperinus and the amplifiability of the DNA extracts obtained from this fish feed was checked with a 18S rDNA target (46,47).
The four mixes of fish feed containing 0.1% of different processed A. diaperinus meals were tested and gave positive results with the cadherin PCR test ( Table 9). The 10-fold dilutions provided evidence that there was no inhibitory effect of the feed matrix on the amplification of the A. diaperinus target.  36.14 ± 0.50 (60) For efficiency, each concentration was analyzed in six replicates and on four PCR plates (n = 24). For LOD 95% , the concentration at 10 copies was analyzed in 60 replicates on 1 PCR plate (n = 60).
Concerning the transferability, the PCR efficiency was evaluated at 90.0% on genomic DNA in a second laboratory and therefore met the acceptance criterion proposed by Broeders et al. (52). Table 10 indicates the mean Cq values obtained with the different copy numbers tested (from 5,000 to 10), and no outliers were encountered. When calculated per plate, the efficiency was higher than 90% for three plates and slightly below for the fourth one with an efficiency at 87.7%. The LOD 6 was estimated at five copies following the former AFNOR XP V03-020-2 standard approach (51) and at 10 copies for the LOD 95% with 60/60 positive signals. The mean Cq value at 10 copies is of 37.08 cycles. For efficiency, each concentration was analyzed in six replicates and on four PCR plates (n = 24). For LOD 95% , the concentration at 10 copies was analyzed in 60 replicates on 1 PCR plate (n = 60). The PCR test easily reached the recommended performance criterion (≤ 20 copies). The transferability of the method was therefore demonstrated.

DISCUSSION
The study describes a specific, sensitive, and robust test to detect A. diaperinus. With the recent authorization in the EU legislation to use eight insect species, among which lesser mealworm, in aquafeed, pig, and poultry feed, the interest of such a PCR test is increasing. For efficiency, each concentration was analyzed in 6 replicates and on 4 PCR plates (n = 24). For LOD 95% , the concentration at 10 copies was analyzed in 60 replicates on 1 PCR plate (n = 60).
During the search for the target and the validation study, special care was taken to be able to distinguish lesser mealworm from other representatives of the Tenebrionidae family (Supplementary File 1). Blasting the sequences of the developed primers and of the targeted fragment against the NCBI nucleotide database showed that the PCR test is specific to the target species A. diaperinus. Unfortunately, cadherin gene sequences for other Alphitobius species are not available in the NCBI and DNA Data Bank of Japan databases. Sequence alignments were therefore not possible. The other Alphitobius species are, however, not produced nor marketed, with the exception of A. laevigatus.
This confusion seems also to occur at the level of the sequences available in the NCBI database. Two COI sequences published by (72,73)  This probably also explains that the ordered samples of A. laevigatus were wrongly labeled by their providers. The Sanger sequencing revealed that these samples were in reality A. diaperinus samples. No real commercial sample of A. laevigatus was found, and even if we did not check with all possible sources, it seems that it is not that easy to find marketed A. laevigatus larvae. That is the reason why from a merely practical point of view we consider that the 134-bp cadherin target published in this study is specific to A. diaperinus. Nevertheless, from a scientific viewpoint, it is not impossible that the target is specific only at genus level. It could not be checked because of a lack of appropriate samples to do so.
The applicability of the PCR method was also successfully tested on real-world samples. This validates the implementation of the method on industrial samples.
Finally, it should be stressed that the target is present only once per haploid genome, which makes it a suitable target for possible quantification purposes. This would require to have a general insect target, present only once per haploid genome. To our knowledge, such a target has unfortunately not yet been identified.

CONCLUSIONS
The developed PCR method based on the cadherin gene fits for the purpose of detection of A. diaperinus in feed. Indeed, the efficiency met the required criterion, the specificity gave good results, and only A. diaperinus was detected with respect to the insect, animal, and plant species tested. The acceptance criteria were also reached for sensitivity (LOD 6 and LOD 95% ) and robustness. The PCR method was applicable on real-life samples from industry even when A. diaperinus was present at 0.1% in mass fraction in a fish feed. Finally, the transferability of the method in a second laboratory was also demonstrated by testing the efficiency and LOD.
The developed method was primarily aimed for application on feed. However, it is not excluded that several of the PCR tests developed in the study might be helpful to taxonomists in a better delineation of the several taxons within the genus Alphitobius.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
AM, HS, BD, and FD contributed to the design and implementation of the research. AM, HS, BD, JM, and FD designed and performed the experiments. SG and J-FM participated to the transferability study. AM, JM, BD, SG, J-FM, and FD analyzed the data and interpreted the results. AM wrote the manuscript with the help of FD. GB, BD, FF, SG, J-FM, and OF provided valuable comments to improve the quality of the manuscript. All authors contributed to the article and approved the submitted version.

FUNDING
This research was financially supported by the European Commission in the frame of Horizon 2020 Public-Private Partnership Bio-Based Industries Joint Undertaking (topic BBI.2018.F2 -Large-scale production of proteins for food and feed applications from alternative, sustainable sources) through the FARMYNG project.