Systematic Exploration of SARS-CoV-2 Adaptation to Vero E6, Vero E6/TMPRSS2, and Calu-3 Cells

Abstract Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to spread globally, and scientists around the world are currently studying the virus intensively in order to fight against the on-going pandemic of the virus. To do so, SARS-CoV-2 is typically grown in the lab to generate viral stocks for various kinds of experimental investigations. However, accumulating evidence suggests that such viruses often undergo cell culture adaptation. Here, we systematically explored cell culture adaptation of two SARS-CoV-2 variants, namely the B.1.36.16 variant and the AY.30 variant, a sub lineage of the B.1.617.2 (Delta) variant, propagated in three different cell lines, including Vero E6, Vero E6/TMPRSS2, and Calu-3 cells. Our analyses detected numerous potential cell culture adaptation changes scattering across the entire virus genome, many of which could be found in naturally circulating isolates. Notable ones included mutations around the spike glycoprotein's multibasic cleavage site, and the Omicron-defining H655Y mutation on the spike glycoprotein, as well as mutations in the nucleocapsid protein's linker region, all of which were found to be Vero E6-specific. Our analyses also identified deletion mutations on the non-structural protein 1 and membrane glycoprotein as potential Calu-3-specific adaptation changes. S848C mutation on the non-structural protein 3, located to the protein's papain-like protease domain, was also identified as a potential adaptation change, found in viruses propagated in all three cell lines. Our results highlight SARS-CoV-2 high adaptability, emphasize the need to deep-sequence cultured viral samples when used in intricate and sensitive biological experiments, and illustrate the power of experimental evolutionary study in shedding lights on the virus evolutionary landscape.

[1]  A. Alam,et al.  Unique mutations in SARS-CoV-2 Omicron subvariants' non-spike proteins: Potential impacts on viral pathogenesis and host immune evasion , 2022, Microbial Pathogenesis.

[2]  F. De Maio,et al.  Rapid Detection of the Omicron (B.1.1.529) SARS-CoV-2 Variant Using a COVID-19 Diagnostic PCR Assay , 2022, Microbiology spectrum.

[3]  D. Guo,et al.  Evidence of Infection of Human Embryonic Stem Cells by SARS-CoV-2 , 2022, Frontiers in Cellular and Infection Microbiology.

[4]  E. Criscuolo,et al.  Proper Selection of In Vitro Cell Model Affects the Characterization of the Neutralizing Antibody Response against SARS-CoV-2 , 2022, Viruses.

[5]  I. Terenin,et al.  Clinically observed deletions in SARS‐CoV‐2 Nsp1 affect its stability and ability to inhibit translation , 2022, FEBS letters.

[6]  J. Gohda,et al.  SARS-CoV-2 Omicron spike H655Y mutation is responsible for enhancement of the endosomal entry pathway and reduction of cell surface entry pathways , 2022, bioRxiv.

[7]  M. Bathe,et al.  Secondary structural ensembles of the SARS-CoV-2 RNA genome in infected cells , 2022, Nature Communications.

[8]  F. Kirchhoff,et al.  Omicron: What Makes the Latest SARS-CoV-2 Variant of Concern So Concerning? , 2022, Journal of virology.

[9]  Haihai Jiang,et al.  Potential Inhibitors Targeting Papain-Like Protease of SARS-CoV-2: Two Birds With One Stone , 2022, Frontiers in Chemistry.

[10]  C. Agoti,et al.  Optimization of the SARS-CoV-2 ARTIC Network V4 Primers and Whole Genome Sequencing Protocol , 2022, Frontiers in Medicine.

[11]  B. La Scola,et al.  Emerging SARS-CoV-2 Genotypes Show Different Replication Patterns in Human Pulmonary and Intestinal Epithelial Cells , 2021, Viruses.

[12]  Kimberly J. Stemple,et al.  Propagation of SARS-CoV-2 in Calu-3 Cells to Eliminate Mutations in the Furin Cleavage Site of Spike , 2021, Viruses.

[13]  Y. Kawaoka,et al.  Enhanced fusogenicity and pathogenicity of SARS-CoV-2 Delta P681R mutation , 2021, Nature.

