key: cord-0950025-7bzvvizi
authors: Lumata, Jenica L.; Ball, Darby; Shahrivarkevishahi, Arezoo; Luzuriaga, Michael A.; Herbert, Fabian C.; Brohlin, Olivia; Lee, Hamilton; Hagge, Laurel M.; D’Arcy, Sheena; Gassensmith, Jeremiah J.
title: Identification and physical characterization of a spontaneous mutation of the tobacco mosaic virus in the laboratory environment
date: 2021-07-23
journal: Sci Rep
DOI: 10.1038/s41598-021-94561-2
sha: bdd27ca95688bf4adff4a3ee4e617bec6f449be6
doc_id: 950025
cord_uid: 7bzvvizi

Virus-like particles are an emerging class of nano-biotechnology with the Tobacco Mosaic Virus (TMV) having found a wide range of applications in imaging, drug delivery, and vaccine development. TMV is typically produced in planta, and, as an RNA virus, is highly susceptible to natural mutation that may impact its properties. Over the course of 2 years, from 2018 until 2020, our laboratory followed a spontaneous point mutation in the TMV coat protein—first observed as a 30 Da difference in electrospray ionization mass spectrometry (ESI–MS). The mutation would have been difficult to notice by electrophoretic mobility in agarose or SDS-PAGE and does not alter viral morphology as assessed by transmission electron microscopy. The mutation responsible for the 30 Da difference between the wild-type (wTMV) and mutant (mTMV) coat proteins was identified by a bottom-up proteomic approach as a change from glycine to serine at position 155 based on collision-induced dissociation data. Since residue 155 is located on the outer surface of the TMV rod, it is feasible that the mutation alters TMV surface chemistry. However, enzyme-linked immunosorbent assays found no difference in binding between mTMV and wTMV. Functionalization of a nearby residue, tyrosine 139, with diazonium salt, also appears unaffected. Overall, this study highlights the necessity of standard workflows to quality-control viral stocks. We suggest that ESI–MS is a straightforward and low-cost way to identify emerging mutants in coat proteins.

www.nature.com/scientificreports/ nanoparticles, use an extremophile strain of the Australian plant Nicotiana benthamiana for large-scale production of TMV. It is the most widely used host for plant virus replication as it is very susceptible to infection. An RNA silencing gene in N. benthamiana (NbRdRP1m) is mutated and less active, contributing to the species' susceptibility of RNA-virus infection [25] [26] [27] . Single-stranded RNA viruses, like TMV, are very susceptible to mutation 28 . Though the mechanism is not established yet, Sanjuán et al. (2016) suggested that the genetic material in viruses with single-stranded genomes are more exposed to oxidative deamination and other chemical stress 28 . In fact, mutant populations are ubiquitous in plant RNA viruses and their presence is rarely considered significant until the mutant phenotype dominates 29 . RNA viruses, except Coronavirus, contain RNA polymerases that lack 3′ → 5′ exonuclease proofreading mechanisms, which makes them more prone to error than DNA viruses 28, 30 . One study showed that infecting tomato plants containing the Tm-2 Gene with TMV strain Ltbl, resulted in point mutations in the 30 kDa movement protein 31 . TMV is well-known to be susceptible to environmental and evolutionary stressors and the emergence of mutants should not be surprising; however, there is little discussion in the literature of how coat protein mutations might alter the viruses' physical properties vis-à-vis chemical functionalization and antibody binding. Best practices for monitoring TMV mutation and the impacts of TMV mutation on virus usability are not well established. Knowing what to look for, the relative time scale, and how to identify mutation in plant-sourced viruses would be helpful to the broader community.

Here, we report the emergence of a mutant strain of TMV over the course of two years in the absence of any purposeful chemical or environmental stressor. The strain contains a single point mutation in the coat protein and was identified by electrospray ionization mass spectrometry (ESI-MS) of TMV stocks archived from January 2018 to January 2020. Over this period, the mutant strain largely displaced the wild-type strain, suggesting a competitive advantage. While we were following best practices for purifying and testing isolated TMV at the time, we have since identified additional practices that should be followed to identify emerging mutants. In many assays, the mutant was indistinguishable from wild-type. This work highlights the need for standard workflows to quality-control viral stocks in the chemical and bioengineering laboratory, particularly if they are being considered for translational purposes.

