Tese de mestrado. Bioinformática e Biologia Computacional. Universidade de Lisboa, Faculdade de Ciências, 2010Streptococcus pneumoniae, also known as pneumococcus, is a commensal bacterium very common in the nasopharynx of young children but that can also be found in older children and adults. Carriage may lead to infection. Although this is a rare event, it has a significant impact on human health. Indeed, diseases caused by pneumococcus include infections as common as otitis media and as dangerous as pneumonia or meningitis. It is therefore important to gain more knowledge of its transmission dynamics as modulated by environmental factors and of how it is affected by host population specificities. To characterize pneumococcus transmission dynamics in a Portuguese day-care center, data was used from a one-year longitudinal study on the state of colonization by S. pneumoniae in children attending a day-care center in Lisboa, Portugal [Sá-Leão et al.2008]. The data refer to 1998, before the introduction of the seven-valent pneumococcal conjugate vaccine. A conceptual model for pneumococcus transmission was built, which considered genotype colonizations and clearances as dependent on the number of carriers, the number of non-carriers and the values of four parameters: the clearance rate m, the within-group transmission parameter b, the community rate of acquisition k and the between-genotypes competition parameter f. Bayesian inference was used to estimate these parameters. Colonizations and clearances were modeled as Poisson processes and the joint posterior probability distributions of the model’s parameters were estimated by Markov Chain Monte Carlo sampling. The number of transitions that occurred in each sampling interval was counted directly from the sampled states, assuming that children did not undergo more than one transition per sampling interval. The posterior mean for the transmission parameters were 0.5974for b, 0.0107for k, 0.6280 for f and 0.3059 for m. Data was simulated using the posterior estimates for these parameters from a study of Finnish DCCs [Hoti et al.2009]. Sampling this data monthly, the method was found to give biased estimations, since the assumption that children did not undergo more than one transition per sampling interval did not hold. The precision could be improved by sampling for a longer period, 30 months were used. To significantly improve the accuracy, the sampling interval needed to be extremely short, daily samples were taken from the simulation. The estimation model used was found to be impractical. Another estimation method should be used that infers the possible carriage histories consistent with the observed states.Streptococcus pneumoniae, também denominado pneumococo, é uma bactéria comensal que só coloniza o homem. A aquisição desta bactéria começa pela colonização da nasofaringe, o indivíduo torna-se portador e pode transmitir a bactéria a outras pessoas através de contacto directo. Devido aos contactos frequentes e próximos entre crianças a frequentar infantários e à imaturidade do seu sistema imunitário, as crianças até aos seis anos de idade representam um grupo onde a prevalência de portadores de pneumococo é elevada [Bogaert et al.2004]. Estudos transversais realizados em crianças a frequentar infantários, mostraram que a prevalência de pneumococo nestas crianças pode ser superior a 70% [Mato et al. 2005]. Embora a colonização seja geralmente assintomática, pode, por vezes, evoluír para infecção. Na verdade, o pneumococo é uma causa comum de otite média, pneumonia e meningite, sendo um agente principal de doenças graves em crianças jovens [Centers for Disease Control and Prevention2000]. É, por isso, especialmente importante estudar a dinâmica de transmissão de pneumococo entre as crianças que frequentam infantários, bem como tentar perceber de que forma esta dinâmica pode ser afectada por factores ambientais e características populacionais. Existe uma grande variedade de estirpes de pneumococo, tendo sido descritos mais de noventa de cápsulas, ou serótipos, a que correspondem diferentes propriedades, tais como, um potencial patogénico variável [Brueggemann et al. 2003]. Um estudo recente, usando testes in vitro sobre amostras clínicas, encontrou uma possível relação, com significado biológico, entre o tipo de cápsula e a capacidade de resistir ao sistema imunitário do hospedeiro [Weinberger et al. 2009]. é possível que várias estirpes de pneumococo compitam na nasofaringe. Diversos tipos de competição foram estudados a nível teórico [Lipsitch 1997, Zhang et al. 2004] e os diferentes potenciais de colonização ou resistência a colonização foram testados para alguns serótipos em ratos [Lipsitch et al. 2000]. Um número limitado de serótipos foi estudado e poucos deram diferenças estatisticamente significativas. O primeiro tratamento utilizado contra pneumococo foi a soroterapia em estados iniciais de doença [Klugman 2008]. A partir de uma amostra de expectoração, determinava-se o serótipo e injectava-se no paciente soro animal contendo anticorpos contra esse serótipo. Em 1940 esta terapia estava disponível para cinco serótipos. Entretanto, foram descobertos antibióticos eficazes, sendo de destacar a penicilina. Embora a mortalidade por doença pneumocócica nos muito jovens e idosos tenha continuado alta, a penicilina parecia uma droga milagrosa por ser eficaz contra inúmeras doenças, algumas das quais quase sempre fatais, sem necessidade determinar o serótipo. Com o aumento do uso de penicilina no mundo, os tipos epidémicos mais invasivos, geralmente transmitidos por indivíduos doentes a pessoas saudáveis, mas que raramente se encontravam em indivíduos saudáveis, tornaram-se menos frequentes. Por outro lado, serótipos pediátricos menos virulentos, geralmente colonizadores por períodos prolongados e causa comum de otite média, começaram a substituir os tipos epidémicos nas infecç ões graves e a ganhar resistência aos antibióticos. Finalmente, na década de 80, uma vacina polissacárida eficaz contra as doenças pneumocócicas chegou ao mercado. Infelizmente, esta vacina não era eficiente em crianças. A primeira vacina conjugada a proteger contra a doença e o transporte de sete serótipos (PCV-7), também eficaz em crianças de menos de dois anos, só foi licenciada no ano 2000. Após o uso disseminado da vacina, a doença invasiva diminuiu, embora o transporte de pneumococo se tivesse mantido [Revai et al. 2006]. A última vacina conjugada aprovada protege contra 13 serótipos (PCV-13). O objectivo deste trabalho é estudar a dinâmica de transmissão de S. pneumoniae usando dados de um estudo longitudinal do estado de colonização, por esta bactéria, de crianças que frequentavam um infantário de Lisboa, Portugal [Sá-Leão et al. 2008]. Estes dados foram obtidos em 1998, antes da introdução da vacina PCV7. Ao todo, 21 clones de pneumococo foram identificados correspondentes a 13 serótipos. Neste trabalho, o termo genótipo foi usado para designar o clone. A análise exploratória permitiu perceber que os dados eram insuficientes para modelar a transmissão para genótipos diferentes de forma diferencial. Considerou-se, por isso, que todos os genótipos tinham igual capacidade de ser transmitidos e de resistir à eliminação. Para caracterizar a dinâmica de transmissão de pneumococo no infantário, construiu-se um modelo conceitual da transmissão de pneumococo que pressupõe que todas as crianças amostradas estão em contacto entre si e se misturam, que a colonização de uma criança portadora de pneumococo pode ser diferente da colonização de uma criança não portadora e que as taxas de transmissão e de clearance se mantêm constantes. Embora tendo em conta o número de portadores de cada genótipo, considera-se que todos têm características iguais. Este modelo depende de quatro parâmetros: uma taxa de transmissão entre as crianças do infantário b, uma taxa de aquisição de pneumococo pela comunidade k, um factor de competição entre genótipos f e uma taxa de clearance m. O valor destes parâmetros determina a dinâmica de transmissão própria ao infantário e foi estimado por inferência bayesiana. O modelo de verosimilhança individuais a partir dos estados observados como método alternativo que, por não depender deste pressuposto, permite estimar os parâmetros usando dados amostrados com um intervalo mais viável
