key: cord-0010751-vhdnwtek
authors: Zhang, Yina; Yu, Tianqi; Ding, Yingnan; Li, Yahui; Lei, Jing; Hu, Boli; Zhou, Jiyong
title: Analysis of Expression Profiles of Long Noncoding RNAs and mRNAs in A549 Cells Infected with H3N2 Swine Influenza Virus by RNA Sequencing
date: 2019-11-27
journal: Virol Sin
DOI: 10.1007/s12250-019-00170-9
sha: a54724f4753c7401fb85000a254b88ef78786d55
doc_id: 10751
cord_uid: vhdnwtek

Long noncoding RNAs (lncRNAs) participate in regulating many biological processes. However, their roles in influenza A virus (IAV) pathogenicity are largely unknown. Here, we analyzed the expression profiles of lncRNAs and mRNAs in H3N2-infected cells and mock-infected cells by high-throughput sequencing. The results showed that 6129 lncRNAs and 50,031 mRNA transcripts in A549 cells displayed differential expression after H3N2 infection compared with mock infection. Among the differentially expressed lncRNAs, 4963 were upregulated, and 1166 were downregulated. Functional annotation and enrichment analysis using gene ontology and Kyoto Encyclopedia of Genes and Genomes databases (KEGG) suggested that target genes of the differentially expressed lncRNAs were enriched in some biological processes, such as cellular metabolism and autophagy. The up- or downregulated lncRNAs were selected and further verified by quantitative real-time polymerase chain reaction (RT-qPCR) and reverse transcription PCR (RT-PCR). To the best of our knowledge, this is the first report of a comparative expression analysis of lncRNAs in A549 cells infected with H3N2. Our results support the need for further analyses of the functions of differentially expressed lncRNAs during H3N2 infection. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s12250-019-00170-9) contains supplementary material, which is available to authorized users.

