Introduction

Flies of family calliphoridae (blow flies) are strong flying insects, highly mobile and typically one of the first flies reaching corpses within minutes after death (Goff et al. 1993; Kabadaia 2015). Accordingly, these flies may provide a useful solution in determining the minimum time since death which is known as the minimum post mortem interval (PMIMIN). Usually, forensic entomologists precisely calculate PMIMIN depending on species identification, duration of different stages at different temperatures, and other certain abiotic factors as geographic location, climate, latitude and so on (Turchetto and Vanin 2004; Sharif and Qamar 2021).

Usually, in death investigations, forensic professionals collect all insect specimens (adults, eggs, maggots and pupae) on corpses or around them. Sometimes, fresh insect samples are absent and only puparial exuviae (cases) are typically found. The exuvial identification of forensically important flies is so problematic, as they are often destroyed by the mechanical activity of adult emergence. Accordingly, traditional taxonomical identification of these deteriorated exuviae is very difficult. Also, natural degradation of DNA, enzymes and proteins during aging process deeply compromise molecular analysis of forensic samples (Gibbs and Crockettj 1998; Ye et al. 2007; Moore et al. 2021). One alternative method to identify and age insect species is cuticular hydrocarbon analysis (CHC). The external layer of insect exoskeleton is the cuticle which acts as a mechanical support, prevents desiccation, protects against different microorganisms and serves as a contact or close-range pheromone. It is composed of a mixture of esters, alcohols, ketones, aldehydes, fatty acids and hydrocarbons (Suarez et al. 2011; Sharif et al. 2023). In many insect species; hydrocarbons predominate the cuticle (Blomquist and Ginzel 2021) and proven to be very stable (Drijfhout 2010; Braga et al. 2016; Moore et al. 2022). Cuticular hydrocarbons are composed mainly of n-alkanes, branched methyl-alkanes, and alkenes (Blomquist and Bagnères 2010; Drijfhout 2010). The constituents of CHCs profile differ among various insect taxa; in the number of compounds, their proportions, chemical compositions, and chain lengths (Howard and Blomquist 2005; Sprenger and Menzel 2020; Sharif et al. 2023). These differences, along with certain other properties allowed the cuticular hydrocarbon content to be utilized for precise calculations of the weathering time of exuviae as well as the post mortem interval (Zhu et al. 2007; Sharif et al. 2023).

Many authors reported the uniqueness of CHCs and successfully identified many species (Anyanwu et al. 2000, 2001; Horne and Priestmann 2002; Bejarano et al. 2003; Shaalan et al. 2019; Moore et al. 2022). Also, studies confirmed that several CHCs exhibit changes in their relative abundance with chronological age of insect samples in different life stages (Urech et al. 2005; Zhu et al. 2006; Braga et al. 2016; Moore et al. 2021). In many insect species; CHCs function as sexual attractant pheromones or clues for species discrimination hence, very useful in speciation (Rundle et al. 2005). Changes in the composition of these clues alter mating preferences and pre-mating isolation and can be produced by changing the diet (Stennett and Etges 1997) or temperature (Buckley et al. 2003). This explains the divergence in the cuticular hydrocarbons of geographically isolated populations due to differences in food and/or temperature which may lead to reproductive isolation. Hence, CHCs represent better indicators of recent speciation events and reproductive isolation than other genetic and morphological characters, that require more time to be expressed after speciation events.

As far as we know this is the first study on Egyptian calliphorids investigating their cuticular hydrocarbon composition. Only few studies were done on other insect taxa as Hymenoptera (Surtasi et al. 2016; Elshaier 2021), Mantodea (Mohammad et al. 2009), or other dipteran species (Galhoum 2017, 2018; Shaalan et al. 2019). So, the aim of this preliminary study is to use the technique of gas chromatography/mass spectrometry (GC–MS) to analyze the chemical composition of the puparial exuviae of three widely distributed Egyptian blow flies of forensic relevance.

Materials and methods

Flies collection and identification

Stock colonies of Lucilia sericata, Chrysomya albiceps and Chrysomya marginalis were established from flies initially collected during May, June & July 2019 from El-Mansuryia, Giza Governorate and Cairo Governorate, Egypt. Collected adults were transferred to be reared in the Entomology laboratory, Zoology department, Zagazig University where they were maintained in rearing cages under laboratory conditions at (27˚C ± 2) and (55–70%) relative humidity.

