Estimativa de entropia de Muscina stabulans (Fallén) (Diptera, Muscidae) em condições artificiais

BIOLOGIA, ECOLOGIA E DIVERSIDADE

Estimativa de entropia de Muscina stabulans (Fallén) (Diptera, Muscidae) em condições artificiais¹

Rodrigo Ferreira Krüger^I,II; Simone Giehl Erthal^III

^IDepto. de Zoologia e Genética, IB, UFPel, Caixa Postal 354, 96010-900, Pelotas, RS, Brasil. [e-mail].

^IIPrograma de Pós-Graduação em Ciências Biológicas/Entomologia, Depto. de Zoologia, UFPR.

^IIICentro Educacional de Jovens e Adultos Donizete Romualdo da Silva, Rua Serra Azul, 2428, Centro, 78983-000 Espigão do Oeste-RO, Brasil.

RESUMO

Estimativa de entropia de Muscina stabulans (Fallén) (Diptera, Muscidae) em condições artificiais. O conceito de entropia (H) foi adaptado da mecânica estatística para a demografia para quantificar o impacto da mortalidade na expectativa de vida e demonstrar quantitativamente a tendência da mortalidade em populações experimentais. Isto foi verificado para 160 casais de Muscina stabulans (Fallén, 1817) mantidos em câmara climatizada a 24,8ºC ± 0,6ºC, umidade relativa do ar entre 70 e 80% e fotofase de 12 horas. Nestas condições, machos e fêmeas apresentaram valores de H intermediários aos valores teóricos de H = 0 e H = 0,5 demonstrando que para esta espécie, a curva de sobrevivência é do tipo retangular. A distribuição da mortalidade por idade específica indicou que a força desse parâmetro age de dois modos sobre os adultos desta espécie. Em um, a mortalidade tem maior força nos intervalos compreendidos entre a emergência dos adultos e o 10º dia após este processo. No segundo modo, a força de mortalidade é maior entre o 20º e 30º dias após a emergência, sendo que pequenas variações na mortalidade causam maior impacto na sobrevivência das fêmeas do que nos machos.

Palavras-Chave: Expectativa de vida; dinâmica populacional; mortalidade; sobrevivência.

ABSTRACT

Entropy estimation in Muscina stabulans (Fallén) (Diptera, Muscidae) under laboratory conditions. Entropy (H) as a concept in demography was adapted from that of physics to quantify the impact of mortality on life expectancy and to quantitatively demonstrate mortality tendencies in experimental populations. Entropy was estimated for 160 pairs of Muscina stabulans (Fallén, 1817) in a climate-controlled chamber (24.8ºC ± 0.6ºC, relative humidity 70 – 80%, 12 h photophase). Both sexes had H values intermediate to those theoretically expected (0.0 – 0.5), showing that the survival curve is rectangular. The age-specific mortality distribution shows that mortality affects adults of this species in two ways. First, mortality is higher for 10 days from pupal emergence. Second, mortality again increases between 20-30 days after emergence. Small differences in mortality rates cause a greater overall impact on the survival of females.

Keywords: Life expectancy; mortality; population dynamics; survival.

Muscina stabulans (Fallén, 1817) (Diptera, Muscidae) é uma espécie anautógena (Tirone et al. 1996), cosmopolita, com maior freqüência na Região Neotropical (Carvalho 2002). No Brasil foi encontrada na Bahia, Paraná e São Paulo (Carvalho et al. 2002) com preferência por ambientes urbanos (Linhares 1981; Carvalho et al. 1984), importante para a entomologia forense em países subtropicais (Liu & Greenberg 1989; Centeno et al. 2002) e capaz de transmitir diversos patógenos (Greenberg 1971). Em granjas avícolas, podem reduzir drasticamente as populações de Musca domestica L., 1758 pois suas larvas são predadoras facultativas de terceiro instar (Skidmore 1985).

