Genetic relationship between reproductive traits in Nellore cattle

Pós-graduação em Genética e Melhoramento Animal - FCAV

Estimativas de componentes de (CO) variância de características de crescimentos na raça Nelore, utilizando inferência bayesiana

Pós-graduação em Genética e Melhoramento Animal - FCAV

Diversidade genética e estrutura populacional da anta (Tapirus terrestris) na Serra do Mar, São Paulo, Brasil

Pós-graduação em Ciências Biológicas (Zoologia) - IBRC

Caracterização molecular de cepas do Vírus da raiva (Lyssavirus; Rhabdoviridae) isoladas no Estado do Pará

A raiva é uma das zoonoses mais antigas e temidas pelo homem devido a seu desfecho fatal. Os cães ainda são considerados os principais responsáveis pela manutenção e transmissão da raiva para o homem. Porém, nos últimos anos os morcegos hematófagos têm ganhado destaque como potenciais transmissores de raiva para animais e humanos nas Américas. Recentemente, várias epidemias de raiva humana transmitida por morcegos hematófagos foram relatados no estado do Pará, o que mostra uma grande alteração no ambiente natural destes animais. A amplificação parcial do gene N pela técnica de RT-PCR foi aplicada em 62 amostras positivas para o Vírus da raiva, pela imunofluorescência direta e prova biológica. As seqüências nucleotídicas obtidas foram comparadas entre si e com outras amostras de vírus rábico isoladas no Brasil, utilizando os métodos de análise filogenética máxima verossimilhança e Bayesiano. Estas análises permitiram traçar o perfil epidemiológico molecular das variantes virais circulantes no estado do Pará, observando a emergência da transmissão de casos associados à variante antigênica 3 (VAg3), comumente encontrada em morcegos hematófagos Desmodus rotundus em detrimento dos casos relacionados à variante antigênica 2 (VAg2) associada a cães domésticos, bem como a identificação de três linhagens genéticas relacionadas a VAg3 e uma relacionada a VAg2 e uma possível nova variante isolada de morcego frugívoro Uroderma bilobatum.

Pesagem em movimento e caracterização do tráfego ferroviário com uso da técnica B-WIM

Nesta dissertação é apresentado o desenvolvimento de algoritmos para aplicação do método Bridge-Weigh In Motion (B-WIM) para a pesagem em movimento de trens e para a caracterização do tráfego ferroviário, permitindo-se obter informações sobre a velocidade de passagem dos trens, número e espaçamento entre eixos. Os sistemas B-WIM a partir de uma simples instrumentação permitem determinar as cargas por eixo de veículos em movimento, eliminando o efeito dinâmico. Foram implementados os algoritmos para a determinação dos valores referentes a geometria do trem e das cargas, que foi validado a partir de um exemplo teórico, onde se simulou a passagem de um trem de características conhecidas sobre a ponte e as cargas por eixos foram determinadas com 100% de exatidão. Além disso, foi feito um exemplo numérico em elementos finitos, de um viaduto em concreto armado para aplicação do método, onde foi feita a determinação das cargas por eixo para diferentes velocidades de passagem do trem. A fim de reduzir o tempo de processamento nas análises do exemplo numérico, foi desenvolvido um algoritmo para a geração de cargas nodais no modelo numérico que reduziram o tempo de processamento em até 96% quando comparado com a análise de múltiplos passos (“Multi-Step”), que simula automaticamente a passagem do trem sobre a estrutura. Finalmente, o método foi testado em um caso real a partir de monitorações realizadas em um viaduto de concreto armado da Estrada de Ferro Carajás. Apesar de não ter sido possível a determinação das cargas por eixo da locomotiva, foi possível medir precisamente o peso bruto total da locomotiva quando se utilizou o modelo constitutivo de Collins & Mitchell (1991) para o concreto.

Análise Espacial dos Anos Potenciais de Vida Perdidos por Causas Externas no Estado de São Paulo 2000 e 2010

Pós-graduação em Saúde Coletiva - FMB