Adults were provided with water, sugar and meat as oviposition media. Meat was supplied in a clear plastic cup with damp cotton piece to prevent drying out of meat and checked daily for oviposition. After that, each deposited egg batch was transferred to a new plastic jar containing fresh meat. Newly hatching larvae were transferred to new jars containing fresh meat, covered with muslin and fastened with rubber bands. Dry autoclaved sieved sawdust was used as a medium for pupation. The pupae were sieved from the sawdust and transferred in petri dishes to the rearing cages for adult emergence. After adult emergence, puparial exuviae were collected for cuticular hydrocarbon analysis.

Morphological identifications were done using the identification key of adult Calliphoridae (Lutz et al. 2018) at Entomology Department, Faculty of Science, Ein Shams University, Egypt.

Cuticular hydrocarbon analysis

The extraction procedures of Ye et al. (2007) were slightly modified. Three replicates were analyzed for each species. Eight puparial cases of L. sericata and Ch. albiceps and only six puparial cases of Ch. marginalis (as large size) were used for each replicate. Puparial exuviae were washed with distilled water, cleaned by tip of fine paint brush and then dried at filter papers. Puparial exuviae of each replicate were immersed in 5mL n-hexane in glass vial and gentle swirl for 10 to 15 min at room temperature. After that, puparial cases were removed from the extracts. The extracts were then filtrated and collected in clean glass vials and stored at -20 ˚C till GC-MS analysis.

Gas chromatography–mass spectrometry analysis (GC-MS)

The GC-MS system (Agilent Technologies) was equipped with gas chromatograph (7890B) and mass spectrometer detector (5977 A) at Central Laboratories Network, National Research Centre, Cairo, Egypt. The GC was equipped with HP-5MS column (30 m x 0.25 mm internal diameter and 0.25 μm film thickness). Analyses were carried out using helium as the carrier gas at a flow rate of 1.0 ml/min at a splitless, injection volume of 1 µl and the following temperature program: 45 °C for 2 min; rising at 10 °C /min to 300 °C and held for 10 min. The injector and detector were held at 280 and 300 °C, respectively. Mass spectra were obtained by electron ionization (EI) at 70 eV; using a spectral range of m/z 25–700. Identification of different constituents was determined by comparing the spectrum fragmentation pattern with those stored in Wiley and NIST Mass Spectral Library data.

Statistical analysis

In the present study, the data were analyzed using SPSS version 22. The peaks were selected for analysis based on their retention time (RT). Peaks with RT > 18 was chosen and annotated according to their retention times. The relative abundance of each peak was calculated based on their computed area under curve for each hydrocarbon. Discriminant analysis was executed, and discriminant functions were conducted using Wilk’s lambda method. One-way ANOVA was applied to study the statistical effect of fly species on the percent composition of hydrocarbons. Least significant differences (LSD) test was used to illustrate the statistical differences in the studied variables among the different species.

Results

Twelve classes of chemical compounds were identified from the empty puparial cases of the three fly species at retention time 18.78 to 35.03. The percentages and the types of the extracted compounds were listed in Table 1. The profiles for Chrysomya albiceps and Lucilia sericata were very similar. Alcohol represented the highest percentage of compounds with 28.6% in L. sericata. However, in the three fly species L. sericata, Ch. albiceps and Chrysomya marginalis, alkanes (n-alkanes, branched alkanes and cycloalkanes) constitute the major component of cuticular hydrocarbons with 28.5, 50 and 89.4%, respectively. The chromatographs in Fig. 1 showed the characteristic peaks for each fly. Among studied species, the CHCs abundance in Ch. albiceps was lower than that in the others (18 compounds) as shown in Table 2.

Fig. 1
figure 1

Representative gas chromatographs of cuticular compounds of Lucilia sericata, Chrysomya albiceps and Chrysomya marginalis