As larvas de M. stabulans utilizam recursos efêmeros e temporalmente irregulares, com maior velocidade de desenvolvimento a 26ºC (Mascarini & Prado 2002). Espécies que utilizam recursos desse tipo servem como modelos para estudos de estrutura e dinâmica populacional (Moura & Bonatto 1999), que podem ser realizados através das tabelas de vida, amplamente utilizadas pois determinam a probabilidade de mortalidade e sobrevivência por idade específica. Tais modelos podem prever um ponto de partida para questões fundamentais sobre processos populacionais, identificando dois parâmetros demográficos importantes, a distribuição da mortalidade dos indivíduos por idade específica e a taxa de crescimento da população, sugerindo ainda, uma direção para esse crescimento (Carey 1993).

A taxa intrínseca de aumento natural (r_m) e a expectativa de vida (e_x) são parâmetros cruciais que descrevem a idade média de morte dos indivíduos em uma população (Demetrius 1979) já que a mortalidade não é igual em todas as idades (Krebs 1994) o que leva a uma distribuição de mortes por idade que poderá assumir dois padrões extremos (Demetrius 1978), sendo possível que todos os indivíduos morram com a mesma idade ou que todos os indivíduos possuam a mesma probabilidade de morrer em cada idade (Carey 1993).

A mortalidade pode ser influenciada pela estrutura genética de uma população e as condições ambientais em que essa população vive, podendo alterar a convexidade da curva de sobrevivência (Demetrius 1979) que tem sido geralmente analisada de forma qualitativa (Moura & Bonatto 1999).

A quantificação dos valores de convexidade desta curva fornece informações úteis sobre a população, como a tendência da mortalidade e a força de mortalidade por idade específica (Demetrius 1979). Isso pode ser obtido por meio da utilização do conceito de entropia (H) em demografia (Demetrius, 1978, 1979; Goldman & Lord 1986).

O conceito de entropia ou informação pode ser interpretado como a medida da heterogeneidade de uma população com respeito à mortalidade em diferentes idades. Este conceito foi adaptado a partir da mecânica estatística e da teoria ergódica que possibilita ainda descrever o desvio da linearidade da variação do número de sobreviventes em função da idade, quantificando assim o impacto da mortalidade sobre a e_x,sendo considerado um sumário estatístico tão significativo quanto r (Demetrius 1978), sintetizando a informação contida na curva de sobrevivência (Tuljapurkar 1982).

Quando H é igual a zero, obtém-se a concavidade mais extrema possível (Carey 1993), uma vez que todos os indivíduos morrem exatamente ao mesmo tempo, em estágios mais avançados de vida, o que dá uma forma retangular à curva de sobrevivência (Vaupel 1986). Já para H igual a um, ocorre um decréscimo geométrico na sobrevivência em relação à idade (Carey 1993), o que corresponderia a uma probabilidade constante de mortalidade em cada idade. Para H igual a 0,5 há uma relação linear entre mortalidade e idade (Demetrius 1978).

A força de mortalidade por idade específica pode ser medida através do produto e_xd_x. Esse produto representa a média ponderada da e_x em idades específicas ou estágios de desenvolvimento. A mortalidade pode ser relacionada com a e_x em idades específicas por intermédio dos valores deste produto que significa os dias que seriam ganhos na e_x desta espécie caso essas mortes fossem evitadas (Carey 1993; Goldman & Lord 1986).

Como observado por Demetrius (1979), os valores de entropia e da força de mortalidade por idade específica podem mudar de acordo com a variabilidade genética da população nos períodos de vida de diferentes indivíduos, mascarando os resultados finais devido à heterogeneidade populacional. Para avaliar a heterogeneidade dos dados relacionados à quantificação das curvas de sobrevivência são necessárias criações de adultos que formem populações independentes.

