key: cord-0923529-xqokfg3l
authors: Wang, Anqi; Sun, Lipei; Wang, Mingshu; Jia, Renyong; Zhu, Dekang; Liu, Mafeng; Sun, Kunfeng; Yang, Qiao; Wu, Ying; Chen, Xiaoyue; Cheng, Anchun; Chen, Shun
title: Identification of IFITM1 and IFITM3 in Goose: Gene Structure, Expression Patterns, and Immune Reponses against Tembusu Virus Infection
date: 2017-03-13
journal: Biomed Res Int
DOI: 10.1155/2017/5149062
sha: c96d197b59d78b7d5dba68fb9a4e5ef67bd066d1
doc_id: 923529
cord_uid: xqokfg3l

As interferon-stimulated genes (ISGs), interferon-inducible transmembrane proteins 1 and 3 (IFITM1 and IFITM3) can effectively inhibit the replication of multiple viruses. Here, goose IFITM1 and IFITM3 were cloned and identified for the first time. The two proteins share the same topological structure and several important sites critical for the antiviral functions in other species are conserved in the goose. Goose IFITM1 and IFITM3 are most closely related to their respective orthologs in ducks; these proteins exhibited high mRNA transcript levels in immune-related tissues, including the thymus, bursa of Fabricius, and Harderian gland, compared to other tissues. Moreover, goose IFITM1 was highly constitutively expressed in gastrointestinal tract tissues, while goose IFITM3 was expressed in respiratory organs. Furthermore, goose IFITM3 was activated in goose peripheral blood mononuclear cells (PBMCs) infected with Tembusu virus (TMUV) or treated with Toll-like receptors (TLRs) agonists, while only the R848 and Poly (I:C) agonists induced significant upregulation of goose IFITM1. Furthermore, goose IFITM1 and IFITM3 were upregulated in the sampled tissues, to some extent, after TMUV infection. Notably, significant upregulation of goose IFITM1 and IFITM3 was detected in the cecum and cecal tonsil, where TMUV was primarily distributed. These data provide new insights into the immune effectors in geese and promote our understanding of the role of IFITM1 and IFITM3 in the defense against TMUV.

Types I and II interferons (IFNs) are critical for establishing an antiviral state, which is mediated by downstream IFNstimulated genes (ISGs) [1] . Viruses have evolved diverse strategies to escape immune defenses [2] . However, virus evasion of the interferon-inducible transmembrane proteins (IFITMs) restriction is not apparent.

Recently, IFITMs have become a popular research topic with the discovery that the immune-related IFITMs (IFITM1, IFITM2, and IFITM3) inhibit the early replication of multiple viruses, including influenza A virus (IAV), dengue virus (DENV), West Nile virus (WNV), severe acute respiratory syndrome coronavirus (SARS CoV), hepatitis C virus (HCV), vesicular stomatitis virus (VSV), Ebola virus (EBOV), and human immunodeficiency virus type 1 (HIV-1) [3] [4] [5] [6] [7] [8] . IFITMs are a subfamily of a larger transmembrane protein family called Dispanins [9] , which are generally considered to be comprised of IFITM1, IFITM2, IFITM3, IFITM5, and IFITM10 [10] . IFITM4P, a pseudogene, is located on mouse chromosome as IFITM6 and IFITM7 [11] . Moreover, IFITM1, IFITM2, and IFITM3, as viral restriction factors, are known as immune-related IFITMs, the antiviral activity of which is conserved from prokaryotes to vertebrates. In addition, the immune-related IFITMs exhibit high constitutive expression in target cells and are strongly induced by types I and II IFNs [12] . IFITMs have common domains that consist of N-and C-termini, two transmembrane domains, and a conserved cytoplasmic domain and share common properties: (1) the proteins contain a CD225 domain composed of transmembrane domain 1 (TM1) and a cytoplasmic domain and (2) two exons encode transmembrane polypeptides [13] . Among these domains, the N-terminal region plays an important role in the correct cellular localization of IFITM3. Removing the N-terminal 21 amino acids of IFITM3 results in a loss of association with the endosomal compartments and relocation to the cell surface, thereby abrogating its antiviral function [14] . IFITM1 possesses a shorter N-terminal region and predominately localizes at the cell periphery and in early endosomes [4] . The different cellular localization of IFITM1 and IFITM3 likely underlies their diverse antiviral efficiency against various viruses [15] . IFITM1 exerts more resistance to viruses that fuse at the plasma membrane or early endosome, such as HIV-1 [8] . However, IFITM3 is more proficient than IFITM1 at preventing infection by the late endosomal-or lysosomal-entering viruses, including IAV and DENV [3] . At present, IFITMs have been widely identified in vertebrates, and homologs of IFITMs are even found in bacteria [16] . However, little attention has been paid to IFITMs in birds. Only a small number of bird IFITM sequences have been deposited in public databases, including Gallus gallus IFITM1-like (GenBank: XM_001233949), Anas platyrhynchos IFITM1 (GenBank: KF584226), Gallus gallus IFITM3-like (GenBank: XM_420925), Anas platyrhynchos IFITM3 (GenBank: KF584228), Serinus canaria IFITM3 (GenBank: XM_009102512), and Nipponia nippon IFITM3 (GenBank: XM_009463652). In addition, it was reported that chicken IFITM3 [17] and duck IFITM3 [18] could restrict influenza viruses, but information regarding goose IFITMs was unavailable.

