Investigation of the Dominant Microbiota in Ready-to-Eat Grasshoppers and Mealworms and Quantification of Carbapenem Resistance Genes by qPCR

In this study, 30 samples of processed edible mealworms (Tenebrio molitor L.) and 30 samples of grasshoppers (Locusta migratoria migratorioides) were obtained from producers located in Europe (Belgium and the Netherlands) and Asia (Thailand) and subjected to PCR-DGGE analyses. The PCR-DGGE analyses showed that species in the genus Staphylococcus were predominant in the samples of mealworms from Belgium and grasshoppers from the Netherlands; species in the genus Bacillus were detected in the samples of mealworms and grasshoppers from Thailand. Moreover, Weissella cibaria/confusa/spp. was found in grasshoppers from Belgium. Since data concerning the role of novel foods such as edible insects in the dissemination of carbapenem resistance are currently lacking, the quantification of five carbapenemase encoding genes (bla_NDM−1, bla_VIM, bla_GES, bla_OXA−48, and bla_KPC) by qPCR was also carried out in all the samples under study. The genes coding for GES and KPC were not detected in the analyzed samples. A very low frequency of bla_OXA−48 (3%) and bla_NDM−1 (10%) genes was detected among mealworms. In contrast, grasshoppers were characterized by a high incidence of the genes for OXA-48 and NDM-1, accounting for 57 and 27% of the overall grasshopper samples, respectively. The bla_VIM gene was detected exclusively in two grasshopper samples from Thailand, showing only 7% positivity. The analysis of variance showed that all the effects (producers, species, and producers × species) were statistically significant for bla_NDM−1, whereas for bla_OXA−48 and bla_VIM, no significant effects were detected for the same source of variation. Further studies are necessary to assess the possible role of edible insects as reservoirs for the resistance to carbapenems and to understand the correlation with the insect microbiota. Furthermore, an intensified surveillance plan examining the occurrence of carbapenemase encoding genes in the food chain and in environmental compartments is needed for a proper risk assessment. In such a context, the appropriate use of antimicrobials represents the main preventive action that should always be applied.

Introduction

The use of insects for human consumption is a frequent practice worldwide, mainly in Thailand and other Asian countries, Africa, America, and Australia (van Huis et al., 2013; Schlüter et al., 2017). In Europe, edible insects represent an innovative and uncommon protein source, although in some European countries, especially in the Netherlands and Belgium, the rearing and the industrial production of edible insects is gradually increasing (ANSES, 2014; Schlüter et al., 2017). Indeed, several promising aspects are associated with insect consumption since insects (i) are generally characterized by a positive nutrient profile in terms of high-quality proteins and amino acids, good lipids, vitamins, minerals, and fiber, (ii) are easy to breed, and (iii) cause lower emissions of greenhouse gases and ammonia than traditional livestock (Klunder et al., 2012; van Huis et al., 2013; Schlüter et al., 2017). Due to these numerous nutritional, social, and environmental benefits, edible insects are considered the “food of the future” and categorized as novel foods by Regulation (EU) No 2015/2283 of the European Parliament and of the Council.

Among the edible insects, mealworms and grasshoppers are included within those that are already commercialized as food in EU countries (Schlüter et al., 2017) and that have been partially investigated for the presence of relevant pathogens or potentially pathogenic microorganisms (Ali et al., 2010; Klunder et al., 2012; Stoops et al., 2016; Garofalo et al., 2017; Osimani et al., 2017a, 2018a) as well as for the presence of transferable resistances to antibiotics (Milanović et al., 2016; Osimani et al., 2017b,c, 2018b).

Carbapenems are broad-spectrum β-lactam antibiotics, currently considered the last-line antibiotics for the treatment of severe human infections caused by multidrug-resistant Gram-negative bacteria (EFSA BIOHAZ Panel, 2013; Guerra et al., 2014; Woodford et al., 2014). The production of carbapenemases, which are β-lactamases capable of hydrolyzing carbapenems and almost all β-lactams, represents the main mechanism of resistance to carbapenems (Tzouvelekis et al., 2012; Doi and Paterson, 2015; Fischer et al., 2017). Among the carbapenemases, the plasmid-acquired class A serine-β-lactamases KPC (Klebsiella pneumoniae carbapenemase, 17 variants) and GES (Guiana extended spectrum, 9 variants), class B metallo-β-lactamases VIM (Verona integron-encoded metallo-beta-lactamase, 40 variants) and NDM (New Delhi metallo-beta-lactamase, 10 variants), and class D serine-β-lactamases including OXA carbapenemases (Carbapenem-hydrolyzing oxacillinase) such as OXA-48 are among the most common and important from an epidemiological point of view and are thus the main clinical concern (Queenan and Bush, 2007; Grundmann et al., 2010; Miriagou et al., 2010; Walsh, 2010; Cantón et al., 2012; Monteiro et al., 2012; Pfeifer et al., 2012; Guerra et al., 2014; Woodford et al., 2014; Fischer et al., 2017). Carbapenemases are encoded by genes that are easily transferable among bacteria by horizontal gene transfer events since they are located on mobile genetic elements, thus increasing their worldwide spread among bacteria in different reservoirs (Tzouvelekis et al., 2012; Woodford et al., 2014; Doi and Paterson, 2015; Fischer et al., 2017). Indeed, during the last few years, the spread of carbapenem resistance has increased, especially in the Enterobacteriaceae, as well as in non-fermenters such as Pseudomonas spp. and Acinetobacter spp. and in non-pathogenic bacteria such as Stenotrophomonas spp. and Myroides spp. (Miriagou et al., 2003; Grundmann et al., 2010; Nordmann et al., 2011, 2012; Cantón et al., 2012; Tzouvelekis et al., 2012; Guerra et al., 2014; Doi and Paterson, 2015; Morrison and Rubin, 2015). In the last decade, the rapid and global dissemination of infections caused by carbapenemase-producing Enterobacteriaceae (CPE) in hospitals and healthcare institutions is of great concern since these outbreaks are often associated with high mortality rates due to the limited and inadequate alternative treatments (Nordmann et al., 2011; Doi and Paterson, 2015; Rossolini, 2015; Grundmann et al., 2017; Zhang et al., 2017). From all of these data, the current opinion is that the acquired carbapenemases are a primary pressing public health threat related to antibiotic resistance (AR) (Tzouvelekis et al., 2012; EFSA BIOHAZ Panel, 2013; Woodford et al., 2014; Doi and Paterson, 2015).