[14]  J. Barrett,et al.  Variation at Spike position 142 in SARS-CoV-2 Delta genomes is a technical artifact caused by dropout of a sequencing amplicon , 2021, medRxiv.

[15]  M. Farzan,et al.  Mechanisms of SARS-CoV-2 entry into cells , 2021, Nature reviews. Molecular cell biology.

[16]  M. Takeda Proteolytic activation of SARS‐CoV‐2 spike protein , 2021, Microbiology and immunology.

[17]  E. Domingo,et al.  Mutation Rates, Mutation Frequencies, and Proofreading-Repair Activities in RNA Virus Genetics , 2021, Viruses.

[18]  Nicholas T. Ingolia,et al.  The N-terminal domain of SARS-CoV-2 nsp1 plays key roles in suppression of cellular gene expression and preservation of viral gene expression , 2021, Cell Reports.

[19]  S. Ovchinnikov,et al.  ColabFold: making protein folding accessible to all , 2022, Nature Methods.

[20]  Haibo Wu,et al.  Nucleocapsid mutations R203K/G204R increase the infectivity, fitness, and virulence of SARS-CoV-2 , 2021, Cell Host & Microbe.

[21]  Oriol Vinyals,et al.  Highly accurate protein structure prediction with AlphaFold , 2021, Nature.

[22]  O. Pybus,et al.  Assignment of epidemiological lineages in an emerging pandemic using the pangolin tool , 2021, Virus evolution.

[23]  M. Cotten,et al.  Spike Protein Cleavage-Activation Mediated by the SARS-CoV-2 P681R Mutation: A Case-Study From Its First Appearance in Variant of Interest (VOI) A.23.1 Identified in Uganda , 2021, SSRN Electronic Journal.

[24]  S. Boulant,et al.  TMPRSS2 expression dictates the entry route used by SARS‐CoV‐2 to infect host cells , 2021, The EMBO journal.

[25]  M. Galiano,et al.  A Sanger sequencing protocol for SARS‐CoV‐2 S‐gene , 2021, Influenza and other respiratory viruses.

[26]  Ryan V. Moriarty,et al.  A cautionary perspective regarding the isolation and serial propagation of SARS-CoV-2 in Vero cells , 2021, NPJ vaccines.

[27]  L. Hoang,et al.  Mutations in emerging variant of concern lineages disrupt genomic sequencing of SARS-CoV-2 clinical specimens , 2021, International Journal of Infectious Diseases.

[28]  R. Guiomar,et al.  Mutation rate of SARS-CoV-2 and emergence of mutators during experimental evolution , 2021, bioRxiv.

[29]  M. Giacca,et al.  The furin cleavage site in the SARS-CoV-2 spike protein is required for transmission in ferrets , 2021, Nature Microbiology.

[30]  W. P. Duprex,et al.  Recurrent deletions in the SARS-CoV-2 spike glycoprotein drive antibody escape , 2021, Science.

[31]  Natacha S. Ogando,et al.  Genomic monitoring of SARS-CoV-2 uncovers an Nsp1 deletion variant that modulates type I interferon response , 2021, Cell Host & Microbe.

[32]  M. Saier,et al.  The SARS-Coronavirus Infection Cycle: A Survey of Viral Membrane Proteins, Their Functional Interactions and Pathogenesis , 2021, International journal of molecular sciences.

[33]  N. Wu,et al.  Human airway cells prevent SARS-CoV-2 multibasic cleavage site cell culture adaptation , 2021, bioRxiv.

[34]  Lijun Rong,et al.  Targeting SARS‐CoV‐2 viral proteases as a therapeutic strategy to treat COVID‐19 , 2021, Journal of medical virology.

[35]  R. Rottier,et al.  SARS-CoV-2 entry into human airway organoids is serine protease-mediated and facilitated by the multibasic cleavage site , 2021, eLife.

[36]  Xiangxi Wang,et al.  The architecture of the SARS-CoV-2 RNA genome inside virion , 2020, Nature communications.

[37]  M. Islam,et al.  Structure and dynamics of membrane protein in SARS-CoV-2 , 2020, Journal of biomolecular structure & dynamics.

