key: cord-0031063-cfofst5d
authors: Sui, Nana; Zhang, Ruihua; Jiang, Yue; Yu, Honglei; Xu, Guige; Wang, Jingyu; Zhu, Yanli; Xie, Zhijing; Hu, Jiaqing; Jiang, Shijin
title: Long Noncoding RNA Expression Rofiles Elucidate the Potential Roles of lncRNA- XR_003496198 in Duck Hepatitis A Virus Type 1 Infection
date: 2022-04-22
journal: Front Cell Infect Microbiol
DOI: 10.3389/fcimb.2022.858537
sha: 05315da00963aa57e51870889ad0a6121e96e2c3
doc_id: 31063
cord_uid: cfofst5d

Duck hepatitis A virus type 1 (DHAV-1) is a highly lethal virus that severely affects the duck industry worldwide. Long noncoding RNAs (lncRNAs) exert crucial roles in pathogen attacks. Here, we conducted deep transcriptome analysis to investigate the dynamic changes of host lncRNAs profiles in DHAV-1-infected duck embryo fibroblasts. We identified 16,589 lncRNAs in total and characterized their genomic features. Moreover, 772 and 616 differentially expressed lncRNAs (DELs) were screened at 12 and 24 h post-infection. Additionally, we predicted the DELs’ cis- and trans-target genes and constructed lncRNA-target genes regulatory networks. Functional annotation analyses indicated that the putative target genes of DELs participated in diverse vital biological processed, including immune responses, cellular metabolism, and autophagy. For example, we confirmed the dysregulation of pattern recognition receptors (TLR3, RIG-I, MDA5, LGP2, cGAS), signal transducers (STAT1), transcription factors (IRF7), immune response mediators (IL6, IL10, TRIM25, TRIM35, TRIM60, IFITM1, IFITM3, IFITM5), and autophagy-related genes (ULK1, ULK2, EIF4EBP2) using RT-qPCR. Finally, we confirmed that one DHAV-1 induced lncRNA-XR_003496198 is likely to inhibit DHAV-1 replication in DEFs. Our study comprehensively analyzed the lncRNA profiles upon DHAV-1 infection and screened the target genes involved in the innate immune response and autophagy signaling pathway, thereby revealing the essential roles of duck lncRNAs and broadening our understanding of host-virus interactions.

Duck virus hepatitis (DVH) is a highly lethal and fast-spreading disease affecting ducklings, which is caused by the duck hepatitis A virus (DHAV) (Chalmers and Woolcock, 1984; Yang et al., 2012; Zhang et al., 2020a) . As the only Avihepatovirus genus in the Picornaviridae family, the DHAV genus is genetically classified into three types: DHAV-1, DHAV-2, and DHAV-3 (Tseng et al., 2007; Gao et al., 2012) . Among the causative agents of DVH, DHAV-1 is the most common and widely distributed and is responsible for acute hepatitis characterized by neurological symptoms and ecchymotic hemorrhages on liver surfaces (Song et al., 2014; Ou et al., 2017) . In detail, DHAV-1 transmission occurs in the digestive and respiratory tracts via breath and direct contact between ducklings with a morbidity rate of up to 90% (Ding and Zhang, 2007; Jin et al., 2008; Wen et al., 2018) . Currently, DHAV-1 has become one of the major concerns in the duck industry (Fu et al., 2008; Zhang et al., 2020b) .

High throughput sequencing allows us to accurately detect and identify both protein-coding and noncoding RNAs (ncRNAs) (Stark et al., 2019) . Advanced genome-and transcriptome-wide analysis indicated that only 1%-2% of the genome encodes proteins, while over 75% is actively transcribed into non-protein-coding RNAs (Carninci et al., 2005) . According to the sequence length cutoff, ncRNAs are classified into two subclasses: small ncRNA (sncRNA, < 200 nucleotides) and long ncRNAs (lncRNAs, > 200 nucleotides). These ncRNAs exert central modulatory roles in gene expression regulation, RNA editing, and epigenetic activities (Hombach and Kretz, 2016) . To date, transcriptome sequencing combined with bioinformatics analysis has identified many new mRNAs and ncRNAs in animals coping with multiple pathogen attacks (Hu et al., 2016; Agliano et al., 2019) , which could help uncover novel pathogenic mechanisms of host-pathogen interactions.

Based on their genome location and direction relative to nearby protein-coding genes, lncRNAs can be designated as five types: sense, antisense, intergenic, intronic, and bidirectional (Iyer et al., 2015) . LncRNAs frequently act as molecular signals, scaffolds, decoys, or competing endogenous RNAs to affect gene expression and perform regulatory roles Yao et al., 2019) . Accumulating evidence suggests that lncRNAs serve as versatile regulators during diverse biological and physiological processes, especially in hostpathogen interactions against various infectious agents (Agliano et al., 2019) . For example, lnc-THRIL is a nuclear lncRNA that forms complex with heterogeneous nuclear ribonucleoprotein L (hnRNPL) to trans-regulate tumor necrosis factor a expression by binding to its promoter (Li et al., 2014) . Long intergenic ncRNA (lincRNA)-COX2 is another trans-acting lncRNA that regulates both the activation and repression of the inflammatory gene cyclooxygenase 2 by interacting with hnRNP-A/B and hnRNP-A2/B1 (Carpenter et al., 2013) . Meanwhile, lnc-LUCAT1 inhibits type I interferon (IFN-I) production by interacting with STAT1 in the nucleus, thereby suppressing the antiviral immune response (Agarwal et al., 2020) . Recent reports have emphasized the importance of lncRNAs, particularly in modulating host-pathogen interactions (Agliano et al., 2019; Hao et al., 2021) . However, limited information on lncRNAs' expression and function in ducks during DHAV-1 attack.

