Transcriptomic analysis of the dialogue between Pseudorabies virus and porcine epithelial cells during infection

Advances in the immunoescape mechanisms exploited by alphaherpesviruses

Alphaherpesviruses, categorized as viruses with linear DNA composed of two complementary strands, can potentially to induce diseases in both humans and animals as pathogens. Mature viral particles comprise of a core, capsid, tegument, and envelope. While herpesvirus infection can elicit robust immune and inflammatory reactions in the host, its persistence stems from its prolonged interaction with the host, fostering a diverse array of immunoescape mechanisms. In recent years, significant advancements have been achieved in comprehending the immunoescape tactics employed by alphaherpesviruses, including pseudorabies virus (PRV), herpes simplex virus (HSV), varicella-zoster virus (VZV), feline herpesvirus (FeHV), equine herpesvirus (EHV), and caprine herpesvirus type I (CpHV-1). Researchers have unveiled the intricate adaptive mechanisms existing between viruses and their natural hosts. This review endeavors to illuminate the research advancements concerning the immunoescape mechanisms of alphaherpesviruses by delineating the pertinent proteins and genes involved in virus immunity. It aims to furnish valuable insights for further research on related mechanisms and vaccine development, ultimately contributing to virus control and containment efforts.

Acute transcriptomic changes in murine RAW 264.7 cells following pseudorabies virus infection

Next-generation sequencing enables the evaluation of gene expression changes resulting from virus-host interactions at the RNA level. Pseudorabies virus (PRV) causes substantial economic loss in the swine industry. Recent research has revealed that PRV can be transmitted to and infect humans as well. To identify physiopathological and pathological responses post-PRV infection, we characterized transcriptomic changes in the murine RAW 264.7 cell line over the course of 36 h. In total, 156, 153, and 190 differentially expressed genes were identified at 2 h, 12 h, and 36 h, respectively. Seven differentially expressed genes (Trim27, Ccdc117, Mrps12, Ccl4, Cerkl, Ubald1, and Hmga1-rs1) were present across all treatment groups. Our findings expand our knowledge of gene regulation and immune response following PRV infection. These differentially expressed genes can subsequently improve our understanding of PRV pathogenesis.

Pseudorabies Virus Infection Results in a Broad Inhibition of Host Gene Transcription

Pseudorabies virus (PRV) is a porcine alphaherpesvirus that belongs to the Herpesviridae family. We showed earlier that infection of porcine epithelial cells with PRV triggers activation of the nuclear factor κB (NF-κB) pathway, a pivotal signaling axis in the early immune response. However, PRV-induced NF-κB activation does not lead to NF-κB-dependent gene expression. Here, using electrophoretic mobility shift assays (EMSAs), we show that PRV does not disrupt the ability of NF-κB to interact with its κB target sites. Assessing basal cellular transcriptional activity in PRV-infected cells by quantitation of prespliced transcripts of constitutively expressed genes uncovered a broad suppression of cellular transcription by PRV, which also affects the inducible expression of NF-κB target genes. Host cell transcription inhibition was rescued when viral genome replication was blocked using phosphonoacetic acid (PAA). Remarkably, we found that host gene expression shutoff in PRV-infected cells correlated with a substantial retention of the NF-κB subunit p65, the TATA box binding protein, and RNA polymerase II-essential factors required for (NF-κB-dependent) gene transcription-in expanding PRV replication centers in the nucleus and thereby away from the host chromatin. This study reveals a potent mechanism used by the alphaherpesvirus PRV to steer the protein production capacity of infected cells to viral proteins by preventing expression of host genes, including inducible genes involved in mounting antiviral responses. IMPORTANCE Herpesviruses are highly successful pathogens that cause lifelong persistent infections of their host. Modulation of the intracellular environment of infected cells is imperative for the success of virus infections. We reported earlier that a DNA damage response in epithelial cells infected with the alphaherpesvirus pseudorabies virus (PRV) results in activation of the hallmark proinflammatory NF-κB signaling axis but, remarkably, that this activation does not lead to NF-κB-induced (proinflammatory) gene expression. Here, we report that PRV-mediated inhibition of host gene expression stretches beyond NF-κB-dependent gene expression and in fact reflects a broad inhibition of host gene transcription, which correlates with a substantial recruitment of essential host transcription factors in viral replication compartments in the nucleus, away from the host chromatin. These data uncover a potent alphaherpesvirus mechanism to interfere with production of host proteins, including proteins involved in antiviral responses.

Phospho-proteomics identifies a critical role of ATF2 in pseudorabies virus replication

Pseudorabies virus (PRV), an etiological agent of pseudorabies in livestock, has negatively affected the porcine industry all over the world. Epithelial cells are reported as the first site of PRV infection. However, the role of host proteins and its related signaling pathways in PRV replication is largely unclear. In this study, we performed a quantitative phosphoproteomics screening on PRV-infected porcine kidney (PK-15) epithelial cells. Totally 5723 phosphopeptides, corresponding to 2180 proteins, were obtained, and the phosphorylated states of 810 proteins were significantly different in PRV-infected cells compared with mock-infected cells (P ​< ​0.05). GO and KEGG analysis revealed that these differentially expressed phosphorylated proteins were predominantly related to RNA transport and MAPK signaling pathways. Further functional studies of NF-κB, transcription activator factor-2 (ATF2), MAX and SOS genes in MAPK signaling pathway were analyzed using RNA interference (RNAi) knockdown. It showed that only ATF2-knockdown reduces both PRV titer and viral genome copy number. JNK pathway inhibition and CRISPR/Cas9 gene knockout showed that ATF2 was required for the effective replication of PRV, especially during the biogenesis of viral genome DNA. Subsequently, by overexpression of the ATF2 gene and point mutation of the amino acid positions 69/71 of ATF2, it was further demonstrated that the phosphorylation of ATF2 promoted PRV replication. These findings suggest that ATF2 may provide potential therapeutic target for inhibiting PRV infection.

•

Phosphoproteomic profiling of PRV-infected PK-15 ​cells with iTRAQ-quantification.

•

JNK pathway regulates ATF2 phosphorylation and PRV replication.

•

Phosphorylation of ATF2 promotes PRV replication.

Overexpression of Sedoheptulose-1,7-Bisphosphatase Enhances Photosynthesis in Chlamydomonas reinhardtii and Has No Effect on the Abundance of Other Calvin-Benson Cycle Enzymes

The productivity of plants and microalgae needs to be increased to feed the growing world population and to promote the development of a low-carbon economy. This goal can be achieved by improving photosynthesis via genetic engineering. In this study, we have employed the Modular Cloning strategy to overexpress the Calvin-Benson cycle (CBC) enzyme sedoheptulose-1,7-bisphosphatase (SBP1) up to threefold in the unicellular green alga Chlamydomonas reinhardtii. The protein derived from the nuclear transgene represented ∼0.3% of total cell protein. Photosynthetic rate and growth were significantly increased in SBP1-overexpressing lines under high-light and elevated CO2 conditions. Absolute quantification of the abundance of all other CBC enzymes by the QconCAT approach revealed no consistent differences between SBP1-overexpressing lines and the recipient strain. This analysis also revealed that the 11 CBC enzymes represent 11.9% of total cell protein in Chlamydomonas. Here, the range of concentrations of CBC enzymes turned out to be much larger than estimated earlier, with a 128-fold difference between the most abundant CBC protein (rbcL) and the least abundant (triose phosphate isomerase). Accordingly, the concentrations of the CBC intermediates are often but not always higher than the binding site concentrations of the enzymes for which they act as substrates. The enzymes with highest substrate to binding site ratios might represent good candidates for overexpression in subsequent engineering steps.

Functional analysis of prv-miR-LLT11a encoded by pseudorabies virus

