# Proteogenomics uncovers critical elements of host response in bovine soft palate epithelial cells following in vitro infection with foot-and-mouth disease Virus

Florian Pfaff, Sara Hagglund, Martina Zoli, Sandra Blaise-Boisseau, Eve Laloy, Susanne Koethe, Daniela Zuehlke, Katharina Riedel, Stéphan Zientara, Labib Bakkali, et al.

► To cite this version:

Florian Pfaff, Sara Hagglund, Martina Zoli, Sandra Blaise-Boisseau, Eve Laloy, et al. Proteogenomics uncovers critical elements of host response in bovine soft palate epithelial cells following in vitro infection with foot-and-mouth disease Virus. Viruses, 2019, 11 (1), pp.1-19. 10.3390/v11010053. hal-02628114

HAL Id: hal-02628114

https://hal.inrae.fr/hal-02628114

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

Article

# Proteogenomics Uncovers Critical Elements of Host Response in Bovine Soft Palate Epithelial Cells Following In Vitro Infection with Foot-And-Mouth Disease Virus

Florian Pfaff <fnref n="1" />,*, Sara Hägglund <fnref n="2" />, Martina Zoli <fnref n="1" />,*, Sandra Blaise-Boisseau <fnref n="3" />, Eve Laloy <fnref n="3" />,<fnref n="4" />, Susanne Koethe <fnref n="1" />, Daniela Zühlke <fnref n="5" />, Katharina Riedel <fnref n="5" />, Stephan Zientara <fnref n="3" />, Labib Bakkali-Kassimi <fnref n="3" />, Jean-François Valarcher <fnref n="2" />,*, Dirk Höper <fnref n="1" />,*, Martin Beer <fnref n="1" /> and Michael Eschbaumer <fnref n="1" />,*

Received: 20 December 2018; Accepted: 11 January 2019; Published: 12 January 2019

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Foot-and-mouth disease (FMD) is the most devastating disease of cloven-hoofed livestock, with a crippling economic burden in endemic areas and immense costs associated with outbreaks in free countries. Foot-and-mouth disease virus (FMDV), a picornavirus, will spread rapidly in naïve populations, reaching morbidity rates of up to 100% in cattle. Even after recovery, over 50% of cattle remain subclinically infected and infectious virus can be recovered from the nasopharynx. The pathogen and host factors that contribute to FMDV persistence are currently not understood. Using for the first time primary bovine soft palate multilayers in combination with proteogenomics, we analyzed the transcriptional responses during acute and persistent FMDV infection. During the acute phase viral RNA and protein was detectable in large quantities and in response hundreds of interferon-stimulated genes (ISG) were overexpressed, mediating antiviral activity and apoptosis. Although the number of pro-apoptotic ISGs and the extent of their regulation decreased during persistence, some ISGs with antiviral activity were still highly expressed at that stage. This indicates a long-lasting but ultimately ineffective stimulation of ISGs during FMDV persistence. Furthermore, downregulation of relevant genes suggests an interference with the extracellular matrix that may contribute to the skewed virus-host equilibrium in soft palate epithelial cells.

**Keywords:** foot-and-mouth disease virus (FMDV); bovine soft palate; nasopharynx; transcriptomics; proteomics; bioinformatics; virus-host interaction; innate immune system; interferon-stimulated genes (ISG)

<footnote n="1">¹ Institute of Diagnostic Virology, Friedrich-Loeffler-Institut, Federal Research Institute for Animal Health, 17493 Greifswald, Germany; Florian.Pfaff@fli.de (F.P.); martina.zoli2811@gmail.com (M.Z.); Susanne.Koethe@fli.de (S.K.); Dirk.Hoeper@fli.de (D.H.); Martin.Beer@fli.de (M.B.)</footnote>

<footnote n="2">² Swedish University of Agricultural Sciences, Host-pathogen interaction group, Division of Ruminant Medicine, 75007 Uppsala, Sweden; sara.hagglund@slu.se (S.H.); jean-francois.valarcher@slu.se (J.-F.V.)</footnote>

<footnote n="3">³ Laboratoire de Santé Animale de Maisons-Alfort, UMR 1161 virologie, INRA, Ecole Nationale Vétérinaire d'Alfort, ANSES, Université Paris-Est, 94700 Maisons-Alfort, France; sandra.blaise-boisseau@anses.fr (S.B.-B.); eve.laloy@vet-alfort.fr (E.L.); stephan.zientara@vet-alfort.fr (S.Z.); Labib.BAKKALI-KASSIMI@anses.fr (L.B.-K.)</footnote>

<footnote n="4">⁴ Biopôle EnvA, Ecole Nationale Vétérinaire d'Alfort, Université Paris-Est, 94700 Maisons-Alfort, France</footnote>

<footnote n="5">⁵ Institute of Microbiology, Department for Microbial Physiology and Molecular Biology, University of Greifswald, 17489 Greifswald, Germany; daniela.zuehlke@uni-greifswald.de (D.Z.); riedela@uni-greifswald.de (K.R.)</footnote>

<footnote>* Correspondence: Michael.Eschbaumer@fli.de; Tel.: +49-38351-7-1211</footnote>
# 1. Introduction

Foot-and-mouth disease (FMD) is an acute and severe systemic vesicular disease of cloven-hoofed animals (*Artiodactyla*) with tremendous economic impact. The last FMD epizootic in the European Union in the United Kingdom, Ireland, France and the Netherlands in 2001 culminated in the slaughter of more than 6.5 million animals and an economic toll of over €5 billion [1]. The etiological agent is foot-and-mouth disease virus (FMDV), the type species of the genus *Aphthovirus* in the family *Picornaviridae* [2]. FMDV particles comprise a non-enveloped icosahedral capsid that surrounds a single-stranded positive-sense RNA genome with an approximate length of 8.4 kilobases [3]. FMD mainly affects livestock such as cattle, buffalo, pigs, goats, and sheep, but can also be transmitted to deer and wild boar. It is endemic in wild buffalo in Southern Africa [4,5]. Although more than 70 species are known to be susceptible to FMDV, its primary host seem to be buffalo and cattle, in which the disease causes very high morbidity, but only low mortality in adults [6]. During the onset of acute infection, cattle are highly febrile and small vesicles develop on the mucosal membranes of the muzzle, lips, and oral cavity, as well as on the coronary band and interdigital space, and the teats of the udder.

Usually, cattle clinically recover within 2-3 weeks if no secondary infection occurs. At this time, many have completely cleared the virus, however, about 50% of animals may remain subclinically infected for up to three years, depending on the species [7]. The World Organisation for Animal Health (OIE) defines an animal from which infectious FMDV can be recovered by probang sampling later than 28 days post infection (dpi) as persistently infected or a so-called "carrier" [7]. The fear of contagion from carrier animals has severe consequences for trade in live animals and animal products [8]. In vivo studies have shown that in cattle more than 50% of animals will become persistently infected [9], even if they had been vaccinated against FMDV and did not develop clinical disease [10-12]. The exact anatomical sites of persistence are still debated, but different tissues of the upper respiratory tract including the nasopharynx have been suggested. A recent study of the tissue-specific localization of FMDV in persistently infected steers identified the contiguous epithelia of the dorsal soft palate (SP) and the dorsal nasopharynx as the most likely sites of persistence [10]. Evidence for FMDV persistence in lymph nodes and germinal centers was also put forward although no viral replication could be detected [13].

The cellular factors that promote establishment of FMDV persistence and the viral strategies of immune avoidance in the bovine nasopharynx remain currently poorly understood, making it impossible to predict which animals will develop into carriers and which will clear the virus. Previous in vitro work that aimed to decipher cellular responses during persistent FMDV infection hardly reflected persistence in vivo as the used models were either based on immortalized non-bovine cell lines, such as hamster kidney (BHK) cells [14], or bovine cell lines, e.g., bovine kidney (MDBK, EBK) cells [15,16], that are not from the primary site of persistence. O'Donnell et al. [17] used a persistently infected bovine pharynx cell line to examine the gene expression of selected cellular cytokines by RT-qPCR, showing differences in the expression of key antiviral cytokines between acute and persistent infection.

The rationale of this study was therefore to use a novel air-liquid interphase cell culture model based on primary SP cells from cattle [18] that closely resembles the situation in vivo, together with an proteogenomics approach that combined transcriptomics by high-throughput sequencing, RT-qPCR, proteomics and bioinformatics. Our results provide detailed insights into the transcriptional responses of the SP in reaction to acute FMDV infection and reveal long-lasting changes of gene activation throughout persistence. The identified pathways and genes may give rise to further investigations leading to early detection of persistence in cattle and novel vaccines that prevent the carrier state.

# 2. Materials and Methods

A schematic visualization of the sample processing and analysis workflow can be found in the Supplementary Figure S1.
## 2.1. *Ethics Statement*

The tissues used for the study were collected from animals slaughtered for food production. The animals were being processed as part of the normal work of the abattoir, therefore no ethics approval was required.

## 2.2. *Bovine Epithelial Cultures from Soft Palate*

Multilayers of bovine dorsal soft palate cells were propagated at the air-liquid interface for 5 weeks before FMDV infection, as described previously [18]. Briefly, bovine dorsal soft palate tissue, collected immediately after slaughter, was dissected and digested at 4 °C overnight in incubation medium supplemented with protease XIV (Sigma-Aldrich, St. Louis, MO, USA). Epithelial cells were thereafter scraped off the underlying tissue, filtered, and incubated in cell culture flasks for 4 h at 37 °C and 5% CO₂. Cells that did not adhere to the plastic were centrifuged at 200× g for 10 min at room temperature, frozen, thawed, and propagated for three to five passages in cell culture flasks before being seeded in 12 mm diameter Corning® Transwell-COL collagen-coated PTFE membrane inserts with 3.0 µm pores (Sigma-Aldrich). The cell culture medium was removed from the upper compartment after five days of culture and changed in the lower compartment every two or three days. The average number of cells that constituted the upper layer of the multilayer was estimated at 750,000.

