Discovery and functional interrogation of SARS-CoV-2 protein-RNA interactions
This article has been Reviewed by the following groups
Listed in
- Evaluated articles (ScreenIT)
Abstract
The COVID-19 pandemic is caused by severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The betacoronvirus has a positive sense RNA genome which encodes for several RNA binding proteins. Here, we use enhanced crosslinking and immunoprecipitation to investigate SARS-CoV-2 protein interactions with viral and host RNAs in authentic virus-infected cells. SARS-CoV-2 proteins, NSP8, NSP12, and nucleocapsid display distinct preferences to specific regions in the RNA viral genome, providing evidence for their shared and separate roles in replication, transcription, and viral packaging. SARS-CoV-2 proteins expressed in human lung epithelial cells bind to 4773 unique host coding RNAs. Nine SARS-CoV-2 proteins upregulate target gene expression, including NSP12 and ORF9c, whose RNA substrates are associated with pathways in protein N-linked glycosylation ER processing and mitochondrial processes. Furthermore, siRNA knockdown of host genes targeted by viral proteins in human lung organoid cells identify potential antiviral host targets across different SARS-CoV-2 variants. Conversely, NSP9 inhibits host gene expression by blocking mRNA export and dampens cytokine productions, including interleukin-1α/β. Our viral protein-RNA interactome provides a catalog of potential therapeutic targets and offers insight into the etiology of COVID-19 as a safeguard against future pandemics.
Article activity feed
-
SciScore for 10.1101/2022.02.21.481223: (What is this?)
Please note, not all rigor criteria are appropriate for all manuscripts.
Table 1: Rigor
Ethics IRB: This study protocol was approved by the Institutional Review Board of UCSD’s Human Research Protections Program (181180). Sex as a biological variable not detected. Randomization The foreground was a bed file of significant IDR peaks; the background was randomly defined peaks within the same annotated region as the foreground peaks. Blinding not detected. Power Analysis not detected. Cell Line Authentication Authentication: Cell culture and cell line generation: BEAS-2B, HEK293T and Vero E6 cells were purchased from the American Type Culture Collection and were not further authenticated.
Contamination: Cells were routinely tested for mycoplasma contamination with a MycoAlert mycoplasma test kit …SciScore for 10.1101/2022.02.21.481223: (What is this?)
Please note, not all rigor criteria are appropriate for all manuscripts.
Table 1: Rigor
Ethics IRB: This study protocol was approved by the Institutional Review Board of UCSD’s Human Research Protections Program (181180). Sex as a biological variable not detected. Randomization The foreground was a bed file of significant IDR peaks; the background was randomly defined peaks within the same annotated region as the foreground peaks. Blinding not detected. Power Analysis not detected. Cell Line Authentication Authentication: Cell culture and cell line generation: BEAS-2B, HEK293T and Vero E6 cells were purchased from the American Type Culture Collection and were not further authenticated.
Contamination: Cells were routinely tested for mycoplasma contamination with a MycoAlert mycoplasma test kit (Lonza) and were found negative for mycoplasma.Table 2: Resources
Antibodies Sentences Resources The remaining lysates were immunoprecipitated using 15 µl anti-Strep or 10 µl anti-FLAG antibody (depending on the epitope tag of the construct; Supplementary Table 4) on Sheep Anti-Mouse IgG Dynabeads M-280 (ThermoFisher) overnight at 4°C. anti-FLAGsuggested: NoneAnti-Mouse IgGsuggested: (Thermo Fisher Scientific Cat# 11201D, RRID:AB_2783640)Negative control samples are wild type (WT) BEAS-2B cells, and performed using both anti-Strep and anti-FLAG antibodies (separately). anti-Strepsuggested: NoneCells were fixed with 4% paraformaldehyde for 30 min, which inactivates the virus, before transferring from BSL3 to BSL2, and proceeding with immunofluorescence staining using anti-Nucleocapsid antibody (40143-R019, Sino Biological). BSL2suggested: Noneanti-Nucleocapsidsuggested: NoneNext, cells were incubated with primary antibodies (Supplementary Table 4) at 1:250-2000 dilutions in blocking buffer for 16 h at 4 °C, washed with PBS+0.01% Triton X-100 three times for 5 min each at room temperature, and then incubated with secondary antibody (goat anti-rabbit secondary IgG (H+L) anti-rabbitsuggested: NoneExperimental Models: Cell Lines Sentences Resources Cell culture and cell line generation: BEAS-2B, HEK293T and Vero E6 cells were purchased from the American Type Culture Collection and were not further authenticated. HEK293Tsuggested: KCB Cat# KCB 200744YJ, RRID:CVCL_0063)The ACE2-overexpressing A549 cell line (A549-ACE2) was clonally generated and a gift from Benjamin tenOever52. A549suggested: NoneBEAS-2B cells were cultured on Matrigel (Corning) coated plates and maintained in the PneumaCult-Ex Plus Medium (Stem Cell Technologies), supplemented with 33 µg/ml hydrocortisone (Stem Cell Technologies). BEAS-2Bsuggested: NoneHEK293T, Vero E6 and A549-ACE2 cells were cultured in DMEM (ThermoFisher) supplemented with 10% FBS (ThermoFisher) and passaged every three days. Vero E6suggested: NoneA549-ACE2suggested: NoneSARS-CoV-2 virus infection: SARS-CoV-2 isolates USA-WA1/2020 (BEI Resources, #NR-52281), hCoV-19/USA/CA_UCSD_5574/2020 (lineage B.1.1.7) and hCoV-19/South Africa/KRISP-K005325/2020 (lineage B.1.351 BEI Resources NR-54009) were propagated and infectious units quantified by plaque assay using Vero E6 cells for the WA1 variant, and Vero-TMPRSS2 cells for variants B.1.1.7 and B.1.351. Vero-TMPRSS2suggested: JCRB Cat# JCRB1818, RRID:CVCL_YQ48)Experimental Models: Organisms/Strains Sentences Resources Cell culture and cell line generation: BEAS-2B, HEK293T and Vero E6 cells were purchased from the American Type Culture Collection and were not further authenticated. BEAS-2Bsuggested: NoneRecombinant DNA Sentences Resources Plasmid construction: 2xStrep-tagged plasmids in a pLVX vector expressing SARS-CoV-2 proteins were a gift from Nevan Krogan1. pLVXsuggested: RRID:Addgene_174088)Cloning into the pcDNA3.4 was performed using FastDigest restriction enzymes EcoRI and BshT1 (Invitrogen) and Gibson assembly (NEB). pcDNA3.4suggested: RRID:Addgene_131198)SARS-CoV-2 ORFs were amplified by PCR (KAPA HiFi HotStart ReadyMix, Roche) from the 2xStrep-tagged or 3xFLAG-tagged plasmids with oligonucleotide primers containing attB recombination sites and recombined into pDONR221 using BP clonase II (ThermoFisher) (Supplementary Table 6). pDONR221suggested: RRID:Addgene_63743)Software and Algorithms Sentences Resources This study protocol was approved by the Institutional Review Board of UCSD’s Human Research Protections Program (181180). Human Research Protections Programsuggested: NoneFor reads mapped to the SARS-CoV-2 genome, bedgraph densities were generated using SAM tools v1.9 to obtain read densities at each nucleotide position. SAMsuggested: (SAM, RRID:SCR_010951)For reads mapped to the African Green Monkey genome, peaks were called on the usable reads by CLIPper54 and assigned to gene annotations in Ensembl ChlSab1.1 release 102 and Sars_cov_2 ASM985889 v3.101 were used to annotate peaks mapped to the African Green Monkey and SARS-CoV-2 genome. Ensemblsuggested: (Ensembl, RRID:SCR_002344)Gene Ontology analysis of eCLIP target genes was performed using ENRICHR (https://maayanlab.cloud/Enrichr/https://maayanlab.cloud/Enrichr/). ENRICHRsuggested: (Enrichr, RRID:SCR_001575)Cluster maps were visualized using Cytoscape version 3.8.1 Cytoscapesuggested: (Cytoscape, RRID:SCR_003032)GENCODE v19 gene annotations and featureCounts (v.1.5.30) were used to create read count matrices. featureCountssuggested: (featureCounts, RRID:SCR_012919)Sequencing reads are first processed as RNA-seq libraries, where RNA-seq reads were trimmed of adaptor sequences using cutadapt (v1.4.0) and mapped to repetitive elements (RepBase v18.04) using STAR (v2.4.0i). STARsuggested: (STAR, RRID:SCR_004463)Multiple sequence alignment was performed using MAFFT v7.453 and default parameters, and sequence alignment was visualized using Jalview (version 1.0). MAFFTsuggested: (MAFFT, RRID:SCR_011811)Jalviewsuggested: (Jalview, RRID:SCR_006459)De novo motif analysis: HOMER was used to identify de novo motifs using reads from IDR peaks. HOMERsuggested: (HOMER, RRID:SCR_010881)The length of each region within the metagene was then scaled to 8%, 62% and 30%, corresponding to the average length of regions from the most highly expressed transcripts in ENCODE HepG2 RNA-seq control datasets. ENCODEsuggested: (Encode, RRID:SCR_015482)Interval features such as introns and exons were extracted from the Gencode v19 coordinates. Gencodesuggested: (GENCODE, RRID:SCR_014966)Images were collected via Zeiss ZEN software and converted to tiff for downstream analysis. Zeiss ZENsuggested: NoneImages were analyzed using a custom-developed pipeline in CellProfiler (v.3.1.09). CellProfilersuggested: NoneResults from OddPub: Thank you for sharing your code and data.
Results from LimitationRecognizer: An explicit section about the limitations of the techniques employed in this study was not found. We encourage authors to address study limitations.Results from TrialIdentifier: No clinical trial numbers were referenced.
Results from Barzooka: We did not find any issues relating to the usage of bar graphs.
Results from JetFighter: Please consider improving the rainbow (“jet”) colormap(s) used on page 64. At least one figure is not accessible to readers with colorblindness and/or is not true to the data, i.e. not perceptually uniform.
Results from rtransparent:- Thank you for including a conflict of interest statement. Authors are encouraged to include this statement when submitting to a journal.
- Thank you for including a funding statement. Authors are encouraged to include this statement when submitting to a journal.
- No protocol registration statement was detected.
Results from scite Reference Check: We found no unreliable references.
-