Bayesian phylogenetic modelling of lateral gene transfers

PhD ThesisPhylogenetic trees represent the evolutionary relationships between a set of species. Inferring these trees from data is particularly challenging sometimes since the transfer of genetic material can occur not only from parents to their o spring but also between organisms via lateral gene transfers (LGTs). Thus, the presence of LGTs means that genes in a genome can each have di erent evolutionary histories, represented by di erent gene trees. A few statistical approaches have been introduced to explore non-vertical evolution through collections of Markov-dependent gene trees. In 2005 Suchard described a Bayesian hierarchical model for joint inference of gene trees and an underlying species tree, where a layer in the model linked gene trees to the species tree via a sequence of unknown lateral gene transfers. In his model LGT was modeled via a random walk in the tree space derived from the subtree prune and regraft (SPR) operator on unrooted trees. However, the use of SPR moves to represent LGT in an unrooted tree is problematic, since the transference of DNA between two organisms implies the contemporaneity of both organisms and therefore it can allow unrealistic LGTs. This thesis describes a related hierarchical Bayesian phylogenetic model for reconstructing phylogenetic trees which imposes a temporal constraint on LGTs, namely that they can only occur between species which exist concurrently. This is achieved by taking into account possible time orderings of divergence events in trees, without explicitly modelling divergence times. An extended version of the SPR operator is introduced as a more adequate mechanism to represent the LGT e ect in a tree. The extended SPR operation respects the time ordering. It additionaly di ers from regular SPR as it maintains a 1-to-1 correspondence between points on the species tree and points on each gene tree. Each point on a gene tree represents the existence of a population containing that gene at some point in time. Hierarchical phylogenetic models were used in the reconstruction of each gene tree from its corresponding gene alignment, enabling the pooling of information across genes. In addition to Suchard's approach, we assume variation in the rate of evolution between di erent sites. The species tree is assumed to be xed. A Markov Chain Monte Carlo (MCMC) algorithm was developed to t the model in a Bayesian framework. A novel MCMC proposal mechanism for jointly proposing the gene tree topology and branch lengths, LGT distance and LGT history has been developed as well as a novel graphical tool to represent LGT history, the LGT Biplot. Our model was applied to simulated and experimental datasets. More speci cally we analysed LGT/reassortment presence in the evolution of 2009 Swine-Origin In uenza Type A virus. Future improvements of our model and algorithm should include joint inference of the species tree, improving the computational e ciency of the MCMC algorithm and better consideration of other factors that can cause discordance of gene trees and species trees such as gene loss

Epidemic history of hepatitis C virus genotypes and subtypes in Portugal

Any successful strategy to prevent and control HCV infection requires an understanding of the epidemic behaviour among the different genotypes. Here, we performed the first characterization of the epidemic history and transmission dynamics of HCV subtypes in Portugal. Direct sequencing of NS5B was performed on 230 direct-acting antiviral drugs (DAA)-treatment naïve patients in Lisbon. Phylogenetic analysis was used for subtyping and transmission cluster identification. Bayesian methods were used to reconstruct the epidemic history of HCV subtypes. Sequences were analysed for resistance-associated substitutions (RAS). The majority of strains were HCV-GT1 (62.6%), GT3 (18.3%, all subtype 3a) and GT4 (16.1%). Among GT1, the most frequent were subtypes 1a (75.5%) and 1b (24.5%). Polyphyletic patterns were found in all but 12 lineages suggesting multiple introductions of the different subtypes in this population. Five distinct epidemics were identified. The first significant HCV epidemic in Portugal occurred between 1930s and 1960s, was caused almost exclusively by GT1b and was likely associated with blood transfusions. Rapid expansion of GT3a occurred in the 1960s and GT1a in the 1980s, associated with intravenous drug use. The most recent epidemics were caused by GT4a and GT4d and seem to be associated with the resurgence of opioid use. The C316N substitution was found in 31.4% of GT1b-patients. Close surveillance of patients bearing this mutation and undergoing dasabuvir-based regimens will be important to determine its impact on treatment outcome.publishersversionpublishe

