Movatterモバイル変換

[0]ホーム
URL:

Skip to main content
NCBI home page
Search in PMCSearch
As a library, NLM provides access to scientific literature. Inclusion in an NLM database does not imply endorsement of, or agreement with, the contents by NLM or the National Institutes of Health.
Learn more:PMC Disclaimer | PMC Copyright Notice
Scientific Reports logo

Temperature models of development forNecrodes littoralis L. (Coleoptera: Silphidae), a carrion beetle of forensic importance in the Palearctic region

Joanna Gruszka1,2,3,✉,Szymon Matuszewski1,2
1Laboratory of Criminalistics, Adam Mickiewicz University, Św. Marcin 90, 61-809 Poznań, Poland
2Center for Advanced Technologies, Adam Mickiewicz University, Uniwersytetu Poznańskiego 10, 61-614 Poznań, Poland
3Department of Animal Taxonomy and Ecology, Adam Mickiewicz University, Uniwersytetu Poznańskiego 6, 61-614 Poznań, Poland
✉

Corresponding author.

Received 2022 Feb 25; Accepted 2022 May 30; Collection date 2022.

© The Author(s) 2022

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visithttp://creativecommons.org/licenses/by/4.0/.

PMCID: PMC9188545  PMID:35690667

Abstract

Analysis of insects can provide evidence in death cases, for example, by answering the question about the time of death. Apart from flies, beetles are the second most useful insect group in forensic entomology. To elucidate the time of death based on insect evidence, developmental models of a given species are necessary. In this study, we developed such models forNecrodes littoralis, a necrophagous beetle, which is common in the Palearctic region and has great potential in forensic entomology. We monitored the development at 10 constant temperatures (14–30 °C). Larvae were reared in aggregations. Thermal summation models, isomorphen and isomegalen diagrams and growth curves were derived using the data. Depending on the temperature, development lasted between about 23 and 89 days. Mortality was high at the extremes of the temperature range. The thermal summation constant for the total development was 434.7 ± 28.86 accumulated degree-days above a developmental threshold of 9.04 ± 0.55 °C. This is the first comprehensive dataset on the development ofN. littoralis. Implications for its use in forensic casework are discussed.

Subject terms: Entomology, Ecology

Introduction

Forensic entomology uses insects and other arthropods as evidence in legal investigations. Arthropods are used mainly in homicide, suicide and mysterious death cases, when insect evidence allows to establish some circumstances of death1,2. For instance, due to the ecological specialization of specific insect taxa, it is possible to determine the environment where death occurred and to answer the question whether the body was relocated after death3,4. Most often, however, insects are used to elucidate the time of death, usually through the estimation of minimum postmortem interval (PMImin)5. Two main methods can be used for this purpose. The first one uses the regularity of insect succession on carcasses6,7. Using models of insect succession on carrion for given environmental conditions and a set of insect taxa found on the death scene, we can try to estimate how much time has passed since death8. The second method uses the regularity of insect development, in particular a close relationship between temperature and the rate of development1,9. Developmental studies of necrophagous insects allow for the creation of developmental models that can be used to estimate PMImin.

Forensic entomologists created various kinds of developmental models. Some of them, such as the isomegalen diagram or growth curve, represent the relationship between development time and size of the larvae at a given temperature. The isomorphen diagram shows the duration of developmental stages at given temperatures. The thermal summation model (TSM) assumes that within the certain temperature range, there is a linear relationship between the rate of development and temperature and that development stops below a certain temperature10. Accordingly, there is some constant amount of heat that needs to be accumulated by insects of a given species to reach certain developmental landmark (e.g. hatching, pupation or eclosion). This value is called a thermal summation constant (K) and is expressed in accumulated degree hours (ADH) or accumulated degree days (ADD). The temperature below which the insect species ceases to develop is termed a lower developmental threshold (D0) and is expressed in temperature units10. One of the most frequently used methods to derive TSM, is the one proposed by Ikemoto and Takai11. It is a modification of the classic thermal summation method. Authors proposed an equation:

DT=K+D0∗D

whereD is the duration of development,D0 is the lower developmental threshold,K is the thermal summation constant andT is the environmental temperature. The authors also suggested using the Reduced Major Axis (RMA) regression to derive TSM. Slope of the RMA model is the lower developmental threshold (D0) and they intercept is the thermal summation constant (K). Although TSMs simplify complexity of insect development and ignore its substantial intraspecific variation (e.g. variation in development time between the sexes12,13 or between insects of different sizes14,15), their use enables satisfactorily accurate estimation of insect age and eventually PMImin8,16,17.

Flies are the most frequently used group of insects in forensic entomology. They usually appear on cadavers as first colonizers and age of their immature stages allow to estimate the PMImin that is close to the true PMI18. Typically, the first colonizers are blow flies (Calliphoridae). They appear on the corpse within the first hours or even minutes after death10,19. During feeding, their larvae form large aggregations20–22. Beetles appear later on carcasses, and are often present there until the remains stage. Hence, they can be an important tool in estimating the PMI in the advanced stages of cadaver decomposition23. Numerous studies of insect succession on carrion and descriptions of specific forensic cases demonstrated that beetles may be useful for the estimation of PMI6,17,24–34. Among the forensically important families of beetles, the most frequently mentioned ones are carrion beetles (Silphidae), rove beetles (Staphylinidae), checkered beetles (Cleridae), skin beetles (Dermestidae), clown beetles (Histeridae) and sap beetles (Nitidulidae). The greatest limitation for the use of beetles in forensic entomology is the lack of developmental models.

Developmental models for only 15 species of forensically important beetles have been published so far. First models were created forThanatophilus micans (Silphidae), population from South Africa35. The other carrion beetles with published models comprise:Oxelytrum discicolle population from northern part of South America36,Thanatophilus mutilatus population from South Africa37,Necrodes littoralis15 population from Central Europe,Necrophila (Calosilpha) brunnicollis population from East Asia38 and populations from Central Europe ofThanatophilus sinuatus39 andThanatophilus rugosus40. Models were created also for four species ofDermestidae: populations from south-western Europe ofDermestes frischi, D. undulatus andD. maculatus and Chinese population ofD. tessellatocollis41,42. Developmental models were also published for central European population ofSciodrepoides watsoni (Leiodidae: Cholevinae)43, Chinese44 and central European45 populations ofCreophilus maxillosus (Staphylinidae) and Chinese populations ofNecrobia rufipes (Cleridae)46 andOmosita colon (Nitidulidae)47. Although there are beetle taxa that were extensively studied (e.g.Thanatophilus orCreophilus), still many forensically useful beetle species lack developmental models.

