Cellular response to persistent foot-and-mouth disease virus infection is linked to specific types of alterations in the host cell transcriptome

Figure 1. Emergence of FMDV-resistant cells during persistent infection. (a) BHK-21 cells and BHK-VECs were infected with FMDV at 2.5 × 10−4 PFU/cell for 18 h, and whole-cell extracts were analyzed by western blotting using rabbit polyclonal anti-3D antibody and mouse monoclonal anti-GAPDH. (b) BHK-VECs were infected with FMDV at 5 × 10−2 PFU/cell. Infected cells were subcultured, and intracellular viral RNA of each generation was measured using qRT-PCR. (c) Expression of FMDV integrin receptors (ITGB6, ITGAV, ITGB1) in BHK-VEC, BHK-Op, and uninfected BHK-21 cell cultures were detected by western blotting using specific antibodies. (d) Adsorption of FMDV to BHK-21 cells and BHK-VECs. BHK21 cells and BHK-VECs were infected with FMDV at 2.8 PFU/cell and maintained at 0–4 °C. At the indicated time point, supernatant was removed rapidly, and monolayers were washed with cold PBS three times. One mL serum-free MEM was added to the monolayers. Cell suspensions were then subjected to three cycles of freezing and thawing and centrifuged for 5 min at 12,000 × g at 4 °C. Supernatant was used to quantify the bound virus titer by TCID50 assay. (e) Detection of intracellular viral RNA replication in…

Figure 2. BHK-VECs resisted infection of FMDV-Op. (a) BHK-21 cells and BHK-VECs were infected with FMDV-Op at 2 × 10−4 PFU/cell for 24 h. Whole-cell extracts were analyzed by western blotting using rabbit polyclonal anti-3D-specific or mouse monoclonal anti-GAPDH antibodies. (b) BHK-VECs were infected with FMDV-Op at 2 × 10−4 PFU/cell. Infected cells were subcultured, and the secreted infectious virus in the supernatant of each generation was measured by titration assay. (c) BHK-21 cells were infected with FMDV or FMDV-Op at 2 × 10−4 PFU/cell. Intracellular RNA was isolated for quantification of intracellular virus RNA by qRT-PCR at the indicated time points. ***p < 0.001.

Figure 3. Analysis of DEGs in BHK-VECs. (a) GO enrichment analysis of DEGs in BHK-VECs. We included the top 30 most significant GO categories of DEGs (*p < 0.05). KEGG pathway enrichment analysis was performed for differentially expressed (b) up-regulated and (c) down-regulated genes in BHK-VECs. The graph depicts the most significant 20 enriched pathway entries. (d) Validation of DEGs in BHK-VECs by (i) qRT-PCR or (ii) western blotting. For all qRT-PCR and western blot validations, GAPDH served as the internal reference gene to normalize data.

Figure 4. AS induced in BHK-VECs. (a) Summary of significant differential AS events in BHK-VECs compared with BHK-21 cells. The intron-exon structure involved in each splicing pattern is shown. (b–e) Verification of AS events by RT-PCR. Visualization of RNA-seq reads was first performed by IGV. Then, RT-PCR was performed using specific primers located upstream and downstream of SE events. The y-axis shows the number of mapped reads, and the red arrows at the bottom show the target area for primer amplification. Hnrnpa2b1, Pde4dip, and Hypk showed increased SEs in BHK-VECs as shown by a larger number of products of 324, 390, and 144 bp. Syt12 showed a decrease in SEs in BHK-VECs as shown by a larger number of products of 227 bp.

Figure 5. Differential expression of immunity-related genes in BHK-VECs. (a) Differential expression of innate immunity-related genes between BHK-21 cells and BHK-VECs were artificially chosen for expression analysis. The x-axis is the periphery of the circle showing the names of genes, and the y-axis is the radius of the circle showing the corresponding FPKM values. The farther away from the center of the circle, the higher the expression of the gene. (b) Comparison of adaptive immunity-related gene expression between BHK-21 cells and BHK-VECs. (c) Hierarchical clustering heatmap of differentially expressed immunity-related genes between BHK-VECs and BHK-21 cells. The color range from green to red corresponds to a low to high abundance of gene expression, respectively. (d,e) Validation of differentially expressed immunity-related genes in BHK-VECs by qRT-PCR or western blotting. GAPDH served as the internal reference gene to normalize data.

Figure 6. Differential expression of MAPK signaling-related genes in BHK-VECs. (a) qRT-PCR validation of DEGs related to MAPK pathways in BHK-VECs. GAPDH served as the internal reference gene to normalize data, and three biologically independent replicates were performed. (b,c) Effect of MAPK/ERK or p38/MAPK inhibition on replication of FMDV in BHK-21 cells. BHK-21 cells were pre-incubated (1 h) with DMSO, (b) 20 or 50 mM U0126, or (c) 20 or 50 mM SB202190 and then infected with FMDV at 2.5 × 10−4 PFU/cell for 24 h in the presence of DMSO, U0126, or SB202190. Protein extracts were examined using western blotting with FMDV 3D-specific or the indicated antibodies; the effectiveness of inhibition was monitored by detecting the phosphorylation of inhibitor-specific target protein (P-MAPK/ERK and pS15-Hsp27). (d) (i) p38δ (MAPK13) was down-regulated in BHK-VECs compared with BHK-21 cells, and (ii) overexpression of MAPK13 genes in BHK-21 cells promoted the replication of FMDV.