Estimação Bayesiana dos parâmetros da distribuição Exponencial Generalizada Bivariada

Pós-graduação em Matematica Aplicada e Computacional - FCT

Análise da relação entre a distribuição espacial das morbidades por obesidade e hipertensão arterial para o estado de São Paulo, Brasil, de 2000 a 2010

A expansão da obesidade em diversos países do mundo na última década tem resultado no aumento da morbidade e mortalidade por hipertensão arterial e suas complicações. O objetivo deste trabalho é analisar a distribuição espacial da obesidade e hipertensão arterial no estado de São Paulo no período de 2000 a 2010, a partir de registros hospitalares e internação do Sistema de Informações Hospitalares do Sistema Único de Saúde (SIH - SUS). Foram utilizados coeficientes de prevalência das doenças em cada município suavizadas pelo método bayesiano empírico, permitindo uma visualização do padrão espacial dessas morbidades no Estado. Foi explorada a dependência espacial destes padrões verificando a autocorrelação entre os indicadores por meio do cálculo do Índice de Autocorrelação Espacial de Moran. Além disso, estudou-se a correlação positiva (Pearson) entre obesidade e hipertensão. Os dados e os mapas mostraram clusters de 87 municípios onde há maior e menor prevalência de hipertensão e obesidade no espaço com forte autocorrelação entre os municípios vizinhos. O coeficiente correlação de Pearson encontrado para esses municípios foi de 0,404 e sugere associação entre as morbidades. As técnicas de análise espacial mostraram-se úteis para o planejamento de ações de saúde pública.

Ecology, stable isotopes, and management of grassland songbirds at National Park Service properties on the Great Plains

Grassland ecosystems have been severely reduced and grassland bird populations have experienced consistent declines. National Park Service (NPS) properties on the Great Plains provide breeding habitat for grassland songbirds, though little is known about the quality of this habitat. A short-term study on songbirds at three NPS properties complemented current monitoring, providing an among park comparison addressing grassland bird productivity and fidelity relative to NPS property size. During 2008-2009, I assessed avian species richness, and estimated bird density and grassland songbird nest success. Bird species richness was greatest at small and medium sites, while number of nesting obligate species was greatest at the large site. Species-specific densities varied among sites, with few grassland obligates occurring at all three sites. Nest success estimates for grassland obligates were highest at the small site and lower at the large site. Another method to quantify habitat quality is assessment of breeding site fidelity. Current extrinsic markers used in monitoring site fidelity are inadequate for small birds; stable isotope analyses provide an alternative. I compared two techniques for assigning stable isotope tissue origin and measured grassland songbird site fidelity. My method of assigning origin provided site-specific variances of expected stable isotope values, an improvement over the most commonly used method. Fidelity tended to be higher at the large site, which may indicate a more robust breeding community of grassland birds. The small size of two of my sites precluded large sample sizes and made strong inferences difficult. To quantify how scientists cope with weak inference, I conducted a literature review. Strong inference was rarely observed, and most authors of weak-inference papers provided specific management recommendations. I suggest that adaptive management is an ideal method to resolve uncertainty from weak inference. Managers should consider my results within the context of regional and global management and the extent to which their unit might aide songbird conservation.

Testing allele homogeneity: the problem of nested hypotheses

Background: The evaluation of associations between genotypes and diseases in a case-control framework plays an important role in genetic epidemiology. This paper focuses on the evaluation of the homogeneity of both genotypic and allelic frequencies. The traditional test that is used to check allelic homogeneity is known to be valid only under Hardy-Weinberg equilibrium, a property that may not hold in practice. Results: We first describe the flaws of the traditional (chi-squared) tests for both allelic and genotypic homogeneity. Besides the known problem of the allelic procedure, we show that whenever these tests are used, an incoherence may arise: sometimes the genotypic homogeneity hypothesis is not rejected, but the allelic hypothesis is. As we argue, this is logically impossible. Some methods that were recently proposed implicitly rely on the idea that this does not happen. In an attempt to correct this incoherence, we describe an alternative frequentist approach that is appropriate even when Hardy-Weinberg equilibrium does not hold. It is then shown that the problem remains and is intrinsic of frequentist procedures. Finally, we introduce the Full Bayesian Significance Test to test both hypotheses and prove that the incoherence cannot happen with these new tests. To illustrate this, all five tests are applied to real and simulated datasets. Using the celebrated power analysis, we show that the Bayesian method is comparable to the frequentist one and has the advantage of being coherent. Conclusions: Contrary to more traditional approaches, the Full Bayesian Significance Test for association studies provides a simple, coherent and powerful tool for detecting associations.

