Luca Freschi

Clinical utilization of genomics data produced by the international Pseudomonas aeruginosa consortium

The International Pseudomonas aeruginosa Consortium is sequencing over 1000 genomes and building an analysis pipeline for the study of Pseudomonas genome evolution, antibiotic resistance and virulence genes. Metadata, including genomic and phenotypic data for each isolate of the collection, are available through the International Pseudomonas Consortium Database (http://ipcd.ibis.ulaval.ca/). Here, we present our strategy and the results that emerged from the analysis of the first 389 genomes. With as yet unmatched resolution, our results confirm that P. aeruginosa strains can be divided into three major groups that are further divided into subgroups, some not previously reported in the literature. We also provide the first snapshot of P. aeruginosa strain diversity with respect to antibiotic resistance. Our approach will allow us to draw potential links between environmental strains and those implicated in human and animal infections, understand how patients become infected and how the infection evolves over time as well as identify prognostic markers for better evidence-based decisions on patient care.

Phosphorylation network rewiring by gene duplication

Elucidating how complex regulatory networks have assembled during evolution requires a detailed understanding of the evolutionary dynamics that follow gene duplication events, including changes in post‐translational modifications. We compared the phosphorylation profiles of paralogous proteins in the budding yeast Saccharomyces cerevisiae to that of a species that diverged from the budding yeast before the duplication of those genes. We found that 100 million years of post‐duplication divergence are sufficient for the majority of phosphorylation sites to be lost or gained in one paralog or the other, with a strong bias toward losses. However, some losses may be partly compensated for by the evolution of other phosphosites, as paralogous proteins tend to preserve similar numbers of phosphosites over time. We also found that up to 50% of kinase–substrate relationships may have been rewired during this period. Our results suggest that after gene duplication, proteins tend to subfunctionalize at the level of post‐translational regulation and that even when phosphosites are preserved, there is a turnover of the kinases that phosphorylate them.

Characterization of Spindle Checkpoint Kinase Mps1 Reveals Domain with Functional and Structural Similarities to Tetratricopeptide Repeat Motifs of Bub1 and BubR1 Checkpoint Kinases

Background: The N terminus is required for localization and functions of Mps1, Bub1, and BubR1 kinases. Results: A novel Bub1/BubR1-related TPR motif is identified in Mps1 and is required for kinase activity. Conclusion: TPR domain of Mps1 regulates kinase activity, Mps1 chromosome alignment, and checkpoint functions. Significance: Identification of a novel domain in Mps1 enhances our understanding of its contribution to maintaining genome integrity. Kinetochore targeting of the mitotic kinases Bub1, BubR1, and Mps1 has been implicated in efficient execution of their functions in the spindle checkpoint, the self-monitoring system of the eukaryotic cell cycle that ensures chromosome segregation occurs with high fidelity. In all three kinases, kinetochore docking is mediated by the N-terminal region of the protein. Deletions within this region result in checkpoint failure and chromosome segregation defects. Here, we use an interdisciplinary approach that includes biophysical, biochemical, cell biological, and bioinformatics methods to study the N-terminal region of human Mps1. We report the identification of a tandem repeat of the tetratricopeptide repeat (TPR) motif in the N-terminal kinetochore binding region of Mps1, with close homology to the tandem TPR motif of Bub1 and BubR1. Phylogenetic analysis indicates that TPR Mps1 was acquired after the split between deutorostomes and protostomes, as it is distinguishable in chordates and echinoderms. Overexpression of TPR Mps1 resulted in decreased efficiency of both chromosome alignment and mitotic arrest, likely through displacement of endogenous Mps1 from the kinetochore and decreased Mps1 catalytic activity. Taken together, our multidisciplinary strategy provides new insights into the evolution, structural organization, and function of Mps1 N-terminal region.

Functional divergence and evolutionary turnover in mammalian phosphoproteomes.

Protein phosphorylation is a key mechanism to regulate protein functions. However, the contribution of this protein modification to species divergence is still largely unknown. Here, we studied the evolution of mammalian phosphoregulation by comparing the human and mouse phosphoproteomes. We found that 84% of the positions that are phosphorylated in one species or the other are conserved at the residue level. Twenty percent of these conserved sites are phosphorylated in both species. This proportion is 2.5 times more than expected by chance alone, suggesting that purifying selection is preserving phosphoregulation. However, we show that the majority of the sites that are conserved at the residue level are differentially phosphorylated between species. These sites likely result from false-negative identifications due to incomplete experimental coverage, false-positive identifications and non-functional sites. In addition, our results suggest that at least 5% of them are likely to be true differentially phosphorylated sites and may thus contribute to the divergence in phosphorylation networks between mouse and humans and this, despite residue conservation between orthologous proteins. We also showed that evolutionary turnover of phosphosites at adjacent positions (in a distance range of up to 40 amino acids) in human or mouse leads to an over estimation of the divergence in phosphoregulation between these two species. These sites tend to be phosphorylated by the same kinases, supporting the hypothesis that they are functionally redundant. Our results support the hypothesis that the evolutionary turnover of phosphorylation sites contributes to the divergence in phosphorylation profiles while preserving phosphoregulation. Overall, our study provides advanced analyses of mammalian phosphoproteomes and a framework for the study of their contribution to phenotypic evolution.

Turnover of protein phosphorylation evolving under stabilizing selection

Most proteins are regulated by posttranslational modifications and changes in these modifications contribute to evolutionary changes as well as to human diseases. Phosphorylation of serines, threonines, and tyrosines are the most common modifications identified to date in eukaryotic proteomes. While the mode of action and the function of most phosphorylation sites remain unknown, functional studies have shown that phosphorylation affects protein stability, localization and ability to interact. Two broad modes of action have been described for protein phosphorylation. The first mode corresponds to the canonical and qualitative view whereby single phosphorylation sites act as molecular switches that either turn on or off specific protein functions through direct or allosteric effects. The second mode is more akin to a rheostat than a switch. In this case, a group of phosphorylation sites in a given protein region contributes collectively to the modification of the protein, irrespective of the precise position of individual sites, through an aggregate property. Here we discuss these two types of regulation and examine how they affect the rate and patterns of protein phosphorylation evolution. We describe how the evolution of clusters of phosphorylation sites can be studied under the framework of complex traits evolution and stabilizing selection.

The Widespread Multidrug-Resistant Serotype O12 Pseudomonas aeruginosa Clone Emerged through Concomitant Horizontal Transfer of Serotype Antigen and Antibiotic Resistance Gene Clusters

ABSTRACT The O-specific antigen (OSA) in Pseudomonas aeruginosa lipopolysaccharide is highly varied by sugar identity, side chains, and bond between O-repeats. These differences classified P. aeruginosa into 20 distinct serotypes. In the past few decades, O12 has emerged as the predominant serotype in clinical settings and outbreaks. These serotype O12 isolates exhibit high levels of resistance to various classes of antibiotics. Here, we explore how the P. aeruginosa OSA biosynthesis gene clusters evolve in the population by investigating the association between the phylogenetic relationships among 83 P. aeruginosa strains and their serotypes. While most serotypes were closely linked to the core genome phylogeny, we observed horizontal exchange of OSA biosynthesis genes among phylogenetically distinct P. aeruginosa strains. Specifically, we identified a “serotype island” ranging from 62 kb to 185 kb containing the P. aeruginosa O12 OSA gene cluster, an antibiotic resistance determinant (gyrAC248T), and other genes that have been transferred between P. aeruginosa strains with distinct core genome architectures. We showed that these genes were likely acquired from an O12 serotype strain that is closely related to P. aeruginosa PA7. Acquisition and recombination of the “serotype island” resulted in displacement of the native OSA gene cluster and expression of the O12 serotype in the recipients. Serotype switching by recombination has apparently occurred multiple times involving bacteria of various genomic backgrounds. In conclusion, serotype switching in combination with acquisition of an antibiotic resistance determinant most likely contributed to the dissemination of the O12 serotype in clinical settings. IMPORTANCE Infection rates in hospital settings by multidrug-resistant (MDR) Pseudomonas aeruginosa clones have increased during the past decades, and serotype O12 is predominant among these epidemic strains. It is not known why the MDR phenotype is associated with serotype O12 and how this clone type has emerged. This study shows that evolution of MDR O12 strains involved a switch from an ancestral O4 serotype to O12. Serotype switching was the result of horizontal transfer and genetic recombination of lipopolysaccharide (LPS) biosynthesis genes originating from an MDR taxonomic outlier P. aeruginosa strain. Moreover, the recombination event also resulted in acquisition of antibiotic resistance genes. These results impact on our understanding of MDR outbreak strain and serotype evolution and can potentially assist in better monitoring and prevention. Infection rates in hospital settings by multidrug-resistant (MDR) Pseudomonas aeruginosa clones have increased during the past decades, and serotype O12 is predominant among these epidemic strains. It is not known why the MDR phenotype is associated with serotype O12 and how this clone type has emerged. This study shows that evolution of MDR O12 strains involved a switch from an ancestral O4 serotype to O12. Serotype switching was the result of horizontal transfer and genetic recombination of lipopolysaccharide (LPS) biosynthesis genes originating from an MDR taxonomic outlier P. aeruginosa strain. Moreover, the recombination event also resulted in acquisition of antibiotic resistance genes. These results impact on our understanding of MDR outbreak strain and serotype evolution and can potentially assist in better monitoring and prevention.

A Syst-OMICS Approach to Ensuring Food Safety and Reducing the Economic Burden of Salmonellosis

The Salmonella Syst-OMICS consortium is sequencing 4,500 Salmonella genomes and building an analysis pipeline for the study of Salmonella genome evolution, antibiotic resistance and virulence genes. Metadata, including phenotypic as well as genomic data, for isolates of the collection are provided through the Salmonella Foodborne Syst-OMICS database (SalFoS), at https://salfos.ibis.ulaval.ca/. Here, we present our strategy and the analysis of the first 3,377 genomes. Our data will be used to draw potential links between strains found in fresh produce, humans, animals and the environment. The ultimate goals are to understand how Salmonella evolves over time, improve the accuracy of diagnostic methods, develop control methods in the field, and identify prognostic markers for evidence-based decisions in epidemiology and surveillance.

Integrative avenues for exploring the dynamics and evolution of protein interaction networks.

Over the past decade, the study of protein interaction networks (PINs) has shed light on the organizing principles of living cells. However, PINs have been mostly mapped in one single condition. We outline three of the most promising avenues of investigation in this field, namely the study of first, how PINs are rewired by mutations and environmental perturbations; secondly, how inter-species interactions affect PIN achitectures; thirdly, what mechanisms and forces drive PIN evolution. These investigations will unravel the dynamics and condition dependence of PINs and will thus lead to a better functional annotation of network architecture. One major challenge to reach these goals is the integration of PINs with other cellular regulatory networks in the context of complex cellular phenotypes.

Genomics of antibiotic-resistance prediction in Pseudomonas aeruginosa

Antibiotic resistance is a worldwide health issue spreading quickly among human and animal pathogens, as well as environmental bacteria. Misuse of antibiotics has an impact on the selection of resistant bacteria, thus contributing to an increase in the occurrence of resistant genotypes that emerge via spontaneous mutation or are acquired by horizontal gene transfer. There is a specific and urgent need not only to detect antimicrobial resistance but also to predict antibiotic resistance in silico. We now have the capability to sequence hundreds of bacterial genomes per week, including assembly and annotation. Novel and forthcoming bioinformatics tools can predict the resistome and the mobilome with a level of sophistication not previously possible. Coupled with bacterial strain collections and databases containing strain metadata, prediction of antibiotic resistance and the potential for virulence are moving rapidly toward a novel approach in molecular epidemiology. Here, we present a model system in antibiotic‐resistance prediction, along with its promises and limitations. As it is commonly multidrug resistant, Pseudomonas aeruginosa causes infections that are often difficult to eradicate. We review novel approaches for genotype prediction of antibiotic resistance. We discuss the generation of microbial sequence data for real‐time patient management and the prediction of antimicrobial resistance.

A Multi-Species TaqMan PCR Assay for the Identification of Asian Gypsy Moths (Lymantria spp.) and Other Invasive Lymantriines of Biosecurity Concern to North America.

Preventing the introduction and establishment of forest invasive alien species (FIAS) such as the Asian gypsy moth (AGM) is a high-priority goal for countries with extensive forest resources such as Canada. The name AGM designates a group of closely related Lymantria species (Lepidoptera: Erebidae: Lymantriinae) comprising two L. dispar subspecies (L. dispar asiatica, L. dispar japonica) and three closely related Lymantria species (L. umbrosa, L. albescens, L. postalba), all considered potential FIAS in North America. Ships entering Canadian ports are inspected for the presence of suspicious gypsy moth eggs, but those of AGM are impossible to distinguish from eggs of innocuous Lymantria species. To assist regulatory agencies in their identification of these insects, we designed a suite of TaqMan® assays that provide significant improvements over existing molecular assays targeting AGM. The assays presented here can identify all three L. dispar subspecies (including the European gypsy moth, L. dispar dispar), the three other Lymantria species comprising the AGM complex, plus five additional Lymantria species that pose a threat to forests in North America. The suite of assays is built as a “molecular key” (analogous to a taxonomic key) and involves several parallel singleplex and multiplex qPCR reactions. Each reaction uses a combination of primers and probes designed to separate taxa through discriminatory annealing. The success of these assays is based on the presence of single nucleotide polymorphisms (SNPs) in the 5’ region of mitochondrial cytochrome c oxidase I (COI) or in its longer, 3’ region, as well as on the presence of an indel in the “FS1” nuclear marker, generating North American and Asian alleles, used here to assess Asian introgression into L. dispar dispar. These assays have the advantage of providing rapid and accurate identification of ten Lymantria species and subspecies considered potential FIAS.