Intrapuparial Development and Age Estimation of Calliphora grahami (Diptera: Calliphoridae) for Postmortem Interval Estimation

Development and age estimation of the intrapuparial period of Phormia regina (Diptera: Calliphoridae) for postmortem interval estimation

Phormia regina (Meigen, 1826; Diptera: Calliphoridae) is a Holarctic species that rapidly colonizes carcasses and has been used as an indicator for determining the minimum postmortem interval. However, studies using morphological methods to estimate the intrapuparial age of P. regina are lacking. In this study, morphological changes within the puparium were observed under a stereomicroscope at 7 constant temperatures ranging from 16 °C to 34 °C. The intrapuparial period was categorized into 12 substages. Morphological indicators, including compound eyes, mouthparts, antennae, thorax, legs, wings, and abdomen, were recorded in detail. The observed morphological changes were divided into 6-10 substages, and the duration of each substage was also recorded in detail. The results of this study provide primary data for using the intrapuparial morphology of P. regina when pupae are collected at a crime scene and estimating the minimum postmortem interval.

Research trends on forensic entomology for five decades worldwide

Forensic entomology, the study of insects and arthropods in criminal investigations, is crucial for estimating postmortem intervals, determining changes in corpse positioning, and identifying causes of death. This study analyzes global research trends in forensic entomology using data from the Scopus database spanning 1970 to 2024, with data visualized through VOSviewer. A total of 2,261 articles were identified, with an average productivity of 42 papers per year. The leading countries in forensic entomology research are the United States (n = 444), Brazil (n = 266), China (n = 198), the United Kingdom (n = 194), and Germany (n = 156). Current hot topics in the field include pupae, feeding behavior, beetles, and genetics. These findings underscore the ongoing interest and advancements in forensic entomology, highlighting its significance and diverse applications in criminal investigations.

Intrapuparial stage aging and PMI estimation based on the developmental transcriptomes of forensically important Aldrichina grahami (Diptera: Calliphoridae) gene expression

Background

The expression profiles of differentially expressed genes (DEGs) during pupal development have been demonstrated to be vital in age estimation of forensic entomological study. Here, using forensically important Aldrichina grahami (Diptera: Calliphoridae), we aimed to explore the potential of intrapuparial stage aging and postmortem interval (PMI) estimation based on characterization of successive developmental transcriptomes and gene expression patterns.

Methods

We collected A. grahami pupae at 11 successive intrapuparial stages at 20 °C and used the RNA-seq technique to build the transcriptome profiles of their intrapuparial stages. The DEGs were identified during the different intrapuparial stages using comparative transcriptome analysis. The selected marker DEGs were classified and clustered for intrapuparial stage aging and PMI estimation and then further verified for transcriptome data validation. Ultimately, we categorized the overall gene expression levels as the dependent variable and the age of intrapuparial A. grahami as the independent variable to conduct nonlinear regression analysis.

Results

We redefined the intrapuparial stages of A. grahami into five key successive substages (I, II, III, IV, and V), based on the overall gene expression patterns during pupal development. We screened 99 specific time-dependent expressed genes (stage-specific DEGs) to determine the different intrapuparial stages based on comparison of the gene expression levels during the 11 developmental intrapuparial stages of A. grahami. We observed that 55 DEGs showed persistent upregulation during the development of intrapuparial A. grahami. We then selected four DEGs (act79b, act88f, up and ninac) which presented consistent upregulation using RT-qPCR (quantitative real-time PCR) analysis, along with consideration of the maximum fold changes during the pupal development. We conducted nonlinear regression analysis to simulate the calculations of the relationships between the expression levels of the four selected DEGs and the developmental time of intrapuparial A. grahami and constructed fitting curves. The curves demonstrated that act79b and ninac showed continuous relatively increasing levels.

Conclusions