Although the occurrence of carbapenemases was first discovered and mainly investigated at hospitals and healthcare facilities, scientific studies reporting carbapenemase-producers, and carbapenemase-encoding genes (CEG) in animals, the environment and food are increasingly frequent. Specifically, carbapenem resistance has been detected in livestock and in their environments in France, Germany, Switzerland, the USA, and China (Guerra et al., 2014; Webb et al., 2016; Zurfluh et al., 2016; Fischer et al., 2017); companion animals and wildlife (Guerra et al., 2014; Woodford et al., 2014); aquatic environments (Zurfluh et al., 2013; Guerra et al., 2014; Woodford et al., 2014; Fernando et al., 2016); retail chicken meat from Egypt (Abdallah et al., 2015); vegetables and seafood from Asia, India and Brazil (Guerra et al., 2014; Morrison and Rubin, 2015; Zurfluh et al., 2015, 2016). This suggests that non-human sources may be reservoirs of CPE and CEG. Since it is widely recognized that the food chain is one of the main routes for the introduction of antibiotic-resistant bacteria and their genes into the human digestive tract, and for the diffusion and spread of AR in human pathogens (Clementi and Aquilanti, 2011; Rolain, 2013; Milanović et al., 2017), there is a need for intensified surveillance of the occurrence of CEG in the food chain and in different environmental compartments. This is also underscored by the European Food Safety Authority, which has recently recognized the need to improve European legislation to ensure the monitoring of carbapenem resistance in animals and food (EFSA BIOHAZ Panel, 2013).

Based on all these premises, edible insects deserve great attention in terms of safety, including the assessment of the microbiota and of the incidence of AR genes, in particular CEG.

Therefore, in order to obtain an overview of the predominant bacterial species in samples of processed edible mealworms (Tenebrio molitor L.) and grasshoppers (Locusta migratoria migratorioides) obtained from producers in Europe (Belgium and the Netherlands) and Asia (Thailand), the total microbial DNA was analyzed by culture-independent PCR-DGGE.

Concerning the AR issue, it is worth noting that currently, only a few scientific studies are available on the occurrence of transferable AR genes in edible insects (Milanović et al., 2016; Osimani et al., 2017b,c, 2018b; Vandeweyer et al., 2018), and to the authors' knowledge, none have been published on the detection of CEG.

To address this gap, the edible insect samples under study were subjected to screening by quantitative PCR (qPCR) of five among the most common carbapenem resistance genes (bla_NDM−1, bla_VIM, bla_GES, bla_OXA−48, and bla_KPC) (Monteiro et al., 2012). Statistical analyses were performed to determine if edible insect species (mealworms and grasshoppers) or geographical location correlated with the occurrence of carbapenem resistance genes in this study.

Materials and Methods

Sampling

Thirty samples of edible mealworms and 30 samples of grasshoppers (boiled, dried, and salted) were purchased via the internet from dealers located in Europe (Belgium and the Netherlands) and Asia (Thailand). Ten mealworm and ten grasshopper samples from each country were collected and marked as follows: TN1-TN10 (mealworms from the Netherlands, Producer 1), TB1-TB10 (mealworms from Belgium, Producer 2), TT1-TT10 (mealworms from Thailand, Producer 3), GN1-GN10 (grasshoppers from the Netherlands, Producer 1), GB1-GB10 (grasshoppers from Belgium, Producer 2), and GT1-GT10 (grasshoppers from Thailand, Producer 3). All the samples were provided in sealed plastic containers and delivered at ambient temperature via international shipping. No information was available on the rearing and hygiene conditions of processing, transport and storage applied to these edible insects before marketing.

Bacterial DNA Extraction

Total microbial DNA was extracted directly from the insect samples using PowerFood Microbial DNA Isolation Kit (Mo Bio Laboratories, Carlsbad, CA, USA) as described by Osimani et al. (2017a). The extracted DNA was quantified and checked for the purity using a NanoDrop ND 1000 (Thermo Fisher Scientific, Wilmington, DE, USA) and then standardized to 2 ng μL⁻¹ for qPCR assays and to 25 ng μL⁻¹ for PCR-DGGE analysis. The effective extraction of bacterial DNA was confirmed by conventional PCR amplification of 2 μL (50 ng) of extracted DNA suspensions in a My Cycler Thermal Cycler (BioRad Laboratories, Hercules, CA, USA) using universal prokaryotic primers 27F and 1495R as described by Osimani et al. (2015).

PCR-DGGE Analysis

The equal portions of DNAs extracted from insects were mixed together with the goal of obtaining six pooled samples (TB, TN, TT, GB, GN, and GT), each representing an insect type (mealworms and grasshoppers) and country of origin (Belgium, the Netherlands, and Thailand). The amplification products obtained from the 27F-1495R primer pair as described above were purified using the Illustra GFX PCR DNA and Gel Band Purification Kit (GE Healthcare Life Sciences, Buckinghamshire, UK). Subsequently, 2 μL of the purified PCR products was reamplified using the universal prokaryotic primers U968GC (added with GC clamp) and L1401 (Muyzer et al., 1993; Randazzo et al., 2002) following the PCR conditions previously described by Aquilanti et al. (2013). Following amplification, 5 μL of the PCR reaction was loaded on a 1.5% agarose gel together with a 100 bp molecular weight marker (HyperLadder™ 100 bp) to check for the expected product size of 480 bp. Twenty microliters of these PCR amplicons were analyzed by DGGE (30–60% urea-formamide denaturing gradient; 4 h at 130 V) using DCode Universal Mutation Detection System (Bio-Rad Laboratories) as described by Garofalo et al. (2015). All DGGE bands visible under UV light were excised from the gel, and the DNA was eluted overnight at 4°C in 50 μL of molecular biology grade water (Garofalo et al., 2008) and reamplified via PCR as described above, but with the forward primer U968 lacking the GC clamp. The PCR products were sent to Genewiz (Takeley, UK) for purification and sequencing, and the obtained sequences were identified at species level as described above by Osimani et al. (2018c).

