Browsing by Subject "Sequence alignment"

Now showing 1 - 18 of 18
  • Thumbnail Image
    Item
    Open Access
    Avoidance of Protein Fold Disruption in Natural Virus Recombinants
    (Public Library of Science, 2007) Lefeuvre, Pierre; Lett, Jean-Michel; Reynaud, Bernard; Martin, Darren P
    Author Summary The exchange of genetic material between different virus species, called inter-species recombination, has the potential to generate, within a single genome replication cycle, an almost unimaginable number of genetically distinct virus strains, including many that might cause deadly new human, animal, or plant diseases. Many fear that inter-species recombination could provide viruses with quick access to evolutionary innovations such as broader host ranges, altered tissue tropisms, or increased severities. However, mounting evidence suggests that recombination is not an unconstrained process and that most inter-species recombinants that occur in nature are probably defective. It is suspected that networks of coevolved interactions between different parts of virus genomes and their encoded proteins must be kept intact for newly formed inter-species recombinants to have any chance of out-competing their parents. One category of coevolved interaction is that between contacting amino acids within the 3-D structures of folded proteins. Here we examine the distributions of recombination events across the genomes of a group of rampantly recombining plant viruses and find very good evidence that this class of interaction tends to be preserved amongst recombinant sequences sampled from nature. This indicates that selection against misfolded proteins strongly influences the survival of natural recombinants.
  • Thumbnail Image
    Item
    Open Access
    Bat accelerated regions identify a bat forelimb specific enhancer in the HoxD locus
    (Public Library of Science, 2016) Booker, Betty M; Friedrich, Tara; Mason, Mandy K; VanderMeer, Julia E; Zhao, Jingjing; Eckalbar, Walter L; Logan, Malcolm; Illing, Nicola; Pollard, Katherine S; Ahituv, Nadav
    Author Summary: The limb is a classic example of vertebrate homology and is represented by a large range of morphological structures such as fins, legs and wings. The evolution of these structures could be driven by alterations in gene regulatory elements that have critical roles during development. To identify elements that may contribute to bat wing development, we characterized sequences that are conserved between vertebrates, but changed significantly in the bat lineage. We then overlapped these sequences with predicted developing limb enhancers as determined by ChIP-seq, finding 166 bat accelerated sequences (BARs). Five BARs that were tested for enhancer activity in mice all drove expression in the limb. Testing the mouse orthologous sequence showed that three had differences in their limb enhancer activity as compared to the bat sequence. Of these, BAR116 was of particular interest as it is located near the HoxD locus, an essential gene complex required for proper spatiotemporal patterning of the developing limb. The bat BAR116 sequence drove robust forelimb expression but the mouse BAR116 sequence did not show enhancer activity. These experiments correspond to analyses of HoxD gene expressions in developing bat limbs, which had strong forelimb versus weak hindlimb expression for Hoxd10 - 11 . Combined, our studies highlight specific genomic regions that could be important in shaping the morphological differences that led to the development of the bat wing.
  • Thumbnail Image
    Item
    Open Access
    Comparative genomic and phylogenetic approaches to characterize the role of genetic recombination in mycobacterial evolution
    (Public Library of Science, 2012) Smith, Silvia E; Showers-Corneli, Patrice; Dardenne, Caitlin N; Harpending, Henry H; Martin, Darren P; Beiko, Robert G
    The genus Mycobacterium encompasses over one hundred named species of environmental and pathogenic organisms, including the causative agents of devastating human diseases such as tuberculosis and leprosy. The success of these human pathogens is due in part to their ability to rapidly adapt to their changing environment and host. Recombination is the fastest way for bacterial genomes to acquire genetic material, but conflicting results about the extent of recombination in the genus Mycobacterium have been reported. We examined a data set comprising 18 distinct strains from 13 named species for evidence of recombination. Genomic regions common to all strains (accounting for 10% to 22% of the full genomes of all examined species) were aligned and concatenated in the chromosomal order of one mycobacterial reference species. The concatenated sequence was screened for evidence of recombination using a variety of statistical methods, with each proposed event evaluated by comparing maximum-likelihood phylogenies of the recombinant section with the non-recombinant portion of the dataset. Incongruent phylogenies were identified by comparing the site-wise log-likelihoods of each tree using multiple tests. We also used a phylogenomic approach to identify genes that may have been acquired through horizontal transfer from non-mycobacterial sources. The most frequent associated lineages (and potential gene transfer partners) in the Mycobacterium lineage-restricted gene trees are other members of suborder Corynebacterinae, but more-distant partners were identified as well. In two examined cases of potentially frequent and habitat-directed transfer ( M. abscessus to Segniliparus and M. smegmatis to Streptomyces ), observed sequence distances were small and consistent with a hypothesis of transfer, while in a third case ( M. vanbaalenii to Streptomyces ) distances were larger. The analyses described here indicate that whereas evidence of recombination in core regions within the genus is relatively sparse, the acquisition of genes from non-mycobacterial lineages is a significant feature of mycobacterial evolution.
  • Thumbnail Image
    Item
    Open Access
    Differences in genotype and virulence among four multidrug-resistant Streptococcus pneumoniae isolates belonging to the PMEN1 clone
    (Public Library of Science, 2011) Hiller, N Luisa; Eutsey, Rory A; Powell, Evan; Earl, Joshua P; Janto, Benjamin; Martin, Darren P; Dawid, Suzanne; Ahmed, Azad; Longwell, Mark J; Dahlgren, Margaret E
    We report on the comparative genomics and characterization of the virulence phenotypes of four S. pneumoniae strains that belong to the multidrug resistant clone PMEN1 (Spain 23F ST81). Strains SV35-T23 and SV36-T3 were recovered in 1996 from the nasopharynx of patients at an AIDS hospice in New York. Strain SV36-T3 expressed capsule type 3 which is unusual for this clone and represents the product of an in vivo capsular switch event. A third PMEN1 isolate - PN4595-T23 - was recovered in 1996 from the nasopharynx of a child attending day care in Portugal, and a fourth strain - ATCC700669 - was originally isolated from a patient with pneumococcal disease in Spain in 1984. We compared the genomes among four PMEN1 strains and 47 previously sequenced pneumococcal isolates for gene possession differences and allelic variations within core genes. In contrast to the 47 strains - representing a variety of clonal types - the four PMEN1 strains grouped closely together, demonstrating high genomic conservation within this lineage relative to the rest of the species. In the four PMEN1 strains allelic and gene possession differences were clustered into 18 genomic regions including the capsule, the blp bacteriocins, erythromycin resistance, the MM1-2008 prophage and multiple cell wall anchored proteins. In spite of their genomic similarity, the high resolution chinchilla model was able to detect variations in virulence properties of the PMEN1 strains highlighting how small genic or allelic variation can lead to significant changes in pathogenicity and making this set of strains ideal for the identification of novel virulence determinants.
  • Thumbnail Image
    Item
    Open Access
    Evolutionary time-scale of the begomoviruses: evidence from integrated sequences in the Nicotiana genome
    (Public Library of Science, 2011) Lefeuvre, Pierre; Harkins, Gordon W; Lett, Jean-Michel; Briddon, Rob W; Chase, Mark W; Moury, Benoit; Martin, Darren P
    Despite having single stranded DNA genomes that are replicated by host DNA polymerases, viruses in the family Geminiviridae are apparently evolving as rapidly as some RNA viruses. The observed substitution rates of geminiviruses in the genera Begomovirus and Mastrevirus are so high that the entire family could conceivably have originated less than a million years ago (MYA). However, the existence of geminivirus related DNA (GRD) integrated within the genomes of various Nicotiana species suggests that the geminiviruses probably originated >10 MYA. Some have even suggested that a distinct New-World (NW) lineage of begomoviruses may have arisen following the separation by continental drift of African and American proto-begomoviruses ∼110 MYA. We evaluate these various geminivirus origin hypotheses using Bayesian coalescent-based approaches to date firstly the Nicotiana GRD integration events, and then the divergence of the NW and Old-World (OW) begomoviruses. Besides rejecting the possibility of a<2 MYA OW-NW begomovirus split, we could also discount that it may have occurred concomitantly with the breakup of Gondwanaland 110 MYA. Although we could only confidently narrow the date of the split down to between 2 and 80 MYA, the most plausible (and best supported) date for the split is between 20 and 30 MYA - a time when global cooling ended the dispersal of temperate species between Asia and North America via the Beringian land bridge.
  • Thumbnail Image
    Item
    Open Access
    Frequent toggling between alternative amino acids is driven by selection in HIV-1
    (Public Library of Science, 2008) Delport, Wayne; Scheffler, Konrad; Seoighe, Cathal
    Author Summary Viruses, such as HIV, are able to evade host immune responses through escape mutations, yet sometimes they do so at a cost. This cost is the reduction in the ability of the virus to replicate, and thus selective pressure exists for a virus to revert to its original state in the absence of the host immune response that caused the initial escape mutation. This pattern of escape and reversion typically occurs when viruses are transmitted between individuals with different immune responses. We develop a phylogenetic model of immune escape and reversion and provide evidence that it outperforms existing models for the detection of selective pressure associated with host immune responses. Finally, we demonstrate that amino acid toggling is a pervasive process in HIV-1 evolution, such that many of the positions in the virus that evolve rapidly, under the influence of positive Darwinian selection, nonetheless display quite low sequence diversity. This highlights the limitations of HIV-1 evolution, and sites such as these are potentially good targets for HIV-1 vaccines.
  • Thumbnail Image
    Item
    Open Access
    Generation of genic diversity among Streptococcus pneumoniae strains via horizontal gene transfer during a chronic polyclonal pediatric infection
    (Public Library of Science, 2010) Hiller, N Luisa; Ahmed, Azad; Powell, Evan; Martin, Darren P; Eutsey, Rory; Earl, Josh; Janto, Benjamin; Boissy, Robert J; Hogg, Justin; Barbadora, Karen
    Although there is tremendous interest in understanding the evolutionary roles of horizontal gene transfer (HGT) processes that occur during chronic polyclonal infections, to date there have been few studies that directly address this topic. We have characterized multiple HGT events that most likely occurred during polyclonal infection among nasopharyngeal strains of Streptococcus pneumoniae recovered from a child suffering from chronic upper respiratory and middle-ear infections. Whole genome sequencing and comparative genomics were performed on six isolates collected during symptomatic episodes over a period of seven months. From these comparisons we determined that five of the isolates were genetically highly similar and likely represented a dominant lineage. We analyzed all genic and allelic differences among all six isolates and found that all differences tended to occur within contiguous genomic blocks, suggestive of strain evolution by homologous recombination. From these analyses we identified three strains (two of which were recovered on two different occasions) that appear to have been derived sequentially, one from the next, each by multiple recombination events. We also identified a fourth strain that contains many of the genomic segments that differentiate the three highly related strains from one another, and have hypothesized that this fourth strain may have served as a donor multiple times in the evolution of the dominant strain line. The variations among the parent, daughter, and grand-daughter recombinant strains collectively cover greater than seven percent of the genome and are grouped into 23 chromosomal clusters. While capturing in vivo HGT, these data support the distributed genome hypothesis and suggest that a single competence event in pneumococci can result in the replacement of DNA at multiple non-adjacent loci.
  • Thumbnail Image
    Item
    Open Access
    Genetic variability among complete human respiratory syncytial virus subgroup A genomes: bridging molecular evolutionary dynamics and epidemiology
    (Public Library of Science, 2012) Tan, Lydia; Lemey, Philippe; Houspie, Lieselot; Viveen, Marco C; Jansen, Nicolaas J G; Loon, Anton M van; Wiertz, Emmanuel; Bleek, Grada M van; Martin, Darren P; Coenjaerts, Frank E
    Human respiratory syncytial virus (RSV) is an important cause of severe lower respiratory tract infections in infants and the elderly. In the vast majority of cases, however, RSV infections run mild and symptoms resemble those of a common cold. The immunological, clinical, and epidemiological profile of severe RSV infections suggests a disease caused by a virus with typical seasonal transmission behavior, lacking clear-cut virulence factors, but instead causing disease by modifying the host's immune response in a way that stimulates pathogenesis. Yet, the interplay between RSV-evoked immune responses and epidemic behavior, and how this affects the genomic evolutionary dynamics of the virus, remains poorly understood. Here, we present a comprehensive collection of 33 novel RSV subgroup A genomes from strains sampled over the last decade, and provide the first measurement of RSV-A genomic diversity through time in a phylodynamic framework. In addition, we map amino acid substitutions per protein to determine mutational hotspots in specific domains. Using Bayesian genealogical inference, we estimated the genomic evolutionary rate to be 6.47×10 −4 (credible interval: 5.56×10 −4 , 7.38×10 −4 ) substitutions/site/year, considerably slower than previous estimates based on G gene sequences only. The G gene is however marked by elevated substitution rates compared to other RSV genes, which can be attributed to relaxed selective constraints. In line with this, site-specific selection analyses identify the G gene as the major target of diversifying selection. Importantly, statistical analysis demonstrates that the immune driven positive selection does not leave a measurable imprint on the genome phylogeny, implying that RSV lineage replacement mainly follows nonselective epidemiological processes. The roughly 50 years of RSV-A genomic evolution are characterized by a constant population size through time and general co-circulation of lineages over many epidemic seasons - a conclusion that might be taken into account when developing future therapeutic and preventive strategies.
  • Thumbnail Image
    Item
    Open Access
    Genetic variation and recombination of RdRp and HSP 70h genes of Citrus tristeza virus isolates from orange trees showing symptoms of citrus sudden death disease
    (BioMed Central Ltd, 2008) Gomes, Clarissa; Nagata, Tatsuya; de Jesus, Waldir; Neto, Carlos; Pappas, Georgios; Martin, Darren
    BACKGROUND: Citrus sudden death (CSD), a disease that rapidly kills orange trees, is an emerging threat to the Brazilian citrus industry. Although the causal agent of CSD has not been definitively determined, based on the disease's distribution and symptomatology it is suspected that the agent may be a new strain of Citrus tristeza virus (CTV). CTV genetic variation was therefore assessed in two Brazilian orange trees displaying CSD symptoms and a third with more conventional CTV symptoms. RESULTS: A total of 286 RNA-dependent-RNA polymerase (RdRp) and 284 heat shock protein 70 homolog (HSP70h) gene fragments were determined for CTV variants infecting the three trees. It was discovered that, despite differences in symptomatology, the trees were all apparently coinfected with similar populations of divergent CTV variants. While mixed CTV infections are common, the genetic distance between the most divergent population members observed (24.1% for RdRp and 11.0% for HSP70h) was far greater than that in previously described mixed infections. Recombinants of five distinct RdRp lineages and three distinct HSP70h lineages were easily detectable but respectively accounted for only 5.9 and 11.9% of the RdRp and HSP70h gene fragments analysed and there was no evidence of an association between particular recombinant mosaics and CSD. Also, comparisons of CTV population structures indicated that the two most similar CTV populations were those of one of the trees with CSD and the tree without CSD. CONCLUSION: We suggest that if CTV is the causal agent of CSD, it is most likely a subtle feature of population structures within mixed infections and not merely the presence (or absence) of a single CTV variant within these populations that triggers the disease.
  • Thumbnail Image
    Item
    Open Access
    Genotype at the P554L variant of the hexose-6 phosphate dehydrogenase gene is associated with carotid intima-medial thickness
    (Public Library of Science, 2011) Rahman, Thahira J; Walker, Elizabeth A; Mayosi, Bongani M; Hall, Darroch H; Avery, Peter J; Connell, John M C; Watkins, Hugh; Stewart, Paul M; Keavney, Bernard
    Objective The combined thickness of the intima and media of the carotid artery (carotid intima-medial thickness, CIMT) is associated with cardiovascular disease and stroke. Previous studies indicate that carotid intima-medial thickness is a significantly heritable phenotype, but the responsible genes are largely unknown. Hexose-6 phosphate dehydrogenase (H6PDH) is a microsomal enzyme whose activity regulates corticosteroid metabolism in the liver and adipose tissue; variability in measures of corticosteroid metabolism within the normal range have been associated with risk factors for cardiovascular disease. We performed a genetic association study in 854 members of 224 families to assess the relationship between polymorphisms in the gene coding for hexose-6 phosphate dehydrogenase (H6PD) and carotid intima-medial thickness. METHODS: Families were ascertained via a hypertensive proband. CIMT was measured using B-mode ultrasound. Single nucleotide polymorphisms (SNPs) tagging common variation in the H6PD gene were genotyped. Association was assessed following adjustment for significant covariates including "classical" cardiovascular risk factors. Functional studies to determine the effect of particular SNPs on H6PDH were performed. RESULTS: There was evidence of association between the single nucleotide polymorphism rs17368528 in exon five of the H6PD gene, which encodes an amino-acid change from proline to leucine in the H6PDH protein, and mean carotid intima-medial thickness (p = 0.00065). Genotype was associated with a 5% (or 0.04 mm) higher mean carotid intima-medial thickness measurement per allele, and determined 2% of the population variability in the phenotype. CONCLUSIONS: Our results suggest a novel role for the H6PD gene in atherosclerosis susceptibility.
  • Thumbnail Image
    Item
    Open Access
    Genotypic diversity is associated with clinical outcome and phenotype in cryptococcal meningitis across Southern Africa
    (Public Library of Science, 2015) Beale, Mathew A; Sabiiti, Wilber; Robertson, Emma J; Fuentes-Cabrejo, Karen M; O'Hanlon, Simon J; Jarvis, Joseph N; Loyse, Angela; Meintjes, Graeme; Harrison, Thomas S; May, Robin C
    Cryptococcal meningitis is a major cause of mortality throughout the developing world, yet little is known about the genetic markers underlying Cryptococcal virulence and patient outcome. We studied a cohort of 230 Cryptococcus neoformans (Cn) isolates from HIV-positive South African clinical trial patients with detailed clinical follow-up using multi-locus sequence typing and in vitro phenotypic virulence assays, correlating these data with clinical and fungal markers of disease in the patient. South African Cn displayed high levels of genetic diversity and locus variability compared to globally distributed types, and we identified 50 sequence types grouped within the main molecular types VNI, VNII and VNB, with 72% of isolates typed into one of seven 'high frequency' sequence types. Spatial analysis of patients' cryptococcal genotype was not shown to be clustered geographically, which might argue against recent local acquisition and in favour of reactivation of latent infection. Through comparison of MLST genotyping data with clinical parameters, we found a relationship between genetic lineage and clinical outcome, with patients infected with the VNB lineage having significantly worse survival (n=8, HR 3.35, CI 1.51-7.20, p=0.003), and this was maintained even after adjustment for known prognostic indicators and treatment regimen. Comparison of fungal genotype with in vitro phenotype (phagocytosis, laccase activity and CSF survival) performed on a subset of 89 isolates revealed evidence of lineage-associated virulence phenotype, with the VNII lineage displaying increased laccase activity (p=0.001) and ex vivo CSF survival (p=0.0001). These findings show that Cryptococcus neoformans is a phenotypically heterogeneous pathogen, and that lineage plays an important role in cryptococcal virulence during human infection. Furthermore, a detailed understanding of the genetic diversity in Southern Africa will support further investigation into how genetic diversity is structured across African environments, allowing assessment of the risks different ecotypes pose to infection.
  • Thumbnail Image
    Item
    Open Access
    HIV evolution in early infection: selection pressures, patterns of insertion and deletion, and the impact of APOBEC
    (Public Library of Science, 2009) Wood, Natasha; Bhattacharya, Tanmoy; Keele, Brandon F; Giorgi, Elena; Liu, Michael; Gaschen, Brian; Daniels, Marcus; Ferrari, Guido; Haynes, Barton F; McMichael, Andrew
    Author Summary HIV is a rapidly evolving virus, displaying enormous genetic diversity between and even within infected individuals, with implications for vaccine design and drug treatment. Yet, recent research has shown that most new infections result from transmission of a single virus resulting in a homogeneous viral population in early infection. The process of diversification from the transmitted virus provides information about the selection pressures experienced by the virus during the establishment of a new infection. In this paper, we studied early diversification of the envelope gene in a cohort of 81 subjects acutely infected with HIV-1 subtype B and found evidence of adaptive evolution, with a proportion of sites that tended to diversify more rapidly than expected under a model of neutral evolution. Several of these rapidly diversifying sites facilitate escape from early cytotoxic immune responses. Interestingly, hypermutation of the virus, brought about by host proteins as a strategy to restrict infection, appeared to be associated with early immune escape. In addition to single base substitutions, insertions and deletions are an important aspect of HIV evolution. We show that insertion and deletion mutations occur evenly across the gene, but are preferentially fixed in the variable loop regions.
  • Thumbnail Image
    Item
    Open Access
    A multilocus phylogeny of the world Sycoecinae fig wasps (Chalcidoidea: Pteromalidae)
    (Public Library of Science, 2013) Cruaud, Astrid; Underhill, Jenny G; Huguin, Maïlis; Genson, Gwenaëlle; Jabbour-Zahab, Roula; Tolley, Krystal A; Rasplus, Jean-Yves; van Noort, Simon
    The Sycoecinae is one of five chalcid subfamilies of fig wasps that are mostly dependent on Ficus inflorescences for reproduction. Here, we analysed two mitochondrial ( COI , Cyt b ) and four nuclear genes (ITS2, EF-1α, RpL27a, mago nashi ) from a worldwide sample of 56 sycoecine species. Various alignment and partitioning strategies were used to test the stability of major clades. All topologies estimated using maximum likelihood and Bayesian methods were similar and well resolved but did not support the existing classification. A high degree of morphological convergence was highlighted and several species appeared best described as species complexes. We therefore proposed a new classification for the subfamily. Our analyses revealed several cases of probable speciation on the same host trees (up to 8 closely related species on one single tree of F. sumatrana ), which raises the question of how resource partitioning occurs to avoid competitive exclusion. Comparisons of our results with fig phylogenies showed that, despite sycoecines being internally ovipositing wasps host-switches are common incidents in their evolutionary history. Finally, by studying the evolutionary properties of the markers we used and profiling their phylogenetic informativeness, we predicted their utility for resolving phylogenetic relationships of Chalcidoidea at various taxonomic levels.
  • Thumbnail Image
    Item
    Open Access
    Non-Negative Matrix Factorization for Learning Alignment-Specific Models of Protein Evolution
    (Public Library of Science, 2011) Murrell, Ben; Weighill, Thomas; Buys, Jan; Ketteringham, Robert; Moola, Sasha; Benade, Gerdus; Buisson, Lise du; Kaliski, Daniel; Hands, Tristan; Scheffler, Konrad
    Models of protein evolution currently come in two flavors: generalist and specialist. Generalist models (e.g. PAM, JTT, WAG) adopt a one-size-fits-all approach, where a single model is estimated from a number of different protein alignments. Specialist models (e.g. mtREV, rtREV, HIVbetween) can be estimated when a large quantity of data are available for a single organism or gene, and are intended for use on that organism or gene only. Unsurprisingly, specialist models outperform generalist models, but in most instances there simply are not enough data available to estimate them. We propose a method for estimating alignment-specific models of protein evolution in which the complexity of the model is adapted to suit the richness of the data. Our method uses non-negative matrix factorization (NNMF) to learn a set of basis matrices from a general dataset containing a large number of alignments of different proteins, thus capturing the dimensions of important variation. It then learns a set of weights that are specific to the organism or gene of interest and for which only a smaller dataset is available. Thus the alignment-specific model is obtained as a weighted sum of the basis matrices. Having been constrained to vary along only as many dimensions as the data justify, the model has far fewer parameters than would be required to estimate a specialist model. We show that our NNMF procedure produces models that outperform existing methods on all but one of 50 test alignments. The basis matrices we obtain confirm the expectation that amino acid properties tend to be conserved, and allow us to quantify, on specific alignments, how the strength of conservation varies across different properties. We also apply our new models to phylogeny inference and show that the resulting phylogenies are different from, and have improved likelihood over, those inferred under standard models.
  • Thumbnail Image
    Item
    Open Access
    A novel diagnostic target in the hepatitis C virus genome
    (Public Library of Science, 2009) Drexler, Jan Felix; Kupfer, Bernd; Petersen, Nadine; Grotto, Rejane Maria Tommasini; Rodrigues, Silvia Maria Corvino; Grywna, Klaus; Panning, Marcus; Annan, Augustina; Silva, Giovanni Faria; Douglas, Jill
    Christian Drosten and colleagues develop, validate, and make openly available a prototype hepatitis C virus assay based on the conserved 3' X-tail element, with potential for clinical use in developing countries.
  • Thumbnail Image
    Item
    Open Access
    Phylogenetic Evidence That Two Distinct Trichuris Genotypes Infect both Humans and Non-Human Primates
    (Public Library of Science, 2012) Ravasi, Damiana F.; O’Riain, Mannus J.; Davids, Faezah; Illing, Nicola
    Although there has been extensive debate about whether Trichuris suis and Trichuris trichiura are separate species, only one species of the whipworm T. trichiura has been considered to infect humans and non-human primates. In order to investigate potential cross infection of Trichuris sp. between baboons and humans in the Cape Peninsula, South Africa, we sequenced the ITS1-5.8S-ITS2 region of adult Trichuris sp. worms isolated from five baboons from three different troops, namely the Cape Peninsula troop, Groot Olifantsbos troop and Da Gama Park troop. This region was also sequenced from T. trichiura isolated from a human patient from central Africa (Cameroon) for comparison. By combining this dataset with Genbank records for Trichuris isolated from other humans, non-human primates and pigs from several different countries in Europe, Asia, and Africa, we confirmed the identification of two distinct Trichuris genotypes that infect primates. Trichuris sp. isolated from the Peninsula baboons fell into two distinct clades that were found to also infect human patients from Cameroon, Uganda and Jamaica (named the CP-GOB clade) and China, Thailand, the Czech Republic, and Uganda (named the DG clade), respectively. The divergence of these Trichuris clades is ancient and precedes the diversification of T. suis which clustered closely to the CP-GOB clade. The identification of two distinct Trichuris genotypes infecting both humans and non-human primates is important for the ongoing treatment of Trichuris which is estimated to infect 600 million people worldwide. Currently baboons in the Cape Peninsula, which visit urban areas, provide a constant risk of infection to local communities. A reduction in spatial overlap between humans and baboons is thus an important measure to reduce both cross-transmission and zoonoses of helminthes in Southern Africa.
  • Thumbnail Image
    Item
    Open Access
    Pleistocene aridification cycles shaped the contemporary genetic architecture of southern african baboons
    (Public Library of Science, 2015) Sithaldeen, Riashna; Ackermann, Rebecca Rogers; Bishop, Jacqueline M
    Plio-Pleistocene environmental change influenced the evolutionary history of many animal lineages in Africa, highlighting key roles for both climate and tectonics in the evolution of Africa’s faunal diversity. Here, we explore diversification in the southern African chacma baboon Papio ursinus sensu lato and reveal a dominant role for increasingly arid landscapes during past glacial cycles in shaping contemporary genetic structure. Recent work on baboons ( Papio spp.) supports complex lineage structuring with a dominant pulse of diversification occurring 1-2Ma, and yet the link to palaeoenvironmental change remains largely untested. Phylogeographic reconstruction based on mitochondrial DNA sequence data supports a scenario where chacma baboon populations were likely restricted to refugia during periods of regional cooling and drying through the Late Pleistocene. The two lineages of chacma baboon, ursinus and griseipes , are strongly geographically structured, and demographic reconstruction together with spatial analysis of genetic variation point to possible climate-driven isolating events where baboons may have retreated to more optimum conditions during cooler, drier periods. Our analysis highlights a period of continuous population growth beginning in the Middle to Late Pleistocene in both the ursinus and the PG2 griseipes lineages. All three clades identified in the study then enter a state of declining population size (Ne f ) through to the Holocene; this is particularly marked in the last 20,000 years, most likely coincident with the Last Glacial Maximum. The pattern recovered here conforms to expectations based on the dynamic regional climate trends in southern Africa through the Pleistocene and provides further support for complex patterns of diversification in the region’s biodiversity.
  • Thumbnail Image
    Item
    Open Access
    Scoring protein relationships in functional interaction networks predicted from sequence data
    (Public Library of Science, 2011) Mazandu, Gaston K; Mulder, Nicola J
    The abundance of diverse biological data from various sources constitutes a rich source of knowledge, which has the power to advance our understanding of organisms. This requires computational methods in order to integrate and exploit these data effectively and elucidate local and genome wide functional connections between protein pairs, thus enabling functional inferences for uncharacterized proteins. These biological data are primarily in the form of sequences, which determine functions, although functional properties of a protein can often be predicted from just the domains it contains. Thus, protein sequences and domains can be used to predict protein pair-wise functional relationships, and thus contribute to the function prediction process of uncharacterized proteins in order to ensure that knowledge is gained from sequencing efforts. In this work, we introduce information-theoretic based approaches to score protein-protein functional interaction pairs predicted from protein sequence similarity and conserved protein signature matches. The proposed schemes are effective for data-driven scoring of connections between protein pairs. We applied these schemes to the Mycobacterium tuberculosis proteome to produce a homology-based functional network of the organism with a high confidence and coverage. We use the network for predicting functions of uncharacterised proteins. Availability Protein pair-wise functional relationship scores for Mycobacterium tuberculosis strain CDC1551 sequence data and python scripts to compute these scores are available at http://web.cbio.uct.ac.za/~gmazandu/scoringschemes .