Transforming Growth Factor Beta 2 and Heme Oxygenase 1 Genes Are Risk Factors for the Cerebral Malaria Syndrome in Angolan Children

BACKGROUND: Cerebral malaria (CM) represents a severe outcome of the Plasmodium falciparum infection. Recent genetic studies have correlated human genes with severe malaria susceptibility, but there is little data on genetic variants that increase the risk of developing specific malaria clinical complications. Nevertheless, susceptibility to experimental CM in the mouse has been linked to host genes including Transforming Growth Factor Beta 2 (TGFB2) and Heme oxygenase-1 (HMOX1). Here, we tested whether those genes were governing the risk of progressing to CM in patients with severe malaria syndromes. METHODOLOGY/PRINCIPAL FINDINGS: We report that the clinical outcome of P. falciparum infection in a cohort of Angolan children (n = 430) correlated with nine TGFB2 SNPs that modify the risk of progression to CM as compared to other severe forms of malaria. This genetic effect was explained by two haplotypes harboring the CM-associated SNPs (Pcorrec. = 0.035 and 0.036). In addition, one HMOX1 haplotype composed of five CM-associated SNPs increased the risk of developing the CM syndrome (Pcorrec. = 0.002) and was under-transmitted to children with uncomplicated malaria (P = 0.036). Notably, the HMOX1-associated haplotype conferred increased HMOX1 mRNA expression in peripheral blood cells of CM patients (P = 0.012). CONCLUSIONS/SIGNIFICANCE: These results represent the first report on CM genetic risk factors in Angolan children and suggest the novel hypothesis that genetic variants of the TGFB2 and HMOX1 genes may contribute to confer a specific risk of developing the CM syndrome in patients with severe P. falciparum malaria. This work may provide motivation for future studies aiming to replicate our findings in larger populations and to confirm a role for these genes in determining the clinical course of malaria

Transforming Growth Factor Beta 2 and Heme Oxygenase 1 Genes Are Risk Factors for the Cerebral Malaria Syndrome in Angolan Children

BACKGROUND: Cerebral malaria (CM) represents a severe outcome of the Plasmodium falciparum infection. Recent genetic studies have correlated human genes with severe malaria susceptibility, but there is little data on genetic variants that increase the risk of developing specific malaria clinical complications. Nevertheless, susceptibility to experimental CM in the mouse has been linked to host genes including Transforming Growth Factor Beta 2 (TGFB2) and Heme oxygenase-1 (HMOX1). Here, we tested whether those genes were governing the risk of progressing to CM in patients with severe malaria syndromes. METHODOLOGY/PRINCIPAL FINDINGS: We report that the clinical outcome of P. falciparum infection in a cohort of Angolan children (n = 430) correlated with nine TGFB2 SNPs that modify the risk of progression to CM as compared to other severe forms of malaria. This genetic effect was explained by two haplotypes harboring the CM-associated SNPs (Pcorrec. = 0.035 and 0.036). In addition, one HMOX1 haplotype composed of five CM-associated SNPs increased the risk of developing the CM syndrome (Pcorrec. = 0.002) and was under-transmitted to children with uncomplicated malaria (P = 0.036). Notably, the HMOX1-associated haplotype conferred increased HMOX1 mRNA expression in peripheral blood cells of CM patients (P = 0.012). CONCLUSIONS/SIGNIFICANCE: These results represent the first report on CM genetic risk factors in Angolan children and suggest the novel hypothesis that genetic variants of the TGFB2 and HMOX1 genes may contribute to confer a specific risk of developing the CM syndrome in patients with severe P. falciparum malaria. This work may provide motivation for future studies aiming to replicate our findings in larger populations and to confirm a role for these genes in determining the clinical course of malaria

Epidemic history of hepatitis C virus genotypes and subtypes in Portugal

