Identification of epistatic mutations and insights into the evolution of the influenza virus using a mass-based protein phylogenetic approach

Structural Phylogenetics with Protein Mass Spectrometry: A Proof-of-Concept

It is demonstrated, for the first time, that a mass spectrometry approach (known as phylonumerics) can be successfully implemented for structural phylogenetics investigations to chart the evolution of a protein's structure and function. Illustrated for the compact globular protein myoglobin, peptide masses produced from the proteolytic digestion of the protein across animal species generate trees congruent to the sequence tree counterparts. Single point mutations calculated during the same mass tree building step can be followed along interconnected branches of the tree and represent a viable structural metric. A mass tree built for 15 diverse animal species, easily resolve the birds from mammal species, and the ruminant mammals from the remainder of the animals. Mutations within helix-spanning peptide segments alter both the mass and structure of the protein in these segments. Greater evolution is found in the B-helix over the A, E, F, G and H helices. A further mass tree study, of six more closely related primate species, resolves gorilla from the other primates based on a P22S mutation within the B-helix. The remaining five primates are resolved into two groups based on whether they contain a glycine or serine at position 23 in the same helix. The orangutan is resolved from the gibbon and siamang by its G-helix C110S mutation, while homo sapiens are resolved from chimpanzee based on the Q116H mutation. All are associated with structural perturbations in such helices. These structure altering mutations can be tracked along interconnecting branches of a mass tree, to follow the protein's structure and evolution, and ultimately the evolution of the species in which the proteins are expressed. Those that have the greatest impact on a protein's structure, its function, and ultimately the evolution of the species, can be selectively tracked or monitored.

25 Years Responding to Respiratory and Other Viruses with Mass Spectrometry

This review article presents the development and application of mass spectrometry (MS) approaches, developed in the author's laboratory over the past 25 years, to detect; characterise, type and subtype; and distinguish major variants and subvariants of respiratory viruses such as influenza and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). All features make use of matrix-assisted laser desorption ionisation (MALDI) mass maps, recorded for individual viral proteins or whole virus digests. A MALDI-based immunoassay in which antibody-peptide complexes were preserved on conventional MALDI targets without their immobilisation led to an approach that enabled their indirect detection. The site of binding, and thus the molecular antigenicity of viruses, could be determined. The same approach was employed to study antivirals bound to their target viral protein, the nature of the binding residues, and relative binding affinities. The benefits of high-resolution MS were exploited to detect sequence-conserved signature peptides of unique mass within whole virus and single protein digests. These enabled viruses to be typed, subtyped, their lineage determined, and variants and subvariants to be distinguished. Their detection using selected ion monitoring improved analytical sensitivity limits to aid the identification of viruses in clinical specimens. The same high-resolution mass map data, for a wide range of viral strains, were input into a purpose-built algorithm (MassTree) in order to both chart and interrogate viral evolution. Without the need for gene or protein sequences, or any sequence alignment, this phylonumerics approach also determines and displays single-point mutations associated with viral protein evolution in a single-tree building step.

The putative metal-binding proteome of the Coronaviridae family

Metalloproteins are well-known for playing various physicochemical processes in all life forms, including viruses. Some life-threatening viruses (such as some members of the Coronaviridae family of viruses) are emerged and remerged frequently and are rapidly transmitted throughout the globe. This study aims to identify and characterize the metal-binding proteins (MBPs) of the Coronaviridae family of viruses and further provides insight into the MBP's role in sustaining and propagating viruses inside a host cell and in the outer environment. In this study, the available proteome of the Coronaviridae family was exploited. Identified potential MBPs were analyzed for their functional domains, structural aspects, and subcellular localization. We also demonstrate phylogenetic aspects of all predicted MBPs among other Coronaviridae family members to understand the evolutionary trend among their respective hosts. A total of 256 proteins from 51 different species of coronaviruses are predicted as MBPs. These MBPs perform various key roles in the replication and survival of viruses within the host cell. Cysteine, aspartic acid, threonine, and glutamine are key amino acid residues interacting with respective metal ions. Our observations also indicate that the metalloproteins of this family of viruses circulated and evolved in different hosts, which supports the zoonotic nature of coronaviruses. The comprehensive information on MBPs of the Coronaviridae family may be further helpful in designing novel therapeutic metalloprotein targets. Moreover, the study of viral MBPs can also help to understand the roles of MBPs in virus pathogenesis and virus-host interactions.

SEQUENCE-FREE PHYLOGENETICS WITH MASS SPECTROMETRY

An alternative, more rapid, sequence-free approach to build phylogenetic trees has been conceived and implemented. Molecular phylogenetics has continued to mostly focus on improvement in tree construction based on gene sequence alignments. Here protein-based phylogenies are constructed using numerical data sets ("phylonumerics") representing the masses of peptide segments recorded in a mass mapping experiment. This truly sequence-free method requires no gene sequences, nor their alignment, to build the trees affording a considerable time and cost-saving to conventional phylogenetics methods. The approach also calculates single point amino acid mutations from a comparison of mass pairs from different maps in the data set and displays these at branch nodes across the tree together with their frequency. Studies of the consecutive, and near-consecutive, ancestral and descendant mutations across interconnected branches of a mass tree allow putative adaptive, epistatic, and compensatory mutations to be identified in order to investigate mechanisms associated with evolutionary processes and pathways. A side-by-side comparison of this sequence-free approach and conventional gene sequence phylogenetics is discussed.

Protein phylogenetics with mass spectrometry. A comparison of methods

