Infants with cystic fibrosis have altered fecal functional capacities with potential clinical and metabolic consequences

Background

Infants with cystic fibrosis (CF) suffer from gastrointestinal (GI) complications, including pancreatic insufficiency and intestinal inflammation, which have been associated with impaired nutrition and growth. Recent evidence identified altered fecal microbiota taxonomic compositions in infants with CF relative to healthy infants that were characterized by differences in the abundances of taxa associated with GI health and nutrition. Furthermore, these taxonomic differences were more pronounced in low length infants with CF, suggesting a potential link to linear growth failure. We hypothesized that these differences would entail shifts in the microbiome’s functional capacities that could contribute to inflammation and nutritional failure in infants with CF.

Results

To test this hypothesis, we compared fecal microbial metagenomic content between healthy infants and infants with CF, supplemented with an analysis of fecal metabolomes in infants with CF. We identified notable differences in CF fecal microbial functional capacities, including metabolic and environmental response functions, compared to healthy infants that intensified during the first year of life. A machine learning-based longitudinal metagenomic age analysis of healthy and CF fecal metagenomic functional profiles further demonstrated that these differences are characterized by a CF-associated delay in the development of these functional capacities. Moreover, we found metagenomic differences in functions related to metabolism among infants with CF that were associated with diet and antibiotic exposure, and identified several taxa as potential drivers of these functional differences. An integrated metagenomic and metabolomic analysis further revealed that abundances of several fecal GI metabolites important for nutrient absorption, including three bile acids, correlated with specific microbes in infants with CF.

Conclusions

Our results highlight several metagenomic and metabolomic factors, including bile acids and other microbial metabolites, that may impact nutrition, growth, and GI health in infants with CF. These factors could serve as promising avenues for novel microbiome-based therapeutics to improve health outcomes in these infants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02305-z.

Maintenance tobramycin primarily affects untargeted bacteria in the CF sputum microbiome

RATIONALE

The most common antibiotic used to treat people with cystic fibrosis (PWCF) is inhaled tobramycin, administered as maintenance therapy for chronic Pseudomonas aeruginosa lung infections. While the effects of inhaled tobramycin on P. aeruginosa abundance and lung function diminish with continued therapy, this maintenance treatment is known to improve long-term outcomes, underscoring how little is known about why antibiotics work in CF infections, what their effects are on complex CF sputum microbiomes and how to improve these treatments.

OBJECTIVES

To rigorously define the effect of maintenance tobramycin on CF sputum microbiome characteristics.

METHODS AND MEASUREMENTS

We collected sputum from 30 PWCF at standardised times before, during and after a single month-long course of maintenance inhaled tobramycin. We used traditional culture, quantitative PCR and metagenomic sequencing to define the dynamic effects of this treatment on sputum microbiomes, including abundance changes in both clinically targeted and untargeted bacteria, as well as functional gene categories.

MAIN RESULTS

CF sputum microbiota changed most markedly by 1 week of antibiotic therapy and plateaued thereafter, and this shift was largely driven by changes in non-dominant taxa. The genetically conferred functional capacities (ie, metagenomes) of subjects' sputum communities changed little with antibiotic perturbation, despite taxonomic shifts, suggesting functional redundancy within the CF sputum microbiome.

CONCLUSIONS

Maintenance treatment with inhaled tobramycin, an antibiotic with demonstrated long-term mortality benefit, primarily impacted clinically untargeted bacteria in CF sputum, highlighting the importance of monitoring the non-canonical effects of antibiotics and other treatments to accurately define and improve their clinical impact.

Human and Extracellular DNA Depletion for Metagenomic Analysis of Complex Clinical Infection Samples Yields Optimized Viable Microbiome Profiles