Viral-encoded microRNAs (miRNAs) have vital roles in the regulation of virus replications and host immune responses. The results of previous studies have indicated that miRNA clusters are involved in the replication and virulence of the pseudorabies virus (PRV), which may potentially lead to immune escape or facilitation of PRV replication. This study's previous research revealed that prv-miR-LLT11a was differentially expressed during PRV infection. The present study's results have demonstrated that prv-miR-LLT11a could significantly inhibit PRV replication. It was further determined that SLA-1 was the target gene of prv-miR-LLT11a, and simultaneously, that overexpression of prv-miR-LLT11a could downregulate the mRNA and protein levels of SLA-1 in a dose-independent manner. Furthermore, the present study also observed that prv-miR-LLT11a can downregulate TAP1 expression. Our findings provide a better understanding of the molecular mechanism involved in the effects of prv-miR-LLT11a on SLA-1 and TAP1 as well as its involvement in immune system evasion of PRV.

Unified feature association networks through integration of transcriptomic and proteomic data

High-throughput multi-omics studies and corresponding network analyses of multi-omic data have rapidly expanded their impact over the last 10 years. As biological features of different types (e.g. transcripts, proteins, metabolites) interact within cellular systems, the greatest amount of knowledge can be gained from networks that incorporate multiple types of -omic data. However, biological and technical sources of variation diminish the ability to detect cross-type associations, yielding networks dominated by communities comprised of nodes of the same type. We describe here network building methods that can maximize edges between nodes of different data types leading to integrated networks, networks that have a large number of edges that link nodes of different-omic types (transcripts, proteins, lipids etc). We systematically rank several network inference methods and demonstrate that, in many cases, using a random forest method, GENIE3, produces the most integrated networks. This increase in integration does not come at the cost of accuracy as GENIE3 produces networks of approximately the same quality as the other network inference methods tested here. Using GENIE3, we also infer networks representing antibody-mediated Dengue virus cell invasion and receptor-mediated Dengue virus invasion. A number of functional pathways showed centrality differences between the two networks including genes responding to both GM-CSF and IL-4, which had a higher centrality value in an antibody-mediated vs. receptor-mediated Dengue network. Because a biological system involves the interplay of many different types of molecules, incorporating multiple data types into networks will improve their use as models of biological systems. The methods explored here are some of the first to specifically highlight and address the challenges associated with how such multi-omic networks can be assembled and how the greatest number of interactions can be inferred from different data types. The resulting networks can lead to the discovery of new host response patterns and interactions during viral infection, generate new hypotheses of pathogenic mechanisms and confirm mechanisms of disease.

Novel nucleotide variants in SLA-DOB and CD4 are associated with immune traits in pregnant sows

Reinforcing the immunity of pregnant sows can not only improve their own health condition but also increase the survival rate and healthy status of their piglets. This study aims to find single-nucleotide polymorphism (SNP) and molecular markers that are associated with the immune traits of pregnant sows. SLA-DOB and CD4 were selected as candidate genes, and blood samples were randomly collected from pregnant Landrace sows and used to detect T-lymphocyte subsets, interferon alpha, interleukin 6, Toll-like receptor 3, serum antibody immunoglobulin G, and porcine reproductive and respiratory syndrome virus-specific antibody. Then, association analyses were conducted for the polymorphic sites of candidate genes with immune traits. We found 12 mutations in the two genes and conducted an association study with eight of them. Our results indicated that among the eight mutations, SNP1, SNP2, and SNP3 of the SLA-DOB gene and Ins9, SNP10, and SNP11 in the CD4 gene are newly discovered mutations. Except for SNP1, SNP3, and SNP11, the other five SNPs are associated with at least one immune trait tested. Especially, SNP2 and Ins9 are significantly associated with at least one of the T-lymphocyte subgroups and at least one antibody. These novel mutations have potential important effects on the polymorphic loci of the above immune traits in pregnant sows. The results suggest that the SLA-DOB and CD4 genes and their genetic mutations can be considered as important candidate genes and mutations for the immunity of pregnant sows.

Interactions between the transcription and replication machineries regulate the RNA and DNA synthesis in the herpesviruses

The temporal coordination of viral gene expression is imperative for the regulation of the herpesvirus replication cycle. While the main factors of this transcriptional coordination are known, the subtler control mechanisms of gene expression remain elusive. Recent long read sequencing-based approached have revealed an intricate meshwork of overlaps between the herpesvirus transcripts and the overlap of the replication origins with noncoding RNAs. It has been shown that the transcriptional apparatuses can physically interfere with one another while transcribing overlapping regions. We hypothesize that transcriptional interference regulates the global gene expression across the herpesvirus genome. Additionally, an overall decrease in transcriptional activity in individual viral genes has been observed following the onset of DNA replication. An overlap of the replication origins with specific transcripts has also been described in several herpesviruses. The genome-wide interactions between the transcriptional apparatuses and between the replication and transcriptional machineries suggest the existence of novel layers of genetic regulation.

Identification of DNA methylation regulated novel host genes relevant to inhibition of virus replication in porcine PK15 cell using double stranded RNA mimics and DNA methyltransferase inhibitor

During RNA viruses's replication, double-stranded RNA (dsRNA) is normally produced and induce host innate immune response. Most of gene activation due cytokine mediated but which are due to methylation mediated is still unknown. In the study, DNA methylome was integrated with our previous transcriptome data to investigate the differentially methylated regions and genes using MeDIP-chip technology. We found that the transcriptional expressions of 15, 37 and 18 genes were negatively related with their promoter DNA methylation levels in the cells treated by PolyI:C, Aza-CdR, as well as PolyI:C plus Aza-CdR, respectively, compared with the untreated cells. GO analysis revealed hypo-methylated genes (BNIP3L and CDK9) and a hyper-methylated gene (ZC3HAV1) involved in the host response to viral replication. Our results suggest that these novel genes targeted by DNA methylation can be potential markers relevant to virus replication and host innate immune response to set up a medical model of infectious diseases.

Transcriptomics Sequencing Provides Insights into Understanding the Mechanism of Grass Carp Reovirus Infection

Grass carp is an important aquaculture fish species in China that is affected by severe diseases, especially haemorrhagic disease caused by grass carp reovirus (GCRV). However, the mechanisms of GCRV invasion and infection remain to be elucidated. In the present study, Ctenopharyngodon idellus kidney (CIK) cells were infected with GCRV, harvested at 0, 8, 24, and 72 h post infection, respectively, and then subjected to transcriptomics sequencing. Each sample yielded more than 6 Gb of clean data and 40 million clean reads. To better understand GCRV infection, the process was divided into three phases: the early (0-8 h post infection), middle (8-24 h post infection), and late (24-72 h) stages of infection. A total of 76 (35 up-regulated, 41 down-regulated), 553 (463 up-regulated, 90 down-regulated), and 284 (150 up-regulated, 134 down-regulated) differently expressed genes (DEGs) were identified during the early, middle, and late stages of infection, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that DEGs were mainly involved in carbohydrate biosynthesis, transport, and endocytosis in the early stage, phagocytosis and lysosome pathways were mainly enriched in the middle stage, and programmed cell death, apoptosis, and inflammation were largely associated with the late stage. These results suggest GCRV infection is a gradual process involving adsorption on the cell surface, followed by endocytosis into cells, transport by lysosomes, and eventually resulted in cell necrosis and/or apoptosis. Our findings provide insight into the mechanisms of grass carp reovirus infection.

Evaluation of the impact of ul54 gene-deletion on the global transcription and DNA replication of pseudorabies virus