Necrodes littoralis (Silphidae) is widely distributed in the Palearctic region. It prefers open and forest natural habitats, but has also been recorded in urban open and quasi-indoor habitats24,31,34,48,49. Larvae feed on carrion mainly in spring and summer, adult beetles have also been recorded in the fall24,48. The species is associated with cadavers at active and advanced decay7,31,48,50.N. littoralis prefers large carrion, on which its larvae frequently form aggregations that allow them to drive active decay similarly to blow flies50–52.N. littoralis has been reported from many forensic cases, however it has infrequently been used to estimate PMI17,31,34,49,53–55. Comprehensive data on its occurrence on human cadavers in France was provided by Charabidze et al.31. According to their analysis,N. littoralis (larvae or adult beetles) were present in 154 cases (1 in 8 cases examined), with 91.6% of outdoor cases (mainly forests and bushes). Most of them occurred during spring or summer. In more than 85% of the casesN. littoralis was observed from the early to advanced decomposition stages. Moreover, the pre-appearance interval (PAI) of adult and larvalN. littoralis was found to be strongly related to the preceding ambient temperature56 and for this reason, it may easily be estimated using temperature methods for PAI57. Furthermore, the size of adultN. littoralis was found to be negatively correlated with its physiological age at maturity. Thus, such features as length or weight of adultN. littoralis can be useful to calibrate developmental constants for this species and to improve the accuracy of age estimation15. Some partial developmental datasets forN. littoralis were already published15,58. However, there is no comprehensive dataset providing all development models and based on the full temperature range for this species. This deficiency partially explains low frequency ofN. littoralis use for the estimation of PMI. Current study aims to create the first comprehensive set of robust development models forN. littoralis. Some of the data used in this work was also used for previously published analyzes15.

Results

Life cycle

AdultNecrodes littoralis is mainly active after dark. Mating occurs usually during the night. The female lays eggs in the soil, in batches, usually 50–70 eggs in each. The first instar larvae hatch creamy white and migrate quickly in search for food. It is the period when they are very fragile and particularly vulnerable to injuries. They darken, and their cuticle hardens with time. Likewise, the second and third instar larvae are white and non-sclerotized shortly after ecdysis. When feeding is complete, third instar larvae burrow into the ground. Then they form pupal chambers by thrashing the abdomen and thus compacting the soil around them. They go through the prepupal, pupal and teneral adult stages inside the chambers. The pupa is creamy white at the beginning and with time it gradually sclerotizes and darkens. The metamorphosis ends with the appearance of tenerals—creamy white and non-sclerotized adult beetles. Over time, their cuticle hardens and darkens, giving the insects their final black color. The beetles dig out of the pupal chamber after they became fully sclerotized and colored.

Influence of in vivo measurements

In general, non-measured beetles developed longer (Fig. 1), with significant differences recorded at 15 °C (Mann–WhitneyU test:Z = 6.02,p < 0.001), 16 °C (Mann–WhitneyU test:Z = 5.92,p < 0.001), 19 °C (Mann–WhitneyU test:Z = 3.56,p < 0.001) and 20 °C (Mann–WhitneyU test:Z = 3.20,p = 0.001). The development time of non-measured beetles was shorter only at 14 °C (Mann–WhitneyU test:Z = − 3.16,p = 0.001). The largest difference was observed in 16 °C, where non-measured beetles developed 5.39% longer than measured beetles. However, in higher temperatures, starting from 17 °C, the differences between measured and non-measured beetles were very small (Fig. 1).

Figure 1.

Figure 1

Differences in duration of total immature development between measured and non-measured specimens ofN. littoralis at ten constant temperatures. *Statistically significant difference in Mann–WhitneyU test at α = 0.001.

Developmental models calculated using measured and non-measured specimens were only slightly different (Table1). There were, however no significant differences in the relative error of age estimation between the tested models (Wilcoxon signed-rank test:Z = 0.96,p = 0.34,N = 110). Both models yielded estimates of age with the average relative error below 0.08 (Fig. 2).

Table 1.

Thermal summation models for the total immature development ofN. littoralis, calculated using measured or non-measured specimens.

ModelTemperature range [°C]Thermal summation constant—K (SE) [degree-days]Developmental threshold—D0 (SE) [°C]r2Np
Measured14–30421.09 (29.06)9.251 (0.552)0.9689 < 0.001
Non-measured14–30436.90 (29.05)9.051 (0.548)0.9679 < 0.001

Figure 2.

Figure 2

Relative error of age estimation using models calculated with measured or non-measured beetles.

Mortality and development

The highest mortality was observed at extreme temperatures: 86.5% at 14 °C and 64.75% at 30 °C. At the other temperatures mortality was below 25%, except for 17 °C where it was 53% (Fig. 3). There were differences between the temperatures in the distribution of deaths across the life stages (Fig. 3.) At the lowest and the highest temperatures deaths occured during all developmental stages, while at the optimal temperatures beetles died mostly at postfeeding larval or pupal stages. At 17 °C deaths were also distributed across all developmental stages. The hypothesis as to why and consequences will be presented in the discussion section.

Figure 3.

Figure 3

Mortality ofN. littoralis at immature developmental stages in different rearing temperatures.

The development from egg to adult stage took between 22.84 days at 30 °C and 89.11 days at 14 °C (Table2, Fig. 4.). Results from 17 °C deviated from the general pattern of decrease in development time with an increase in temperature. The first and second instar larvae at 17 °C developed much longer than larvae at 16 °C and 15 °C, and thus the total development time at 17 °C lengthened on average to 60.98 days. Due to this inconsistency and the exceptionally high mortality of larvae at 17 °C, we decided not to use the data from this temperature when developing the models (except for the model for hatching).

Table. 2.