Reference Strains

DNA extracted from five reference strains (Table 1), each carrying one of the carbapenem resistance genes under study, was used as positive control in the qPCR reactions and for the construction of qPCR standard curves.

TABLE 1

Table 1. Bacterial reference strains carrying carbapenems resistance genes, used as positive controls in the qPCR reactions.

qPCR

Absolute quantification of each carbapenemase gene (bla_NDM−1, bla_VIM, bla_GES, bla_OXA−48, and bla_KPC) in the insect samples was performed by qPCR in a Mastercycler^® ep realplex machine (Eppendorf, Hamburg, Germany) using the qPCR primers and cycling conditions described by Monteiro et al. (2012). To check for product specificity, all cycles were followed by a melt curve step analysis with temperature gradually increasing from 65 to 95°C by 0.2°C/s. Each qPCR reaction consisted of 4 μL (8 ng) of the extracted DNA; 5 μL of Type-it 2X HRM PCR Master Mix (Qiagen, Hilden, Germany) containing HotStarTaq Plus DNA Polymerase, EvaGreen Dye, an optimized concentration of Q-solution, dNTPs and MgCl₂; 0.2 μM of forward and reverse primers for each gene; and nuclease-free molecular biology grade water to a final reaction volume of 10 μL. The exogenous standards for each gene were prepared by qPCR amplification of the DNA extracted from the reference strains (Table 1) as described above but in a final reaction volume of 25 μL. The correct melting temperatures (Tm) and sizes of the obtained PCR products were checked by melting curve analysis and electrophoresis on a 1.5% agarose gel, respectively. The Illustra GFX PCR DNA and Gel Band Purification Kit (GE Healthcare Life Sciences) was used for the purification of the amplicons following the manufacturer's instructions. The quantity and the purity of the products were determined using the NanoDrop ND 1000 (Thermo Fisher Scientific). The calculation of each gene copy number was performed using an online calculator (www.idtdna.com) based on the mass and size of the purified qPCR products. The standard curves were created by qPCR amplification of 10-fold serial dilutions of exogenous standards. The amplification efficiencies were estimated from the slopes of the standard curves, and the correlation coefficients (R²) (Stolovitzky and Cecchi, 1996) were calculated automatically by Mastercycler^® ep realplex software. To determine the qPCR detection limit for each gene, the standard curves were generated in the range from ~1 to 10⁷ gene copies per reaction.

