Shuai Le

Mapping the tail fiber as the receptor binding protein responsible for differential host specificity of Pseudomonas aeruginosa bacteriophages PaP1 and JG004.

The first step in bacteriophage infection is recognition and binding to the host receptor, which is mediated by the phage receptor binding protein (RBP). Different RBPs can lead to differential host specificity. In many bacteriophages, such as Escherichia coli and Lactococcal phages, RBPs have been identified as the tail fiber or protruding baseplate proteins. However, the tail fiber-dependent host specificity in Pseudomonas aeruginosa phages has not been well studied. This study aimed to identify and investigate the binding specificity of the RBP of P. aeruginosa phages PaP1 and JG004. These two phages share high DNA sequence homology but exhibit different host specificities. A spontaneous mutant phage was isolated and exhibited broader host range compared with the parental phage JG004. Sequencing of its putative tail fiber and baseplate region indicated a single point mutation in ORF84 (a putative tail fiber gene), which resulted in the replacement of a positively charged lysine (K) by an uncharged asparagine (N). We further demonstrated that the replacement of the tail fiber gene (ORF69) of PaP1 with the corresponding gene from phage JG004 resulted in a recombinant phage that displayed altered host specificity. Our study revealed the tail fiber-dependent host specificity in P. aeruginosa phages and provided an effective tool for its alteration. These contributions may have potential value in phage therapy.

Genomic and Proteomic Analyses of the Terminally Redundant Genome of the Pseudomonas aeruginosa Phage PaP1: Establishment of Genus PaP1-Like Phages

We isolated and characterized a new Pseudomonas aeruginosa myovirus named PaP1. The morphology of this phage was visualized by electron microscopy and its genome sequence and ends were determined. Finally, genomic and proteomic analyses were performed. PaP1 has an icosahedral head with an apex diameter of 68–70 nm and a contractile tail with a length of 138–140 nm. The PaP1 genome is a linear dsDNA molecule containing 91,715 base pairs (bp) with a G+C content of 49.36% and 12 tRNA genes. A strategy to identify the genome ends of PaP1 was designed. The genome has a 1190 bp terminal redundancy. PaP1 has 157 open reading frames (ORFs). Of these, 143 proteins are homologs of known proteins, but only 38 could be functionally identified. Sodium dodecyl sulfate-polyacrylamide gel electrophoresis and high-performance liquid chromatography-mass spectrometry allowed identification of 12 ORFs as structural protein coding genes within the PaP1 genome. Comparative genomic analysis indicated that the Pseudomonas aeruginosa phage PaP1, JG004, PAK_P1 and vB_PaeM_C2-10_Ab1 share great similarity. Besides their similar biological characteristics, the phages contain 123 core genes and have very close phylogenetic relationships, which distinguish them from other known phage genera. We therefore propose that these four phages be classified as PaP1-like phages, a new phage genus of Myoviridae that infects Pseudomonas aeruginosa.

Chromosomal DNA deletion confers phage resistance to Pseudomonas aeruginosa

Bacteria develop a broad range of phage resistance mechanisms, such as prevention of phage adsorption and CRISPR/Cas system, to survive phage predation. In this study, Pseudomonas aeruginosa PA1 strain was infected with lytic phage PaP1, and phage-resistant mutants were selected. A high percentage (~30%) of these mutants displayed red pigmentation phenotype (Red mutant). Through comparative genomic analysis, one Red mutant PA1r was found to have a 219.6 kb genomic fragment deletion, which contains two key genes hmgA and galU related to the observed phenotypes. Deletion of hmgA resulted in the accumulation of a red compound homogentisic acid; while A galU mutant is devoid of O-antigen, which is required for phage adsorption. Intriguingly, while the loss of galU conferred phage resistance, it significantly attenuated PA1r in a mouse infection experiment. Our study revealed a novel phage resistance mechanism via chromosomal DNA deletion in P. aeruginosa.

Functional Constraint Profiling of a Viral Protein Reveals Discordance of Evolutionary Conservation and Functionality

Viruses often encode proteins with multiple functions due to their compact genomes. Existing approaches to identify functional residues largely rely on sequence conservation analysis. Inferring functional residues from sequence conservation can produce false positives, in which the conserved residues are functionally silent, or false negatives, where functional residues are not identified since they are species-specific and therefore non-conserved. Furthermore, the tedious process of constructing and analyzing individual mutations limits the number of residues that can be examined in a single study. Here, we developed a systematic approach to identify the functional residues of a viral protein by coupling experimental fitness profiling with protein stability prediction using the influenza virus polymerase PA subunit as the target protein. We identified a significant number of functional residues that were influenza type-specific and were evolutionarily non-conserved among different influenza types. Our results indicate that type-specific functional residues are prevalent and may not otherwise be identified by sequence conservation analysis alone. More importantly, this technique can be adapted to any viral (and potentially non-viral) protein where structural information is available.

HIV-1 Quasispecies Delineation by Tag Linkage Deep Sequencing

Trade-offs between throughput, read length, and error rates in high-throughput sequencing limit certain applications such as monitoring viral quasispecies. Here, we describe a molecular-based tag linkage method that allows assemblage of short sequence reads into long DNA fragments. It enables haplotype phasing with high accuracy and sensitivity to interrogate individual viral sequences in a quasispecies. This approach is demonstrated to deduce ∼2000 unique 1.3 kb viral sequences from HIV-1 quasispecies in vivo and after passaging ex vivo with a detection limit of ∼0.005% to ∼0.001%. Reproducibility of the method is validated quantitatively and qualitatively by a technical replicate. This approach can improve monitoring of the genetic architecture and evolution dynamics in any quasispecies population.

Anti-obesity Effect of Capsaicin in Mice Fed with High-Fat Diet Is Associated with an Increase in Population of the Gut Bacterium Akkermansia muciniphila

