Accepted Manuscript Title: Larval development rates of Chrysomya rufifacies Macquart, 1842 (Diptera: Calliphoridae) within its native range in South-East Asia Author: Surasuk Yanmanee Martin Husemann Mark Eric Benbow Guntima Suwannapong PII: DOI: Reference:
S0379-0738(16)30188-8 http://dx.doi.org/doi:10.1016/j.forsciint.2016.04.033 FSI 8449
To appear in:
FSI
Received date: Revised date: Accepted date:
25-3-2016 21-4-2016 27-4-2016
Please cite this article as: S. Yanmanee, M. Husemann, M.E. Benbow, G. Suwannapong, Technical Note, Forensic Science International (2016), http://dx.doi.org/10.1016/j.forsciint.2016.04.033 This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
1 2
Highlights
3
Chrysomya rufifacies represents an indicator species in forensic entomology
5
we investigated the development of the species from the east of Thailand
6
we present isomorphen and isomegalen diagrams for nine temperatures
7
the lower threshold temperature for total development was 9.5°C
cr
ip t
4
8
us
9 10
an
11
Ac ce
pt
ed
M
12
1
Page 1 of 23
Technical Note:
13
Larval development rates of Chrysomya rufifacies Macquart,
14
1842 (Diptera: Calliphoridae) within its native range in South-
15
East Asia
ip t
12
16
Surasuk Yanmanee1, Martin Husemann2, Mark Eric Benbow3, Guntima Suwannapong4*
cr
17
19 20
1
21 22
2
23 24
3
25 26
4
us
18
an
Biological Science Program, Faculty of Science, Burapha University, Chon Buri 20131, Thailand General Zoology, Department of Biology, Martin-Luther University Halle-Wittenberg, 06120 Halle (Saale), Germany
M
Departments of Entomology and of Osteopathic Medical Specialty, Michigan State University, East Lansing, MI 48824, USA
ed
Department of Biology, Faculty of Science, Burapha University, Chon Buri 20131, Thailand
pt
27
*Corresponding author:
29
Guntima Suwannapong
30
Department of Biology
31
Faculty of Science
32
Burapha University
33
Chon Buri 20131
34
Thailand
35
E-mail:
[email protected]
36
Tel. 038-102222 Ext.3088
Ac ce
28
37 38
2
Page 2 of 23
Technical Note:
39
Larval development rates of Chrysomya rufifacies Macquart,
40
1842 (Diptera: Calliphoridae) within its native range in South-
41
East Asia
ip t
38
cr
42
Abstract
44
Chrysomya rufifacies represents an important indicator species in forensic entomology
45
that is often used to estimate the minimum postmortem interval (PMImin) in crime scene
46
investigation. However, developmental rates differ locally, so that estimates should be
47
based on regionally generated development data. Therefore, we determined the
48
developmental rates of C. rufifacies within its native range in Thailand under nine constant
49
temperature regimes: 15, 18, 21, 24, 27, 30, 33, 36 and 39°C. Developmental times from
50
egg to adult varied among the temperatures and were longest at 15°C (618 h) and shortest
51
at 33°C (168 h). No pupae emerged at 39°C. We used linear regression models to estimate
52
the minimum development threshold temperatures for each life stage: egg stage = 9.5°C,
53
first to second instar = 10.8°C, second to third instar = 11.5°C, third instar to pupariation =
54
11.4°C, pupariation to adults = 5.0°C; the minimum threshold to complete all larvae stages
55
was 11.1°C and to complete all life stages from eggs to adult was 9.5°C. We further
56
generated isomorphen and isomegalen diagrams that can be used to quickly estimate the
57
PMImin for forensic applications.
Ac ce
pt
ed
M
an
us
43
58 59
Keywords: Development time, isomegalen diagram, isomorphen diagram, minimum
60
threshold temperature 3
Page 3 of 23
Introduction
62
Forensic entomology is based on the use of developmental and distribution information of
63
necrophagous insects, such as blow flies (Diptera: Calliphoridae), to support crime scene
64
investigations. As insects are poikilotherms, their development rates are strongly
65
temperature dependent and are usually predictable [1]. Further, insects require a certain
66
amount of heat unit-energy units (degree-days) to develop from one life stage to the other
67
[2] and do not develop below or above certain threshold temperatures [3,4]. The lower
68
developmental threshold temperature is most important for PMImin estimation in forensic
69
entomology. Forensic entomologists often use 0°C as lower threshold for estimates based
70
on blow fly species; yet, some reports suggest that temperatures between 6 and 10°C are
71
more appropriate for specific species [5].
M
an
us
cr
ip t
61
72
Generally, the temperature thresholds are species, or even population specific, and should
74
be estimated for each taxon or regional population separately as insects adapt to the local
75
conditions [6]. For instance, it has been reported that Phormia regina (Diptera:
76
Calliphoridae) does not successfully complete development at temperatures below 10°C
77
[7], whereas for Calliphora vicina (Diptera: Calliphoridae) a lower developmental
78
threshold of 1°C was reported [8]. Marchenko (2001) provided lower developmental
79
threshold temperatures for several blow flies and some other fly species further
80
demonstrating interspecific variability: 2.0°C for C. vicina, 3.0°C for Calliphora
81
vomitoria, 7.8°C for Protophormia terraenovae, 7.8°C for Lucilia sericata, 10.2°C for
82
Chrysomya albiceps, 11.4°C for P. regina, 7.2°C for Muscina stabulans, 7.9°C for
83
Muscina assimilis, 7.8°C for Boettcherisca septenrionalis and 6.4°C for Piophila
84
foveolata [9]. Only few studies addressed variability of thermal thresholds within blow
85
flies species; yet, such local estimates may be important, especially for geographically
Ac ce
pt
ed
73
4
Page 4 of 23
86
distant populations [7,10] of widespread taxa. This is specifically critical if widespread
87
species are employed in forensic applications. Hence, studies are needed to better
88
understand regional variability in developmental rates for individual species and their local
89
populations.
ip t
90
Chrysomya rufifacies has a wide distribution and is a blow fly species of special interest
92
for forensic science in many countries of the world, including Thailand [11-16]. The
93
larvae of C. rufifacies are commonly used in forensics [17], but often based on
94
developmental data generated in distant geographic regions. However, in order to more
95
accurately estimate PMImin, local or regional developmental data are needed as
96
populations in different latitudes and climates may vary in development due to plasticity
97
or local adaptation [18]. Although developmental rates of Ch. rufifacies have been studied
98
under ambient temperature in Chiang Mai, northern Thailand [13], the variability of
99
developmental parameters under different temperatures has not been explored in Thailand.
100
Therefore, the objective of this study was to determine the developmental variation of a
101
local population of Ch. rufifacies from central Thailand at a range of temperatures to
102
generate one of the first developmental data sets for this species within its native range in
103
South East Asia.
us
an
M
ed
pt
Ac ce
104
cr
91
105
Material and Methods
106
Study species, sampling and colony establishment and maintenance
107
Adults and larvae of Ch. rufifacies were collected from chicken carcasses in Chon Buri,
108
Thailand between 2012-2015. Adults were identified using the key by Marshall et al. [19]
109
and third instars were identified according to Sukontasan et al. [13] and Smith [20]. One
5
Page 5 of 23
hundred Ch. rufifacies adults (50:50 male and female) were caged (6 x 18 x 13 cm3) under
111
nine different temperature regimes (15, 18, 21, 24, 27, 30, 33, 36 and 39°C) with
112
approximately 60-70% RH and a photoperiod of 12:12 (L.D). Each cage was supplied
113
with water, 50% sucrose and granular sugar ad libitum. Flies were fed with blended swine
114
liver to facilitate ovary development, and on the fourth day were provided ad libitum with
115
beef muscle tissue.
cr
ip t
110
us
116
Egg deposition and larval development
118
To study developmental time from oviposition to hatching at different temperatures (15,
119
18, 21, 24, 27, 30, 33, 36 and 39°C), 100 adult individuals (sex ratio of 1:1) were reared in
120
plastic cages (6 x 18 x 13 cm3). Egg batches were collected within 30 min of oviposition
121
and transferred to fresh beef muscle placed in a plastic container (7.0 x 9.5 x 5 cm3).
122
Specimens were kept at the same condition as described above. For the first 100 hatched
123
larvae the time of emergence was recorded.
M
ed
pt
124
an
117
Another 100 freshly hatched larvae were transferred to 20 g beef muscle placed on an
126
aluminum boat (4 x 5 x 1 cm3) within the plastic container. To evaluate the developmental
127
time from the first to the second instar, from the second to the third instar and from the
128
early third instar to the post feeding stage, 10 larvae from each instar were parboiled
129
(100°C) for 30 s before determining the larval stage under a stereomicroscope [21]. To
130
minimize any potential effects of changing densities, larvae at the exact same stage were
131
replaced from an extra cage keeping the total number of larvae constant at 100 individuals
132
throughout experiment. 20 g of fresh beef muscle were replaced daily. For length
133
measurements, the four largest larvae in each cage were removed from the plastic
Ac ce
125
6
Page 6 of 23
container every three hours following the methods of Grassberger and Reiter [10] and
135
Byrd and Butler [22]. Larvae were parboiled for 30 s and measured [21] immediately after
136
cooling under a binocular microscope in 0.1 mm units using a Vernier caliper [10]. When
137
the third instars entered the post feeding stage, they were moved into 200 g of sterile soil
138
sawdust (50% w/w) in a new plastic container (13 x 18 x 6 cm3) and the time of larval post
139
feeding migration and emergence of each adult was recorded. Three replicates were
140
performed for each rearing temperature.
cr
ip t
134
us
141
Statistical Analyses
143
The relationship between the developmental rate and temperature was modeled according
144
to a linear developmental rate model (ADH) [23]. The lower developmental threshold for
145
each stage was estimated by the x-intercept method [24]. The reciprocal (1/b) of the slope
146
of the regression line represented the sum of degree hours (ADH) defined as the K-value
147
that insects required above the lower developmental threshold [25]. Developmental times
148
and larval body sizes were compared across temperatures using an analysis of variance
149
(ANOVA), and Duncan’s Multiple Rang test (α = 0.05). Statistical analyses were
150
performed with PAST [26].
M
ed
pt
Ac ce
151
an
142
152
Results
153
Development curves from constant temperature regimes
154
The mean maximum body lengths of all instars were plotted against time at each
155
temperature (from hatching until the end of the third instar) (Fig. 1). The growth curve had
7
Page 7 of 23
156
a sigmoidal in form. As expected, larvae grew faster at higher temperatures. The periods
157
of each developmental stage under all temperature regimes are presented in Table 1.
158
Isomegalen- and Isomorphen-diagrams
160
We plotted all developmental data from egg hatching to eclosion at nine constant
161
temperature regimes as isomegalen and isomorphen diagrams (Fig. 2, 3). Isomegalen-
162
diagrams plot the time from egg hatching to post feeding against temperature. Each line
163
represents an identical larval body length (Fig. 2). Further, time of all stages from
164
oviposition to eclosion was plotted against different constant temperatures, where each
165
line represents morphological changes such has molting from one instar to the next. Areas
166
between lines define developmental stages (Fig. 3).
M
an
us
cr
ip t
159
ed
167
Developmental thresholds and thermal constants
169
The developmental rates of all stages increased with temperature (Table S1). The highest
170
developmental rate was found when the eggs were reared at 39°C (0.163 1/hours), while
171
the lowest developmental rate was found for pupariation at 15°C (0.003 1/hours). The
172
developmental rate from oviposition to eclosion was lowest at 15°C and highest at 36°C.
Ac ce
173
pt
168
174
Minimum development thresholds were extrapolated from regression lines for each model
175
(Supplemental Figs. 1, 2, 3). The linear regression of the development from oviposition to
176
pupariation was described as y = 0.0007x - 0.0078 (R2 = 0.99) with a temperature
177
threshold of 11.14 ± 0.58. The rate from oviposition to eclosion was described by the
178
equation y = 0.0002x - 0.0019 (R2 = 0.96) with a minimum thermal threshold of 9.50 ±
8
Page 8 of 23
179
1.04. All minimum threshold temperatures and thermal constants are presented in
180
Supplemental Table S2.
181
Accumulated Degree Hours (ADH)
183
The ADH of Ch. rufifacies from egg to eclosion was highest at 36°C. Different base
184
temperatures showed differences in ADH. All values are provided in Supplemental Tables
185
S3-5.
us
cr
ip t
182
an
186
Discussion
188
In this study we investigated several developmental parameters at different temperatures
189
for the forensically important blow fly species Ch. rufifacies from Thailand. We generated
190
isomegalen and isomorphen diagrams for the use in forensics with the motivation to
191
provide much needed developmental data for this forensically important species from a
192
tropical region.
ed
pt
Ac ce
193
M
187
194
Comparison to other studies of the same species
195
The development times and sizes of peak feeding larvae of Ch. rufifacies we generated
196
here were relatively similar to those reported by Sukontason et. al. [13] at ambient
197
temperatures suggesting consistency of rates within a geographic region and supporting
198
the quality of our data. Byrd and Butler [27] studied the development of eggs, larvae and
199
pupae of Ch. rufifacies under 4 mean cyclic temperatures (15.6, 21.1, 26.7, and 35.0°C)
200
and one constant temperature (25.0°C). The development from egg to adult under all
9
Page 9 of 23
regimes ranged from 190 to 598 h. The constant temperature of 25°C produced a range of
202
pupariation times from 134 to 162 h, with adult emergence ranging from 237 to 289 h.
203
These results are similar to what we report in the present study (see Table 2). No pupae
204
survived at the highest temperature (39°C), indicating that the upper threshold temperature
205
for Ch. rufifacies is at or below this temperature. The upper threshold temperature reported
206
here is lower than the previous upper thermal limits documented for this species, which
207
was 43°C [28]. The developmental time from egg to eclosion at 27°C was similar to the
208
findings of Byrd and Butler [27], but was faster compared to findings of Swiger et al. [29].
209
Whether these differences have a genetic basis or are plastic responses to different
210
environmental conditions needs further evaluation.
an
us
cr
ip t
201
211
Comparison to other taxa
213
Our developmental data for Ch. rufifacies slightly differ from that of other blow fly
214
species; Protophormia terraenovae, for example, has a lower threshold temperature of
215
10.3°C for egg hatching, of 10.7°C for the first instar larval stage, of 10.7-11.0°C for the
216
second instar larval stage, of 11.0°C for the third instar larval stage, of 11.5-11.7°C for
217
pupariation, and of 11.7°C for adult emergence [30]. Further, the lower threshold
218
temperature from oviposition to eclosion in this species differed among geographic
219
regions: for a population from Austria it was estimated at approximately 9°C [10,31],
220
while in a Russian population it was estimated to be 7.8°C [9] and in a population from
221
British Columbia it was 9.8°C [4]. Similar differences between local populations were also
222
found in other species. Davies and Ratcliffe found a threshold of 3.5°C for Calliphora
223
vicina in England [32], which differs from the estimate by Donovan et al. [33], who
224
reported a threshold of 1°C from the same country, whereas Marchenko recorded a base
225
temperature of 2°C for the same species from the Soviet Union [9]. These examples show
Ac ce
pt
ed
M
212
10
Page 10 of 23
226
that temperature threshold may differ even within a country, further supporting the need to
227
explore regional variation and to develop local data for forensically important taxa.
228
Conclusions
230
The developmental rate of Ch. rufifacies was significantly and positively associated with
231
temperature. The longest development times of egg stage, larval stage and pupariation
232
were found at the lowest rearing temperature of 15°C, while the shortest developmental
233
times were found at the higher temperatures. However, no pupal survival was found at the
234
highest temperature (39°C), indicating that the upper threshold temperature for this species
235
in its native range must be below this temperature. Furthermore, larval body size was
236
significantly influenced by temperature where the minimum threshold temperature was
237
estimated to be 11.14°C for the transition from oviposition to pupariation, while it was
238
9.5°C from oviposition to eclosion. The comparison to literature data showed consistency,
239
especially with data from the same region, but also revealed local differences in this and in
240
other species suggesting the need for locally adapted data sets for forensic studies.
pt
ed
M
an
us
cr
ip t
229
Ac ce
241 242
References
243
[1] A. Campbell, B.D. Frazer, N. Gilbert, A.P. Gutierrez, M. Mackauer, Temperature
244
requirements of some aphids and their parasites, J. Appl. Ecol. 11 (1974) 431.
245 246 247 248
[2] K. Bergant, S. Trdan, How reliable are thermal constants for insect development when estimated from laboratory experiments, Entomol. Exp. Appl. 120 (2006) 251–256. [3] B. Greenberg, J.C. Kunich, Entomology and the law: flies as forensic indicators, Cambridge Univ. Press, Cambridge, 2002, pp. 102–103.
11
Page 11 of 23
249
[4] J. Warren, The development of Protophormia terraenovae (Robineau-Desvoidy)
250
(Diptera: Calliphoridae) at constant and fluctuating temperatures, MA. Thesis, Simon
251
Fraser Univ., Burnaby, British Columbia, 2006. [5] L.G. Higley, N.H. Haskell, Insect development and forensic entomology, in: J.H. Byrd,
253
J.L. Castner eds., Forensic entomology: the utility of arthropods in legal investigations,
254
CRC Press, Boca Raton : United States, 2001, pp. 287–302.
257 258
cr
forensically useful flies on carcasses, Int. J. Legal Med. 128 (2014) 1013–1020.
us
256
[6] S. Matuszewski, M. Szafałowicz, A. Grzywacz, Temperature-dependent appearance of
[7] J.H. Byrd, J.C. Allen, The development of the black blow fly, Phormia regina (Meigen), Forensic Sci. Int. 120 (2001) 79–88.
an
255
ip t
252
[8] S.E. Donovan, M.J.R. Hall, B.D. Turner, C.B. Moncrieff, Larval growth rates of the
260
blowfly, Calliphora vicina, over a range of temperatures, Med. Vet. Entomol. 20
261
(2006)
264
ed
263
[9] M.I. Marchenko, Medicolegal relevance of cadaver entomofauna for the determination of the time of death, Forensic Sci. Int. 120 (2001) 89–109. [10] M. Grassberger, C. Reiter, Effect of temperature on development of the forensically
pt
262
M
259
important holarctic blow fly Protophormia terraenovae (Robineau-Desvoidy) (Diptera:
266
Calliphoridae), Forensic Sci. Int. 128 (2002) 177–182.
267
Ac ce
265
[11] K. Sukontason, K.L. Sukontason, S. Piangjai, T. Chaiwong, N. Boonchu, H.
268
Kurahashi, Hairy maggot of Chrysomya villeneuvi (Diptera: Calliphoridae), a fly
269
species of forensic importance, J. Med. Entomol. 40 (2003) 983–984.
270
[12] K.L. Sukontason, K. Sukontason, P. Narongchai, S. Lertthamnongtham, S. Piangjai,
271
J.K. Olson, Chrysomya rufifacies (Macquart) as a forensically-important fly species in
272
Thailand: a case report, J. Vector Ecol. J. Soc. Vector Ecol. 26 (2001) 162–164.
12
Page 12 of 23
273
[13] K. Sukontason, S. Piangjai, S. Siriwattanarungsee, K.L. Sukontason, Morphology and developmental rate of blowflies Chrysomya megacephala and Chrysomya rufifacies in
275
Thailand: application in forensic entomology, Parasitol. Res. 102 (2008) 1207–1216.
276
[14] K.L. Sukontason, K. Sukontason, S. Lertthamnongtham, B. Kuntalue, N. Thijuk, R.C.
277
Vogtsberger, et al., Surface ultrastructure of Chrysomya rufifacies (Macquart) larvae
278
(Diptera: Calliphoridae), J. Med. Entomol. 40 (2003) 259–267.
281
cr
280
[15] J.D. Wells, F.A. Sperling, Molecular phylogeny of Chrysomya albiceps and C. rufifacies (Diptera: Calliphoridae), J. Med. Entomol. 36 (1999) 222–226.
us
279
ip t
274
[16] J.Y. Rosati, S.L. VanLaerhoven, New record of Chrysomya rufifacies (Diptera: Calliphoridae) in Canada: predicted range expansion and potential effects on native
283
species, Can. Entomol. 139 (2007) 670–677.
[17] R.D. Hall, N.H. Haskell, Forensic entomology—applications in medicolegal
M
284
an
282
investigations, in: C.H. Wecht ed., Forensic sciences, Mathew Bender: United States,
286
1995.
287
ed
285
[18] K. Sukontason, P. Narongchai, C. Kanchai, K. Vichairat, P. Sribanditmongkol, T. Bhoopat, et al., Forensic entomology cases in Thailand: a review of cases from 2000 to
289
2006, Parasitol. Res. 101 (2007) 1417–1423.
Ac ce
290
pt
288
[19] S.A. Marshall, T. Whitworth, L. Roscoe, Blow flies (Diptera; Calliphoridae) of
291
eastern Canada with a key to Calliphoridae subfamilies and genera of eastern North
292
America, and a key to the eastern Canadian species of Calliphorinae, Luciliinae and
293
Chrysomyiinae. Can. J. Arthropod Identif. 11 (2011) 1–93.
294 295
[20] K.G.V. Smith, A manual of forensic entomology, Cornell Univ. Press, Ithaca: United States, 1986.
13
Page 13 of 23
296
[21] Z.J.O. Adams, M.J.R. Hall, Methods used for the killing and preservation of blowfly
297
larvae, and their effect on post-mortem larval length, Forensic Sci. Int. 138 (2003) 50–
298
61. [22] J.H. Byrd, J.F. Butler, Effects of temperature on Sarcophaga haemorrhoidalis
300
(Diptera: Sarcophagidae) development, J. Med. Entomol. 35 (1998) 694–698.
303
vicina, Z. Für Rechtsmed. 91 (1984) 295–308.
cr
302
[23] C. Reiter, Zum Wachstumsverhalten der Maden der blauen Schmeißfliege Calliphora
[24] G.D. Broufas, D.S. Koveos, Threshold temperature for post-diapause development
us
301
ip t
299
and degree-days to hatching of winter eggs of the European red mite (Acari:
305
Tetranychidae) in Northern Greece, Environ. Entomol. 29 (2000) 710–713.
306
an
304
[25] B. Bourel, B. Callet, V. Hédouin, D. Gosset, Flies eggs: a new method for the estimation of short-term post-mortem interval? Forensic Sci. Int. 135 (2003) 27–34.
308
[26] Ø. Hammer, D.A.T. Harper, P.D. Ryan, PAST: Paleontological statistics software
312 313 314
ed
311
[27] J.H. Byrd, J.F. Butler, Effects of temperature on Chrysomya rufifacies (Diptera: Calliphoridae) development, J. Med. Entomol. 34 (1997) 353–358.
pt
310
package for education and data analysis, Palaeontol. Electron. 4 (2001) 1–9.
[28] A.P. Johnson, J.F. Wallman, Effect of massing on larval growth rate, Forensic Sci.
Ac ce
309
M
307
Int. 241 (2014) 141–149.
[29] S.L. Swiger, J.A. Hogsette, J.F. Butler, Larval distribution and behavior of
315
Chrysomya rufifacies (Macquart) (Diptera: Calliphoridae) relative to other species on
316
Florida Black Bear (Carnivora: Ursidae) decomposing carcasses, Neotrop. Entomol. 43
317
(2014) 21–26.
318
[30] G.S. Anderson, J.A. Warren, Establishing lower developmental thresholds for a
319
common blowfly: for use in estimating elapsed time since death using entomological
320
methods. Contract Report: Defense Research and Development Canada, 2011.
14
Page 14 of 23
321
[31] M. Grassberger, C. Reiter, Effect of temperature on Lucilia sericata (Diptera:
322
Calliphoridae) development with special reference to the isomegalen- and isomorphen-
323
diagram, Forensic Sci. Int. 120 (2001) 32–36.
325
[32] L. Davies, G.G. Ratcliffe, Development rates of some pre-adult stages in blowflies with reference to low temperatures, Med. Vet. Entomol. 8 (1994) 245–254.
ip t
324
[33] S.E. Donovan, M.J. Hall, B.D. Turner, C.B. Moncrieff, Larval growth rates of the
327
blow fly, Calliphora vicina, over a range of temperatures. Med. Vet. Entomol. 20
328
(2006) 106-114.
Ac ce
pt
ed
M
an
us
cr
326
15
Page 15 of 23
ip t cr
Table 1. Developmental times (mean ± SD) of different life stages of Ch. rufifacies at nine constant temperatures: egg stage, F8 = 100.83,
330
P