The retention times, names, and the frequency of each hydrocarbon in the three flies were listed (Table 2). Forty-two compounds were identified with chain lengths ranging from C12 to C45. Heptacosane (7.6%) is the n-alkane dominated the chemical profile of Ch. marginalis, while Dodecane is the major one found in L. sericata (2.36%) and Ch. albiceps (1.88%). The predominant methyl branched alkane is 2-Methyltetracosane in L. sericata and Ch. albiceps and accounted for 5.53 and 4.92%, respectively. While the most abundant methyl branched alkane in Ch. marginalis is 2,6,10,14-Tetramethylhexadecane, (2.65%). The three species shared one compound in common which is the cycloalkane, 1-(2-Octyldecyl)octahydropentalene. Halogen branched hydrocarbons were detected in the chemical profiles of L. sericata and Ch. albiceps, but none was found in Ch. marginalis profile. Also, alkenes with different function groups as acid anhydride, alcohol and ester were detected in L. sericata and/or Ch. albiceps as illustrated in Table 2 peaks number 29, (26, 18) and (11, 7), respectively.

Alkadienes were represented in L. sericata with peak 30 (aldehyde) and in L. sericata and Ch. albiceps profiles with peak 35 (alcohol). The only cycloalkadienes observed in the chromatogram is the ketone compound, 3-(Dodecenyl)dihydro-2,5-furandione (peak 36) in the profiles of L. sericata and Ch. marginalis. The later species revealed several specific compounds demonstrated by peaks 34, 33, 31, 28, 27, 25, 23, 20, 17, 15, 13, 9, 8, 6, 5, 4 and 2. The chromatogram of Ch. marginalis shows more alkanes than those found in other species. According to test equality of group means, most peaks differed significantly (P < 0.01) among all species.

Discriminant analysis

According to multiple regression analyses using Fisher discriminant method, twelve spectral peaks identified as characteristic variables among the three species of flies. They include peak 1 (Hydroxymethylcyclododecane), peak 2 (Heneicosane), peak 3 (Dodecane), peak 4 (Triacontane), peak 5 (Tetracosamethyl-cyclododecasiloxane), peak 6 (Tricosane), peak 7 (Oxalic acid, allyl pentadecyl ester), peak 10 (9-t-Butyl-4-iodo-2,2-dimethyladamantane), peak 11 (Oxalic acid, allyl octadecyl ester), peak 12 (2-Ethyl-1-decanol), peak 14 (2-Butyl-1-octanol), and peak 18 (Phytol).

Peaks 1, 3, 11, 12, 14 and 18 were found in Lucilia sericata. Peaks 1, 3, 7, 10, 14, and 18 were identified in Chrysomya albiceps. However, in Chrysomya marginalis, peaks 2, 4, 5 and 6 were only presented. According to one-way ANOVA, all the peaks showed significant differences among the studied species.

Two canonical standardized functions were obtained by discriminant analysis (Table 3). Function 1 explained 83.8% of the variations in the dependent variables (Fly species) and function 2 interpreted 16.2% of variable rate. To validate this result, species were plotted according to their scores on these two functions. Each individual was correctly assigned to its species as observed in Fig. 2 permitting establishment of a confident identification.

Fig. 2
figure 2

Cuticular hydrocarbon composition of Lucilia sericata, Chrysomya albiceps and Chrysomya marginalis distributed in the space of discriminant functions 1 and 2

In Table 4, the unstandardized coefficients of the canonical discriminant function are displayed. The higher the value of the coefficient, the higher the ability to predict the change in the dependent variable. The canonical discriminant function for the three species is discernible.

Table 1 The frequency and percentage of each hydrocarbon classes in the three flies
Table 2 Classes and percent composition of the compounds isolated from the cuticle of the three flies
Table 3 Summary of Canonical Discriminant Functions
Table 4 Canonical discriminant function coefficients (unstandardized)

Discussion

Cuticular hydrocarbons are proven to be species-specific in many insect taxa including Diptera (Braga et al. 2013; Moore et al. 2022). It is expected to be a promising tool when comes into the field of forensic entomology especially in cases where only empty puparia are available in a scene. Analysis of CHCs provides very helpful information in identifying ambiguous specimens due to physical damage, degradation of the genetic material or even in case of sexually dimorphic or morphologically similar species (Braga et al. 2013; Moore et al. 2021). Moreover, various studies confirmed the reliability of this technique in assigning individuals to certain geographic population (Charabidze et al. 2017; Moore et al. 2022) which in turn can reveal the presence of non- native population on a cadaver, hence cadaver movement from original death location. The main aim of this study was to establish if a distinction could be made among the empty puparial cases of the three blow fly species (Lucilia. sericata, Chrysomya albiceps and Chrysomya marginalis) using cuticular hydrocarbon analysis. As far as we know, this is the first study that deals with the cuticular chemical composition of some Egyptian flies of forensic importance. More investigation should be done on the cuticle of necrophagous flies as it could greatly facilitate species identification and accelerate solving forensic cases without the need to rear larvae or pupae to adult stage (Paula et al. 2017).