Metagenomic sequencing is a promising approach for identifying and characterizing organisms and their functional characteristics in complex, polymicrobial infections, such as airway infections in people with cystic fibrosis. These analyses are often hampered, however, by overwhelming quantities of human DNA, yielding only a small proportion of microbial reads for analysis. In addition, many abundant microbes in respiratory samples can produce large quantities of extracellular bacterial DNA originating either from biofilms or dead cells. We describe a method for simultaneously depleting DNA from intact human cells and extracellular DNA (human and bacterial) in sputum, using selective lysis of eukaryotic cells and endonuclease digestion. We show that this method increases microbial sequencing depth and, consequently, both the number of taxa detected and coverage of individual genes such as those involved in antibiotic resistance. This finding underscores the substantial impact of DNA from sources other than live bacteria in micro-biological analyses of complex, chronic infection specimens.

Nelson et al. describe a method for reducing both human cellular DNA and extracellular DNA (human and bacterial) in a complex respiratory sample using hypotonic lysis and endonuclease digestion. This method increases effective microbial sequencing depth and minimizes bias introduced into subsequent phylogenetic analysis by bacterial extracellular DNA.

CFTR dysregulation drives active selection of the gut microbiome

Patients with cystic fibrosis (CF) have altered fecal microbiomes compared to those of healthy controls. The magnitude of this dysbiosis correlates with measures of CF gastrointestinal (GI) disease, including GI inflammation and nutrient malabsorption. However, whether this dysbiosis is caused by mutations in the CFTR gene, the underlying defect in CF, or whether CF-associated dysbiosis augments GI disease was not clear. To test the relationships between CFTR dysfunction, microbes, and intestinal health, we established a germ-free (GF) CF mouse model and demonstrated that CFTR gene mutations are sufficient to alter the GI microbiome. Furthermore, flow cytometric analysis demonstrated that colonized CF mice have increased mesenteric lymph node and spleen TH17+ cells compared with non-CF mice, suggesting that CFTR defects alter adaptive immune responses. Our findings demonstrate that CFTR mutations modulate both the host adaptive immune response and the intestinal microbiome.

It has been difficult to establish causal relationships between host genetics and the selection of the vast multitude of micro-organisms that live in and on us (i.e. the microbiota). Cystic fibrosis has been shown to be associated with changes in the fecal microbiome (the genetic constitution of the microbiota) although it was not evident whether mutation of CFTR, the gene mutated in CF, could drive this selection or whether the frequent use of antibiotics in this population was at fault. Here, by using a germfree (i.e. sterile, lacking all microbiota) mouse model of CF we clearly demonstrate that mutated CFTR alone can alter the microbiome. We also show an increase in an adaptive immune cell type (TH17 cells) in the mesenteric lymph nodes and spleens of CF mice compared to control mice. Our study provides new insights into the dominant role that CFTR plays in microbiome determination and suggests that therapies restoring CFTR function could also correct the microbial dysbiosis observed in CF.

In vivo protein interaction network analysis reveals porin-localized antibiotic inactivation in Acinetobacter baumannii strain AB5075

The nosocomial pathogen Acinetobacter baumannii is a frequent cause of hospital-acquired infections worldwide and is a challenge for treatment due to its evolved resistance to antibiotics, including carbapenems. Here, to gain insight on A. baumannii antibiotic resistance mechanisms, we analyse the protein interaction network of a multidrug-resistant A. baumannii clinical strain (AB5075). Using in vivo chemical cross-linking and mass spectrometry, we identify 2,068 non-redundant cross-linked peptide pairs containing 245 intra- and 398 inter-molecular interactions. Outer membrane proteins OmpA and YiaD, and carbapenemase Oxa-23 are hubs of the identified interaction network. Eighteen novel interactors of Oxa-23 are identified. Interactions of Oxa-23 with outer membrane porins OmpA and CarO are verified with co-immunoprecipitation analysis. Furthermore, transposon mutagenesis of oxa-23 or interactors of Oxa-23 demonstrates changes in meropenem or imipenem sensitivity in strain AB5075. These results provide a view of porin-localized antibiotic inactivation and increase understanding of bacterial antibiotic resistance mechanisms.

The bacterial pathogen Acinetobacter baumannii has evolved resistance to many antibiotics, including carbapenems. Here, Wu et al. show that the carbapenemase Oxa-23 interacts with the outer membrane porin CarO in an A. baumannii isolate, indicative of porin-localised antibiotic inactivation.