Our TMV stock is grown in N. benthamiana eight weeks after germination while the plants are around 11.5 cm tall. N. benthamiana is grown in a purpose-built plant growth room with constant temperature (22 °C) and humidity (67%). TMV (15 mg) from a prior harvest is added to 1 g of silicon carbide as an abrasive and rubbed into the leaves of N. benthamiana. Once the leaves become discolored, two weeks after infection, the leaves are collected and stored at − 80 °C until the virus is extracted as previously described 12 .

A small amount of TMV is always set aside for the next infection. In the literature, the TMV stock is typically characterized by electrophoretic mobility (e.g. SDS-PAGE), transmission electron microscopy (TEM), size exclusion chromatography, and intact protein mass spectrometry 6, 9, 12, 14, 24 . RNA sequencing has been done in some cases 32, 33 but is not routinely reported. If mutant TMV strains have been detected by sequencing, they have not been widely reported in the literature to date.

We monitor our TMV purifications using liquid chromatography mass spectrometry (LC-MS). We conducted LC-MS using ESI-MS on denatured viral samples. 20 µL of 10 mg/mL TMV was mixed with 40 µL of glacial acetic acid and the RNA removed via centrifugation. Over two years, from January 2018 till January 2020, we observed the appearance of a new mutant TMV strain (Fig. 1A) . In Jan 2018, we observed the anticipated coat protein mass of 17,534 Da, the theoretical mass of wild-type TMV coat protein (wTMV) acetylated at its N-terminus 6, 9, 12, 34 . A second peak was also observed at approximately 17,564 Da. This peak was a relatively minor peak but was consistent in every sample of extracted TMV. By Jan 2020, the peak at 17,564 Da, which we now know is a mutant TMV coat protein (mTMV), had become dominant.

Aside from its mass, the mTMV is indistinguishable from wTMV. Characterization approaches that follow best practice per the literature, do not indicate the heterogeneity of the TMV preparations. When the wTMV and mTMV samples are analyzed by electromobility in SDS-PAGE and agarose gels (Fig. 1B) , for example, there is no difference in protein migration. SDS-PAGE shows the typical single protein mass. Agarose, which measures the migration of intact viral nanoparticles, can be affected by both changes in mass, as well as charge, and it also showed no difference between wTMV and mTMV. The intact viral nanoparticles of both mTMV and wTMV ( Fig. 1C ) also show identical rod morphology in TEM. Had we not routinely conducted ESI-MS analyses of our isolates, the mutant would have gone undetected. Matrix assisted laser desorption ionization-time of flight mass spectrometry (MALDI-MS) is widely used to measure protein mass. From our assessment of the literature, most mass analyses for TMV coat proteins are conducted via MALDI-MS, including many from our group 6, 14, 32, 35 . Nominal resolution MALDI is attractive because the instruments are low cost, sample preparation and analyses are straightforward, and the gentle ionization conveniently produces the + 1 m/z peak. However, nominal resolution MALDI-TOF typically lacks the resolving power to identify small changes in mass on a large protein 36 .

To identify the precise mutation(s) responsible for the 30 Da difference between wTMV and mTMV coat protein, we adopted a bottom-up approach. We compared a sample containing 60% wTMV and 40% mTMV, to one almost entirely composed of mTMV ( Fig. 2A) . Tryptic digests of both samples produced seven shared peptides that covered 90% of the coat protein sequence. Of these seven peptides, only one showed over an order of magnitude difference in intensity between the two samples ( Fig. 2B left) . This peptide covered residues 142-158, suggesting the TMV mutation was located within these C-terminal residues of the protein. For further confirmation, peptic digests of both samples were also performed. These similarly identified a change in intensity for a peptide covering residues 151-158 ( Fig. 2B right) . These data show that the mutation is in one of the last 8 residues of the coat protein (Fig. 3B ). www.nature.com/scientificreports/ We next generated a list of all possible single point mutations that would cause a mass increase of 30 ± 1 Da in the last 8 residues of the coat protein. The list contained nine possible sequences and each sequence was used to interrogate the tryptic and peptic digests of the mTMV sample. Only one sequence gave a passing result and identified the C-terminal peptide of the mTMV coat protein. This sequence contained a G155S mutation. In the two C-terminal peptides (142-158 and 151-158), collision-induced dissociation data show a serine at position 155 with consecutive fragmentation products recovered at the precise site of mutation (Fig. 2D ). When using the sequence with G155S, the C-terminal peptides also had more intensity in the mTMV sample compared to the mixed sample (Fig. 2C) . The difference between wTMV and mTMV is thus a G155S point mutation in the coat protein. At the nucleic acid level, the most likely mutation is a single guanine to adenine which changes the codon of residue 155 from GGU (glycine) to AGU (serine).