[38]  D. Lavillette,et al.  The SARS-CoV-2 envelope and membrane proteins modulate maturation and retention of the spike protein, allowing assembly of virus-like particles , 2020, Journal of Biological Chemistry.

[39]  Weijin Huang,et al.  Cathepsin L plays a key role in SARS-CoV-2 infection in humans and humanized mice and is a promising target for new drug development , 2020, Signal Transduction and Targeted Therapy.

[40]  S. Olsen,et al.  Activity profiling and crystal structures of inhibitor-bound SARS-CoV-2 papain-like protease: A framework for anti–COVID-19 drug design , 2020, Science Advances.

[41]  Tsuyoshi Sekizuka,et al.  Disentangling primer interactions improves SARS-CoV-2 genome sequencing by multiplex tiling PCR , 2020, PloS one.

[42]  Meitian Wang,et al.  Crystal structure of SARS-CoV-2 papain-like protease , 2020, Acta Pharmaceutica Sinica B.

[43]  Gerbrand J. van der Heden van Noort,et al.  Mechanism and inhibition of the papain‐like protease, PLpro, of SARS‐CoV‐2 , 2020, The EMBO journal.

[44]  F. Klawonn,et al.  Rapid SARS-CoV-2 Adaptation to Available Cellular Proteases , 2020, Journal of virology.

[45]  Y. Xiong,et al.  Nonstructural Protein 1 of SARS-CoV-2 Is a Potent Pathogenicity Factor Redirecting Host Protein Synthesis Machinery toward Viral RNA , 2020, bioRxiv.

[46]  P. Rocchi,et al.  Coronavirus RNA Proofreading: Molecular Basis and Therapeutic Targeting , 2020, Molecular Cell.

[47]  G. Hummer,et al.  Papain-like protease regulates SARS-CoV-2 viral spread and innate immunity , 2020, Nature.

[48]  D. Matthews,et al.  Characterisation of the transcriptome and proteome of SARS-CoV-2 reveals a cell passage induced in-frame deletion of the furin-like cleavage site from the spike glycoprotein , 2020, Genome Medicine.

[49]  N. Ban,et al.  SARS-CoV-2 Nsp1 binds ribosomal mRNA channel to inhibit translation , 2020, bioRxiv.

[50]  N. Loman,et al.  Identification of Common Deletions in the Spike Protein of Severe Acute Respiratory Syndrome Coronavirus 2 , 2020, Journal of Virology.

[51]  A. McElroy,et al.  SARS-CoV-2 growth, furin-cleavage-site adaptation and neutralization using serum from acutely infected, hospitalized COVID-19 patients , 2020, bioRxiv.

[52]  P. Tuchinda,et al.  High-content screening of Thai medicinal plants reveals Boesenbergia rotunda extract and its component Panduratin A as anti-SARS-CoV-2 agents , 2020, Scientific Reports.

[53]  S. Pascarella,et al.  Sars-CoV-2 Envelope and Membrane Proteins: Structural Differences Linked to Virus Characteristics? , 2020, BioMed research international.

[54]  Thomas Becker,et al.  Structural basis for translational shutdown and immune evasion by the Nsp1 protein of SARS-CoV-2 , 2020, Science.

[55]  Sunil Thomas The Structure of the Membrane Protein of SARS-CoV-2 Resembles the Sugar Transporter SemiSWEET , 2020, Pathogens & immunity.

[56]  M. Hoffmann,et al.  A Multibasic Cleavage Site in the Spike Protein of SARS-CoV-2 Is Essential for Infection of Human Lung Cells , 2020, Molecular Cell.

[57]  A. Pruijssers,et al.  The coronavirus proofreading exoribonuclease mediates extensive viral recombination , 2020, bioRxiv.

[58]  Natacha S. Ogando,et al.  SARS-coronavirus-2 replication in Vero E6 cells: replication kinetics, rapid adaptation and cytopathology , 2020, bioRxiv.

[59]  Paul J. Hanson,et al.  Gene expression and in situ protein profiling of candidate SARS-CoV-2 receptors in human airway epithelial cells and lung tissue , 2020, European Respiratory Journal.