GUTSS: An Alignment-Free Sequence Comparison Method for Use in Human Intestinal Microbiome and Fecal Microbiota Transplantation Analysis

BACKGROUND

Comparative analysis of gut microbiomes in clinical studies of human diseases typically rely on identification and quantification of species or genes. In addition to exploring specific functional characteristics of the microbiome and potential significance of species diversity or expansion, microbiome similarity is also calculated to study change in response to therapies directed at altering the microbiome. Established ecological measures of similarity can be constructed from species abundances, however methods for calculating these commonly used ecological measures of similarity directly from whole genome shotgun (WGS) metagenomic sequence are lacking.

RESULTS

We present an alignment-free method for calculating similarity of WGS metagenomic sequences that is analogous to the Bray-Curtis index for species, implemented by the General Utility for Testing Sequence Similarity (GUTSS) software application. This method was applied to intestinal microbiomes of healthy young children to measure developmental changes toward an adult microbiome during the first 3 years of life. We also calculate similarity of donor and recipient microbiomes to measure establishment, or engraftment, of donor microbiota in fecal microbiota transplantation (FMT) studies focused on mild to moderate Crohn's disease. We show how a relative index of similarity to donor can be calculated as a measure of change in a patient's microbiome toward that of the donor in response to FMT.

CONCLUSION

Because clinical efficacy of the transplant procedure cannot be fully evaluated without analysis methods to quantify actual FMT engraftment, we developed a method for detecting change in the gut microbiome that is independent of species identification and database bias, sensitive to changes in relative abundance of the microbial constituents, and can be formulated as an index for correlating engraftment success with clinical measures of disease. More generally, this method may be applied to clinical evaluation of human microbiomes and provide potential diagnostic determination of individuals who may be candidates for specific therapies directed at alteration of the microbiome.

Low Level Engraftment and Improvement following a Single Colonoscopic Administration of Fecal Microbiota to Patients with Ulcerative Colitis

Objective

Fecal microbiota transplantation (FMT) is an investigational treatment for diseases thought to involve alterations in the intestinal microbiota including ulcerative colitis (UC). Case reports have described therapeutic benefit of FMT in patients with UC, possibly due to changes in the microbiota. We measured the degree to which the transplanted microbiota engraft following FMT in patients with UC using a donor similarity index (DSI).

Methods

Seven patients with mild to moderate UC (UC disease activity index scores 3–10) received a single colonoscopic administration of FMT. Metagenomic sequence data from stool were analyzed using an alignment-free comparison tool, to measure the DSI, and a phylogenetic analysis tool, to characterize taxonomic changes. Clinical, endoscopic, histologic, and fecal calprotectin outcome measures were also collected.

Results

One of 5 patients from whom sequencing data were available achieved the primary endpoint of 50% donor similarity at week 4; an additional 2 patients achieved 40% donor similarity. One patient with 40% donor similarity achieved clinical and histologic remission 1 month after FMT. However, these were lost by 2−3 months, and loss correlated with a decrease in DSI. The remaining patients did not demonstrate clinical response or remission. Histology scores improved in all but 1 patient. No patients remained in remission at 3 months after FMT.

Conclusions

Following a single colonoscopic fecal transplant, a DSI of 40-50% is achieved in about two-thirds of recipients. This level of engraftment correlated with a temporary clinical improvement in only 1/5 patients. Larger sample sizes could further validate this method for measuring engraftment, and changes in transplant frequency or method might improve microbiota engraftment and efficacy.

Trial Registration

ClinicalTrials.gov NCT01742754

An underlying mechanism for the increased mutagenesis of lagging-strand genes in Bacillus subtilis