A dinâmica populacional de M. stabulans tem sido pouco explorada, apesar de sua importância epidemiológica e ecológica. O conceito de entropia foi utilizado neste trabalho para responder três questionamentos relativos a demografia desta espécie. Qual a tendência da curva de sobrevivência de machos e fêmeas? Durante o período de vida de machos e fêmeas, qual ou quais intervalos etários sofrem maior impacto da mortalidade? E, de que forma os resultados obtidos podem contribuir com o entendimento da performance populacional desta espécie?

MATERIAL E MÉTODOS

Foi utilizada uma colônia de M. stabulans mantida em laboratório, em câmara climatizada, a partir de adultos capturados no Campus da Universidade Federal de Pelotas (UFPel) na cidade de Capão do Leão, RS (31º45'48''S e 52º29'02''W), para obtenção de ovos, larvas e pupas. Os adultos foram acondicionados em gaiolas com 30x30x30 cm e alimentados com leite em pó, açúcar refinado e farinha de peixe, numa proporção de 2:2:1, além de água em pequenos frascos com espuma de poliestireno cobrindo a superfície do líquido. A dieta e água foram oferecidas ad libitum conforme estabelecido para outras espécies de Muscidae (Ribeiro et al. 2000; Krüger et al. 2003).

Para obtenção das posturas foi oferecida uma dieta composta de farinha de peixe e serragem, numa proporção de 3:2, levando em consideração dietas utilizadas por outros muscídeos (Ribeiro et al., 2000; Krüger et al. 2003), adicionando-se água até torná-la pastosa. Os ovos foram transferidos para recipientes contendo a mesma dieta utilizada para sua obtenção e colocados em funil para coleta de larvas de terceiro instar na fase pós-alimentar. Neste meio, após a eclosão, as larvas foram alimentadas até o estágio de L_3, abandonando o funil e caindo em recipiente contendo serragem úmida. As L₃ na fase pós-alimentar e as pupas foram transferidas para recipientes de vidro e mantidas até a emergência dos adultos.

A colônia foi mantida durante todo experimento em câmara climatizada com umidade relativa do ar entre 70% e 80% e fotofase de 12 horas. A temperatura foi estimada na câmara da seguinte maneira. As temperaturas mínima e máxima foram estimadas às 7h de cada dia e a temperatura média da câmara durante o período de fotofase foi estimada em três ocasiões diárias (10h, 13h e 18h), gerando uma média de temperatura para este intervalo. A temperatura média do período experimental foi de 24,8ºC ± 0,6ºC, conforme a média diária das três temperaturas, mínima, máxima e média do período de fotofase, com variação total de 22ºC a 28ºC (Tabela I).

Foram mantidas quatro gaiolas, com as mesmas dimensões descritas anteriormente, com 80 indivíduos em cada, sendo 40 casais por gaiola. Para formar os casais, as pupas foram individualizadas em tubos de ensaio com serragem úmida. Após a emergência foi identificado o sexo e posteriormente os casais foram acondicionados em gaiolas, observados diariamente em intervalos de 24 horas, removendo-se as posturas e insetos mortos.

Para o cálculo da estimativa da entropia (H) foram utilizados os seguintes parâmetros, calculados a partir dos dados de sobrevivência em tabela de vida (Carey 1993): x, é a classe etária compreendida em intervalos de cinco dias; l_x, a proporção de sobreviventes em relação ao intervalo x de idade; d_x = l_x l_x+1 fornece a proporção de indivíduos que morreram entre x e x + 1; q_x = d_x/l_x, a probabilidade de morte na idade x e, p_x = 1-q_x, a probabilidade de sobrevivência. A probabilidade média de sobrevivência foi calculada por L_x = l_x (d_x/2) e a esperança de vida (ex) estimada por T_x/l_x.

O número total de dias restantes aos sobreviventes até a idade máxima w, T_x, foi estimado por:

Os valores de entropia (H) foram calculados seguindo-se a expressão: onde e_x representa a esperança de vida na idade x; d_x a proporção de indivíduos que morreram entre as idades x e x+1; e₀a esperança de vida na idade inicial e; w idade máxima alcançada pelo último indivíduo sobrevivente (Carey 1993).

