The phage shock protein (PSP) envelope stress response: discovery of novel partners and evolutionary history.

Bacterial phage shock protein (PSP) systems stabilize the bacterial cell membrane and protect against envelope stress. These systems have been associated with virulence, but despite their critical roles, PSP components are not well characterized outside proteobacteria. Using comparative genomics and protein sequence-structure-function analyses, we systematically identified and analyzed PSP homologs, phyletic patterns, domain architectures, and gene neighborhoods. This approach underscored the evolutionary significance of the system, revealing that its core protein PspA (Snf7 in ESCRT outside bacteria) was present in the last universal common ancestor and that this ancestral functionality has since diversified into multiple novel, distinct PSP systems across life. Several novel partners of the PSP system were identified: (i) the Toastrack domain, likely facilitating assembly of sub-membrane stress-sensing and signaling complexes, (ii) the newly defined HTH-associated α-helical signaling domain-PadR-like transcriptional regulator pair system, and (iii) multiple independent associations with ATPase, CesT/Tir-like chaperone, and Band-7 domains in proteins thought to mediate sub-membrane dynamics. Our work also uncovered links between the PSP components and other domains, such as novel variants of SHOCT-like domains, suggesting roles in assembling membrane-associated complexes of proteins with disparate biochemical functions. Results are available at our interactive web app, https://jravilab.org/psp.IMPORTANCEPhage shock proteins (PSP) are virulence-associated, cell membrane stress-protective systems. They have mostly been characterized in Proteobacteria and Firmicutes. We now show that a minimal PSP system was present in the last universal common ancestor that evolved and diversified into newly identified functional contexts. Recognizing the conservation and evolution of PSP systems across bacterial phyla contributes to our understanding of stress response mechanisms in prokaryotes. Moreover, the newly discovered PSP modularity will likely prompt new studies of lineage-specific cell envelope structures, lifestyles, and adaptation mechanisms. Finally, our results validate the use of domain architecture and genetic context for discovery in comparative genomics.

Global-scale GWAS associates a subset of SNPs with animal-adapted variants in M. tuberculosis complex.

BACKGROUND: While Mycobacterium tuberculosis complex (MTBC) variants are clonal, variant tuberculosis is a human-adapted pathogen, and variant bovis infects many hosts. Despite nucleotide identity between MTBC variants exceeding 99.95%, it remains unclear what drives these differences. Markers of adaptation into variants were sought by bacterial genome-wide association study of single nucleotide polymorphisms extracted from 6,362 MTBC members from varied hosts and countries. RESULTS: The search identified 120 genetic loci associated with MTBC variant classification and certain hosts. In many cases, these changes are uniformly fixed in certain variants while absent in others in this dataset, providing good discriminatory power in distinguishing variants by polymorphisms. Multiple changes were seen in genes for cholesterol and fatty acid metabolism, pathways previously proposed to be important for host adaptation, including Mce4F (part of the fundamental cholesterol intake Mce4 pathway), 4 FadD and FadE genes (playing roles in cholesterol and fatty acid utilization), and other targets like Rv3548c and PTPB, genes shown essential for growth on cholesterol by transposon studies. CONCLUSIONS: These findings provide a robust set of genetic loci associated with the split of variant bovis and variant tuberculosis, and suggest that adaptation to new hosts could involve adjustments in uptake and catabolism of cholesterol and fatty acids, like the proposed specialization to different populations in MTB lineages by alterations to host lipid composition. Future studies are required to elucidate how the associations between cholesterol profiles and pathogen utilization differences between hosts and MTBC variants, as well as the investigation of uncharacterized genes discovered in this study. This information will likely provide an understanding on the diversification of MBO away from humans and specialization towards a broad host range.

Attenuated but immunostimulatory Mycobacterium tuberculosis variant bovis strain Ravenel shows variation in T cell epitopes.

Tuberculosis, caused by Mycobacterium tuberculosis complex (MTBC) organisms, affects a range of humans and animals globally. Mycobacterial pathogenesis involves manipulation of the host immune system, partially through antigen presentation. Epitope sequences across the MTBC are evolutionarily hyperconserved, suggesting their recognition is advantageous for the bacterium. Mycobacterium tuberculosis var. bovis (MBO) strain Ravenel is an isolate known to provoke a robust immune response in cattle, but typically fails to produce lesions and persist. Unlike attenuated MBO BCG strains that lack the critical RD1 genomic region, Ravenel is classic-type MBO structurally, suggesting genetic variation is responsible for defective pathogenesis. This work explores variation in epitope sequences in MBO Ravenel by whole genome sequencing, and contrasts such variation against a fully virulent clinical isolate, MBO strain 10-7428. Validated MTBC epitopes (n = 4818) from the Immune Epitope Database were compared to their sequences in MBO Ravenel and MBO 10-7428. Ravenel yielded 3 modified T cell epitopes, in genes rpfB, argC, and rpoA. These modifications were predicted to have little effect on protein stability. In contrast, T cells epitopes in 10-7428 were all WT. Considering T cell epitope hyperconservation across MTBC variants, these altered MBO Ravenel epitopes support their potential contribution to overall strain attenuation. The affected genes may provide clues on basic pathogenesis, and if so, be feasible targets for reverse vaccinology.