We previously reported that lagging-strand genes accumulate mutations faster than those encoded on the leading strand in Bacillus subtilis. Although we proposed that orientation-specific encounters between replication and transcription underlie this phenomenon, the mechanism leading to the increased mutagenesis of lagging-strand genes remained unknown. Here, we report that the transcription-dependent and orientation-specific differences in mutation rates of genes require the B. subtilis Y-family polymerase, PolY1 (yqjH). We find that without PolY1, association of the replicative helicase, DnaC, and the recombination protein, RecA, with lagging-strand genes increases in a transcription-dependent manner. These data suggest that PolY1 promotes efficient replisome progression through lagging-strand genes, thereby reducing potentially detrimental breaks and single-stranded DNA at these loci. Y-family polymerases can alleviate potential obstacles to replisome progression by facilitating DNA lesion bypass, extension of D-loops, or excision repair. We find that the nucleotide excision repair (NER) proteins UvrA, UvrB, and UvrC, but not RecA, are required for transcription-dependent asymmetry in mutation rates of genes in the two orientations. Furthermore, we find that the transcription-coupling repair factor Mfd functions in the same pathway as PolY1 and is also required for increased mutagenesis of lagging-strand genes. Experimental and SNP analyses of B. subtilis genomes show mutational footprints consistent with these findings. We propose that the interplay between replication and transcription increases lesion susceptibility of, specifically, lagging-strand genes, activating an Mfd-dependent error-prone NER mechanism. We propose that this process, at least partially, underlies the accelerated evolution of lagging-strand genes.

Genomic analysis of the emergence of 20th century epidemic dysentery

Background

Shigella dysenteriae type 1 (Sd1) causes recurrent epidemics of dysentery associated with high mortality in many regions of the world. Sd1 infects humans at very low infectious doses (10 CFU), and treatment is complicated by the rapid emergence of antibiotic resistant Sd1 strains. Sd1 is only detected in the context of human infections, and the circumstances under which epidemics emerge and regress remain unknown.

Results

Phylogenomic analyses of 56 isolates collected worldwide over the past 60 years indicate that the Sd1 clone responsible for the recent pandemics emerged at the turn of the 20^th century, and that the two world wars likely played a pivotal role for its dissemination. Several lineages remain ubiquitous and their phylogeny indicates several recent intercontinental transfers. Our comparative genomics analysis reveals that isolates responsible for separate outbreaks, though closely related to one another, have independently accumulated antibiotic resistance genes, suggesting that there is little or no selection to retain these genes in-between outbreaks. The genomes appear to be subjected to genetic drift that affects a number of functions currently used by diagnostic tools to identify Sd1, which could lead to the potential failure of such tools.

Conclusions

Taken together, the Sd1 population structure and pattern of evolution suggest a recent emergence and a possible human carrier state that could play an important role in the epidemic pattern of infections of this human-specific pathogen. This analysis highlights the important role of whole-genome sequencing in studying pathogens for which epidemiological or laboratory investigations are particularly challenging.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-355) contains supplementary material, which is available to authorized users.

Gastrointestinal pathology in juvenile and adult CFTR-knockout ferrets

Cystic fibrosis (CF) is a multiorgan disease caused by loss of a functional cystic fibrosis transmembrane conductance regulator (CFTR) chloride channel in many epithelia of the body. Here we report the pathology observed in the gastrointestinal organs of juvenile to adult CFTR-knockout ferrets. CF gastrointestinal manifestations included gastric ulceration, intestinal bacterial overgrowth with villous atrophy, and rectal prolapse. Metagenomic phylogenetic analysis of fecal microbiota by deep sequencing revealed considerable genotype-independent microbial diversity between animals, with the majority of taxa overlapping between CF and non-CF pairs. CF hepatic manifestations were variable, but included steatosis, necrosis, biliary hyperplasia, and biliary fibrosis. Gallbladder cystic mucosal hyperplasia was commonly found in 67% of CF animals. The majority of CF animals (85%) had pancreatic abnormalities, including extensive fibrosis, loss of exocrine pancreas, and islet disorganization. Interestingly, 2 of 13 CF animals retained predominantly normal pancreatic histology (84% to 94%) at time of death. Fecal elastase-1 levels from these CF animals were similar to non-CF controls, whereas all other CF animals evaluated were pancreatic insufficient (<2 μg elastase-1 per gram of feces). These findings suggest that genetic factors likely influence the extent of exocrine pancreas disease in CF ferrets and have implications for the etiology of pancreatic sufficiency in CF patients. In summary, these studies demonstrate that the CF ferret model develops gastrointestinal pathology similar to CF patients.