## 2.3. *Experimental Design and FMDV Infection*

After 5 weeks of culture on inserts without passage, cells were infected with a twice-plaque-purified viral clone (FMDV O Clone 2.2, "Cl 2.2") derived from the O/FRA/1/2001 strain that was further propagated on BHK-21 cells (four passages) [16], or negative cell lysate, as described previously [18]. Two experiments (experiment 1 and 2) were performed with SP cells that originated from two different animals (a male and female, respectively) and that were infected at a multiplicity of infection (MOI) of 0.01 (compare Supplementary Table S1). Briefly, the inserts were incubated for one hour with 500 µL of clarified cell lysate from infected or uninfected cell cultures. Following infection and thereafter at a maximum interval of 3 days, the upper compartments were washed with 500 µL cell culture medium containing 10% FCS and, similarly, but only from 2 dpi, the medium in the lower compartments was changed. For each experiment, at days 0, 1 and 28, SP cells from 2 inserts were lysed with 750 µL of TRIzol Reagent (Life Technologies, Carlsbad, CA, USA) and frozen separately at -80 °C for transcriptomic and proteomic analyses.

## 2.4. *RNA and Protein Isolation*

For isolation of high-quality total RNA and proteins, 150 µL of trichlormethane (Carl Roth, Karlsruhe, Germany) was added to the lysed cells in 750 µL TRIzol Reagent and the mixture was centrifuged in order to separate the RNA-containing aqueous phase from the DNA- and protein-containing organic phase. The aqueous phase was then mixed with an equal amount of 100% ethanol (Carl Roth) and total RNA was extracted using the RNeasy Mini Kit (Qiagen, Hilden, Germany) with on-column DNase digestion with the RNase-Free DNase Set (Qiagen), following the manufacturer's instructions. The quantity and quality of total RNA was subsequently analyzed using a NanoDrop 1000 spectrophotometer (Peqlab, Erlangen, Germany) and RNA 6000 Pico chips on an Agilent 2100 Bioanalyzer (Agilent Technologies, Böblingen, Germany). All samples were checked for contamination using the 260/280 and 260/230 nm ratios, as well as for RNA degradation using the RNA Integrity Number (RIN).

Polyadenylated mRNA was subsequently isolated from 1–3 µg of high-quality total RNA using the Dynabeads mRNA DIRECT Micro kit (Invitrogen, Carlsbad, CA, USA) following the manufacturer's instructions. Prior to isolation, the ERCC ExFold RNA Spike-In mix 1 (Invitrogen) was supplemented and used as an internal control for all following steps. The quality of the mRNA and the extent
of ribosomal RNA contamination was assessed using RNA 6000 Pico chips on the Agilent 2100 Bioanalyzer (Agilent Technologies).

Proteins were extracted from the organic TRIzol phase following the manufacturer's instructions. Briefly, 0.3 mL of 100% ethanol (Carl Roth) were added and DNA was pelleted by centrifugation. To the supernatant, 1.5 mL of isopropanol (Carl Roth) were added and proteins were pelleted by centrifugation. Protein pellets were washed three times with 0.3 M guanidine hydrochloride (Carl Roth) in 95% ethanol (Carl Roth). After a final washing step using 95% ethanol, proteins were air-dried and the pellet was resuspended in freshly prepared 1% SDS (Carl Roth) by ultrasonication. Quality and quantity of the isolated proteins were checked using a 12% polyacrylamide gel (SDS-PAGE) and a colorimetric bicinchoninic acid (BCA) assay.

## 2.5. *Library Preparation and Sequencing*

For preparation of whole-transcriptome libraries the Ion Total RNA-Seq Kit v2 (Life Technologies) was used, following the manufacturer's instructions. Briefly, between 1 and 100 ng of the mRNA containing the ERCC spike-in control was treated with RNase III at 37 °C for 10 min. The fragmented mRNA was subsequently purified using the Magnetic Bead Cleanup Module (Life Technologies) and the resulting size distribution was assessed with the Agilent 2100 Bioanalyzer as described above. After hybridization and ligation of appropriate adapters, the fragmented mRNA was reverse transcribed into cDNA using SuperScript III enzyme. The cDNA was purified as described above and amplified for 14 cycles using Platinum PCR SuperMix High Fidelity along with appropriate primers for the generation of barcoded libraries. The resulting libraries were again purified using the method described above and the size distribution was assessed using the Agilent 2100 Bioanalyzer together with the DNA 7500 kit and chip (Agilent Technologies). All libraries were subsequently quantified using the KAPA Library Quantification Kit Ion Torrent (Kapa Biosystems, Wilmington, MA, USA) on a CFX96 Real-Time PCR Detection System (Bio-Rad Laboratories, München, Germany) and pooled at an equimolar ratio. For sequencing an Ion S5XL sequencing system (Life Technologies) along with the Ion 540 OT2 and Chip kit (Life Technologies) for the generation of up to 200 bp reads was used. Each library was sequenced in at least two independent sequencing runs.

## 2.6. *Statistical Analysis of Differential Gene Expression*

In order to detect problems and biases during mRNA isolation and library preparation, we used the ERCC_Analysis Plugin (version 5.8.0.1) provided in the Torrent Suite software (version 5.8.0). Only forward strand reads were selected for mapping and the minimum transcript count was set to 50. The raw reads from each sequencing library were quality checked using FastQC (version 0.11.7; Babraham Institute) with a focus on the read length distribution and adapter contamination. In order to quantify the expression of known bovine transcripts in these datasets we used Salmon (version 0.9.1) that uses a lightweight alignment method (quasi-mapping) for rapid transcript abundance estimation [19]. Briefly, we obtained the transcript reference of cattle (GCF000003055.6 Bos Taurus UMD 3.1.1) from NCBI and selected only transcripts that were featured as mRNAs, non-coding RNAs and miscellaneous RNAs (for details see Supplementary Table S2). A Salmon index was created using the option for a perfect hash, rather than a dense hash. Each sequencing library was then used as input for the major Salmon function "quant" using appropriate options for the library type (stranded single-end protocol with reads coming from the forward strand) and the mean read length. The number of bootstraps was set to 100 replicates. Subsequently a "tx2gene" table was prepared using the accessions of each RNA transcript from the aforementioned reference in the first column and the corresponding gene symbol in the second column (see also Supplementary Table S3). Using this table, the transcript abundance datasets from Salmon "quant" were imported into R workspace (version 3.4.1; [20]) using the "tximport" package (version 1.6.0; [21]) as implemented in the Bioconductor library. For handling and manipulating R scripts, the software RStudio (version 1.0.153) was used. Technical replicates of samples (multiple sequencing of the same library from a single sample) were
combined and the datasets were pre-filtered using only genes were more than four samples had raw gene counts greater than or equal to 100. In order to check the dataset for influence of treatment (infection, time and animal) and repeatability (replicates of same treatment) a principal components analysis (PCA) as well as a heatmap clustering was conducted using the regularized log transformed (rlog) read counts. DESeq2 (version 1.18.1; [22]) was then used to identify differentially expressed genes between the treatments based on the negative binomial distribution. The resulting *p*-values were adjusted with the Benjamini–Hochberg procedure and only genes with an adjusted *p*-value below 0.001 and an absolute fold change of >1 were considered significant. The logarithmic fold changes were further shrunken as recommended and described by Love et al. 2014 [22] in order to account for genes with low read counts. Significant differently expressed genes were annotated using the AnnotationDbi package (version 1.40.0) and used for further pathway enrichment analysis. Briefly, gene sets were analyzed using the “enrichPathway” function of the ReactomePA package (version 1.22.0; [23]) and the “enrichKEGG” function of the clusterProfiler package (version 3.6.0; [24]).

## 2.7. *Quantitative Reverse Transcription PCR (RT-qPCR)*

In order to confirm the results from the RNA sequencing experiment and to include samples from additional time points, a subset of six target genes (*ANKRD1*, *CASP7*, *IDO1*, *IFIH1*, *NCAM1*, *OAS2*) and two reference genes (*ACTB*, *GAPDH*) was selected for quantitative reverse transcription PCR (RT-qPCR) analysis. For each gene, two intron-spanning primer sets were designed using Primer3 (version 0.4.0; [25]) (for primer sequences see Supplementary Table S4). For RT-qPCR the QuantiTect Probe RT-PCR Kit (Qiagen) was used together with LightCycler 480 ResoLight Dye (Roche, Mannheim, Germany) according to the manufacturers' instructions. The following temperature profile was used on a Bio-Rad CFX96: 50 °C for 30 min, 95 °C for 15 min and 45 cycles of 94 °C for 15 s and 60 °C for 1 min. After each cycle, the fluorescence in the SYBR channel was detected and threshold cycle (Ct) values were deduced after each run. For each newly designed primer set the PCR efficiency was determined using an appropriate dilution series from an independent bovine RNA control. Using the ΔΔCt method, Ct values of target genes were first normalized by subtracting the Ct value of reference genes (ΔCt = Ct<sub>Target</sub> − Ct<sub>Reference</sub>). Mean and standard deviation of ΔCt were calculated from two technical and two biological replicates for each treatment group. Subsequently, the ΔΔCt was calculated by subtracting the ΔCT of the control group from the ΔCt of the treatment group at different time points (ΔΔCt = ΔCt<sub>Control</sub> − ΔCt<sub>Treatment</sub>) and used for calculation of the fold change of each gene (2<sup>−ΔΔCT</sup>).