Cold Cas: reevaluating the occurrence of CRISPR/Cas systems in Mycobacteriaceae.

Bacterial CRISPR/Cas systems target foreign genetic elements such as phages and regulate gene expression by some pathogens, even in the host. The system is a marker for evolutionary history and has been used for inferences in Mycobacterium tuberculosis for 30 years. However, knowledge about mycobacterial CRISPR/Cas systems remains limited. It is believed that Type III-A Cas systems are exclusive to Mycobacterium canettii and the M. tuberculosis complex (MTBC) of organisms and that very few of the >200 diverse species of non-tuberculous mycobacteria (NTM) possess any CRISPR/Cas system. This study sought unreported CRISPR/Cas loci across NTM to better understand mycobacterial evolution, particularly in species phylogenetically near the MTBC. An analysis of available mycobacterial genomes revealed that Cas systems are widespread across Mycobacteriaceae and that some species contain multiple types. The phylogeny of Cas loci shows scattered presence in many NTM, with variation even within species, suggesting gains/losses of these loci occur frequently. Cas Type III-A systems were identified in pathogenic Mycobacterium heckeshornense and the geological environmental isolate Mycobacterium SM1. In summary, mycobacterial CRISPR/Cas systems are numerous, Type III-A systems are unreliable as markers for MTBC evolution, and mycobacterial horizontal gene transfer appears to be a frequent source of genetic variation.

Mycobacterium bovis Strain Ravenel Is Attenuated in Cattle.

Mycobacterium tuberculosis variant bovis (MBO) has one of the widest known mammalian host ranges, including humans. Despite the characterization of this pathogen in the 1800s and whole genome sequencing of a UK strain (AF2122) nearly two decades ago, the basis of its host specificity and pathogenicity remains poorly understood. Recent experimental calf infection studies show that MBO strain Ravenel (MBO Ravenel) is attenuated in the cattle host compared to other pathogenic strains of MBO. In the present study, experimental infections were performed to define attenuation. Whole genome sequencing was completed to identify regions of differences (RD) and single nucleotide polymorphisms (SNPs) to explain the observed attenuation. Comparative genomic analysis of MBO Ravenel against three pathogenic strains of MBO (strains AF2122-97, 10-7428, and 95-1315) was performed. Experimental infection studies on five calves each, with either MBO Ravenel or 95-1315, revealed no visible lesions in all five animals in the Ravenel group despite robust IFN-γ responses. Out of 486 polymorphisms in the present analysis, 173 were unique to MBO Ravenel among the strains compared. A high-confidence subset of nine unique SNPs were missense mutations in genes with annotated functions impacting two major MBO survival and virulence pathways: (1) Cell wall synthesis & transport [espH (A103T), mmpL8 (V888I), aftB (H484Y), eccC5 (T507M), rpfB (E263G)], and (2) Lipid metabolism & respiration [mycP1(T125I), pks5 (G455S), fadD29 (N231S), fadE29 (V360G)]. These substitutions likely contribute to the observed attenuation. Results from experimental calf infections and the functional attributions of polymorphic loci on the genome of MBO Ravenel provide new insights into the strain's genotype-disease phenotype associations.

Genome Sequences of Two Mycobacterium tuberculosis Isolates from Asian Elephants in Nepal.

This report describes the genome sequences of two Mycobacterium tuberculosis isolates, S1 and S3, recovered from Asian elephants in Nepal. These genome sequences will enhance our understanding of the genomic epidemiology of Mycobacterium tuberculosis in Asian elephants.

Genome Sequences of Mycobacterium Strains Recovered from Captive Elephants with Tuberculosis.

Members of the Mycobacterium tuberculosis complex cause tuberculosis, infamous for enormous impacts on human health. As zoonoses, they also threaten endangered species like African/Asian elephants. We report the whole-genome sequences of Mycobacterium tuberculosis bv. tuberculosis and Mycobacterium tuberculosis bv. bovis from two zoo elephants in the United States.

Genome Sequences of Mycobacterium tuberculosis Biovar bovis Strains Ravenel and 10-7428.

We report the draft genomes of two Mycobacterium tuberculosis biovar bovis strains. Strain Ravenel was isolated in the 1900s and has been shown to be attenuated in cattle. Strain 10-7428 is considered highly pathogenic in cattle and was isolated from a bovine tuberculosis outbreak.