Morphological differentiation can be noticed among the adults of those flies (Lutz et al. 2018), while identification of larvae is time consuming and challenging specially in early instars (Szpila et al. 2014). When comes into pupae, usual morphological distinction is very difficult or even impossible for scientists other than taxonomists due to deformation or weathering conditions (Ye et al. 2007; Moore et al. 2022). Despite being known in many insect species, the chemical composition of the cuticle of many Egyptian species is still unknown and requires a thorough investigation. Our results showed that the three fly species have a distinct fingerprint profile. Their CHCs are like those of other insects and consisted of alkanes, methylalkanes, halogenated alkanes and cyclic hydrocarbons (Ye et al. 2007; Braga et al. 2013; Galhoum 2018; Moore et al. 2022). We also included all compounds obtained from the chromatogram like alcohols, ketones, aldehydes, esters and acids into our analysis. As detected previously (Frederickx et al. 2012; Kranz et al. 2017); those compounds yielded distinct peaks that can be used to distinguish between the three species. The classes of the chemical compounds obtained from the chromatogram of L. sericata and Ch. albiceps included hydrocarbons and alcohols, ketones, esters, ethers, acid anhydrides, epoxides and an aldehyde. Similar results were recorded by many authors as (Al-Dawsary 2014) who found that the most prevalent chemical groups in the cuticle of the red palm weevil Rhynchophorus ferrugineus (Olivier) are alcohols then hydrocarbons, carboxylic acid, esters, aldehydes and ketones respectively. While, (Alnajim et al. 2019) found the most abundant classes in the cuticle of Tribolium castaneum (Herbst) and Rhyzopertha dominica (Fabricius) are hydrocarbons, fatty acids and a sterole. Same was obtained by (Elshaier 2021) who found the cuticle of Anthidium amabile (Alfken) dominated by fatty acids then hydrocarbons sterols, glycerides, one ketone and one alcohol. Although the number of hydrocarbons in the cuticle of L. sericata (6 compounds) and Ch. albiceps (7 compounds) recorded in this work was significantly smaller than those obtained from previous studies on the same species (Ye et al. 2007; Braga et al. 2013; Moore et al. 2014, 2022; Paula et al. 2017), some authors like (Elshaier 2021) recorded only five hydrocarbons from the cuticle of the wool-carder bees Anthidium amabile from Egypt. However, (Drijfhout et al. 2009) estimated the total number of hydrocarbons in the cuticle of insects as ranging from five to fifty compounds. Similar to many other dipteran flies, the chemical profile of the cuticle of Ch. marginalis composed mainly of n-alkanes, with the most abundant compound is heptacosane (C27: H56) (Goodrich 1970; Ye et al. 2007; Moore et al. 2022; Kula et al., 2023). Squalene was found only in the profile of Ch. marginalis and most likely was ingested during feeding as insects don’t produce this compound (Braga et al. 2013).

Our results showed that, the only shared compound between the three flies is 1-(2-Octyldecyl)octahydropentalene (C26: H50). This compound was encountered in essential oils extracted from medicinal plants for cytotoxic, antimicrobial and insecticidal activities (Mohamed et al. 2015; Al-Mazroa et al. 2015; Hamada et al. 2018; Sadiq et al. 2018; Mamza et al. 2021; Kewlani et al. 2022). Also, was detected in ground water samples used for drinking and irrigation in Egypt (Abd-Elgawad et al. 2022). So, the presence of this substance may be due to the feeding habits of the three fly species. Kranz et al. (2017), found that diet outmost the impact of any other abiotic factors on the structure of insects cuticle, resulting in significant influence on their profiles. Until now, there is no study reported the presence of such compound in insect cuticle and the exact role of it is still unknown.

In conclusion, the use of GC-MS chemical analysis of puparial cases can accurately distinguish between the studied blow fly species without the need for specialized taxonomists for identification. This method has a lot of potential to be exploited in criminal investigations and post mortem interval estimation. Further research is needed to confirm these findings and to investigate the impact of factors such as temperature, diet, and location on cuticular components.