## 2.8. *Protein Identification and Quantification, and Statistical Analysis of Differential Protein Expression*

Thirty µg of the extracted proteins were separated by SDS-PAGE, stained with Colloidal Coomassie Brilliant Blue G-250 and afterwards lanes were cut into ten equidistant pieces. In-gel digestion using trypsin and purification of tryptic peptides using Ziptips (C18, Millipore) prior to MS analysis were done as described previously [26]. LC-MS/MS were conducted using an EASY-nLC II coupled to a LTQ Orbitrap-Velos mass spectrometer (Thermo Fisher Scientific, Waltham, MA, USA). Peptides were separated at a constant flow rate of 300 nL/min using a binary 76 min gradient from 5% B to 75% B 99.9% ACN, 0.1% acetic acid). Survey scans in the Orbitrap were recorded with a resolution of 60.000 in a m/z range of 300–1700. The 20 most intense peaks per scan cycle were selected for CID fragmentation in the LTQ. Ions with unknown charge state, as well as singly charged ions were excluded from fragmentation. Dynamic exclusion of precursor ions for 30 s was enabled. Internal calibration (lock mass 445.120025) was enabled as well. For protein identification, resulting spectra were searched against a database containing sequences of *Bos taurus* including reverse sequences and common laboratory contaminants (44,376 entries). Database searches using Sorcerer-SEQUEST (version v.27, rev.11; Sage-N Research, Inc., Milpitas, CA, USA) and Scaffold (version v.4.8.4; Proteome Software, Portland, OR, USA) were done as described earlier [27].

In order to identify differentially expressed proteins we imported the raw NSAF values into R workspace and created an ExpressionSet using the Biobase package (version 2.38.0; [28]). We then
used the “Power Law Global Error Model” (PLGEM) as implemented in the “plgem” package (version 1.50.0; [29]) and first fitted it to the dataset using default settings and “FMDV infection” (samples from 24 h post infection (hpi) and 28 dpi) as fitting condition. After computation of observed and resampled signal-to-noise ratios, *p*-values for each detected protein were calculated.

## 2.9. *Data Availability*

The raw sequencing data along with deduced Salmon read count tables and substantial metadata are available at ArrayExpress (http://www.ebi.ac.uk/arrayexpress) under the accession number E-MTAB-7605. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://www.proteomexchange.org) via the PRIDE [30] partner repository with the dataset identifier PXD012242

# 3. Results

## 3.1. *A Cell Culture Model for FMDV Persistence*

The establishment and full characterization of the primary SP cell culture model are presented by Hägglund et al. [18]. In summary, the primary SP cells formed into multilayers and showed typical features of stratified squamous epithelia, such as tight junctions and impermeability to cell culture media. The cultures were inoculated either with a viral clone derived from the O/FRA/1/2001 strain or negative cell lysate for control. Infection with FMDV at low MOI resulted in limited cytopathic effect, high viral loads and presence of detectable viral antigen. After 28 dpi, the multilayers were still intact, but FMDV antigen and genome remained detectable at very low levels and viable virus could be isolated—giving a clear indication of FMDV persistence.

## 3.2. *RNA Sequencing and Exploratory Data Analysis*

The polyadenylated RNA fraction from 21 primary bovine SP cell culture samples was sequenced. From these, 10 and 11 samples originated from a female and male bovine, respectively. Cells for analysis were harvested immediately before inoculation (0 hpi), 24 hpi and 28 dpi (Figure 1A). For each time point (0 hpi, 24 hpi, 28 dpi) and treatment (FMDV, control), a minimum of four biological replicates (two from each animal) was sequenced. The number of reads for each sample ranged from 21.8 to 28.0 million, with an average of 24.7 million. In total, 518.4 million reads were included in the analysis (Supplementary Table S1).

These reads were assigned to a cattle transcript reference and the raw gene count data was transformed with respect to library size and transcript length using an appropriate model (Supplementary Figure S2). A principal components analysis (PCA) based on the normalized gene counts revealed that the samples from non-infected cultures remained closely together independent of the time in culture, while FMDV-infected samples showed a clear time-dependent grouping (Figure 1B). Accordingly, the first principal component (50% of variance) was assumed to represent the transcriptional differences between the samples caused by FMDV infection. Furthermore, transcriptional differences between both animals were clearly visible and represented by principal component 2 (29% of variance). The PCA showed no batch-to-batch effects and biological replicates from the same animal were grouped tightly.

Unsupervised cluster analysis of the 45 most variable gene transcripts clearly separated the FMDV-infected samples from the non-infected controls (Figure 1C). In detail, the infected samples showed a strong positive deviation from the per gene mean count (Figure 1C, green and blue transcript cluster). Interestingly, samples from the acute and persistent phases of infection were clearly separated, because the activation of the aforementioned genes was reduced at 28 dpi (Figure 1C, blue cluster). As in the PCA, the non-infected samples were divided by the differential expression of another set of transcripts (Figure 1C, red transcript cluster) and grouped according to the donor animal they
originated from. In summary, the explorative data analysis showed that the FMDV infection has a
clear effect on gene expression that is distinct from the animal-dependent effects.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Experimental setup of RNA sequencing and exploratory data analysis. (A) Primary soft palate (SP) cell cultures were obtained from two animals. Baseline samples for RNA sequencing were collected immediately before inoculation (circle). Subsequently, the cell cultures were inoculated with foot-and-mouth disease virus (FMDV) (orange symbols) or mock-infected for use as controls (blue symbols). Cells were then harvested for sequencing at 24 h post infection (hpi) (triangles) or 28 days post infection (dpi) (squares), representing acute and persistent infection, respectively. (B) Principal components analysis based on normalized gene counts of the 1000 most variable genes. (C) The variance of normalized gene counts was calculated for each gene and the 45 genes with the highest variance were selected and visualized in a heat map. The color of the cells indicates the difference from the mean normalized gene count of the corresponding gene. Samples and genes were clustered according to these differences (trees). The transcript cluster highlighted in green is similarly activated by acute and persistent FMDV infection, while the blue cluster highlights transcripts that are only active during the acute phase of infection. The red cluster comprises transcripts whose abundance differs between donor animals used in this experiment.</figcaption></figure>

3.3. *Differential Expression during Acute and Persistent FMDV Infection and associated Pathways*

In order to confirm the trends from the explorative data analysis and to address the transcriptional
changes during acute and persistent infection in more detail, a differential expression analysis
was conducted. The transcriptional response of cells from bovine SP tissue to FMDV infection
was assessed at 24 hpi and 28 dpi, in contrast to non-infected controls from the same time points
(Figure 2). A total of 312 and 73 gene transcripts were differentially expressed at 24 hpi and 28 dpi,
respectively. With 305/312 (97.7%) for 24 hpi and 65/73 (89%) for 28 dpi, the majority of differentially
expressed transcripts were up-regulated, while only 8/312 (2.5%) and 8/73 (11%) were down-regulated
(Figure 2A,B). From a total of 324 differentially expressed transcripts at both time points, 249 and 10
transcripts were solely regulated during either acute or persistent phase, respectively. The expression
of 63 transcripts was significantly up-regulated during both stages of infection (Figure 2C). The log2 fold change of these transcripts was much lower at 28 dpi when compared to 24 hpi (Figure 2D), with highest differences observed for *IFIT2*, *LOC101907799* and *CMPK2*. In contrast, the log2 fold change for *IFI27*, *PARP9*, *IFI6*, and *OAS1Z* was comparable or nearly equal during both phases of infection. *PLAC8* was slightly stronger regulated at 28 dpi than at 24 hpi. A full list of differentially expressed transcripts can be found in the supplementary material (Supplementary Table S5).

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Differential gene expression of SP cells during FMDV infection. (A) Volcano plot showing the log2 fold change (x-axis) and the adjusted *p*-value (y-axis) for all differentially expressed genes during the acute (24 hpi, green) and persistent (28 dpi, orange) phase of FMDV infection. The log2 fold change and adjusted *p*-value are calculated relative to a non-infected control from the same time point. The grey dotted lines indicate the cutoff values: adjusted *p*-value < 0.001 and |log2 fold change| > 1. (B) The histogram summarizes the number of genes that have a certain log2 fold change. (C) The total number of differentially expressed genes is visualized in a Venn diagram. The overlap indicates the 63 genes that are differentially expressed during both phases of infection. (D) The log2 fold change of these 63 genes is compared for acute and persistent infection.</figcaption></figure>

The significantly differently expressed genes during acute and persistent infection were further matched to specific metabolic and signaling pathways using the Reactome database [31] (Figure 3). During both infection stages, most of the regulated genes are associated with interferon signaling, in particular interferon α, β and γ signaling (Figure 3A and Table 1). The induction of these interferons
appears to be mediated by the DDX58/IFIH1 pathway. Accordingly, a number of interferon-stimulated genes (ISG), such as the ubiquitin-like family member ISG15, that induce strong antiviral mechanisms of the innate immune system are highly activated. During acute infection, transcripts involved in the interleukin-1 (IL1) and -10 signaling pathways as well as the programmed cell death (apoptosis) pathway were enriched. However, it is of note that the IL1 and IL10 genes were not significantly differently expressed. The MHC I antigen presentation pathway, a part of both the innate and adaptive immune system, was also enriched at 24 hpi. During persistence, most of the genes belonging to the interleukin-1, -10, apoptosis and MHC I antigen presentation pathways were not significantly differentially expressed in comparison to control samples.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Pathway enrichment analyses of significantly differentially expressed genes during acute and persistent FMDV infection of bovine soft palate cells. (A) The number of genes that contribute to a significantly enriched Reactome pathway is shown for 24 hpi (blue dots) and 28 dpi (orange dots). Black lines connect pathways that are enriched in both datasets and the size of the circles represents the ratio of genes found for a certain pathway to the overall gene number. (B) The log2 fold change of genes that contribute to selected enriched Reactome pathways are highlighted (24 hpi—blue dots) and 28 dpi—orange dots). The grey dotted lines indicate a |log2 fold change| > 1. The numbers on the right correspond to the number of genes associated with each pathway.</figcaption></figure>
<table><caption>**Table 1.** Selected enriched pathways and corresponding genes during acute and persistent FMDV infection.</caption><thead><tr><th>Metabolic Complex</th><th>Pathway</th><th>24 hpi</th><th>28 dpi</th></tr></thead><tbody><tr><td rowspan="2">Innate immune system</td><td>DDX58/IFIH1-mediated induction of interferon-alpha/beta</td><td>DDX58, DHX58, HERC5, IFIH1, IFNB1, IRF1, IRF3, IRF7, ISG15, NFKBIA, NFKBIB, NLRC5, TNFAIP3, TRIM25, UBA7</td><td>DDX58, DHX58, HERC5, IFIH1, IRF7, ISG15, UBA7</td></tr><tr><td>Cytokine signaling in immune system</td><td>Interferon alpha, beta signaling<br>ADAR, GBP2, IFI27, IFI35, IFI6, IFIT2, IFIT3, IFITM1, IFITM2, IFNB1, IRF1, IRF3, IRF5, IRF7, IRF9, ISG15, MX1, MX2, OAS2, PSMB8, RNASEL, RSAD2, SOCS1, STAT2, USP18, XAF1</td><td>GBP2, IFI27, IFI35, IFI6, IFIT2, IFITM1, IRF7, IRF9, ISG15, MX1, MX2, OAS2, RSAD2, USP18, XAF1</td></tr><tr><td rowspan="3">Adaptive immune system</td><td>Interferon gamma signaling</td><td>GBP2, IFI27, IFI35, IFI6, IFIT2, IFITM1, IRF7, IRF9, ISG15, MX1, MX2, OAS2, RSAD2, USP18, XAF1</td><td>GBP2, GBP4, GBP5, IRF7, IRF9, NCAM1, OAS2, SP100</td></tr><tr><td>Interleukin-1 family signaling</td><td>IL18, IL18BP, MAP3K8, NFKBIA, NFKBIB, PEL11, PSMB10, PSMB8, PSMB9, PSME2, PSMF1</td><td></td></tr><tr><td>Antiviral mechanism by IFN-stimulated genes/ISG15 antiviral mechanism</td><td>CCL2, CCL5, CXCL10, CXCL2, CXCL8, IL18, IL6</td><td>CCL5</td></tr><tr><td rowspan="2">Programmed cell death</td><td>Class I MHC mediated antigen processing and presentation</td><td>DDX58, EIF2AK2, HERC5, IRF3, ISG15, MX1, MX2, TRIM25, UBA7, USP18</td><td>DDX58, HERC5, ISG15, MX1, MX2, UBA7, USP18</td></tr><tr><td>Activation of NF-kappaB in B cells</td><td>AREL1, CTSS, DTX3L, ERAP2, HECTD2, HERC5, HERC6, PSMB10, PSMB8, PSMB9, PSME2, PSMF1, RBCK1, RNF114, RNF19B, SOCS1, TAP1, TRIM21, UBA7</td><td>DTX3L, HERC5, HERC6, UBA7</td></tr><tr><td></td><td>Programmed cell death</td><td>NFKBIA, NFKBIB, NFKBIE, PSMB10, PSMB8, PSMB9, PSME2, PSMF1</td><td></td></tr><tr><td></td><td></td><td>CASP7, CFLAR, PMAIP1, PSMB10, PSMB8, PSMB9, PSME2, PSMF1, RIPK3, TNFSF10</td><td>TNFSF10</td></tr></tbody></table>