In this research, we first analyzed the lncRNA expression profiles in DHAV-1-infected duck embryo fibroblasts (DEFs), explored the lncRNAs' genomic structure, and screened the differentially expressed lncRNAs (DELs) transcripts for target genes predicting. Additionally, DELs and trans-target genes were obtained to establish co-expression networks. Gene Ontology items (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that many target genes participated in immune-related and autophagy-related pathways. Finally, we demonstrated the duck lncRNA-XR_003496198 (lnc-XR_003496198) as an antiviral factor against DHAV-1. Thus, our findings could help uncover the detailed pathogenesis of the DHAV-1 attack.

Fresh DEFs were isolated from 10-day-old SPF duck embryos and cultured as described previously (Lan et al., 2019) . When cells reached 80%-90% confluence, the DEFs were stimulated with pathogenic DHAV-1 LY0801 strain (3 multiplicity of infection (MOI)). At 2 h post-infection (hpi), the inoculum was discarded and replaced with a fresh maintenance medium containing 2% chicken serum (CS, Solarbio, Beijing, China) and 1% penicillin and streptomycin. The cell lysates were harvested in triplicates at 12, 24, and 36 hpi to detect the viral load.

All animal experiments were conducted following the guidelines issued by the Animal Care and Use Committee of Shandong Agricultural University (Approval Number: # SDAUA-2018-045).

DEFs' total RNA was isolated following the standard protocol using the Trizol ® reagent (Ambion). RNA quantification and quality were evaluated using a Bioanalyzer 2100 system (Agilent Technologies, CA, USA). Stranded mRNA sample preparation and next-generation sequencing were conducted by Novogene (Beijing, PR China). A total of nine libraries (which contained three biological replicates) were constructed and labeled as N1, N2, N3, D1, D2, D3, H1, H2, and H3 (N: non-infected DEFs harvest at 0 hpi, D: infected DEFs harvested at 12 hpi, H: infected DEFs harvested at 24 hpi). Transcriptome sequencing was performed on an Illumina Hiseq 2500 platform and PE150 (paired-end) runs. Clean reads were obtained according to the manufacturer's instruction and were then mapped to the genome of Anas platyrhychos using HISAT2 . Finally, the aligned reads were transferred to StringTie (Pertea et al., 2016) for transcript assembly.

The resulting transcripts were merged using Cuffmerge (Ghosh and Chan, 2016) and then identified lncRNAs from the assembled transcripts. Only transcripts with length > 200 nt and exon > 2 were retained, and the transcripts were then matched with a reference annotation using Cuffcompare (Ghosh and Chan, 2016) . The Coding potential was predicted using CNCI (Sun et al., 2013) , CPC (Kong et al., 2007) , and Pfam-scan (Finn et al., 2014; El-Gebali et al., 2019) . The transcripts that failed to encode any proteins were considered as candidate lncRNAs. Moreover, class code clustered candidates lncRNAs into lincRNA, intronic lncRNA, and antisense lncRNA. The transcripts' genomic features (length, exon length, and exon number) between different lncRNA types were compared.

The DESeq R package analyzed DELs of the control and DHAV-1 infection groups. The mapped outputs were counted using StringTie software, and expression values were normalized using Reads Per kilobase of transcript per Million mapped reads (RPKM) (Pertea et al., 2016) . DELs with | log2 (Fold Change) | > 0 & padj < 0.05 were filtered for further target genes prediction.

The DELs' target genes were predicted by applying two algorithms: the cis-(co-location) and trans-acting (coexpression) modes. The protein-coding genes near the position of lncRNA within 10/100kb were considered potential cis-targets.

Functional annotation enrichment analyses (GO and KEGG) for the cis-and trans-target genes were implemented using the cluster Profiler R package. The Go terms and KEGG pathways were considered significant when the corrected P-value was < 0.05.

Total RNA was isolated by the method mentioned above, cDNA was then generated using a reverse transcription kit (UE, China), and RT-qPCR was performed on a Light Cycler 480II instrument (Roche, Basel, Switzerland) using SYBER Green Real-Time PCR Kit (TaKaRa, Dalian, China) following the manufacturer's protocols.

Primers were either synthesized based on the reference sequences using NCBI primer BLAST or from published literature ( Table 1) . The DHAV-1 virus copies were identified following the method established by our laboratory (Lin et al., 2016) . The data were normalized via the comparative 2 -DDCT method using endogenous b-actin (DCt) as a control (Livak and Schmittgen, 2001) .

The selected lncRNAs were amplified from DEFs by traditional PCR, and the PCR products were subcloned into the PCAGGSflag plasmid with a homologous recombination kit (New Cell & Molecular Biotech Co, Ltd, China) and sequenced by Sangon (Shanghai, China). The small interference RNA (siRNA) targeting lnc-XR_003496198 and negative control siRNA were obtained from Genepharma (shanghai, China). Knockdown efficiency of the target lncRNAs was verified by RT-qPCR analysis.

DEFs seeded in six-well plates were transfected with the abovementioned plasmids, then stimulated with DHAV-1 for 24 h. The cell lysates were used for immunoblotting, and the supernatant was prepared for viral titers detection.

DEFs grown in 96-well plates were used for viral titration in triplicate with 10-fold serial dilutions. The virus titers were 

Data were presented as means ± SD. Significant differences were performed using Student's t-tests. Data analyses were evaluated using SAS version 9.0 (Cary, NC) and GraphPad Prism Software version 8.0 (San Diego, USA). * indicated P < 0.05 and ** indicated P < 0.01.

FCS (fetal calf serum) nonspecifically inhibits DHAV-1, whereas CS exerted no or minimal inhibitory influence (Chalmers and Woolcock, 1984) . Therefore, according to a previous report, we tried to culture DHAV-1 in a maintenance medium with 2% CS . To identify DHAV-1's replication efficiency in DEFs, we established the growth kinetics curve of the virulent LY0801 strain. The infection was monitored by observing cytopathic effects and virus copies. Compared to mock-infected cells without no cytopathic signs, the cytopathic effect (CPE) in infected DEFs started within 12 hpi, became more visible at 24 hpi, and the cells were crumbling and falling at 36 hpi ( Figure 1A ). DHAV-1 copies gradually increased before 36 hpi and decreased at 48 hpi ( Figure 1B ). Due to the severe shedding of infected DEFs at 36 hpi, we collected the DEFs at 12 and 24 hpi for transcriptome sequencing.

The Illumina Hiseq Xten platform was applied to profile lncRNA expression during DHAV-1 infection. More than 78 million raw reads were produced from each library (Accession number PRJNA715039). After initial processing, the clean reads comprised > 94% of the raw data, indicating the sequencing data's high accuracy. Subsequently, the clean reads mapped to the reference genome were highlighted, demonstrating a high fraction of reading alignment accuracy (average alignment rate ≥ 86.49%) ( Table 2) . Moreover, Pearson's correlation coefficients of all pairwise comparisons were > 0.90, attesting to the high consistency of the data ( Figure 2 ). Taken together, the sequencing data and biological replicates indicated sufficient data quality for further analysis.

To analyze the features of lncRNAs, we investigated the characteristics and transcription patterns based on the lncRNA profiles. A total of 11941 annotated lncRNAs and 4648 novel lncRNAs were identified. All detected lncRNAs were classified into four types, including 8840 intergenic, 1216 sense-intron, 1858 antisense, and 27 sense-overlapping lncRNAs ( Figure 3A ). Nearly half of duck lncRNAs were medium-length lncRNA (1000 -5000 bp: 49.8%), 31.4% were short lncRNA (200 -1000 bp), and only 18.1% were long lncRNAs (> 5 Kb) ( Figure 3B) . Moreover, the exons revealed that most lncRNAs had fewer than 5 exons ( Figure 3C ). The mapping of the lncRNAs to the duck chromosomes showed that all lncRNAs were distributed in each chromosome without a chromosome location preference ( Figure 3D ).

To identify the high-confidence DELs, the cutoff P < 0.05 and | log 2 FC | > 0 were used as a strict screening criterion. The D-vs-N comparison revealed 436 up-regulated and 336 down-regulated DELs ( Figure 4A , Table S1 ), while the H-vs-N comparison revealed 290 up-regulated and 326 down-regulated DELs ( Figure 4B and Table S1 ). Notably, the two lncRNA sets shared 268 DELs ( Figure 4C ). The clustering analysis of both DELs showed each group clustered together, demonstrating the low lncRNA profile variability within the same treatment ( Figure 4D ).

The target genes of DELs were predicted using cis and transregulatory relationships to explore interactive duck lncRNA-mRNA pairs. Regarding cis-action, 772 DELs (D-vs-N comparison) corresponded to 2404 protein-coding genes (Table  S2) , whereas 616 DELs (H-vs-N comparison) neighbored 1970 protein-coding regions (Table S2) . We also identified 2866 potential target genes for 772 trans-acting lncRNAs (D-vs-N comparison) ( Table S2 ) and 1807 potential target genes for 616 trans-acting lncRNAs (H-vs-N comparison) ( Table S2) . 

To further characterize the biological functions of the target genes, GO and KEGG analyses were performed. Only the significantly enriched GO terms (P < 0.05) were listed derived from cis-and trans-targets (D-vs-N, H-vs-N) ( Table S3) . The GO analysis showed that DELs' cis-target genes (D-vs-N comparison) were mostly engaged in stimulus-response, signal transduction, phosphorylation, and immune system processes ( Figure 5A ), whereas those of the H-vs-N comparison was mainly engaged in metabolic processes ( Figure 5B) . Moreover, DELs' trans-target genes (D-vs-N comparison) were mostly involved in biological, cellular, metabolic, single organism processes ( Figure 5C ).

Meanwhile, those of the H-vs-N comparison was mainly involved in metabolic processes ( Figure 5D ). KEGG enrichment analysis predicted the enriched pathways based on target genes ( Table S4) . The four groups' top 20 significantly enriched pathways were presented based on a Qvalue < 0.05 ( Figure 6) . Furthermore, we noted that the target genes were closely related to the immune system and signal transduction, suggesting that these pathways play an essential role in DHAV-1 progression.

Multiple trans-acting lncRNAs can simultaneously target some genes. We constructed an interaction network diagram to elucidate their co-expression patterns more clearly (Figure 7) . The results showed that the top six target genes (which have the most connections, including KCNV2, BORCS8, LOC106019909, LOC113841698, KCTD19, and COLEC10) co-expressed with 38, 35, 33, 33, 31, and 30 DELs, respectively, implying that they play essential roles in DHAV-1 infection.

We screened crucial immune-related signaling pathways and target genes upon DHAV-1 infection to investigate the virusinduced host innate immune response ( Table S5) (CCL20, IL8, TNFRSF8, IL10, IL6) . Furthermore, the number of target genes involved in these immune-related pathways was listed ( Table 3) . It was notable that most cis-target genes were related to the MAPK signaling pathway, whereas most trans-target genes were related to cytokine-cytokine receptor interaction. Notably, the immunerelated pathways had more cis-and trans-target genes enriched at 12 hpi (146/213) than at 24 hpi (100/148), indicating that DHAV-1 infection induced a stronger immune response at an early stage.

IFN-stimulated genes (ISGs) participate in antiviral defense. This study identified 100 and 57 potential ISGs among the cisand trans-target genes in the D-vs-N comparison, while 73 and 39 potential ISGs were found in the cis-and trans-target genes D-vs-H comparison, respectively ( Table 4 ). All these genes were previously identified in chickens and ducks (Dai et al., 2020; Yang et al., 2020a) . These results revealed that DHAV-1 infection could induce an intensive immune response and activate ISGs expression. Autophagy is a highly conserved self-eating process in eukaryotic cells. In our study, DELs' target genes were also involved in important autophagy pathways, including the mTOR signaling pathway (AKT1, HIF1A, RRAGD, IGF1) and the autophagy signaling pathway (ULK1, ULK2, IFNG).

The crucial immune-related and autophagy-related DELs were quantified using RT-qPCR to verify the relevance of the transcriptomic data. In general, the exhibited expression patterns of RT-qPCR were relatively matched the RNA-Seq data ( Figure 8A ) with a correlation coefficient of 0.86, confirming that the RNA-Seq data were reliable ( Figure 8B ). We also quantified the 15 immune-related and three autophagy-related target genes' expression ( Figure 8C ). In the infected group, the 15 immune-related genes were significantly up-regulated than those in the N group, whereas the autophagyrelated target genes ULK1, ULK2, and EIF4EBP2 exhibited low expression compared to the control group.

Based on the opinion that lncRNAs accumulated could somehow be involved in regulating virus replication, we screened and cloned five of those mentioned above strongly up-regulated lncRNAs plasmids and transfected them to DEFs to test if lncRNAs regulated DHAV-1. Firstly, we detected the transfection efficiency of the overexpression plasmids, and the results confirmed that the corresponding plasmids were significantly increased the lncRNA expression ( Figure 9A) . Subsequently, we transfected the above plasmids to DEFs, then followed infected with pathogenic DHAV-1. Surprisingly, among theses-lnc-XR_003496198, it was found to inhibit DHAV-1 infection in DEFs ( Figure 9B ). This transcript had no obvious annotation in the NCBI, but after checking its sequence using tools available with the CPC yielded a score of 0.40266 confirmed its noncoding characteristics. CNCI and PFAM also indicated it has almost no coding potential. We prepared two RNAi pairs and determined their efficiency on the knockdown to further evaluate the antiviral role of endogenous lnc-XR_003496198 during DHAV-1 infection lnc-XR_003496198 using RT-qPCR. The results confirmed that lnc-XR_003496198-RNAi #2 transfection suppressed lnc-XR_003496198 to 10% of control, while the RNAi #1 transfection knockdown lnc-XR_003496198 to 60% of the control ( Figure 9C ). As shown in Figures 9D, F , knockdown of lnc-XR_003496198 with siRNA significantly promoted the DHAV-1 virus titers. Meanwhile, immunoblotting demonstrated that the expression of DHAV-1 protein was inhibited in DEFs with overexpression of lnc-XR_003496198, whereas increased the DHAV-1 protein in DEFs with reduced expression of lnc-XR_003496198 compared to control cells ( Figure 9E) . Altogether, these results suggest that lnc-XR_003496198 elicited an inhibitory effect on DHAV-1 replication in DEFs.

DVH, mainly caused by DHAV, affects the duck industry by decreasing production efficiency (Tseng et al., 2007) . DHAV-1 causes severe acute liver injuries and is the most prevalent serotype in ducklings (Wen et al., 2018) . Thus, clarifying the DHAV-1-host interaction molecular mechanisms is urgent. LncRNAs exert vital roles in multiple biological processes such as metabolism and immune modulation (Hombach and Kretz, 2016) . Many studies have shown that lncRNAs are critical antiviral effectors during virus-host interactions (Zhu et al., 2020; Wang et al., 2021) . However, the role of lncRNAs in DHAV-1-exposed host cells remains unknown. Here, we profiled the expression patterns of lncRNAs in DHAV-1infected DEFs to elucidate the DHAV-1-induced host immune response. Among all the transcripts, a total of 16589 lncRNAs were found, and 4648 of them were novel (Table S1 ). In DHAV-1-infected DEFs, the most represented lncRNA category was lincRNA, and the majority of the lncRNA lengths ranged from 1000 to 5000 bp. The DEFs lncRNA characteristics were consistent with those of the previous investigation (Lin et al., 2020) .

Numerous pathogens can induce and regulate DELs upon infection (Lu et al., 2019; Ma et al., 2019; You et al., 2019) . Our transcriptome analysis revealed that DHAV-1 greatly dysregulated lncRNAs, which confirmed that these DELs might be potential DHAV-1 development regulators. Notably, more DELs were detected in the 12 hpi vs. control group (D-vs-N: 772 and H-vs-N: 616) , reflecting more intensive changes in the early stages of DHAV-1 infection.

The most important pattern for lncRNAs to function on neighboring genes is mainly in cis-or trans-manners (Bai et al., 2015; Yamada, 2017; Gao et al., 2020) . Consequently, we predicted the DELs target genes and performed GO and KEGG analysis to further explore the regulatory pathways potentially involved in DHAV-1 infection. According to functional annotation, infection enriched numerous innate immune response target genes. Therefore, most follow-up studies have focused on the immune-related target genes in response to DHAV-1, including the key pattern recognition receptors (PRRs), downstream pathways of PRRs, and effector molecules, such as cytokines and ISGs.

Innate PRRs play important defensive roles against viral infections (Zhao et al., 2016) . We found KEGG pathway annotations for various target genes involved in the Toll-like receptor (TLR) signaling, NOD-like receptor signaling, RIG-Ilike receptor signaling, and cytosolic DNA sensing pathways. For example, a panel of TLRs including TLR3, TLR4, and TLR7 (targeted by lnc-XR_001188764.3, lnc-TCONS_00161164; lnc-T C O N S _ 0 0 1 2 4 2 4 1 ; l n c -T C O N S _ 0 0 1 4 2 3 3 2 , l n c -TCONS_00002813, lnc-XR_002398570.2) was involved in the TLR signaling pathway. TLR3 is an essential PRR for monitoring viral double-stranded RNA (dsRNA), and lncRNAs are potential regulators in TLR3-driven antiviral responses (Murphy and Medvedev, 2016) . For example, lnc-NEAT1 facilitates TLR3stimulated CXCL8 expression in HeLa and A549 cells by binding to SFPQ (a transcriptional CXCL8 inhibitor), leading to SFPQ transposition away from the CXCL8 promotor (Boo and Yang, 2010) . Additionally, the up-regulation of TLR3 occurs during multiple virus infections such as Muscovy Duck Reovirus (MDRV), Duck Tembusu virus (DTMUV), and Zika virus (ZIKV) (Bendelja et al., 2010; Zhang et al., 2015; Dang et al., 2016) . Similarly, our results showed that DHAV-1 greatly increases TLR3 expression. The retinoic acid-inducible gene-I-like receptor (RLR) family consists of RIG-I, MDA5, and LGP2. RIG-I and MDA5 primarily recognize cytosolic ssRNA and dsRNA (Kato et al., 2006) , and they share similar protein structures, notably an N-terminal tandem caspase activation and recruitment domain (CARD), a central DExD/H box RNA helicase domain, and a C-terminal repressor domain (Loo and Gale, 2011). LGP2 lacks an Nterminal CARD, and its exact role is still controversial (Saito et al., 2007) . Upon ligand recognition, RIG-I and MDA5 recruit the MAVS adaptor protein via CARD-CARD interaction and activate the downstream TBK1 and inducible IKKs. Signal transduction activates NF-kB and IRF7, leading to their translocation to the nucleus and promoting IFN-I expression (Loo and Gale, 2011) . In this study, several components of the RLR signaling pathway, including DDX58 (RIG-I, targeted by six DELs), IFIH1 (MDA5, targeted by eight DELs), IRF7 (targeted by two DELs), TRIM25 (targeted by two DELs), and DHX58 (LGP2, targeted by five DELs) were all predicted based on the KEGG analysis and these genes expression were significantly upregulated from RT-qPCR analysis. Moreover, accumulating evidence suggests that lncRNAs play an important role in the RLR signaling pathway. For instance, lnc-ITPRIP-1 upregulation can positively regulate the RLR signaling pathway by targeting MDA5 and suppressing HCV replication (Xie et al., 2018) . Similarly, lnc-NEAT1 inhibits HIV infection by positively regulating RIG-I-DDX60-mediated IFN responses (Ma et al., 2017) . Moreover, lnc-zc3h7a and TRIM25 synergize to block virus attacks by enhancing TRIM25-mediated K63-linked ubiquitination of RIG-I, which is required for downstream antiviral signaling activation (Lin et al., 2019) . Therefore, we speculate that these DELs could target crucial RLR signaling pathway genes, thereby defending against DHAV-1 infection. Interestingly, DHAV-1 treatment up-regulated cGAS (targeted by lnc-XR_001193465.3) expression. cGAS plays a role in the cytosolic DNA receptor signaling pathway, and it is generally accepted that DNA viruses activate this pathway. Indeed, RNA viruses can also manipulate DNA sensors (Liu et al., 2017a; Abe et al., 2019) . It has been proven that DTMUVencoded NS2A protein disrupts du-STING-dependent antiviral cellular defenses by binding to du-STING . The potential interaction between the viral protein and cytosolic DNA sensing pathway facilitates the virus' escaping strategy. However, interactions between this pathway and DHAV-1 remain unclear and require in-depth research.

PRRs independently activated the JAK-STAT and the MAPK signaling pathways. The D-vs-N and H-vs-N groups contain 27 and 22 trans-target genes participating in the JAK-STAT signaling pathway, respectively, and 47 and 27 trans-target genes involved in the MAPK signaling pathway. STAT1 is a critical transcription factor involved in signal transduction (Mariéet al., 2018) targeted by the lnc-TCONS_00113417 in our study, and its expression was significantly up-regulated. However, many viruses associate with STAT1 to antagonize or inhibit the IFN-triggered JAK-STAT signal transduction cascade from evading the host's innate immune responses. For instance, TCBV NS4A inhibits STAT1/STAT2 phosphorylation and dimerization to block nuclear translocation and transcriptional regulation . Besides, DEV evades innate immunity clearance thanks to its UL47 protein, which interacts with STAT1 to suppress signal transduction (Dar et al., 2005) . Therefore, we speculated that DHAV-1 infection could activate the JAK-STAT signaling pathway efficiently.

Cytokines are soluble extracellular mediators such as chemokines, IFN, interleukins (IL), lymphokines, and tumor necrosis factors, which have important roles in immune response activation. IL-6 is a pleiotropic cytokine that regulates hematopoiesis and immune and inflammatory responses and contributes to T and B cell activation, proliferation, and differentiation (Chen et al., 2020) . IL-8, also known as CXCL8, is an important immune response mediator produced by macrophages and peripheral blood mononuclear cells (Zou et al., 2016) . In contrast, IL-10 can inhibit the synthesis and release of immune mediators and block antigen presentation to T cells and macrophages, thereby reducing the inflammatory response (Liu et al., 2013) . The KEGG analysis predicted that DHAV-1 would affect many important cytokines and cytokine receptor genes, including IL-6, IL-8, and IL-10. Besides, the lnc-TCONS_00113417, whose target gene was IL-8, exhibited the second highest expression level change (12.4-fold) among all DELs, and the expression of both IL-6 and IL-10 were upregulated, indicating that the protective response occurred during DHAV-1 infection. Compared with cytokine changes in DHAV-3-infected duck livers, IL-6 and IL-10 expression increased significantly (Zhang et al., 2018) . Therefore, we assumed that DHAV-1 and DHAV-3 induce similar changes in immunological characteristics.

Type "I" signaling pathway activation increases ISG expression, which plays an important role in antiviral defense, antiproliferative activities, and adaptive immunity stimulation (Sen and Sarkar, 2007) . In our study, we respectively identified 100 and 57 potential ISGs in the cis-and trans-target genes in the D-vs-N group and 73 and 39 in the D-vs-H group (Table S6) . These include interferon-induced transmembrane (IFITM) proteins, such as IFITM1, IFITM3, and IFITM5, which generally play a part in antiviral immunity. IFITM proteins restrict virus replication by blocking viral entry and preventing virus fusion with the late endosome or lysosome and penetration into the cytoplasm (Huang et al., 2011; Li et al., 2017) . Recently, many tripartite motif family proteins have been confirmed to be ISGs, which regulate the innate immune response. For example, TRIM25 has a role in TRAF6-mediated NF-kB activation and positively regulates the RLR signaling pathway (Lee et al., 2015) . Meanwhile, TRIM35 defends the host against the influenza A virus by forming circuits among TRIM35, TRAF3, and influenza A virus' PB2 . We confirmed the up-regulation of TRIM25, TRIM35, and TRIM66 (targeted by lnc-TCONS_00182498, lnc-XR_003500542.1, lnc-TCONS_00012417, lnc-TCONS_00191851; lnc-TCONS_00049845, lnc-XR_002401149.2; lnc-XR_003498594.1) using RT-qPCR. Hence, these genes may be essential to the response against DHAV-1 infection.

Besides the target genes participating in the innate immune response during virus attack, our transcriptome sequencing data also proved that some DELs targeted ULK1, ULK2, ATG3, and ATG12, which have been reported to be involved in autophagy (Levine and Kroemer, 2019) . Moreover, multiple DELs could simultaneously target genes, implying important roles in DHAV-1 infection. In our study, six target genes (KCNV2, LOC106019909, 113841698, BORCS8, COLE10, and KCTD19) were co-expressed with 38, 33, 33, 35, 30, and 31 DELs, respectively. Thus, the lncRNA-mRNA crosstalk provided more considerable information and confirmed the complexity of the gene regulatory network.

LncRNAs are the critical regulators between virus and host and are involved in the biological process ranging from control of viral and host genes, regulation of stability and translation of mRNAs, and manuscript of host antiviral response (Carpenter et al., 2013; Liu and Ding, 2017) . Our study found that lnc-XR_003496198, which has no obvious coding potential, is located at chromosome 3. Increased the expression of lnc-XR_003496198 can repress the replication of DHAV-1 at transcript and protein levels, while knockdown of it can remarkably promote DHAV-1 proliferation. However, the exact regulatory mechanisms of lnc-XR_003496198 have not been investigated and require further research.

In summary, we reported the first comprehensive description of lncRNA profiles upon DHAV-1 infection and screened multiple lncRNAs and putative target genes involved in innate immunity and autophagy. Moreover, seven lncRNAs were significantly induced upon DHAV-1 infection were determined, and the results showed that XR_003496198 overexpression in DEFs led to a reduced DHAV-1 transcription and translation. These findings provide valuable transcriptomic information for better understanding the immune response of host cells. Further investigations will focus on explaining specific regulation mechanisms of these lncRNAs and searching for a potential therapeutic target against Picornavirusrelated diseases.

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

The animal study was reviewed and approved by the Animal Care and Use Committee of Shandong Agricultural University (# SDAUA-2018-045).

Conceptualization, JH and NS. Methodology, JH. Software, NS and YJ. Validation, NS, YJ, and HY. Formal analysis, JW.

Investigation, GX and YZ. Resources, ZX. Data curation, NS. Writing-original draft preparation, NS and JH. Writingreview and editing, NS, JH and SJ. Visualization, NS, YZ. Supervision, NS, JH and SJ. All authors contributed to the article and approved the submitted version.

This study was supported by grants from the National Natural Science Foundation of China (31772754), Shandong Modern Agricultural Technology & Industry System, China , and Funds of Shandong. "Double Tops" Program, China (SYL2017YSTD11).

Cytosolic DNA-Sensing Immune Response and Viral Infection

The Long Noncoding RNA LUCAT1 is a Negative Feedback Regulator of Interferon Responses in Humans

Long Noncoding RNAs in Host-Pathogen Interactions

RNA Regulatory Networks in Animals and Plants: A Long Noncoding RNA Perspective

Decreased Toll-Like Receptor 8 Expression and Lower TNF-Alpha Synthesis in Infants With Acute RSV Infection

Intrinsic Cellular Defenses Against Virus Infection by Antiviral Type I Interferon

The Transcriptional Landscape of the Mammalian Genome

A Long Noncoding RNA Mediates Both Activation and Repression of Immune Response Genes

The Effect of Animal Sera on Duck Hepatitis Virus

Functional Characterization of Duck STING in IFN-b Induction and Anti-H9N2 Avian Influenza Viruses Infections

Gene Regulation in the Immune System by Long Noncoding RNAs

Silencing CDC25A Inhibits the Proliferation of Liver Cancer Cells by Downregulating IL −6 In Vitro and In Vivo

Duck Interferon Regulatory Factor 7 (IRF7) can Control Duck Tembusu Virus (DTMUV) Infection by Triggering Type I Interferon Production and its Signal Transduction Pathway

Systematic Identification of Chicken Type I, II and III Interferon-Stimulated Genes

Zika Virus Depletes Neural Progenitors in Human Cerebral Organoids Through Activation of the Innate Immune Receptor Tlr3

Transcriptional Analysis of Avian Embryonic Tissues Following Infection With Avian Infectious Bronchitis Virus

Molecular Analysis of Duck Hepatitis Virus Type 1

Pfam: The Protein Families Database

Molecular Detection and Typing of Duck Hepatitis A Virus Directly From Clinical Specimens

Genetic Variation of the VP1 Gene of the Virulent Duck Hepatitis A Virus Type 1 (DHAV-1) Isolates in Shandong Province of China

Comparative Transcriptome Analysis Uncovers Regulatory Roles of Long Noncoding RNAs Involved in Resistance to Powdery Mildew in Melon

Analysis of RNA-Seq Data Using TopHat and Cufflinks

Identification and Characterization of lncRNA AP000253 in Occult Hepatitis B Virus Infection

Distinct Patterns of IFITM-Mediated Restriction of Filoviruses, SARS Coronavirus, and Influenza A Virus

CpG Oligodeoxynucleotides Induce Differential Cytokine and Chemokine Gene Expression Profiles in Dapulian and Landrace Pigs

The Landscape of Long Noncoding RNAs in the Human Transcriptome

Identification and Molecular Analysis of the Highly Pathogenic Duck Hepatitis Virus Type 1 in Hubei Province of China

Differential Roles of MDA5 and RIG-I Helicases in the Recognition of RNA Viruses

HISAT: A Fast Spliced Aligner With Low Memory Requirements

CPC: Assess the Protein-Coding Potential of Transcripts Using Sequence Features and Support Vector Machine

Quantitative Proteomic Analysis Uncovers the Mediation of Endoplasmic Reticulum Stress-Induced Autophagy in DHAV-1-Infected DEF Cells

Regulation of MDA5-MAVS Antiviral Signaling Axis by TRIM25 Through TRAF6-Mediated NF-kb Activation

Biological Functions of Autophagy Genes: A Disease Perspective

The Long Noncoding RNA THRIL Regulates TNFalpha Expression Through its Interaction With hnRNPL

Circulation and In Vivo Distribution of Duck Hepatitis A Virus Types 1 and 3 in Infected Ducklings

The Long Noncoding RNA Lnczc3h7a Promotes a TRIM25-Mediated RIG-I Antiviral Innate Immune Response

Differently Expression Analysis and Function Prediction of Long Non-Coding RNAs in Duck Embryo Fibroblast Cells Infected by Duck Tembusu Virus

Coordinated Regulation of IFITM1, 2 and 3 Genes by an IFN-Responsive Enhancer Through Long-Range Chromatin Interactions

Roles of LncRNAs in Viral Infections

Higher Serum Levels of Interleukin 10 Occur at Onset of Acute Otitis Media Caused by Streptococcus Pneumoniae Compared to Haemophilus Influenzae and Moraxella Catarrhalis

RIGulation of STING Expression: At the Crossroads of Viral RNA and DNA Sensing Pathways

Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2-DDct Method

Immune Responses of Ducks Infected With Duck Tembusu Virus

Immune Signaling by RIG-I-Like Receptors

Identification and Analysis of Long Noncoding RNAs in Response to H5N1 Influenza Viruses in Duck (Anas Platyrhynchos)

The Long Noncoding RNA NEAT1 Exerts Antihantaviral Effects by Acting as Positive Feedback for RIG-I Signaling

HDAC Stimulates Gene Expression Through BRD4 Availability in Response to IFN and in Interferonopathies

Identification and Analysis of Long Noncoding RNAs That are Involved in Inflammatory Process in Response to Transmissible Gastroenteritis Virus Infection

Long Noncoding RNAs as Regulators of Toll-Like Receptor Signaling and Innate Immunity

Comparative Analysis of Virus-Host Interactions Caused by a Virulent and an Attenuated Duck Hepatitis A Virus Genotype 1

Transcript-Level Expression Analysis of RNA-Seq Experiments With HISAT, StringTie and Ballgown

Duck Interferon Regulatory Factor 1 Acts as a Positive Regulator in Duck Innate Antiviral Response

Regulation of Innate Antiviral Defenses Through a Shared Repressor Domain in RIG-I and LGP2

The Interferon-Stimulated Genes: Targets of Direct Signaling by Interferons, Double-Stranded RNA, and Viruses

Effect of Age on the Pathogenesis of DHV-1 in Pekin Ducks and on the Innate Immune Responses of Ducks to Infection

RNA Sequencing: The Teenage Years

TRIM35 Mediates Protection Against Influenza Infection by Activating TRAF3 and Degrading Viral PB2

Utilizing Sequence Intrinsic Composition to Classify Protein-Coding and Long Noncoding Transcripts

Molecular Analysis of Duck Hepatitis Virus Type 1 Indicates That it Should be Assigned to a New Genus

Fetal Calf Serum Exerts an Inhibitory Effect on Replication of Duck Hepatitis A Virus Genotype 1 in Duck Embryo Fibroblast Cells

Positive Feedback Loop of Long Noncoding RNA OASL-IT1 and Innate Immune Response Restricts the Replication of Zika Virus in Epithelial A549 Cells

Molecular Epidemiology of Duck Hepatitis a Virus Types 1 and 3 in Chin-2015

Long Noncoding RNA ITPRIP-1 Positively Regulates the Innate Immune Response Through Promotion of Oligomerization and Activation of MDA5

DHAV-1 Inhibits Type I Interferon Signaling to Assist Viral Adaption by Increasing the Expression of SOCS3

Functions of Long Intergenic Noncoding (Linc) RNAs in Plants

Development and Application of a Reverse Transcription Loop-Mediated Isothermal Amplification Method for Rapid Detection of Duck Hepatitis A Virus Type 1

Quantitative Proteomic Analysis of Duck Embryo Fibroblasts Infected With Novel Duck Reovirus

Tick-Borne Encephalitis Virus NS4A Ubiquitination Antagonizes Type I Interferon-Stimulated STAT1/2 Signalling Pathway

Cellular Functions of Long Noncoding RNAs

Integrated Analysis of lncRNA and mRNA Repertoires in Marek's Disease Infected Spleens Identifies Genes Relevant to Resistance

Comparative Liver Transcriptome Analysis in Ducklings Infected With Duck Hepatitis A Virus 3 (DHAV-3) at 12 and 48 Hours Post-Infection Through RNA-Seq

Binding of Duck Tembusu Virus Nonstructural Protein 2A to Duck STING Disrupts Induction of Its Signal Transduction Cascade To Inhibit Beta Interferon Induction

Molecular Cloning of Peking Duck Toll-Like Receptor 3 (Dutlr3) Gene and its Responses to Reovirus Infection

A Novel Peptide Isolated From a Phage Display Peptide Library Modeling Antigenic Epitope of DHAV-1 and DHAV-3

Evidence of Possible Vertical Transmission of Duck Hepatitis A Virus Type 1 in Ducks

Molecular and Genetic Inflammation Networks in Major Human Diseases

Long Noncoding RNAs are Associated With Seneca Valley Virus Infection

Identification of Downstream Target Genes Regulated by the Nitric Oxide-Soluble Guanylate Cyclase-Cyclic Guanosine Monophosphate Signal Pathway in Pulmonary Hypertension

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.