Comparative genomics of three strains of Ehrlichia ruminantium: a review

Bioinformatic Assessment of Factors Affecting the Correlation between Protein Abundance and Elongation Efficiency in Prokaryotes

Protein abundance is crucial for the majority of genetically regulated cell functions to act properly in prokaryotic organisms. Therefore, developing bioinformatic methods for assessing the efficiency of different stages of gene expression is of great importance for predicting the actual protein abundance. One of these steps is the evaluation of translation elongation efficiency based on mRNA sequence features, such as codon usage bias and mRNA secondary structure properties. In this study, we have evaluated correlation coefficients between experimentally measured protein abundance and predicted elongation efficiency characteristics for 26 prokaryotes, including non-model organisms, belonging to diverse taxonomic groups The algorithm for assessing elongation efficiency takes into account not only codon bias, but also number and energy of secondary structures in mRNA if those demonstrate an impact on predicted elongation efficiency of the ribosomal protein genes. The results show that, for a number of organisms, secondary structures are a better predictor of protein abundance than codon usage bias. The bioinformatic analysis has revealed several factors associated with the value of the correlation coefficient. The first factor is the elongation efficiency optimization type-the organisms whose genomes are optimized for codon usage only have significantly higher correlation coefficients. The second factor is taxonomical identity-bacteria that belong to the class Bacilli tend to have higher correlation coefficients among the analyzed set. The third is growth rate, which is shown to be higher for the organisms with higher correlation coefficients between protein abundance and predicted translation elongation efficiency. The obtained results can be useful for further improvement of methods for protein abundance prediction.

Intragenic codon usage in proteobacteria: Translational selection, IS expansion and genomic shrinkage

This manuscript presents a method to systematically study intragenic variations in codon usage using correspondence analysis and the effective number of codons. The method is applied to >1100 proteobacteria. Codon usage biases (measured as inertia) increases with genome size, the same is true for the percentage of inertia explained by the first axis. It is shown that there is often a relaxed or more uniform codon usage near the gene termini. Ithis is not seen n small genomes, notably those of intracellular organisms like Buchnera aphidicola or Rickettsia prowazekii where translational selection plays less of a role. When genes from E. coli, for which translational selection is well described, are split into low, intermediate and high expression, respectively, it is shown that the intragenic codon usage pattern with more uniform usage at termini exist across all three expression groups. Furthermore, the correspondence analysis reveals a unique pattern in Bordetella pertussis due to IS expansion. This study thus shows that translational selection, genome shrinkage and IS expansion result in characteristic patterns in intragenic codon usage.

Revisiting Ehrlichia ruminantium Replication Cycle Using Proteomics: The Host and the Bacterium Perspectives

The Rickettsiales Ehrlichia ruminantium, the causal agent of the fatal tick-borne disease Heartwater, induces severe damage to the vascular endothelium in ruminants. Nevertheless, E. ruminantium-induced pathobiology remains largely unknown. Our work paves the way for understanding this phenomenon by using quantitative proteomic analyses (2D-DIGE-MS/MS, 1DE-nanoLC-MS/MS and biotin-nanoUPLC-MS/MS) of host bovine aorta endothelial cells (BAE) during the in vitro bacterium intracellular replication cycle. We detect 265 bacterial proteins (including virulence factors), at all time-points of the E. ruminantium replication cycle, highlighting a dynamic bacterium–host interaction. We show that E. ruminantium infection modulates the expression of 433 host proteins: 98 being over-expressed, 161 under-expressed, 140 detected only in infected BAE cells and 34 exclusively detected in non-infected cells. Cystoscape integrated data analysis shows that these proteins lead to major changes in host cell immune responses, host cell metabolism and vesicle trafficking, with a clear involvement of inflammation-related proteins in this process. Our findings led to the first model of E. ruminantium infection in host cells in vitro, and we highlight potential biomarkers of E. ruminantium infection in endothelial cells (such as ROCK1, TMEM16K, Albumin and PTPN1), which may be important to further combat Heartwater, namely by developing non-antibiotic-based strategies.

Sweet and Sour Ehrlichia: Glycoproteomics and Phosphoproteomics Reveal New Players in Ehrlichia ruminantium Physiology and Pathogenesis