All of the genes contributing to the aforementioned pathways were up-regulated during the acute and persistent phases of infection, with a single exception for the nuclear cell adhesion molecule 1 (NCAM1) that is related to interferon γ signaling and was down-regulated during persistence (Figure 3B, and Table 2). In addition to NCAM1, other proteins that are associated with extracellular membrane (ECM) pathways, such as type I and III collagens (COL1A1 and COL3A1) were also down-regulated during persistence.

<table><caption>**Table 2.** Genes that are differentially expressed only during the persistent phase of FMDV infection.</caption><thead><tr><th>Gene</th><th>Description</th><th>LFC<sup>†</sup></th><th>Adjusted<br>p-value</th><th>Enriched DAVID Terms</th></tr></thead><tbody><tr><td>NCAM1</td><td>neural cell adhesion molecule 1</td><td>-4.48</td><td>3.35 × 10<sup>-4</sup></td><td>signal peptide, secreted</td></tr><tr><td>ANKRD1</td><td>ankyrin repeat domain 1</td><td>-4.19</td><td>6.61 × 10<sup>-5</sup></td><td>positive regulation of apoptotic process</td></tr><tr><td>SFRP2</td><td>secreted frizzled related protein 2</td><td>-4.16</td><td>3.31 × 10<sup>-5</sup></td><td>signal peptide, secreted, positive<br>regulation of apoptotic process</td></tr><tr><td>COL1A1</td><td>collagen type I alpha 1 chain</td><td>-3.50</td><td>4.16 × 10<sup>-4</sup></td><td>signal peptide, secreted</td></tr><tr><td>COL3A1</td><td>collagen type III alpha 1 chain</td><td>-3.28</td><td>3.23 × 10<sup>-4</sup></td><td>signal peptide, secreted</td></tr><tr><td>MYLK</td><td>myosin light chain kinase</td><td>-3.26</td><td>9.61 × 10<sup>-7</sup></td><td></td></tr><tr><td>HTRA3</td><td>HtrA serine peptidase 3</td><td>-3.21</td><td>2.93 × 10<sup>-7</sup></td><td>signal peptide, secreted</td></tr><tr><td>ALDH1A2</td><td>aldehyde dehydrogenase 1 family member A2</td><td>-3.09</td><td>7.83 × 10<sup>-5</sup></td><td>positive regulation of apoptotic process</td></tr><tr><td>LY6E</td><td>lymphocyte antigen 6 family member E</td><td>3.45</td><td>2.59 × 10<sup>-6</sup></td><td>signal peptide</td></tr><tr><td>FBP1</td><td>fructose-bisphosphatase 1</td><td>4.09</td><td>2.94 × 10<sup>-5</sup></td><td></td></tr></tbody></table>

### 3.4. RT-qPCR and Quantitative Proteomics

Results from the differential gene expression analysis based on RNA sequencing were confirmed using RT-qPCR. A panel of two reference (ACTB, GAPDH) and six target genes (OAS2, IFIH1, NCAM1, ANKRD1, IDO1, CASP7) were selected for analysis. These genes were chosen according to their apparent regulation: OAS2, IFIH1 and IDO1 were up-regulated during acute and persistent infection, while NCAM1 and ANKRD1 were only down-regulated during persistence. CASP7 was only up-regulated during the acute phase of FMDV infection. The selected genes were also involved in relevant signaling pathways: IFIH1 is part of the interferon induction cascade, while OAS2 is an

<footnote>†: log2 fold change.</footnote>
interferon-induced gene with antiviral activity. *CASP7* plays a key role in apoptosis and *IDO1* encodes a metabolically active protein that supports the immune system. *ANKRD1* and *NCAM1* were chosen as interesting target genes that were only regulated during persistence. The RT-qPCR analysis included samples from additional time points (2, 24, and 48 hpi; 7 and 28 dpi) for which log2 fold changes were then calculated using the delta-delta Ct method. For each time point and cell type four biological replicates (two replicates per animal) were used for analysis with the exception for 48 hpi, from which only two replicates derived from a single animal were available. The log2 fold-changes obtained by RT-qPCR of SP cells at the time points 24 hpi and 28 dpi confirm the results of the RNA sequencing (Figure 4A). All of the analyzed genes showed a time-dependent regulation during infection with FMDV (Figure 4B). The genes *IDO1*, *OAS2*, and *IFIH1* were strongly up-regulated from 2 to 24 hpi and were still up-regulated after 28 dpi. In contrast, *CASP7* was only up-regulated during acute infection. *ANKRD1* and *NCAM1* were down regulated as early as 48 hpi and throughout 28 dpi.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Expression analysis for selected genes by RT-qPCR at different time points. (A) Gene expression in FMDV infected soft palate cells compared to uninfected controls, at 24 hpi (blue panel) and 28 dpi (orange panel). The change in gene expression was analyzed for six selected target genes and two reference genes (not shown) using RT-qPCR. Log2 fold changes were calculated in comparison to a non-infected control from the same time points. Grey and white bars indicate the log2 fold change of gene expression using results from RT-qPCR and RNA sequencing (RNASeq), respectively. The grey dotted lines indicate a |log2 fold change| > 1. (B) The gene expression of the six selected genes was traced in a time course 2, 24, and 48 hpi; and 7 and 28 dpi. Statistically significant changes in comparison to control samples are highlighted with asterisks (*p*-value < 0.05: *; *p*-value < 0.01: **). Time points 24 hpi and 28 dpi are illustrated in blue and orange, respectively.</figcaption></figure>