Escherichia coli dysbiosis correlates with gastrointestinal dysfunction in children with cystic fibrosis

Cystic fibrosis gastrointestinal disease includes nutrient malabsorption and intestinal inflammation. We show that the abundances of Escherichia coli in fecal microbiota were significantly higher in young children with cystic fibrosis than in controls and correlated with fecal measures of nutrient malabsorption and inflammation, suggesting that E. coli could contribute to cystic fibrosis gastrointestinal dysfunction.

Evolution of Burkholderia pseudomallei in recurrent melioidosis

Burkholderia pseudomallei, the etiologic agent of human melioidosis, is capable of causing severe acute infection with overwhelming septicemia leading to death. A high rate of recurrent disease occurs in adult patients, most often due to recrudescence of the initial infecting strain. Pathogen persistence and evolution during such relapsing infections are not well understood. Bacterial cells present in the primary inoculum and in late infections may differ greatly, as has been observed in chronic disease, or they may be genetically similar. To test these alternative models, we conducted whole-genome comparisons of clonal primary and relapse B. pseudomallei isolates recovered six months to six years apart from four adult Thai patients. We found differences within each of the four pairs, and some, including a 330 Kb deletion, affected substantial portions of the genome. Many of the changes were associated with increased antibiotic resistance. We also found evidence of positive selection for deleterious mutations in a TetR family transcriptional regulator from a set of 107 additional B. pseudomallei strains. As part of the study, we sequenced to base-pair accuracy the genome of B. pseudomallei strain 1026b, the model used for genetic studies of B. pseudomallei pathogenesis and antibiotic resistance. Our findings provide new insights into pathogen evolution during long-term infections and have important implications for the development of intervention strategies to combat recurrent melioidosis.

PGAT: a multistrain analysis resource for microbial genomes

Motivation: The Prokaryotic-genome Analysis Tool (PGAT) is a web-based database application for comparing gene content and sequence across multiple microbial genomes facilitating the discovery of genetic differences that may explain observed phenotypes. PGAT supports database queries to identify genes that are present or absent in user-selected genomes, comparison of sequence polymorphisms in sets of orthologous genes, multigenome display of regions surrounding a query gene, comparison of the distribution of genes in metabolic pathways and manual community annotation.

Availability and Implementation:The PGAT website may be accessed at http://nwrce.org/pgat.

Contact:mbrittna@uw.edu

A genome-wide in vitro bacterial-infection screen reveals human variation in the host response associated with inflammatory disease

Recent progress in cataloguing common genetic variation has made possible genome-wide studies that are beginning to elucidate the causes and consequences of our genetic differences. Approaches that provide a mechanistic understanding of how genetic variants function to alter disease susceptibility and why they were substrates of natural selection would complement other approaches to human-genome analysis. Here we use a novel cell-based screen of bacterial infection to identify human variation in Salmonella-induced cell death. A loss-of-function allele of CARD8, a reported inhibitor of the proinflammatory protease caspase-1, was associated with increased cell death in vitro (p = 0.013). The validity of this association was demonstrated through overexpression of alternative alleles and RNA interference in cells of varying genotype. Comparison of mammalian CARD8 orthologs and examination of variation among different human populations suggest that the increase in infectious-disease burden associated with larger animal groups (i.e., herds and colonies), and possibly human population expansion, may have naturally selected for loss of CARD8. We also find that the loss-of-function CARD8 allele shows a modest association with an increased risk of systemic inflammatory response syndrome in a small study (p = 0.05). Therefore, a by-product of the selected benefit of loss of CARD8 could be increased inflammatory diseases. These results demonstrate the utility of genome-wide cell-based association screens with microbes in the identification of naturally selected variants that can impact human health.

A Francisella mutant in lipid A carbohydrate modification elicits protective immunity