Clusters of leprosy transmission and of late diagnosis in a highly endemic area in Brazil: focus on different spatial analysis approaches

Objective The Brazilian National Hansens Disease Control Program recently identified clusters with high disease transmission. Herein, we present different spatial analytical approaches to define highly vulnerable areas in one of these clusters. Method The study area included 373 municipalities in the four Brazilian states Maranha o, Para ', Tocantins and Piaui '. Spatial analysis was based on municipalities as the observation unit, considering the following disease indicators: (i) rate of new cases / 100 000 population, (ii) rate of cases < 15 years / 100 000 population, (iii) new cases with grade-2 disability / 100 000 population and (iv) proportion of new cases with grade-2 disabilities. We performed descriptive spatial analysis, local empirical Bayesian analysis and spatial scan statistic. Results A total of 254 (68.0%) municipalities were classified as hyperendemic (mean annual detection rates > 40 cases / 100 000 inhabitants). There was a concentration of municipalities with higher detection rates in Para ' and in the center of Maranha o. Spatial scan statistic identified 23 likely clusters of new leprosy case detection rates, most of them localized in these two states. These clusters included only 32% of the total population, but 55.4% of new leprosy cases. We also identified 16 significant clusters for the detection rate < 15 years and 11 likely clusters of new cases with grade-2. Several clusters of new cases with grade-2 / population overlap with those of new cases detection and detection of children < 15 years of age. The proportion of new cases with grade-2 did not reveal any significant clusters. Conclusions Several municipality clusters for high leprosy transmission and late diagnosis were identified in an endemic area using different statistical approaches. Spatial scan statistic is adequate to validate and confirm high-risk leprosy areas for transmission and late diagnosis, identified using descriptive spatial analysis and using local empirical Bayesian method. National and State leprosy control programs urgently need to intensify control actions in these highly vulnerable municipalities.

Automatische objektive Identifikation und Berechnung der Eigenschaften von Rossbywellenzügen