Additionally, a GeLC-MS/MS-based comparative proteome analysis was conducted to confirm the results of the RNA sequencing on the protein level, using two biological replicates for each time point and treatment from animal 2. Proteins with significantly different abundance in both experiments were identified and quantified by spectral-counting (Table 3). At the given thresholds, 11 and 8 proteins were differentially expressed when comparing FMDV-infected and control samples at 24 hpi and
28 dpi, respectively. Of these, 6 were differentially expressed at both time points. The Normalized Spectral Abundance Factor (NSAF), measured as the proportion of a single protein in relation to the detected proteins overall revealed an accumulation of the proteins *HERC6*, *IFI44*, *IFI44L*, *ISG15*, *MX1*, *MX2*, *OAS1X*, and *OAS1Y* in cells during persistent infection. In contrast, the proportion of *ATAD1*, *IFIT1*, and *IFIT2* to all other detected proteins was higher during acute than in persistent infection. FMDV polyprotein, *IFIT3* and *RSAD2* were not detected at the protein level at 28 dpi.

<table><caption>**Table 3.** Proteins with significantly different abundance at 24 hpi or 28 dpi. Only proteins that were also differentially expressed in the RNA sequencing experiment at the given time point are shown.</caption><thead><tr><th>Protein</th><th>Description</th><th>24 hpi<br>Mean NSAF</th><th>28 dpi<br>Mean NSAF</th><th>NSAF<sub>28 dpi</sub>/<br>NSAF<sub>24 hpi</sub></th></tr></thead><tbody><tr><td>ATAD1</td><td>ATPase Family, AAA Domain Containing 1</td><td>2.50 × 10<sup>-4</sup> **</td><td>1.56 × 10<sup>-4</sup></td><td>0.6</td></tr><tr><td>CAG23917.14</td><td>FMDV polyprotein</td><td>3.16 × 10<sup>-4</sup> **</td><td>0<sup>†</sup></td><td>n.a.</td></tr><tr><td>HERC6</td><td>HECT And RLD Domain Containing E3 Ubiquitin Protein Ligase<br>Family Member 6</td><td>2.46 × 10<sup>-5</sup></td><td>2.74 × 10<sup>-4</sup> **</td><td>11.2</td></tr><tr><td>IFI44</td><td>Interferon Induced Protein 44</td><td>1.12 × 10<sup>-4</sup> **</td><td>4.58 × 10<sup>-4</sup> **</td><td>4.1</td></tr><tr><td>IFI44L</td><td>Interferon Induced Protein 44 like</td><td>8.32 × 10<sup>-5</sup> **</td><td>6.98 × 10<sup>-4</sup> **</td><td>8.4</td></tr><tr><td>IFIT1</td><td>Interferon-Induced Protein with Tetratricopeptide Repeats 1</td><td>5.41 × 10<sup>-4</sup> **</td><td>7.62 × 10<sup>-5</sup> **</td><td>0.1</td></tr><tr><td>IFIT2</td><td>Interferon Induced Protein With Tetratricopeptide Repeats 2</td><td>2.31 × 10<sup>-4</sup> **</td><td>9.34 × 10<sup>-5</sup> **</td><td>0.4</td></tr><tr><td>IFIT3</td><td>Interferon Induced Protein With Tetratricopeptide Repeats 3</td><td>3.78 × 10<sup>-4</sup> **</td><td>0<sup>†</sup></td><td>n.a.</td></tr><tr><td>ISG15</td><td>ISG15 Ubiquitin-Like Modifier</td><td>3.53 × 10<sup>-3</sup> **</td><td>7.52 × 10<sup>-3</sup> **</td><td>2.1</td></tr><tr><td>MX1</td><td>MX Dynamin Like GTPase 1</td><td>8.42 × 10<sup>-4</sup> **</td><td>2.27 × 10<sup>-3</sup> **</td><td>2.7</td></tr><tr><td>MX2</td><td>MX Dynamin Like GTPase 2</td><td>2.44 × 10<sup>-4</sup> **</td><td>5.04 × 10<sup>-4</sup> **</td><td>2.1</td></tr><tr><td>OAS1X</td><td>2'-5'-Oligoadenylate Synthetase 1 X</td><td>2.47 × 10<sup>-4</sup> **</td><td>4.85 × 10<sup>-4</sup> **</td><td>2.0</td></tr><tr><td>OAS1Y</td><td>2'-5'-Oligoadenylate Synthetase 1 Y</td><td>2.31 × 10<sup>-4</sup> **</td><td>3.68 × 10<sup>-4</sup> **</td><td>1.6</td></tr><tr><td>RSAD2</td><td>Radical S-Adenosyl Methionine Domain Containing 2</td><td>5.69 × 10<sup>-4</sup> **</td><td>0<sup>†</sup></td><td>n.a.</td></tr></tbody></table>

## 4. Discussion

The failure of the host response to clear virus from the nasopharynx is one of the most important features of FMDV infection, and has been studied previously using monolayers of primary cells and non-primary cell lines from hamsters (BHK-21) [14], swine (SK6) [32], and cattle (EBK, MDBK and pharynx cells) [15–17]. These cells needed regular cell culture passage and were, except for the pharynx cells, not from anatomical locations relevant to natural FMDV infection and persistence, which resulted in a suboptimal reflection of the in vivo situation.

In this study we used multilayers of primary bovine dorsal soft palate (SP) epithelial cells in an air-liquid interface cell culture model [18], as it has been shown that this tissue, along with the dorsal nasopharynx, is the most likely site of primary FMDV infection and persistence in cattle [10,33]. This model mimics several properties of the squamous epithelium that are observed in vivo and allows detailed studies of the target cells of primary FMDV infection. Specialized immune cells, such as B cells, T cells, and natural killer cells, which contribute to host responses in vivo, were depleted by the establishment of the cell culture in air [18], leading to improved standardization and a focus on the primary target cell. In vivo, neutralizing antibodies secreted by B cells are essential for resolving systemic FMDV infection, but have no effect on persistent infection in the nasopharynx [9]. Recent studies of gene expression in nasopharyngeal tissues of cattle have suggested an important role of cytotoxic T cells in the clearance of persistent FMDV infection [34,35]. The transcriptional changes analyzed in the present study, however, are limited to the response of SP epithelial cells themselves to infection with FMDV.

This response was analyzed using RNA sequencing, RT-qPCR and quantitative proteomics. Based on the results from exploratory RNA sequencing alone, infection with FMDV resulted in detectable changes in the transcriptome of the cells and the induced changes were distinct during the acute and persistent phases of infection. Despite that only a small proportion of cells were persistently infected, these influenced the transcriptome significantly. The non-infected controls did not show any comparable, time-dependent changes, indicating stable culture conditions without extensive cell differentiation or degradation. The observed changes induced by FMDV infection were also

<footnote>*: *p*-value < 0.01; **: *p*-value < 0.001; †: not significant; n.a. not applicable.</footnote>
independent of the donor animals used for the preparation of the primary cell cultures, although transcriptional differences between the animals were detectable, e.g., sex-associated differences in XIST and keratin expression [36]. In summary, we observed genes that were highly regulated during both stages of infection (Figure 1C, green cluster), genes exclusively regulated during acute infection (Figure 1C, blue cluster) and animal-specific gene expression (Figure 1C, red cluster). During acute and persistent infection, comparably low numbers of genes were significantly differentially expressed and most of these were up-regulated. Interestingly, 63 genes were differentially expressed during both infection phases, indicating that FMDV infection induces long-lasting changes in the soft palate transcriptome. Furthermore, these genes were generally more up-regulated during acute infection. The transcriptional changes observed by RNA sequencing were supported at the protein level using quantitative proteomics. Although the number of detected differentially expressed proteins was much lower, it was confirmed that overexpression of ISGs leads to detectable levels of these proteins in the cells.

We found a strong activation of the innate immune response at 24 hpi that appears to be triggered by sensing of viral dsRNA over the cytosolic RNA sensors *IFIH1/MDA5* and *RIG-I/DDX58*, as the expression of *MDA5* and *RIG-I* was significantly increased during acute infection. While *RIG-I* is sensing minus-strand RNA viruses by specific binding to 5'-triphosphate uncapped RNA genomes, *MDA5* specifically senses plus-strand RNA viruses, such as FMDV, by binding to their dsRNA replication intermediates [37,38]. Although *MDA5* is thought to be the main cellular detector of FMDV, it has been shown that *RIG-I* transcription is also elevated during FMDV infection [39]. After specific binding of *MDA5* or *RIG-I* to their ligands they interact with *TRIM25* and the mitochondrial protein *MAVS* that activate the transcription factors *IRF3/IRF7* [40] as well as *NF-κB* [41], which is in accordance with the observed up-regulation of *IRF3*, *IRF7*, *NFKBIA*, *NFKBIB* and *TRIM25* (Table 1). Activated *IRF3/IRF7* and *NF-κB* complexes thereafter induce the expression of type I/III interferons (IFNs) and proinflammatory cytokines, respectively [42]. We found that during acute infection especially *IFN-β* was expressed at elevated levels with significant overexpression of *IFNB1*, while *IFN-α* and *IFN-λ* were only expressed at very low levels (Supplementary Figure S3). However, we observed expression of several ISGs with antiviral activity, such as *IRF1*, *ISG15*, *MX1/2*, *OAS1/2*, and *RSAD2* (see Table 1). This indicates that while *IFN-α* transcription may be limited, *IFN-β* alone is able to induce a potent antiviral response in bovine SP cells. It is known that the viral leader proteinase L^pro^ of FMDV can inhibit *IFN-β* transcription and protein translation, thereby blocking the cells innate immune response [43]. However, at least in this model L^pro^ does not induce a full blocking of *IFN-β* transcription. A specific blocking of *IFN-α* transcription has only been shown for swine dendritic cell populations from blood and skin during acute FMDV infection [44]. The induction of *IFN-β* and associated ISGs was furthermore coincident with high viral genome copy numbers and presence of viral proteins (Table 3) as has been described before [45].