Duration of immature developmental stages ofN. littoralis at ten constant temperatures.

Temp. [°C]Mean duration of development [days] (SE;N)
Egg1st instar larva2nd instar larva3rd instar larvaPupaTotal development
149.63 (0.15; 8)10.33 (0.23; 8)10.31 (0.41; 8)41.35 (0.07; 63)18.19 (0.42; 54)89.11 (0.45; 54)
156.55 (0.22; 8)7.65 (0.17; 8)7.15 (0.35 ;8)33.89 (0.15; 319)16.34 (0.12; 302)71.29 (0.17; 302)
165.95 (0.12; 8)5.27 (0.11; 8)4.83 (0.11; 8)28.97 (0.16; 373)13.35 (0.08; 370)57.78 (0.14; 370)
175.53 (0.08; 8)9.54 (1.31; 8)8.55 (0.52; 8)26.90 (0.16; 240)12.10 (0.10; 188)60.98 (0.21; 188)
184.93 (0.15; 8)4.13 (0.16; 8)4.01 (0.16; 8)23.97 (0.11; 354)11.85 (0.10; 341)48.81 (0.20; 341)
194.02 (0.12; 8)4.04 (0.18; 8)3.92 (0.21; 8)22.61 (0.13; 380)11.80 (0.06; 364)45.14 (0.13; 364)
203.67 (0.06; 8)3.96 (0.10; 8)3.77 (0.15; 8)22.27 (0.11; 400)10.35 (0.07; 357)43.75 (0.08; 357)
222.73 (0.09; 8)2.35 (0.08; 8)2.29 (0.08; 8)13.46 (0.05; 390)7.43 (0.02; 385)28.28 (0.06; 385)
262.66 (0.04; 8)1.55 (0.03; 8)1.68 (0.08; 8)12.82 (0.07; 338)7.03 (0.06; 327)25.65 (0.08; 327)
302.38 (0.09; 8)1.39 (0.07; 8)1.45 (0.06; 8)11.89 (0.05; 187)5.93 (0.07; 141)22.84 (0.10; 141)

Figure 4.

Figure 4

Isomorphen diagram forN. littoralis based on median times to reach developmental events at each of the nine temperatures tested (Supplementary Table1). Horizontal bars represent interquartile ranges. Areas between lines represent developmental stages. Different symbols represent developmental events.

Immediately after hatching, the larvae had an average length of 6.91 ± 0.79 mm (Isomegalen diagram: Supplementary Fig. 1, Growth curves: Supplementary Figs. 2–10). At the growth peak, larvae were on average 25.14 ± 1.82 mm in length. Moulting usually occurred at a similar length of the larvae: 10.71 ± 0.79 mm for the first ecdysis and 16.10 ± 1.30 mm for the second ecdysis. Growth curves were distinctly sigmoidal (Supplementary Figs. 2–10).

All included temperature points were within 95% confidence interval for all reduced major axis regression models (Fig. 5). The coefficient of determination for each thermal summation model exceeded 0.96, indicating high fit of the models. Lower developmental thresholds (D0) ranged from 9.04 °C for eclosion to 10.86 °C for second ecdysis. Thermal summation constants (K) were between 39.97 ADD for hatching and 434.705 ADD for eclosion (Table3).

Figure 5.

Figure 5

Reduced major axis (RMA) regression models for five developmental events ofN. littoralis. DT is the time in days to reach developmental event multiplied by the constant rearing temperature. Temperature values are noted next to the points. Dashed lines represent 95% confidence intervals.

Table 3.

Thermal summation models for five developmental events ofN. littoralis calculated using Ikemoto and Takai method.

Developmental eventTemperature range [°C]Model equationThermal summation constant—K (SE) [degree days]Developmental threshold—D0 (SE) [°C]Nr2p
Hatching14–30DT = 39.977 + 9.627 × D39.977 (2.902)9.627 (0.539)100.968 < 0.001
First ecdysis14–30DT = 65.747 + 10.583 × D65.747 (3.235)10.583 (0.315)90.992 < 0.001
Second ecdysis14–30DT = 92.060 + 10.861 × D92.060 (5.096)10.861 (0.326)90.992 < 0.001
Pupation14–30DT = 303.294 + 9.701 × D303.294 (21.680)9.701 (0.531)90.973 < 0.001
Eclosion14–30DT = 434.705 + 9.044 × D434.705 (28.862)9.044 (0.548)90.967 < 0.001

Discussion

Current results on the development time are consistent with the development times ofN. littoralis obtained by Dekeirsschieter58. She studied the development at two temperatures only. The development (without the egg stage) took on average 42.79 days at 18 °C and 23.60 days at 23 °C. By subtracting the egg stage from the mean development time at equivalent temperature, in the current study, we get consistent values. TSM for the eclosion from this paper is different compared to the general model for the total immature development given in our previous paper (Table1 in Gruszka & Matuszewski15). The previous model was created without data from 17 °C and 22 °C, as regression analysis showed that they were outside the 95% confidence interval. In the current study, we excluded the data from 17 °C before the analysis and therefore they were not used while calculating initial TSMs (we give the reasons for this in the results section and discuss proximate and distant causes later in this section). Omitting data from 17 °C changed the results of the regression analysis, and the data point from 22 °C was now within the 95% confidence interval, so it was not omitted while deriving the current final TSMs. Therefore, current article gives more accurate thermal summation values for the eclosion ofN. littoralis and we encourage using the current TSM for the eclosion instead of the less accurate previous one.

The value of thermal summation constant (K) for the eclosion ofN. littoralis is larger compared to those obtained for other silphid beetles (Table4). Larger accumulation of thermal units needed to reach maturity by an insect species indicates its longer development time. The development ofN. littoralis is the longest among studied silphid beetles. Furthermore, this species usually colonizes cadavers later than for instanceThanatophilus species59,60. Therefore,Necrodes littoralis expands the timeframe in which PMI may be estimated using carrion beetles. In other beetle families, the thermal summation values are larger than in carrion beetles (Table4).

Table 4.

Comparison of thermal summation models for the eclosion created for forensically important beetles.

