Phylogeography of Francisella tularensis subspecies holarctica from the country of Georgia

<p>Abstract</p> <p>Background</p> <p><it>Francisella tularensis</it>, the causative agent of tularemia, displays subspecies-specific differences in virulence, geographic distribution, and genetic diversity. <it>F. tularensis </it>subsp. <it>holarctica </it>is widely distributed throughout the Northern Hemisphere. In Europe, <it>F. tularensis </it>subsp. <it>holarctica </it>isolates have largely been assigned to two phylogenetic groups that have specific geographic distributions. Most isolates from Western Europe are assigned to the B.Br.FTNF002-00 group, whereas most isolates from Eastern Europe are assigned to numerous lineages within the B.Br.013 group. The eastern geographic extent of the B.Br.013 group is currently unknown due to a lack of phylogenetic knowledge about populations at the European/Asian juncture and in Asia. In this study, we address this knowledge gap by describing the phylogenetic structure of <it>F. tularensis </it>subsp. <it>holarctica </it>isolates from the country of Georgia, and by placing these isolates into a global phylogeographic context.</p> <p>Results</p> <p>We identified a new genetic lineage of <it>F. tularensis </it>subsp. <it>holarctica </it>from Georgia that belongs to the B.Br.013 group. This new lineage is genetically and geographically distinct from lineages previously described from the B.Br.013 group from Central-Eastern Europe. Importantly, this new lineage is basal within the B.Br.013 group, indicating the Georgian lineage diverged before the diversification of the other known B.Br.013 lineages. Although two isolates from the Georgian lineage were collected nearby in the Ukrainian region of Crimea, all other global isolates assigned to this lineage were collected in Georgia. This restricted geographic distribution, as well as the high levels of genetic diversity within the lineage, is consistent with a relatively older origin and localized differentiation.</p> <p>Conclusions</p> <p>We identified a new lineage of <it>F. tularensis </it>subsp. <it>holarctica </it>from Georgia that appears to have an older origin than any other diversified lineages previously described from the B.Br.013 group. This finding suggests that additional phylogenetic studies of <it>F. tularensis </it>subsp. <it>holarctica </it>populations in Eastern Europe and Asia have the potential to yield important new insights into the evolutionary history and phylogeography of this broadly dispersed <it>F. tularensis </it>subspecies.</p

Host-Pathogen O-Methyltransferase Similarity and Its Specific Presence in Highly Virulent Strains of Francisella tularensis Suggests Molecular Mimicry

Whole genome comparative studies of many bacterial pathogens have shown an overall high similarity of gene content (>95%) between phylogenetically distinct subspecies. In highly clonal species that share the bulk of their genomes subtle changes in gene content and small-scale polymorphisms, especially those that may alter gene expression and protein-protein interactions, are more likely to have a significant effect on the pathogen's biology. In order to better understand molecular attributes that may mediate the adaptation of virulence in infectious bacteria, a comparative study was done to further analyze the evolution of a gene encoding an o-methyltransferase that was previously identified as a candidate virulence factor due to its conservation specifically in highly pathogenic Francisella tularensis subsp. tularensis strains. The o-methyltransferase gene is located in the genomic neighborhood of a known pathogenicity island and predicted site of rearrangement. Distinct o-methyltransferase subtypes are present in different Francisella tularensis subspecies. Related protein families were identified in several host species as well as species of pathogenic bacteria that are otherwise very distant phylogenetically from Francisella, including species of Mycobacterium. A conserved sequence motif profile is present in the mammalian host and pathogen protein sequences, and sites of non-synonymous variation conserved in Francisella subspecies specific o-methyltransferases map proximally to the predicted active site of the orthologous human protein structure. Altogether, evidence suggests a role of the F. t. subsp. tularensis protein in a mechanism of molecular mimicry, similar perhaps to Legionella and Coxiella. These findings therefore provide insights into the evolution of niche-restriction and virulence in Francisella, and have broader implications regarding the molecular mechanisms that mediate host-pathogen relationships