Capsaicin (CAP) reduces body weight mainly through activation of transient receptor potential vanilloid 1 (TRPV1) cation channel. However, recent evidence indicates that the gut microbiota influences many physiological processes in host and might provoke obesity. This study determined whether the anti-obesity effect of CAP is related to the changes in gut microbiota. C57BL/6 mice were fed either with high-fat diet (HFD) or HFD with CAP (HFD-CAP) for 9 weeks. We observed a significantly reduced weight gain and improved glucose tolerance in HFD-CAP-fed mice compared with HFD-fed mice. 16S rRNA gene sequencing results showed a decrease of phylum Proteobacteria in HFD-CAP-fed mice. In addition, HFD-CAP-fed mice showed a higher abundance of Akkermansia muciniphila, a mucin-degrading bacterium with beneficial effects on host metabolism. Further studies found that CAP directly up-regulates the expression of Mucin 2 gene Muc2 and antimicrobial protein gene Reg3g in the intestine. These data suggest that the anti-obesity effect of CAP is associated with a modest modulation of the gut microbiota.

Global Transcriptomic Analysis of Interactions between Pseudomonas aeruginosa and Bacteriophage PaP3.

The interactions between Bacteriophage (phage) and host bacteria are widespread in nature and influences of phage replication on the host cells are complex and extensive. Here, we investigate genome-wide interactions of Pseudomonas aeruginosa (P. aeruginosa) and its temperate phage PaP3 at five time points during phage infection. Compared to the uninfected host, 38% (2160/5633) genes of phage-infected host were identified as differentially expressed genes (DEGs). Functional analysis of the repressed DEGs revealed infection-stage-dependent pathway communications. Based on gene co-expression analysis, most PaP3 middle genes were predicted to have negative impact on host transcriptional regulators. Sub-network enrichment analysis revealed that adjacent genes of PaP3 interacted with the same host genes and might possess similar functions. Finally, our results suggested that during the whole infection stage, the early genes of PaP3 had stronger regulatory role in host gene expression than middle and late genes, while the host genes involved amino acid metabolism were the most “vulnerable” targets of these phage genes. This work provides the basis for understanding survival mechanisms of parasites and host, and seeking phage gene products that could potentially be used in anti-bacterial infection.

SeqKit: A Cross-Platform and Ultrafast Toolkit for FASTA/Q File Manipulation

FASTA and FASTQ are basic and ubiquitous formats for storing nucleotide and protein sequences. Common manipulations of FASTA/Q file include converting, searching, filtering, deduplication, splitting, shuffling, and sampling. Existing tools only implement some of these manipulations, and not particularly efficiently, and some are only available for certain operating systems. Furthermore, the complicated installation process of required packages and running environments can render these programs less user friendly. This paper describes a cross-platform ultrafast comprehensive toolkit for FASTA/Q processing. SeqKit provides executable binary files for all major operating systems, including Windows, Linux, and Mac OSX, and can be directly used without any dependencies or pre-configurations. SeqKit demonstrates competitive performance in execution time and memory usage compared to similar tools. The efficiency and usability of SeqKit enable researchers to rapidly accomplish common FASTA/Q file manipulations. SeqKit is open source and available on Github at https://github.com/shenwei356/seqkit.

Identification and Characterization of the HicAB Toxin-Antitoxin System in the Opportunistic Pathogen Pseudomonas aeruginosa

Toxin-antitoxin (TA) systems are small genetic modules that are widely distributed in the genomes of bacteria and archaea and have been proposed to fulfill numerous functions. Here, we describe the identification and characterization of a type II TA system, comprising the hicAB locus in the human opportunistic pathogen Pseudomonas aeruginosa. The hicAB locus consists of genes hicA and hicB encoding a toxin and its cognate antitoxin, respectively. BLAST analysis revealed that hicAB is prevalent in approximately 36% of P. aeruginosa strains and locates in the same genomic region. RT-PCR demonstrated that hicAB forms a bicistronic operon that is cotranscribed under normal growth conditions. Overproduction of HicA inhibited the growth of Escherichia coli, and this effect could be counteracted by co-expression of HicB. The Escherichia coli kill/rescue assay showed that the effect of HicA is bacteriostatic, rather than bactericidal. Deletion of hicAB had no effect on the biofilm formation and virulence of P. aeruginosa in a mice infection model. Collectively, this study presents the first characterization of the HicAB system in the opportunistic pathogen P. aeruginosa.

Coupling high-throughput genetics with phylogenetic information reveals an epistatic interaction on the influenza A virus M segment

BackgroundEpistasis is one of the central themes in viral evolution due to its importance in drug resistance, immune escape, and interspecies transmission. However, there is a lack of experimental approach to systematically probe for epistatic residues.ResultsBy utilizing the information from natural occurring sequences and high-throughput genetics, this study established a novel strategy to identify epistatic residues. The rationale is that a substitution that is deleterious in one strain may be prevalent in nature due to the presence of a naturally occurring compensatory substitution. Here, high-throughput genetics was applied to influenza A virus M segment to systematically identify deleterious substitutions. Comparison with natural sequence variation showed that a deleterious substitution M1 Q214H was prevalent in circulating strains. A coevolution analysis was then performed and indicated that M1 residues 121, 207, 209, and 214 naturally coevolved as a group. Subsequently, we experimentally validated that M1 A209T was a compensatory substitution for M1 Q214H.ConclusionsThis work provided a proof-of-concept to identify epistatic residues by coupling high-throughput genetics with phylogenetic information. In particular, we were able to identify an epistatic interaction between M1 substitutions A209T and Q214H. This analytic strategy can potentially be adapted to study any protein of interest, provided that the information on natural sequence variants is available.