The contribution of BvgR, RisA, and RisS to global gene regulation, intracellular cyclic-di-GMP levels, motility, and biofilm formation in Bordetella bronchiseptica.

Bordetella bronchiseptica is a highly contagious respiratory bacterial veterinary pathogen. In this study the contribution of the transcriptional regulators BvgR, RisA, RisS, and the phosphorylation of RisA to global gene regulation, intracellular cyclic-di-GMP levels, motility, and biofilm formation were evaluated. Next Generation Sequencing (RNASeq) was used to differentiate the global gene regulation of both virulence-activated and virulence-repressed genes by each of these factors. The BvgAS system, along with BvgR, RisA, and the phosphorylation of RisA served in cyclic-di-GMP degradation. BvgR and unphosphorylated RisA were found to temporally regulate motility. Additionally, BvgR, RisA, and RisS were found to be required for biofilm formation.

Differential Correlation of Transcriptome Data Reveals Gene Pairs and Pathways Involved in Treatment of Citrobacter rodentium Infection with Bioactive Punicalagin.

This study is part of the work investigating bioactive fruit enzymes as sustainable alternatives to parasite anthelmintics that can help reverse the trend of lost efficacy. The study looked to define biological and molecular interactions that demonstrate the ability of the pomegranate extract punicalagin against intracellular parasites. The study compared transcriptomic reads of two distinct conditions. Condition A was treated with punicalagin (PA) and challenged with Citrobacter rodentium, while condition B (CM) consisted of a group that was challenged and given mock treatment of PBS. To understand the effect of punicalagin on transcriptomic changes between conditions, a differential correlation analysis was conducted. The analysis examined the regulatory connections of genes expressed between different treatment conditions by statistically querying the relationship between correlated gene pairs and modules in differing conditions. The results indicated that punicalagin treatment had strong positive correlations with the over-enriched gene ontology (GO) terms related to oxidoreductase activity and lipid metabolism. However, the GO terms for immune and cytokine responses were strongly correlated with no punicalagin treatment. The results matched previous studies that showed punicalagin to have potent antioxidant and antiparasitic effects when used to treat parasitic infections in mice and livestock. Overall, the results indicated that punicalagin enhanced the effect of tissue-resident genes.

Transcriptomic Analysis of Liver Indicates Novel Vaccine to Porcine Reproductive and Respiratory Virus Promotes Homeostasis in T-Cell and Inflammatory Immune Responses Compared to a Commercial Vaccine in Pigs.

One of the largest impediments for commercial swine production is the presence of Porcine Reproductive and Respiratory Syndrome Virus (PRRSV), a devastating RNA viral infection that is responsible for over $1 billion in loss in the U.S. annually. The challenge with combating PRRSV is a combination of the effect of an extraordinary rate of mutation, the ability to infect macrophages, and subversion of host immune response through a series of actions leading to both immunomodulation and immune evasion. Currently there are a handful of commercial vaccines on the market that have been shown to be effective against homologous infections, but struggle against heterologous or mixed strain infections. However, vaccination is the current best strategy for combating PRRSV, making research into new vaccine technology key. To address these issues with PRRSV and host antiviral functions a novel modified-live vaccine (MLV) able to stimulate known antiviral interferons was created and examined for its ability to potentiate effective immunity and better protection. Here, we examine gene expression in the liver of pigs vaccinated with our novel vaccine, given the liver's large role in antiviral responses and vaccine metabolism. Our study indicated that pigs administered the novel vaccine experience homeostatic gene expression consistent with less inflammation and T-cell depletion risk than pigs administered the commercial vaccine.

Impact of Porcine Arterivirus, Influenza B, and Their Coinfection on Antiviral Response in the Porcine Lung.