Pseudorabies virus (PRV) is an animal alphaherpesvirus with a wide host range. PRV has 67 protein-coding genes and several non-coding RNA molecules, which can be classified into three temporal groups, immediate early, early and late classes. The ul54 gene of PRV and its homolog icp27 of herpes simplex virus have a multitude of functions, including the regulation of viral DNA synthesis and the control of the gene expression. Therefore, abrogation of PRV ul54 function was expected to exert a significant effect on the global transcriptome and on DNA replication. Real-time PCR and real-time RT-PCR platforms were used to investigate these presumed effects. Our analyses revealed a drastic impact of the ul54 mutation on the genome-wide expression of PRV genes, especially on the transcription of the true late genes. A more than two hour delay was observed in the onset of DNA replication, and the amount of synthesized DNA molecules was significantly decreased in comparison to the wild-type virus. Furthermore, in this work, we were able to successfully demonstrate the utility of long-read SMRT sequencing for genotyping of mutant viruses.

Validation of a novel Mho microarray for a comprehensive characterisation of the Mycoplasma hominis action in HeLa cell infection

Mycoplasma hominis is the second smallest facultative pathogen of the human urogenital tract. With less than 600 protein-encoding genes, it represents an ideal model organism for the study of host-pathogen interactions. For a comprehensive characterisation of the M. hominis action in infection a customized Mho microarray, which was based on two genome sequences (PG21 and LBD-4), was designed to analyze the dynamics of the mycoplasma transcriptome during infection and validated for M. hominis strain FBG. RNA preparation was evaluated and adapted to ensure the highest recovery of mycoplasmal mRNAs from in vitro HeLa cell infection assays. Following cRNA hybridization, the read-out strategy of the hybridization results was optimized and confirmed by RT-PCR. A statistically robust infection assay with M. hominis strain FBG enabled the identification of differentially regulated key effector molecules such as critical cytoadhesins (4 h post infection (pI)), invasins (48 h pI) and proteins associated with establishing chronic infection of the host (336 h pI). Of the 294 differentially regulated genes (>2-fold) 128 (43.5%) encoded hypothetical proteins, including lipoproteins that seem to play a central role as virulence factors at each stage of infection: P75 as a novel cytoadhesin candidate, which is also differentially upregulated in chronic infection; the MHO_2100 protein, a postulated invasin and the MHO_730-protein, a novel ecto-nuclease and domain of an ABC transporter, the function of which in chronic infection has still to be elucidated. Implementation of the M. hominis microarray strategy led to a comprehensive identification of to date unknown candidates for virulence factors at relevant stages of host cell infection.

iTRAQ-based Proteomic Analysis of Porcine Kidney Epithelial PK15 cells Infected with Pseudorabies virus

Pseudorabies virus (PRV) is one of the most important pathogens of swine, resulting in severe economic losses to the pig industry. To improve our understanding of the host responses to PRV infection, we applied isobaric tags for relative and absolute quantification (iTRAQ) labeling coupled with liquid chromatography-tandem mass spectrometry to quantitatively identify the differentially expressed cellular proteins in PRV-infected PK15 cells. In total, relative quantitative data were identified for 4333 proteins in PRV and mock- infected PK15 cells, among which 466 cellular proteins were differentially expressed, including 234 upregulated proteins and 232 downregulated proteins. Bioinformatics analysis disclosed that most of these differentially expressed proteins were involved in metabolic processes, cellular growth and proliferation, endoplasmic reticulum (ER) stress response, cell adhesion and cytoskeleton. Moreover, expression levels of four representative proteins, beta-catenin, STAT1, GRB2 and PCNA, were further confirmed by western blot analysis. This is the first attempt to analyze the protein profile of PRV-infected PK15 cells using iTRAQ technology, and our findings may provide valuable information to help understand the host response to PRV infection.

Gene Expression Profiling with Cre-Conditional Pseudorabies Virus Reveals a Subset of Midbrain Neurons That Participate in Reward Circuitry

The mesolimbic dopamine pathway receives inputs from numerous regions of the brain as part of a neural system that detects rewarding stimuli and coordinates a behavioral response. The capacity to simultaneously map and molecularly define the components of this complex multisynaptic circuit would thus advance our understanding of the determinants of motivated behavior. To accomplish this, we have constructed pseudorabies virus (PRV) strains in which viral propagation and fluorophore expression are activated only after exposure to Cre recombinase. Once activated in Cre-expressing neurons, the virus serially labels chains of presynaptic neurons. Dual injection of GFP and mCherry tracing viruses simultaneously illuminates nigrostriatal and mesolimbic circuitry and shows no overlap, demonstrating that PRV transmission is confined to synaptically connected neurons. To molecularly profile mesolimbic dopamine neurons and their presynaptic inputs, we injected Cre-conditional GFP virus into the NAc of (anti-GFP) nanobody-L10 transgenic mice and immunoprecipitated translating ribosomes from neurons infected after retrograde tracing. Analysis of purified RNA revealed an enrichment of transcripts expressed in neurons of the dorsal raphe nuclei and lateral hypothalamus that project to the mesolimbic dopamine circuit. These studies identify important inputs to the mesolimbic dopamine pathway and further show that PRV circuit-directed translating ribosome affinity purification can be broadly applied to identify molecularly defined neurons comprising complex, multisynaptic circuits.SIGNIFICANCE STATEMENT The mesolimbic dopamine circuit integrates signals from key brain regions to detect and respond to rewarding stimuli. To further define this complex multisynaptic circuit, we constructed a panel of Cre recombinase-activated pseudorabies viruses (PRVs) that enabled retrograde tracing of neural inputs that terminate on Cre-expressing neurons. Using these viruses and Retro-TRAP (translating ribosome affinity purification), a previously reported molecular profiling method, we developed a novel technique that provides anatomic as well as molecular information about the neural components of polysynaptic circuits. We refer to this new method as PRV-Circuit-TRAP (PRV circuit-directed TRAP). Using it, we have identified major projections to the mesolimbic dopamine circuit from the lateral hypothalamus and dorsal raphe nucleus and defined a discrete subset of transcripts expressed in these projecting neurons, which will allow further characterization of this important pathway. Moreover, the method we report is general and can be applied to the study of other neural circuits.

Characterization of the Dynamic Transcriptome of a Herpesvirus with Long-read Single Molecule Real-Time Sequencing

Herpesvirus gene expression is co-ordinately regulated and sequentially ordered during productive infection. The viral genes can be classified into three distinct kinetic groups: immediate-early, early, and late classes. In this study, a massively parallel sequencing technique that is based on PacBio Single Molecule Real-time sequencing platform, was used for quantifying the poly(A) fraction of the lytic transcriptome of pseudorabies virus (PRV) throughout a 12-hour interval of productive infection on PK-15 cells. Other approaches, including microarray, real-time RT-PCR and Illumina sequencing are capable of detecting only the aggregate transcriptional activity of particular genomic regions, but not individual herpesvirus transcripts. However, SMRT sequencing allows for a distinction between transcript isoforms, including length- and splice variants, as well as between overlapping polycistronic RNA molecules. The non-amplified Isoform Sequencing (Iso-Seq) method was used to analyse the kinetic properties of the lytic PRV transcripts and to then classify them accordingly. Additionally, the present study demonstrates the general utility of long-read sequencing for the time-course analysis of global gene expression in practically any organism.

gD-Independent Superinfection Exclusion of Alphaherpesviruses

Many viruses have the capacity to prevent a cell from being infected by a second virus, often termed superinfection exclusion. Alphaherpesviruses, including the human pathogen herpes simplex virus 1 (HSV-1) and the animal herpesvirus pseudorabies virus (PRV), encode a membrane-bound glycoprotein, gD, that can interfere with subsequent virion entry. We sought to characterize the timing and mechanism of superinfection exclusion during HSV-1 and PRV infection. To this end, we utilized recombinant viruses expressing fluorescent protein (FP) markers of infection that allowed the visualization of viral infections by microscopy and flow cytometry as well as the differentiation of viral progeny. Our results demonstrated the majority of HSV-1- and PRV-infected cells establish superinfection exclusion by 2 h postinfection. The modification of viral infections by virion inactivation and phosphonoacetic acid, cycloheximide, and actinomycin D treatments indicated new protein synthesis is needed to establish superinfection exclusion. Primary infection with gene deletion PRV recombinants identified that new gD expression is not required to establish superinfection exclusion of a secondary viral inoculum. We also identified the timing of coinfection events during axon-to-cell spread, with most occurring within a 2-h window, suggesting a role for cellular superinfection exclusion during neuroinvasive spread of infection. In summary, we have characterized a gD-independent mechanism of superinfection exclusion established by two members of the alphaherpesvirus family and identified a potential role of exclusion during the pathogenic spread of infection.

