Arun Kommadath

Increasing gene discovery and coverage using RNA-seq of globin RNA reduced porcine blood samples.

BackgroundTranscriptome analysis of porcine whole blood has several applications, which include deciphering genetic mechanisms for host responses to viral infection and vaccination. The abundance of alpha- and beta-globin transcripts in blood, however, impedes the ability to cost-effectively detect transcripts of low abundance. Although protocols exist for reduction of globin transcripts from human and mouse/rat blood, preliminary work demonstrated these are not useful for porcine blood Globin Reduction (GR). Our objectives were to develop a porcine specific GR protocol and to evaluate the GR effects on gene discovery and sequence read coverage in RNA-sequencing (RNA-seq) experiments.ResultsA GR protocol for porcine blood samples was developed using RNase H with antisense oligonucleotides specifically targeting porcine hemoglobin alpha (HBA) and beta (HBB) mRNAs. Whole blood samples (n = 12) collected in Tempus tubes were used for evaluating the efficacy and effects of GR on RNA-seq. The HBA and HBB mRNA transcripts comprised an average of 46.1% of the mapped reads in pre-GR samples, but those reads reduced to an average of 8.9% in post-GR samples. Differential gene expression analysis showed that the expression level of 11,046 genes were increased, whereas 34 genes, excluding HBA and HBB, showed decreased expression after GR (FDR <0.05). An additional 815 genes were detected only in post-GR samples.ConclusionsOur porcine specific GR primers and protocol minimize the number of reads of globin transcripts in whole blood samples and provides increased coverage as well as accuracy and reproducibility of transcriptome analysis. Increased detection of low abundance mRNAs will ensure that studies relying on transcriptome analyses do not miss information that may be vital to the success of the study.

Expansion of ruminant-specific microRNAs shapes target gene expression divergence between ruminant and non-ruminant species

BackgroundUnderstanding how species-specific microRNAs (miRNAs) contribute to species-specific phenotypes is a central topic in biology. This study aimed to elucidate the role of ruminant-specific miRNAs in shaping mRNA expression divergence between ruminant and non-ruminant species.ResultsWe analyzed miRNA and mRNA transcriptomes generated by Illumina sequencing from whole blood samples of cattle and a closely related non-ruminant species, pig. We found evidence of expansion of cattle-specific miRNAs by analyzing miRNA conservation among 57 vertebrate species. The emergence of cattle-specific miRNAs was accompanied by accelerated sequence evolution at their target sites. Further, the target genes of cattle-specific miRNAs show markedly reduced expression compared to their pig and human orthologues. We found that target genes with conserved or non-conserved target sites of cattle-specific miRNAs exhibit reduced expression. One of the significantly enriched KEGG pathway terms for the target genes of the cattle-specific miRNAs is the insulin signalling pathway, raising the possibility that some of these miRNAs may modulate insulin resistance in ruminants.ConclusionsWe provide evidence of rapid miRNA-mediated regulatory evolution in the ruminant lineage. Cattle-specific miRNAs play an important role in shaping gene expression divergence between ruminant and non-ruminant species, by influencing the expression of targets genes through both conserved and cattle-specific target sites.

Genome-wide whole blood microRNAome and transcriptome analyses reveal miRNA-mRNA regulated host response to foodborne pathogen Salmonella infection in swine.

To understand the role of miRNAs in regulating genes involved in host response to bacterial infection and shedding of foodborne pathogens, a systematic profiling of miRNAs and mRNAs from the whole blood of pigs upon Salmonella challenge was performed. A total of 62 miRNAs were differentially expressed post infection (false discovery rate <0.1). An integrative analysis of both the differentially expressed miRNAs and mRNAs using sequence-based miRNA target prediction and negative correlation of miRNA-mRNA profiles helped identify miRNA-mRNA networks that may potentially regulate host response to Salmonella infection. From these networks, miR-214 and miR-331-3p were identified as new candidates potentially associated with Salmonella infection. An miRNA seed sequence analysis suggested that these miRNAs regulate several critical immune-related genes including SLC11A1, PIGE-108A11.3 and VAV2. We showed that challenged pigs had reduced miR-214 expression and increased miR-331-3p expression in the whole blood. Furthermore, the expression of the proposed targets of miR-214 (SLC11A1 and PIGE-108A11.3) increased while that of the proposed target of miR-331-3p (VAV2) decreased following challenge (expression changes confirmed by in vitro assays). Based on these observations, we propose potential roles for miR-214 and miR-331-3p in regulation of immune responses to Salmonella infection.

