Comparative genomics and transcriptomics of trait-gene association

Identification of the active mechanism of aminoglycoside entry in V. cholerae through characterization of sRNA ctrR, regulating carbohydrate utilization and transport

The possible active entry of aminoglycosides in bacterial cells has been debated since the development of this antibiotic family. Here we report the identification of their active transport mechanism in Vibrio species. We combined genome-wide transcriptional analysis and fitness screens to identify alterations driven by treatment of V. cholerae with sub-minimum inhibitory concentrations (sub-MIC) of the aminoglycoside tobramycin. RNA-seq data showed downregulation of the small non-coding RNA ncRNA586 during such treatment, while Tn-seq revealed that inactivation of this sRNA was associated with improved fitness in the presence of tobramycin. This sRNA is located near sugar transport genes and previous work on a homologous region in Vibrio tasmaniensis suggested that this sRNA stabilizes gene transcripts for carbohydrate transport and utilization, as well as phage receptors. The role for ncRNA586, hereafter named ctrR, in the transport of both carbohydrates and aminoglycosides, was further investigated. Flow cytometry on cells treated with a fluorescent aminoglycoside confirmed the role of ctrR and of carbohydrate transporters in differential aminoglycoside entry. Despite sequence diversity, ctrR showed functional conservation across the Vibrionales. This system in directly modulated by carbon sources, suggesting regulation by carbon catabolite repression, a widely conserved mechanism in Gram-negative bacteria, priming future research on aminoglycoside uptake by sugar transporters in other bacterial species.

In vitro propagation and genome sequencing of three 'atypical' Ehrlichia ruminantium isolates

Three isolates of Ehrlichia ruminantium (Kümm 2, Omatjenne and Riverside), the causative agent of heartwater in domestic ruminants, were isolated in Ixodes scapularis (IDE8) tick cell cultures using the leukocyte fraction of infected sheep blood. All stocks were successfully propagated in IDE8 cells, whereas initiation attempts using endothelial cell cultures were unsuccessful. Therefore, the new technique should be included in any attempt to isolate field strains of E. ruminantium to enhance the probability of getting E. ruminantium isolates which might not be initiated in endothelial cells. Draft genome sequences of all three isolates were generated and compared with published genomes. The data confirmed previous phylogenetic studies that these three isolates are genetically very close to each other, but distinct from previously characterised E. ruminantium isolates. Genome comparisons indicated that the gene content and genomic synteny were highly conserved, with the exception of the membrane protein families. These findings expand our understanding of the genetic diversity of E. ruminantium and confirm the distinct phenotypic and genetic characteristics shared by these three isolates.

Macromolecular crowding links ribosomal protein gene dosage to growth rate in Vibrio cholerae

BACKGROUND

In fast-growing bacteria, the genomic location of ribosomal protein (RP) genes is biased towards the replication origin (oriC). This trait allows optimizing their expression during exponential phase since oriC neighboring regions are in higher dose due to multifork replication. Relocation of s10-spc-α locus (S10), which codes for most of the RP, to ectopic genomic positions shows that its relative distance to the oriC correlates to a reduction on its dosage, its expression, and bacterial growth rate. However, a mechanism linking S10 dosage to cell physiology has still not been determined.

RESULTS

We hypothesized that S10 dosage perturbations impact protein synthesis capacity. Strikingly, we observed that in Vibrio cholerae, protein production capacity was independent of S10 position. Deep sequencing revealed that S10 relocation altered chromosomal replication dynamics and genome-wide transcription. Such changes increased as a function of oriC-S10 distance. Since RP constitutes a large proportion of cell mass, lower S10 dosage could lead to changes in macromolecular crowding, impacting cell physiology. Accordingly, cytoplasm fluidity was higher in mutants where S10 is most distant from oriC. In hyperosmotic conditions, when crowding differences are minimized, the growth rate and replication dynamics were highly alleviated in these strains.

CONCLUSIONS

The genomic location of RP genes ensures its optimal dosage. However, besides of its essential function in translation, their genomic position sustains an optimal macromolecular crowding essential for maximizing growth. Hence, this could be another mechanism coordinating DNA replication to bacterial growth.

Dynamics of repeat-associated plasticity in the aaap gene family in Anaplasma marginale