IMPORTANCE

Superinfection exclusion is a widely observed phenomenon initiated by a primary viral infection to prevent further viruses from infecting the same cell. The capacity for alphaherpesviruses to infect the same cell impacts rates of interviral recombination and disease. Interviral recombination allows genome diversification, facilitating the development of resistance to antiviral therapeutics and evasion of vaccine-mediated immune responses. Our results demonstrate superinfection exclusion occurs early, through a gD-independent process, and is important in the directed spread of infection. Identifying when and where in an infected host viral genomes are more likely to coinfect the same cell and generate viral recombinants will enhance the development of effective antiviral therapies and interventions.

Effects of Pseudorabies Virus Infection on the Tracheobronchial Lymph Node Transcriptome

This study represents the first swine transcriptome hive plots created from gene set enrichment analysis (GSEA) data and provides a novel insight into the global transcriptome changes occurring in tracheobronchial lymph nodes (TBLN) and spanning the swine genome. RNA isolated from draining TBLN from 5-week-old pigs, either clinically infected with a feral isolate of Pseudorabies virus or uninfected, was interrogated using Illumina Digital Gene Expression Tag Profiling. More than 100 million tag sequences were observed, representing 4,064,189 unique 21-base sequences collected from TBLN at time points 1, 3, 6, and 14 days post-inoculation (dpi). Multidimensional statistical tests were applied to determine the significant changes in tag abundance, and then the tags were annotated. Hive plots were created to visualize the differential expression within the swine transcriptome defined by the Broad Institute’s GSEA reference datasets between infected and uninfected animals, allowing us to directly compare different conditions.

Current transcriptomics in pig immunity research

Swine performance in the face of disease challenge is becoming progressively more important. To improve the pig’s robustness and resilience against pathogens through selection, a better understanding of the genetic and epigenetic factors in the immune response is required. This review highlights results from the most recent transcriptome research, and the meta-analyses performed, in the context of pig immunity. A technological overview is given including wholegenome microarrays, immune-specific arrays, small-scale high-throughput expression methods, high-density tiling arrays, and next generation sequencing (NGS). Although whole genome microarray techniques will remain complementary to NGS for some time in domestic species, research will transition to sequencing-based methods due to cost-effectiveness and the extra information that such methods provide. Furthermore, upcoming high-throughput epigenomic studies, which will add greatly to our knowledge concerning the impact of epigenetic modifications on pig immune response, are listed in this review. With emphasis on the insights obtained from transcriptomic analyses for porcine immunity, we also discuss the experimental design in pig immunity research and the value of the newly published porcine genome assembly in using the pig as a model for human immune response. We conclude by discussing the importance of establishing community standards to maximize the possibility of integrative computational analyses, such as was clearly beneficial for the human ENCODE project.

Exposure to zearalenone mycotoxin alters in vitro porcine intestinal epithelial cells by differential gene expression

The gut represents the main route of intoxication with mycotoxins. To evaluate the effect and the underlying molecular changes that occurred when the intestine is exposed to zearalenone, a Fusarium sp mycotoxin, porcine epithelial cells (IPEC-1) were treated with 10μM of ZEA for 24h and analysed by microarray using Gene Spring GX v.11.5. Our results showed that 10μM of ZEA did not affect cell viability, but can increase the expression of toll like receptors (TLR1-10) and of certain cytokines involved in inflammation (TNF-α, IL-1β, IL-6, IL-8, MCP-1, IL-12p40, CCL20) or responsible for the recruitment of immune cells (IL-10, IL-18). Microarray results identified 190 genes significantly and differentially expressed, of which 70% were up-regulated. ZEA determined the over expression of ITGB5 gene, essential against the attachment and adhesion of ETEC to porcine jejunal cells and of TFF2 implicated in mucosal protection. An up-regulation of glutathione peroxidase enzymes (GPx6, GPx2, GPx1) was also observed. Upon ZEA challenge, genes like GTF3C4 responsible for the recruitment of polymerase III and initiation of tRNA transcription in eukaryotes and STAT5B were significantly higher induced. The up-regulation of CD97 gene and the down-regulation of tumour suppressor genes (DKK-1, PCDH11X and TC531386) demonstrates the carcinogenic potential of ZEA.

The peripheral blood transcriptome reflects variations in immunity traits in swine: towards the identification of biomarkers

BACKGROUND

Immune traits (ITs) are potentially relevant criteria to characterize an individual's immune response. Our aim was to investigate whether the peripheral blood transcriptome can provide a significant and comprehensive view of IT variations in pig.

RESULTS

Sixty-day-old Large White pigs classified as extreme for in vitro production of IL2, IL10, IFNγ and TNFα, phagocytosis activity, in vivo CD4⁻/CD8⁺ or TCRγδ + cell counts, and anti-Mycoplasma antibody levels were chosen to perform a blood transcriptome analysis with a porcine generic array enriched with immunity-related genes. Differentially expressed (DE) genes for in vitro production of IL2 and IL10, phagocytosis activity and CD4⁻/CD8⁺ cell counts were identified. Gene set enrichment analysis revealed a significant over-representation of immune response functions. To validate the microarray-based results, a subset of DE genes was confirmed by RT-qPCR. An independent set of 74 animals was used to validate the covariation between gene expression levels and ITs. Five potential gene biomarkers were found for prediction of IL2 (RALGDS), phagocytosis (ALOX12) or CD4⁻/CD8⁺ cell count (GNLY, KLRG1 and CX3CR1). On average, these biomarkers performed with a sensitivity of 79% and a specificity of 86%.

CONCLUSIONS

Our results confirmed that gene expression profiling in blood represents a relevant molecular phenotype to refine ITs in pig and to identify potential biomarkers that can provide new insights into immune response analysis.

Genome-wide effects of DNA methyltransferase inhibitor on gene expression in double-stranded RNA transfected porcine PK15 cells

Double-stranded RNA (dsRNA) is produced in host cells during viral replication. The effects of DNA demethylation on gene expression in dsRNA transfected swine cells are unclear. The study aims to profile the transcriptome changes which are induced by DNA methyltransferase inhibitor (Aza-CdR) in porcine PK15 cells transfected with viral-like dsRNA (Poly(I:C)). A total of 44, 76 and 952 differentially expressed genes (DEGs) were detected in the cells treated by Poly(I:C) plus Aza-CdR (P+A), Poly(I:C) (P) or Aza-CdR (A) alone compared to the controls (C). Immune response-related pathways are observed in the comparison of A vs. C and P vs. C, and the genes in the pathways were recovered in the comparison of (P+A) vs. C. GO analysis indicated that Aza-CdR has negative regulatory effects on viral reproduction. The results suggest that the stimulant of Poly(I:C) could be regressed by Aza-CdR. These observations provide new insights into the epigenetic regulatory effects on viral replication.

Physiological change under OsHV-1 contamination in Pacific oyster Crassostrea gigas through massive mortality events on fields

BACKGROUND

Massive mortalities have been observed in France since 2008 on spat and juvenile Pacific oysters, Crassostrea gigas. A herpes virus called OsHV-1, easily detectable by PCR, has been implicated in the mortalities as demonstrated by the results of numerous field studies linking mortality with OsHV-1 prevalence. Moreover, experimental infections using viral particles have documented the pathogenicity of OsHV-1 but the physiological responses of host to pathogen are not well known.