In contrast, no interferon expression was observed during the persistent phase of infection, although the aforementioned virus sensors *MDA5/RIG-I* and their associated transcription factor *IRF7* were still highly expressed (see *IFIH1* in Figure 3B). A possible factor reducing the transcription of interferons during persistence is the down-regulation of *ANKRD1*, as *ANKRD1* is directly involved in the signal transduction of *IRF3* and *IRF7* by binding *IRF7/IRF3* complexes and thereby enhances the expression of type I/III interferons [46]. *ANKRD1* was initially observed in human dermal endothelial cells, where it is induced by inflammatory cytokines [47] and is substantially involved in the fibroblast mediated wound healing [48]. In agreement with our results, the expression of *ANKRD1* has previously been shown to be down-regulated in vivo during persistent FMDV infection of bovine nasopharyngeal tissue [35]. Gene silencing of *ANKRD1* in cells infected with herpes simplex virus resulted in increased viral load and reduced *IFN-β* (*IFNB1*) and *IFN-λ* (IL29) expression [46]. Therefore, the strong downregulation of *ANKRD1* as early as 48 hpi may be involved in the decrease of *IFN-β* expression during persistence.
Although no interferon expression was detectable at 28 dpi, the expression of many ISGs with known antiviral activity, such as *MX1/2* and *OAS1/2*, was significantly increased (see Figures 2D and 3B). Furthermore, based on the proteome data, the translation of these genes was stable, as their encoded proteins accumulated in the cells (Table 3). This indicates long-lasting induction of ISGs during persistent infection, with viral protein below the detection limit and low FMDV genome copy numbers. A comparable pattern of highly expressed ISGs as observed during FMDV persistence (*ISG15*, *MX1*, *OAS1/2*, and *USP18*) has been identified in liver samples from chronic infections with hepatitis C virus (HCV) [49,50]. Interestingly, the elevated expression of *ISG15/USP18* during chronic HCV infection, as we also observed during FMDV persistence, correlated with decreased responses to IFN-α treatment. Furthermore, the strong expression of *OAS1/2*, *MX1/2* and *IRF7*, but not of *IRF3* at 28 dpi is in accordance with increased interferon receptor (IFNR) signaling [51,52]. Observations from persistent lymphocytic choriomeningitis virus (LCMV) infections revealed that persistence is driven by chronic IFNR signaling, characterized by similar ISG activation patterns as observed in the present study [52]. The blockade of IFNR by antibodies abolished the expression of *IL-10* and *PD-L1*, two immunosuppressive T-cell exhaustion factors expressed by dendritic cells, and ultimately led to clearance of persistent LCMV infection by activation of IFN-γ expressing CD4+ T cells [52]. In this study, we did not observe upregulation of *IL-10* or *PD-L1*, as no dendritic cells were present in the SP cultures; however, an overexpression of T-cell exhaustion factors was previously found in nasopharyngeal tissues of FMDV carrier animals [35] and therefore may play a role in FMDV persistence in vivo.

The immunosuppressive factors identified in that study also included transforming growth factor β (TGF-β). While an overexpression of TGF-β itself was not evident in our data, lymphocyte antigen 6 family member E (*LY6E*) was one of the two genes that were up-regulated during the persistent phase only. *LY6E*, an ISG, has been implicated in the TGF-β-mediated escape from immune surveillance in many forms of cancer [53]. *LY6E* has been previously shown to promote viral infection [54] and it is essential for clathrin-mediated endocytosis of virus particles [55], a pathway that is also used by FMDV [56].

During acute infection, the interferon-mediated induction of ISGs appears to trigger apoptosis, as indicated by pathway analysis. Interestingly, genes associated with apoptosis, such as *CASP7*, were only expressed during acute infection and their expression waned with increasing time post infection (Figure 3B). The apparent absence of apoptotic processes during persistent infection is in accordance with recent findings that indicated an inhibition of apoptotic pathways in nasopharyngeal tissues from FMDV carrier animals [34,35]. Furthermore, we observed a specific down-regulation of pro-apoptotic genes during persistence, such as *ALDH1A2*, *ANKRD1*, and *SFRP2* (Table 2). Similarly, overexpression of fructose-bisphosphatase 1, the other gene that was only up-regulated during the persistent phase, inhibits many forms of apoptosis by increasing total cellular glutathione [57].

Another interesting observation was the downregulation of *NCAM1* (CD56), a member of the immunoglobulin superfamily, as early as 48 hpi and its further decrease until 28 dpi. Besides its role as a differentiation marker for natural killer cells, *NCAM1* is involved in cell binding and migration [58]. Downregulation of *NCAM1* is associated with decreased cell adhesion capacity, enhanced tumor cell invasiveness and is triggered in other viral infections [59]. The expression of *HTRA3*, a serine protease involved in remodeling of the extracellular matrix (ECM), as well as the expression of the ECM components *COL1A1* and *COL3A1* were significantly down-regulated during persistent infection. The role of these changes during FMDV persistence is currently unknown but lends itself to some speculation. The interaction of cells with the ECM plays a key role in epithelial maturation [60], which in turn is critical for the life cycle of some persistent viruses such as papillomaviruses [61]. Even though FMDV and papillomaviruses are biologically very different, certain features of FMDV persistence in vivo are conspicuously similar to what is observed in papillomavirus infection—particularly the difference between the distribution of viral genome, which is concentrated in the basal stratum germinativum [62], and viral antigen, which is concentrated in the superficial
layers of the soft palate epithelium [63]. This segregation helps papillomaviruses evade the immune
system as high levels of viral replication and protein synthesis occur only in terminally differentiated
cells that are not subject to immune surveillance [64]. Whether a similar mechanism is involved in
the maintenance of FMDV persistence remains to be investigated. However, it is difficult to faithfully
recreate the complex epithelial structure and ECM interactions in vitro and investigations of the role of
the ECM in persistent FMDV infection have to be performed with ex vivo tissue samples.

**5. Conclusions**

In conclusion, our study independently confirms earlier findings of a polygenic inhibition of apoptosis during persistent FMDV infection, which has been put forward as one of the principal mechanisms for the maintenance of persistence [34,35]. Another proposed mechanism, Th2 polarization and T-cell exhaustion, was not directly represented in our data, because the SP culture model does not include specialized immune cells. Nevertheless, we demonstrated the utility of state-of-the-art proteogenomics for the analysis of transcriptional signatures of acute and persistent FMDV infection in a near-natural in vitro system. We will proceed to apply this technology to tissue samples collected from carrier animals to obtain the first comprehensive picture of the transcriptomic and proteomic alterations associated with FMDV persistence in the natural host. Unraveling the cellular mechanisms of FMDV persistence may ultimately give rise to improved diagnostics and prevention of the FMDV carrier state.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/1999-4915/11/1/53/s1, Figure S1: Schematic visualization of the proteogenomic workflow, Figure S2: Comparison of different transformation methods for the raw read counts per gene, Figure S3: Comparison of interferon expression, Table S1: Summary of primary soft palate cell culture samples, Table S2: Summary of excluded and included RNA reference transcripts, Table S3: Replaced gene aliases of the *Bos taurus* reference genome, Table S4: Sequences of primers used for the RT-qPCR analysis, Table S5: Differentially expressed genes, File S1: Raw read counts for all significantly differently expressed genes.

**Author Contributions:** Writing—Original Draft, F.P. and M.E.; Methodology, S.H., F.P. and M.Z.; Investigation, F.P., S.H., M.Z., S.B., E.L., S.K., D.Z. and J.V.; Formal Analysis, Visualization, and Data Curation, F.P.; Supervision, S.Z., L.B., J.V., K.R., D.H., M.B. and M.E.; Conceptualization, S.Z., L.B., J.V., D.H., M.B. and M.E.; Project Administration, S.Z.

**Funding:** This work was supported by the Swedish Research Council (Formas, Sweden), L'Agence Nationale de la Recherche (ANR, France) and the Federal Office for Agriculture and Food (BLE, Germany) through the ANIHWA ERA-Net project TRANSCRIPTOVAC.

**Acknowledgments:** We thank Cells for Life Platform, partly funded by the Infrastructure Committee at SLU, Sweden, for providing facilities and equipment. We acknowledge Patrick Zitzow and Jenny Lorke of the Institute of Diagnostic Virology at the Friedrich-Loeffler-Institut for their excellent technical assistance. We are grateful to Dörte Becher and Claudia Hirschfeld from the Institute of Microbiology at the University of Greifswald for helping with the mass spectrometry analyses.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. Knight-Jones, T.J.; Rushton, J. The economic impacts of foot and mouth disease—What are they, how big are they and where do they occur? *Prev. Vet. Med.* **2013**, *112*, 161–173. [CrossRef] [PubMed]

2. Zell, R.; Delwart, E.; Gorbalenya, A.E.; Hovi, T.; King, A.M.Q.; Knowles, N.J.; Lindberg, A.M.; Pallansch, M.A.; Palmenberg, A.C.; Reuter, G.; et al. ICTV Virus Taxonomy Profile: Picornaviridae. *J. Gen. Virol.* **2017**, *98*, 2421–2422. [CrossRef] [PubMed]

3. Forss, S.; Strebel, K.; Beck, E.; Schaller, H. Nucleotide sequence and genome organization of foot-and-mouth disease virus. *Nucleic Acids Res.* **1984**, *12*, 6587–6601. [CrossRef] [PubMed]

4. Alexandersen, S.; Mowat, N. Foot-and-mouth disease: Host range and pathogenesis. *Curr. Top. Microbiol. Immunol.* **2005**, *288*, 9–42. [CrossRef] [PubMed]