MicroRNA Buffering and Altered Variance of Gene Expression in Response to Salmonella Infection

One potential role of miRNAs is to buffer variation in gene expression, although conflicting results have been reported. To investigate the buffering role of miRNAs in response to Salmonella infection in pigs, we sequenced miRNA and mRNA in whole blood from 15 pig samples before and after Salmonella challenge. By analyzing inter-individual variation in gene expression patterns, we found that for moderately and lowly expressed genes, putative miRNA targets showed significantly lower expression variance compared with non-miRNA-targets. Expression variance between highly expressed miRNA targets and non-miRNA-targets was not significantly different. Further, miRNA targets demonstrated significantly reduced variance after challenge whereas non-miRNA-targets did not. RNA binding proteins (RBPs) are significantly enriched among the miRNA targets with dramatically reduced variance of expression after Salmonella challenge. Moreover, we found evidence that targets of young (less-conserved) miRNAs showed lower expression variance compared with targets of old (evolutionarily conserved) miRNAs. These findings point to the importance of a buffering effect of miRNAs for relatively lowly expressed genes, and suggest that the reduced expression variation of RBPs may play an important role in response to Salmonella infection.

Genetic architecture of gene expression underlying variation in host response to porcine reproductive and respiratory syndrome virus infection

It has been shown that inter-individual variation in host response to porcine reproductive and respiratory syndrome (PRRS) has a heritable component, yet little is known about the underlying genetic architecture of gene expression in response to PRRS virus (PRRSV) infection. Here, we integrated genome-wide genotype, gene expression, viremia level, and weight gain data to identify genetic polymorphisms that are associated with variation in inter-individual gene expression and response to PRRSV infection in pigs. RNA-seq analysis of peripheral blood samples collected just prior to experimental challenge (day 0) and at 4, 7, 11 and 14 days post infection from 44 pigs revealed 6,430 differentially expressed genes at one or more time points post infection compared to the day 0 baseline. We mapped genetic polymorphisms that were associated with inter-individual differences in expression at each day and found evidence of cis-acting expression quantitative trait loci (cis-eQTL) for 869 expressed genes (qval < 0.05). Associations between cis-eQTL markers and host response phenotypes using 383 pigs suggest that host genotype-dependent differences in expression of GBP5, GBP6, CCHCR1 and CMPK2 affect viremia levels or weight gain in response to PRRSV infection.

Differences in Whole Blood Gene Expression Associated with Infection Time-Course and Extent of Fetal Mortality in a Reproductive Model of Type 2 Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) Infection.

Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) infection of pregnant females causes fetal death and increased piglet mortality, but there is substantial variation in the extent of reproductive pathology between individual dams. This study used RNA-sequencing to characterize the whole blood transcriptional response to type 2 PRRSV in pregnant gilts during the first week of infection (at 0, 2, and 6 days post-inoculation), and attempted to identify gene expression signatures associated with a low or high level of fetal mortality rates (LFM and HFM; n = 8/group) at necropsy, 21 days post-inoculation. The initial response to infection measured at 2 days post-inoculation saw an upregulation of genes involved in innate immunity, such as interferon-stimulated antiviral genes and inflammatory markers, and apoptosis. A concomitant decrease in expression of protein synthesis and T lymphocyte markers was observed. By day 6 the pattern had reversed, with a drop in innate immune signaling and an increase in the expression of genes involved in cell division and T cell signaling. Differentially expressed genes (DEGs) associated with extremes of litter mortality rate were identified at all three time-points. Among the 15 DEGs upregulated in LFM gilts on all three days were several genes involved in platelet function, including integrins ITGA2B and ITGB3, and the chemokine PF4 (CXCL4). LFM gilts exhibited a higher baseline expression of interferon-stimulated and pro-inflammatory genes prior to infection, and of T cell markers two days post-infection, indicative of a more rapid progression of the immune response to PRRSV. This study has increased our knowledge of the early response to PRRSV in the blood of pregnant gilts, and could ultimately lead to the development of a biomarker panel that can be used to predict PRRSV-associated reproductive pathology.