In dieser Arbeit wird eine Methode entwickelt, die Rossbywellenzüge automatisch identifiziert und und deren Eigenschaften quantifiziert. Mit dieser Methode wird der Wellenzug als eine Einheit in Raum und Zeit interpretiert. Einheit im Raum heißt, dass nicht die einzelnen Tröge und Rücken eines Wellenzugs betrachtet werden, sondern deren Einhüllende. Einheit in der Zeit bedeutet, dass der Wellenzug nicht nur zu einem Zeitpunkt betrachtet wird, sondern über seine gesamte Lebensdauer hinweg. Um den Wellenzug als räumliche und zeitliche Einheit zu erhalten, werden die Einhüllenden der Wellenzüge in Längengrad-Zeit Diagrammen, sogenannten Hovmöllerdiagrammen, betrachtet. Dort werden zusammenhängende Regionen als Objekte, die jeweils einen Wellenzug repräsentieren, identifiziert. Deren Eigenschaften werden dann automatisch berechnet. Diese Eigenschaften können nun direkt dem zugrunde liegenden Rossbywellenzug zugeordnet werden.rnDie neue Methode wird in zwei verschiedenen Szenarien angewendet: erstens zur Beurteilung der Vorhersagequalität eines einzelnen Rossbywellenzugs und zweitens für die klimatologische Betrachtung von Rossbywellenzügen im ERA-40 Reanalysedatensatz. Sie wurde weiterhin mit bisher verwendeten Methoden zur Identifikation und Quantifizierung von Rossbywellenzügen verglichen.rnDie Untersuchung der Vorhersagequalität ergab, dass in dem betrachteten Fall die Übereinstimmung der Vorhersage mit der Analyse des Wellenzugs gering war, sofern das Modell initialisiert wurde, bevor der Rossbywellenzug eingesetzt hatte. Im Gegensatz dazu nahm die Vorhersagequalität deutlich zu, wenn der Wellenzug bereits in den Vorhersagedaten enthalten war. Dies deutet darauf hin, dass es in dem vorliegenden Fall problematisch ist, mit dem Modell den Auslösemechanismus korrekt voherzusagen. Für die weitere Untersuchung der Vorhersagequalität wurde eine spezielle Art der Darstellung der Daten verwendet, mit deren Hilfe deutlich wurde, dass das verwendete Modell in der Lage ist, diesen Wellenzug ungefähr sechs Tage im Voraus vorherzusagen. Diese Zeitspanne ist deutlich kürzer als die Lebensdauer des Wellenzugs, die etwa 10 Tage beträgt.rnIm Rahmen der klimatologischen Studie ergab sich eine positive Korrelation zwischen der Lebensdauer eines Rossbywellenzugs und des Bereichs den dieser Wellenzug während seiner gesamten Existenz in zonaler Richtung überstreicht. Für Wellenzüge mit einer kurzen Lebensdauer ergab sich eine ebenfalls positive Korrelation zwischen der mittleren Amplitude und der Dauer des Wellenzugs. Für eine längere Lebensdauer geht diese Korrelation aber in eine Sättigung über und die mittlere Amplitude steigt nicht mehr weiter an. Als eine mögliche Erklärung für dieses Verhalten wird angeführt, dass eine gewisse Stärke der Amplitude benötigt wird um stromabwärtige Entwicklung zu erhalten aber zu große Amplituden im Allgemeinen zum Brechen der Welle führen. Das Brechen leitet den letzten Abschnitt im Lebenszyklus eines Rossbywellenzuges ein, welcher im Anschluss meist zerfällt. Ein weiteres Ergebnis der klimatologischen Untersuchung ist das Auffinden bevorzugter Regionen der Entstehung und des Abklingens von Rossbywellenzügen. Diese Regionen unterscheiden sich erheblich für Rossbywellenzüge unterschiedlicher minimaler Lebensdauer. Langlebige Rossbywellenzüge entstehen demnach hauptsächlich über Ostasien und dem Westpazifik und vergehen dann über Europa.rnSchließlich wurde die entwickelte Methode in einen systematischen Vergleich anderer Methoden zur Identifikation und Quantifizierung von Rossbywellenzügen eingereiht. Die betrachteten Methoden beinhalten verschiedene Trog-und-Rücken Hovmöllerdiagramme des Meridionalwindes, Methoden die Rossbywellenzüge als eine Einheit identifizieren und Methoden die den Beitrag verschiedener physikalischer Aspekte zu der Entwicklung von Rossbywellenenzügen quantifizieren. Der Vergleich macht deutlich, dass jede Methode ihre individuellen Stärken und Schwächen hat. Dies bedeutet insbesondere, dass die Eignung der Methode von dem Stadium des Lebenszyklus, in dem sich der Rossbywellenzug befindet und dem Fokus, den man bei der Betrachtung hat, abhängt. Ideal ist eine Kombination mehrerer Methoden, da dies ein vollständigeres Bild eines Rossbywellenzuges ergibt als einzelne Methoden es zu liefern vermögen. Obwohl alle Methoden für die Anwendungen, für die sie jeweils konzipiert wurden, geeignet sind, ergeben sich bei der Diagnose der Rossbywellenzüge beträchtliche Unterschiede. Letztendlich stellt sich heraus, dass sogar die Defintion eines Rossbywellenzugs bis zu einem gewissen Grad von der zu seiner Identifizierung verwendeten Methode abhängt.

Estimating the Number of Essential Genes in a Genome by Random Transposon Mutagenesis