Unraveling which proteins and post-translational modifications (PTMs) affect bacterial pathogenesis and physiology in diverse environments is a tough challenge. Herein, we used mass spectrometry-based assays to study protein phosphorylation and glycosylation in Ehrlichia ruminantium Gardel virulent (ERGvir) and attenuated (ERGatt) variants and, how they can modulate Ehrlichia biological processes. The characterization of the S/T/Y phosphoproteome revealed that both strains share the same set of phosphoproteins (n = 58), 36% being overexpressed in ERGvir. The percentage of tyrosine phosphorylation is high (23%) and 66% of the identified peptides are multi-phosphorylated. Glycoproteomics revealed a high percentage of glycoproteins (67% in ERGvir) with a subset of glycoproteins being specific to ERGvir (n = 64/371) and ERGatt (n = 36/343). These glycoproteins are involved in key biological processes such as protein, amino-acid and purine biosynthesis, translation, virulence, DNA repair, and replication. Label-free quantitative analysis revealed over-expression in 31 proteins in ERGvir and 8 in ERGatt. While further PNGase digestion confidently localized 2 and 5 N-glycoproteins in ERGvir and ERGatt, respectively, western blotting suggests that many glycoproteins are O-GlcNAcylated. Twenty-three proteins were detected in both the phospho- and glycoproteome, for the two variants. This work represents the first comprehensive assessment of PTMs on Ehrlichia biology, rising interesting questions regarding ER–host interactions. Phosphoproteome characterization demonstrates an increased versatility of ER phosphoproteins to participate in different mechanisms. The high number of glycoproteins and the lack of glycosyltransferases-coding genes highlight ER dependence on the host and/or vector cellular machinery for its own protein glycosylation. Moreover, these glycoproteins could be crucial to interact and respond to changes in ER environment. PTMs crosstalk between of O-GlcNAcylation and phosphorylation could be used as a major cellular signaling mechanism in ER. As little is known about the Ehrlichia proteins/proteome and its signaling biology, the results presented herein provide a useful resource for further hypothesis-driven exploration of Ehrlichia protein regulation by phosphorylation and glycosylation events. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium with the data set identifier PXD012589.

Comparative Proteomic Profiling of Ehrlichia ruminantium Pathogenic Strain and Its High-Passaged Attenuated Strain Reveals Virulence and Attenuation-Associated Proteins

The obligate intracellular bacterium Ehrlichia ruminantium (ER) causes heartwater, a fatal tick-borne disease in livestock. In the field, ER strains present different levels of virulence, limiting vaccine efficacy, for which the molecular basis remains unknown. Moreover, there are no genetic tools currently available for ER manipulation, thus limiting the knowledge of the genes/proteins that are essential for ER pathogenesis and biology. As such, to identify proteins and/or mechanisms involved in ER virulence, we performed the first exhaustive comparative proteomic analysis between a virulent strain (ERGvir) and its high-passaged attenuated strain (ERGatt). Despite their different behaviors in vivo and in vitro, our results from 1DE-nanoLC-MS/MS showed that ERGvir and ERGatt share 80% of their proteins; this core proteome includes chaperones, proteins involved in metabolism, protein-DNA-RNA biosynthesis and processing, and bacterial effectors. Conventional 2DE revealed that 85% of the identified proteins are proteoforms, suggesting that post-translational modifications (namely glycosylation) are important in ER biology. Strain-specific proteins were also identified: while ERGatt has an increased number and overexpression of proteins involved in cell division, metabolism, transport and protein processing, ERGvir shows an overexpression of proteins and proteoforms (DIGE experiments) involved in pathogenesis such as Lpd, AnkA, VirB9 and B10, providing molecular evidence for its increased virulence in vivo and in vitro. Overall, our work reveals that ERGvir and ERGatt proteomes are streamlined to fulfill their biological function (maximum virulence for ERGvir and replicative capacity for ERGatt), and we provide both pioneering data and novel insights into the pathogenesis of this obligate intracellular bacterium.

Ehrlichia's molecular tricks to manipulate their host cells

Ehrlichia is a large genus of obligate intracellular Gram-negative bacteria transmitted by ticks that cause several emerging infectious diseases in humans and are pathogenic on rodents, ruminants, and dogs. Ehrlichia spp. invade and replicate either in endothelial cells, white blood cells, or within midgut cells and salivary glands of their vector ticks. In this review, we discuss the insights that functional studies are providing on how this group of bacteria exploits their host by subverting host innate immunity and hijacking cellular processes.

Comparative genomics of first available bovine Anaplasma phagocytophilum genome obtained with targeted sequence capture

BACKGROUND

Anaplasma phagocytophilum is a zoonotic and obligate intracellular bacterium transmitted by ticks. In domestic ruminants, it is the causative agent of tick-borne fever, which causes significant economic losses in Europe. As A. phagocytophilum is difficult to isolate and cultivate, only nine genome sequences have been published to date, none of which originate from a bovine strain.Our goals were to; 1/ develop a sequencing methodology which efficiently circumvents the difficulties associated with A. phagocytophilum isolation and culture; 2/ describe the first genome of a bovine strain; and 3/ compare it with available genomes, in order to both explore key genomic features at the species level, and to identify candidate genes that could be specific to bovine strains.

RESULTS

DNA was extracted from a bovine blood sample infected by A. phagocytophilum. Following a whole genome capture approach, A. phagocytophilum DNA was enriched 197-fold in the sample and then sequenced using Illumina technology. In total, 58.9% of obtained reads corresponded to the A. phagocytophilum genome, covering 85.3% of the HZ genome. Then by performing comparisons with nine previously-sequenced A. phagocytophilum genomes, we determined the core genome of these ten strains. Following analysis, 1281 coding DNA sequences, including 1001 complete sequences, were detected in the A. phagocytophilum bovine genome, of which four appeared to be unique to the bovine isolate. These four coding DNA sequences coded for "hypothetical proteins of unknown function" and require further analysis. We also identified nine proteins common to both European domestic ruminants tested.