This study redefined the intrapuparial stages of A. grahami based on expression profiles of developmental transcriptomes for the first time. The stage-specific DEGs and those with consistent tendencies of expression were found to have potential in age estimation of intrapuparial A. grahami and could be supplementary to a more accurate prediction of PMI.

Age determination of Chrysomya megacephala (Diptera: Calliphoridae) using lifespan patterns, gene expression, and pteridine concentration under constant and variable temperatures

Chrysomya megacephala (Fabricius, 1794) (Diptera: Calliphoridae), is a blowfly species widely studied in medical, veterinary, and entomological research. Our study examined the impact of constant (15, 20, 25, 30, and 35 °C) and variable (ranging from 21.0 to 25.4 °C, with an average of 23.31 °C) temperatures on the development and larval body length of C. megacephala. Additionally, we analyzed the age of the adult C. megacephala through pteridine content and related metabolic genes analysis. Our findings revealed three distinct growth patterns: isomorphen diagram, isomegalen diagram, and thermal accumulated models. At constant temperatures of 15, 20, 25, 30, and 35 °C, egg-hatching times were 44.5 ± 8.9, 26.7 ± 4.6, 12.6 ± 1.1, 11.0 ± 1.0, and 9.9 ± 1.9 h, respectively, while it was 15.3 ± 5.9 h at variable temperatures. The total development times from oviposition to adult eclosion in C. megacephala required 858.1 ± 69.2, 362.3 ± 5.9, 289.6 ± 17.8, 207.3 ± 9.3, and 184.7 ± 12.1 h at constant temperatures of 15, 20, 25, 30, and 35 °C, respectively. This duration was extended to 282.0 ± 64.1 h under variable temperatures. However, no significant differences were found in hatching times and the total developmental durations between 25 °C and variable temperatures. A developmental threshold temperature (D0) of 9.90 ± 0.77 °C and a thermal summation constant (K) of 4244.0 ± 347.0° hours were ascertained. Pteridine content patterns varied significantly across constant temperatures, but not between 25 °C and variable temperatures. Sex and temperature were identified as the primary factors influencing pteridine levels in the head of C. megacephala. Gene expression associated with pteridine metabolism decreased following adult eclosion, matching with increased pteridine concentration. Further investigations are needed to explore the use of pteridine cofactors for age-grading adult necrophagous flies. These findings provide valuable insights into the lifespan of C. megacephala, thereby offering valuable groundwork for forthcoming investigations and PMImin determination.

Estimating the intra-puparial period of Hydrotaea spinigera (Stein,1910) (Diptera: Muscidae) with morphological and gene expression changes

Hydrotaea spinigera (Stein, 1910) (Diptera: Muscidae) is a forensically important sarcosaprophagous species widely distributed throughout the Oriental and Australasian regions. At the advanced decomposition stage or the skeletonize stage, the immature stages of H. spinigera, especially the pupae, can still be found in large quantities and could be used as important indicators to estimate the minimum postmortem interval (PMI_min). However, there have been no studies on the intra-puparial period of this species. Herein, we studied morphological and differential gene expression changes during the intra-puparial development of H. spinigera, aiming to accurately estimate the intra-puparial age of H. spinigera. The intra-puparial morphological changes of H. spinigera were observed at seven constant temperatures ranging from 16°C to 34°C and divided into 12 sub-stages. Structures that could be used to estimate the intra-puparial age, such as compound eyes, mouthparts, antennae, thorax, legs, wings, and abdomen, were observed in detail, and the developmental process of each structure was divided into 5 to 10 stages. The time range of each sub-stage, or when a structure appeared, was recorded. For the gene expression section, the most suitable reference genes were screened by geNorm, NormFinder, BestKeeper, and ΔCt methods. Based on the selected reference genes, real-time quantitative PCR (RT qPCR) was used to detect the expression changes of the white and hsp90 genes with developmental time at 19°C, 25°C, and 31°C. Results showed that the trend of Hsp90 gene expression under different temperatures was not consistent, while white genes exhibited regular changes during development, and could thus be used for age estimation of H. spinigera. This study provides an important basis for forensic entomology to use morphological and differential gene expression for estimating the age of H. spinigera during the intra-puparial period. Moreover, the combination of the two methods can produce a more accurate minimum postmortem interval (PMI_min) estimate compared to when each method is used separately.