Comparative Genomic Characterization of Francisella tularensis Strains Belonging to Low and High Virulence Subspecies

Tularemia is a geographically widespread, severely debilitating, and occasionally lethal disease in humans. It is caused by infection by a gram-negative bacterium, Francisella tularensis. In order to better understand its potency as an etiological agent as well as its potential as a biological weapon, we have completed draft assemblies and report the first complete genomic characterization of five strains belonging to the following different Francisella subspecies (subsp.): the F. tularensis subsp. tularensis FSC033, F. tularensis subsp. holarctica FSC257 and FSC022, and F. tularensis subsp. novicida GA99-3548 and GA99-3549 strains. Here, we report the sequencing of these strains and comparative genomic analysis with recently available public Francisella sequences, including the rare F. tularensis subsp. mediasiatica FSC147 strain isolate from the Central Asian Region. We report evidence for the occurrence of large-scale rearrangement events in strains of the holarctica subspecies, supporting previous proposals that further phylogenetic subdivisions of the Type B clade are likely. We also find a significant enrichment of disrupted or absent ORFs proximal to predicted breakpoints in the FSC022 strain, including a genetic component of the Type I restriction-modification defense system. Many of the pseudogenes identified are also disrupted in the closely related rarely human pathogenic F. tularensis subsp. mediasiatica FSC147 strain, including modulator of drug activity B (mdaB) (FTT0961), which encodes a known NADPH quinone reductase involved in oxidative stress resistance. We have also identified genes exhibiting sequence similarity to effectors of the Type III (T3SS) and components of the Type IV secretion systems (T4SS). One of the genes, msrA2 (FTT1797c), is disrupted in F. tularensis subsp. mediasiatica and has recently been shown to mediate bacterial pathogen survival in host organisms. Our findings suggest that in addition to the duplication of the Francisella Pathogenicity Island, and acquisition of individual loci, adaptation by gene loss in the more recently emerged tularensis, holarctica, and mediasiatica subspecies occurred and was distinct from evolutionary events that differentiated these subspecies, and the novicida subspecies, from a common ancestor. Our findings are applicable to future studies focused on variations in Francisella subspecies pathogenesis, and of broader interest to studies of genomic pathoadaptation in bacteria

Comparative Genomic Analysis of Human Fungal Pathogens Causing Paracoccidioidomycosis

Paracoccidioides is a fungal pathogen and the cause of paracoccidioidomycosis, a health-threatening human systemic mycosis endemic to Latin America. Infection by Paracoccidioides, a dimorphic fungus in the order Onygenales, is coupled with a thermally regulated transition from a soil-dwelling filamentous form to a yeast-like pathogenic form. To better understand the genetic basis of growth and pathogenicity in Paracoccidioides, we sequenced the genomes of two strains of Paracoccidioides brasiliensis (Pb03 and Pb18) and one strain of Paracoccidioides lutzii (Pb01). These genomes range in size from 29.1 Mb to 32.9 Mb and encode 7,610 to 8,130 genes. To enable genetic studies, we mapped 94% of the P. brasiliensis Pb18 assembly onto five chromosomes. We characterized gene family content across Onygenales and related fungi, and within Paracoccidioides we found expansions of the fungal-specific kinase family FunK1. Additionally, the Onygenales have lost many genes involved in carbohydrate metabolism and fewer genes involved in protein metabolism, resulting in a higher ratio of proteases to carbohydrate active enzymes in the Onygenales than their relatives. To determine if gene content correlated with growth on different substrates, we screened the non-pathogenic onygenale Uncinocarpus reesii, which has orthologs for 91% of Paracoccidioides metabolic genes, for growth on 190 carbon sources. U. reesii showed growth on a limited range of carbohydrates, primarily basic plant sugars and cell wall components; this suggests that Onygenales, including dimorphic fungi, can degrade cellulosic plant material in the soil. In addition, U. reesii grew on gelatin and a wide range of dipeptides and amino acids, indicating a preference for proteinaceous growth substrates over carbohydrates, which may enable these fungi to also degrade animal biomass. These capabilities for degrading plant and animal substrates suggest a duality in lifestyle that could enable pathogenic species of Onygenales to transfer from soil to animal hosts.National Institute of Allergy and Infectious Diseases (U.S.)National Institutes of Health. Department of Health and Human Services (contract HHSN266200400001C)National Institutes of Health. Department of Health and Human Services(contract HHSN2722009000018C)Brazil. National Council for Scientific and Technological Developmen