RESULTS

The aim of this study was to understand mechanisms brought into play against the virus during infection in the field. A microarray assay has been developed for a major part of the oyster genome and used for studying the host transcriptome across mortality on field. Spat with and without detectable OsHV-1 infection presenting or not mortality respectively were compared by microarray during mortality episodes. In this study, a number of genes are regulated in the response to pathogen infection on field and seems to argue to an implication of the virus in the observed mortality. The result allowed establishment of a hypothetic scheme of the host cell's infection by, and response to, the pathogen.

CONCLUSIONS

This response shows a "sensu stricto" innate immunity through genic regulation of the virus OsHV-1 life cycle, but also others biological processes resulting to complex interactions between host and pathogens in general.

Structural and functional annotation of the porcine immunome

BACKGROUND

The domestic pig is known as an excellent model for human immunology and the two species share many pathogens. Susceptibility to infectious disease is one of the major constraints on swine performance, yet the structure and function of genes comprising the pig immunome are not well-characterized. The completion of the pig genome provides the opportunity to annotate the pig immunome, and compare and contrast pig and human immune systems.

RESULTS

The Immune Response Annotation Group (IRAG) used computational curation and manual annotation of the swine genome assembly 10.2 (Sscrofa10.2) to refine the currently available automated annotation of 1,369 immunity-related genes through sequence-based comparison to genes in other species. Within these genes, we annotated 3,472 transcripts. Annotation provided evidence for gene expansions in several immune response families, and identified artiodactyl-specific expansions in the cathelicidin and type 1 Interferon families. We found gene duplications for 18 genes, including 13 immune response genes and five non-immune response genes discovered in the annotation process. Manual annotation provided evidence for many new alternative splice variants and 8 gene duplications. Over 1,100 transcripts without porcine sequence evidence were detected using cross-species annotation. We used a functional approach to discover and accurately annotate porcine immune response genes. A co-expression clustering analysis of transcriptomic data from selected experimental infections or immune stimulations of blood, macrophages or lymph nodes identified a large cluster of genes that exhibited a correlated positive response upon infection across multiple pathogens or immune stimuli. Interestingly, this gene cluster (cluster 4) is enriched for known general human immune response genes, yet contains many un-annotated porcine genes. A phylogenetic analysis of the encoded proteins of cluster 4 genes showed that 15% exhibited an accelerated evolution as compared to 4.1% across the entire genome.

CONCLUSIONS

This extensive annotation dramatically extends the genome-based knowledge of the molecular genetics and structure of a major portion of the porcine immunome. Our complementary functional approach using co-expression during immune response has provided new putative immune response annotation for over 500 porcine genes. Our phylogenetic analysis of this core immunome cluster confirms rapid evolutionary change in this set of genes, and that, as in other species, such genes are important components of the pig's adaptation to pathogen challenge over evolutionary time. These comprehensive and integrated analyses increase the value of the porcine genome sequence and provide important tools for global analyses and data-mining of the porcine immune response.

An assessment of opportunities to dissect host genetic variation in resistance to infectious diseases in livestock

This paper reviews the evidence for host genetic variation in resistance to infectious diseases for a wide variety of diseases of economic importance in poultry, cattle, pig, sheep and Atlantic salmon. Further, it develops a method of ranking each disease in terms of its overall impact, and combines this ranking with published evidence for host genetic variation and information on the current state of genomic tools in each host species. The outcome is an overall ranking of the amenability of each disease to genomic studies that dissect host genetic variation in resistance. Six disease-based assessment criteria were defined: industry concern, economic impact, public concern, threat to food safety or zoonotic potential, impact on animal welfare and threat to international trade barriers. For each category, a subjective score was assigned to each disease according to the relative strength of evidence, impact, concern or threat posed by that particular disease, and the scores were summed across categories. Evidence for host genetic variation in resistance was determined from available published data, including breed comparison, heritability studies, quantitative trait loci (QTL) studies, evidence of candidate genes with significant effects, data on pathogen sequence and on host gene expression analyses. In total, 16 poultry diseases, 13 cattle diseases, nine pig diseases, 11 sheep diseases and three Atlantic salmon diseases were assessed. The top-ranking diseases or pathogens, i.e. those most amenable to studies dissecting host genetic variation, were Salmonella in poultry, bovine mastitis, Marek's disease and coccidiosis, both in poultry. The top-ranking diseases or pathogens in pigs, sheep and Atlantic salmon were Escherichia coli, mastitis and infectious pancreatic necrosis, respectively. These rankings summarise the current state of knowledge for each disease and broadly, although not entirely, reflect current international research efforts. They will alter as more information becomes available and as genome tools become more sophisticated for each species. It is suggested that this approach could be used to rank diseases from other perspectives as well, e.g. in terms of disease control strategies.

Viral infection resistance conferred on mice by siRNA transgenesis

RNA interference is an attractive strategy to fight against viral diseases by targeting the mRNA of viral genes. Most studies have reported the transient delivery of small interfering RNA or small hairpin (shRNA) expression constructs. Here, we present the production of transgenic mice stably expressing shRNA or miRNA targeting the IE180 mRNA (immediate early gene) of the pseudorabies virus (PRV) which infects mice and farm animals. We firstly designed non-retroviral shRNA or miRNA expression vectors. Secondly, we selected the most efficient shRNA construct that targeted either the 5'part or 3'UTR of the IE mRNA and was able to knockdown the target gene expression in cultured cells, by measuring systematically the shRNA content and comparing this with the interfering effects. We then produced four lines of transgenic mice expressing different amounts of shRNA or miRNA in the brain but without signs of stimulation of innate immunity. Lastly, we tested their resistance to PRV infection. In all transgenic lines, we observed a significant resistance to viral challenge, the best being achieved with the shRNA construct targeting the 3'UTR of the IE gene. Viral DNA levels in the brains of infected mice were always lower in transgenic mice, even in animals that did not survive. Finally, this work reports an effective strategy to generate transgenic animals producing shRNA from non-retroviral expression vectors. Moreover, these mice are the first transgenic animal models producing shRNA with a significant antiviral effect but without any apparent shRNA toxicity.

Microarray analysis of mediastinal lymph node of pigs naturally affected by postweaning multisystemic wasting syndrome

Postweaning multisystemic wasting syndrome (PMWS) is one of the pig diseases with major economic impact worldwide. Clinical, pathologic and some immunologic aspects of this disease are relatively well-known, but the molecular mechanisms underlying the pathogenic mechanisms of the disease are still poorly understood. The objective of the present study was to investigate the global transcriptome changes in the mediastinal lymph nodes from pigs naturally affected by PMWS, as well as healthy counterparts, using the Affymetrix Porcine Genechip(®). From 366 transcripts showing significant differential abundance in the PMWS group of pigs relative to healthy animals, 229 showed higher and 137 lower abundance. A relative increased abundance of mRNAs coded by a large set of genes involved in the inflammatory responses (e.g. cytokines, acute phase proteins, and respiratory burst) was observed in PMWS affected pigs. The Gpnmb and Lgals3 genes, which have antagonistic functions in regulation of inflammatory processes, showed high mRNA levels in diseased pigs. The complement system was altered by PMWS, notably by the lower levels of Cr1 mRNA, which might favour both complement deposition and secondary infections by impairing phagocytosis. Decreased mRNA abundance of several genes involved in lymphocyte activation/differentiation, such as Cd79b, Cd19, Cd21 and MybL1, and the high level of Vsig4 mRNA, which can compromise the activation of residing T-cells, pointed towards a defective adaptive immunity. This is the first study on gene expression in pigs naturally affected by PMWS. The present results allowed identifying potential mechanisms underlying the inflammation and lymphocyte depletion in lymphoid tissues by complement mediated damage and immunosuppression, which are key features of PMWS.