5. Thomson, G.R.; Vosloo, W.; Bastos, A.D. Foot and mouth disease in wildlife. *Virus Res.* **2003**, *91*, 145–161. [CrossRef]

6. Grubman, M.J.; Baxt, B. Foot-and-mouth disease. *Clin. Microbiol. Rev.* **2004**, *17*, 465–493. [CrossRef] [PubMed]
7. World Organisation for Animal Health (OIE). Foot and mouth disease (infection with foot and mouth disease virus). In *Manual of Diagnostic Tests and Vaccines for Terrestrial Animals (Terrestrial Manual)*; OIE: Paris, France, 2018; 2.1.8.

8. Arzt, J.; Belsham, G.J.; Lohse, L.; Botner, A.; Stenfeldt, C. Transmission of Foot-and-Mouth Disease from Persistently Infected Carrier Cattle to Naive Cattle via Transfer of Oropharyngeal Fluid. *mSphere* **2018**, *3*. [CrossRef]

9. Alexandersen, S.; Zhang, Z.; Donaldson, A.I. Aspects of the persistence of foot-and-mouth disease virus in animals—The carrier problem. *Microbes Infect.* **2002**, *4*, 1099–1110. [CrossRef]

10. Pacheco, J.M.; Smoliga, G.R.; O'Donnell, V.; Brito, B.P.; Stenfeldt, C.; Rodriguez, L.L.; Arzt, J. Persistent Foot-and-Mouth Disease Virus Infection in the Nasopharynx of Cattle; Tissue-Specific Distribution and Local Cytokine Expression. *PLoS ONE* **2015**, *10*, e0125698. [CrossRef]

11. Cox, S.J.; Voyce, C.; Parida, S.; Reid, S.M.; Hamblin, P.A.; Paton, D.J.; Barnett, P.V. Protection against direct-contact challenge following emergency FMD vaccination of cattle and the effect on virus excretion from the oropharynx. *Vaccine* **2005**, *23*, 1106–1113. [CrossRef]

12. Cox, S.J.; Voyce, C.; Parida, S.; Reid, S.M.; Hamblin, P.A.; Hutchings, G.; Paton, D.J.; Barnett, P.V. Effect of emergency FMD vaccine antigen payload on protection, sub-clinical infection and persistence following direct contact challenge of cattle. *Vaccine* **2006**, *24*, 3184–3190. [CrossRef] [PubMed]

13. Juleff, N.; Windsor, M.; Reid, E.; Seago, J.; Zhang, Z.; Monaghan, P.; Morrison, I.W.; Charleston, B. Foot-and-mouth disease virus persists in the light zone of germinal centres. *PLoS ONE* **2008**, *3*, e3434. [CrossRef] [PubMed]

14. Han, L.; Xin, X.; Wang, H.; Li, J.; Hao, Y.; Wang, M.; Zheng, C.; Shen, C. Cellular response to persistent foot-and-mouth disease virus infection is linked to specific types of alterations in the host cell transcriptome. *Sci. Rep.* **2018**, *8*, 5074. [CrossRef] [PubMed]

15. Zhu, J.; Weiss, M.; Grubman, M.J.; de los Santos, T. Differential gene expression in bovine cells infected with wild type and leaderless foot-and-mouth disease virus. *Virology* **2010**, *404*, 32–40. [CrossRef]

16. Kopliku, L.; Relmy, A.; Romey, A.; Gorna, K.; Zientara, S.; Bakkali-Kassimi, L.; Blaise-Boisseau, S. Establishment of persistent foot-and-mouth disease virus (FMDV) infection in MDBK cells. *Arch. Virol.* **2015**, *160*, 2503–2516. [CrossRef] [PubMed]

17. O'Donnell, V.; Pacheco, J.M.; Larocco, M.; Gladue, D.P.; Pauszek, S.J.; Smoliga, G.; Krug, P.W.; Baxt, B.; Borca, M.V.; Rodriguez, L. Virus-host interactions in persistently FMDV-infected cells derived from bovine pharynx. *Virology* **2014**, *468–470*, 185–196. [CrossRef] [PubMed]

18. Hägglund, S.; Laloy, E.; Näslund, K.; Pfaff, F.; Eschbaumer, M.; Romey, A.; Relmy, A.; Rikberg, A.; Svensson, A.; Huet, H.; et al. Model of persistent foot-and-mouth disease virus infection in multilayered cells derived from bovine dorsal soft palate. *Transbound. Emerg. Dis.* **2018**, under review.

19. Patro, R.; Duggal, G.; Love, M.I.; Irizarry, R.A.; Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. *Nat. Methods* **2017**, *14*, 417–419. [CrossRef]

20. R Core Team. *R: A Language and Environment for Statistical Computing*; R Foundation for Statistical Computing: Vienna, Austria, 2017.

21. Soneson, C.; Love, M.I.; Robinson, M.D. Differential analyses for RNA-seq: Transcript-level estimates improve gene-level inferences. *F1000Research* **2015**, *4*, 1521. [CrossRef]

22. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biol.* **2014**, *15*, 550. [CrossRef]

23. Yu, G.; He, Q.Y. ReactomePA: An R/Bioconductor package for reactome pathway analysis and visualization. *Mol. Biosyst.* **2016**, *12*, 477–479. [CrossRef] [PubMed]

24. Yu, G.; Wang, L.G.; Han, Y.; He, Q.Y. clusterProfiler: An R package for comparing biological themes among gene clusters. *Omics* **2012**, *16*, 284–287. [CrossRef] [PubMed]

25. Untergasser, A.; Cutcutache, I.; Koressaar, T.; Ye, J.; Faircloth, B.C.; Remm, M.; Rozen, S.G. Primer3—New capabilities and interfaces. *Nucleic Acids Res.* **2012**, *40*, e115. [CrossRef] [PubMed]

26. Eymann, C.; Lassek, C.; Wegner, U.; Bernhardt, J.; Fritsch, O.A.; Fuchs, S.; Otto, A.; Albrecht, D.; Schiefelbein, U.; Cernava, T.; et al. Symbiotic Interplay of Fungi, Algae, and Bacteria within the Lung Lichen *Lobaria pulmonaria* L. Hoffm. as Assessed by State-of-the-Art Metaproteomics. *J. Proteome Res.* **2017**, *16*, 2160–2173. [CrossRef] [PubMed]
27. López-Mondéjar, R.; Zühlke, D.; Větrovský, T.; Becher, D.; Riedel, K.; Baldrian, P. Decoding the complete arsenal for cellulose and hemicellulose deconstruction in the highly efficient cellulose decomposer Paenibacillus O199. *Biotechnol. Biofuels* **2016**, *9*, 104. [CrossRef] [PubMed]

28. Huber, W.; Carey, V.J.; Gentleman, R.; Anders, S.; Carlson, M.; Carvalho, B.S.; Bravo, H.C.; Davis, S.; Gatto, L.; Girke, T.; et al. Orchestrating high-throughput genomic analysis with Bioconductor. *Nat. Methods* **2015**, *12*, 115–121. [CrossRef] [PubMed]

29. Pavelka, N.; Pelizzola, M.; Vizzardelli, C.; Capozzoli, M.; Splendiani, A.; Granucci, F.; Ricciardi-Castagnoli, P. A power law global error model for the identification of differentially expressed genes in microarray data. *BMC Bioinform.* **2004**, *5*, 203. [CrossRef]

30. Vizcaino, J.A.; Csordas, A.; del-Toro, N.; Dianes, J.A.; Griss, J.; Lavidas, I.; Mayer, G.; Perez-Riverol, Y.; Reisinger, F.; Ternent, T.; et al. 2016 update of the PRIDE database and its related tools. *Nucleic Acids Res.* **2016**, *44*, D447–D456. [CrossRef]

31. Fabregat, A.; Jupe, S.; Matthews, L.; Sidiropoulos, K.; Gillespie, M.; Garapati, P.; Haw, R.; Jassal, B.; Korninger, F.; May, B.; et al. The Reactome Pathway Knowledgebase. *Nucleic Acids Res.* **2018**, *46*, D649–D655. [CrossRef]

32. Ni, Z.; Yang, F.; Cao, W.; Zhang, X.; Jin, Y.; Mao, R.; Du, X.; Li, W.; Guo, J.; Liu, X.; et al. Differential gene expression in porcine SK6 cells infected with wild-type and SAP domain-mutant foot-and-mouth disease virus. *Virol. Sin.* **2016**, *31*, 249–257. [CrossRef]

33. Arzt, J.; Juleff, N.; Zhang, Z.; Rodriguez, L.L. The pathogenesis of foot-and-mouth disease I: Viral pathways in cattle. *Transbound. Emerg. Dis.* **2011**, *58*, 291–304. [CrossRef] [PubMed]

34. Stenfeldt, C.; Eschbaumer, M.; Smoliga, G.R.; Rodriguez, L.L.; Zhu, J.; Arzt, J. Clearance of a persistent picornavirus infection is associated with enhanced pro-apoptotic and cellular immune responses. *Sci. Rep.* **2017**, *7*, 17800. [CrossRef] [PubMed]

35. Eschbaumer, M.; Stenfeldt, C.; Smoliga, G.R.; Pacheco, J.M.; Rodriguez, L.L.; Li, R.W.; Zhu, J.; Arzt, J. Transcriptomic Analysis of Persistent Infection with Foot-and-Mouth Disease Virus in Cattle Suggests Impairment of Apoptosis and Cell-Mediated Immunity in the Nasopharynx. *PLoS ONE* **2016**, *11*, e0162750. [CrossRef] [PubMed]

