Comparative genomics, pangenomics, and phenomic studies of Pectobacterium betavasculorum strains isolated from sugar beet, potato, sunflower, and artichoke: insights into pathogenicity, virulence determinants, and adaptation to the host plant

Introduction

Bacteria of genus Pectobacterium, encompassing economically significant pathogens affecting various plants, includes the species P. betavasculorum, initially associated with beetroot infection. However, its host range is much broader. It causes diseases of sunflower, potato, tomato, carrots, sweet potato, radish, squash, cucumber, and chrysanthemum. To explain this phenomenon, a comprehensive pathogenomic and phenomic characterisation of P. betavasculorum species was performed.

Methods

Genomes of P. betavasculorum strains isolated from potato, sunflower, and artichoke were sequenced and compared with those from sugar beet isolates. Metabolic profiling and pathogenomic analyses were conducted to assess virulence determinants and adaptation potential. Pathogenicity assays were performed on potato tubers and chicory leaves to confirm in silico predictions of disease symptoms. Phenotypic assays were also conducted to assess the strains ability to synthesise homoserine lactones and siderophores.

Results

The genome size ranged from 4.675 to 4.931 kbp, and GC % was between 51.0% and 51.2%. The pangenome of P. betavasculorum is open and comprises, on average, 4,220 gene families. Of these, 83% of genes are the core genome, and 2% of the entire pangenome are unique genes. Strains isolated from sugar beet have a smaller pangenome size and a higher number of unique genes than those from other plants. Interestingly, genomes of strains from artichoke and sunflower share 391 common CDS that are not present in the genomes of other strains from sugar beet or potato. Those strains have only one unique gene. All strains could use numerous sugars as building materials and energy sources and possessed a high repertoire of virulence determinants in the genomes. P. betavasculorum strains were able to cause disease symptoms on potato tubers and chicory leaves. They were also able to synthesise homoserine lactones and siderophores.

Discussion

The findings underscore the adaptability of P. betavasculorum to diverse hosts and environments. Strains adapted to plants with high sugar content in tissues have a different composition of fatty acids in membranes and a different mechanism of replenishing nitrogen in case of deficiency of this compound than strains derived from other plant species. Extensive phenomics and genomic analyses performed in this study have shown that P. betavasculorum species is an agronomically relevant pathogen.

Clonal reproduction of Moniliophthora roreri and the emergence of unique lineages with distinct genomes during range expansion

The basidiomycete Moniliophthora roreri causes frosty pod rot of cacao (Theobroma cacao) in the western hemisphere. Moniliophthora roreri is considered asexual and haploid throughout its hemibiotrophic life cycle. To understand the processes driving genome modification, using long-read sequencing technology, we sequenced and assembled 5 high-quality M. roreri genomes out of a collection of 99 isolates collected throughout the pathogen's range. We obtained chromosome-scale assemblies composed of 11 scaffolds. We used short-read technology to sequence the genomes of 22 similarly chosen isolates. Alignments among the 5 reference assemblies revealed inversions, translocations, and duplications between and within scaffolds. Isolates at the front of the pathogens' expanding range tend to share lineage-specific structural variants, as confirmed by short-read sequencing. We identified, for the first time, 3 new mating type A locus alleles (5 in total) and 1 new potential mating type B locus allele (3 in total). Currently, only 2 mating type combinations, A1B1 and A2B2, are known to exist outside of Colombia. A systematic survey of the M. roreri transcriptome across 2 isolates identified an expanded candidate effector pool and provided evidence that effector candidate genes unique to the Moniliophthoras are preferentially expressed during the biotrophic phase of disease. Notably, M. roreri isolates in Costa Rica carry a chromosome segment duplication that has doubled the associated gene complement and includes secreted proteins and candidate effectors. Clonal reproduction of the haploid M. roreri genome has allowed lineages with unique genome structures and compositions to dominate as it expands its range, displaying a significant founder effect.

Mycobacterium bovis Strain Ravenel Is Attenuated in Cattle