Multimethod combination for age estimation of Sarcophaga peregrina (Diptera: Sarcophagidae) with implications for estimation of the postmortem interval

Sarcophaga peregrina (Robineau-Desvoidy, 1830) (Diptera: Sarcophagidae) is a forensically important flesh fly with potential value for estimating the minimum postmortem interval (PMI_min). Here, the developmental patterns of S. peregrina were investigated at 5 constant temperatures (15-35 °C). Morphological changes at different developmental stages of the pupa were observed at 4 constant temperatures (15-30 °C) by removing the puparium and staining the pupa with hematoxylin and eosin. Furthermore, differentially expressed genes (DEGs) were analyzed at 25 °C in the intrapuparial period to estimate the age of S. peregrina during the intrapuparial stage. S. peregrina completed development from larviposition to adult eclosion at 15 °C, 20 °C, 25 °C, and 30 °C; the developmental durations were 1090.3 ± 30.6 h, 566.6 ± 21.9 h, 404.6 ± 13.01 h, and 280.3 ± 4.5 h, respectively, while the development could not be completed at 35 °C. The intrapuparial period of S. peregrina was divided into 12 sub-stages on the basis of the overall external morphological changes; 6 sub-stages on the basis of individual morphological structures such as the compound eyes, antennae, thorax, legs, wings, and abdomen; and 10 sub-stages on the basis of internal morphological changes detected using histological analysis. The period of each sub-stage or structure that appeared was determined. Moreover, we found that 6 genes (NDUFS2, CPAMD8, NDUFV2, Hsp27, Hsp23, and TPP) with differential expression can be used for the precise age estimation of S. peregrina during the intrapuparial period. This study provided basic developmental data for the use of S. peregrina in PMI_min estimation, and we successfully estimated PMI_min in a real forensic case by using a multimethod combination.

Pupal Age Estimation of Sarcophaga peregrina (Diptera: Sarcophagidae) at Different Constant Temperatures Utilizing ATR-FTIR Spectroscopy and Cuticular Hydrocarbons

Sarcophaga peregrina (Robineau-Desvoidy, 1830) (Diptera: Sarcophagidae) is a forensically important flesh fly that has potential value in estimating the PMImin. The precise pupal age estimation has great implications for PMImin estimation. During larval development, the age determination is straightforward by the morphological changes and variation of length and weight, however, the pupal age estimation is more difficult due to anatomical and morphological changes not being visible. Thus, it is necessary to find new techniques and methods that can be implemented by standard experiments for accurate pupal age estimation. In this study, we first investigated the potential of attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectroscopy and cuticular hydrocarbons (CHCs) for the age estimations of S. peregrina pupae at different constant temperatures (20 °C, 25 °C, and 30 °C). The orthogonal projections latent structure discrimination analysis (OPLS-DA) classification model was used to distinguish the pupae samples of different developmental ages. Then, a multivariate statistical regression model, partial least squares (PLS), was established with the spectroscopic and hydrocarbon data for pupal age estimations. We identified 37 CHCs with a carbon chain length between 11 and 35 in the pupae of S. peregrina. The results of the OPLS-DA model show a significant separation between different developmental ages of pupae (R2X > 0.928, R2Y > 0.899, Q2 > 0.863). The PLS model had a satisfactory prediction with a good fit between the actual and predicted ages of the pupae (R2 > 0.927, RMSECV < 1.268). The results demonstrate that the variation tendencies of spectroscopy and hydrocarbons were time-dependent, and ATR-FTIR and CHCs may be optimal for the age estimations of pupae of forensically important flies with implications for PMImin estimation in forensic practice.