Anaplasmosis, the most prevalent tick-transmitted disease of cattle, is caused by the rickettsial intracellular parasite Anaplasma marginale. The pathogen replicates within a parasitophorous vacuole formed from the invagination of the erythrocyte membrane. Several strains of A. marginale form "tails" or "appendages" which are attached to, and extend out from, the cytoplasmic side of the parasitophorous vacuole. Genomic analysis of the parasite antigen distributed along the appendage led to the discovery of the aaap (Anaplasma appendage associated protein) gene family located within a highly plastic region in the genome. The aaap gene family consists of aaap and several alps (for aaap-like proteins), depending on the strain. These genes/proteins are characterized by repeat sequences. To investigate locus plasticity, different versions of the locus were cloned from the same strain as well as from different strains, sequenced and aligned to identify changes. Our findings show that repeat sequences both within and between genes facilitated rearrangement events within the locus. Structural variation of the locus in the St. Maries strain was further investigated during infection of different cellular environments, i.e., bovine erythrocytes and tick cells, with a reduction in subpopulations of the aaap locus within the tick as compared to erythrocytes. Interestingly, subpopulations bearing alternative locus structures began to arise again when the pathogen was transferred from the tick environment into a naïve calf. Additionally, the Aaap protein expression profile between blood and tick samples showed a regulatory shift, indicating a host-specific response. Alignment of the protein sequences from different species of Anaplasma reveals six similar repeating motifs that appear to be unique to a few species of Anaplasma. The role the aaap locus may play in the pathogenesis of the bovine host or in tick infection/transmission remains unknown; however, the changes in aaap locus subpopulations, locus structure, and protein expression indicate that these genes have a role in strain diversification.

Impact of Three Different Mutations in Ehrlichia chaffeensis in Altering the Global Gene Expression Patterns

The rickettsial pathogen Ehrlichia chaffeensis causes a tick-borne disease, human monocytic ehrlichiosis. Mutations within certain genomic locations of the pathogen aid in understanding the pathogenesis and in developing attenuated vaccines. Our previous studies demonstrated that mutations in different genomic sites in E. chaffeensis caused variable impacts on their growth and attenuation in vertebrate and tick hosts. Here, we assessed the effect of three mutations on transcriptional changes using RNA deep-sequencing technology. RNA sequencing aided in detecting 66-80% of the transcripts of wildtype and mutant E. chaffeensis. Mutation in an antiporter gene (ECH_0379) causing attenuated growth in vertebrate hosts resulted in the down regulation of many transcribed genes. Similarly, a mutation downstream to the ECH_0490 coding sequence resulted in minimal impact on the pathogen's in vivo growth, but caused major changes in its transcriptome. This mutation caused enhanced expression of several host stress response genes. Even though the ECH_0660 gene mutation caused the pathogen's rapid clearance in vertebrate hosts and aids in generating a protective response, there was minimal impact on the transcriptome. The transcriptomic data offer novel insights about the impact of mutations on global gene expression and how they may contribute to the pathogen's resistance and/or clearance from the host.

Transcriptomic Determinants of Scrapie Prion Propagation in Cultured Ovine Microglia

Susceptibility to infection by prions is highly dependent on the amino acid sequence and host expression of the cellular prion protein (PrP^C); however, cellular expression of a genetically susceptible PrP^C is insufficient. As an example, it has been shown in cultured cells that permissive and resistant sublines derived from the same parental population often have similar expression levels of PrP^C. Thus, additional cellular factors must influence susceptibility to prion infection. The aim of this study was to elucidate the factors associated with relative permissiveness and resistance to scrapie prions in cultured cells derived from a naturally affected species. Two closely related ovine microglia clones with different prion susceptibility, but no detectable differences in PrP^C expression levels, were inoculated with either scrapie-positive or scrapie-negative sheep brainstem homogenates. Five passages post-inoculation, the transcriptional profiles of mock and infected clones were sequenced using Illumina technology. Comparative transcriptional analyses identified twenty-two differentially transcribed genes, most of which were upregulated in poorly permissive microglia. This included genes encoding for selenoprotein P, endolysosomal proteases, and proteins involved in extracellular matrix remodeling. Furthermore, in highly permissive microglia, transforming growth factor β–induced, retinoic acid receptor response 1, and phosphoserine aminotranspherase 1 gene transcripts were upregulated. Gene Set Enrichment Analysis identified proteolysis, translation, and mitosis as the most affected pathways and supported the upregulation trend of several genes encoding for intracellular proteases and ribosomal proteins in poorly permissive microglia. This study identifies new genes potentially involved in scrapie prion propagation, corroborates results from other studies, and extends those results into another cell culture model.

ESCCAR international congress on Rickettsia and other intracellular bacteria

The European Society for the study of Chlamydia, Coxiella, Anaplasma and Rickettsia (ESCCAR) held his triennial international meeting in Lausanne. This meeting gathered 165 scientists from 28 countries and all 5 continents, allowing efficient networking and major scientific exchanges. Topics covered include molecular and cellular microbiology, genomics, as well as epidemiology, veterinary and human medicine. Several breakthroughs have been revealed at the meeting, such as (i) the presence of CRISPR (the "prokaryotic immune system") in chlamydiae, (ii) an Anaplasma effector involved in host chromatin remodelling, (iii) the polarity of the type III secretion system of chlamydiae during the entry process revealed by cryo-electron tomography. Moreover, the ESCCAR meeting was a unique opportunity to be exposed to cutting-edge science and to listen to comprehensive talks on current hot topics.