Mycobacterium tuberculosis variant bovis (MBO) has one of the widest known mammalian host ranges, including humans. Despite the characterization of this pathogen in the 1800s and whole genome sequencing of a UK strain (AF2122) nearly two decades ago, the basis of its host specificity and pathogenicity remains poorly understood. Recent experimental calf infection studies show that MBO strain Ravenel (MBO Ravenel) is attenuated in the cattle host compared to other pathogenic strains of MBO. In the present study, experimental infections were performed to define attenuation. Whole genome sequencing was completed to identify regions of differences (RD) and single nucleotide polymorphisms (SNPs) to explain the observed attenuation. Comparative genomic analysis of MBO Ravenel against three pathogenic strains of MBO (strains AF2122-97, 10-7428, and 95-1315) was performed. Experimental infection studies on five calves each, with either MBO Ravenel or 95-1315, revealed no visible lesions in all five animals in the Ravenel group despite robust IFN-γ responses. Out of 486 polymorphisms in the present analysis, 173 were unique to MBO Ravenel among the strains compared. A high-confidence subset of nine unique SNPs were missense mutations in genes with annotated functions impacting two major MBO survival and virulence pathways: (1) Cell wall synthesis & transport [espH (A103T), mmpL8 (V888I), aftB (H484Y), eccC₅ (T507M), rpfB (E263G)], and (2) Lipid metabolism & respiration [mycP₁(T125I), pks5 (G455S), fadD29 (N231S), fadE29 (V360G)]. These substitutions likely contribute to the observed attenuation. Results from experimental calf infections and the functional attributions of polymorphic loci on the genome of MBO Ravenel provide new insights into the strain's genotype-disease phenotype associations.

Natural language processing approach to model the secretion signal of type III effectors

Type III effectors are proteins injected by Gram-negative bacteria into eukaryotic hosts. In many plant and animal pathogens, these effectors manipulate host cellular processes to the benefit of the bacteria. Type III effectors are secreted by a type III secretion system that must "classify" each bacterial protein into one of two categories, either the protein should be translocated or not. It was previously shown that type III effectors have a secretion signal within their N-terminus, however, despite numerous efforts, the exact biochemical identity of this secretion signal is generally unknown. Computational characterization of the secretion signal is important for the identification of novel effectors and for better understanding the molecular translocation mechanism. In this work we developed novel machine-learning algorithms for characterizing the secretion signal in both plant and animal pathogens. Specifically, we represented each protein as a vector in high-dimensional space using Facebook's protein language model. Classification algorithms were next used to separate effectors from non-effector proteins. We subsequently curated a benchmark dataset of hundreds of effectors and thousands of non-effector proteins. We showed that on this curated dataset, our novel approach yielded substantially better classification accuracy compared to previously developed methodologies. We have also tested the hypothesis that plant and animal pathogen effectors are characterized by different secretion signals. Finally, we integrated the novel approach in Effectidor, a web-server for predicting type III effector proteins, leading to a more accurate classification of effectors from non-effectors.

Wheat stem rust recorded for the first time in decades in Ireland

Wheat stem rust, caused by the fungus Puccinia graminis f. sp. tritici (Pgt), occurs in most wheat-growing areas worldwide, and, in western Europe since 2013, has started to re-emerge after many decades of absence. Following this trend across western Europe, in 2020, we also detected and recorded wheat stem rust for the first time in five decades in experimental plots across five locations in Ireland. To examine the potential origin of the Irish Pgt infection in 2020, we carried out transcriptome sequencing on 12 Pgt-infected wheat samples collected across Ireland and compared these to 76 global P. graminis isolates. This analysis identified a close genetic relationship between the Irish Pgt isolates and those from Ethiopia collected in 2015 after a severe stem rust epidemic caused by the TKTTF Pgt race, and with the UK-01 Pgt isolate that was previously assigned to the TKTTF race. Subsequent pathology-based race profiling designated two Irish isolates and recent UK and French Pgt isolates to the TKTTF Pgt race group. This suggests that the Irish Pgt occurrence most probably originated from recent long-distance windborne dispersal of Pgt urediniospores from neighbouring countries in Europe where we confirmed the Pgt TKTTF race continues to be prevalent. The identification of wheat stem rust in Ireland at multiple locations in 2020 illustrates that the disease can occur in Ireland and emphasizes the need to re-initiate local monitoring for this re-emergent threat to wheat production across western Europe.

Antibiotic Resistance and Pathogenomics of Staphylococci Circulating in Novosibirsk, Russia