SpeciesFamilyGeographic populationTemp. range [°C]Thermal summation constant [ADD]Developmental threshold [°C]References
Necrodes littoralisSilphidaeCentral European14–30434.70 ± 28.869.04 ± 0.55This study
Thanatophilus micansSilphidaeSouth African17–20197.97 ± 19.7413.26 ± 0.58Ridgeway et al.37
Thanatophilus mutilatusSilphidaeSouth African15–27.5384.11 ± 16.999.04 ± 0.36Ridgeway et al.37
Thanatophilus sinuatusSilphidaeCentral European14–26360.46 ± 10.759.85 ± 0.23Montoya-Molina et al.39
Thanatophilus rugosusSilphidaeCentral European12–22362.76 ± 4.978.53 ± 0.08Montoya-Molina et al.40
Dermestes tessellatocollisDermestidaeChinese16–25664.39 ± 55.8712.07 ± 0.53Wang et al.42
Creophilus maxillosusStaphylinidaeChinese17.5–30492.06 ± 23.619.60 ± 0.58Wang et al.44
Central European15–30405.16 ± 14.6311.66 ± 0.24Frątczak-Łagiewska et al.45
Necrobia rufipesCleridaeChinese22–36591.00 ± 39.5316.62 ± 0.63Hu et al.46
Omosita colonNitidulidaeChinese16–31514.1 ± 8.710.65 ± 0.16Wang et al.47

The pattern of mortality across temperatures is consistent with the patterns revealed for other forensically important beetles. InT. micans,T sinuatus,D. tessellatocollis,C. maxillosus andN. rufipes the highest mortality was observed at extreme temperatures35,39,42,45,46. The same pattern is reported in this study. At intermediate temperatures, mortality ofN. littoralis was below 25%, which is also in line with mortality patterns reported for the other beetle species42,45. In the study of Dekeirsschieter at 23 °C, mortality ofN. littoralis was surprisingly high. It amounted to 30%58, whereas only 3.75% of the beetles died in 22 °C in the current study. Probably the higher mortality reported by Dekeirsschieter resulted from individual rearing of larvae, since mortality rate of individually reared larvae ofN littoralis was found to be higher compared to the larvae that were reared in aggregations61. However, other effects might have been important, as well.

Mortality was exceptionally high at 17 °C in this study. More than half of the beetles died. Moreover, the development at this temperature deviated significantly from the general pattern, with the first and second instar larvae developing longer than at lower temperatures (16 °C and 15 °C). For this reason, we decided not to use data from 17 °C to calculate TSMs. Development extension and high mortality were probably caused here by rearing problems and were not directly related to the temperature itself. At 17 °C large numbers of nematodes were recorded in all containers. They were sticking to the larvae, which probably interrupted larval development and finally led to death of some larvae.

In order to create an isomegalen diagram and growth curves, it is necessary to measure larvae throughout their development. For this purpose, larvae can be killed and then measured. However, killing and storage methods may deform larvae62. Therefore, we decided to measure larvae in vivo. This method also has some drawbacks. First, the temperature changes when larvae are removed from the temperature chambers for measurement. Second, by taking a measurement, stress is induced that in turn can cause an increase in juvenile hormone levels, which interfere with the development63. Previous study64 revealed differences in development time between measured and non-measuredCreophilus maxillosus beetles. However, after TSMs were built, it turned out that the estimation error did not differ significantly between the two models. The authors concluded that differences in development time were an effect of repeated stress rather than temperature changes during the measurements. Moreover, lack of significant differences between measured and non-measured TSMs was due to a small size of the differences and moderate size of the validation sample. In the current study, differences in the estimation error between the models were also insignificant. Based on this finding it can be concluded that in vivo measurements of forensically useful beetles have negligible impact on the resultant development models.

Third instar larvae in this study had longer maximum body length than those in previous studies ofN. littoralis58,65. In the studies of Dekeirsschieter58, larvae were reared individually on small pieces of meat that could result in their smaller size. In the second study, larvae were collected from pig carcasses on the 17th, 20th and 24th days of decomposition65. The largest third instar larvae used in that study were therefore not the largest (longest) larvae present during decomposition, since after 24th day of decay larvae continued to feed and grow. In the current study, a substantial increase in larval length (above 22 mm) took place at the very end of the feeding phase (Supplementary Figs. 3–10). Fully mature Silphinae larvae range in size from 12 to 40 mm66, current results on larval length are therefore consistent with this range.

When rearing larvae in aggregations, it is difficult to ensure that manual samples are taken at random. We tried to sample larvae irrespective of their position inside the aggregation and their size. However owing to a natural tendency to pick larger larvae, there may be a size bias in our analysis.

A transition of the third instar larva to the post-feeding phase is a critical moment inN. littoralis development. This is when the larva buries itself into the soil to form a pupal chamber. The transition was impossible in the containers in which larvae were reared in aggregations. Therefore, it was necessary to transfer them to separate, smaller containers after they ceased feeding. This transfer might slightly interfere with development. Identification of the exact moment of transition to the post-feeding phase is, however, difficult due to rather long duration and lack of unambiguous markers of transition. The same problem was reported in the study ofD. tessellatocollis development, in case of which durations of the last larval stage and the prepupal stage were summed up due to the difficulties in identifying the transition point42. In the study onT. sinuatus, T. rugosus andN. brunnicollis, the post-feeding stage was distinguished38–40. Similarly, in the studies ofN. littoralis by Dekeirsschieter, the post-feeding stage was identified, supposedly due to the rearing of larvae individually in Petri dishes that facilitated monitoring of the larvae58.

Necrodes littoralis is a species with great forensic potential. However, the lack of development models has largely limited its use in forensic cases. By providing the first comprehensive developmental dataset for this species, current study makes a significant contribution to the advancement of forensic entomology, particularly in central and northern Europe.

Materials and methods

Laboratory colony

Adult beetles came from our laboratory colony established in 2017 using insects collected in the alder forest of Biedrusko military range (Western Poland, Central Europe, 52° 31ʹ N, 16° 54ʹ E). They were kept at room temperature in medium-sized terrariums, 20–30 adult beetles per box (2–3 boxes maintained simultaneously). Boxes contained soil and cotton wool soaked with water. Beetles were fed with pork ad libitum.

Experimental rearing