Effects of deletion of the early protein 0 gene of pseudorabies virus on the overall viral gene expression

Real-time RT-PCR analysis was applied to evaluate the impact of deletion of the early protein 0 (EP0) gene of pseudorabies virus (PRV) on the global expression of the viral transcripts during lytic infection in cultured porcine kidney cells. Our analysis showed that EP0 exerted an inhibitory effect on the transcription of the PRV genes in the early stage of infection, and alternating stimulatory and inhibitory effects on the viral gene expressions in the late stage of infection. The data also suggested that a general function of EP0 might be to reverse the kinetics of expression of early viral genes. We also observed that EP0 facilitated the development of correlations in the transcription kinetics between the immediate early 180 gene and the PRV transcripts, indicating that a major function of EP0 could be to modify the effects of the IE180 protein on the PRV transcriptome.

Investigation of the relationship between SLA-1 and SLA-3 gene expression and susceptibility to Escherichia coli F18 in post-weaning pigs

Porcine post-weaning diarrhea and edema disease are principally caused by Escherichia coli strains that produce F18 adhesin. FUT1 genotyping and receptor binding studies divided piglets into E. coli F18-resistant and -sensitive groups, and the roles of SLA-1 and SLA-3 were investigated. SLA-1 and SLA-3 expression was detected in 11 pig tissues, with higher levels of SLA-1 in lung, immune tissues and gastrointestinal tract, and higher levels of SLA-3 also in lung and lymphoid tissues. Both genes were expressed higher in F18-resistant piglets, and their expression was positively correlated in different tissues; a negative correlation was observed in some tissues of F18-sensitive group, particularly in lung and lymphatic samples. Gene ontology and pathway analyses showed that SLA-1 and SLA-3 were involved in 37 biological processes, including nine pathways related to immune functions. These observations help to elucidate the relationship between SLA class I genes and E. coli F18-related porcine gastrointestinal tract diseases.

Localization, expression change in PRRSV infection and association analysis of the porcine TAP1 gene

