skip to main content


Title: Comparative Phylogenomic Analysis Reveals Evolutionary Genomic Changes and Novel Toxin Families in Endophytic Liberibacter Pathogens
ABSTRACT Liberibacter pathogens are the causative agents of several severe crop diseases worldwide, including citrus Huanglongbing and potato zebra chip. These bacteria are endophytic and nonculturable, which makes experimental approaches challenging and highlights the need for bioinformatic analysis in advancing our understanding about Liberibacter pathogenesis. Here, we performed an in-depth comparative phylogenomic analysis of the Liberibacter pathogens and their free-living, nonpathogenic, ancestral species, aiming to identify major genomic changes and determinants associated with their evolutionary transitions in living habitats and pathogenicity. Using gene neighborhood analysis and phylogenetic classification, we systematically uncovered, annotated, and classified all prophage loci into four types, including one previously unrecognized group. We showed that these prophages originated through independent gene transfers at different evolutionary stages of Liberibacter and only the SC-type prophage was associated with the emergence of the pathogens. Using ortholog clustering, we vigorously identified two additional sets of genomic genes, which were either lost or gained in the ancestor of the pathogens. Consistent with the habitat change, the lost genes were enriched for biosynthesis of cellular building blocks. Importantly, among the gained genes, we uncovered several previously unrecognized toxins, including new toxins homologous to the EspG/VirA effectors, a YdjM phospholipase toxin, and a secreted endonuclease/exonuclease/phosphatase (EEP) protein. Our results substantially extend the knowledge of the evolutionary events and potential determinants leading to the emergence of endophytic, pathogenic Liberibacter species, which will facilitate the design of functional experiments and the development of new methods for detection and blockage of these pathogens. IMPORTANCE Liberibacter pathogens are associated with several severe crop diseases, including citrus Huanglongbing, the most destructive disease to the citrus industry. Currently, no effective cure or treatments are available, and no resistant citrus variety has been found. The fact that these obligate endophytic pathogens are not culturable has made it extremely challenging to experimentally uncover the genes/proteins important to Liberibacter pathogenesis. Further, earlier bioinformatics studies failed to identify key genomic determinants, such as toxins and effector proteins, that underlie the pathogenicity of the bacteria. In this study, an in-depth comparative genomic analysis of Liberibacter pathogens along with their ancestral nonpathogenic species identified the prophage loci and several novel toxins that are evolutionarily associated with the emergence of the pathogens. These results shed new light on the disease mechanism of Liberibacter pathogens and will facilitate the development of new detection and blockage methods targeting the toxins.  more » « less
Award ID(s):
1928770
NSF-PAR ID:
10308495
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  
Editor(s):
Burbank, Lindsey Price
Date Published:
Journal Name:
Microbiology Spectrum
Volume:
9
Issue:
2
ISSN:
2165-0497
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Miller, Samuel I. (Ed.)
    ABSTRACT The emergence of fluoroquinolone resistance in nosocomial pathogens has restricted the clinical efficacy of this antibiotic class. In Acinetobacter baumannii , the majority of clinical isolates now show high-level resistance due to mutations in gyrA (DNA gyrase) and parC (topoisomerase IV [topo IV]). To investigate the molecular basis for fluoroquinolone resistance, an exhaustive mutation analysis was performed in both drug-sensitive and -resistant strains to identify loci that alter ciprofloxacin sensitivity. To this end, parallel fitness tests of over 60,000 unique insertion mutations were performed in strains with various alleles in genes encoding the drug targets. The spectra of mutations that altered drug sensitivity were found to be similar in the drug-sensitive and gyrA parC double-mutant backgrounds, having resistance alleles in both genes. In contrast, the introduction of a single gyrA resistance allele, resulting in preferential poisoning of topo IV by ciprofloxacin, led to extreme alterations in the insertion mutation fitness landscape. The distinguishing feature of preferential topo IV poisoning was enhanced induction of DNA synthesis in the region of two endogenous prophages, with DNA synthesis associated with excision and circularization of the phages. Induction of the selective DNA synthesis in the gyrA background was also linked to heightened prophage gene transcription and enhanced activation of the mutagenic SOS response relative to that observed in either the wild-type (WT) or gyrA parC double mutant. Therefore, the accumulation of mutations that result in the stepwise evolution of high ciprofloxacin resistance is tightly connected to modulation of the SOS response and endogenous prophage DNA synthesis. IMPORTANCE Fluoroquinolones have been extremely successful antibiotics due to their ability to target multiple bacterial enzymes critical to DNA replication, the topoisomerases DNA gyrase and topo IV. Unfortunately, mutations lowering drug affinity for both enzymes are now widespread, rendering these drugs ineffective for many pathogens. To undermine this form of resistance, we examined how bacteria with target alterations differentially cope with fluoroquinolone exposures. We studied this problem in the nosocomial pathogen A. baumannii , which causes drug-resistant life-threatening infections. Employing genome-wide approaches, we uncovered numerous pathways that could be exploited to raise fluoroquinolone sensitivity independently of target alteration. Remarkably, fluoroquinolone targeting of topo IV in specific mutants caused dramatic hyperinduction of prophage replication and enhanced the mutagenic DNA damage response, but these responses were muted in strains with DNA gyrase as the primary target. This work demonstrates that resistance evolution via target modification can profoundly modulate the antibiotic stress response, revealing potential resistance-associated liabilities. 
    more » « less
  2. Gilbert, Jack A. (Ed.)
    ABSTRACT Host association—the selective adaptation of pathogens to specific host species—evolves through constant interactions between host and pathogens, leaving a lot yet to be discovered on immunological mechanisms and genomic determinants. The causative agents of Lyme disease (LD) are spirochete bacteria composed of multiple species of the Borrelia burgdorferi sensu lato complex, including B. burgdorferi ( Bb ), the main LD pathogen in North America—a useful model for the study of mechanisms underlying host-pathogen association. Host adaptation requires pathogens’ ability to evade host immune responses, such as complement, the first-line innate immune defense mechanism. We tested the hypothesis that different host-adapted phenotypes among Bb strains are linked to polymorphic loci that confer complement evasion traits in a host-specific manner. We first examined the survivability of 20 Bb strains in sera in vitro and/or bloodstream and tissues in vivo from rodent and avian LD models. Three groups of complement-dependent host-association phenotypes emerged. We analyzed complement-evasion genes, identified a priori among all strains and sequenced and compared genomes for individual strains representing each phenotype. The evolutionary history of ospC loci is correlated with host-specific complement-evasion phenotypes, while comparative genomics suggests that several gene families and loci are potentially involved in host association. This multidisciplinary work provides novel insights into the functional evolution of host-adapted phenotypes, building a foundation for further investigation of the immunological and genomic determinants of host association. IMPORTANCE Host association is the phenotype that is commonly found in many pathogens that preferential survive in particular hosts. The Lyme disease (LD)-causing agent, B. burgdorferi ( Bb ), is an ideal model to study host association, as Bb is mainly maintained in nature through rodent and avian hosts. A widespread yet untested concept posits that host association in Bb strains is linked to Bb functional genetic variation conferring evasion to complement, an innate defense mechanism in vertebrate sera. Here, we tested this concept by grouping 20 Bb strains into three complement-dependent host-association phenotypes based on their survivability in sera and/or bloodstream and distal tissues in rodent and avian LD models. Phylogenomic analysis of these strains further correlated several gene families and loci, including ospC , with host-specific complement-evasion phenotypes. Such multifaceted studies thus pave the road to further identify the determinants of host association, providing mechanistic insights into host-pathogen interaction. 
    more » « less
  3. ABSTRACT The Fusarium oxysporum species complex (FOSC) is a group of soilborne pathogens causing severe disease in more than 100 plant hosts, while individual strains exhibit strong host specificity. Both chromosome transfer and comparative genomics experiments have demonstrated that lineage-specific (LS) chromosomes contribute to the host-specific pathogenicity. However, little is known about the functional importance of genes encoded in these LS chromosomes. Focusing on signaling transduction, this study compared the kinomes of 12 F. oxysporum isolates, including both plant and human pathogens and 1 nonpathogenic biocontrol strain, with 7 additional publicly available ascomycete genomes. Overall, F. oxysporum kinomes are the largest, facilitated in part by the acquisitions of the LS chromosomes. The comparative study identified 99 kinases that are present in almost all examined fungal genomes, forming the core signaling network of ascomycete fungi. Compared to the conserved ascomycete kinome, the expansion of the F. oxysporum kinome occurs in several kinase families such as histidine kinases that are involved in environmental signal sensing and target of rapamycin (TOR) kinase that mediates cellular responses. Comparative kinome analysis suggests a convergent evolution that shapes individual F. oxysporum isolates with an enhanced and unique capacity for environmental perception and associated downstream responses. IMPORTANCE Isolates of Fusarium oxysporum are adapted to survive a wide range of host and nonhost conditions. In addition, F. oxysporum was recently recognized as the top emerging opportunistic fungal pathogen infecting immunocompromised humans. The sensory and response networks of these fungi undoubtedly play a fundamental role in establishing the adaptability of this group. We have examined the kinomes of 12 F. oxysporum isolates and highlighted kinase families that distinguish F. oxysporum from other fungi, as well as different isolates from one another. The amplification of kinases involved in environmental signal relay and regulating downstream cellular responses clearly sets Fusarium apart from other Ascomycetes . Although the functions of many of these kinases are still unclear, their specific proliferation highlights them as a result of the evolutionary forces that have shaped this species complex and clearly marks them as targets for exploitation in order to combat disease. 
    more » « less
  4. Alexandre, Gladys (Ed.)
    ABSTRACT Xylella fastidiosa infects several economically important crops in the Americas, and it also recently emerged in Europe. Here, using a set of Xylella genomes reflective of the genus-wide diversity, we performed a pan-genome analysis based on both core and accessory genes for two purposes: (i) to test associations between genetic divergence and plant host species and (ii) to identify positively selected genes that are potentially involved in arms-race dynamics. For the former, tests yielded significant evidence for the specialization of X. fastidiosa to plant host species. This observation contributes to a growing literature suggesting that the phylogenetic history of X. fastidiosa lineages affects the host range. For the latter, our analyses uncovered evidence of positive selection across codons for 5.3% (67 of 1,257) of the core genes and 5.4% (201 of 3,691) of the accessory genes. These genes are candidates to encode interacting factors with plant and insect hosts. Most of these genes had unknown functions, but we did identify some tractable candidates, including nagZ_2 , which encodes a beta-glucosidase that is important for Neisseria gonorrhoeae biofilm formation; cya , which modulates gene expression in pathogenic bacteria, and barA , a membrane associated histidine kinase that has roles in cell division, metabolism, and pili formation. IMPORTANCE Xylella fastidiosa causes devasting diseases to several critical crops. Because X. fastidiosa colonizes and infects many plant species, it is important to understand whether the genome of X. fastidiosa has genetic determinants that underlie specialization to specific host plants. We analyzed genome sequences of X. fastidiosa to investigate evolutionary relationships and to test for evidence of positive selection on specific genes. We found a significant signal between genome diversity and host plants, consistent with bacterial specialization to specific plant hosts. By screening for positive selection, we identified both core and accessory genes that may affect pathogenicity, including genes involved in biofilm formation. 
    more » « less
  5. Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome ofBothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2(PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option.

     
    more » « less