Reduced Infectivity in cattle for an outer membrane protein mutant of Anaplasma marginale

Anaplasma marginale is the causative agent of anaplasmosis in cattle. Transposon mutagenesis of this pathogen using the Himar1 system resulted in the isolation of an omp10 operon insertional mutant referred to as the omp10::himar1 mutant. The work presented here evaluated if this mutant had morphological and/or growth rate defects compared to wild-type A. marginale. Results showed that the morphology, developmental cycle, and growth in tick and mammalian cell cultures are similar for the mutant and the wild type. Tick transmission experiments established that tick infection levels with the mutant were similar to those with wild-type A. marginale and that infected ticks successfully infected cattle. However, this mutant exhibited reduced infectivity and growth in cattle. The possibility of transforming A. marginale by transposon mutagenesis coupled with in vitro and in vivo assessment of altered phenotypes can aid in the identification of genes associated with virulence. The isolation of deliberately attenuated organisms that can be evaluated in their natural biological system is an important advance for the rational design of vaccines against this species.

Choice of next-generation sequencing pipelines

The next-generation sequencing (NGS) technologies are revolutionary tools which have made possible achieving remarkable advances in genetics since the beginning of the twenty-first century. Thanks to the possibility to produce large amount of sequence data, these tools are going to completely substitute other high-throughput technologies. Moreover, the large applications of NGS protocols are increasing the genetic decoding of biological systems through studies of genome anatomy and gene mapping, coupled to the transcriptome pictures. The application of NGS pipelines such as (1) de-novo genomic sequencing by mate-paired and whole-genome shotgun strategies; (2) specific gene sequencing on large bacterial communities; and (3) RNA-seq methods including whole transcriptome sequencing and Serial Analysis of Gene Expression (Sage-analysis) are fundamental in the genome-wide fields like metagenomics. Recently, the availability of these advanced protocols has allowed to overcome the usual sequencing technical issues related to the mapping specificity over standard shotgun library sequencing, the detection of large structural genomes variations and bridging sequencing gaps, as well as more precise gene annotation. In this chapter we will discuss how to manage a successful NGS pipeline from the planning of sequencing projects through the choice of the platforms up to the data analysis management.

Global transcriptional analysis reveals surface remodeling of Anaplasma marginale in the tick vector

BACKGROUND

Pathogens dependent upon vectors for transmission to new hosts undergo environment specific changes in gene transcription dependent on whether they are replicating in the vector or the mammalian host. Differential gene transcription, especially of potential vaccine candidates, is of interest in Anaplasma marginale, the tick-borne causative agent of bovine anaplasmosis.

METHODS

RNA-seq technology allowed a comprehensive analysis of the transcriptional status of A. marginale genes in two conditions: bovine host blood and tick derived cell culture, a model for the tick vector. Quantitative PCR was used to assess transcription of a set of genes in A. marginale infected tick midguts and salivary glands at two time points during the transmission cycle.

RESULTS

Genes belonging to fourteen pathways or component groups were found to be differentially transcribed in A. marginale in the bovine host versus the tick vector. One of the most significantly altered groups was composed of surface proteins. Of the 56 genes included in the surface protein group, eight were up regulated and 26 were down regulated. The down regulated surface protein encoding genes include several that are well studied due to their immunogenicity and function. Quantitative PCR of a set of genes demonstrated that transcription in tick cell culture most closely approximates transcription in salivary glands of recently infected ticks.

CONCLUSIONS

The ISE6 tick cell culture line is an acceptable model for early infection in tick salivary glands, and reveals disproportionate down regulation of surface protein genes in the tick. Transcriptional profiling in other cell lines may help us simulate additional microenvironments. Understanding vector-specific alteration of gene transcription, especially of surface protein encoding genes, may aid in the development of vaccines or transmission blocking therapies.

Knockout of an outer membrane protein operon of Anaplasma marginale by transposon mutagenesis

Background

The large amounts of data generated by genomics, transcriptomics and proteomics have increased our understanding of the biology of Anaplasma marginale. However, these data have also led to new assumptions that require testing, ideally through classical genetic mutation. One example is the definition of genes associated with virulence. Here we describe the molecular characterization of a red fluorescent and spectinomycin and streptomycin resistant A. marginale mutant generated by Himar1 transposon mutagenesis.

Results