A systematic review of the psychometric properties of self-report research utilization measures used in healthcare

<p>Abstract</p> <p>Background</p> <p>In healthcare, a gap exists between what is known from research and what is practiced. Understanding this gap depends upon our ability to robustly measure research utilization.</p> <p>Objectives</p> <p>The objectives of this systematic review were: to identify self-report measures of research utilization used in healthcare, and to assess the psychometric properties (acceptability, reliability, and validity) of these measures.</p> <p>Methods</p> <p>We conducted a systematic review of literature reporting use or development of self-report research utilization measures. Our search included: multiple databases, ancestry searches, and a hand search. Acceptability was assessed by examining time to complete the measure and missing data rates. Our approach to reliability and validity assessment followed that outlined in the <it>Standards for Educational and Psychological Testing</it>.</p> <p>Results</p> <p>Of 42,770 titles screened, 97 original studies (108 articles) were included in this review. The 97 studies reported on the use or development of 60 unique self-report research utilization measures. Seven of the measures were assessed in more than one study. Study samples consisted of healthcare providers (92 studies) and healthcare decision makers (5 studies). No studies reported data on acceptability of the measures. Reliability was reported in 32 (33%) of the studies, representing 13 of the 60 measures. Internal consistency (Cronbach's Alpha) reliability was reported in 31 studies; values exceeded 0.70 in 29 studies. Test-retest reliability was reported in 3 studies with Pearson's <it>r </it>coefficients > 0.80. No validity information was reported for 12 of the 60 measures. The remaining 48 measures were classified into a three-level validity hierarchy according to the number of validity sources reported in 50% or more of the studies using the measure. Level one measures (n = 6) reported evidence from any three (out of four possible) <it>Standards </it>validity sources (which, in the case of single item measures, was all applicable validity sources). Level two measures (n = 16) had evidence from any two validity sources, and level three measures (n = 26) from only one validity source.</p> <p>Conclusions</p> <p>This review reveals significant underdevelopment in the measurement of research utilization. Substantial methodological advances with respect to construct clarity, use of research utilization and related theory, use of measurement theory, and psychometric assessment are required. Also needed are improved reporting practices and the adoption of a more contemporary view of validity (<it>i.e.</it>, the <it>Standards</it>) in future research utilization measurement studies.</p

The evolutionary history of amino acid variations mediating increased resistance of S. aureus identifies reversion mutations in metabolic regulators.

The evolution of resistance in Staphylococcus aureus occurs rapidly, and in response to all known antimicrobial treatments. Numerous studies of model species describe compensatory roles of mutations in mediating competitive fitness, and there is growing evidence that these mutation types also drive adaptation of S. aureus strains. However, few studies have tracked amino acid changes during the complete evolutionary trajectory of antibiotic adaptation or been able to predict their functional relevance. Here, we have assessed the efficacy of computational methods to predict biological resistance of a collection of clinically known Resistance Associated Mutations (RAMs). We have found that >90% of known RAMs are incorrectly predicted to be functionally neutral by at least one of the prediction methods used. By tracing the evolutionary histories of all of the false negative RAMs, we have discovered that a significant number are reversion mutations to ancestral alleles also carried in the MSSA476 methicillin-sensitive isolate. These genetic reversions are most prevalent in strains following daptomycin treatment and show a tendency to accumulate in biological pathway reactions that are distinct from those accumulating non-reversion mutations. Our studies therefore show that in addition to non-reversion mutations, reversion mutations arise in isolates exposed to new antibiotic treatments. It is possible that acquisition of reversion mutations in the genome may prevent substantial fitness costs during the progression of resistance. Our findings pose an interesting question to be addressed by further clinical studies regarding whether or not these reversion mutations lead to a renewed vulnerability of a vancomycin or daptomycin resistant strain to antibiotics administered at an earlier stage of infection

Prevalence of non-reversion and reversion mutations in complex biosynthetic metabolic pathways.

<p>Reversion (A) and non-reversion mutations (B) identified in protein sequences from all MRSA isolates exposed to drug treatment were mapped by homology to the <i>Staphylococcus aureus</i> metabolome <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0056466#pone.0056466-AubryDamon1" target="_blank">[31]</a> (mutations in compounds acting in Nucleoside/nucleotide Biosynthesis (A, orange) and Amino Acid Biosynthesis (B, red) are shown). Compounds accumulating mutations interact in complex pathway networks that involve highly conserved proteins (solid lines) as well as “pathway-holes” (light grey lines) where identified sequence or functional homologs have not yet been identified in the <i>Staphylococcus aureus</i> genomes. The propensity of different mutant types to accumulate in distinct biological pathways, and significant correlation of reversions arising in strains exposed to daptomycin, is suggestive of how exposure to different drugs likely presents evolutionary adaptation biases in metabolic regulation in order to minimize fitness cost while promoting the progression of antibiotic resistance.</p

Estimation of the accuracy of predicting RAMs using three methods.

<p>The accuracy statistics of three methods in predicting StaphRAMs (A–C) is compared to the known accuracy of the same methods for predicting Human DAMs <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0056466#pone.0056466-Lopes1" target="_blank">[21]</a> (D) Maximum likelihood estimation of a binormal ROC curve with an asymmetric 95% confidence interval using 3 methods of prediction (B), with CAROL predictions representing the strongest evidence of a damaging prediction (C). The area of the ROC curve is 0.63 with a standard deviation of 0.06. Accuracy, sensitivity, and specificity were estimated to be 62.6%, 94.0% and 38.5% respectively. Red and blue represent the fitted ROC curve (B). Grey lines denote the 95% confidence interval of the fitted ROC curve (B). A boxplot of the spread estimate for all three methods is shown and ROCR curves for each method are colorized according to the prediction cutoffs (C).</p

IL‐33 facilitates oncogene‐induced cholangiocarcinoma in mice by an interleukin‐6‐sensitive mechanism

UnlabelledCholangiocarcinoma (CCA) is a lethal hepatobiliary neoplasm originating from the biliary apparatus. In humans, CCA risk factors include hepatobiliary inflammation and fibrosis. The recently identified interleukin (IL)-1 family member, IL-33, has been shown to be a biliary mitogen which also promotes liver inflammation and fibrosis. Our aim was to generate a mouse model of CCA mimicking the human disease. Ectopic oncogene expression in the biliary tract was accomplished by the Sleeping Beauty transposon transfection system with transduction of constitutively active AKT (myr-AKT) and Yes-associated protein. Intrabiliary instillation of the transposon-transposase complex was coupled with lobar bile duct ligation in C57BL/6 mice, followed by administration of IL-33 for 3 consecutive days. Tumors developed in 72% of the male mice receiving both oncogenes plus IL-33 by 10 weeks but in only 20% of the male mice transduced with the oncogenes alone. Tumors expressed SOX9 and pancytokeratin (features of CCA) but were negative for HepPar1 (a marker of hepatocellular carcinoma). Substantive overlap with human CCA specimens was revealed by RNA profiling. Not only did IL-33 induce IL-6 expression by human cholangiocytes but it likely facilitated tumor development in vivo by an IL-6-sensitive process as tumor development was significantly attenuated in Il-6(-/-) male animals. Furthermore, tumor formation occurred at a similar rate when IL-6 was substituted for IL-33 in this model.ConclusionThe transposase-mediated transduction of constitutively active AKT and Yes-associated protein in the biliary epithelium coupled with lobar obstruction and IL-33 administration results in the development of CCA with morphological and biochemical features of the human disease; this model highlights the role of inflammatory cytokines in CCA oncogenesis