The mutation of glycine to serine at position 155 would not be expected to drastically impact the folding or stability of either the coat protein or the intact TMV. Glycine-155 is in a loop that is exposed to solvent on the exterior-side of the TMV rod (Fig. 3A) . It is not involved in contacts between the many copies of the coat protein. A glycine to serine change is also somewhat conservative with both having small side chains. Differences may occur, however, as glycine allows for greater backbone flexibility and serine brings an additional hydroxyl www.nature.com/scientificreports/ group for hydrogen bonding. Algorithms predict a slight reduction in protein stability with the G155S mutation (ΔΔG = -0.41) 37 .

Given the surface localization of the G155S mutation, we sought to test its impact on TMV antibody binding and chemical modification. With identical concentrations of TMV, changes in the relative binding of polyclonal antibodies would be obvious by comparing a dilution series in an enzyme-linked immunosorbent assay (ELISA). ELISA run on both mTMV and wTMV bind TMV antibodies at the same concentrations to produce linear curves with a correlation coefficient of 0.998, indicating the binding is almost identical (Fig. 4A ). Therefore, it seems unlikely this mutation would have impacted any immunological or cell studies. So far, the ELISA result shows no significant difference between the binding of the epitopes of mutant and wild-type TMVs to TMV antibodies, but if there is any change in the epitope, it might likely be missed by the method, and further study is needed. TMV is also widely used for its facile chemical modifiability at the external tyrosine residue via a diazonium coupling reaction (Fig. 4B) . We were able to show that the tyrosine residue can be functionalized through diazonium coupling quantitatively on the mTMV. Figure 4C shows the deconvoluted ESI mass spectra of wTMV-Alk at 17,662 Da-which is a 128 Da increase from wTMV mass from the attached azo-alkyne, although very small wTMV peak can still be observed on the spectrum. A preceding peak can be also observed (17,645 Da), which has been previously reported 6, 9 . For mTMV-Alk, a complete conjugation is observed with the major peak at 17,692 Da, a 30 Da increase from wTMV-Alk. The presence of wTMV in the m TMV sample, generated four peaks in total. Based on these data, we conclude that the mutation does not have an obvious effect on the ability to functionalize Y130 via diazonium coupling. The morphology between wTMV-Alk and mTMV-Alk remains the same (Fig. 4D) . These data show that the G155S mutation had no obvious effects on the physical properties of the TMV.

In a span of 2 years, the wild-type TMV has spontaneously mutated in the laboratory environment as observed in the ESI-MS data collected in that period. It should be noted that similar data collected on intact proteins using nominal resolution MALDI-TOF, where the detected ion is the + 1 m/z, would likely have insufficient resolution to detect the presence of single point mutant strains. Further RNA sequencing methods that can detect mutants when they are not the dominant species are complicated and expensive. Bottom-up analysis and peptide sequencing with collision-induced dissociation identified the residue 155 as the site of a glycine to serine point mutation. Dramatic physical effects of this mutation are not predicted but given its surface localization there was the possibility of altered binding and surface chemistry. However, ELISA shows that both mutant and wild-type TMV have identical binding to TMV antibodies and can be equally functionalized at tyrosine 139. It is unclear precisely why the G155S mutant emerged and overtook the wild-type in our laboratory environment. www.nature.com/scientificreports/ Our discovery, however, highlights the necessity to monitor every batch of purified TMV as spontaneous point mutation can occur even in the absence of stress and mutagens. Instrumentation. LC-MS were obtained using an Agilent 1100 HPLC with a PLRP-S column for separation and an AB Sciex 4000 QTRAP system for detection; also, a Waters SYNAPT G2-Si Q-TOF with an M-class UPLC. TEM was conducted using a JEOL JEM-1400Plus transmission electron microscope. Bio-rad ChemiDoc MP was used for gel imaging. ELISA data were obtained using a BioTek Synergy H4 Hybrid microplate reader.