36. Brown, C.J.; Hendrich, B.D.; Rupert, J.L.; Lafreniere, R.G.; Xing, Y.; Lawrence, J.; Willard, H.F. The human XIST gene: Analysis of a 17 kb inactive X-specific RNA that contains conserved repeats and is highly localized within the nucleus. *Cell* **1992**, *71*, 527–542. [CrossRef]

37. Feng, Q.; Hato, S.V.; Langereis, M.A.; Zoll, J.; Virgen-Slane, R.; Peisley, A.; Hur, S.; Semler, B.L.; van Rij, R.P.; van Kuppeveld, F.J. MDA5 detects the double-stranded RNA replicative form in picornavirus-infected cells. *Cell Rep.* **2012**, *2*, 1187–1196. [CrossRef]

38. Kato, H.; Takeuchi, O.; Sato, S.; Yoneyama, M.; Yamamoto, M.; Matsui, K.; Uematsu, S.; Jung, A.; Kawai, T.; Ishii, K.J.; et al. Differential roles of MDA5 and RIG-I helicases in the recognition of RNA viruses. *Nature* **2006**, *441*, 101–105. [CrossRef]

39. Zhu, Z.; Wang, G.; Yang, F.; Cao, W.; Mao, R.; Du, X.; Zhang, X.; Li, C.; Li, D.; Zhang, K.; et al. Foot-and-Mouth Disease Virus Viroporin 2B Antagonizes RIG-I-Mediated Antiviral Effects by Inhibition of Its Protein Expression. *J. Virol.* **2016**, *90*, 11106–11121. [CrossRef]

40. Lin, R.; Heylbroeck, C.; Pitha, P.M.; Hiscott, J. Virus-dependent phosphorylation of the IRF-3 transcription factor regulates nuclear translocation, transactivation potential, and proteasome-mediated degradation. *Mol. Cell. Biol.* **1998**, *18*, 2986–2996. [CrossRef]

41. Liu, S.; Chen, J.; Cai, X.; Wu, J.; Chen, X.; Wu, Y.T.; Sun, L.; Chen, Z.J. MAVS recruits multiple ubiquitin E3 ligases to activate antiviral signaling cascades. *Elife* **2013**, *2*, e00785. [CrossRef]

42. Belgnaoui, S.M.; Paz, S.; Hiscott, J. Orchestrating the interferon antiviral response through the mitochondrial antiviral signaling (MAVS) adapter. *Curr. Opin. Immunol.* **2011**, *23*, 564–572. [CrossRef]

43. De Los Santos, T.; de Avila Botton, S.; Weiblen, R.; Grubman, M.J. The leader proteinase of foot-and-mouth disease virus inhibits the induction of beta interferon mRNA and blocks the host innate immune response. *J. Virol.* **2006**, *80*, 1906–1914. [CrossRef] [PubMed]

44. Nfon, C.K.; Ferman, G.S.; Toka, F.N.; Gregg, D.A.; Golde, W.T. Interferon-alpha production by swine dendritic cells is inhibited during acute infection with foot-and-mouth disease virus. *Viral. Immunol.* **2008**, *21*, 68–77. [CrossRef] [PubMed]
45. Eschbaumer, M.; Stenfeldt, C.; Rekant, S.I.; Pacheco, J.M.; Hartwig, E.J.; Smoliga, G.R.; Kenney, M.A.; Golde, W.T.; Rodriguez, L.L.; Arzt, J. Systemic immune response and virus persistence after foot-and-mouth disease virus infection of naive cattle and cattle vaccinated with a homologous adenovirus-vectored vaccine. *BMC Vet. Res.* **2016**, *12*, 205. [CrossRef] [PubMed]

46. Bin, L.; Li, X.; Richers, B.; Streib, J.E.; Hu, J.W.; Taylor, P.; Leung, D.Y.M. Ankyrin repeat domain 1 regulates innate immune responses against herpes simplex virus 1: A potential role in eczema herpeticum. *J. Allergy Clin. Immunol.* **2018**. [CrossRef] [PubMed]

47. Chu, W.; Burns, D.K.; Swerlick, R.A.; Presky, D.H. Identification and characterization of a novel cytokine-inducible nuclear protein from human endothelial cells. *J. Biol. Chem.* **1995**, *270*, 10236-10245. [CrossRef]

48. Samaras, S.E.; Almodovar-Garcia, K.; Wu, N.; Yu, F.; Davidson, J.M. Global deletion of Ankrd1 results in a wound-healing phenotype associated with dermal fibroblast dysfunction. *Am. J. Pathol.* **2015**, *185*, 96-109. [CrossRef] [PubMed]

49. Chen, L.; Borozan, I.; Feld, J.; Sun, J.; Tannis, L.L.; Coltescu, C.; Heathcote, J.; Edwards, A.M.; McGilvray, I.D. Hepatic gene expression discriminates responders and nonresponders in treatment of chronic hepatitis C viral infection. *Gastroenterology* **2005**, *128*, 1437-1444. [CrossRef] [PubMed]

50. MacQuillan, G.C.; Mamotte, C.; Reed, W.D.; Jeffrey, G.P.; Allan, J.E. Upregulation of endogenous intrahepatic interferon stimulated genes during chronic hepatitis C virus infection. *J. Med. Virol.* **2003**, *70*, 219-227. [CrossRef]

51. Au, W.C.; Moore, P.A.; Lowther, W.; Juang, Y.T.; Pitha, P.M. Identification of a member of the interferon regulatory factor family that binds to the interferon-stimulated response element and activates expression of interferon-induced genes. *Proc. Natl. Acad. Sci. USA* **1995**, *92*, 11657-11661. [CrossRef]

52. Wilson, E.B.; Yamada, D.H.; Elsaesser, H.; Herskovitz, J.; Deng, J.; Cheng, G.; Aronow, B.J.; Karp, C.L.; Brooks, D.G. Blockade of chronic type I interferon signaling to control persistent LCMV infection. *Science* **2013**, *340*, 202-207. [CrossRef]

53. AlHossiny, M.; Luo, L.; Frazier, W.R.; Steiner, N.; Gusev, Y.; Kallakury, B.; Glasgow, E.; Creswell, K.; Madhavan, S.; Kumar, R.; et al. Ly6E/K Signaling to TGFbeta Promotes Breast Cancer Progression, Immune Escape, and Drug Resistance. *Cancer Res.* **2016**, *76*, 3376-3386. [CrossRef] [PubMed]

54. Schoggins, J.W.; Wilson, S.J.; Panis, M.; Murphy, M.Y.; Jones, C.T.; Bieniasz, P.; Rice, C.M. A diverse range of gene products are effectors of the type I interferon antiviral response. *Nature* **2011**, *472*, 481-485. [CrossRef] [PubMed]

55. Hackett, B.A.; Cherry, S. Flavivirus internalization is regulated by a size-dependent endocytic pathway. *Proc. Natl. Acad. Sci. USA* **2018**, *115*, 4246-4251. [CrossRef] [PubMed]

56. Han, S.C.; Guo, H.C.; Sun, S.Q.; Jin, Y.; Wei, Y.Q.; Feng, X.; Yao, X.P.; Cao, S.Z.; Liu, D.X.; Liu, X.T. Productive Entry of Foot-and-Mouth Disease Virus via Macropinocytosis Independent of Phosphatidylinositol 3-Kinase. *Sci. Rep.* **2016**, *6*, 19294. [CrossRef] [PubMed]

57. Voehringer, D.W.; Hirschberg, D.L.; Xiao, J.; Lu, Q.; Roederer, M.; Lock, C.B.; Herzenberg, L.A.; Steinman, L.; Herzenberg, L.A. Gene microarray identification of redox and mitochondrial elements that control resistance or sensitivity to apoptosis. *Proc. Natl. Acad. Sci. USA* **2000**, *97*, 2680-2685. [CrossRef] [PubMed]

58. Moos, M.; Tacke, R.; Scherer, H.; Teplow, D.; Früh, K.; Schachner, M. Neural adhesion molecule L1 as a member of the immunoglobulin superfamily with binding domains similar to fibronectin. *Nature* **1988**, *334*, 701-703. [CrossRef] [PubMed]

59. Blaheta, R.A.; Beecken, W.D.; Engl, T.; Jonas, D.; Oppermann, E.; Hundemer, M.; Doerr, H.W.; Scholz, M.; Cinatl, J. Human cytomegalovirus infection of tumor cells downregulates NCAM (CD56): A novel mechanism for virus-induced tumor invasiveness. *Neoplasia* **2004**, *6*, 323-331. [CrossRef]

60. Bonnans, C.; Chou, J.; Werb, Z. Remodelling the extracellular matrix in development and disease. *Nat. Rev. Mol. Cell Biol.* **2014**, *15*, 786-801. [CrossRef]

61. Feller, L.; Khammissa, R.A.; Wood, N.H.; Lemmer, J. Epithelial maturation and molecular biology of oral HPV. *Infect. Agent Cancer* **2009**, *4*, 16. [CrossRef]

62. Prato Murphy, M.L.; Forsyth, M.A.; Belsham, G.J.; Salt, J.S. Localization of foot-and-mouth disease virus RNA by in situ hybridization within bovine tissues. *Virus Res.* **1999**, *62*, 67-76. [CrossRef]
63. Stenfeldt, C.; Eschbaumer, M.; Rekant, S.I.; Pacheco, J.M.; Smoliga, G.R.; Hartwig, E.J.; Rodriguez, L.L.; Arzt, J. The Foot-and-Mouth Disease Carrier State Divergence in Cattle. *J. Virol.* **2016**, *90*, 6344–6364. [CrossRef] [PubMed]

64. McBride, A.A. Mechanisms and strategies of papillomavirus replication. *Biol. Chem.* **2017**, *398*, 919–927. [CrossRef] [PubMed]

<figure><img src="image_8.png" /></figure>

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).