High throughput genome sequencing to determine the Himar1-A. marginale genome junctions established that the transposon sequences were integrated within the coding region of the omp10 gene. This gene is arranged within an operon with AM1225 at the 5’ end and with omp9, omp8, omp7 and omp6 arranged in tandem at the 3’ end. RNA analysis to determine the effects of the transposon insertion on the expression of omp10 and downstream genes revealed that the Himar1 insertion not only reduced the expression of omp10 but also that of downstream genes. Transcript expression from omp9, and omp8 dropped by more than 90% in comparison with their counterparts in wild-type A. marginale. Immunoblot analysis showed a reduction in the production of Omp9 protein in these mutants compared to wild-type A. marginale.

Conclusions

These results demonstrate that transposon mutagenesis in A. marginale is possible and that this technology can be used for the creation of insertional gene knockouts that can be evaluated in natural host-vector systems.

Genetic Diversity of Tick-Borne Rickettsial Pathogens; Insights Gained from Distant Strains

The ability to capture genetic variation with unprecedented resolution improves our understanding of bacterial populations and their ability to cause disease. The goal of the pathogenomics era is to define genetic diversity that results in disease. Despite the economic losses caused by vector-borne bacteria in the Order Rickettsiales, little is known about the genetic variants responsible for observed phenotypes. The tick-transmitted rickettsial pathogen Anaplasma marginale infects cattle in tropical and subtropical regions worldwide, including Australia. Genomic analysis of North American A. marginale strains reveals a closed core genome defined by high levels of Single Nucleotide Polymorphisms (SNPs). Here we report the first genome sequences and comparative analysis for Australian strains that differ in virulence and transmissibility. A list of genetic differences that segregate with phenotype was evaluated for the ability to distinguish the attenuated strain from virulent field strains. Phylogenetic analyses of the Australian strains revealed a marked evolutionary distance from all previously sequenced strains. SNP analysis showed a strikingly reduced genetic diversity between these strains, with the smallest number of SNPs detected between any two A. marginale strains. The low diversity between these phenotypically distinct bacteria presents a unique opportunity to identify the genetic determinants of virulence and transmission.

Whole-genome sequencing in bacteriology: state of the art

Over the last ten years, genome sequencing capabilities have expanded exponentially. There have been tremendous advances in sequencing technology, DNA sample preparation, genome assembly, and data analysis. This has led to advances in a number of facets of bacterial genomics, including metagenomics, clinical medicine, bacterial archaeology, and bacterial evolution. This review examines the strengths and weaknesses of techniques in bacterial genome sequencing, upcoming technologies, and assembly techniques, as well as highlighting recent studies that highlight new applications for bacterial genomics.

Transcriptional pathways associated with the slow growth phenotype of transformed Anaplasma marginale

Background

The ability to genetically manipulate bacteria has been fundamentally important for both basic biological discovery and translational research to develop new vaccines and antibiotics. Experimental alteration of the genetic content of prokaryotic pathogens has revealed both expected functional relationships and unexpected phenotypic consequences. Slow growth phenotypes have been reported for multiple transformed bacterial species, including extracellular and intracellular pathogens. Understanding the genes and pathways responsible for the slow growth phenotype provides the opportunity to develop attenuated vaccines as well as bacteriostatic antibiotics. Transformed Anaplasma marginale, a rickettsial pathogen, exhibits slow growth in vitro and in vivo as compared to the parent wild type strain, providing the opportunity to identify the underlying genes and pathways associated with this phenotype.

Results

Whole genome transcriptional profiling allowed for identification of specific genes and pathways altered in transformed A. marginale. Genes found immediately upstream and downstream of the insertion site, including a four gene operon encoding key outer membrane proteins, were not differentially transcribed between wild type and transformed A. marginale. This lack of significant difference in transcription of flanking genes and the large size of the insert relative to the genome were consistent with a trans rather than a cis effect. Transcriptional profiling across the complete genome identified the most differentially transcribed genes, including an iron transporter, an RNA cleaving enzyme and several genes involved in translation. In order to confirm the trend seen in translation-related genes, K-means clustering and Gene Set Enrichment Analysis (GSEA) were applied. These algorithms allowed evaluation of the behavior of genes as groups that share transcriptional status or biological function. Clustering and GSEA confirmed the initial observations and found additional pathways altered in transformed A. marginale. Three pathways were significantly altered as compared to the wild type: translation, translation elongation, and purine biosynthesis.

Conclusions

Identification of perturbed genes and networks through genome wide transcriptional profiling highlights the relevance of pathways such as nucleotide biosynthesis, translation, and translation elongation in the growth phenotype of obligate intracellular bacteria. These genes and pathways provide specific targets for development of slow growing attenuated vaccines and for bacteriostatic antibiotics.