Francisella tularensis (Ft) is a highly infectious Gram-negative bacterium and the causative agent of the human disease tularemia. Ft is designated a class A select agent by the Centers for Disease Control and Prevention. Human clinical isolates of Ft produce lipid A of similar structure to Ft subspecies novicida (Fn), a pathogen of mice. We identified three enzymes required for Fn lipid A carbohydrate modifications, specifically the presence of mannose (flmF1), galactosamine (flmF2), or both carbohydrates (flmK). Mutants lacking either galactosamine (flmF2) or galactosamine/mannose (flmK) addition to their lipid A were attenuated in mice by both pulmonary and subcutaneous routes of infection. In addition, aerosolization of the mutants (flmF2 and flmK) provided protection against challenge with wild-type (WT) Fn, whereas subcutaneous administration of only the flmK mutant provided protection from challenge with WT Fn. Furthermore, infection of an alveolar macrophage cell line by the flmK mutant induced higher levels of tumor necrosis factor-α (TNF-α) and macrophage inhibitory protein-2 (MIP-2) when compared to infection with WT Fn. Bone marrow–derived macrophages (BMMø) from Toll-like receptor 4 (TLR4) and TLR2/4 knockout mice infected with the flmK mutant also produced significantly higher amounts of interleukin-6 (IL-6) and MIP-2 than BMMø infected with WT Fn. However, production of IL-6 and MIP-2 was undetectable in BMMø from MyD88^−/− mice infected with either strain. MyD88^−/− mice were also susceptible to flmK mutant infection. We hypothesize that the ability of the flmK mutant to activate pro-inflammatory cytokine/chemokine production and innate immune responses mediated by the MyD88 signaling pathway may be responsible for its attenuation, leading to the induction of protective immunity by this mutant.

Bacterial pathogens modify outer membrane components, such as lipid A or endotoxin, the lipid anchor of lipopolysaccharide, to enhance the ability to colonize, spread to different tissues, and/or avoid the host's immune defenses. Lipopolysaccharide also plays an essential role in maintaining membrane integrity and is a key factor in host innate immune recognition of Gram-negative bacterial infections. Francisella tularensis is the causative agent of the human disease tularemia and is classified as a category A select agent. Francisella novicida (Fn) is the murine counterpart of F. tularensis. The structure of Francisella spp. lipid A is unique in that it is modified by various carbohydrates that play a role in virulence and altered endotoxicity. In our study, we identified and defined the role of three genes involved in the carbohydrate modification of the base Fn lipid A structure. We showed that the lack of specific modification(s) of the Fn lipid A molecule lead to bacterial attenuation and activation of a protective immune response against a lethal wild-type infection. Therefore, alteration of Francisella lipid A structure may represent a pathogenesis strategy common to the Francisella species, and specific lipid A mutant strains may be candidates for inclusion in future vaccine studies.

PSAT: a web tool to compare genomic neighborhoods of multiple prokaryotic genomes

Background

The conservation of gene order among prokaryotic genomes can provide valuable insight into gene function, protein interactions, or events by which genomes have evolved. Although some tools are available for visualizing and comparing the order of genes between genomes of study, few support an efficient and organized analysis between large numbers of genomes. The Prokaryotic Sequence homology Analysis Tool (PSAT) is a web tool for comparing gene neighborhoods among multiple prokaryotic genomes.

Results

PSAT utilizes a database that is preloaded with gene annotation, BLAST hit results, and gene-clustering scores designed to help identify regions of conserved gene order. Researchers use the PSAT web interface to find a gene of interest in a reference genome and efficiently retrieve the sequence homologs found in other bacterial genomes. The tool generates a graphic of the genomic neighborhood surrounding the selected gene and the corresponding regions for its homologs in each comparison genome. Homologs in each region are color coded to assist users with analyzing gene order among various genomes. In contrast to common comparative analysis methods that filter sequence homolog data based on alignment score cutoffs, PSAT leverages gene context information for homologs, including those with weak alignment scores, enabling a more sensitive analysis. Features for constraining or ordering results are designed to help researchers browse results from large numbers of comparison genomes in an organized manner. PSAT has been demonstrated to be useful for helping to identify gene orthologs and potential functional gene clusters, and detecting genome modifications that may result in loss of function.

Conclusion

PSAT allows researchers to investigate the order of genes within local genomic neighborhoods of multiple genomes. A PSAT web server for public use is available for performing analyses on a growing set of reference genomes through any web browser with no client side software setup or installation required. Source code is freely available to researchers interested in setting up a local version of PSAT for analysis of genomes not available through the public server. Access to the public web server and instructions for obtaining source code can be found at .