Laboratory rearing was conducted at ten constant temperatures: 14, 15, 16, 17, 18, 19, 20, 22, 26 and 30 °C. We followed the same protocol for each of the tested temperatures. Rearing was carried out inside the temperature chambers (type ST 1/1 BASIC or ST 1/ 1 + , POL-EKO, Poland). We controlled the chambers using temperature recorders. Deviations from the settings were consistent with the manufacturer's assurances (they did not exceed ± 0.5 °C). In order to induce oviposition adult beetles from the main colony were matched in pairs and placed in 0.5-L containers filled with soil. They had constant access to water and pork ad libitum. Containers were inspected every four hours for the presence of fresh eggs. If eggs were recorded, adult beetles were taken out from the container and eggs were left for hatching. Further inspections were carried at intervals representing no more than 10% of the egg stage duration. Inspection times were calculated based on the results of pilot studies. After hatching, first instar larvae were counted, and 50 larvae were placed in a small terrarium (18 cm × 11 cm × 14 cm) for further rearing. There were eight replicates per temperature (8 containers with 50 larvae each). Rearing boxes contained soil, cotton wool soaked with water, and pork meat ad libitum. Meat was covered with aluminum foil to avoid drying out. After feeding was completed and larvae started to bury themselves, they were transferred to new 0.5-L containers (8–10 larvae per container), filled with soil to allow them to form pupal chambers and complete their development.

Inspections and measurements

Immature beetles from every terrarium were checked for transition to the next developmental stage: second instar larva, third instar larva, post-feeding larva, pupa and adult beetle. Inspections were carried out at time intervals that were no longer than 10% of the stage duration. In the case of the second instar, third instar and post-feeding larvae data were collected per container. Transitions were identified when more than half of the larvae passed to the next stage. Data for pupae and adult beetles were collected per insect. Mortality was measured for each stage of development. For this purpose, live individuals were counted after transition to the next stage. In addition, 24 sampled larvae from four out of eight containers (six larvae per container) were measured in vivo during each inspection. Length of their body was measured from the anterior margin of clypeus to the posterior end of the last abdominal segment. Measurements were made using a geometrical micrometer67. Larvae in the other four containers were inspected for transitions and mortality only.

Data analysis

Influence of in vivo measurements

To test the influence of multiple in vivo measurements on the development ofN. littoralis, we compared the total development time of measured and non-measured beetles. 100 randomly selected individuals were used per temperature (54 beetles at 14 °C due to high mortality). Non-parametric Mann–Whitney U test was used in this comparison. To account for the multiple tests we used Bonferonni correction (since there were 10 tests, we used 0.005 level of significance). Then, using data for measured and non-measured beetles, we created separate thermal summation models with the Ikemoto and Takai method11. In order to validate these models we used 110 non-measured beetles (15 beetles per temperature; insects from 17 °C were omitted; due to the high mortality, we had no specimens from 14 °C and only 5 beetles from 30 °C). The true physiological age of the insects (known based on laboratory rearing) was compared with the age estimated using the models. The differences between true and estimated age were used to quantify the estimation errors for measured and non-measured TSM. These errors were compared using the Wilcoxon signed-rank test. Analyzes were performed with Statistica 13 (TIBCO Software Inc.).

Mortality and development

We calculated the percentage mortality for each developmental stage and each of the temperatures tested. Due to the difficulty in accurate separation of the active feeding and post-feeding third instar larvae, we summed up the mortality for these stages. They were similarly combined for the purpose of all the subsequent analyzes. Using data for all individuals, we calculated the mean duration for each developmental stage at each temperature. Median times to reach developmental landmarks (i.e. hatching, first ecdysis, second ecdysis, pupation and eclosion) were calculated and used to build the isomorphen diagram. Using larval measurements, we built isomegalen diagram and growth curves. Thermal summation models were derived using the Ikemoto and Takai method11. For this purpose, we used measured and non-measured beetles. For hatching, first ecdysis and second ecdysis, the sample comprised eight observations per temperature, since we had eight containers for each temperature and could use only “container data” in these analyzes. For pupation and eclosion, the sample comprised 100 observations per temperature, since we could use “individual data” in these analyzes (for 14 °C the sample comprised 63 beetles upon pupation and 54 beetles upon eclosion). Growth curves and thermal summation models were created in R 3.5.2. The other analyzes were performed using Statistica 13 (TIBCO Software Inc.).

Ethical approval

The study comprised laboratory experiments using insect speciesNecrodes littoralis (Coleoptera: Silphidae). The species is not under protection. No permission or approval from Ethic Commission were needed.

Supplementary Information

Acknowledgements

The study was funded by the National Science Center of Poland (grant no. 2016/21/B/NZ8/00788).

Author contributions

J.G. performed laboratory rearing and data collection, analyzed the results and prepared figures. S.M. obtained funds and supervised the research. Both authors designed the study and wrote the manuscript.

Data availability

The datasets generated and/or analyzed during the study are available from the corresponding author on a reasonable request.

Competing interests

The authors declare no competing interests.

Footnotes

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-022-13901-y.