As testemunhas que emergiram da criação em laboratório foram depositadas na Coleção de Insetos de Importância Médica e Veterinária, Departamento de Microbiologia e Parasitologia (CDEMP), Universidade Federal de Pelotas.

RESULTADOS E DISCUSSÃO

As criações tiveram resultados diferentes sob as mesmas condições experimentais devido à variabilidade genética da população de M. stabulans. Esta variação foi sentida nas curvas de sobrevivência para machos e fêmeas das quatro gaiolas bem como nos valores de entropia que estão demonstrados na figura 1 e nos parâmetros necessários para o cálculo da entropia na Tabela II.

Os valores de entropia para machos (0,251; 0,260; 0,214 e 0;280) e fêmeas (0,199; 0,276; 0,225 e 0,284) para as respectivas gaiolas de A a D (Fig. 1) foram intermediários aos valores teóricos de H = 0 e H = 0,5. Com exceção a gaiola A, os valores de H para machos e fêmeas das outras gaiolas foram muito semelhantes, apesar das fêmeas apresentarem valores maiores do que os machos e da e₀ e w diferir bastante nestas condições (Tabela II).

Nos machos e fêmeas a mortalidade em idades avançadas seguiu um padrão verificado para outras espécies de Diptera. Os valores de entropia (fig. 1) foram semelhantes aos obtidos para Chrysomya megacephala (Fabricius, 1794) (Diptera, Calliphoridae) (Von Zuben et al. 1996), Sarconesia chlorogaster (Wiedemann, 1830) (Calliphoridae) (Moura & Bonatto 1999) e Lucilia cuprina (Calliphoridae) (Fernandes et al. 2003).

Como a concavidade da curva l_x pode ser obtida por intermédio de H (Demetrius 1978), os machos de M. stabulans apresentaram l_x levemente mais côncava que a l_x das fêmeas (Fig. 1 B,^C,^D), com exceção da gaiola A, onde ocorreu o inverso (fig. 1A). Para S. chlorogaster e L. cuprina há tendência dos machos possuírem l_x levemente mais côncava (Moura & Bonatto 1999; Fernandes et al. 2003) diferindo do que foi encontrado para C. megacephala (Von Zuben et al. 1996).

Conforme Vaupel (1986), H fornece indícios sobre o efeito da mortalidade por idade específica na e_x através da mudança percentual produzida por uma redução de 1% na força de mortalidade em todas as idades. Com isso, uma redução de 1% na mortalidade de M. stabulans, iria mudar a e_x nos machos de 0,214 a 0,280 dias e nas fêmeas de 0,199 a 0,284 dias sendo semelhante ao que foi obtido para S. chlorogaster (Moura & Bonatto 1999).

O efeito da mortalidade em diferentes idades para esta espécie está demonstrado nas Tabelas III e ^IV através da análise do produto e_xd_x. Os valores desse produto para machos são maiores do 20º ao 30º, 15º ao 25º, 20º ao 30º e do 5º ao 15º dias respectivamente para as gaiolas A, B, C e D (Tabela III). Caso as mortes neste período fossem evitadas haveria adição de quatro dias para as gaiolas A e B e cerca de dois a três dias para as gaiolas C e D. Para as fêmeas estes valores foram de quatro, cinco, três e dois dias para as gaiolas A, B, C e D (^{Tabela IV}). Nessas gaiolas, para fêmeas, estes dias que seriam adicionados à e_x, caso estas mortes fossem evitadas, foram compreendidos pelos intervalos do 25º ao 35º, 20º ao 30º, 5º ao 15º e de zero ao 10º dias. Estes períodos são semelhantes entre machos e fêmeas com exceção da gaiola C (Tabelas III e ^IV).