Tembusu virus (TMUV) is a newly emerging pathogenic virus that, together with DENV, WNV, Japanese encephalitis virus (JEV), and yellow fever virus (YFV), belongs to the genus Flavivirus [19] . The TMUV genome is a single-stranded RNA with an open reading frame that encodes three structural proteins, namely, the core (C), membrane (prM), and envelope (E) proteins, and seven nonstructural proteins [20] . In China, the virus was initially isolated from ducks and can cause severe egg drop syndrome in laying ducks, resulting in huge economic losses. However, recent studies showed that chickens, geese, and house sparrows are also susceptible to this virus [21] [22] [23] . Most importantly, serum samples from workers in the duck industry contained TMUV antibodies, which may indicate potential zoonotic transmission [24] . In view of the significant antiviral activity of IFITMs in the IFN-mediated innate antiviral defense and the lack of studies focusing on avian IFITMs, goose IFITM1 and IFITM3 were cloned for the first time. In addition, the amino acid sequences and topological structures of these proteins were predicted. To understand the evolutionary relationship of goose IFITM1 and IFITM3, a phylogenetic tree was constructed. Moreover, the tissue distribution of goose IFITM1 and IFITM3 in healthy geese was identified. It has been demonstrated that IFITMs play an important role in blocking the entry of pathogenic flaviviruses such as DENV and WNV.

As TMUV is a novel member of Flavivirus, the interaction between TMUV and IFITM-mediated immune responses is unknown. Interestingly, we observed that goose IFITM1 and IFITM3, to different degrees, positively responded to TMUV infection both in vivo and in vitro. Collectively, our results promote the research of avian IFITMs and may provide the molecular foundation for further research regarding the contribution of goose IFITM1 and IFITM3 to the innate immune responses to TMUV infection.

2.1. Animals, Viruses, and Agonists. All geese were purchased from the waterfowl breeding center of Sichuan Agriculture University and provided with sufficient water, fodder, and vegetables prior to and during experiments. All experiments were approved by the Institutional Animal Care and Use Committee of Sichuan Agriculture University (Number XF2014-18), Sichuan, China.

TMUV was a kind gift of the Avian Diseases Research Center of Sichuan Agricultural University. The agonists R848 (InvivoGen, USA), Poly (I:C) (Sigma, USA), ODN2006 (InvivoGen, USA), and LPS (InvivoGen, USA) were used to mimic treatment with single-stranded RNA, double-stranded RNA, synthetic oligonucleotides, and lipopolysaccharides, respectively.