A total of 394 strains of staphylococci found in humans and pets in Novosibirsk, Siberian Russia, were characterized in terms of antibiotic resistance and corresponding genes. Two coagulase-positive and 17 coagulase-negative species were identified. The majority of isolates, with the exception of S. haemolyticus and hospital S. epidermidis isolates, were sensitive to most of the tested antibiotics, and isolates from pets displayed the lowest level of resistance. Nevertheless, methicillin-resistant (MRS) and/or multidrug-resistant (MDR) isolates were found in all prevailed species, including coagulase-negative. A set of genes corresponding to the detected resistance was identified: mecA (beta-lactam resistance), aac(6')-Ie-aph(2″)-Ia, aph(3')-IIIa, ant(4')-Ia (aminoglycoside-modifying enzymes), ermA/ermC, and msrA (macrolide resistance). Complete genome analysis for ten MDR S. epidermidis and five MDR S. haemolyticus isolates revealed additional antibiotic resistance genes mphC, qacA/qacB, norA, dfrC/dfrG, lnuA, BseSR, and fosB. NorA, dfrC, and fosB were present in all S. epidermidis genomes, whereas mphC and msrA were identified in all S. haemolyticus ones. All investigated MDR S. epidermidis and four of five S. haemolyticus strains were moderate or strong biofilm producers, whereas multiple genes responsible for this function and for virulence and pathogenicity were identified mostly in S. epidermidis, but were less frequently represented in S. haemolyticus.

Comparative Pathogenomics of Escherichia coli: Polyvalent Vaccine Target Identification through Virulome Analysis

Comparative genomics of bacterial pathogens has been useful for revealing potential virulence factors. Escherichia coli is a significant cause of human morbidity and mortality worldwide but can also exist as a commensal in the human gastrointestinal tract. With many sequenced genomes, it has served as a model organism for comparative genomic studies to understand the link between genetic content and potential for virulence. To date, however, no comprehensive analysis of its complete “virulome” has been performed for the purpose of identifying universal or pathotype-specific targets for vaccine development. Here, we describe the construction of a pathotype database of 107 well-characterized completely sequenced pathogenic and nonpathogenic E. coli strains, which we annotated for major virulence factors (VFs). The data are cross referenced for patterns against pathotype, phylogroup, and sequence type, and the results were verified against all 1,348 complete E. coli chromosomes in the NCBI RefSeq database. Our results demonstrate that phylogroup drives many of the “pathotype-associated” VFs, and ExPEC-associated VFs are found predominantly within the B2/D/F/G phylogenetic clade, suggesting that these phylogroups are better adapted to infect human hosts. Finally, we used this information to propose polyvalent vaccine targets with specificity toward extraintestinal strains, targeting key invasive strategies, including immune evasion (group 2 capsule), iron acquisition (FyuA, IutA, and Sit), adherence (SinH, Afa, Pap, Sfa, and Iha), and toxins (Usp, Sat, Vat, Cdt, Cnf1, and HlyA). While many of these targets have been proposed before, this work is the first to examine their pathotype and phylogroup distribution and how they may be targeted together to prevent disease.

Comparative and Evolutionary Genomics of Isolates Provide Insight into the Pathoadaptation of Aeromonas

Aeromonads are ubiquitous aquatic bacteria that cause opportunistic infections in humans, but their pathogenesis remains poorly understood. A pathogenomic approach was undertaken to provide insights into the emergence and evolution of pathogenic traits in aeromonads. The genomes of 64 Aeromonas strains representative of the whole genus were analyzed to study the distribution, phylogeny, and synteny of the flanking sequences of 13 virulence-associated genes. The reconstructed evolutionary histories varied markedly depending on the gene analyzed and ranged from vertical evolution, which followed the core genome evolution (alt and colAh), to complex evolution, involving gene loss by insertion sequence-driven gene disruption, horizontal gene transfer, and paraphyly with some virulence genes associated with a phylogroup (aer, ser, and type 3 secretion system components) or no phylogroup (type 3 secretion system effectors, Ast, ExoA, and RtxA toxins). The general pathogenomic overview of aeromonads showed great complexity with diverse evolution modes and gene organization and uneven distribution of virulence genes in the genus; the results provided insights into aeromonad pathoadaptation or the ability of members of this group to emerge as pathogens. Finally, these findings suggest that aeromonad virulence-associated genes should be examined at the population level and that studies performed on type or model strains at the species level cannot be generalized to the whole species.

Pathogenomics of Uterine Fibroids Development

We review recent studies dealing with the molecular genetics and basic results of omics analysis of uterine leiomyoma (LM)-a common benign muscle tumor of the uterus. Whole genome studies of LM resulted in the discovery of many new gene nets and biological pathways, including its origin, transcriptomic, and epigenetic profiles, as well as the impact of the inter-cell matrix in LM growth and involvement of microRNA in its regulation. New data on somatic cell mutations ultimately involved in the origin, distribution and growth of LM are reviewed. Putative identification of LM progenitor SC (stem cells) giving rise to maternal fibroid nodes and junctional zones provide a new clue for hypotheses on the pathogenomics of LM. The reviewed data are consistent with at least two different but probably intimately interacted molecular mechanisms of LM. One of them (the genetic hypothesis) is focused primarily on the MED12 gene mutations and suggests its onset in the side population of embryonic myoblasts of the female reproductive system, which later gave rise to multiple small and medium fibroids. The single and usually large-size fibroids are induced by predominantly epigenetic disorders in LM SC, provoked by enhanced expression of the HMGA2 gene caused by its hypomethylation and epigenetic deregulation enhanced by hypoxia, muscle tension, or chromosome instability/aberrations. The pathogenomics of both genetic and epigenetic programs of LM with many peculiarities at the beginning later became rather similar and partly overlapped due to the proximity of their gene nets and epigenetic landscape. Pathogenomic studies of LM open ways for elaboration of novel strategies of prevention and treatment of this common disease.

Pathogenomics of Emerging Campylobacter Species

Campylobacter is among the four main causes of gastroenteritis worldwide and has increased in both developed and developing countries over the last 10 years. The vast majority of reported Campylobacter infections are caused by Campylobacter jejuni and, to a lesser extent, C. coli; however, the increasing recognition of other emerging Campylobacter pathogens is urgently demanding a better understanding of how these underestimated species cause disease, transmit, and evolve. In parallel to the enhanced clinical awareness of campylobacteriosis due to improved diagnostic protocols, the application of high-throughput sequencing has increased the number of whole-genome sequences available to dozens of strains of many emerging campylobacters. This has allowed for comprehensive comparative pathogenomic analyses for several species, such as C. fetus and C. concisus These studies have started to reveal the evolutionary forces shaping their genomes and have brought to light many genomic features related to pathogenicity in these neglected species, promoting the development of new tools and approaches relevant for clinical microbiology. Despite the need for additional characterization of genomic diversity in emerging campylobacters, the increasing body of literature describing pathogenomic studies on these species deserves to be discussed from an integrative perspective. This review compiles the current knowledge and highlights future work toward deepening our understanding about genome dynamics and the mechanisms governing the evolution of pathogenicity in emerging Campylobacter species, which is urgently needed to develop strategies to prevent or control the spread of these pathogens.

Coagulase-Negative Staphylococci Pathogenomics

Coagulase-negative Staphylococci (CoNS) are skin commensal bacteria. Besides their role in maintaining homeostasis, CoNS have emerged as major pathogens in nosocomial settings. Several studies have investigated the molecular basis for this emergence and identified multiple putative virulence factors with regards to Staphylococcus aureus pathogenicity. In the last decade, numerous CoNS whole-genome sequences have been released, leading to the identification of numerous putative virulence factors. Koch’s postulates and the molecular rendition of these postulates, established by Stanley Falkow in 1988, do not explain the microbial pathogenicity of CoNS. However, whole-genome sequence data has shed new light on CoNS pathogenicity. In this review, we analyzed the contribution of genomics in defining CoNS virulence, focusing on the most frequent and pathogenic CoNS species: S. epidermidis, S. haemolyticus, S. saprophyticus, S. capitis, and S. lugdunensis.

Advances in Wheat and Pathogen Genomics: Implications for Disease Control

The gene pool of wheat and its wild and domesticated relatives contains a plethora of resistance genes that can be exploited to make wheat more resilient to pathogens. Only a few of these genes have been isolated and studied at the molecular level. In recent years, we have seen a shift from classical breeding to genomics-assisted breeding, which makes use of the enormous advancements in DNA sequencing and high-throughput molecular marker technologies for wheat improvement. These genomic advancements have the potential to transform wheat breeding in the near future and to significantly increase the speed and precision at which new cultivars can be bred. This review highlights the genomic improvements that have been made in wheat and its pathogens over the past years and discusses their implications for disease-resistance breeding.

Genomic Characterization and Copy Number Variation of Bacillus anthracis Plasmids pXO1 and pXO2 in a Historical Collection of 412 Strains

Bacillus anthracis microorganisms are of historical and epidemiological importance and are among the most homogenous bacterial groups known, even though the B. anthracis genome is rich in mobile elements. Mobile elements can trigger the diversification of lineages; therefore, characterizing the extent of genomic variation in a large collection of strains is critical for a complete understanding of the diversity and evolution of the species. Here, we sequenced a large collection of B. anthracis strains (>400) that were recovered from human, animal, and environmental sources around the world. Our results confirmed the remarkable stability of gene content and synteny of the anthrax plasmids and revealed no signal of plasmid exchange between B. anthracis and pathogenic B. cereus isolates but rather predominantly vertical descent. These findings advance our understanding of the biology and pathogenomic evolution of B. anthracis and its plasmids.

Bacillus anthracis plasmids pXO1 and pXO2 carry the main virulence factors responsible for anthrax. However, the extent of copy number variation within the species and how the plasmids are related to pXO1/pXO2-like plasmids in other species of the Bacillus cereussensu lato group remain unclear. To gain new insights into these issues, we sequenced 412 B. anthracis strains representing the total phylogenetic and ecological diversity of the species. Our results revealed that B. anthracis genomes carried, on average, 3.86 and 2.29 copies of pXO1 and pXO2, respectively, and also revealed a positive linear correlation between the copy numbers of pXO1 and pXO2. No correlation between the plasmid copy number and the phylogenetic relatedness of the strains was observed. However, genomes of strains isolated from animal tissues generally maintained a higher plasmid copy number than genomes of strains from environmental sources (P < 0.05 [Welch two-sample t test]). Comparisons against B. cereus genomes carrying complete or partial pXO1-like and pXO2-like plasmids showed that the plasmid-based phylogeny recapitulated that of the main chromosome, indicating limited plasmid horizontal transfer between or within these species. Comparisons of gene content revealed a closed pXO1 and pXO2 pangenome; e.g., plasmids encode <8 unique genes, on average, and a single large fragment deletion of pXO1 in one B. anthracis strain (2000031682) was detected. Collectively, our results provide a more complete view of the genomic diversity of B. anthracis plasmids, their copy number variation, and the virulence potential of other Bacillus species carrying pXO1/pXO2-like plasmids.

IMPORTANCEBacillus anthracis microorganisms are of historical and epidemiological importance and are among the most homogenous bacterial groups known, even though the B. anthracis genome is rich in mobile elements. Mobile elements can trigger the diversification of lineages; therefore, characterizing the extent of genomic variation in a large collection of strains is critical for a complete understanding of the diversity and evolution of the species. Here, we sequenced a large collection of B. anthracis strains (>400) that were recovered from human, animal, and environmental sources around the world. Our results confirmed the remarkable stability of gene content and synteny of the anthrax plasmids and revealed no signal of plasmid exchange between B. anthracis and pathogenic B. cereus isolates but rather predominantly vertical descent. These findings advance our understanding of the biology and pathogenomic evolution of B. anthracis and its plasmids.

Whole-Genome Resequencing and Pan-Transcriptome Reconstruction Highlight the Impact of Genomic Structural Variation on Secondary Metabolite Gene Clusters in the Grapevine Esca Pathogen Phaeoacremonium minimum

The Ascomycete fungus Phaeoacremonium minimum is one of the primary causal agents of Esca, a widespread and damaging grapevine trunk disease. Variation in virulence among Pm. minimum isolates has been reported, but the underlying genetic basis of the phenotypic variability remains unknown. The goal of this study was to characterize intraspecific genetic diversity and explore its potential impact on virulence functions associated with secondary metabolism, cellular transport, and cell wall decomposition. We generated a chromosome-scale genome assembly, using single molecule real-time sequencing, and resequenced the genomes and transcriptomes of multiple isolates to identify sequence and structural polymorphisms. Numerous insertion and deletion events were found for a total of about 1 Mbp in each isolate. Structural variation in this extremely gene dense genome frequently caused presence/absence polymorphisms of multiple adjacent genes, mostly belonging to biosynthetic clusters associated with secondary metabolism. Because of the observed intraspecific diversity in gene content due to structural variation we concluded that a transcriptome reference developed from a single isolate is insufficient to represent the virulence factor repertoire of the species. We therefore compiled a pan-transcriptome reference of Pm. minimum comprising a non-redundant set of 15,245 protein-coding sequences. Using naturally infected field samples expressing Esca symptoms, we demonstrated that mapping of meta-transcriptomics data on a multi-species reference that included the Pm. minimum pan-transcriptome allows the profiling of an expanded set of virulence factors, including variable genes associated with secondary metabolism and cellular transport.

Pathogenomics of Endometriosis Development

For over 100 years, endometriosis, as a chronic, estrogen-dependent, inflammatory, heritable disease affecting approximately 5–10% of women in reproductive age has been the focus of clinicians and scientists. In spite of numerous environmental, genetic, epigenetic, endocrine, and immunological studies, our knowledge of endometriosis is still fragmentary, and its precise pathophysiology and pathogenomics remain a mystery. The implementation of new technologies has provided tremendous progress in understanding the many intrinsic molecular mechanisms in the development of endometriosis, with progenitor and stem cells (SCs) of the eutopic endometrium as the starting players and endometriotic lesions as the final pathomorphological trait. Novel data on the molecular, genetic, and epigenetic mechanisms of the disease are briefly outlined. We hypothesize the existence of an endometriosis development genetic program (EMDP) that governs the origin of endometrium stem cells programmed for endometriosis (1), their transition (metaplasia) into mesenchymal SCs (2), and their invasion of the peritoneum and progression to endometriotic lesions (3). The pros and cons of the recent unifying theory of endometriosis are also discussed. Complex genomic and epigenetic interactions at different stages of the endometriosis process result in different forms of the disease, with specific features and clinical manifestations. The significance of the EMDP in elaborating a new strategy for endometriosis prediction, prevention, and treatment is discussed.

Pathogenomic Analysis of Wheat Yellow Rust Lineages Detects Seasonal Variation and Host Specificity

Recent disease outbreaks caused by (re-)emerging plant pathogens have been associated with expansions in pathogen geographic distribution and increased virulence. For example, in the past two decades' wheat yellow (stripe) rust, Puccinia striiformis f. sp. tritici, has seen the emergence of new races that are adapted to warmer temperatures, have expanded virulence profiles, and are more aggressive than previous races, leading to wide-scale epidemics. Here, we used field-based genotyping to generate high-resolution data on P. striiformis genetics and carried out global population analysis. We also undertook comparative analysis of the 2014 and 2013 UK populations and assessed the temporal dynamics and host specificity of distinct pathogen genotypes. Our analysis revealed that P. striiformis lineages recently detected in Europe are extremely diverse and in fact similar to globally dispersed populations. In addition, we identified a considerable shift in the UK P. striiformis population structure including the first identification of one infamous race known as Kranich. Next, by establishing the genotype of both the pathogen and host within a single infected field sample, we uncovered evidence for varietal specificity for genetic groups of P. striiformis. Finally, we found potential seasonal specificity for certain genotypes of the pathogen with several lineages identified only in samples collected in late spring and into the summer, whereas one lineage was identified throughout the wheat growing season. Our discovery of which wheat varieties are susceptible to which specific P. striiformis isolates, and when those isolates are prevalent throughout the year, represents a powerful tool for disease management.

Insights from the Genome Sequence of Acidovorax citrulli M6, a Group I Strain of the Causal Agent of Bacterial Fruit Blotch of Cucurbits

Acidovorax citrulli is a seedborne bacterium that causes bacterial fruit blotch of cucurbit plants including watermelon and melon. A. citrulli strains can be divided into two major groups based on DNA fingerprint analyses and biochemical properties. Group I strains have been generally isolated from non-watermelon cucurbits, while group II strains are closely associated with watermelon. In the present study, we report the genome sequence of M6, a group I model A. citrulli strain, isolated from melon. We used comparative genome analysis to investigate differences between the genome of strain M6 and the genome of the group II model strain AAC00-1. The draft genome sequence of A. citrulli M6 harbors 139 contigs, with an overall approximate size of 4.85 Mb. The genome of M6 is ∼500 Kb shorter than that of strain AAC00-1. Comparative analysis revealed that this size difference is mainly explained by eight fragments, ranging from ∼35-120 Kb and distributed throughout the AAC00-1 genome, which are absent in the M6 genome. In agreement with this finding, while AAC00-1 was found to possess 532 open reading frames (ORFs) that are absent in strain M6, only 123 ORFs in M6 were absent in AAC00-1. Most of these M6 ORFs are hypothetical proteins and most of them were also detected in two group I strains that were recently sequenced, tw6 and pslb65. Further analyses by PCR assays and coverage analyses with other A. citrulli strains support the notion that some of these fragments or significant portions of them are discriminative between groups I and II strains of A. citrulli. Moreover, GC content, effective number of codon values and cluster of orthologs' analyses indicate that these fragments were introduced into group II strains by horizontal gene transfer events. Our study reports the genome sequence of a model group I strain of A. citrulli, one of the most important pathogens of cucurbits. It also provides the first comprehensive comparison at the genomic level between the two major groups of strains of this pathogen.

Identification of novel Xanthomonas euvesicatoria type III effector proteins by a machine-learning approach

The Gram-negative bacterium Xanthomonas euvesicatoria (Xcv) is the causal agent of bacterial spot disease in pepper and tomato. Xcv pathogenicity depends on a type III secretion (T3S) system that delivers effector proteins into host cells to suppress plant immunity and promote disease. The pool of known Xcv effectors includes approximately 30 proteins, most identified in the 85-10 strain by various experimental and computational techniques. To identify additional Xcv 85-10 effectors, we applied a genome-wide machine-learning approach, in which all open reading frames (ORFs) were scored according to their propensity to encode effectors. Scoring was based on a large set of features, including genomic organization, taxonomic dispersion, hypersensitive response and pathogenicity (hrp)-dependent expression, 5' regulatory sequences, amino acid composition bias and GC content. Thirty-six predicted effectors were tested for translocation into plant cells using the hypersensitive response (HR)-inducing domain of AvrBs2 as a reporter. Seven proteins (XopAU, XopAV, XopAW, XopAP, XopAX, XopAK and XopAD) harboured a functional translocation signal and their translocation relied on the HrpF translocon, indicating that they are bona fide T3S effectors. Remarkably, four belong to novel effector families. Inactivation of the xopAP gene reduced the severity of disease symptoms in infected plants. A decrease in cell death and chlorophyll content was observed in pepper leaves inoculated with the xopAP mutant when compared with the wild-type strain. However, populations of the xopAP mutant in infected leaves were similar in size to those of wild-type bacteria, suggesting that the reduction in virulence was not caused by impaired bacterial growth.

Computational analyses of ancient pathogen DNA from herbarium samples: challenges and prospects

The application of DNA sequencing technology to the study of ancient DNA has enabled the reconstruction of past epidemics from genomes of historically important plant-associated microbes. Recently, the genome sequences of the potato late blight pathogen Phytophthora infestans were analyzed from 19th century herbarium specimens. These herbarium samples originated from infected potatoes collected during and after the Irish potato famine. Herbaria have therefore great potential to help elucidate past epidemics of crops, date the emergence of pathogens, and inform about past pathogen population dynamics. DNA preservation in herbarium samples was unexpectedly good, raising the possibility of a whole new research area in plant and microbial genomics. However, the recovered DNA can be extremely fragmented resulting in specific challenges in reconstructing genome sequences. Here we review some of the challenges in computational analyses of ancient DNA from herbarium samples. We also applied the recently developed linkage method to haplotype reconstruction of diploid or polyploid genomes from fragmented ancient DNA.

Leptospiral pathogenomics

Leptospirosis, caused by pathogenic spirochetes belonging to the genus Leptospira, is a zoonosis with important impacts on human and animal health worldwide. Research on the mechanisms of Leptospira pathogenesis has been hindered due to slow growth of infectious strains, poor transformability, and a paucity of genetic tools. As a result of second generation sequencing technologies, there has been an acceleration of leptospiral genome sequencing efforts in the past decade, which has enabled a concomitant increase in functional genomics analyses of Leptospira pathogenesis. A pathogenomics approach, by coupling of pan-genomic analysis of multiple isolates with sequencing of experimentally attenuated highly pathogenic Leptospira, has resulted in the functional inference of virulence factors. The global Leptospira Genome Project supported by the U.S. National Institute of Allergy and Infectious Diseases to which key scientific contributions have been made from the international leptospirosis research community has provided a new roadmap for comprehensive studies of Leptospira and leptospirosis well into the future. This review describes functional genomics approaches to apply the data generated by the Leptospira Genome Project towards deepening our knowledge of virulence factors of Leptospira using the emerging discipline of pathogenomics.

Prediction of molecular mimicry candidates in human pathogenic bacteria

Molecular mimicry of host proteins is a common strategy adopted by bacterial pathogens to interfere with and exploit host processes. Despite the availability of pathogen genomes, few studies have attempted to predict virulence-associated mimicry relationships directly from genomic sequences. Here, we analyzed the proteomes of 62 pathogenic and 66 non-pathogenic bacterial species, and screened for the top pathogen-specific or pathogen-enriched sequence similarities to human proteins. The screen identified approximately 100 potential mimicry relationships including well-characterized examples among the top-scoring hits (e.g., RalF, internalin, yopH, and others), with about 1/3 of predicted relationships supported by existing literature. Examination of homology to virulence factors, statistically enriched functions, and comparison with literature indicated that the detected mimics target key host structures (e.g., extracellular matrix, ECM) and pathways (e.g., cell adhesion, lipid metabolism, and immune signaling). The top-scoring and most widespread mimicry pattern detected among pathogens consisted of elevated sequence similarities to ECM proteins including collagens and leucine-rich repeat proteins. Unexpectedly, analysis of the pathogen counterparts of these proteins revealed that they have evolved independently in different species of bacterial pathogens from separate repeat amplifications. Thus, our analysis provides evidence for two classes of mimics: complex proteins such as enzymes that have been acquired by eukaryote-to-pathogen horizontal transfer, and simpler repeat proteins that have independently evolved to mimic the host ECM. Ultimately, computational detection of pathogen-specific and pathogen-enriched similarities to host proteins provides insights into potentially novel mimicry-mediated virulence mechanisms of pathogenic bacteria.

Hunger for iron: the alternative siderophore iron scavenging systems in highly virulent Yersinia

Low molecular weight siderophores are used by many living organisms to scavenge scarcely available ferric iron. Presence of at least a single siderophore-based iron acquisition system is usually acknowledged as a virulence-associated trait and a pre-requisite to become an efficient and successful pathogen. Currently, it is assumed that yersiniabactin (Ybt) is the solely functional endogenous siderophore iron uptake system in highly virulent Yersinia (Yersinia pestis, Y. pseudotuberculosis, and Y. enterocolitica biotype 1B). Genes responsible for biosynthesis, transport, and regulation of the yersiniabactin (ybt) production are clustered on a mobile genetic element, the High-Pathogenicity Island (HPI) that is responsible for broad dissemination of the ybt genes in Enterobacteriaceae. However, the ybt gene cluster is absent from nearly half of Y. pseudotuberculosis O3 isolates and epidemic Y. pseudotuberculosis O1 isolates responsible for the Far East Scarlet-like Fever. Several potential siderophore-mediated iron uptake gene clusters are documented in Yersinia genomes, however, neither of them have been proven to be functional. It has been suggested that at least two siderophores alternative to Ybt may operate in the highly virulent Yersinia pestis/Y. pseudotuberculosis group, and are referred to as pseudochelin (Pch) and yersiniachelin (Ych). Furthermore, most sporadic Y. pseudotuberculosis O1 strains possess gene clusters encoding all three iron scavenging systems. Thus, the Ybt system appears not to be the sole endogenous siderophore iron uptake system in the highly virulent yersiniae and may be efficiently substituted and/or supplemented by alternative iron siderophore scavenging systems.