Any successful strategy to prevent and control HCV infection requires an understanding of the epidemic behaviour among the different genotypes. Here, we performed the first characterization of the epidemic history and transmission dynamics of HCV subtypes in Portugal. Direct sequencing of NS5B was performed on 230 direct-acting antiviral drugs (DAA)-treatment naïve patients in Lisbon. Phylogenetic analysis was used for subtyping and transmission cluster identification. Bayesian methods were used to reconstruct the epidemic history of HCV subtypes. Sequences were analysed for resistance-associated substitutions (RAS). The majority of strains were HCV-GT1 (62.6%), GT3 (18.3%, all subtype 3a) and GT4 (16.1%). Among GT1, the most frequent were subtypes 1a (75.5%) and 1b (24.5%). Polyphyletic patterns were found in all but 12 lineages suggesting multiple introductions of the different subtypes in this population. Five distinct epidemics were identified. The first significant HCV epidemic in Portugal occurred between 1930s and 1960s, was caused almost exclusively by GT1b and was likely associated with blood transfusions. Rapid expansion of GT3a occurred in the 1960s and GT1a in the 1980s, associated with intravenous drug use. The most recent epidemics were caused by GT4a and GT4d and seem to be associated with the resurgence of opioid use. The C316N substitution was found in 31.4% of GT1b-patients. Close surveillance of patients bearing this mutation and undergoing dasabuvir-based regimens will be important to determine its impact on treatment outcome.The authors would like to thank all the patients for their participation. C.P. is supported by the Portuguese ‘Fundação para a Ciência e Tecnologia’ (FCT) (grant number SFRH/BPD/77448/2011, part of the EDCTP2 program supported by the European Union). R.M. was supported by a Ph.D. scholarship by the Portuguese ‘Fundação para a Ciência e Tecnologia’ (FCT). V.B. is supported by the Miguel Servet program run by the ‘Fondo de Investigación Sanitaria’ (ISCIII) (grant number CP13/00098). I.M. is supported by the Portuguese ‘Fundação para a Ciência e Tecnologia’ (FCT) (grant number SFRH/BD/131062/2017).S

Epidemic history of hepatitis C virus genotypes and subtypes in Portugal

Any successful strategy to prevent and control HCV infection requires an understanding of the epidemic behaviour among the different genotypes. Here, we performed the first characterization of the epidemic history and transmission dynamics of HCV subtypes in Portugal. Direct sequencing of NS5B was performed on 230 direct-acting antiviral drugs (DAA)-treatment naïve patients in Lisbon. Phylogenetic analysis was used for subtyping and transmission cluster identification. Bayesian methods were used to reconstruct the epidemic history of HCV subtypes. Sequences were analysed for resistance-associated substitutions (RAS). The majority of strains were HCV-GT1 (62.6%), GT3 (18.3%, all subtype 3a) and GT4 (16.1%). Among GT1, the most frequent were subtypes 1a (75.5%) and 1b (24.5%). Polyphyletic patterns were found in all but 12 lineages suggesting multiple introductions of the different subtypes in this population. Five distinct epidemics were identified. The first significant HCV epidemic in Portugal occurred between 1930s and 1960s, was caused almost exclusively by GT1b and was likely associated with blood transfusions. Rapid expansion of GT3a occurred in the 1960s and GT1a in the 1980s, associated with intravenous drug use. The most recent epidemics were caused by GT4a and GT4d and seem to be associated with the resurgence of opioid use. The C316N substitution was found in 31.4% of GT1b-patients. Close surveillance of patients bearing this mutation and undergoing dasabuvir-based regimens will be important to determine its impact on treatment outcome.The authors would like to thank all the patients for their participation. C.P. is supported by the Portuguese ‘Fundação para a Ciência e Tecnologia’ (FCT) (grant number SFRH/BPD/77448/2011, part of the EDCTP2 program supported by the European Union). R.M. was supported by a Ph.D. scholarship by the Portuguese ‘Fundação para a Ciência e Tecnologia’ (FCT). V.B. is supported by the Miguel Servet program run by the ‘Fondo de Investigación Sanitaria’ (ISCIII) (grant number CP13/00098). I.M. is supported by the Portuguese ‘Fundação para a Ciência e Tecnologia’ (FCT) (grant number SFRH/BD/131062/2017).info:eu-repo/semantics/publishedVersio