Advances in protein mass spectrometry have provided the ability to identify and sequence proteins with unprecedented speed, sensitivity and accuracy. These benefits now offer advantages for studies of protein evolution and phylogeny avoiding the need to generate and align DNA sequences which can prove time consuming, costly and difficult in the case of large genomes and for highly diverse organisms. The methods of phylogenetic analysis using protein mass spectrometry can be classified into three categories: (1) de novo protein sequencing followed by multiple sequence alignment for classical phylogenetic reconstruction, (2) direct phylogenetic reconstruction using expressed protein mass profiles exploited in microbial biotyping applications, and (3) the construction of trees using proteolytic peptide mass map or fingerprint data. This review describes the three approaches together with the relevant tools and algorithms required to implement them. It also compares each of these alternative protein based methods alongside conventional gene sequence based phylogenetics.

Episodic evolution of coadapted sets of amino acid sites in mitochondrial proteins

The rate of evolution differs between protein sites and changes with time. However, the link between these two phenomena remains poorly understood. Here, we design a phylogenetic approach for distinguishing pairs of amino acid sites that evolve concordantly, i.e., such that substitutions at one site trigger subsequent substitutions at the other; and also pairs of sites that evolve discordantly, so that substitutions at one site impede subsequent substitutions at the other. We distinguish groups of amino acid sites that undergo coordinated evolution and evolve discordantly from other such groups. In mitochondrion-encoded proteins of metazoans and fungi, we show that concordantly evolving sites are clustered in protein structures. By analysing the phylogenetic patterns of substitutions at concordantly and discordantly evolving site pairs, we find that concordant evolution has two distinct causes: epistatic interactions between amino acid substitutions and episodes of selection independently affecting substitutions at different sites. The rate of substitutions at concordantly evolving groups of protein sites changes in the course of evolution, indicating episodes of selection limited to some of the lineages. The phylogenetic positions of these changes are consistent between proteins, suggesting common selective forces underlying them.

The mode and rate of evolution of a protein site depends on the effect of its mutations on protein fitness. The fitness effect of a mutation itself can change in the course of evolution for at least two reasons. First, it can be modulated by substitutions occurring at other sites, a phenomenon called epistasis. Second, changes in selection can be non-epistatic, affecting sites independently of one another. Here, we analyse substitutions accumulated by the evolving lineages of the five proteins encoded by the mitochondrial genomes of thousands of species of metazoans and fungi. We show that substitutions at different amino acid sites occur in a coordinated fashion, and this coordination is caused both by epistasis and by episodes of selection affecting groups of sites. We partition each protein into several groups of concordantly evolving sites such that evolution of sites from different groups is discordant, and show that the proteins encoded by the mitochondrial genome consist of coevolving structural blocks. Some of these blocks have a clear functional specialization, e.g. are associated with interfaces between proteins composing respiratory complexes. Together, our results reveal a previously unrecognized complexity in the causes of variation in evolutionary rates between protein sites.

Existing Host Range Mutations Constrain Further Emergence of RNA Viruses

RNA viruses are capable of rapid host shifting, typically due to a point mutation that confers expanded host range. As additional point mutations are necessary for further expansions, epistasis among host range mutations can potentially affect the mutational neighborhood and frequency of niche expansion. We mapped the mutational neighborhood of host range expansion using three genotypes of the double-stranded RNA (dsRNA) bacteriophage φ6 (wild type and two isogenic host range mutants) on the novel host Pseudomonas syringae pv. atrofaciens. Both Sanger sequencing of 50 P. syringae pv. atrofaciens mutant clones for each genotype and population Illumina sequencing revealed the same high-frequency mutations allowing infection of P. syringae pv. atrofaciens. Wild-type φ6 had at least nine different ways of mutating to enter the novel host, eight of which are in p3 (host attachment protein gene), and 13/50 clones had unchanged p3 genes. However, the two isogenic mutants had dramatically restricted neighborhoods: only one or two mutations, all in p3. Deep sequencing revealed that wild-type clones without mutations in p3 likely had changes in p12 (morphogenic protein), a region that was not polymorphic for the two isogenic host range mutants. Sanger sequencing confirmed that 10/13 of the wild-type φ6 clones had nonsynonymous mutations in p12, and 2 others had point mutations in p9 and p5. None of these genes had previously been associated with host range expansion in φ6. We demonstrate, for the first time, epistatic constraint in an RNA virus due to host range mutations themselves, which has implications for models of serial host range expansion.IMPORTANCE RNA viruses mutate rapidly and frequently expand their host ranges to infect novel hosts, leading to serial host shifts. Using an RNA bacteriophage model system (Pseudomonas phage φ6), we studied the impact of preexisting host range mutations on another host range expansion. Results from both clonal Sanger and Illumina sequencing show that extant host range mutations dramatically narrow the neighborhood of potential host range mutations compared to that of wild-type φ6. This research suggests that serial host-shifting viruses may follow a small number of molecular paths to enter additional novel hosts. We also identified new genes involved in φ6 host range expansion, expanding our knowledge of this important model system in experimental evolution.

Mutation and Epistasis in Influenza Virus Evolution

Influenza remains a persistent public health challenge, because the rapid evolution of influenza viruses has led to marginal vaccine efficacy, antiviral resistance, and the annual emergence of novel strains. This evolvability is driven, in part, by the virus's capacity to generate diversity through mutation and reassortment. Because many new traits require multiple mutations and mutations are frequently combined by reassortment, epistatic interactions between mutations play an important role in influenza virus evolution. While mutation and epistasis are fundamental to the adaptability of influenza viruses, they also constrain the evolutionary process in important ways. Here, we review recent work on mutational effects and epistasis in influenza viruses.