Em geral, uma história de vida ótima é representada pelo balanceamento entre os custos e benefícios de uma mudança na l_x estando relacionada ao esforço reprodutivo empreendido pelas fêmeas em uma dada idade específica (Reznick 1985). Conforme haja variação na força de mortalidade por idade específica, poderá haver influência significativa na fecundidade das fêmeas de uma população, caso esta força seja maior em intervalos etários onde haja mais oviposições. Isto foi verificado para Ophyra albuquerquei Lopes, 1986 (Diptera, Muscidae), através de análise qualitativa dos dados de sobrevivência e fecundidade e foi constatado que após sucessivos picos de oviposição ocorreram picos de mortalidade de fêmeas, o que poderia indicar um alto custo fisiológico (Krüger et al. 2004). Este custo se dá, geralmente, nos primeiros dias de uma coorte de uma espécie que tenha alta velocidade de desenvolvimento como é o caso de M. stabulans.

Respondendo as questões propostas para este trabalho, machos e fêmeas de M. stabulans apresentam curva de sobrevivência com tendência retangular. Neste caso, a mortalidade tem maior impacto nos primeiros intervalos etários, logo após a emergência dos adultos. Esta força possui dois modos distintos de manifestação nas populações desta espécie nos primeiros dias de vida dos adultos, sendo a variabilidade genética a única resposta possível para este resultado, já que as condições de criação foram iguais para todas as gaiolas.

Mesmo apresentando similaridades quanto à forma da curva de sobrevivência, o impacto da mortalidade é maior nas fêmeas do que nos machos, o que pode alterar os sistemas onde esta espécie é encontrada, conforme preconizado para outras espécies de Diptera (Von Zuben et al. 1993; Moura & Bonatto 1999).

Estes resultados auxiliam no entendimento da performance populacional de M. stabulans, na medida que são as fêmeas os indivíduos responsáveis pela distribuição dos ovos nos substratos de criação das larvas. Logo, se estes indivíduos sofrem maior influência da mortalidade, a probabilidade de localização de novos sítios de criação poderá ser alterada (Roff 1977), já que entre a emergência dos adultos e os períodos de maior influência da mortalidade, decorre pouco tempo, ressaltando que a mortalidade pode ser influenciada pela própria fecundidade através de custos reprodutivos.

Agradecimentos. A CAPES pela bolsa de doutorado ao primeiro autor e a FAPERGS pela bolsa de iniciação científica concedida ao segundo autor. Ao Dr. Paulo Bretanha Ribeiro pelo empréstimo das instalações do Laboratório de Biologia de Insetos da UFPel. Ao Dr. James Roper pelas correções do abstract.

Recebido em 04/07/2005; aceito em 24/04/2006