CONCLUSION

Using a whole genome capture approach, we have sequenced the first A. phagocytophilum genome isolated from a cow. To the best of our knowledge, this is the first time that this method has been used to selectively enrich pathogenic bacterial DNA from samples also containing host DNA. The four proteins unique to the A. phagocytophilum bovine genome could be involved in host tropism, therefore their functions need to be explored.

Progress and obstacles in vaccine development for the ehrlichioses

Ehrlichia are tick-borne obligately intracellular bacteria that cause significant diseases in veterinary natural hosts, including livestock and companion animals, and are now considered important zoonotic pathogens in humans. Vaccines are needed for these veterinary and zoonotic human pathogens, but many obstacles exist that have impeded their development. These obstacles include understanding genetic and antigenic variability, influence of the host on the pathogen phenotype and immunogenicity, identification of the ehrlichial antigens that stimulate protective immunity and those that elicit immunopathology, development of animal models that faithfully reflect the immune responses of the hosts and understanding molecular host-pathogen interactions involved in immune evasion or that may be blocked by the host immune response. We review the obstacles and progress in addressing barriers associated with vaccine development to protect livestock, companion animals and humans against these host defense-evasive and cell function-manipulative, vector-transmitted pathogens.

Genomic acquisition of a capsular polysaccharide virulence cluster by non-pathogenic Burkholderia isolates

Background

Burkholderia thailandensis is a non-pathogenic environmental saprophyte closely related to Burkholderia pseudomallei, the causative agent of the often fatal animal and human disease melioidosis. To study B. thailandensis genomic variation, we profiled 50 isolates using a pan-genome microarray comprising genomic elements from 28 Burkholderia strains and species.

Results

Of 39 genomic regions variably present across the B. thailandensis strains, 13 regions corresponded to known genomic islands, while 26 regions were novel. Variant B. thailandensis isolates exhibited isolated acquisition of a capsular polysaccharide biosynthesis gene cluster (B. pseudomallei-like capsular polysaccharide) closely resembling a similar cluster in B. pseudomallei that is essential for virulence in mammals; presence of this cluster was confirmed by whole genome sequencing of a representative variant strain (B. thailandensis E555). Both whole-genome microarray and multi-locus sequence typing analysis revealed that the variant strains formed part of a phylogenetic subgroup distinct from the ancestral B. thailandensis population and were associated with atypical isolation sources when compared to the majority of previously described B. thailandensis strains. In functional assays, B. thailandensis E555 exhibited several B. pseudomallei-like phenotypes, including colony wrinkling, resistance to human complement binding, and intracellular macrophage survival. However, in murine infection assays, B. thailandensis E555 did not exhibit enhanced virulence relative to other B. thailandensis strains, suggesting that additional factors are required to successfully colonize and infect mammals.

Conclusions

The discovery of such novel variant strains demonstrates how unbiased genomic surveys of non-pathogenic isolates can reveal insights into the development and emergence of new pathogenic species.

Innovative approach for transcriptomic analysis of obligate intracellular pathogen: selective capture of transcribed sequences of Ehrlichia ruminantium

Background

Whole genome transcriptomic analysis is a powerful approach to elucidate the molecular mechanisms controlling the pathogenesis of obligate intracellular bacteria. However, the major hurdle resides in the low quantity of prokaryotic mRNAs extracted from host cells. Our model Ehrlichia ruminantium (ER), the causative agent of heartwater, is transmitted by tick Amblyomma variegatum. This bacterium affects wild and domestic ruminants and is present in Sub-Saharan Africa and the Caribbean islands. Because of its strictly intracellular location, which constitutes a limitation for its extensive study, the molecular mechanisms involved in its pathogenicity are still poorly understood.

Results

We successfully adapted the SCOTS method (Selective Capture of Transcribed Sequences) on the model Rickettsiales ER to capture mRNAs. Southern Blots and RT-PCR revealed an enrichment of ER's cDNAs and a diminution of ribosomal contaminants after three rounds of capture. qRT-PCR and whole-genome ER microarrays hybridizations demonstrated that SCOTS method introduced only a limited bias on gene expression. Indeed, we confirmed the differential gene expression between poorly and highly expressed genes before and after SCOTS captures. The comparative gene expression obtained from ER microarrays data, on samples before and after SCOTS at 96 hpi was significantly correlated (R² = 0.7). Moreover, SCOTS method is crucial for microarrays analysis of ER, especially for early time points post-infection. There was low detection of transcripts for untreated samples whereas 24% and 70.7% were revealed for SCOTS samples at 24 and 96 hpi respectively.

Conclusions

We conclude that this SCOTS method has a key importance for the transcriptomic analysis of ER and can be potentially used for other Rickettsiales. This study constitutes the first step for further gene expression analyses that will lead to a better understanding of both ER pathogenicity and the adaptation of obligate intracellular bacteria to their environment.