The animals were euthanized, and then the various tissues were rapidly removed and frozen in liquid nitrogen. Total RNA was extracted from tissues using TRIzol (Invitrogen, Carlsbad, CA, USA) following the manufacturer's protocol and then reverse transcribed using 5x All-in-One RT MasterMix (Abm, Canada), which could remove the genomic DNA. Finally, the synthetic cDNA was stored at −80 ∘ C until needed for cloning and quantitative real-time PCR (qRT-PCR).

Coding Sequences (CDs). The primers IFITM1-F, IFITM1-R, IFITM3-F, and IFITM3-R (all sequences of the optimal primers used in this study are listed in Table 1 ) were designed based on the conserved regions of their respective counterparts in chicken and duck to amplify the partial sequences of goose IFITM1 and IFITM3. Subsequently, the full-length cDNA of IFITM3 and the 5 -untranslated region (UTR) of IFITM1 were obtained by rapid amplification of cDNA ends (RACE). Briefly, the gene-specific primers (GSPs), including 5 GSP0, 5 GSP1, 5 GSP2, 3 GSP1, and 3 GSP2, were designed based on the obtained partial sequences. For 5 -RACE, a homopolymeric tail was added to the 5 -end of the synthetic cDNA using terminal deoxynucleotidyl transferase and dCTP (TaKaRa, Japan), and then the first-strand cDNA was synthesized using the primer 5 GSP0 and a HiScript 1st Strand cDNA Synthesis Kit (Vazyme, USA). The 5 -end of the target gene was obtained following two rounds of nested PCR using the primer 5 GSP1 with the Abridged Anchor Primer (AAP) and 5 GSP2 with the Abridged Universal Amplification Primer (AUAP). For 3 -RACE, the Adapter Primer (AP) was used to amplify the first-strand cDNA. The 3 -end of IFITM3 was also amplified by nested PCR with the primer pairs 3 GSP1/AP1 and 3 GSP2/AP2. The primer IFITM3-3 R, which was used for amplifying the 3 -UTR of IFITM3, was designed according to the predicted IFITM3 sequence from goose transcriptome data [25] . Finally, the complete CDs of IFITM1 and IFITM3 were amplified using Primer STAR Max DNA polymerase (TaKaRa, Japan) and were cloned into the pMD19-T vector (TaKaRa, Japan) for sequencing.

Healthy Geese. To evaluate the tissue expression profiles of IFITM1 and IFITM3, various tissues, including the brain, bursa of Fabricius, cecum, cecal tonsil, gizzard, heart, Harderian gland, kidney, liver, lung, muscle, pancreas, proventriculus, small intestine, skin, spleen, thymus, and trachea, were collected from two-week-old gosling and adult goose. The mRNA expression levels of IFITM1 and IFITM3 in the various tissues were detected by the Bio-Rad CFX96 Real-Time PCR Detection System (Bio-Rad, USA) and EvaGreen 2x qPCR MasterMix-No Dye Kit (Abm, Canada) according to the manufacturer's protocol. The relative expression of target genes was calculated using the Livak and Schmittgen 2 −ΔΔCT method [26] and was normalized to goose GAPDH.

PBMCs were separated from goose blood using a goose Peripheral Blood Lymphocyte Separation Kit (TBD Sciences, China) according to the manufacturer's protocol. Then, PBMCs were cultured and maintained in RPMI 1640 medium (Gibco, USA) at 37 ∘ C in 5% CO 2 overnight. Subsequently, the cells were treated with different Toll-like receptors (TLRs) agonists, including R848 (5 g/mL), Poly (I:C) (30 g/mL), ODN2006 (50 g/mL), and LPS (25 g/mL), or challenged with 50 L TMUV. Cells treated with PBS were used as an experimental control, and the treated cells were challenged with each agonist or virus in quadruplicate and incubated for 6 h. Finally, the cells from each treatment group were harvested for qRT-PCR analysis of the expression levels of goose IFITM1 and IFITM3.