For the absolute quantification of CEG, the DNAs extracted from mealworms and grasshoppers were run along with the 10-fold serial dilutions of the standards prepared as described above. The gene copy number of each gene detected in the analyzed insect samples was determined using the slope of the corresponding standard curves. The baseline and threshold calculations were performed automatically by the Mastercycler^® ep realplex software. In addition to melting curve analysis, the correct sizes of the amplification products were checked by electrophoresis on 1.5% agarose gels using 100 bp DNA Ladder (HyperLadder™ 100 bp, Bioline, UK) as a molecular weight marker. Moreover, the accuracy of the amplification reactions was validated by the sequencing (Genewiz) of randomly selected positive samples (TB6, GT7, and GN1 for the bla_OXA−48 gene; GT2 and TT8 for the bla_NDM−1 gene; GT7 and GT8 for the bla_VIM gene). The resulting sequences were compared with those from the GenBank database (http://www.ncbi.nlm.nih.gov) using Basic Local Alignment Search Tool (BLAST) (Altschul et al., 1990). All qPCR reactions were performed in triplicate and the results were expressed as the mean gene copy number per ng of DNA ± standard deviation for each gene.

Statistical Analysis

Descriptive statistics, calculated on 20 samples for each producer and on 30 samples for each insect species were carried out for the bla_OXA−48, bla_NDM−1, and bla_VIM gene copies by computing the means ± standard deviation.

After first checking for conformance to a normal distribution and identification of outliers, analysis of variance (ANOVA) was carried out using JMP statistical software version 11.0.0 (SAS Institute Inc., Cary, NC, USA) to test the following main effects: producers (Belgium, Thailand, the Netherlands), insect species (mealworm, grasshoppers) and producer × species.

Principal component analysis (PCA) was applied to discriminate among mealworms and grasshoppers coming from different producers (Belgium, the Netherlands, and Thailand), and the presence of genes related to resistance to carbapenems (bla_NDM−1, bla_VIM, bla_GES, bla_OXA−48, and bla_KPC). PCA was carried out using the Unscrambler 7.5 software (CAMO ASA, Oslo, Norway). The mean data were normalized to neutralize any influence of hidden factors. PCA provides a graphical representation of the overall differences in terms of distribution of genes between the insects coming from different producers.

Results and Discussion

In the present study, the microbiota of commercialized read-to-eat grasshoppers and mealworms from different countries was investigated via PCR-DGGE, as well as the quantification and distribution of five common CEG within the same matrices, in order to have a more complete picture of some safety aspects related to edible insects.

Determination of Microbial Diversity

To have an overview of the predominant bacterial species found in the edible insects considered in this study, the total microbial DNA was extracted from the samples, the DNAs were mixed in order to obtain six pooled samples, each representing an insect type (mealworms, grasshoppers) and country of origin (Belgium, the Netherlands, Thailand) and then analyzed by a culture-independent PCR-DGGE method. The results obtained are reported in Table 2.

TABLE 2

Table 2. Sequencing results of the bands excised from the DGGE gel obtained from the amplified fragments of pooled bacterial DNA extracted directly from mealworms and grasshoppers.

The dominant species found in mealworms from Belgium and grasshoppers from the Netherlands belonged to the genus Staphylococcus, and species in the genus Bacillus were found in mealworms and grasshoppers from Thailand. Grasshoppers from Belgium were positive for Weissella cibaria/confusa/spp., while bacterial species with a percentage identity below 97% were found in mealworms from the Netherlands. It is interesting to note that the predominance of lactic acid bacteria (LAB), and in particular the Weissella spp., has been reported for processed and fresh grasshoppers from Belgium and the Netherlands (Stoops et al., 2016; Garofalo et al., 2017; Osimani et al., 2017a), thus suggesting that the specific rearing conditions may have selected for this microbial group or that this bacterial species is intrinsically associated with this edible insect.

The genera Staphylococcus and Bacillus identified among the other pooled samples of mealworms and grasshoppers may contain pathogenic species such as Staphylococcus aureus and Bacillus cereus, and these data are in agreement with other studies on the microbiota in fresh and processed mealworms and grasshoppers (Stoops et al., 2016; Garofalo et al., 2017). The presence of Staphylococcus members may be due to an environmental contamination occurring during human handling or processing. Indeed, these insects were boiled and salted, and since Staphylococcus spp. is a halophile bacterium usually predominating in environments with low microbial competition, it could have found conditions suitable for growth.

The lack of detection of bacterial DNA belonging to the Enterobacteriaceae has been already reported by Osimani et al. (2018c), although it is generally reported that Enterobacteriaceae represents a predominant bacterial group in the edible insect gut microbiota (Stoops et al., 2016; Garofalo et al., 2017; Osimani et al., 2017b). These data suggest that good manufacturing practices were applied during rearing and processing of the insects and/or that a degutting step may have been applied. A further explanation is that members of the Enterobacteriaceae family could be present within the processed insect samples but with a lower abundance in respect with other microbial groups. This latter hypotesis is supported by the fact that these microorganisms had previously found below the detection limit of microbial counts (<1 Log cfu g⁻¹) in the same samples (Osimani et al., 2017b,c).

Additionally, the lack of detection of bacterial DNA belonging to Pseudomonadaceae is unusual, but it is possible, as suggested for Enterobacteriaceae, that PCR-DGGE was not able to detect members of this microbial group if they were in the minority relative to the others.

Quantification and Distribution of carbapenem Resistance Genes

This study represents the first report on the screening of five carbapenemase encoding genes (bla_NDM−1, bla_VIM, bla_GES, bla_OXA−48, and bla_KPC) in processed edible mealworms (T. molitor L.) and edible grasshoppers (L. migratoria migratorioides) from producers located in Europe (Belgium and the Netherlands) and Asia (Thailand). The identification of the genes coding for these carbapenemases in the samples of edible insects under investigation was conducted by using qPCR. As previously underlined by Monteiro et al. (2012), molecular assays are considered the best solutions for the rapid detection of carbapenem resistance genes and for the identification of the resistance mechanism. The detection limit, defined as the lowest concentration at which 95% of the positive samples are detected was <10 gene copies per reaction for all the genes. The efficiencies of the qPCR reactions were 1.00 for the genes bla_OXA−48, bla_VIM, and bla_GES; 1.01 for the gene bla_KPC and 0.96 for the gene bla_NDM−1. The R² was 1.000 for the genes bla_OXA−48 and bla_NDM−1; 0.999 for the gene bla_GES, 0.996 for the gene bla_KPC and 0.995 for the gene bla_VIM. Moreover, the specificity of the primers used for the amplification of carbapenem resistance genes was confirmed by the results of the sequencing of randomly selected positive samples, which showed >97% similarity with the corresponding gene sequences deposited in the GenBank database. In more detail, the results of the BLAST analysis for the bla_OXA−48 gene showed 98% of the similarity with the sequences deposited in GenBank such as K. pneumoniae (CP031374), Citrobacter freundii (MG430338), Enterobacter ludwigii (MG436907), Pantoea agglomerans (MG436898), Escherichia coli (NG_055499), Shewanella sp. (NG_055475), and Proteus mirabilis (KT175900); for the bla_NDM−1 gene 99% similarity with the Enterobacter hormaechei (AP018835), K. pneumoniae (AP018834), Raoultella planticola (MH257689), E. coli (CP021206), Pseudomonas aeruginosa (KT364224), and Acinetobacter baumannii (KU180703); and for the bla_VIM gene 98% similarity with the K. pneumoniae (MH071811), C. freundii (NG_061412), Paenibacillus sp. (KR822172), E. coli (MF169879), E. hormaechei (LT991955), Kluyvera cryocrescens (MG228427), Alcaligenes faecalis (KY623659), Klebsiella oxytoca (NG_050362), and Enterobacter cloacae (CP030081).

The results of the qPCR quantification of carbapenem resistance genes in samples of ready-to-eat edible mealworms and grasshoppers produced in the Netherlands, Belgium and Thailand are reported in Tables 3, 4.

TABLE 3

Table 3. Results of qPCR quantification of carbapenemase genes in samples of ready-to-eat edible mealworms produced in the Netherlands (TN1-TN10), Belgium (TB1-TB10) and Thailand (TT1-TT10).

TABLE 4

Table 4. Results of qPCR quantification of carbapenemase genes in samples of ready-to-eat edible grasshoppers produced in the Netherlands (GN1-GN10), Belgium (GB1-GB10) and Thailand (GT1-GT10).

Regarding mealworms, none of the samples were positive for the genes bla_GES, bla_KPC, and bla_VIM, while only sample TB6 was found positive for bla_OXA−48 (3% positivity) and samples TN2, TN8, and TT8 were positive for bla_NDM−1 (10% positivity) (Table 3).

Regarding grasshoppers, the genes bla_GES and bla_KPC were not detected in any of the analyzed samples while only two samples from Thailand (GT7 and GT8) were positive for the bla_VIM gene (7% positivity). Interestingly, a high prevalence of bla_OXA−48 was noted (57% positivity), followed by bla_NDM−1 (27% positivity) (Table 4). Specifically, bla_OXA−48 was prevalent in 80% of the samples from Belgium, in 50% of the samples from the Netherlands and in 40% of the samples from Thailand. The highest frequency of bla_NDM−1 was found among samples from Thailand (40%), followed by samples from Belgium (30%) and the Netherlands (10%) (Table 4).

All the insect samples analyzed in this study were previously screened for the presence of 12 selected genes coding for resistance to tetracyclines [tet(M), tet(O), tet(S), and tet(K)], macrolide-lincosamide-streptogramin B (MLS_B) [erm(A), erm(B), erm(C)], vancomycin (vanA and vanB), beta-lactams (blaZ and mecA) and aminoglycosides [aac(6′)-Ie aph(2″)-Ia referred as aac-aph] through PCR and nested PCR assays (Osimani et al., 2017b,c). It is interesting to note that the mealworms samples TN2 and TN8 from the Netherlands were also found to be positive for the presence of genes coding for resistance to tetracyclines [tet(M), tet(K), tet(S)] and erythromycin [erm(B), erm(C)], while sample TB6 was positive for genes coding for resistance to tetracyclines [tet(M), tet(K)], and sample TT8 from Thailand was positive for genes coding for resistance to tetracyclines [tet(K)], erythromycin [erm(B)], and aminoglycosides (aac-aph) (Osimani et al., 2017b). Moreover, among grasshoppers, almost all of the samples that were found to be positive for bla_OXA−48, bla_NDM−1, and bla_VIM previously showed positivity for AR genes coding for resistance to tetracyclines [tet(M), tet(S), tet(K)], erythromycin [erm(B), erm(C)], aminoglycosides (aac-aph) and beta-lactams (blaZ) (Osimani et al., 2017c).

The average levels of gene copies ng⁻¹ in the 60 samples of edible insects were as follows: 0.59 ± 2.39 with a minimum value of 0 and a maximum value of 16.56 for bla_OXA−48, 2.94 ± 11.53 with a minimum value of 0 and a maximum value of 70.37 for bla_NDM−1, 7.25 ± 50.85 with a minimum value of 0 and a maximum value of 392.25 for bla_VIM.

Descriptive statistics on 20 samples from each producer are shown in Table 5. In addition, descriptive statistics on 30 samples from each insect species are reported in Table 6.

TABLE 5

Table 5. Descriptive statistics on 20 samples for each producer.

TABLE 6

Table 6. Descriptive statistics on 30 samples for each insect species.

The analysis of variance (Table 7) showed that all the variables (producers, species, and producers × species) had significant effects (P < 0.05) on the frequency of bla_NDM−1, whereas for bla_OXA−48 and bla_VIM, no significant effects were detected for the same source of variation.

TABLE 7

Table 7. ANOVA results for bla_OXA−48, bla_NDM−1, bla_VIM.

Regarding the distribution of bla_OXA−48 and bla_VIM, multiple comparisons between ACC Least Square Means (LSM) carried out using the Tukey test showed no significant differences among samples from different producers or insect species. Regarding bla_NDM−1, multiple comparisons (Tukey HSD) showed that insect species, but not the origin of the sample, had a significant correlation (P < 0.05) with the frequency of the gene.

PCA did not discriminate between the presence of genes encoding resistance to carbapenems among mealworms and grasshoppers coming from different producers. In contrast, in the previous studies on the occurrence of transferable ARs in ready-to-eat edible insects, PCA showed a differentiation among producers, thus suggesting that different rearing and clinical practices associated with different countries may have played a role in the variability observed (Milanović et al., 2016; Osimani et al., 2017b,c).

As reported by Schlüter et al. (2017), it is presumed that the rearing and processing conditions applied to edible insects will comply with the same food safety regulations as for livestock farming. The use of carbapenems is prohibited in food-producing animals in all countries (OIE, 2015; Webb et al., 2016). Notwithstanding, scientific studies reporting CPE and CEG in livestock and their environment are progressively more frequent (Guerra et al., 2014; Webb et al., 2016; Zurfluh et al., 2016; Fischer et al., 2017). Furthermore, in the last EFSA report on antimicrobial resistance in zoonotic and indicator bacteria from humans, animals and food in 2015, the presumptive extended-spectrum beta-lactamase (ESBL)-/AmpC-/carbapenemase-production in Salmonella and E. coli was monitored in humans, meat (pork and beef), fattening pigs and calves for the first time (EFSA and ECDC, 2017). Varying occurrence/prevalence rates of ESBL-/AmpC-producers were observed between countries, and carbapenemase-producing E. coli were detected in single samples of pig meat and from fattening pigs from two Member States (EFSA and ECDC, 2017). These data indicate that other antimicrobial classes could indirectly select CPE outside the hospital setting and that the rapid dissemination of CPE is also promoted by CEG located on plasmids transmissible by horizontal gene transfer events (Tzouvelekis et al., 2012; Woodford et al., 2014). As reviewed by Caniça et al. (2015), AR comprises a dynamic network that involves several environmental niches (e.g., water, soil, and plants) and different reservoirs (e.g., husbandry, hospitals, wild animal, human settings, human hand, food and global trade in foodstuffs) in which the path of dissemination and dynamics of AR genes has to be taken into consideration in order to understand and prevent the AR transmission and spread. Therefore, it is possible to hypothesize that, irrespective of the use of carbapenems in the edible insect rearing, the CEG may derive from the substrates used for feed or from surfaces and hands of operators or from treatments applied for processing, in addition to transport and storage. It is also interesting to note that grasshoppers and mealworms have different dietary habits since grasshoppers are grass-feeders whereas mealworms are usually reared on cereal-based matrices; therefore, the differences in terms of presence and distribution of CEG among these insect species may derive from different rearing practices and substrates.

Conclusion

Edible insects such as grasshoppers and mealworms represent a novel food that deserves attention in terms of safety, including the assessment of the incidence of AR genes. The investigation of the microbiota of the mealworm and grasshopper samples in this study revealed the presence of potential pathogenic and non-pathogenic species.

Scientific studies reporting carbapenemase-producing microorganisms and CEG in animals, the environment and food are increasingly frequent. The data presented in this study is the first attempt aimed at determining the incidence of CEG among samples of commercialized ready-to-eat grasshoppers and mealworms from Belgium, the Netherlands and Thailand. Although further studies are necessary to understand the correlation of CEG with the insect microbiota and to assess the possible role of edible insects as reservoirs of resistance to carbapenems, an intensified surveillance plan examining the occurrence of CEG in the food chain and in different environmental compartments, along with a prudent use of carbapenems and antimicrobials in general, are primary measures that should be applied.

Author Contributions

VM and AR carried out molecular analyses. MP, ST, LC, and MC performed statistical analyses. CG wrote the manuscript. AO, LA, CV, and FC critically analyzed the results and revised the final manuscript.

Funding

The study was partly supported by funding for the project “PSA2017-Carbapenemase-producing bacteria: from the environment to humans or vice versa?” from the Polytechnic University of Marche.

Conflict of Interest Statement

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.

Acknowledgments

We wish to thank the Esoteric Testing Laboratory, in the Pathology Department at Tampa General Hospital (Tampa, FL, USA) for providing the DNA extracted from reference strains, each carrying one of the carbapenem resistance genes studied in this research.

References

Abdallah, H. M., Reuland, E. A., Wintermans, B. B., Al Naiemi, N., Koek, A., Abdelwahab, A. M., et al. (2015). Extended-spectrum β-Lactamases and/or carbapenemases-producing enterobacteriaceae isolated from retail chicken meat in zagazig, Egypt. PLoS ONE 10:e0136052 doi: 10.1371/journal.pone.0136052

PubMed Abstract | CrossRef Full Text | Google Scholar

Ali, A., Mohamadou, B. A., Saidou, C., Aoudou, Y., and Tchiegang, C. (2010). Physicochemical properties and safety of grasshoppers, important contributors to food security in the far North Region of Cameroon. Res. J. Anim. Sci. 4, 108–111. doi: 10.3923/rjnasci.2010.108.111

CrossRef Full Text | Google Scholar

Altschul, S. F., Gish, W., Miller, W., Myers, E. W., and Lipman, D. J. (1990). Basic local alignment search tool. J. Mol. Biol,. 215, 403–410. doi: 10.1016/S0022-2836(05)80360-2

PubMed Abstract | CrossRef Full Text | Google Scholar

ANSES (2014). Opinion of the French Agency for Food, Environmental and Occupational Health & Safety on the Use of Insects as Food and Feed and the Review of Scientific Knowledge on the Health Risks Related to the Consumption of Insects. Opinion Request No, 2014-SA-0153

Aquilanti, L., Santarelli, S., Babini, V., Osimani, A., and Clementi, F. (2013). Quality evaluation and discrimination of semi-hard and hard cheeses from the Marche region (Central Italy) using chemometric tools. Int. Dairy J. 29, 42–52. doi: 10.1016/j.idairyj.2012.11.001

CrossRef Full Text | Google Scholar

Caniça, M., Manageiro, V., Jones-Dias, D., Clemente, L., Gomes-Neves, E., Poeta, P., et al. (2015). Current perspectives on the dynamics of antibiotic resistance in different reservoirs. Res. Microbiol. 166, 594–600. doi: 10.1016/j.resmic.2015.07.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Cantón R, Akóva M, Carmeli Y, Giske CG, Glupczynski Y, Gniadkowski M. (2012). Rapid evolution and spread of carbapenemases among Enterobacteriaceae in Europe. Clin. Microbiol. Infect. 18, 413–431. doi: 10.1111/j.1469-0691.2012.03821.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Clementi, F., and Aquilanti, L. (2011). Recent investigations and updated criteria for the assessment of antibiotic resistance in food lactic acid bacteria. Anaerobe 6, 394–398. doi: 10.1016/j.anaerobe.2011.03.021

CrossRef Full Text | Google Scholar

Doi, Y., and Paterson, D. L. (2015). Carbapenemase-Producing Enterobacteriaceae. Semin. Respir. Crit. Care Med. 36, 74–84. doi: 10.1055/s-0035-1544208

CrossRef Full Text | Google Scholar

EFSA BIOHAZ Panel (2013). Scientific Opinion on Carbapenem resistance in food animal ecosystems. EFSA J. 11:3501. doi: 10.2903/j.efsa.2013.3501

CrossRef Full Text

EFSA and ECDC (2017). The European Union summary report on antimicrobial resistance inzoonotic and indicator bacteria from humans, animals and food in 2015. EFSA J. 15:4694. doi: 10.2903/j.efsa.2017.4694

CrossRef Full Text

Fernando, D. M., Tun, H. M., Poole, J., Patidar, R., Li, R., Mi, R., et al. (2016). Detection of antibiotic resistance genes in source and drinking water samples from a first nations community in Canada. Appl. Environ. Microbiol. 82, 4767–4775 doi: 10.1128/AEM.00798-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Fischer, J., San José, M., Roschanski, N., Schmoger, S., Baumann, B., Irrgang, A., et al. (2017). Spread and persistence of VIM-1 Carbapenemase-producing Enterobacteriaceae in three German swine farms in 2011 and 2012. Vet. Microbiol. 200, 118–123. doi: 10.1016/j.vetmic.2016.04.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Garofalo, C., Osimani, A., Milanović, V., Aquilanti, L., De Filippis, F., Stellato, G., et al. (2015). Bacteria and yeast microbiota in milk kefir grains from different Italian regions. Food Microbiol. 49, 123–133. doi: 10.1016/j.fm.2015.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Garofalo, C., Osimani, A., Milanović, V., Taccari, M., Cardinali, F., Aquilanti, L., et al. (2017). The microbiota of marketed processed edible insects as revealed by high-throughput sequencing. Food Microbiol. 62, 15–22. doi: 10.1016/j.fm.2016.09.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Garofalo, C., Silvestri, G., Aquilanti, L., and Clementi, F. (2008). PCR-DGGE analysis of lactic acid bacteria and yeast dynamics during the production processes of three varieties of Panettone. J. Appl. Microbiol. 105, 243–254. doi: 10.1111/j.1365-2672.2008.03768.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Grundmann, H., Glasner, C., Albiger, B., Aanensen, D. M., Tomlinson, C. T., Andrasević, A. T., et al. (2017). Occurrence of carbapenemase-producing Klebsiella pneumoniae and Escherichia coli in the European survey of carbapenemase-producing Enterobacteriaceae (EuSCAPE): a prospective, multinational study. Lancet Infect. Dis. 17, 153–163. doi: 10.1016/S1473-3099(16)30257-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Grundmann, H., Livermore, D. M., Giske, C. G., Canton, R., Rossolini, G. M., Campos, J., et al. (2010). Carbapenem-nonsusceptible enterobacteriaceae in Europe: conclusions from a meeting of national experts. Euro Surveill. 15:19711. doi: 10.2807/ese.15.46.19711-en

CrossRef Full Text | Google Scholar

Guerra, B., Fischer, J., and Helmuth, R. (2014). An emerging public health problem: acquired carbapenemase-producing microorganisms are present in food-producing animals, their environment, companion animals and wild birds. Vet. Microbiol. 171, 290–297. doi: 10.1016/j.vetmic.2014.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Klunder, H. C., Wolkers-Rooijackers, J., Korpela, J. M., and Nout, M. J. R. (2012). Microbiological aspects of processing and storage of edible insects. Food Control. 26, 628–631. doi: 10.1016/j.foodcont.2012.02.013

CrossRef Full Text | Google Scholar

Milanović, V., Osimani, A., Aquilanti, L., Tavoletti, S., Garofalo, C., Polverigiani, S., et al. (2017). Occurrence of antibiotic resistance genes in the fecal DNA of healthy omnivores, ovo-lacto vegetarians and vegans. Mol. Nutr. Food Res. 61:1601098. doi: 10.1002/mnfr.201601098

PubMed Abstract | CrossRef Full Text | Google Scholar

Milanović, V., Osimani, A., Pasquini, M., Aquilanti, L., Garofalo, C., Taccari, M., et al. (2016). Getting insight into the prevalence of antibiotic resistance genes in specimens of marketed edible insects. Int. J. Food Microbiol. 227, 22–28. doi: 10.1016/j.ijfoodmicro.2016.03.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Miriagou, V., Cornaglia, G., Edelstein, M., Galani, I., Giske, C. G., Gniadkowski, M., et al. (2010). Acquired carbapenemases in Gram-negative bacterial pathogens: detection and surveillance issues. Clin. Microbiol. Infect. 16, 112–122. doi: 10.1111/j.1469-0691.2009.03116.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Miriagou, V., Tzouvelekis, L. S., Rossiter, S., Tzelepi, E., Angulo, F. J., and Whichard, J. M. (2003). Imipenem resistance in a Salmonella clinical strain due to plasmid-mediated class A carbapenemase KPC-2. Antimicrob. Agents Chemother. 47, 1297–1300.

PubMed Abstract | Google Scholar

Monteiro, J., Widen, R. H., Pignatari, A. C., Kubasek, C., and Silbert, S. (2012). Rapid detection of carbapenemase genes by multiplex real-time, PCR. J. Antimicrob. Chemother. 67, 906–909. doi: 10.1093/jac/dkr563

PubMed Abstract | CrossRef Full Text | Google Scholar

Morrison, B. J., and Rubin, J. E. (2015).Carbapenemase producing bacteria in the food supply escaping detection. PLoS ONE. 10:e0126717. doi: 10.1371/journal.pone.0126717

CrossRef Full Text

Muyzer, G., de Waal, E. C., and Uitterlinden, A. G. (1993). Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Appl. Environ. Microbiol. 59, 695–700.

PubMed Abstract | Google Scholar

Nordmann, P., Gniadkowski, M., Giske, C. G., Poirel, L., Woodford, N., and Miriagou, V. (2012). European network on carbapenemases, identification and screening of carbapenemase-producing enterobacteriaceae. Clin. Microbiol. Infect. 18, 432–438. doi: 10.1111/j.1469-0691.2012.03815.x

CrossRef Full Text | Google Scholar

Nordmann, P., Naas, T., and Poirel, L. (2011). Global spread of carbapenemase producing enterobacteriaceae. Emerg. Infect. Dis. 17, 1791–1798. doi: 10.3201/eid1710.110655

PubMed Abstract | CrossRef Full Text | Google Scholar

OIE (2015). List of Antimicrobial Agents of Veterinary Medicine. World Organization for Animal Health. Available online at: http://www.oie.int/fileadmin/Home/eng/Our_scientific_expertise/docs/pdf/Eng_OIE_List_antimicrobials_May2015.pdf.

Osimani, A., Cardinali, F., Aquilanti, L., Garofalo, C., Roncolini, A., Milanović, V., et al. (2017b). Occurrence of transferable antibiotic resistances in commercialized ready-to-eat mealworms (Tenebrio molitor L.). Int. J. Food Microbiol. 263, 38–46. doi: 10.1016/j.ijfoodmicro.2017.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Osimani, A., Garofalo, C., Aquilanti, L., Milanović, V., and Clementi, F. (2015). Unpasteurised commercial boza as a source of microbial diversity. Int. J. Food Microbiol. 194, 62–70. doi: 10.1016/j.ijfoodmicro.2014.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Osimani, A., Garofalo, C., Aquilanti, L., Milanović, V., Cardinali, F., Taccari, M., et al. (2017c). Transferable antibiotic resistances in marketed edible grasshoppers (Locusta migratoria migratorioides). J. Food Sci. 82, 1184–1192. doi: 10.1111/1750-3841.13700

PubMed Abstract | CrossRef Full Text | Google Scholar

Osimani, A., Garofalo, C., Milanović, V., Taccari, M., Cardinali, F., Aquilanti, L., et al. (2017a). Insight into the proximate composition and microbial diversity of edible insects marketed in the European Union. Eur. Food Res. Technol. 243, 1157–1171. doi: 10.1007/s00217-016-2828-4

CrossRef Full Text | Google Scholar

Osimani, A., Milanović, V., Cardinali, F., Garofalo, C., Clementi, F., Pasquini, M., et al. (2018a). The bacterial biota of laboratory-reared edible mealworms (Tenebrio molitor L.): From feed to frass. Int. J. Food Microbiol. 272, 49–60. doi: 10.1016/j.ijfoodmicro.2018.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Osimani, A., Milanović, V., Garofalo, C., Cardinali, F., Roncolini, A., Sabbatini, R., et al. (2018c). Revealing the microbiota of marketed edible insects through PCR-DGGE, metagenomic sequencing and real-time PCR. Int. J. Food Microbiol. 276, 54–62. doi: 10.1016/j.ijfoodmicro.2018.04.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Osimani, A., Milanović, V., Cardinali, F., Garofalo, C., Clementi, F., Ruschioni, S., et al. (2018b). Distribution of transferable antibiotic resistance genes in laboratory-reared edible mealworms (Tenebrio molitor L.). Front. Microbiol. 9:2702. doi: 10.3389/fmicb.2018.02702

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfeifer, Y., Schlatterer, K., Engelmann, E., Schiller, R. A., Frangenberg, H. R., Stiewe, D., et al. (2012). Emergence of OXA-48- type carbapenemase-producing Enterobacteriaceae in German hospitals. Antimicrob. Agents Chemother. 56, 2125–2128. doi: 10.1128/AAC.05315-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Queenan, A. M., and Bush, K. (2007). Carbapenemases: the Versatile β-Lactamases. Clin. Microbiol. Rev. 20, 440–458. doi: 10.1128/CMR.00001-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Randazzo, C. L., Torriani, S., Akkermans, A. D., de Vos, W. M., and Vaughan, E. E. (2002). Diversity, dynamics, and activity of bacterial communities during production of an artisanal Sicilian cheese as evaluated by 16S rRNA analysis. Appl. Environ. Microbiol. 68, 1882–1892. doi: 10.1128/AEM.68.4.1882-1892.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Rolain, J. M. (2013). Food and human gut as reservoirs of transferable antibiotic resistance encoding genes. Front. Microbiol. 4, 1–10:173. doi: 10.3389/fmicb.2013.00173

PubMed Abstract | CrossRef Full Text | Google Scholar

Rossolini, G. M. (2015). Extensively drug-resistant carbapenemase-producing Enterobacteriaceae: an emerging challenge for clinicians and healthcare systems. J. Int. Med. 277, 528–531. doi: 10.1111/joim.12350

PubMed Abstract | CrossRef Full Text | Google Scholar

Schlüter, O., Rumpold, B., Holzhauser, T., Roth, A., Vogel, R. F., Quasigroch, W., et al. (2017). Safety aspects of the production of foods and food ingredients from insects. Mol. Nutr. Food Res. 61:1600520. doi: 10.1002/mnfr.201600520

PubMed Abstract | CrossRef Full Text | Google Scholar

Stolovitzky, G., and Cecchi, G. (1996). Efficiency of DNA replication in the polymerase chain reaction. Proc. Natl. Acad. Sci. U.S.A. 93, 12947–12952. doi: 10.1073/pnas.93.23.12947

PubMed Abstract | CrossRef Full Text | Google Scholar

Stoops, J., Crauwels, S., Waud, M., Claes, J., Lievens, B., and Van Campenhout, L. (2016). Microbial community assessment of mealworm larvae (Tenebrio molitor) and grasshoppers (Locusta migratoria migratorioides) sold for human consumption. Food Microbiol. 53, 122–127. doi: 10.1016/j.fm.2015.09.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Tzouvelekis, L. S., Markogiannakis, A., Psichogiou, M., Tassios, P. T., and Daikos, G. L. (2012). Carbapenemases in Klebsiella pneumoniae and other enterobacteriaceae: an evolving crisis of global dimensions. Clin. Microbiol. Rev. 25, 682–707. doi: 10.1128/CMR.05035-11

PubMed Abstract | CrossRef Full Text | Google Scholar

van Huis, A., Van Itterbeeck, J., Klunder, H., Mertens, E., Halloran, A., and Vantomme, P. (2013). Edible Insects: Future Prospects for Food and Feed Security. Food and Agriculture Organization of the United Nations. Rome: FAO Forestry Paper, FAO. 187

Google Scholar

Vandeweyer, D., Milanović, V., Garofalo, C., Osimani, A., Clementi, F., Van Campenhout, L., et al. (2018). Real-time PCR detection and quantification of selected transferable antibiotic resistance genes in fresh edible insects from Belgium and the Netherlands. Int. J. Food Microbiol. 290, 288–295. doi: 10.1016/j.ijfoodmicro.2018.10.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Walsh, T. R. (2010). Emerging carbapenemases: a global perspective. Int. J. Antimicrob. Agents 36 S8–S14. doi: 10.1016/S0924-8579(10)70004-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Webb, H. E., Bugarel, M., den Bakker, H. C., Nightingale, K. K., Granier, S. A., Scott, H. M., et al. (2016). Carbapenem-resistant bacteria recovered from faeces of dairy cattle in the high plains region of the USA. PLoS ONE:e0147363. doi: 10.1371/journal.pone.0147363

PubMed Abstract | CrossRef Full Text | Google Scholar

Woodford, N., Wareham, D. W., Guerra, B., and Teale, C. (2014). Carbapenemase-producing enterobacteriaceae and non-enterobacteriaceae from animals and the environment: an emerging public health risk of our own making? J. Antimicrob. Chemother. 69, 287–291. doi: 10.1093/jac/dkt392

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, R., Chan, E. W., Zhou, H., and Chen, S. (2017). Prevalence and genetic characteristics of carbapenem-resistant enterobacteriaceae strains in China. Lancet. 17, 256–257. doi: 10.1016/S1473-3099(17)30072-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Zurfluh, K., Hächler, H., Nüesch-Inderbinen, M., and Stephan, R. (2013). Characteristics of extended-spectrum β-lactamase- and carbapenemase- producing enterobacteriaceae isolates from rivers and lakes in Switzerland. Appl. Environ. Microbiol. 79, 3021–3026. doi: 10.1128/AEM.00054-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Zurfluh, K., Hindermann, D., Nüesch-Inderbinen, M., Poirel, L., Nordmann, P., and Stephan, R. (2016). Occurrence and features of chromosomally encoded carbapenemases in Gram-negative bacteria in farm animals sampled at slaughterhouse level. Kurzmitteilungen 6, 457–460. doi: 10.17236/sat00072

CrossRef Full Text | Google Scholar

Zurfluh, K., Poirel, L., Nordmann, P., Klumpp, J., and Stephan, R. (2015). First detection of Klebsiella variicola producing OXA-181 carbapenemase in fresh vegetable imported from Asia to Switzerland. Antimicrob. Resist. Infect. Control. 4:38. doi: 10.1186/s13756-015-0080-5

PubMed Abstract | CrossRef Full Text | Google Scholar