Interferon (IFN) cytokines induce autonomous antiviral state in cells of the infected site to restrict virus spreading and critically regulate overall antiviral response. The antiviral state leads to host protection through expression of hundreds of IFN-stimulated genes that restrict viral infection through multiple mechanisms, for example, directly in viral genome degradation and indirectly through cellular metabolic inhibition. Young pigs were split into four treatment groups: control, porcine reproductive and respiratory syndrome virus (PRRSV, also known as porcine arterivirus) infected, influenza B virus (IBV) infected, and IBV/PRRSV coinfection. Lung tissue was collected at 3, 5, and 7 days post infection (dpi) for control, PRRSV and IBV/PRRSV coinfection, and at 3 and 5 dpi for IBV. Transcriptomic analysis, using usegalaxy.org tools, was performed against the S.scrofa 11.1 reference genome. Differentially expressed gene (DEG) analysis was carried out using DeSeq2 based on the model treatment + dpi + treatment:dpi + E. Downstream analysis examined the interaction of DEG at each dpi for over-enriched gene ontology (G.O.) terms and pathways. Comparisons of the infected groups vs. the controls yielded a total of (n = 1412) DEGs for the PRRSV group and (n = 1578) for the IBV/PRRSV group across all timepoints. The IBV group had (n = 64) total DEGs across 3 and 5 dpi. Expression data were considered statistically significant based on false discovery rate (FDR) ⫹ 0.1. Venn diagram comparisons of the DEGs across dpi showed that groups shared only 16 DEGs at 3 dpi, no DEGs were shared at 5 dpi, and for 7 dpi, only the PRRSV and IBV/PRRSV groups were compared and shared a total of 43 DEGs. Across the comparisons, differential expression was observed in antiviral genes such as IRF1, MX1, and OAS2. The IBV and IBV/PRRSV groups showed higher expression of antiviral genes at earlier dpi than the PRRSV group. Additionally, downregulated genes from the comparisons clustered around Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways effecting lung development and cellular integrity. Early expression of host IFN and antiviral genes may lead to viral RNA degradation, and assembly and transcription inhibition in the IBV infections. In comparison, expression of antiviral genes in the PRRSV group decreased across time. The decrease may explain why PRRSV infections persist, while IBV clears. Moreover, all infected groups showed prolonged upregulation in neutrophil degranulation pathway activity, possibly exacerbating symptomatic lung lesion pathology seen in these respiratory infections.

Comparison of the Transcriptome Response within the Swine Tracheobronchial Lymphnode Following Infection with PRRSV, PCV-2 or IAV-S.

Porcine reproductive and respiratory syndrome virus (PRRSV) is a major respiratory pathogen of swine that has become extremely costly to the swine industry worldwide, often causing losses in production and animal life due to their ease of spread. However, the intracellular changes that occur in pigs following viral respiratory infections are still scantily understood for PRRSV, as well as other viral respiratory infections. The aim of this study was to acquire a better understanding of the PRRS disease by comparing gene expression changes that occur in tracheobronchial lymph nodes (TBLN) of pigs infected with either porcine reproductive and respiratory syndrome virus (PRRSV), porcine circovirus type 2 (PCV-2), or swine influenza A virus (IAV-S) infections. The study identified and compared gene expression changes in the TBLN of 80 pigs following infection by PRRSV, PCV-2, IAV-S, or sham inoculation. Total RNA was pooled for each group and time-point (1, 3, 6, and 14 dpi) to make 16 libraries-analyses are by Digital Gene Expression Tag Profiling (DGETP). The data underwent standard filtering to generate a list of sequence tag raw counts that were then analyzed using multidimensional and differential expression statistical tests. The results showed that PRRSV, IAV-S and PCV-2 viral infections followed a clinical course in the pigs typical of experimental infection of young pigs with these viruses. Gene expression results echoed this course, as well as uncovered genes related to intersecting and unique host immune responses to the three viruses. By testing and observing the host response to other respiratory viruses, our study has elucidated similarities and differences that can assist in the development of vaccines and therapeutics that shorten or prevent a chronic PRRSV infection.

Differentially Expressed MiRNAs and tRNA Genes Affect Host Homeostasis During Highly Pathogenic Porcine Reproductive and Respiratory Syndrome Virus Infections in Young Pigs.

Background: Porcine respiratory and reproductive syndrome virus (PRRSV) is a single-stranded RNA virus member that infects pigs and causes losses to the commercial industry reaching upward of a billion dollars annually in combined direct and indirect costs. The virus can be separated into etiologies that contain multiple heterologous low and highly pathogenic strains. Recently, the United States has begun to see an increase in heterologous type 2 PRRSV strains of higher virulence (HP-PRRSV). The high pathogenicity of these strains can drastically alter host immune responses and the ability of the animal to maintain homeostasis. Because the loss of host homeostasis can denote underlying changes in gene and regulatory element expression profiles, the study aimed to examine the effect PRRSV infections has on miRNA and tRNA expression and the roles they play in host tolerance or susceptibility. Results: Using transcriptomic analysis of whole blood taken from control and infected pigs at several time points (1, 3, 8 dpi), the analysis returned a total of 149 statistically significant (FDR ⫹ 0.15) miRNAs (n = 89) and tRNAs (n = 60) that were evaluated for possible pro- and anti-viral effects. The tRNA differential expression increased in both magnitude and count as dpi increased, with no statistically significant expression at 1 dpi, but increases at 3 and 8 dpi. The most abundant tRNA amino acid at 3 dpi was alanine, while glycine was the most abundant at 8 dpi. For the miRNAs, focus was put on upregulation that can inhibit gene expression. These results yielded candidates with potential anti- and pro-viral actions such as Ssc-miR-125b, which is predicted to limit PRRSV viral levels, and Ssc-miR-145-5p shown to cause alternative macrophage priming. The results also showed that both the tRNAs and miRNAs displayed expression patterns. Conclusions: The results indicated that the HP-PRRSV infection affects host homeostasis through changes in miRNA and tRNA expression and their subsequent gene interactions that target and influence the function of host immune, metabolic, and structural pathways.