A Mycobacterium avium subsp. paratuberculosis Predicted Serine Protease Is Associated with Acid Stress and Intraphagosomal Survival.

The ability to maintain intra-cellular pH is crucial for bacteria and other microbes to survive in diverse environments, particularly those that undergo fluctuations in pH. Mechanisms of acid resistance remain poorly understood in mycobacteria. Although, studies investigating acid stress in M. tuberculosis are gaining traction, few center on Mycobacterium avium subsp. paratuberculosis (MAP), the etiological agent of chronic enteritis in ruminants. We identified a MAP acid stress response network involved in macrophage infection. The central node of this network was MAP0403, a predicted serine protease that shared an 86% amino acid identity with MarP in M. tuberculosis. Previous studies confirmed MarP as a serine protease integral to maintaining intra-bacterial pH and survival in acid in vitro and in vivo. We show that MAP0403 is upregulated in infected macrophages and MAC-T cells that coincided with phagosome acidification. Treatment of mammalian cells with bafilomcyin A1, a potent inhibitor of phagosomal vATPases, diminished MAP0403 transcription. MAP0403 expression was also noted in acidic medium. A surrogate host, M. smegmatis mc(2) 155, was designed to express MAP0403 and when exposed to either macrophages or in vitro acid stress had increased bacterial cell viability, which corresponds to maintenance of intra-bacterial pH in acidic (pH = 5) conditions, compared to the parent strain. These data suggest that MAP0403 may be the equivalent of MarP in MAP. Future studies confirming MAP0403 as a serine protease and exploring its structure and possible substrates are warranted.

Conservation analysis of the CydX protein yields insights into small protein identification and evolution.

BACKGROUND: The reliable identification of proteins containing 50 or fewer amino acids is difficult due to the limited information content in short sequences. The 37 amino acid CydX protein in Escherichia coli is a member of the cytochrome bd oxidase complex, an enzyme found throughout Eubacteria. To investigate the extent of CydX conservation and prevalence and evaluate different methods of small protein homologue identification, we surveyed 1095 Eubacteria species for the presence of the small protein. RESULTS: Over 300 homologues were identified, including 80 unannotated genes. The ability of both closely-related and divergent homologues to complement the E. coli ΔcydX mutant supports our identification techniques, and suggests that CydX homologues retain similar function among divergent species. However, sequence analysis of these proteins shows a great degree of variability, with only a few highly-conserved residues. An analysis of the co-variation between CydX homologues and their corresponding cydA and cydB genes shows a close synteny of the small protein with the CydA long Q-loop. Phylogenetic analysis suggests that the cydABX operon has undergone horizontal gene transfer, although the cydX gene likely evolved in a progenitor of the Alpha, Beta, and Gammaproteobacteria. Further investigation of cydAB operons identified two additional conserved hypothetical small proteins: CydY encoded in CydAQlong operons that lack cydX, and CydZ encoded in more than 150 CydAQshort operons. CONCLUSIONS: This study provides a systematic analysis of bioinformatics techniques required for the unique challenges present in small protein identification and phylogenetic analyses. These results elucidate the prevalence of CydX throughout the Proteobacteria, provide insight into the selection pressure and sequence requirements for CydX function, and suggest a potential functional interaction between the small protein and the CydA Q-loop, an enigmatic domain of the cytochrome bd oxidase complex. Finally, these results identify other conserved small proteins encoded in cytochrome bd oxidase operons, suggesting that small protein subunits may be a more common component of these enzymes than previously thought.

The Escherichia coli CydX protein is a member of the CydAB cytochrome bd oxidase complex and is required for cytochrome bd oxidase activity.

Cytochrome bd oxidase operons from more than 50 species of bacteria contain a short gene encoding a small protein that ranges from ∼30 to 50 amino acids and is predicted to localize to the cell membrane. Although cytochrome bd oxidases have been studied for more than 70 years, little is known about the role of this small protein, denoted CydX, in oxidase activity. Here we report that Escherichia coli mutants lacking CydX exhibit phenotypes associated with reduced oxidase activity. In addition, cell membrane extracts from ΔcydX mutant strains have reduced oxidase activity in vitro. Consistent with data showing that CydX is required for cytochrome bd oxidase activity, copurification experiments indicate that CydX interacts with the CydAB cytochrome bd oxidase complex. Together, these data support the hypothesis that CydX is a subunit of the CydAB cytochrome bd oxidase complex that is required for complex activity. The results of mutation analysis of CydX suggest that few individual amino acids in the small protein are essential for function, at least in the context of protein overexpression. In addition, the results of analysis of the paralogous small transmembrane protein AppX show that the two proteins could have some overlapping functionality in the cell and that both have the potential to interact with the CydAB complex.