Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The VSports app下载. gov means it’s official. Federal government websites often end in . gov or . mil. Before sharing sensitive information, make sure you’re on a federal government site. .

Https

The site is secure. The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely. V体育官网.

. 2017 Sep 13;22(3):377-386.e5.
doi: 10.1016/j.chom.2017.08.004.

Epitranscriptomic Enhancement of Influenza A Virus Gene Expression and Replication

Affiliations

"V体育官网入口" Epitranscriptomic Enhancement of Influenza A Virus Gene Expression and Replication

David G Courtney et al. Cell Host Microbe. .

Abstract

Many viral RNAs are modified by methylation of the N6 position of adenosine (m6A). m6A is thought to regulate RNA splicing, stability, translation, and secondary structure. Influenza A virus (IAV) expresses m6A-modified RNAs, but the effects of m6A on this segmented RNA virus remain unclear VSports手机版. We demonstrate that global inhibition of m6A addition inhibits IAV gene expression and replication. In contrast, overexpression of the cellular m6A "reader" protein YTHDF2 increases IAV gene expression and replication. To address whether m6A residues modulate IAV RNA function in cis, we mapped m6A residues on the IAV plus (mRNA) and minus (vRNA) strands and used synonymous mutations to ablate m6A on both strands of the hemagglutinin (HA) segment. These mutations inhibited HA mRNA and protein expression while leaving other IAV mRNAs and proteins unaffected, and they also resulted in reduced IAV pathogenicity in mice. Thus, m6A residues in IAV transcripts enhance viral gene expression. .

Keywords: RNA methylation; YTHDF2; epitranscriptomics; influenza A virus; m(6)A; mRNA function; post-transcriptional regulation; viral pathogenesis V体育安卓版. .

PubMed Disclaimer

Figures