We describe a Bayesian method for estimating the number of essential genes in a genome, on the basis of data on viable mutants for which a single transposon was inserted after a random TA site in a genome,potentially disrupting a gene. The prior distribution for the number of essential genes was taken to be uniform. A Gibbs sampler was used to estimate the posterior distribution. The method is illustrated with simulated data. Further simulations were used to study the performance of the procedure.

IDENTIFICATION OF GENES CONTROLLING BIOLOGICAL PROCESSES AND PATHWAYS THROUGH STATISTICAL ANALYSIS AND NETWORK RECONSTRUCTION

The developmental processes and functions of an organism are controlled by the genes and the proteins that are derived from these genes. The identification of key genes and the reconstruction of gene networks can provide a model to help us understand the regulatory mechanisms for the initiation and progression of biological processes or functional abnormalities (e.g. diseases) in living organisms. In this dissertation, I have developed statistical methods to identify the genes and transcription factors (TFs) involved in biological processes, constructed their regulatory networks, and also evaluated some existing association methods to find robust methods for coexpression analyses. Two kinds of data sets were used for this work: genotype data and gene expression microarray data. On the basis of these data sets, this dissertation has two major parts, together forming six chapters. The first part deals with developing association methods for rare variants using genotype data (chapter 4 and 5). The second part deals with developing and/or evaluating statistical methods to identify genes and TFs involved in biological processes, and construction of their regulatory networks using gene expression data (chapter 2, 3, and 6). For the first part, I have developed two methods to find the groupwise association of rare variants with given diseases or traits. The first method is based on kernel machine learning and can be applied to both quantitative as well as qualitative traits. Simulation results showed that the proposed method has improved power over the existing weighted sum method (WS) in most settings. The second method uses multiple phenotypes to select a few top significant genes. It then finds the association of each gene with each phenotype while controlling the population stratification by adjusting the data for ancestry using principal components. This method was applied to GAW 17 data and was able to find several disease risk genes. For the second part, I have worked on three problems. First problem involved evaluation of eight gene association methods. A very comprehensive comparison of these methods with further analysis clearly demonstrates the distinct and common performance of these eight gene association methods. For the second problem, an algorithm named the bottom-up graphical Gaussian model was developed to identify the TFs that regulate pathway genes and reconstruct their hierarchical regulatory networks. This algorithm has produced very significant results and it is the first report to produce such hierarchical networks for these pathways. The third problem dealt with developing another algorithm called the top-down graphical Gaussian model that identifies the network governed by a specific TF. The network produced by the algorithm is proven to be of very high accuracy.

Toward an early diagnosis of lung cancer: an autoantibody signature for squamous cell lung carcinoma

Serum-based diagnosis offers the prospect of early lung carcinoma detection and of differentiation between benign and malignant nodules identified by CT. One major challenge toward a future blood-based diagnostic consists in showing that seroreactivity patterns allow for discriminating lung cancer patients not only from normal controls but also from patients with non-tumor lung pathologies. We addressed this question for squamous cell lung cancer, one of the most common lung tumor types. Using a panel of 82 phage-peptide clones, which express potential autoantigens, we performed serological spot assay. We screened 108 sera, including 39 sera from squamous cell lung cancer patients, 29 sera from patients with other non-tumor lung pathologies, and 40 sera from volunteers without known disease. To classify the serum groups, we employed the standard Naïve Bayesian method combined with a subset selection approach. We were able to separate squamous cell lung carcinoma and normal sera with an accuracy of 93%. Low-grade squamous cell lung carcinoma were separated from normal sera with an accuracy of 92.9%. We were able to distinguish squamous cell lung carcinoma from non-tumor lung pathologies with an accuracy of 83%. Three phage-peptide clones with sequence homology to ROCK1, PRKCB1 and KIAA0376 reacted with more than 15% of the cancer sera, but neither with normal nor with non-tumor lung pathology sera. Our study demonstrates that seroreactivity profiles combined with statistical classification methods have great potential for discriminating patients with squamous cell lung carcinoma not only from normal controls but also from patients with non-tumor lung pathologies.