Comparison of Francisella tularensis genomes reveals evolutionary events associated with the emergence of human pathogenic strains

.Sequencing of the non-pathogenic Francisella tularensis sub-species novicida U112, and comparison with two pathogenic sub-species, provides insights into the evolution of pathogenicity in these species.

Background

Francisella tularensis subspecies tularensis and holarctica are pathogenic to humans, whereas the two other subspecies, novicida and mediasiatica, rarely cause disease. To uncover the factors that allow subspecies tularensis and holarctica to be pathogenic to humans, we compared their genome sequences with the genome sequence of Francisella tularensis subspecies novicida U112, which is nonpathogenic to humans.

Results

Comparison of the genomes of human pathogenic Francisella strains with the genome of U112 identifies genes specific to the human pathogenic strains and reveals pseudogenes that previously were unidentified. In addition, this analysis provides a coarse chronology of the evolutionary events that took place during the emergence of the human pathogenic strains. Genomic rearrangements at the level of insertion sequences (IS elements), point mutations, and small indels took place in the human pathogenic strains during and after differentiation from the nonpathogenic strain, resulting in gene inactivation.

Conclusion

The chronology of events suggests a substantial role for genetic drift in the formation of pseudogenes in Francisella genomes. Mutations that occurred early in the evolution, however, might have been fixed in the population either because of evolutionary bottlenecks or because they were pathoadaptive (beneficial in the context of infection). Because the structure of Francisella genomes is similar to that of the genomes of other emerging or highly pathogenic bacteria, this evolutionary scenario may be shared by pathogens from other species.

MglA regulates Francisella tularensis subsp. novicida (Francisella novicida) response to starvation and oxidative stress

MglA is a transcriptional regulator of genes that contribute to the virulence of Francisella tularensis, a highly infectious pathogen and the causative agent of tularemia. This study used a label-free shotgun proteomics method to determine the F. tularensis subsp. novicida (F. novicida) proteins that are regulated by MglA. The differences in relative protein amounts between wild-type F. novicida and the mglA mutant were derived directly from the average peptide precursor ion intensity values measured with the mass spectrometer by using a suite of mathematical algorithms. Among the proteins whose relative amounts changed in an F. novicida mglA mutant were homologs of oxidative and general stress response proteins. The F. novicida mglA mutant exhibited decreased survival during stationary-phase growth and increased susceptibility to killing by superoxide generated by the redox-cycling agent paraquat. The F. novicida mglA mutant also showed increased survival upon exposure to hydrogen peroxide, likely due to increased amounts of the catalase KatG. Our results suggested that MglA coordinates the stress response of F. tularensis and is likely essential for bacterial survival in harsh environments.

Type IV pili-mediated secretion modulates Francisella virulence

Francisella tularensis are the causative agent of the zoonotic disease, tularaemia. Among four F. tularensis subspecies, ssp. novicida (F. novicida) is pathogenic only for immunocompromised individuals, while all four subspecies are pathogenic for mice. This study utilized proteomic and bioinformatic approaches to identify seven F. novicida secreted proteins and the corresponding Type IV pilus (T4P) secretion system. The secreted proteins were predicted to encode two chitinases, a chitin binding protein, a protease (PepO), and a beta-glucosidase (BglX). The transcription of F. novicida pepO and bglX was regulated by the virulence regulator MglA. Intradermal infection of mice with F. novicida mutants defective in T4P secretion system or PepO resulted in enhanced F. novicida spread to systemic sites. Infection with F. novicida pepO mutants also resulted in increased neutrophil infiltration into the mouse airways. PepO is a zinc protease that is homologous to mammalian endothelin-converting enzyme ECE-1. Therefore, secretion of PepO likely results in increased production of endothelin and increased vasoconstriction at the infection site in skin that limits the F. novicida spread. Francisella human pathogenic strains contain a mutation in pepO predicted to abolish its secretion. Loss of PepO function may have contributed to evolution of highly virulent Francisellae.