[60]  E. Holmes,et al.  The proximal origin of SARS-CoV-2 , 2020, Nature Medicine.

[61]  Fumihiro Kato,et al.  Enhanced isolation of SARS-CoV-2 by TMPRSS2-expressing cells , 2020, Proceedings of the National Academy of Sciences.

[62]  G. Herrler,et al.  SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor , 2020, Cell.

[63]  A. M. Leontovich,et al.  The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2 , 2020, Nature Microbiology.

[64]  B. Graham,et al.  Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation , 2020, Science.

[65]  Kai Zhao,et al.  A pneumonia outbreak associated with a new coronavirus of probable bat origin , 2020, Nature.

[66]  G. Gao,et al.  A Novel Coronavirus from Patients with Pneumonia in China, 2019 , 2020, The New England journal of medicine.

[67]  S. Lo,et al.  A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster , 2020, The Lancet.

[68]  K. To,et al.  Attenuated SARS-CoV-2 variants with deletions at the S1/S2 junction , 2020, Emerging microbes & infections.

[69]  Srinivas Aluru,et al.  Efficient Architecture-Aware Acceleration of BWA-MEM for Multicore Systems , 2019, 2019 IEEE International Parallel and Distributed Processing Symposium (IPDPS).

[70]  Entedar A J Alsaadi,et al.  Membrane binding proteins of coronaviruses , 2019, Future virology.

[71]  Jia Gu,et al.  fastp: an ultra-fast all-in-one FASTQ preprocessor , 2018, bioRxiv.

[72]  Trevor Bedford,et al.  Nextstrain: real-time tracking of pathogen evolution , 2017, bioRxiv.

[73]  Andrey Ziyatdinov,et al.  lme4qtl: linear mixed models with flexible covariance structure for genetic studies of related individuals , 2017, bioRxiv.

[74]  T. Marquès-Bonet,et al.  Bottlenecks and selective sweeps during domestication have increased deleterious genetic variation in dogs , 2015, Proceedings of the National Academy of Sciences.

[75]  A. Wilm,et al.  LoFreq: a sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets , 2012, Nucleic acids research.

[76]  R. Sanjuán From Molecular Genetics to Phylodynamics: Evolutionary Relevance of Mutation Rates Across Viruses , 2012, PLoS pathogens.

[77]  Raul Andino,et al.  Quasispecies Theory and the Behavior of RNA Viruses , 2010, PLoS pathogens.

[78]  Timothy B. Stockwell,et al.  Infidelity of SARS-CoV Nsp14-Exonuclease Mutant Virus Replication Is Revealed by Complete Genome Sequencing , 2010, PLoS pathogens.

[79]  Huanming Yang,et al.  Human Y Chromosome Base-Substitution Mutation Rate Measured by Direct Sequencing in a Deep-Rooting Pedigree , 2009, Current Biology.

[80]  C. Schwegmann-Wessels,et al.  Analysis of ACE2 in polarized epithelial cells: surface expression and function as receptor for severe acute respiratory syndrome-associated coronavirus. , 2006, The Journal of general virology.

[81]  J. Ziebuhr,et al.  Nidovirales: Evolving the largest RNA virus genome , 2006, Virus Research.

[82]  John Bechill,et al.  Identification of Severe Acute Respiratory Syndrome Coronavirus Replicase Products and Characterization of Papain-Like Protease Activity , 2004, Journal of Virology.

[83]  Edward C Holmes,et al.  Error thresholds and the constraints to RNA virus evolution , 2003, Trends in Microbiology.

[84]  E. Domingo,et al.  Lack of evidence for proofreading mechanisms associated with an RNA virus polymerase. , 1992, Gene.

[85]  L. Reed,et al.  A SIMPLE METHOD OF ESTIMATING FIFTY PER CENT ENDPOINTS , 1938 .

[86]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[87]  M. Eigen,et al.  Viral quasispecies. , 1993, Scientific American.

[88]  F. B. Livingstone The founder effect and deleterious genes. , 1969, American journal of physical anthropology.