Figure 1
Figure 1. IAV gene expression is greatly reduced in METTL3 knockout cell lines
(A) Treatment of cells with a non-toxic dose of DAA, an inhibitor of m6A addition, reduced expression of the IAV proteins NS1 and M2 in A549 cells, as determined by Western blot. A549 cells were infected at an MOI of 1.0. (B) Two clonal A549 METTL3 knockout cell lines, M3.1 and M3.2, were established using gene editing and confirmed by Western blot and genomic sequencing. See Fig. S1 for genomic sequences and Table S1 for sgRNA sequences. (C) Two A549 METTL3 knockout lines and two control cell lines, including the parental A549 cell line and a GFP-Flag expressing A549 cell line, were infected with IAV-PR8 at an MOI of 0.01 and NP, NS1 and M2 protein expression determined by Western blot at 24, 48 and 72 hours post infection (hpi). Quantification of the band intensities for NS1 and M2 at 72 hpi revealed expression levels of 0.12 ± 0.01 and 0.12 ± 0.07 for M3.1 and 0.09 ± 0.03 and 0.08 ± 0.03 for M3.2, respectively, when expression in wild type cells was normalized to 1.0. Actin was used as the loading control. (D) Quantitative RT-PCR (qRT-PCR) was performed to determine the levels of the spliced IAV M2 mRNA at the same time points post-infection. See Table S2 for primer sequences. (E) Viral production from the wild type and METTL3 KO cells at 72 hpi was quantified by plaque assay on MDCK cells These data represent the average of three biological replicates with SD indicated (*=p<0.05, **=p<0.01).
Figure 2
Figure 2. Overexpression of YTHDF2 increases all aspects of viral gene expression
(A) A549 cell lines overexpressing YTHDF1-Flag (Y1), YTHDF2-Flag (Y2.1 and Y2.2) or GFP were generated by lentiviral transduction followed by single cell cloning. (B) A549 cell lines were infected with IAV-PR8 at an MOI of 0.01 and expression of the viral proteins NP, NS1 and M2 assessed by Western blot at 24, 48 and 72 hpi. GFP, YTHDF1 and YTHDF2 were detected using anti-Flag. The parental A549 cell line, and A549 cells expressing GFP-Flag, were used as controls. Quantification of the band intensities for NS1 and M2 at 72 hpi revealed expression levels of 6.28 ± 1.28 and 5.19 ± 2.13 for Y2.1 and 5.70 ± 1.65 and 5.84 ± 1.73 for Y2.2, respectively, when the expression in wild type cells was normalized to 1.0. Actin was used as the loading control. Single round infections are described in Fig. S2. (C) qRT-PCR was used to determine the mRNA levels of the spliced M2 IAV mRNA at the same time points post-infection. See Table S2 for primer sequences. (D) The viral titer produced from these cell lines at 72hpi was determined by plaque assay on MDCK cells. These data represent the average of three biological replicates with SD indicated (*=p<0.05, **=p<0.01).
Figure 3
Figure 3. Identification of m6A sites on IAV-PR8 plus sense mRNA
PAR-CLIP and PA-m6A-seq were performed using 293T cell lines 24 hpi with IAV-PR8 at an MOI of 10 (A) Concatenated map of the IAV-PR8 transcriptome that reads were aligned to. (B) Complete transcriptome coverage tracks are shown for PAR-CLIP performed on Flag-GFP, Flag-YTHDF1, Flag-YTHDF2 and Flag-YTHDF3 expressing 293T cells, while PA-m6A-seq was performed using wild type 293T cells. The PA-m6A-seq lane has a Y axis of 0–500 reads, and all others are depicted with Y axes of 0–200 reads. (C) An expanded view of the HA segment of IAV-PR8, with 8 prominent m6A sites numbered. PA-m6A-seq has a Y axis of 0–250 reads, and all others are depicted with Y axes of 0–100 reads. Reverse transcription of crosslinked 4SU residues results in characteristic T>C mutations and the level of T>C conversion at specific residues is indicated by red/blue bars.
Figure 4
Figure 4. Identification of m6A sites on IAV-PR8 minus strand vRNA
PAR-CLIP and PA-m6A-seq were performed on RNA from 293T cell lines 24 hpi with IAV-PR8 at MOI of 10. (A) Concatenated map of the IAV-PR8 genome used for read alignment. (B) Complete genomic coverage tracks are shown, similar to Fig. 2. The PA-m6A-seq lane has a Y axis of 0–2000 reads, and all others are depicted with Y axes of 0–300 reads. (C) An expanded view of the HA segment of IAV-PR8, with 9 prominent m6A sites numbered. PA-m6A-seq has a Y axis of 0–1000 reads, and all others are depicted with Y axes of 0–150 reads. Green/brown bars indicate the level of A>G conversion at specific A residues.
Figure 5
Figure 5. Mutagenesis of m6A motifs on the IAV HA segment inhibits HA mRNA and protein expression
Mutations introduced into the HA segment are described in Fig. S3. (A) Viral RNA extracted from purified IAV-PR8 virions grown in embryonated chicken eggs was separated on a TBE-urea gel and stained. (B) Protein was extracted from purified IAV and Western blotting used to determine HA protein levels in wild type or mutant virions. Band intensity quantification revealed HA0 and HA1 levels of 0.8 ± 0.63 and 0.8 ± 0.54 for the +Mut virions and 0.5 ± 0.53 and 0.6 ± 0.48 for the −Mut virions, respectively, when HA levels in wild type virions was normalized to 1.0. The viral protein M1 was used as a loading control. (C) A multicycle spreading infection was initiated by infection of A549 at an MOI of 0.01 using all 3 IAV-PR8 variants. Total RNA was extracted at 24, 48 and 72 hpi and qRT-PCR then used to determine HA mRNA levels, relative to a GAPDH mRNA internal control. See Table S2 for primer sequences. (D) The same total RNA samples used in panel C were used to quantify M2 mRNA levels, again using GAPDH mRNA as an internal control. (E) Similar to panel C except that Western blotting was performed to evaluate the level of expression of the IAV HA, NS1 and M2 proteins at 24, 48 and 72 hpi. Quantification of HA band intensity at 72 hpi revealed expression levels of 0.69 ± 0.11 for the +Mut virus and 0.32 ± 0.13 for the −Mut virus, when normalized to M2 expression levels. Actin was used as the loading control. Viral spread was also determined by flow cytometry (see Fig. S4). (F) A549 cells were infected at an MOI of 1.0 with the three IAV-PR8 variants and total RNA extracted 12 hpi. qRT-PCR was used to determine the levels of the cellular mRNAs encoding RIG-I, MDA5 and IFNβ1. All data are drawn from three biological replicates with SD indicated (**=p<0.01).
Figure 6
Figure 6. IAV mutants depleted for m6A on the viral HA strand show reduced pathogenicity in vivo
Wildtype C57BL/6 mice were infected with 50 plaque forming units (PFU) of WT PR8, +Mut or −Mut and monitored for (A) body weight with 80% of the starting body weight as a humane endpoint, as indicated by the dashed line and (C) mortality. Similarly, mice were infected with 10 PFU of WT PR8, +Mut or −Mut and monitored for (B) body weight and (D) mortality. Five mice were used per group. Significance is indicated on the graph based on the log-rank Mantel-Cox test.

References

    1. Bader JP, Brown NR, Chiang PK, Cantoni GL. 3-Deazaadenosine, an inhibitor of adenosylhomocysteine hydrolase, inhibits reproduction of Rous sarcoma virus and transformation of chick embryo cells. Virology. 1978;89:494–505. - PubMed
    1. Batista PJ, Molinie B, Wang J, Qu K, Zhang J, Li L, Bouley DM, Lujan E, Haddad B, Daneshvar K, et al. m(6)A RNA modification controls cell fate transition in mammalian embryonic stem cells. Cell Stem Cell. 2014;15:707–719. - "VSports注册入口" PMC - PubMed
    1. Boden D, Pusch O, Lee F, Tucker L, Ramratnam B. Human immunodeficiency virus type 1 escape from RNA interference. J Virol. 2003;77:11531–11535. - PMC - PubMed
    1. Bray M, Driscoll J, Huggins JW. Treatment of lethal Ebola virus infection in mice with a single dose of an S-adenosyl-L-homocysteine hydrolase inhibitor. Antiviral Res. 2000;45:135–147. - "V体育安卓版" PubMed
    1. Chen K, Lu Z, Wang X, Fu Y, Luo GZ, Liu N, Han D, Dominissini D, Dai Q, Pan T, et al. High-resolution N(6) -methyladenosine (m(6) A) map using photo-crosslinking-assisted m(6) A sequencing. Angew Chem Int Ed Engl. 2015;54:1587–1590. - "V体育官网入口" PMC - PubMed

"V体育平台登录" MeSH terms