IFITMs In Vivo. Ten three-day-old goslings were divided into two groups with five goslings in each group. One group of goslings was infected with 500 L of TMUV (6.3 × 10 −6 TCID 50 /mL), and the second group was inoculated with an equivalent amount of normal saline. Five days after infection (dpi), three goslings from each group were euthanized, and eight tissues (bursa of Fabricius, cecum, cecal tonsil, Harderian gland, lung, small intestine, spleen, and thymus) were collected immediately for RNA isolation and cDNA synthesis. The relative expression levels of IFITM1 and IFITM3 in each sample were assessed by qRT-PCR as described above.

Goslings. The full-length E gene from TMUV was amplified using Primer STAR Max DNA polymerase (TaKaRa, Japan) and subcloned into the pET-32a(+) vector. The temperature and primers were optimized with two negative controls (DEPC water) and two positive controls (GPV and H9N2 templates). Three separate dilution series (10 −1 to 10 −8 ) of recombinant plasmid were prepared for the establishment of standard curves. The viral copy numbers (log 10 ) were normalized per 1 g of total RNA [27] . The viral copy number in TMUV-infected gosling tissues was quantified by qRT-PCR.

The potential open reading frames (ORFs) were predicted by ORF finder (https://www .ncbi.nlm.nih.gov/gorf/gorf.html), and the corresponding amino acid sequences of cDNA were translated with DNA-MAN software. The molecular weight was calculated using ExPASy (http://www.expasy.org/tools/). The transmembrane helices prediction was performed by TMHMM 2.0 (http:// www.cbs.dtu.dk/services/TMHMM/). The nucleotide sequences of IFITM1 and IFITM3 from numerous species were compared by pairwise identity analysis using the Species Demarcation Tool [28] . The multiple alignment analysis of amino acid sequences was performed using ClustalW [29] . The phylogenetic evolutionary tree was constructed by MEGA6 using the neighbor-joining method [30] . The qRT-PCR data were analyzed by Bio-Rad CFX Manager Software and GraphPad Prism 5 with an unpaired two-tailed -test.

IFITM3. The full-length cDNA of goose IFITM1 (GenBank: KX594328) is 649 bp, with a 55-bp 5 -UTR, a 390-bp ORF, and a 204-bp 3 -UTR. The identified goose IFITM3 (GenBank: KX594327) contains a 426-bp CD and 19-bp UTRs ( Figure 1 ). The CDs of goose IFITM1 and IFITM3 encode 129 and 141 amino acids (aa) with molecular weights of 14.0 kDa and 15.33 kDa, respectively. The nucleotide sequences from fourteen species were compared by pairwise sequence alignment analysis ( Figure 2 ). The 2D color-coded matrix indicates that goose IFITM1 shares the highest identity with its counterpart in Anas platyrhynchos (GenBank: KF584228.1) (93.3%), followed by 70.6% shared identity with the Gallus gallus gene (GenBank: XM_420925), while it shared less than 60% identity with the sequences of primate genes. In addition, the nucleotide sequence of goose IFITM3 shared 94.5%, 83.4%, and 74.7% identity with the Anas platyrhynchos, Gallus gallus, and Serinus canaria counterparts, respectively. However, only 64.4% identity was shared between goose IFITM1 and IFITM3 (Tables S1 and S2 in Supplementary Material available online at https://doi.org/10.1155/2017/5149062). The pairwise identity scores of IFITM1 and IFITM3 are shown in Tables S1 and S2, respectively.

Proteins. Since these two molecules are transmembrane proteins, the topological structures were predicted by the TMHMM server. Two transmembrane helices were found in both proteins, 1 t a a c t g a g t g g a c c c a a g c a c c g g t g c t c c a a c t c c c c g c t g c t g c c a g g c t g c c t c c t t a a c c c c c c c t g c c c c a a a c c c c a t c a c c c t t c c g a c c c c c t c c g t c c a g c t g a g g g g g t g c a c c c g c c g a g g t g a c c c c g t t g g g a g t g a g c c t c t g g g g t g a g c t c a a g a c g a a t a a a g c a t t c c t c a c t c c c t t t c a a a a a a a a a a a a a a a a a a a Figure 2 : The 2D color-coded matrix of IFITM1 and IFITM3. The nucleotide sequences of IFITM1 (a) and IFITM3 (b) from various species were compared by pairwise alignment using the Species Demarcation Tool. The color spectrum scheme was assigned according to the pairwise identity scores. Goose IFITM1 and IFITM3 shared the highest identity with their respective counterparts in duck.

TM1 CIL TM2 CTD IFITM1  IFITM1  IFITM1  IFITM1  IFITM1  IFITM3  IFITM3  IFITM3  IFITM3  IFITM3   IFITM1  IFITM1  IFITM1  IFITM1  IFITM1  IFITM3  IFITM3  IFITM3  IFITM3  IFITM3   53  52  51  82  70  74  74 predicted to be oriented toward the extracellular space. Multiple amino acid sequence alignments show several conserved residues that are critical for the antiviral function and localization of IFITMs. A highly conserved YxxΦ motif appears in the NTD of IFITM3, which promotes the internalization of IFITM3. However, IFITM1, with a comparatively short NTD, lacks this YxxΦ motif. In addition, three other YxxΦ motifs are conserved in IFITM3 of chickens, ducks, and geese. Two cysteines (C60 and C61) in TM1 are conserved in IFITM3 and undergo palmitoylation [31] . In addition, three lysines (K72, K77 and K93) that are the sites of ubiquitination are also present in IFITM3 [32] (Figure 3 ).

To explore the relationship of IFITM1 and IFITM3, as well as the evolutionary dynamics of goose IFITM1 and IFITM3, a phylogenetic tree was constructed using the amino acid sequences of IFITM1 and IFITM3 from diverse species. The phylogenetic tree indicated that IFITM1 and IFITM3 of diverse species were separated into two clusters. The avian IFITM1 and IFITM3 were clustered into two groups, respectively. Predictably, goose IFITM1 and IFITM3 were closely related to their homologs in ducks (Figure 4) .

The mRNA expression levels of IFITM1 and IFITM3 in eighteen tissues collected from healthy geese were assessed by qRT-PCR. In gosling, goose IFITM1 was highly expressed in the cecum, moderately expressed in the lung, bursa of Fabricius, trachea and cecal tonsil, and minimally expressed in the rest of the tissues. However, goose IFITM3 was widely expressed in all collected tissues. High relative expression levels of goose IFITM3 were detected in the Harderian gland and lung compared to the other tissues ( Figure 5 ). Meanwhile, muscle tissue exhibited the lowest transcription level of goose IFITM3. In adult goose, comparatively high expression levels of goose IFITM1 were observed in the thymus and small intestine. In addition, the production of goose IFITM3 was easily detected in the thymus followed by the lung, skin, and bursa of Fabricius but was only faintly exhibited in the muscle tissue of adult goose ( Figure 5 ).

Challenged PBMCs. IFITMs have been reported to be broadspectrum antiviral molecules [15, 33] . Once host cells are infected, the invasive pathogens are recognized by pattern-recognition receptors (PRRs) [34] , which drive the

Pan troglodytes IFITM3

Capra hircus IFITM3

Xenopus laevis IFITM3

Gallus gallus IFITM3

Anser cygnoides IFITM3

Anser cygnoides IFITM1

Mus musculus IFITM1

Bos taurus IFITM1

Homo sapiens IFITM1 production of IFNs. Subsequently, the interaction of IFNs and IFN receptors triggers the activation of ISGs such as IFITMs [1] . Here, whether goose IFITM1 and IFITM3 were activated in response to exogenous stimulation was monitored by qRT-PCR. As expected, goose IFITM1 and IFITM3 actively responded to treatment with either TMUV or TLRs agonists. TMUV and all the agonists led to extremely significant upregulation of goose IFITM3 in comparison with the PBS group ( Figure 6(b) ). However, only poly(I:C) and R848 effectively caused significant upregulation of goose IFITM1 (Figure 6(a) ).

To investigate the impact of TMUV infection on the transcription levels of goose IFITM1 and IFITM3, the relative expression levels of immune-related and epithelial tissues from goslings treated with normal saline or infected with TMUV were examined. As shown in Figure 7 , goose IFITM1 and IFITM3 were universally upregulated in response to TMUV infection. Compared to the mock-infected group, goose IFITM1 in the cecum, cecal tonsil, Harderian gland, and spleen was significantly upregulated in the TMUVinfected group. In addition, goose IFITM1 in the thymus was extremely significantly upregulated after TMUV infection (Figure 7(a) ). Moreover, extremely significant upregulation of goose IFITM3 was observed in the cecum at 5 dpi with TMUV. Meanwhile, the transcription levels of goose IFITM3 were obviously upregulated in the cecal tonsils and lungs from TMUV-infected animals (Figure 7 (b)). Furthermore, relatively high viral copy numbers were observed in these corresponding tissues, although the viral copy numbers in lung tissues were lower than those of other tissues (Figure 8 ).

Since IFITMs, as cell-intrinsic restriction factors, were first identified to have a role in the restriction of influenza and flavivirus infection in 2009 [3] , the antiviral function of IFITMs has been widely studied in various species, including human [35, 36] , pig [37, 38] , and mouse [39, 40] . However, in avian species, IFITMs were identified in only a small number of birds. The goose IFITM1 and IFITM3 identified here will facilitate a better understanding of the IFITM family in avian species.

As always, the nomenclature of IFITM1 and IFITM3 in chickens and ducks was confusing and controversial. According to the gene synteny with mammals, the predicted chicken Figure 7 were analyzed for the corresponding viral copy numbers by an absolute quantification method in triplicate. The viral copy numbers (log 10 ) were normalized to 1 g of total RNA. IFITM3-like (GenBank: XM_420925) and chicken IFITM1like (GenBank: XM_001233949) were renamed chicken IFITM1 and IFITM3, respectively. In addition, the names of the previously identified duck IFITM3 (GenBank: KF584228) and duck IFITM1 (GenBank: KF584226) were reversed. The pairwise sequence alignment analysis of nucleotide sequences showed that goose IFITM1 shares a much higher identity with its counterparts in duck (GenBank: KF584228) (93.3%) and chicken (GenBank: XM_420925) (70.6%) than with duck IFITM3 (GenBank: KF584226) (63.5%) and chicken IFTM3 (GenBank: XM_001233949) (59.4%). A similar result was observed in the nucleotide sequences alignment of goose IFITM3 with chicken IFITM1 and IFITM3, as well as duck IFITM1 and IFITM3 (Table S2 ). The phylogenetic analysis indicated that the IFITM1 and IFITM3 clusters were divergent subgroups. Goose IFITM1 and the renamed chicken IFITM1 and duck IFITM1, together with the mammalian IFITM1, were classified into the same subgroup. Meanwhile, goose IFITM3 and the renamed duck IFITM3 were located in the same monophyletic group, which was closely related to the renamed chicken IFITM3. In summary, the results revealed that it was necessary to rename the IFITM1 and IFITM3 in chicken and duck and that goose IFITM1 and IFITM3 have been conserved during evolution. The minor IFITM3 allele (SNP rs12252-C), which is equivalent to the N-terminal 21-amino acid-truncated IFITM3, was strikingly enriched in patients who were severely infected with the 2009 pandemic H1N1 virus [35, 41] . In line with the above findings, VSV-G envelope proteinmediated virus entry could not be inhibited by the IFITM3 mutant (the N-terminal 21-amino acid deletion) [14] . These results highlight the important role of the N-terminal region of IFITM3 in antiviral functions [14] . In addition, the 20-YxxΦ-23 sorting signal in the N-terminal region is critical for the internalization of IFITM3 from the cell periphery to endosomal compartments where it acts to impede virus entry [42] . The N-terminal YxxΦ motif was observed in the N-terminal region of goose IFITM3. We predicted that the conserved motif in goose IFITM3 is also involved in the internalization and antiviral actions. Moreover, three other YxxΦ motifs that might be required for the endosomal localization and antiviral functions of duck IFITM3 [18] are also conserved in IFITM3 of chickens and geese. However, whether these YxxΦ motifs in goose IFITM3 are critical for correct cellular localization and inhibition of virus entry has yet to be confirmed.

Analysis of the tissue distribution of goose IFITM1 and IFITM3 in geese showed that both goose IFITM1 and IFITM3 were readily detected in some immune-related tissues, including the thymus, bursa of Fabricius, and Harderian gland. In addition, the expression of goose IFITM1 was restricted and primarily confined to the gastrointestinal tract tissues (cecum, small intestine, and gizzard), which was observed previously in both chickens [17] and ducks [18] . In contrast, goose IFITM3 was constitutively and ubiquitously expressed in various tissues. Notably, high expression levels of goose IFITM3 were observed in respiratory tract tissues (lung and trachea), the target tissues of infection with influenza A viruses, compared to the other tissues, which might contribute to the inhibition of influenza A virus replication [43] . It is notable that goose IFITM1 and IFITM3 both exhibit the lowest levels of expression in muscle tissues, similar to observations in chickens [17] .

In vitro, we explored the immunological characteristics of goose IFITM1 and IFITM3. The mRNA level of IFITM3 in PBMCs sharply increased with TMUV infection and with R848, Poly (I:C), ODN 2006, or LPS treatment, although goose IFITM1 was prominently expressed only in the R848-and Poly (I:C)-stimulated groups. R848, Poly (I:C), ODN 2006, and LPS are known as TLR7, TLR3, TLR21, and TLR4 agonists, respectively. All four agonists and the TMUV genome can be recognized by PRRs, which then activate downstream signal transduction pathways triggering the production of cytokines, such as IFNs, that induce the transcription of ISGs, including IFITMs [1, [44] [45] [46] .

Previous studies demonstrated that immune-related IFITMs were significantly upregulated in highly pathogenic avian influenza-infected ducks [18, 47] . However, the available information on whether IFITMs positively respond to TMUV infection is scant. Obviously, the expression levels of goose IFITM1 and IFITM3 displayed an upward trend in all the tested tissues, especially the cecum (Figure 7) , where TMUV was primarily distributed (Figure 8 ). High viral copy numbers were also observed in the intestinal tract tissues of TMUV-infected ducks [48] . In addition, goose IFITM1 was intensively transcribed in the TMUV-infected spleen and thymus (Figure 7) , where the significant upregulation of the IFNs was previously observed [48] . After TMUV infection, the virus was intensively located in the collected tissues, which reasonably explains why higher expression levels of goose IFITM1 and IFITM3 appeared in the TMUV-infected gosling tissues than those of the mock-infected group. In addition, the lowest viral copy numbers were seen in lung tissues, which might be attributed to the high expression level of IFITM3 in the TMUV-infected goose lungs. In light of the positive response of goose IFITM1 and IFITM3 to TMUV infection, we speculate that goose IFITM1 and IFITM3 might facilitate IFN-mediated defenses against TMUV.

Goose IFITM1 and IFITM3 were characterized for the first time. In addition, high mRNA levels of goose IFITM1 and IFITM3 were presented in some immune-related tissues. Meanwhile, goose IFITM1 and IFITM3 were primarily expressed in intestinal tract tissues and respiratory organs, respectively. Furthermore, goose IFITM1 and IFITM3 were activated in response to TMUV infection in vitro and in vivo. These results indicated that goose IFITM1 and IFITM3 might be potential antiviral effectors that restrict TMUV infection.

A diverse range of gene products are effectors of the type i interferon antiviral response

Inverse interference: how viruses fight the interferon system

The IFITM proteins mediate cellular resistance to influenza A H1N1 virus, west nile virus, and dengue virus

Distinct patterns of IFITM-mediated restriction of filoviruses, SARS coronavirus, and influenza A virus

IFITM1 is a tight junction protein that inhibits hepatitis C virus entry

Interferon-induced cell membrane proteins, IFITM3 and tetherin, inhibit vesicular stomatitis virus infection via distinct mechanisms

Interferon-Induced transmembrane protein-mediated inhibition of host cell entry of ebolaviruses

The IFITM proteins inhibit HIV-1 infection

The Dispanins: a novel gene family of ancient origin that contains 14 human members

Evolution of vertebrate interferon inducible transmembrane proteins

The Fragilis interferon-inducible gene family of transmembrane proteins is associated with germ cell specification in mice

IFITM-family proteins: the cell's first line of antiviral defense

The small interferoninduced transmembrane genes and proteins

The N-terminal region of IFITM3 modulates its antiviral activity by regulating IFITM3 cellular localization

IFITMs restrict the replication of multiple pathogenic viruses

IFITMs from mycobacteria confer resistance to influenza virus when expressed in human cells

Chicken interferoninducible transmembrane protein 3 restricts influenza viruses and lyssaviruses in vitro

Duck interferon-inducible transmembrane protein 3 mediates restriction of influenza viruses

Flaviviruses

Analysis of the complete genome of tembusu virus, a flavivirus isolated from ducks in China

Adapted tembusu-like virus in chickens and geese in China

Identification and molecular characterization of a novel flavivirus isolated from geese in China

Characterization of a Tembusu virus isolated from naturally infected house sparrows (Passer domesticus) in Northern China

Tembusu Virus in Human, China

Immune-related gene expression patterns in GPV-or H9N2-infected goose spleens

Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method

Tembusu virus infection in Cherry Valley ducks: the effect of age at infection

SDT: a virus classification tool based on pairwise sequence alignment and identity calculation

Multiple sequence alignment with the Clustal series of programs

MEGA6: molecular evolutionary genetics analysis version 6.0

Palmitoylome profiling reveals S-palmitoylation-dependent antiviral activity of IFITM3

S-palmitoylation and ubiquitination differentially regulate interferoninduced transmembrane protein 3 (IFITM3)-mediated resistance to influenza virus

The broad-spectrum antiviral functions of IFIT and IFITM proteins

Pathogen recognition and innate immunity

IFITM3 restricts the morbidity and mortality associated with influenza

The Cterminal sequence of IFITM1 regulates its anti-HIV-1 activity

Bat and pig IFN-induced transmembrane protein 3 restrict cell entry by influenza virus and lyssaviruses

Evolutionary characterization of pig interferoninducible transmembrane gene family and member expression dynamics in tracheobronchial lymph nodes of pigs infected with swine respiratory disease viruses

Expression profile and histological distribution of IFITM1 and IFITM3 during H9N2 avian influenza virus infection in BALB/c mice

Human, pig, and mouse interferon-induced transmembrane proteins partially restrict pseudotyped lentiviral vectors

Interferon-induced transmembrane protein-3 genetic variant rs12252-C is associated with severe influenza in Chinese individuals

Identification of an endocytic signal essential for the antiviral action of IFITM3

IFITM3 inhibits influenza a virus infection by preventing cytosolic entry

The history of Toll-like receptors-redefining innate immunity

Toll-like receptor downstream signaling

Interferon response by Toll-like receptor signaling

A comparative analysis of host responses to avian influenza infection in ducks and chickens highlights a role for the interferon-induced transmembrane proteins in viral resistance

The sequential tissue distribution of duck tembusu virus in adult ducks

The authors disclose that there are no financial and personal relationships with other people or organizations that could inappropriately influence this work.

Anqi Wang and Shun Chen contributed equally as co-first authors of this work.