The active site residue V266 of Chlamydial HtrA is critical for substrate binding during both in vitro and in vivo conditions

HtrA is a complex, multimeric chaperone and serine protease important for the virulence and survival of many bacteria. Chlamydia trachomatis is an obligate, intracellular bacterial pathogen that is responsible for severe disease pathology. C. trachomatis HtrA (CtHtrA) has been shown to be highly expressed in laboratory models of disease. In this study, molecular modelling of CtHtrA protein active site structure identified putative S1-S3 subsite residues I242, I265, and V266. These residues were altered by site-directed mutagenesis, and these changes were shown to considerably reduce protease activity on known substrates and resulted in a narrower and distinct range of substrates compared to wild type. Bacterial two-hybrid analysis revealed that CtHtrA is able to interact in vivo with a broad range of protein sequences with high affinity. Notably, however, the interaction was significantly altered in 35 out of 69 clones when residue V266 was mutated, indicating that this residue has an important function during substrate binding.

Prediction of transcriptional regulatory interactions in bacteria : a comparative genomics approach

Exponential growth of genomic data in the last two decades has made manual analyses impractical for all but trial studies. As genomic analyses have become more sophisticated, and move toward comparisons across large datasets, computational approaches have become essential. One of the most important biological questions is to understand the mechanisms underlying gene regulation. Genetic regulation is commonly investigated and modelled through the use of transcriptional regulatory network (TRN) structures. These model the regulatory interactions between two key components: transcription factors (TFs) and the target genes (TGs) they regulate. Transcriptional regulatory networks have proven to be invaluable scientific tools in Bioinformatics. When used in conjunction with comparative genomics, they have provided substantial insights into the evolution of regulatory interactions. Current approaches to regulatory network inference, however, omit two additional key entities: promoters and transcription factor binding sites (TFBSs). In this study, we attempted to explore the relationships among these regulatory components in bacteria. Our primary goal was to identify relationships that can assist in reducing the high false positive rates associated with transcription factor binding site predictions and thereupon enhance the reliability of the inferred transcription regulatory networks. In our preliminary exploration of relationships between the key regulatory components in Escherichia coli transcription, we discovered a number of potentially useful features. The combination of location score and sequence dissimilarity scores increased de novo binding site prediction accuracy by 13.6%. Another important observation made was with regards to the relationship between transcription factors grouped by their regulatory role and corresponding promoter strength. Our study of E.coli ��70 promoters, found support at the 0.1 significance level for our hypothesis | that weak promoters are preferentially associated with activator binding sites to enhance gene expression, whilst strong promoters have more repressor binding sites to repress or inhibit gene transcription. Although the observations were specific to �70, they nevertheless strongly encourage additional investigations when more experimentally confirmed data are available. In our preliminary exploration of relationships between the key regulatory components in E.coli transcription, we discovered a number of potentially useful features { some of which proved successful in reducing the number of false positives when applied to re-evaluate binding site predictions. Of chief interest was the relationship observed between promoter strength and TFs with respect to their regulatory role. Based on the common assumption, where promoter homology positively correlates with transcription rate, we hypothesised that weak promoters would have more transcription factors that enhance gene expression, whilst strong promoters would have more repressor binding sites. The t-tests assessed for E.coli �70 promoters returned a p-value of 0.072, which at 0.1 significance level suggested support for our (alternative) hypothesis; albeit this trend may only be present for promoters where corresponding TFBSs are either all repressors or all activators. Nevertheless, such suggestive results strongly encourage additional investigations when more experimentally confirmed data will become available. Much of the remainder of the thesis concerns a machine learning study of binding site prediction, using the SVM and kernel methods, principally the spectrum kernel. Spectrum kernels have been successfully applied in previous studies of protein classification [91, 92], as well as the related problem of promoter predictions [59], and we have here successfully applied the technique to refining TFBS predictions. The advantages provided by the SVM classifier were best seen in `moderately'-conserved transcription factor binding sites as represented by our E.coli CRP case study. Inclusion of additional position feature attributes further increased accuracy by 9.1% but more notable was the considerable decrease in false positive rate from 0.8 to 0.5 while retaining 0.9 sensitivity. Improved prediction of transcription factor binding sites is in turn extremely valuable in improving inference of regulatory relationships, a problem notoriously prone to false positive predictions. Here, the number of false regulatory interactions inferred using the conventional two-component model was substantially reduced when we integrated de novo transcription factor binding site predictions as an additional criterion for acceptance in a case study of inference in the Fur regulon. This initial work was extended to a comparative study of the iron regulatory system across 20 Yersinia strains. This work revealed interesting, strain-specific difierences, especially between pathogenic and non-pathogenic strains. Such difierences were made clear through interactive visualisations using the TRNDifi software developed as part of this work, and would have remained undetected using conventional methods. This approach led to the nomination of the Yfe iron-uptake system as a candidate for further wet-lab experimentation due to its potential active functionality in non-pathogens and its known participation in full virulence of the bubonic plague strain. Building on this work, we introduced novel structures we have labelled as `regulatory trees', inspired by the phylogenetic tree concept. Instead of using gene or protein sequence similarity, the regulatory trees were constructed based on the number of similar regulatory interactions. While the common phylogentic trees convey information regarding changes in gene repertoire, which we might regard being analogous to `hardware', the regulatory tree informs us of the changes in regulatory circuitry, in some respects analogous to `software'. In this context, we explored the `pan-regulatory network' for the Fur system, the entire set of regulatory interactions found for the Fur transcription factor across a group of genomes. In the pan-regulatory network, emphasis is placed on how the regulatory network for each target genome is inferred from multiple sources instead of a single source, as is the common approach. The benefit of using multiple reference networks, is a more comprehensive survey of the relationships, and increased confidence in the regulatory interactions predicted. In the present study, we distinguish between relationships found across the full set of genomes as the `core-regulatory-set', and interactions found only in a subset of genomes explored as the `sub-regulatory-set'. We found nine Fur target gene clusters present across the four genomes studied, this core set potentially identifying basic regulatory processes essential for survival. Species level difierences are seen at the sub-regulatory-set level; for example the known virulence factors, YbtA and PchR were found in Y.pestis and P.aerguinosa respectively, but were not present in both E.coli and B.subtilis. Such factors and the iron-uptake systems they regulate, are ideal candidates for wet-lab investigation to determine whether or not they are pathogenic specific. In this study, we employed a broad range of approaches to address our goals and assessed these methods using the Fur regulon as our initial case study. We identified a set of promising feature attributes; demonstrated their success in increasing transcription factor binding site prediction specificity while retaining sensitivity, and showed the importance of binding site predictions in enhancing the reliability of regulatory interaction inferences. Most importantly, these outcomes led to the introduction of a range of visualisations and techniques, which are applicable across the entire bacterial spectrum and can be utilised in studies beyond the understanding of transcriptional regulatory networks.

Intra-amniotic ureaplasma infection : adverse outcomes vary depending on gestation

Background: Ureaplasmas are the most prevalent bacteria isolated from preterm deliveries and the prognosis for neonates varies depending on the gestation at delivery. Ureaplasmas vary their surface-exposed antigen (MBA, a virulence mechanism) during chronic intra-amniotic infections, but it is not known when changes first occur during gestation. Method: U. parvum serovar 3 (2x10e7CFU) was injected intra-amniotically (IA) into six experimental cohorts of pregnant ewes (of n=7), 3 days (d) or 7d before delivery at either: 100d, 124d or 140d gestation (term=145d). Control ewes received IA 10B broth. Fetuses were delivered surgically and ureaplasmas cultured from amniotic fluid (AF), chorioamnion, fetal lung (FL) and umbilical cord. Ureaplasmas were tested by western blot to demonstrate MBA variation. Results: The highest number of ureaplasmas were recovered from FL at 100d gestation after 3 days of infection (p<0.03). Six of 7(86%) 100d–3d FL demonstrated an ureaplasma MBA variant, but only 17% and 15% of FL showed an MBA variant after 3d infection at 124d and 140d gestation respectively. Greatest variation of the MBA occurred in AF and FL at 124d gestation after 7d infection. The least MBA variation was observed at 140d; however, at this time the most severe histological chorioamnionitis was observed. Conclusions: After intra-amniotic ureaplasma injections, higher numbers of ureaplasmas gained access to the FL at 100d gestation than observed at later gestations. This may exacerbate the adverse outcomes for neonates delivered early in gestation. In late gestation, ureaplasma MBA variation was minimal, but chorioamnionitis was the most severe. Adverse pregnancy outcomes associated with IA ureaplasma infection may vary depending on the duration of gestation, the number of ureaplasmas isolated from the fetal tissues and the degree of MBA variation.

The role of HtrA as a chaperone and protease in bacterial pathogenesis

HtrA (High Temperature Requirement A) is a critical stress response protease and chaperone for many bacteria. HtrA is a multitasking protein which can degrade unfolded proteins, conduct specific proteolysis of some substrates for correct assembly, interact with substrates to ensure correct folding, assembly or localisation, and chaperone unfolded proteins. These functions are critical for the virulence of a number of bacterial pathogens, in some cases not simply due to the broad activities of HtrA in protection against the protein stress conditions which occur during virulence. But also due to the role of HtrA in either specific proteolysis or assembly of key protein substrates which function directly in virulence. Remarkably, these activities are all conducted without any requirement for ATP. The biochemical mechanism of HtrA relies both on the chymotryptic serine protease active site as well as the presence of two PDZ (protein binding) domains. The mechanism is a unique combination of activation by substrate motifs to alter the confirmation of the active site, and assembly into a multimeric complex which has enhanced degradation and may also act as a protective cage for proteins which are not degraded. The role of this protease in the pathogenesis of a number of bacteria and the details of its distinctive biochemical activation and assembly mechanisms are discussed in this chapter.

The role of HtrA as a chaperone and protease in bacterial pathogenisis

HtrA (High Temperature Requirement A) is a critical stress response protease and chaperone for many bacteria. HtrA is a multitasking protein which can degrade unfolded proteins, conduct specific proteolysis of some substrates for correct assembly, interact with substrates to ensure correct folding, assembly or localisation, and chaperone unfolded proteins. These functions are critical for the virulence of a number of bacterial pathogens, in some cases not simply due to the broad activities of HtrA in protection against the protein stress conditions which occur during virulence. But also due to the role of HtrA in either specific proteolysis or assembly of key protein substrates which function directly in virulence. Remarkably, these activities are all conducted without any requirement for ATP. The biochemical mechanism of HtrA relies both on the chymotryptic serine protease active site as well as the presence of two PDZ (protein binding) domains. The mechanism is a unique combination of activation by substrate motifs to alter the confirmation of the active site, and assembly into a multimeric complex which has enhanced degradation and may also act as a protective cage for proteins which are not degraded. The role of this protease in the pathogenesis of a number of bacteria and the details of its distinctive biochemical activation and assembly mechanisms are discussed in this chapter.

Ureaplasma parvum multiple banded antigen (MBA) size variation - association with fetal inflammation in a sheep model (Published abstract)

Background: Ureaplasma species are the most prevalent isolates from women who deliver preterm. The MBA, a surface exposed lipoprotein, is a key virulence factor of ureaplasmas. We investigated MBA variation after chronic and acute intra-amniotic (IA) ureaplasma infections. Method: U. parvum serovar 3 (2x104 colony-forming-units) was injected IA into pregnant ewes at: 55 days gestation (d, term = 145d) (n=8); 117d (n=8) and 121d (n=8). Fetuses were delivered surgically (124d) and ureaplasmas cultured from amniotic fluid (AF), chorioamnion, fetal lung (FL) and umbilical cord were tested by western blot and PCR assays to demonstrate MBA and mba gene variation respectively. Tissue sections were sectioned and stained by haemotoxylin and eosin and inflammatory cell counts and pathology were reported (blinded to outcome). Results: Numerous MBA/mba variants were generated in vivo after chronic exposure to ureaplasma infection but after acute infection no variants (3d) or very few variants (7d) were generated. Identical MBA variants were detected within the AF and FL but different ureaplasma variants were detected within chorioamnion specimens. The severity of inflammation within chronically infected tissues varied between animals ranging from no inflammation to severe inflammation with/without fibrosis. Chorioamnion, FL and cord from the same animal demonstrated the same degree of inflammation. Conclusions: MBA/mba variation in vivo occurred after the initiation of the host immune response and we propose that ureaplasmas vary the MBA antigen to evade the host immune response. In some animals there was no inflammation despite colonisation with high numbers of ureaplasmas.

Ureaplasma parvum serovar 3 multiple banded antigen size variation after chronic intra-amniotic infection/colonization

Ureaplasma species are the microorganisms most frequently associated with adverse pregnancy outcomes. The multiple banded antigen (MBA), a surface-exposed lipoprotein, is a key virulence factor of ureaplasmas. The MBA demonstrates size variation, which we have shown previously to be correlated with the severity of chorioamnion inflammation. We aimed to investigate U. parvum serovar 3 pathogenesis in vivo, using a sheep model, by investigating: MBA variation after long term (chronic) and short term (acute) durations of in utero ureaplasma infections, and the severity of chorioamnionitis and inflammation in other fetal tissues. Inocula of 2x107 colony-forming-units (CFU) of U. parvum serovar 3 (Up) or media controls (C) were injected intra-amniotically into pregnant ewes at one of three time points: day 55 (69d Up, n=8; C69, n=4); day 117 (7d Up, n=8; C7, n=2); and day 121 (3d Up, n=8; C3, n=2) of gestation (term=145-150d). At day 124, preterm fetuses were delivered surgically. Samples of chorioamnion, fetal lung, and umbilical cord were: (i) snap frozen for subsequent ureaplasma culture, and (ii) fixed, embedded, sectioned and stained by haematoxylin and eosin stain for histological analysis. Selected fetal lung clinical ureaplasma isolates were cloned and filtered to obtain cultures from a single CFU. Passage 1 and clone 2 ureaplasma cultures were tested by western blot to demonstrate MBA variation. In acute durations of ureaplasma infection no MBA variants (3d Up) or very few MBA variants (7d Up) were present when compared to the original inoculum. However, numerous MBA size variants were generated in vivo (alike within contiguous tissues, amniotic fluid and fetal lung, but different variants were present within chorioamnion), during chronic, 69d exposure to ureaplasma infection. For the first time we have shown that the degree of ureaplasma MBA variation in vivo increased with the duration of gestation.

Phylodynamic analysis of the emergence and epidemiological impact of transmissible defective dengue viruses

Intra-host sequence data from RNA viruses have revealed the ubiquity of defective viruses in natural viral populations, sometimes at surprisingly high frequency. Although defective viruses have long been known to laboratory virologists, their relevance in clinical and epidemiological settings has not been established. The discovery of long-term transmission of a defective lineage of dengue virus type 1 (DENV-1) in Myanmar, first seen in 2001, raised important questions about the emergence of transmissible defective viruses and their role in viral epidemiology. By combining phylogenetic analyses and dynamical modelling, we investigate how evolutionary and ecological processes at the intra-host and inter-host scales shaped the emergence and spread of the defective DENV-1 lineage. We show that this lineage of defective viruses emerged between June 1998 and February 2001, and that the defective virus was transmitted primarily through co-transmission with the functional virus to uninfected individuals. We provide evidence that, surprisingly, this co-transmission route has a higher transmission potential than transmission of functional dengue viruses alone. Consequently, we predict that the defective lineage should increase overall incidence of dengue infection, which could account for the historically high dengue incidence reported in Myanmar in 2001-2002. Our results show the unappreciated potential for defective viruses to impact the epidemiology of human pathogens, possibly by modifying the virulence-transmissibility trade-off, or to emerge as circulating infections in their own right. They also demonstrate that interactions between viral variants, such as complementation, can open new pathways to viral emergence.

Yersinia high pathogenicity island genes modify the Escherichia coli primary metabolome independently of siderophore production

Bacterial siderophores may enhance pathogenicity by scavenging iron, but their expression has been proposed to exert a substantial metabolic cost. Here we describe a combined metabolomic-genetic approach to determine how mutations affecting the virulence-associated siderophore yersiniabactin affect the Escherichia coli primary metabolome. Contrary to expectations, we did not find yersiniabactin biosynthesis to correspond to consistent metabolomic shifts. Instead, we found that targeted deletion of ybtU or ybtA, dissimilar genes with similar roles in regulating yersiniabactin expression, were associated with a specific shift in arginine pathway metabolites during growth in minimal media. This interaction was associated with high arginine levels in the model uropathogen Escherichia coli UTI89 compared to its ybtU and ybtA mutants and the K12 strain MG1655, which lacks yersiniabactin-associated genes. Because arginine is not a direct yersiniabactin biosynthetic substrate, these findings show that virulence-associated secondary metabolite systems may shape bacterial primary metabolism independently of substrate consumption

Development status and future prospects for a vaccine against Chlamydia trachomatis infection

Chlamydia trachomatis continues to be the most commonly reported sexually transmitted bacterial infection in many countries with more than 100 million new cases estimated annually. These acute infections translate into significant downstream health care costs, particularly for women, where complications can include pelvic inflammatory disease and other disease sequelae such as tubal factor infertility. Despite years of research, the immunological mechanisms responsible for protective immunity versus immunopathology are still not well understood, although it is widely accepted that T cell driven IFN-g and Th17 responses are critical for clearing infection. While antibodies are able to neutralize infections in vitro, alone they are not protective, indicating that any successful vaccine will need to elicit both arms of the immune response. In recent years, there has been an expansion in the number and types of antigens that have been evaluated as vaccines, and combined with the new array of mucosal adjuvants, this aspect of chlamydial vaccinology is showing promise. Most recently, the opportunities to develop successful vaccines have been given a significant boost with the development of a genetic transformation system for Chlamydia, as well as the identification of the key role of the chlamydial plasmid in virulence. While still remaining a major challenge, the development of a successful C.trachomatis vaccine is starting to look more likely.

Proteolytic activation of Chlamydia trachomatis HTRA is mediated by PDZ1 domain interactions with protease domain loops L3 and LC and beta strand β5

Chlamydia trachomatis is a bacterial pathogen responsible for one of the most prevalent sexually transmitted infections worldwide. Its unique development cycle has limited our understanding of its pathogenic mechanisms. However, CtHtrA has recently been identified as a potential C. trachomatis virulence factor. CtHtrA is a tightly regulated quality control protein with a monomeric structural unit comprised of a chymotrypsin-like protease domain and two PDZ domains. Activation of proteolytic activity relies on the C-terminus of the substrate allosterically binding to the PDZ1 domain, which triggers subsequent conformational change and oligomerization of the protein into 24-mers enabling proteolysis. This activation is mediated by a cascade of precise structural arrangements, but the specific CtHtrA residues and structural elements required to facilitate activation are unknown. Using in vitro analysis guided by homology modeling, we show that the mutation of residues Arg362 and Arg224, predicted to disrupt the interaction between the CtHtrA PDZ1 domain and loop L3, and between loop L3 and loop LD, respectively, are critical for the activation of proteolytic activity. We also demonstrate that mutation to residues Arg299 and Lys160, predicted to disrupt PDZ1 domain interactions with protease loop LC and strand β5, are also able to influence proteolysis, implying their involvement in the CtHtrA mechanism of activation. This is the first investigation of protease loop LC and strand β5 with respect to their potential interactions with the PDZ1 domain. Given their high level of conservation in bacterial HtrA, these structural elements may be equally significant in the activation mechanism of DegP and other HtrA family members.

On the role of RNA silencing in the pathogenicity and evolution of viroids and viral satellites

Viroids and most viral satellites have small, noncoding, and highly structured RNA genomes. How they cause disease symptoms without encoding proteins and why they have characteristic secondary structures are two longstanding questions. Recent studies have shown that both viroids and satellites are capable of inducing RNA silencing, suggesting a possible role of this mechanism in the pathology and evolution of these subviral RNAs. Here we show that preventing RNA silencing in tobacco, using a silencing suppressor, greatly reduces the symptoms caused by the Y satellite of cucumber mosaic virus. Furthermore, tomato plants expressing hairpin RNA, derived from potato spindle tuber viroid, developed symptoms similar to those of potato spindle tuber viroid infection. These results provide evidence suggesting that viroids and satellites cause disease symptoms by directing RNA silencing against physiologically important host genes. We also show that viroid and satellite RNAs are significantly resistant to RNA silencing-mediated degradation, suggesting that RNA silencing is an important selection pressure shaping the evolution of the secondary structures of these pathogens.

Metabolomic analysis of siderophore cheater mutants reveals metabolic costs of expression in uropathogenic escherichia coli

Bacterial siderophores are a group of chemically diverse, virulence-associated secondary metabolites whose expression exerts metabolic costs. A combined bacterial genetic and metabolomic approach revealed differential metabolomic impacts associated with biosynthesis of different siderophore structural families. Despite myriad genetic differences, the metabolome of a cheater mutant lacking a single set of siderophore biosynthetic genes more closely approximate that of a nonpathogenic K12 strain than its isogenic, uropathogen parent strain. Siderophore types associated with greater metabolomic perturbations are less common among human isolates, suggesting that metabolic costs influence success in a human population. Although different siderophores share a common iron acquisition function, our analysis shows how a metabolomic approach can distinguish their relative metabolic impacts in E.coli.

The infectivity and host range of Orgyia anartoides nucleopolyhedrovirus

The painted apple moth (PAM), Teia anartoides (Walker) (Lepidoptera: Lymantriidae) made a recent incursion into New Zealand. A nucleopolyhedrovirus (NPV), Orgyia anartoides NPV (OranNPV), originally isolated from PAM in Australia, was tested for its pathogenicity to PAM and a range of non-target insect species found in New Zealand, to evaluate its suitability as a microbial control for this insect invader. Dosage-mortality tests showed that OranNPV was highly pathogenic to PAM larvae; mean LT50 values for third instars ranged from 17.9 to 8.1 days for doses from 102 to 105 polyhedral inclusion bodies/larva, respectively. The cause of death in infected insects was confirmed as OranNPV. Molecular analysis established that OranNPV can be identified by PCR and restriction digestion, and this process complemented microscopic examination of infected larvae. No lymantriid species occur in New Zealand; however, the virus had no significant effects on species from five other lepidopteran families (Noctuidae, Tortricidae, Geometridae, Nymphalidae and Plutellidae) or on adult honeybees. Thus, all indications from this initial investigation are that OranNPV would be an important tool in the control of PAM in a future incursion of this species into New Zealand.

Global dissemination of a multidrug resistant Escherichia coli clone

Escherichia coli sequence type 131 (ST131) is a globally disseminated, multidrug resistant (MDR) clone responsible for a high proportion of urinary tract and bloodstream infections. The rapid emergence and successful spread of E. coli ST131 is strongly associated with several factors, including resistance to fluoroquinolones, high virulence gene content, the possession of the type 1 fimbriae FimH30 allele, and the production of the CTX-M-15 extended spectrum β-lactamase (ESBL). Here, we used genome sequencing to examine the molecular epidemiology of a collection of E. coli ST131 strains isolated from six distinct geographical locations across the world spanning 2000–2011. The global phylogeny of E. coli ST131, determined from whole-genome sequence data, revealed a single lineage of E. coli ST131 distinct from other extraintestinal E. coli strains within the B2 phylogroup. Three closely related E. coli ST131 sublineages were identified, with little association to geographic origin. The majority of single-nucleotide variants associated with each of the sublineages were due to recombination in regions adjacent to mobile genetic elements (MGEs). The most prevalent sublineage of ST131 strains was characterized by fluoroquinolone resistance, and a distinct virulence factor and MGE profile. Four different variants of the CTX-M ESBL–resistance gene were identified in our ST131 strains, with acquisition of CTX-M-15 representing a defining feature of a discrete but geographically dispersed ST131 sublineage. This study confirms the global dispersal of a single E. coli ST131 clone and demonstrates the role of MGEs and recombination in the evolution of this important MDR pathogen.