The transporter associated with antigen processing (TAP) translocates antigenic peptides from the cytosol into the lumen of the endoplasmic reticular and plays a critical role in the major histocompatibility complex (MHC) class I molecule-mediated antigenic presentation pathway. In this study, the porcine TAP1 gene was mapped to the pig chromosome 7 (SSC7) and was closely linked to the marker SSC2B02 (retention fraction=43%, LOD=15.18). Subcellular localization of TAP1 by transient transfection of PK15 cells indicated that the TAP1 protein might be located in the endoplasmic reticulum (ER) in pig kidney epithelial cells (PK-15). Gene expression analysis by semi-quantitative RT-PCR revealed that TAP1 was selectively expressed in some immune and immune-related tissues. Quantitative real-time PCR (qRT-PCR) analysis revealed that this gene was up-regulated after treatments that mimic viral and bacterial infection (polyriboinosinic-polyribocytidylic acid (poly(I:C)) and lipopolysaccharide (LPS), respectively). In addition, elevated TAP1 expression was detected after porcine reproductive and respiratory syndrome virus (PRRSV) infection in porcine white blood cells (WBCs). One single nucleotide polymorphism (SNP) in exon 3 of TAP1 was detected in a Landrace pig population by Bsp143I restriction enzyme digestion. Different genotypes of this SNP had significant associations (P<0.05) with the red blood cell distribution width (RDW) of 1-day-old (1 d) pigs (P=0.0168), the PRRSV antibody level (PRRSV Ab) (P=0.0445) and the absolute lymphocyte count (LYM#) (P=0.024) of 17 d pigs. Our results showed that the TAP1 gene might have important roles in swine immune responses, and these results provide useful information for further functional studies.

Comparison of gene expression and genome-wide DNA methylation profiling between phenotypically normal cloned pigs and conventionally bred controls

Animal breeding via Somatic Cell Nuclear Transfer (SCNT) has enormous potential in agriculture and biomedicine. However, concerns about whether SCNT animals are as healthy or epigenetically normal as conventionally bred ones are raised as the efficiency of cloning by SCNT is much lower than natural breeding or In-vitro fertilization (IVF). Thus, we have conducted a genome-wide gene expression and DNA methylation profiling between phenotypically normal cloned pigs and control pigs in two tissues (muscle and liver), using Affymetrix Porcine expression array as well as modified methylation-specific digital karyotyping (MMSDK) and Solexa sequencing technology. Typical tissue-specific differences with respect to both gene expression and DNA methylation were observed in muscle and liver from cloned as well as control pigs. Gene expression profiles were highly similar between cloned pigs and controls, though a small set of genes showed altered expression. Cloned pigs presented a more different pattern of DNA methylation in unique sequences in both tissues. Especially a small set of genomic sites had different DNA methylation status with a trend towards slightly increased methylation levels in cloned pigs. Molecular network analysis of the genes that contained such differential methylation loci revealed a significant network related to tissue development. In conclusion, our study showed that phenotypically normal cloned pigs were highly similar with normal breeding pigs in their gene expression, but moderate alteration in DNA methylation aspects still exists, especially in certain unique genomic regions.

Molecular characterization of transcriptome-wide interactions between highly pathogenic porcine reproductive and respiratory syndrome virus and porcine alveolar macrophages in vivo

Porcine reproductive and respiratory syndrome virus (PRRSV) infects mainly the porcine alveolar macrophages (PAMs) and causes porcine reproductive and respiratory syndrome (PRRS). Previous studies have analyzed the global gene expression profiles of lung tissue in vivo and PAMs in vitro following infection with PRRSV, however, transcriptome-wide understanding of the interaction between highly pathogenic PRRSV (HP-PRRSV) and PAMs in vivo has not yet been established. In this study, we employed Affymetrix microarrays to investigate the gene expression patterns of PAMs isolated from Tongcheng piglets (a Chinese indigenous breed) after infection with HP-PRRSV. During the infection, Tongcheng piglets exhibited typical clinical signs, e.g. fever, asthma, coughing, anorexia, lethargy and convulsion, but displayed mild regional lung damage at 5 and 7 dpi. Microarray analysis revealed that HP-PRRSV infection has affected PAMs in expression of the important genes involved in cytoskeleton and exocytosis organization, protein degradation and folding, intracellular calcium and zinc homeostasis. Several potential antiviral strategies might be employed in PAMs, including upregulating IFN-induced genes and increasing intracellular zinc ion concentration. And inhibition of the complement system likely attenuated the lung damage during HP-PRRSV infection. Transcriptomic analysis of PAMs in vivo could lead to a better understanding of the HP-PRRSV-host interaction, and to the identification of novel antiviral therapies and genetic components of swine tolerance/susceptibility to HP-PRRS.

Analysis of porcine MHC using microarrays

The major histocompatibility complex (MHC) in Mammals is one of the most gene dense regions of the genome and contains the polymorphic histocompatibility gene families known to be involved in pathogen response and control of auto-immunity. The MHC is a complex genetic system that provides an interesting model system to study genome expression regulation and genetic diversity at the megabase scale. The pig MHC or SLA (Swine Leucocyte Antigen) complex spans 2.4 megabases and 151 loci have been annotated. We will review key results from previous RNA expression studies using microarrays containing probes specific to annotated loci within SLA and in addition present novel data obtained using high-density tiling arrays encompassing the whole SLA complex. We have focused on transcriptome modifications of porcine peripheral blood mononuclear cells stimulated with a mixture of phorbol myristate acetate and ionomycin known to activate B and T cell proliferation. Our results show that numerous loci mapping to the SLA complex are affected by the treatment. A general decreased level of expression for class I and II genes and an up-regulation of genes involved in peptide processing and transport were observed. Tiling array-based experiments contributed to refined gene annotations as presented for one SLA class I gene referred to as SLA-11. In conclusion, high-density tiling arrays can serve as an excellent tool to draw comprehensive transcription maps, and improve genome annotations for the SLA complex. We are currently studying their relevance to characterize SLA genetic diversity in combination with high throughput next generation sequencing.

Co-expression of host and viral microRNAs in porcine dendritic cells infected by the pseudorabies virus

MicroRNAs are small non-coding RNAs approximately 22 nt long that modulate gene expression in animals and plants. It has been recently demonstrated that herpesviruses encode miRNAs to control the post-transcriptional regulation of expression from their own genomes and possibly that of their host, thus adding an additional layer of complexity to the physiological cross-talk between host and pathogen. The present study focussed on the interactions between porcine dendritic cells (DCs) and the Pseudorabies virus (PRV), an alpha-herpesvirus causing Aujeszky's disease in pigs. A catalogue of porcine and viral miRNAs, expressed eight hours post-infection, was established by deep sequencing. An average of 2 million reads per sample with a size of 21–24 nucleotides was obtained from six libraries representing three biological replicates of infected and mock-infected DCs. Almost 95% of reads mapped to the draft pig genome sequence and pig miRNAs previously annotated in dedicated databases were detected by sequence alignment. In silico prediction allowed the identification of unknown porcine as well as of five miRNAs transcribed by the Large Latency Transcript (LLT) of PRV. The gene target prediction of the viral miRNAs and the Ingenuity Pathway Analysis of differentially expressed pig miRNAs were conducted to contextualize the identified small RNA molecules and functionally characterize their involvement in the post-transcriptional regulation of gene expression. The results support a role for PRV miRNAs in the maintenance of the host cell latency state through the down-regulation of immediate-early viral genes which is similar to other herpesviruses. The differentially expressed swine miRNAs identified a unique network of target genes with highly significant functions in the development and function of the nervous system and in infectious mechanisms, suggesting that the modulation of both host and viral miRNAs is necessary for the establishment of PRV latency.

Transcription specificity of the class Ib genes SLA-6, SLA-7 and SLA-8 of the swine major histocompatibility complex and comparison with class Ia genes

Our aim was to analyse the transcription levels of the three non-classical class Ib genes SLA-6, SLA-7 and SLA-8 of the swine major histocompatibility complex in various tissues and conditions and to compare them to the transcription levels of classical class Ia genes. Twenty-five adult tissues from two pig breeds, pig renal PK15 cells infected with the Pseudorabies virus, and peripheral blood mononuclear cells (PBMCs) stimulated by lipopolysaccharide or a mixture of phorbol myristate acetate and ionomycin were included in our study. Relative transcription was quantified by quantitative real-time PCR. On average, in adult tissues and PBMCs and compared to SLA-6, the transcription level of SLA-Ia genes was 100-1000 times higher, the level of SLA-8 was 10-20 times higher, and that of SLA-7 was five times higher. Thus, SLA-8 is the most transcribed SLA-Ib gene, followed by the SLA-7 and SLA-6 genes. The highest transcription levels of SLA-Ib transcripts were found in the lymphoid organs, followed by the lung and the digestive tract. The tissue variability of expression levels was widest for the SLA-6 gene, with a 1:32 ratio between the lowest and highest levels in contrast to a 1:12 ratio for the SLA-7 and SLA-8 genes and a 1:16 ratio for the SLA-Ia genes. During PK-15 infection and PBMC stimulation, SLA-Ia and SLA-8 genes were downregulated, whereas SLA-6 and SLA-7 were upregulated, downregulated or not significantly modified. Our overall results confirm the tissue-wide transcription of the three SLA-Ib genes and suggest that they have complementary roles.

The effects of viral load on pseudorabies virus gene expression

Background

Herpesvirus genes are classified into distinct kinetic groups on the basis of their expression dynamics during lytic growth of the virus in cultured cells at a high, typically 10 plaque-forming units/cell multiplicity of infection (MOI). It has been shown that both the host response and the success of a pathogen are dependent on the quantity of particles infecting an organism. This work is a continuation of an earlier study [1], in which we characterized the overall expression of PRV genes following low-MOI infection. In the present study, we have addressed the question of whether viral gene expressions are dependent on the multiplicity of infection by comparing gene expressions under low and high-MOI conditions.

Results

In the present study, using a real-time RT-PCR assay, we address the question of whether the expression properties of the pseudorabies virus (PRV) genes are dependent on the number of virion particles infecting a single cell in a culture. Our analysis revealed a significant dependence of the gene expression on the MOI in most of these genes. Specifically, we found that most of the examined viral genes were expressed at a lower level at a low MOI (0.1) than at a high MOI (10) experiment in the early stage of infection; however, this trend reversed by six hour post-infection in more than half of the genes. Furthermore, in the high-MOI infection, several PRV genes substantially declined within the 4 to 6-h infection period, which was not the case in the low-MOI infection. In the low-MOI infection, the level of antisense transcript (AST), transcribed from the antiparallel DNA strand of the immediate-early 180 (ie180) gene, was comparable to that of ie180 mRNA, while in the high-MOI experiment (despite the 10 times higher copy number of the viral genome in the infected cells) the amount of AST dropped by more than two log values at the early phase of infection. Furthermore, our analysis suggests that adjacent PRV genes are under a common regulation. This is the first report on the effect of the multiplicity of infection on genome-wide gene expression of large DNA viruses, including herpesviruses.

Conclusion

Our results show a strong dependence of the global expression of PRV genes on the MOI. Furthermore, our data indicate a strong interrelation between the expressions of ie180 mRNA and AST, which determines the expression properties of the herpesvirus genome and possibly the replication strategy (lytic or latent infection) of the virus in certain cell types.

Methods for transcriptomic analyses of the porcine host immune response: application to Salmonella infection using microarrays

Technological developments in both the collection and analysis of molecular genetic data over the past few years have provided new opportunities for an improved understanding of the global response to pathogen exposure. Such developments are particularly dramatic for scientists studying the pig, where tools to measure the expression of tens of thousands of transcripts, as well as unprecedented data on the porcine genome sequence, have combined to expand our abilities to elucidate the porcine immune system. In this review, we describe these recent developments in the context of our work using primarily microarrays to explore gene expression changes during infection of pigs by Salmonella. Thus while the focus is not a comprehensive review of all possible approaches, we provide links and information on both the tools we use as well as alternatives commonly available for transcriptomic data collection and analysis of porcine immune responses. Through this review, we expect readers will gain an appreciation for the necessary steps to plan, conduct, analyze and interpret the data from transcriptomic analyses directly applicable to their research interests.

Transcriptome analysis of porcine PBMCs after in vitro stimulation by LPS or PMA/ionomycin using an expression array targeting the pig immune response

Background

Designing sustainable animal production systems that better balance productivity and resistance to disease is a major concern. In order to address questions related to immunity and resistance to disease in pig, it is necessary to increase knowledge on its immune system and to produce efficient tools dedicated to this species.

Results

A long-oligonucleotide-based chip referred to as SLA-RI/NRSP8-13K was produced by combining a generic set with a newly designed SLA-RI set that targets all annotated loci of the pig major histocompatibility complex (MHC) region (SLA complex) in both orientations as well as immunity genes outside the SLA complex.

The chip was used to study the immune response of pigs following stimulation of porcine peripheral blood mononuclear cells (PBMCs) with lipopolysaccharide (LPS) or a mixture of phorbol myristate acetate (PMA) and ionomycin for 24 hours. Transcriptome analysis revealed that ten times more genes were differentially expressed after PMA/ionomycin stimulation than after LPS stimulation. LPS stimulation induced a general inflammation response with over-expression of SAA1, pro-inflammatory chemokines IL8, CCL2, CXCL5, CXCL3, CXCL2 and CCL8 as well as genes related to oxidative processes (SOD2) and calcium pathways (S100A9 and S100A12). PMA/ionomycin stimulation induced a stronger up-regulation of T cell activation than of B cell activation with dominance toward a Th1 response, including IL2, CD69 and TNFRSF9 (tumor necrosis factor receptor superfamily, member 9) genes. In addition, a very intense repression of THBS1 (thrombospondin 1) was observed. Repression of MHC class I genes was observed after PMA/ionomycin stimulation despite an up-regulation of the gene cascade involved in peptide processing. Repression of MHC class II genes was observed after both stimulations. Our results provide preliminary data suggesting that antisense transcripts mapping to the SLA complex may have a role during immune response.

Conclusion

The SLA-RI/NRSP8-13K chip was found to accurately decipher two distinct immune response activations of PBMCs indicating that it constitutes a valuable tool to further study immunity and resistance to disease in pig. The transcriptome analysis revealed specific and common features of the immune responses depending on the stimulation agent that increase knowledge on pig immunity.

Global transcriptional response of pig brain and lung to natural infection by Pseudorabies virus

Background

Pseudorabies virus (PRV) is an alphaherpesviruses whose native host is pig. PRV infection mainly causes signs of central nervous system disorder in young pigs, and respiratory system diseases in the adult.

Results

In this report, we have analyzed native host (piglets) gene expression changes in response to acute pseudorabies virus infection of the brain and lung using a printed human oligonucleotide gene set from Illumina. A total of 210 and 1130 out of 23,000 transcript probes displayed differential expression respectively in the brain and lung in piglets after PRV infection (p-value < 0.01), with most genes displaying up-regulation. Biological process and pathways analysis showed that most of the up-regulated genes are involved in cell differentiation, neurodegenerative disorders, the nervous system and immune responses in the infected brain whereas apoptosis, cell cycle control, and the mTOR signaling pathway genes were prevalent in the infected lung. Additionally, a number of differentially expressed genes were found to map in or close to quantitative trait loci for resistance/susceptibility to pseudorabies virus in piglets.

Conclusion

This is the first comprehensive analysis of the global transcriptional response of the native host to acute alphaherpesvirus infection. The differentially regulated genes reported here are likely to be of interest for the further study and understanding of host viral gene interactions.

Whole-genome analysis of pseudorabies virus gene expression by real-time quantitative RT-PCR assay

Background

Pseudorabies virus (PRV), a neurotropic herpesvirus of pigs, serves as an excellent model system with which to investigate the herpesvirus life cycle both in cultured cells and in vivo. Real-time RT-PCR is a very sensitive, accurate and reproducible technique that can be used to detect very small amounts of RNA molecules, and it can therefore be applied for analysis of the expression of herpesvirus genes from the very early period of infection.

Results

In this study, we have developed and applied a quantitative reverse transcriptase-based real-time PCR technique in order to profile transcription from the whole genome of PRV after lytic infection in porcine kidney cells. We calculated the relative expression ratios in a novel way, which allowed us to compare different PRV genes with respect to their expression dynamics, and to divide the PRV genes into distinct kinetic classes. This is the first publication on the whole-genome analysis of the gene expression of an alpha-herpesvirus by qRT2-PCR. We additionally established the kinetic properties of uncharacterized PRV genes and revised or confirmed data on PRV genes earlier examined by traditional methods such as Northern blot analysis. Our investigations revealed that genes with the same expression properties form clusters on the PRV genome: nested overlapping genes belong in the same kinetic class, while most convergent genes belong in different kinetic classes. Further, we detected inverse relationships as concerns the expressions of EP0 and IE180 mRNAs and their antisense partners.

Conclusion

Most (if not all) PRV genes begin to be expressed from the onset of viral expression. No sharp boundary was found between the groups of early and late genes classified on the basis of their requirement for viral DNA synthesis. The expressions of the PRV genes were analyzed, categorized and compared by qRT2-PCR assay, with the average of the minimum cycle threshold used as a control for the calculation of a particular R value. In principle, this new calculation technique is applicable for the analysis of gene expression in all temporally changing genetic systems.

Time course differential gene expression in response to porcine circovirus type 2 subclinical infection

This study was aimed at characterizing the potential differences in gene expression in piglets inoculated with Porcine circovirus type 2 (PCV2), the essential causative agent of postweaning multisystemic wasting syndrome. Seven-day-old caesarean-derived, colostrum-deprived piglets were distributed into two groups: control (n = 8) and pigs inoculated with 10^5.2 TCID₅₀ of the Burgos PCV2 isolate (n = 16). One control and three inoculated pigs were necropsied on days 1, 2, 5, and 8 post-infection (p.i.). The remaining pigs (four of each group) were sequentially bled on days 0, 7, 14, 21, and 29 p.i. (necropsy). Total RNA from the mediastinal lymph node (MLN) and lysed whole blood (LWB) samples were hybridized to Affymetrix Porcine GeneChip^®. Forty-three probes were differentially expressed (DE) in MLN samples (FDR < 0.1, fold change > 2) and were distributed into three clusters: globally down-regulated genes, and up-regulated genes at early (first week p.i.) and late (day 29 p.i.) stages of infection. In LWB samples, maximal differences were observed at day 7 p.i., with 54 probes DE between control and inoculated pigs. Main Gene Ontology biological processes assigned to up-regulated genes were related to the immune response. Six common genes were found in both types of samples, all of which belonged to the interferon signaling antiviral effector pathway. Down-regulated genes were mainly related to cell adhesion and migration in MLN, and cellular organization and biogenesis in LWB. Microarray results were validated by quantitative real-time PCR. This study provides, for the first time, the characterization of the early and late molecular events taking place in response to a subclinical PCV2 infection.

Molecular genetics of the swine major histocompatibility complex, the SLA complex

The swine major histocompatibility complex (MHC) or swine leukocyte antigen (SLA) complex is one of the most gene-dense regions in the swine genome. It consists of three major gene clusters, the SLA class I, class III and class II regions, that span approximately 1.1, 0.7 and 0.5Mb, respectively, making the swine MHC the smallest among mammalian MHC so far examined and the only one known to span the centromere. This review summarizes recent updates to the Immuno Polymorphism Database-MHC (IPD-MHC) website (http://www.ebi.ac.uk/ipd/mhc/sla/) which serves as the repository for maintaining a list of all SLA recognized genes and their allelic sequences. It reviews the expression of SLA proteins on cell subsets and their role in antigen presentation and regulating immune responses. It concludes by discussing the role of SLA genes in swine models of transplantation, xenotransplantation, cancer and allergy and in swine production traits and responses to infectious disease and vaccines.

Quantitative whole-cell proteome analysis of pseudorabies virus-infected cells

A quantitative proteome study using the stable isotope labeling with amino acids in cell culture technique was performed on bovine kidney cells after infection with the alphaherpesvirus pseudorabies virus (PrV), the etiological agent of Aujeszky's disease. To enhance yields of proteins to be identified, raw extracts were fractionated by affinity solid-phase extraction with a combination of a cibacron blue F3G-A and a heparin matrix and with a phosphoprotein-specific matrix. After two-dimensional gel electrophoresis in different pH ranges between pH 3 and pH 10, 2,600 proteins representing 565 genes were identified by mass spectrometry and screened for virus-induced changes in relative protein levels. Four hours after infection, significant quantitative variations were found for constituents of the nuclear lamina, representatives of the heterogeneous nuclear ribonucleoproteins, proteins involved in membrane trafficking and intracellular transport, a ribosomal protein, and heat shock protein 27. Several proteins were present in multiple charge variants that were differentially affected by infection with PrV. As a common pattern for all these proteins, a mass shift in favor of the more acidic isoforms was observed, suggesting the involvement of viral or cellular kinases.