References

  • 1.Smith KGV. Manual of Forensic Entomology. Trustees of the British Museum; 1986. [Google Scholar]
  • 2.Catts EP, Goff ML. Forensic entomology in criminal investigations. Annu. Rev. Entomol. 1992;37:253–272. doi: 10.1146/annurev.en.37.010192.001345. [DOI] [PubMed] [Google Scholar]
  • 3.Charabidze D, Gosselin M, Hedouin V. Use of necrophagous insects as evidence of cadaver relocation: Myth or reality? PeerJ. 2017;5:e3506. doi: 10.7717/peerj.3506. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 4.Matuszewski S, Szafałowicz M, Jarmusz M. Insects colonising carcasses in open and forest habitats of Central Europe: Search for indicators of corpse relocation. Forens. Sci. Int. 2013;231:234–239. doi: 10.1016/j.forsciint.2013.05.018. [DOI] [PubMed] [Google Scholar]
  • 5.Tarone AM, Sanford MR. Is PMI the hypothesis or the null hypothesis? J. Med. Entomol. 2017;54:1109–1115. doi: 10.1093/jme/tjx119. [DOI] [PubMed] [Google Scholar]
  • 6.Grassberger M, Frank C. Initial study of arthropod succession on pig carrion in a central European urban habitat. J. Med. Entomol. 2004;41:511–523. doi: 10.1603/0022-2585-41.3.511. [DOI] [PubMed] [Google Scholar]
  • 7.Matuszewski S, Bajerlein D, Konwerski S, Szpila K. An initial study of insect succession and carrion decomposition in various forest habitats of Central Europe. Forens. Sci. Int. 2008;180:61–69. doi: 10.1016/j.forsciint.2008.06.015. [DOI] [PubMed] [Google Scholar]
  • 8.Amendt J, et al. Best practice in forensic entomology—standards and guidelines. Int. J. Legal Med. 2007;121:90–104. doi: 10.1007/s00414-006-0086-x. [DOI] [PubMed] [Google Scholar]
  • 9.Higley LG, Haskell N. In: Forensic Entomology. The Utility of Arthropods in Legal Investigations. Byrd JH, Castner JL, editors. CRC Press; 2010. pp. 389–407. [Google Scholar]
  • 10.Amendt J, Richards CS, Campobasso CP, Zehner R, Hall MJR. Forensic entomology: Applications and limitations. Forens. Sci. Med. Pathol. 2011;7:379–392. doi: 10.1007/s12024-010-9209-2. [DOI] [PubMed] [Google Scholar]
  • 11.Ikemoto T, Takai K. A new linearized formula for the law of total effective temperature and the evaluation of line-fitting methods with both variables subject to error. Environ. Entomol. 2000;29:671–682. doi: 10.1603/0046-225x-29.4.671. [DOI] [Google Scholar]
  • 12.Picard CJ, et al. Increasing precision in development-based postmortem interval estimates: What's sex got to do with it? J. Med. Entomol. 2013;50:425–431. doi: 10.1603/me12051. [DOI] [PubMed] [Google Scholar]
  • 13.Frątczak-Łagiewska K, Matuszewski S. Sex-specific developmental models for Creophilus maxillosus (L.) (Coleoptera: Staphylinidae): Searching for larger accuracy of insect age estimates. Int. J. Legal Med. 2018;132:887–895. doi: 10.1007/s00414-017-1713-4. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 14.Matuszewski S, Frątczak-Łagiewska K. Size at emergence improves accuracy of age estimates in forensically-usefull beetle Creophilus maxillosus L. (Staphylinidae) Sci. Rep. 2018 doi: 10.1038/s41598-018-20796-1. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 15.Gruszka J, Matuszewski S. Estimation of physiological age at emergence based on traits of the forensically useful adult carrion beetle Necrodes littoralis L. (Silphidae) Forens. Sci. Int. 2020;314:110407. doi: 10.1016/j.forsciint.2020.110407. [DOI] [PubMed] [Google Scholar]
  • 16.Matuszewski S. Post-mortem interval estimation based on insect evidence: Current challenges. Insects. 2021;12:314. doi: 10.3390/insects12040314. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 17.Bajerlein D, Taberski D, Matuszewski S. Estimation of postmortem interval (PMI) based on empty puparia of Phormia regina (Meigen) (Diptera: Calliphoridae) and third larval stage of Necrodes littoralis (L.) (Coleoptera: Silphidae)—advantages of using different PMI indicators. J. Forens. Legal Med. 2018;55:95–98. doi: 10.1016/j.jflm.2018.02.008. [DOI] [PubMed] [Google Scholar]
  • 18.Greenberg B. Flies as forensic indicators. J. Med. Entomol. 1991;28:565–577. doi: 10.1093/jmedent/28.5.565. [DOI] [PubMed] [Google Scholar]
  • 19.Catts EP. Problems in estimating the postmortem interval in death investigations. J. Agric. Entomol. 1992;9:245–255. [Google Scholar]
  • 20.Payne JA. A summer carrion study of the baby pig Sus scrofa Linnaeus. Ecology. 1965;46:592–602. doi: 10.2307/1934999. [DOI] [Google Scholar]
  • 21.Ives AR. Aggregation and coexistence in a carrion fly community. Ecol. Monogr. 1991;61:75–94. doi: 10.2307/1943000. [DOI] [Google Scholar]
  • 22.Slone DH, Gruner SV. Thermoregulation in larval aggregations of carrion-feeding blow flies (Diptera: Calliphoridae) J. Med. Entomol. 2007;44:516–523. doi: 10.1603/0022-2585(2007)44[516:tilaoc]2.0.co;2. [DOI] [PubMed] [Google Scholar]
  • 23.Kulshrestha P, Satpathyb DK. Use of beetles in forensic entomology. Forens. Sci. Int. 2001;120:15–17. doi: 10.1016/s0379-0738(01)00410-8. [DOI] [PubMed] [Google Scholar]
  • 24.Matuszewski S, Bajerlein D, Konwerski S, Szpila K. Insect succession and carrion decomposition in selected forests of Central Europe. Part 2: Composition and residency patterns of carrion fauna. Forens. Sci. Int. 2010;195:42–51. doi: 10.1016/j.forsciint.2009.11.007. [DOI] [PubMed] [Google Scholar]
  • 25.Kadej M, Szleszkowski Ł, Thannhäuser A, Jurek T. A mummified human corpse and associated insects of forensic importance in indoor conditions. Int. J. Legal Med. 2020;134:1963–1971. doi: 10.1007/s00414-020-02373-2. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 26.Mashaly A, Al-Khalifa M, Al-Qahtni A, Alshehri A. Analysis of insects colonised on human corpses during autopsy in Riyadh Saudi Arabia. Entomol. Res. 2020;50:351–360. doi: 10.1111/1748-5967.12441. [DOI] [Google Scholar]
  • 27.Meira LMR, Barbosa TM, Jales JT, Santos AN, Gama RA. Insects associated to crime scenes in the northeast of Brazil: Consolidation of collaboration between entomologists and criminal investigation institutes. J. Med. Entomol. 2020;57:1012–1020. doi: 10.1093/jme/tjaa040. [DOI] [PubMed] [Google Scholar]
  • 28.Wang M, et al. Forensic entomology application in China: Four case reports. J. Forensic Leg. Med. 2019;63:40–47. doi: 10.1016/j.jflm.2019.03.001. [DOI] [PubMed] [Google Scholar]
  • 29.Moemenbellah-Fard MD, Keshavarzi D, Fereidooni M, Soltani A. First survey of forensically important insects from human corpses in Shiraz, Iran. J. Forens. Legal Med. 2018;54:62–68. doi: 10.1016/j.jflm.2017.12.016. [DOI] [PubMed] [Google Scholar]
  • 30.Bonacci T, Vercillo V, Benecke M. Dermestes frischii and D. undulatus (Coleoptera: Dermestidae) on a Human Corpse in Southern Italy: First Report*. Roman. J. Legal Med. 2017;25:180–184. doi: 10.4323/rjlm.2017.180. [DOI] [Google Scholar]
  • 31.Charabidze D, Vincent B, Pasquerault T, Hedouin V. The biology and ecology of Necrodes littoralis, a species of forensic interest in Europe. Int. J. Legal Med. 2016;130:273–280. doi: 10.1007/s00414-015-1253-8. [DOI] [PubMed] [Google Scholar]
  • 32.Charabidze D, Colard T, Vincent B, Pasquerault T, Hedouin V. Involvement of larder beetles (Coleoptera: Dermestidae) on human cadavers: a review of 81 forensic cases. Int. J. Legal Med. 2014;128:1021–1030. doi: 10.1007/s00414-013-0945-1. [DOI] [PubMed] [Google Scholar]
  • 33.Arnaldos MI, Garca MD, Romera E, Presa JJ, Luna A. Estimation of postmortem interval in real cases based on experimentally obtained entomological evidence. Forens. Sci. Int. 2005;149:57–65. doi: 10.1016/j.forsciint.2004.04.087. [DOI] [PubMed] [Google Scholar]
  • 34.Matuszewski S, Mądra-Bielewicz A. Post-mortem interval estimation based on insect evidence in a quasi-indoor habitat. Sci. Justice. 2019;59:109–115. doi: 10.1016/j.scijus.2018.06.004. [DOI] [PubMed] [Google Scholar]
  • 35.Midgley JM, Villet MH. Development of Thanatophilus micans (Fabricius 1794) (Coleoptera: Silphidae) at constant temperatures. Int. J. Legal Med. 2009;123:103–108. doi: 10.1007/s00414-008-0280-0. [DOI] [PubMed] [Google Scholar]
  • 36.Velásquez Y, Viloria AL. Effects of temperature on the development of the Neotropical carrion beetle Oxelytrum discicolle (Brullé, 1840) (Coleoptera: Silphidae) Forens. Sci. Int. 2009;185:107–109. doi: 10.1016/j.forsciint.2008.12.020. [DOI] [PubMed] [Google Scholar]
  • 37.Ridgeway JA, Midgley JM, Collett IJ, Villet MH. Advantages of using developmental models of the carrion beetles Thanatophilus micans (Fabricius), and T. mutilatus (Castelneau) (Coleoptera: Silphidae) for estimating minimum post mortem intervals, verified with case data. Int. J. Legal Med. 2014;128:207–220. doi: 10.1007/s00414-013-0865-0. [DOI] [PubMed] [Google Scholar]
  • 38.Jakubec P, Qubaiová J, Novák M, Růžička J. Developmental Biology of Forensically Important Beetle, Necrophila (Calosilpha) brunnicollis (Coleoptera: Silphidae) J. Med. Entomol. 2020;58:64–70. doi: 10.1093/jme/tjaa170. [DOI] [PubMed] [Google Scholar]
  • 39.Montoya-Molina S, et al. Developmental Models of the Forensically Important Carrion Beetle, Thanatophilus sinuatus (Coleoptera: Silphidae) J. Med. Entomol. 2020 doi: 10.1093/jme/tjaa255. [DOI] [PubMed] [Google Scholar]
  • 40.Montoya-Molina S, et al. Developmental models of the carrion beetle Thanatophilus rugosus (Linnaeus, 1758) (Coleoptera: Silphidae) Sci. Rep. 2021 doi: 10.1038/s41598-021-98833-9. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 41.Martín-Vega D, Díaz-Aranda LM, Baz A, Cifrián B. Effect of temperature on the survival and development of three forensically relevant Dermestes species (Coleoptera: Dermestidae) J. Med. Entomol. 2017;54:1140–1150. doi: 10.1093/jme/tjx110. [DOI] [PubMed] [Google Scholar]
  • 42.Wang Y, et al. Development of Dermestes tessellatocollis Motschulsky under different constant temperatures and its implication in forensic entomology. Forens. Sci. Int. 2021 doi: 10.1016/j.forsciint.2021.110723. [DOI] [PubMed] [Google Scholar]
  • 43.Jakubec P. Thermal summation model and instar determination of all developmental stages of necrophagous beetle, Sciodrepoides watsoni (Spence) (Coleoptera: Leiodidae: Cholevinae) PeerJ. 2016 doi: 10.7717/peerj.1944. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 44.Wang Y, et al. Development of the forensically important beetle Creophilus maxillosus (Coleoptera: Staphylinidae) at constant temperatures. J.Med. Entomol. 2016 doi: 10.1093/jme/tjw193. [DOI] [PubMed] [Google Scholar]
  • 45.Frątczak-Łagiewska K, Grzywacz A, Matuszewski S. Development and validation of forensically useful growth models for Central European population of Creophilus maxillosus L. (Coleoptera: Staphylinidae) Int. J. Legal Med. 2020 doi: 10.1007/s00414-020-02275-3. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 46.Hu G, et al. Development of Necrobia rufipes (De Geer, 1775) (Coleoptera: Cleridae) under constant temperatures and its implication in forensic entomology. Forens. Sci. Int. 2020;311:110275. doi: 10.1016/j.forsciint.2020.110275. [DOI] [PubMed] [Google Scholar]
  • 47.Wang Y, et al. Temperature-dependent development of Omosita colon at constant temperature and its implication for PMI min estimation. J. Forens. Leg. Med. 2020;72:101946. doi: 10.1016/j.jflm.2020.101946. [DOI] [PubMed] [Google Scholar]
  • 48.Anton E, Niederegger S, Beutel RG. Beetles and flies collected on pig carrion in an experimental setting in Thuringia and their forensic implications. Med. Vet. Entomol. 2011;25:353–364. doi: 10.1111/j.1365-2915.2011.00975.x. [DOI] [PubMed] [Google Scholar]
  • 49.Bonacci T, et al. First report of the presence of Necrodes littoralis (L) (Coleoptera: Silphidae) on a human corpse in Italy. J. Forens. Sci. 2021;66:2511–2514. doi: 10.1111/1556-4029.14821. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 50.Matuszewski S, Bajerlein D, Konwerski S, Szpila K. Insect succession and carrion decomposition in selected forests of Central Europe. Part 1: Pattern and rate of decomposition. Forens. Sci. Int. 2010;194:85–93. doi: 10.1016/j.forsciint.2009.10.016. [DOI] [PubMed] [Google Scholar]
  • 51.Matuszewski S, Konwerski S, Frątczak K, Szafałowicz M. Effect of body mass and clothing on decomposition of pig carcasses. Int. J. Legal Med. 2014;128:1039–1048. doi: 10.1007/s00414-014-0965-5. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 52.Gruszka J, et al. Patterns and mechanisms for larval aggregation in carrion beetle Necrodes littoralis (Coleoptera: Silphidae) Anim. Behav. 2020;162:1–10. doi: 10.1016/j.anbehav.2020.01.011. [DOI] [Google Scholar]
  • 53.Saloña-Bordas MI, Perotti MA. First contribution of mites (Acari) to the forensic analysis of hanged corpses: A case study from Spain. Forens. Sci. Int. 2014 doi: 10.1016/j.forsciint.2014.08.005. [DOI] [PubMed] [Google Scholar]
  • 54.Dekeirsschieter J, Frederickx C, Verheggen F, Boxho P, Haubruge E. Forensic entomology investigations from doctor marcel leclercq (1924–2008): A review of cases from 1969 to 2005. J. Med. Entomol. 2013;50:935–954. doi: 10.1603/me12097. [DOI] [PubMed] [Google Scholar]
  • 55.Lutz L, Zehner R, Verhoff MA, Bratzke H, Amendt J. It is all about the insects: A retrospective on 20 years of forensic entomology highlights the importance of insects in legal investigations. Int. J. Legal Med. 2021;135:2637–2651. doi: 10.1007/s00414-021-02628-6. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 56.Matuszewski S, Szafałowicz M. Temperature-dependent appearance of forensically useful beetles on carcasses. Forens. Sci. Int. 2013;229:92–99. doi: 10.1016/j.forsciint.2013.03.034. [DOI] [PubMed] [Google Scholar]
  • 57.Matuszewski S, Mądra-Bielewicz A. Validation of temperature methods for the estimation of pre-appearance interval in carrion insects. Forens. Sci. Med. Pathol. 2016;12:50–57. doi: 10.1007/s12024-015-9735-z. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 58.Dekeirsschieter, J.Etude des interactions entre l’entomofaune et un cadavre:approches biologique,comportementale et chémo-écologique du coléoptère nécrophage,Thanatophilus sinuatus Fabricius (Col., Silphidae) PhD thesis, Gembloux Agro-Bio Tech - University of Liège, (2012).
  • 59.Matuszewski S, Bajerlein D, Konwerski S, Szpila K. Insect succession and carrion decomposition in selected forests of Central Europe. Part 3: Succession of carrion fauna. Forens. Sci. Int. 2011;207:150–163. doi: 10.1016/j.forsciint.2010.09.022. [DOI] [PubMed] [Google Scholar]
  • 60.Matuszewski S, et al. Effect of body mass and clothing on carrion entomofauna. Int. J. Legal Med. 2016;130:221–232. doi: 10.1007/s00414-015-1145-y. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 61.Gruszka J, Matuszewski S. Insect rearing protocols in forensic entomology: benefits from collective rearing of larvae in a carrion beetle Necrodes littoralis L. (Silphidae) PLoS ONE. 2021;16:e0260680. doi: 10.1371/journal.pone.0260680. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 62.Midgley JM, Villet MH. Effect of the killing method on post-mortem change in length of larvae of Thanatophilus micans (Fabricius 1794) (Coleoptera: Silphidae) stored in 70% ethanol. Int. J. Legal Med. 2009;123:103–108. doi: 10.1007/s00414-008-0260-4. [DOI] [PubMed] [Google Scholar]
  • 63.Tauchman SJ, Lorch JM, Orth AP, Goodman WG. Effects of stress on the hemolymph juvenile hormone binding protein titers of Manduca sexta. Insect Biochem. Mol. Biol. 2007;37:847–854. doi: 10.1016/j.ibmb.2007.05.015. [DOI] [PubMed] [Google Scholar]
  • 64.Frątczak-Łagiewska K, Matuszewski S. The quality of developmental reference data in forensic entomology: Detrimental effects of multiple, in vivo measurements in Creophilus maxillosus L. (Coleoptera: Staphylinidae) Forens. Sci. Int. 2019;298:316–322. doi: 10.1016/j.forsciint.2019.02.059. [DOI] [PubMed] [Google Scholar]
  • 65.Frątczak K, Matuszewski S. Instar determination in forensically useful beetles Necrodes littoralis (Silphidae) and Creophilus maxillosus (Staphylinidae) Forens. Sci. Int. 2014;241:20–26. doi: 10.1016/j.forsciint.2014.04.026. [DOI] [PubMed] [Google Scholar]
  • 66.Newton AFJ. In: Immature insects. Stehr FW, editor. Kendall/Hunt; 1991. pp. 339–341. [Google Scholar]
  • 67.Villet MH. An inexpensive geometrical micrometer for measuring small, live insects quickly without harming them. Entomol. Exp. Appl. 2007;122:279–280. doi: 10.1111/j.1570-7458.2006.00520.x. [DOI] [Google Scholar]

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

Data Availability Statement

The datasets generated and/or analyzed during the study are available from the corresponding author on a reasonable request.

Articles from Scientific Reports are provided here courtesy ofNature Publishing Group

ACTIONS

RESOURCES

[8]ページ先頭
©2009-2025 Movatter.jp