Swine influenza A virus (SIV) H1N1 caused the 1918 human influenza pandemic, resulting in approximately 50 million deaths (Lemon and Mahmoud 2005) . Additionally, Hong Kong H3N2 flu killed one million people in 1968 -1969 (Simonsen et al. 1998 Dangi and Jain 2012) . From April 16 to April 29, 2018, 1252 (45.4%) specimens were identified by the World Health Organization Global Influenza Surveillance and Response System laboratories as influenza A (H3N2) (http://www.who.int). Moreover, in recent years, SIV has attracted much attention owing to its ability to infect humans. Pigs have surface receptors for both avian and human influenza viruses (Ito et al. 1998) , and influenza viruses from humans and birds can recombine in pigs (Webster et al. 1992) . Among the various subtypes of SIV, H3N2 and H1N1 strains exist in human populations and cause relatively high morbidity and mortality rates (Horimoto and Kawaoka 2005) . Additionally, compared with H1N1, H3N2 infection is associated with human seasonal influenza with a higher epidemic severity index, and its frequency is increasing (Simonsen et al. 1997 ; Thompson et al. 2003) . Although the pathogenesis of the H3N2 subtype has been extensively studied, the underlying molecular mechanisms mediating abnormal host responses are largely unknown.

Traditional studies of host transcriptional responses to pathogen infection have mostly focused on protein-coding genes. We previously reported mitochondrial proteomic analysis of human pulmonary parenchymal cells (A549 cells) infected with swine influenza virus strain H3N2 (Wu et al. 2013) . We found altered expression of 24 host cell proteins (13 upregulated and 11 downregulated proteins) by two-dimensional gel electrophoresis; however, this approach was limited, and further studies are still needed to elucidate the related mechanisms.

Most mammalian genomes are transcribed into noncoding RNAs (ncRNAs), including small ncRNAs (\ 200 bp) and long ncRNAs (lncRNAs; [ 200 bp) (Costa 2010) . With the development of high-throughput sequencing, genetic manipulation, and other biotechnologies, many functions of lncRNAs have been discovered. Studies have shown that lncRNAs play important roles in gene transcription, splicing, protein translation, protein localization, stem cell pluripotency, cell structural integrity, heat shock response, and human diseases (Szczesniak and Makalowska 2016) . Moreover, lncRNAs have also been shown to have roles in innate immune responses (Carpenter et al. 2013; Ilott et al. 2014; Carpenter 2016; Jiang et al. 2018 ) and host-pathogen interactions (Barriocanal et al. 2014; Ouyang et al. 2014; Nishitsuji et al. 2016; Qiu et al. 2018; Tong et al. 2019) . Although thousands of lncRNAs have been identified in mammalian genomes, only a few are associated with innate immunity against viral pathogenesis and influenza viruses. Importantly, lncRNAs induced by influenza virus affect several host biological processes and can also regulate infection. In addition, lncRNAs influence the feedback mode of viral infection by mediating metabolism (Wang et al. 2017) . However, our understanding of the specific functions of lncRNAs in influenza A virus (IAV) infection is still limited.

Therefore, to explore the function of lncRNAs in cells upon H3N2 infection, lncRNA expression profiles were analyzed in A549 cells infected with mock or H3N2 by high-throughput sequencing technology. Our data provide important insights into the interactions between host lncRNAs and H3N2 infection.

Influenza strain A/swine/Zhejiang/04/H3N2 was stored in our laboratory. A549 cells were cultured in 100-mm cell bottles, grown to approximately 80%-90% confluence, and infected with H3N2 at a multiplicity of infection (MOI) of 10. After 2 h of incubation, the culture medium was removed, and infected cells were maintained in fresh medium containing 2% fetal bovine serum. The viral titer was measured by indirect immunofluorescence assay (IFA) using monoclonal antibodies against IAV M2 protein (Feng et al. 2018) at 12, 18, 24, and 36 h post-infection (hpi) . The 50% tissue culture infectious dose (TCID 50 ) was determined by the Reed-Muench method (Xue et al. 2016) . Uninfected A549 cells were used as a control. H3N2infected and mock-infected cells were harvested at 24 hpi and cryopreserved at -80°C for subsequent total RNA extraction. All samples were analyzed in duplicate. All virus experiments were carried out in a Biosafety Level 2 laboratory.

The total RNA was extracted using TRIzol reagent (Vazyme Biotech Company, China) according to the manufacturer's instructions. The concentration of total RNA was determined using a Nanodrop instrument (Thermo Fisher Scientific). RNA quality was assessed by determining the A260/A280 ratio, with a value of 1.8-2.0 indicating high quality. Ribo-zero-magnetic-kit (Epicentre, USA) was used to remove ribosomal RNA from the samples. RNA libraries were prepared using TruSeq RNA LT Sample Prep Kit v2 (Illumina, San Diego, CA, USA). Library sequencing was performed on an Illumina Hiseq 3000 platform (Illumina) by ShangHai Genergy Biotech (Shanghai, China). The RNA-sequencing data were deposited in the Gene Expression Omnibus database (accession number GSE135482).

Adaptors, low-quality reads, and poly-N were removed using Trim Galore software (http://www.bioinformatics. babraham.ac.uk/projects/). Quality control analysis was performed on clean data using FastQC software (http:// www.bioinformatics.babraham.ac.uk/projects/fastqc/). The filtered reads were then mapped to the human reference genome (version: Homo_sapiens.GRCh38) using STAR software suite (Dobin et al. 2013) . The transcripts were assembled with the mapped reads using StringTie (Pertea et al. 2015) . The assembled transcripts from unknown intergenic transcripts (U), a transfrag falling entirely within a reference intron (I), exonic overlap with reference on the opposite strand (X), and an intron of the transfrag overlaps a reference intron on the opposite strand were identified as novel lncRNAs. The coding ability of lncRNAs was predicted using four tools, including a predictor of lncRNAs and mRNAs based on k-mer scheme (PLEK) (Li et al. 2014) , Coding-Non-Coding-Index (CNCI) , Coding Potential Calculator (CPC) (Kong et al. 2007) , and Coding Potential Assessment Tool (CPAT) . The differentially expressed lncRNAs in H3N2-infected and non-infected groups were analyzed using DESeq2 (Love et al. 2014 ).

Cis-and trans-analyses were implemented to predict the target genes of differentially expressed lncRNAs. The known protein-encoding genes located within a 100-kb window upstream or downstream of lncRNAs were identified as cis target genes. Co-expressed coding genes were classified as trans-regulated target genes of differentially expressed lncRNAs. GO enrichment analysis of lncRNA target genes was performed using the TopGO package (https://www.bioconductor.org/packages/release/bioc/ html/topGO.html). The associated pathways of cis-or trans-regulatory target genes of the dysregulated lncRNAs were predicted by Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The statistical criteria in the analysis were as follows: P \ 0.05 and false discovery rate (FDR) \ 0.05. GO enrichment analysis of dysregulated mRNAs was performed using the TopGO package. The associated pathways of the dysregulated mRNAs were predicted by KEGG database. The statistical criteria in the analysis were as follows: P \ 0.05 and FDR \ 0.05.

In order to validate the results of high-throughput sequencing, RT-qPCR analysis was used to identify differentially expressed lncRNAs. The glyceraldehyde 3-phosphate dehydrogenase (GAPDH) gene was used as an internal standard. The primer sequences are presented in Supplementary Table S1 . RT-qPCR was carried out using AceQ qPCR SYBR Green Master Mix (Vazyme, China) according to the manufacturer's instructions. The relative expression level of each lncRNA was calculated by the 2 -DDct method. All samples were analyzed in triplicate. Statistical analyses were performed using Student's t tests.

To determine the propagation kinetics of H3N2 in A549 cells, we measured the viral titers at various time points after infection. Cells were infected with H3N2 at an MOI of 10 and then monitored by IFA at 12, 18, 24, and 36 hpi. The viral titer gradually increased at 12, 18, and 24 hpi and reached a maximum (10 5.7 TCID 50 /mL) at 24 hpi (Supplementary Fig. S1A and S1B). Based on these results, we selected cells infected with an MOI of 10 after 24 hpi for transcriptome analysis.

High-throughput RNA-seq was performed to determine the expression levels of lncRNAs in A549 cells infected with H3N2 or uninfected (mock). More than 258 billion raw base reads were generated for each sample using an Illumina Hiseq platform. After removing adaptor and lowquality sequences, average size of each clean read was 135 nt, and the clean data Q30 was above 93.45% (Supplementary Table S2 ).

Based on various criteria, including the specific location in genome, multiple exons, length greater than 200 nt, and noncoding characteristics, transcripts were filtered by three steps to identify the annotated and novel lncRNAs. In total, 3031 transcripts were identified as novel lncRNAs (Fig. 1A) . Next, protein-coding or noncoding transcripts were classified using four tools, i.e., CPAT, PLEK, CNCI, and CPC ( Fig. 1B) . Additionally, according to the corresponding genomic locations of transcripts of known protein-coding genes, newly identified lncRNAs were categorized into four groups, i.e., intronic lncRNAs (82%), intergenic lncRNAs (12%), antisense lncRNAs (3%), and bidirectional lncRNAs (3%; Fig. 1C ).

Hierarchical clustering was performed to analyze the lncRNA expression profiles in H3N2-infected or noninfected cells. Obviously, expression levels of lncRNAs were significantly altered after H3N2 infection ( Fig. 2A ). In total, 6129 lncRNAs were differentially expressed, including 4963 upregulated lncRNAs and 1166 downregulated lncRNAs (fold change [FC] C 2, P B 0.05; Fig. 2B ) (Supplementary Table S3 ). Of the differentially regulated lncRNAs, 22 newly identified lncRNAs were altered, with log 2 (FC) values of more than 12, compared with the noninfected group. The lncRNA showing the greatest upregulation was MSTRG.18254.3, with a log 2 (FC) of more than 16; in contrast, the lncRNA showing the greatest The differentially expressed lncRNAs were widely distributed in all chromosomes, although the numbers varied in different chromosomes. Most altered lncRNAs were located on chromosome 1, whereas few altered lncRNAs were located on chromosome Y (Fig. 2C ).

We also detected mRNA expression levels of A549 cells after H3N2 infection. Hierarchical clustering showed that the mRNA expression profiles were significantly altered in the H3N2-infected group compared with that in the noninfected group (Fig. 3A) . In total, 50,031 mRNA transcripts were found to be differentially expressed in A549 cells infected with H3N2 (FC C 2, P B 0.05) (Supplementary Table S4 ). Of the differentially expressed mRNA transcripts, 20,907 were upregulated, and 29,124 were downregulated (Fig. 3B) . Moreover, 833 genes were upregulated with a FC of more than 100 after infection. In our study, the gene showing the greatest upregulation was RPL4, with a log 2 (FC) of 12.3, whereas the gene showing the greatest downregulation was EGFR, with a log 2 (FC) of -14.2 (Table 2) . Similar to the distribution pattern of lncRNAs, the differentially expressed genes in H3N2-infected A549 cells were unevenly distributed among chromosomes. Most altered mRNAs were found on chromosome 1, whereas few altered mRNAs were found on chromosome Y (Fig. 3C ).

Next, we systematically analyzed the basic features of the lncRNAs and compared them with protein-coding genes. The lengths of lncRNA transcripts were longer than those of mRNAs (Fig. 4A) . Additionally, the number of exons of lncRNAs was lower than that of mRNAs (Fig. 4B ).

To better understand the roles of differentially expressed lncRNAs in H3N2-infected cells, GO and KEGG pathway analyses were used to explore the functions of cis-and trans-target genes of H3N2-induced lncRNAs. The results showed that the target genes of these lncRNAs were significantly enriched in biological processes, such as cellular metabolism. The top 20 significant GO terms are listed in Fig. 5A , 5B. The target genes of these differentially expressed lncRNAs participated in various signaling pathways, such as the B cell receptor signaling pathway and autophagy. The top 20 statistically significant enriched KEGG pathways are shown in Fig. 5C, 5D . These findings suggested that the induced lncRNAs regulated cellular metabolic processes, immunity, and autophagy during H3N2 infection. Because many genes were enriched in the cellular metabolic pathway, we further selected 10 of the most dysregulated lncRNAs to generate the cis-or transregulatory networks (shown in Fig. 5E, 5F ).

Dysregulated mRNAs in H3N2-infected cells were used for GO enrichment and KEGG pathway analyses. As shown in Fig. 6A , the majority of the dysregulated mRNAs were significantly enriched in some biological processes, such as cellular metabolic process, organic cyclic compound metabolic process, and cellular macromolecule metabolic process, suggesting that H3N2 infection had a profound effect on cellular metabolism in A549 cells. Furthermore, KEGG functional analyses showed that dysregulated mRNAs were significantly enriched in some pathways (Fig. 6B) , including citrate cycle (tricarboxylic acid [TCA] cycle), DNA replication, biosynthesis of unsaturated fatty acids, and autophagy. Taken together, these results suggested that the differentially expressed mRNAs belonged to multiple pathways. 

RT-qPCR was performed to further detect differentially expressed lncRNAs in sequencing data. Among the 6129 differentially expressed lncRNAs, we validated seven upregulated and two downregulated lncRNAs by RT-qPCR (Fig. 7A ) and RT-PCR (Fig. 7B) . The results showed that changes in lncRNAs expression levels were consistent with RNA-seq data. Moreover, we also tested the reported Regulation of cellular macromolecule biosynthetic process Cellular metabolic process Primary metabolic process Macromolecule metabolic process Cellular macromolecule metabolic process Nitrogen compound metabolic process Cellular nitrogen compound metabolic process Organic cyclic compound metabolic process Heterocycle metabolic process Cellular aromatic compound metabolic process Nucleobase-containing compound metabolic process Gene expression Nucleic acid metabolic process Macromolecule biosynthetic process Cellular macromolecule biosynthetic process Cellular nitrogen compound biosynthetic process Regulation of nitrogen compound metabolic process RNA metabolic process Regulation of macromolecule biosynthetic process lncRNAs, and the results were consistent with the literature (Fig. 7B, 7C ).

Traditional studies of host cell responses to influenza infection have mainly focused on protein-coding genes. Recently, several studies have shown that lncRNAs regulate virus replication by mediating host gene expression. For example, the lncRNA NRAV negatively regulates the initial transcription of multiple critical interferonstimulated genes by affecting histone modification and significantly promotes IAV replication (Ouyang et al. 2014) . Additionally, the lncRNA NEAT1, which is induced by influenza virus, is essential for the formation of nuclear body paraspeckles and thereby facilitates the expression of antiviral genes, including cytokines such as interleukin-8 (IL-8), probably through relocating transcriptional regulator splicing factor proline/glutamine-rich (a NEAT1-binding paraspeckle protein) from the IL-8 promoter to the paraspeckles (Imamura et al. 2014 ). Previous studies have also shown that the lncRNA Bst2/BISPR, which is induced by an influenza virus lacking NS1, is unable to block the interferon response. This lncRNA is located in genome close to BST and positively regulates the expression of BST2, but does not affect the expression of other genes adjacent to BST2 (Barriocanal et al. 2014) . VIN, a large intergenic lncRNA induced by several IAV strains, including H1N1, H3N2, and H7N7, promotes influenza virus replication (Winterling et al. 2014) . In this study, we found that 6129 lncRNAs and 50,031 mRNA transcripts were differentially expressed. Because several lncRNAs have been identified as modulators in virus-induced host responses, we predicted that these lncRNAs may play crucial roles in regulating the host response to IAV, although the mechanisms require further exploration.

Notably, in this study, the functions of lncRNAs were predicted according to their cis-or trans-target genes. GO terms were significantly enriched in biological processes, such as RNA processing, cell cycle, and cellular metabolic process, suggesting that lncRNAs induced by H3N2 infection may regulate metabolism and affect H3N2 replication. To date, several studies have suggested that viruses regulate host metabolism to affect virus replication. Influenza virus manipulates sphingolipid metabolism and increases the level of sphingosine kinase 1 to promote virus replication (Vijayan and Hahm 2014) . Human cytomegalovirus (HCMV) requires glycolysis for replication of the viral genome (McArdle et al. 2011) , and rhinovirus infection results in extensive alterations in cellular metabolism essential for viral replication, including enhancing utilization and uptake of glucose as well as enhancing nucleotide synthesis and lipogenesis (Gualdoni et al. 2018) . HCMV and herpes simplex virus-1 (HSV-1) both induce major metabolic alterations in their host cells, and HCMV profoundly enhances TCA compound levels and alters pathways that generate substrates for lipid metabolism essential for its replication. Moreover, HSV-1 upregulates pyrimidine nucleotide biosynthesis for its replication (Vastag et al. 2011) . LncRNAs regulate metabolic enzymes to affect virus replication. The lncRNA ACOD1 is induced in cells infected with various viruses and can bind the metabolic enzyme glutamic-oxaloacetic transaminase, increased its catalytic activity, and facilitate the production of metabolites that promote viral propagation (Wang et al. 2017) .

In addition to metabolism, KEGG pathway analysis showed that target genes of differentially expressed lncRNAs were enriched in autophagy. IAV infection plays complex roles in regulating autophagy. Indeed, IAV induces the formation of autophagosomes in rapamycintreated Madin-Darby canine kidney cells (Tanida et al. 2008; Zhou et al. 2009 ). Additionally, several viral proteins are involved in the induction of autophagy. For example, M2 protein alone is able to induce the initial steps of autophagosome formation (Gannage et al. 2009; Zhou et al. 2009; Zhirnov and Klenk 2013) , viral HA protein can slightly activate autophagy, and NS1 induces autophagy indirectly by promoting the synthesis of HA and M2 proteins (Zhirnov and Klenk 2013) . Moreover, IAV infection inhibits the degradation of autophagosomes by blocking their fusion with lysosomes (Gannage et al. 2009 ). Therefore, IAV likely regulates autophagy through lncRNAs.

In summary, in this study, we analyzed the expression profiles of lncRNAs in A549 cells infected by H3N2 for the first time. In total, 6129 lncRNAs were differential expressed in H3N2-infected A549 cells compared with that in normal cells. These results indicated that lncRNAs played regulatory roles in metabolism, autophagy, and other related pathways following H3N2 infection. Further studies will be needed to determine the roles of these lncRNAs play in H3N2 infection, and that may lead to novel vaccines developments for preventing IAV.

Long noncoding RNA BST2/BISPR is induced by IFN and regulates the expression of the antiviral factor tetherin

Long noncoding RNA: novel links between gene expression and innate immunity

A long noncoding RNA mediates both activation and repression of immune response genes

Non-coding RNAs: meet thy masters

Influenza virus: a brief overview

STAR: ultrafast universal RNA-seq aligner

Monoclonal antibody against the universal M2 epitope of influenza A virus

Matrix protein 2 of influenza A virus blocks autophagosome fusion with lysosomes

Rhinovirus induces an anabolic reprogramming in host cell metabolism essential for viral replication

Influenza: lessons from past pandemics, warnings from current incidents

Long non-coding RNAs and enhancer RNAs regulate the lipopolysaccharide-induced inflammatory response in human monocytes

Long noncoding RNA NEAT1-dependent SFPQ relocation from promoter region to paraspeckle mediates IL8 expression upon immune stimuli

Molecular basis for the generation in pigs of influenza A viruses with pandemic potential

Self-recognition of an inducible host lncRNA by RIG-I feedback restricts innate immune response

CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine

The threat of pandemic influenza: are we ready?

PLEK: a tool for predicting long noncoding RNAs and messenger RNAs based on an improved k-mer scheme

Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2

Inhibition of calmodulindependent kinase blocks human cytomegalovirus-induced glycolytic activation and severely attenuates production of viral progeny

Long noncoding RNA #32 contributes to antiviral responses by controlling interferon-stimulated gene expression

NRAV, a long noncoding RNA, modulates antiviral responses through suppression of interferon-stimulated gene transcription

StringTie enables improved reconstruction of a transcriptome from RNA-seq reads

Long noncoding RNAs: regulators of viral infection and the interferon antiviral response

The impact of influenza epidemics on mortality: introducing a severity index

Pandemic versus epidemic influenza mortality: a pattern of changing age distribution

Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts

lncRNA-RNA interactions across the human transcriptome

LC3 and autophagy

Mortality associated with influenza and respiratory syncytial virus in the United States

Expression profile and function analysis of long non-coding RNAs in the infection of coxsackievirus B3

Divergent effects of human cytomegalovirus and herpes simplex virus-1 on cellular metabolism

Influenza viral manipulation of sphingolipid metabolism and signaling to modulate host defense system

CPAT: coding-potential assessment tool using an alignment-free logistic regression model

An interferon-independent lncRNA promotes viral replication by modulating cellular metabolism

Evolution and ecology of influenza A viruses

Evidence for a crucial role of a host noncoding RNA in influenza A virus replication

Mitochondrial proteomic analysis of human host cells infected with H3N2 swine influenza virus

Propagation and characterization of influenza virus stocks that lack high levels of defective viral genomes and hemagglutinin mutations

Influenza A virus proteins NS1 and hemagglutinin along with M2 are involved in stimulation of autophagy in infected cells

Autophagy is involved in influenza A virus replication

Conflict of interest All authors declare that they have no conflict of interest.Animal and Human Rights Statement This article does not contain any studies with human or animal subjects.

Author Contributions JZ and YZ designed the experiments and wrote the paper. YZ performed the majority of the experiments. TY and YD analyzed the data. YL prepared the sequencing samples. JL verified sequencing results and BH edited pictures.