TMV solution with 2 weeks of incubation. The harvested infected leaves were stored at − 80 °C. About 100 g of leaves were blended with cold (4 °C) potassium phosphate (KP) buffer (0.1 M, 1000 mL, pH 7.4) and 2-mercaptoethanol (0.2% (v/v)) was added. It was subsequently ground to have an effective extraction. The slurry was filtered and centrifuged at 11,000 × g (4 °C, 20 min). The resulting supernatant was filtered again. The volume was measured and an equal amount of chloroform/1-butanol with 1-to-1 ratio was mixed (4 °C, 30 min). Another centrifugation was done at 4500 × g for 10 min. The collected aqueous phase was mixed with NaCl (final concentration of 0.2 M), PEG 8000 (8% (w/w)), and Triton X-100 surfactant (1% (w/w) Transmission electron microscopy (TEM). TEM imaging was performed on a JEOL JEM-1400 + transmission electron microscope. Samples were prepared by incubating 5 µL of ~ 0.1 mg/mL TMV in water on a 300 mesh formvar-coated copper grid for 30 s. The sample was then stained with 5 μL 2% uranyl acetate for an additional 30 s. The excess liquid was wicked away with a Whatman (#1) filter paper and the grids were left to air dry. Images were taken with an accelerating voltage of 120 kV.

ESI-MS. TMV samples were prepared by denaturing 20 µL of 10 mg/mL in 40 µL of glacial acetic acid. Sample was then centrifuged at 4300 × g for 10 s to separate the precipitated RNA. The supernatant was collected and run on an Agilent 1100 series HPLC system with a PLRP-S column for separation followed by a 4000 QTRAP mass spectrometer. The flow rate is 0.250 mL/min, and the solvent system comprises of Milli Q water, 0.1% formic acid, and pH 7.0 in a 0.1 M sodium phosphate buffer. This system was used for running both TMV and TMV-Alk samples.

Electrophoretic mobility assays. 1% (w/v) Agarose gels were used. The sample was prepared by mixing 3 µg TMV with 5 µL Thermo Scientific 6× DNA Loading Dye. From that mixture, 4 µL was added to each well. The gel was run at 100 eV for 45 min, stained with coomassie brilliant blue, and visualized using Bio-rad ChemiDoc MP gel imager. 10% SDS-PAGE gel was used. The sample was prepared by mixing 3 µg TMV with 5 µL of SDS loading dye (β-Mercaptoethanol (5%), Bromophenol blue (0.02%), Glycerol (30%), SDS (Sodium dodecyl sulfate 10%), Tris-Cl (250 mM, pH 6.8)) and 5 µL of 0.1 M dithiothreitol. The mixture was boiled for 10 min. 4 µL of sample was added to each well and the gel was run at 100 eV for 45 min, stained with coomassie brilliant blue, and visualized using Bio-rad ChemiDoc MP gel imager. www.nature.com/scientificreports/ was 0.1% formic acid. Peptides were separated over a C18 column (Waters Acquity UPLC BEH) and eluted with a linear 3-40% (v/v) acetonitrile gradient for 7 min at 40 μl min −1 and 1 °C. Mass-spectrometry data were acquired using positive-ion mode in HDMS E mode, collecting both low-energy (6 V) and high-energy (ramping 22-44 V) peptide-fragmentation data for peptide identification. All samples were acquired in resolution mode. Capillary voltage was set to 2.8 kV for the sample sprayer. Desolvation gas was set to 650 l h −1 at 175 °C. The source temperature was set to 80 °C. Cone and nebulizer gas was flowed at 90 l h −1 and 6.5 bar, respectively. The sampling cone and source offset were both set to 30 V. Data were acquired at a scan time of 0.4 s with a m/z range of 100-2000. Mass correction was done using [Glu1]-fibrinopeptide B as a reference mass.

Data processing. Raw data were processed by PLGS (Waters Protein Lynx Global Server 3.0.2) using a database containing Sus scrofa pepsin A and native TMV coat protein. In PLGS, the minimum fragment ion matches per peptide was 3, and methionine oxidation and N-terminal acetylation were allowed. Trypsin samples included cysteine carbamidomethylation (CAM) as a fixed modifier. The low and elevated energy thresholds were 250 and 50 counts, respectively, and the overall intensity threshold was 750 counts. Peptides were curated in DynamX 3.0 with thresholds of 0.3 products per amino acid and one consecutive product. www.nature.com/scientificreports/ Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.

Nanoparticles for Biomedical Applications

Viral chemistry: The chemical functionalization of viral architectures to create new technology

Medicago and GSK Announce Start of Phase 2/3 Clinical Trials of Adjuvanted Covid-19 Vaccine Candidate

The emergency response capacity of plant-based biopharmaceutical manufacturing-what it is and what it could be

Advances in Virus Research

Virus Hybrids as Nanomaterials: Methods and Protocols

From stars to stripes: RNA-directed shaping of plant viral protein templates-Structural synthetic virology for smart biohybrid nanostructures

Viruses and Nanotechnology

Nitroxyl modified tobacco mosaic virus as a metal-free high-relaxivity MRI and EPR active superoxide sensor

Plant virus particles with various shapes as potential adjuvants

Plant virus-based materials for biomedical applications: Trends and prospects

Supramolecular and biomacromolecular enhancement of metal-free magnetic resonance imaging contrast agents

Tobacco mosaic virus delivery of mitoxantrone for cancer therapy

Dysprosium-Modified Tobacco Mosaic Virus Nanoparticles for Ultra-High-Field Magnetic Resonance and Near-Infrared Fluorescence Imaging of Prostate Cancer

Bioengineering of tobacco mosaic virus to create a non-infectious positive control for ebola diagnostic assays

Modified TMV particles as beneficial scaffolds to present sensor enzymes

Synthesis of metal-organic frameworks on tobacco mosaic virus templates

Enhanced stability and controlled delivery of MOF-encapsulated vaccines and their immunogenic response in vivo

Investigation of controlled growth of metal-organic frameworks on anisotropic virus particles

Genetically engineered tobacco mosaic virus as nanoparticle vaccines

An improved tobacco mosaic virus (TMV)-conjugated multiantigen subunit vaccine against respiratory tularemia

Development of a multivalent subunit vaccine tularemia using tobacco mosaic virus (TMV) based delivery system

Using tobacco mosaic virus epitope display to enhance papillomavirus L2 antigen presentation

Genomic position affects the expression of tobacco mosaic virus movement and coat protein genes

Differential responses to virus challenge of laboratory and wild accessions of australian species of nicotiana, and comparative analysis of RDR1 gene sequences

A natural variant of a host RNA-dependent RNA polymerase is associated with increased susceptibility to viruses by nicotiana benthamiana

Nicotiana benthamiana: Its history and future as a model for plantpathogen interactions

Mechanisms of viral mutation

Lack of evidence for proofreading mechanisms associated with an RNA virus polymerase

Mutations in the tobacco mosaic virus 30-kD protein gene overcome Tm-2 resistance in tomato

Tobacco mosaic virus-based protein nanoparticles and nanorods for chemotherapy delivery targeting breast cancer

Tobacco mosaic virus-directed reprogramming of auxin/indole acetic acid protein transcriptional responses enhances virus phloem loading

Tobacco mosaic virus assembly and disassembly: Determinants in pathogenicity and resistance

Regulating the uptake of viral nanoparticles in macrophage and cancer cells via a pH switch

General FAQ About Mass Spectrometry

SDM: A server for predicting effects of mutations on protein stability

The authors declare no competing interests.

The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-021-94561-2.Correspondence and requests for materials should be addressed to S.D. or J.J.G.Reprints and permissions information is available at www.nature.com/reprints.