Natural Selection Footprints Among African Chicken Breeds and Village Ecotypes.

Natural selection is likely a major factor in shaping genomic variation of the African indigenous rural chicken, driving the development of genetic footprints. Selection footprints are expected to be associated with adaptation to locally prevailing environmental stressors, which may include diverse factors as high altitude, disease resistance, poor nutrition, oxidative and heat stresses. To determine the existence of a selection footprint, 268 birds were randomly sampled from three indigenous ecotypes from East Africa (Rwanda and Uganda) and North Africa (Baladi), and two registered Egyptian breeds (Dandarawi and Fayoumi). Samples were genotyped using the chicken Affymetrix 600K Axiom® Array. A total of 494,332 SNPs were utilized in the downstream analysis after implementing quality control measures. The intra-population runs of homozygosity (ROH) that occurred in >50% of individuals of an ecotype or in >75% of a breed were studied. To identify inter-population differentiation due to genetic structure, FST was calculated for North- vs. East-African populations and Baladi and Fayoumi vs. Dandarawi for overlapping windows (500 kb with a step-size of 250 kb). The ROH and FST mapping detected several selective sweeps on different autosomes. Results reflected selection footprints of the environmental stresses, breed behavior, and management. Intra-population ROH of the Egyptian chickens showed selection footprints bearing genes for adaptation to heat, solar radiation, ion transport and immunity. The high-altitude-adapted East-African populations' ROH showed a selection signature with genes for angiogenesis, oxygen-heme binding and transport. The neuroglobin gene (GO:0019825 and GO:0015671) was detected on a Chromosome 5 ROH of Rwanda-Uganda ecotypes. The sodium-dependent noradrenaline transporter, SLC6A2 on a Chromosome 11 ROH in Fayoumi breed may reflect its active behavior. Inter-population FST among Egyptian populations reflected genetic mechanisms for the Fayoumi resistance to Newcastle Disease Virus (NDV), while FST between Egyptian and Rwanda-Uganda populations indicated the Secreted frizzled related protein 2, SFRP2, (GO:0009314) on Chromosome 4, that contributes to melanogenic activity and most likely enhances the Dandarawi chicken adaptation to high-intensity of solar radiation in Southern Egypt. These results enhance our understanding of the natural selection forces role in shaping genomic structure for adaptation to the stressful African conditions.

Identification of Coding and Non-coding RNA Classes Expressed in Swine Whole Blood.

The advent of innovative and increasingly powerful next generation sequencing techniques has opened new avenues into the ability to examine the underlying gene expression related to biological processes of interest. These innovations not only allow researchers to observe expression from the mRNA sequences that code for genes that effect cellular function, but also the non-coding RNA (ncRNA) molecules that remain untranslated, but still have regulatory functions. Although researchers have the ability to observe both mRNA and ncRNA expression, it has been customary for a study to focus on one or the other. However, when studies are interested in both mRNA and ncRNA expression, many times they use separate samples to examine either coding or non-coding RNAs due to the difference in library preparations. This can lead to the need for more samples which can increase time, consumables, and animal stress. Additionally, it may cause researchers to decide to prepare samples for only one analysis, usually the mRNA, limiting the number of biological questions that can be investigated. However, ncRNAs span multiple classes with regulatory roles that effect mRNA expression. Because ncRNA are important to fundamental biologic processes and disorder of these processes in during infection, they may, therefore, make attractive targets for therapeutics. This manuscript demonstrates a modified protocol for the generation mRNA and non-coding RNA expression libraries, including viral RNA, from a single sample of whole blood. Optimization of this protocol, improved RNA purity, increased ligation for recovery of methylated RNAs, and omitted size selection, to allow capture of more RNA species.

Identification of small non-coding RNA classes expressed in swine whole blood during HP-PRRSV infection.