• Carey, J. R. 1993. Applied demography for biologists New York, Oxford University Press, XVI+260p.
• Carvalho, C. J. B. (ed.). 2002. Muscidae (Diptera) of the Neotropical Region: Taxonomy Curitiba, Ed. UFPR, 287 p.
• Carvalho, C. J. B.; J. R. Almeida & C. B. Jesus. 1984. Dípteros sinantrópicos de Curitiba e arredores (Paraná, Brasil). I. Muscidae. Revista Brasileira de Entomologia 28: 551560.
• Carvalho, C. J. B.; M. O. Moura & P. B. Ribeiro. 2002. Chave para adultos de dípteros (Muscidae, Fanniidae, Anthomyiidae) associados ao ambiente humano no Brasil. Revista Brasileira de Entomologia 46: 107114.
• Centeno, N.; M. Maldonado & A. Oliva. 2002. Seasonal patterns of arthropods occurring on sheltered and unsheltered pig carcasses in Buenos Aires Province (Argentina). Forensic International Science 126: 6370.
• Demetrius, L. 1978. Adaptative value, entropy and survivorship curve. Nature 275: 213214.
• Demetrius, L. 1979. Relations between demographic parameters. Demography 16: 329338.
• Fernandes, F. M.; D. M. Lapola; R. Neregato; M. C. Carvalho & C. J. Von Zuben. 2003. Curva de sobrevivência e estimativa de entropia em Lucilia cuprina (Diptera, Calliphoridae). Iheringia, Ser. Zoologia 93: 319324.
• Goldman, N. & J. W. Lord. 1986. A new look at entropy and the life table. Demography 23: 275282.
• Greenberg, B. 1971. Flies and disease. Ecology, classification and biotic associations. New Jersey, Princeton University Press, Vol. 1, VII+586p.
• Krebs, C. J. 1994. Ecology: The experimental analysis of distribuition and abundance. New York, Harper Collins College Publishers, XIV+801p.
• Krüger, R. F.; P. B. Ribeiro & C. J. B. Carvalho. 2003. Desenvolvimento de Ophyra albuquerquei Lopes (Diptera, Muscidae) em condições de laboratório. Revista Brasileira de Entomologia 47: 643648.
• Krüger, R. F.; P. B. Ribeiro; C. J. B. Carvalho; F. M. Lambrecht & A. M. Nunes. 2004. Longevidade e oviposição de Ophyra albuquerquei (Diptera, Muscidae) em condições de laboratório. Iheringia, Série Zoologia 94: 211216.
• Linhares, A. X. 1981. Synanthropy of Muscidae, Fanniidae and Anthomyiidae (Diptera) in the city of Campinas, São Paulo, Brazil. Revista Brasileira de Entomologia 25: 231243.
• Mascarini, L. M. & A. P. Prado. 2002. Thermal Constant of an Experimental Population of Muscina stabulans (Fallén 1817) (Diptera:Muscidae) in the Laboratory. Memórias do Instituto Oswaldo Cruz 97: 281283.
• Moura, M. O. & S. R. Bonatto. 1999. Análise de sobrevivência e estimativa de entropia para Sarconesia chlorogaster (Wiedemann) (Diptera, Calliphoridae). Revista Brasileira de Zoologia 16: 221226.
• Reznick, D. 1985. Costs of reproduction: an evaluation of the empirical evidence. Oikos 44: 257267.
• Ribeiro, P. B.; C. J. B. Carvalho; P. R. P. Costa & P. Silveira-Júnior. 2000. Desenvolvimento de Ophyra aenescens Wiedemann, 1830 (Diptera, Muscidae, Azeliinae), em diferentes temperaturas, em condições de laboratório. Revista Brasileira de Agrociência 6: 8087.
• Roff, D. 1977. Dispersal in dipterans: its costs and consequences. Journal of Animal Ecology 46: 443456.
• Skidmore, P. 1985. The biology of the Muscidae of the worldDordrecht. Dordrecht Kunk. 550 p.
• Tirone, G.; P. P. Parise & R. M. P. Avancini. 1996. Ovarian development in Muscina stabulans (Diptera, Muscidae) on differents diets. Revista Brasileira de Entomologia 40: 105106.
• Tuljapurkar, S. D. 1982. Why use population entropy? It determines the rate of convergence. Journal Mathematical of Biology 13: 325337.
• Vaupel, J. W. 1986. How change in age-specific mortality affects life expectancy. Populations Studies 40: 147157.
• Von Zuben, C. J.; S. F. Reis; J. B. R. Do Val; W. A. C. Godoy & O. B. Ribeiro. 1993. Dynamics of a mathematical model of Chrysomya megacephala (Diptera, Calliphoridae). Journal of Medical Entomology 30: 443448.
• Von Zuben, C. J.; W. A. C. Godoy & E. L. A. Monteiro-Filho. 1996. Curva de sobrevivência e estimativa de entropia em Chrysomya megacephala (Diptera, Calliphoridae). Revista Brasileira de Entomologia 40: 221 224

1

Contribuição Nº 1561 do Depto. de Zoologia da Universidade Federal do Paraná (UFPR).

Datas de Publicação

• Publicação nesta coleção
  12 Dez 2006
• Data do Fascículo
  Jun 2006

Histórico

• Recebido
  04 Jul 2005
• Aceito
  24 Abr 2006