It has been established that reduced susceptibility to porcine reproductive and respiratory syndrome virus (PRRSV) has a genetic component. This genetic component may take the form of small non-coding RNAs (sncRNA), which are molecules that function as regulators of gene expression. Various sncRNAs have emerged as having an important role in the immune system in humans. The study uses transcriptomic read counts to profile the type and quantity of both well and lesser characterized sncRNAs, such as microRNAs and small nucleolar RNAs to identify and quantify the classes of sncRNA expressed in whole blood between healthy and highly pathogenic PRRSV-infected pigs. Our results returned evidence on nine classes of sncRNA, four of which were consistently statistically significantly different based on Fisher's Exact Test, that can be detected and possibly interrogated for their effect on host dysregulation during PRRSV infections.

Comparative analysis of signature genes in PRRSV-infected porcine monocyte-derived cells to different stimuli.

Monocyte-derived DCs (mDCs) are major target cells in porcine reproductive and respiratory syndrome virus (PRRSV) pathogenesis; however, the plasticity of mDCs in response to activation stimuli and PRRSV infection remains unstudied. In this study, we polarized mDCs, and applied genome-wide transcriptomic analysis and predicted protein-protein interaction networks to compare signature genes involved in mDCs activation and response to PRRSV infection. Porcine mDCs were polarized with mediators for 30 hours, then mock-infected, infected with PRRSV strain VR2332, or a highly pathogenic PRRSV strain (rJXwn06), for 5 h. Total RNA was extracted and used to construct sequencing libraries for RNA-Seq. Comparisons were made between each polarized and unpolarized group (i.e. mediator vs. PBS), and between PRRSV-infected and uninfected cells stimulated with the same mediator. Differentially expressed genes (DEG) from the comparisons were used for prediction of interaction networks affected by the viruses and mediators. The results showed that PRRSV infection inhibited M1-prone immune activity, downregulated genes, predicted network interactions related to cellular integrity, and inflammatory signaling in favor of M2 activity. Additionally, the number of DEG and predicted network interactions stimulated in HP-PRRSV infected mDCs was superior to the VR-2332 infected mDCs and conformed with HP-PRRSV pathogenicity.

Genomic Comparison of Indigenous African and Northern European Chickens Reveals Putative Mechanisms of Stress Tolerance Related to Environmental Selection Pressure.

Global climate change is increasing the magnitude of environmental stressors, such as temperature, pathogens, and drought, that limit the survivability and sustainability of livestock production. Poultry production and its expansion is dependent upon robust animals that are able to cope with stressors in multiple environments. Understanding the genetic strategies that indigenous, noncommercial breeds have evolved to survive in their environment could help to elucidate molecular mechanisms underlying biological traits of environmental adaptation. We examined poultry from diverse breeds and climates of Africa and Northern Europe for selection signatures that have allowed them to adapt to their indigenous environments. Selection signatures were studied using a combination of population genomic methods that employed FST , integrated haplotype score (iHS), and runs of homozygosity (ROH) procedures. All the analyses indicated differences in environment as a driver of selective pressure in both groups of populations. The analyses revealed unique differences in the genomic regions under selection pressure from the environment for each population. The African chickens showed stronger selection toward stress signaling and angiogenesis, while the Northern European chickens showed more selection pressure toward processes related to energy homeostasis. The results suggest that chromosomes 2 and 27 are the most diverged between populations and the most selected upon within the African (chromosome 27) and Northern European (chromosome 2) birds. Examination of the divergent populations has provided new insight into genes under possible selection related to tolerance of a population's indigenous environment that may be baselines for examining the genomic contribution to tolerance adaptions.

Leading edge analysis of transcriptomic changes during pseudorabies virus infection.

Eight RNA samples taken from the tracheobronchial lymph nodes (TBLN) of pigs that were either infected or non-infected with a feral isolate of porcine pseudorabies virus (PRV) were used to investigate changes in gene expression related to the pathogen. The RNA was processed into fastq files for each library prior to being analyzed using Illumina Digital Gene Expression Tag Profiling sequences (DGETP) which were used as the downstream measure of differential expression. Analyzed tags consisted of 21 base pair sequences taken from time points 1, 3, 6, and 14 days' post infection (dpi) that generated 1,927,547 unique tag sequences. Tag sequences were analyzed for differential transcript expression and gene set enrichment analysis (GSEA) to uncover transcriptomic changes related to PRV pathology progression. In conjunction with the DGETP and GSEA, the study also incorporated use of leading edge analysis to help link the TBLN transcriptome data to clinical progression of PRV at each of the sampled time points. The purpose of this manuscript is to provide useful background on applying the leading edge analysis to GSEA and expression data to help identify genes considered to be of high biological interest. The data in the form of fastq files has been uploaded to the NCBI Gene Expression Omnibus (GEO) (GSE74473) database.