Transcriptome analysis revealed immune responses in the kidney of Atlantic salmon (Salmo salar) co-infected with sea lice (Lepeophtheirus salmonis) and infectious salmon anemia virus.

Sea lice (Lepeophtheirus salmonis) and infectious salmon anemia virus (ISAv) are two of the most important pathogens in Atlantic salmon (Salmo salar) farming and typically cause substantial economic losses to the industry. However, the immune interactions between hosts and these pathogens are still unclear, especially in the scenario of co-infection. In this study, we artificially infected Atlantic salmon with sea lice and ISAv, and investigated the gene expression patterns of Atlantic salmon head kidneys in response to both lice only and co-infection with lice and ISAv by transcriptomic analysis. The challenge experiment indicated that co-infection resulted in a cumulative mortality rate of 47.8 %, while no mortality was observed in the lice alone infection. We identified 240 differentially expressed genes (DEGs) under the lice alone infection, of which 185 were down-regulated and 55 were up-regulated, while a total of 994 DEGs were identified in the co-infection, of which 206 were down-regulated and 788 were significantly up-regulated. The pathway enrichment analysis revealed that single-infection significantly suppressed the innate immune system (e.g., the complement system), whereas co-infection induced a strong immune response, leading to the activation of immune-related signaling pathways such as Toll-like receptors and NOD-like receptors pathways, as well as significant upregulation of genes related to the activation of interferon and MH class I protein complex. Our results provide the first global transcriptomic study of gene expression in the Atlantic salmon head kidney in response to co-infection with sea lice and ISAv, and provided the baseline knowledge for understanding the immune responses during co-infection.

Transcriptome profiling of blood from common bottlenose dolphins (Tursiops truncatus) in the northern Gulf of Mexico to enhance health assessment capabilities.

Following the 2010 Deepwater Horizon disaster and subsequent unusual mortality event, adverse health impacts have been reported in bottlenose dolphins in Barataria Bay, LA including impaired stress response and reproductive, pulmonary, cardiac, and immune function. These conditions were primarily diagnosed through hands-on veterinary examinations and analysis of standard diagnostic panels. In human and veterinary medicine, gene expression profiling has been used to identify molecular mechanisms underlying toxic responses and disease states. Identification of molecular markers of exposure or disease may enable earlier detection of health effects or allow for health evaluation when the use of specialized methodologies is not feasible. To date this powerful tool has not been applied to augment the veterinary data collected concurrently during dolphin health assessments. This study examined transcriptomic profiles of blood from 76 dolphins sampled in health assessments during 2013-2018 in the waters near Barataria Bay, LA and Sarasota Bay, FL. Gene expression was analyzed in conjunction with the substantial suite of health data collected using principal component analysis, differential expression testing, over-representation analysis, and weighted gene co-expression network analysis. Broadly, transcript profiles of Barataria Bay dolphins indicated a shift in immune response, cytoskeletal alterations, and mitochondrial dysfunction, most pronounced in dolphins likely exposed to Deepwater Horizon oiling. While gene expression profiles in Barataria Bay dolphins were altered compared to Sarasota Bay for all years, profiles from 2013 exhibited the greatest alteration in gene expression. Differentially expressed transcripts included genes involved in immunity, inflammation, reproductive failure, and lung or cardiac dysfunction, all of which have been documented in dolphins from Barataria Bay following the Deepwater Horizon oil spill. The genes and pathways identified in this study may, with additional research and validation, prove useful as molecular markers of exposure or disease to assist wildlife veterinarians in evaluating the health of dolphins and other cetaceans.

A novel full-length transcriptome resource from multiple immune-related tissues in turbot (Scophthalmus maximus) using Pacbio SMART sequencing.

Turbot (Scophthalmus maximus) is an important cold-water economic fish. However, the production and development of turbot industry has been constantly hindered by the frequent occurrence of some diseases. Lacking full-length transcriptome for turbot limits immune gene discoveries and gene structures analysis. Therefore, we generated a full-length transcriptome using mixed immune-related tissues of turbot with PacBio Sequel platform. In this study, a total of 31.7 Gb high quality data were generated with the average subreads length of 2618 bp. According to the presence of 5' and 3' primers as well as poly (A) tails, FL (Full-length) and NFL (Non-full-length) isoforms were obtained. Meanwhile, we identified 32,003 non-redundant transcripts, 76.02% of which was novel isoforms of known genes. In addition, 12,176 alternative splicing (AS) events, 6614 polyadenylation (APA) events, 1905 transcription factors, and 2703 lncRNAs were identified. This work is a comprehensive report on the full-length transcriptome of immune-related tissues of turbot, and it also provides valuable molecular resources for future research on the adaptation mechanisms and functional genomics of turbot.

Transcriptome profiling based on larvae at different time points after hatching provides a core set of gene resource for understanding the immune response mechanisms of the egg-protecting behavior against Vibrio anguillarum infection in Amphioctopus fangsiao.

Mollusks have recently received increasing attention because of their unique immune systems. Mollusks such as Amphioctopus fangsiao are economically important cephalopods, and the effects of their egg-protecting behavior on the larval immune response are unclear. Meanwhile, little research has been done on the resistance response of cephalopod larvae infected with pathogenic bacteria such as Vibrio anguillarum. In this study, V. anguillarum was used to infect the primary hatching A. fangsiao larvae under different egg-protecting behaviors for 24 h, and a total of 7156 differentially expressed genes (DEGs) were identified at four time points after hatching based on transcriptome analysis. GO and KEGG enrichment analyses showed that multiple immune-related GO terms and KEGG signaling pathways were enriched. Protein-protein interaction networks (PPI networks) were used to search functional relationships between immune-related DEGs. Finally, 20 hub genes related to multiple gene functions or involved in multiple signaling pathways were identified, and their accuracy was verified using quantitative RT-PCR. PPI networks were first used to study the effects A. fangsiao larvae after infection with V. anguillarum under different egg-protecting behaviors. The results provide significant genetic resources for exploring invertebrate larval immune processes. The data lays a foundation for further study the immune response mechanisms for invertebrates after infection.

Health monitoring in birds using bio-loggers and whole blood transcriptomics.

Monitoring and early detection of emerging infectious diseases in wild animals is of crucial global importance, yet reliable ways to measure immune status and responses are lacking for animals in the wild. Here we assess the usefulness of bio-loggers for detecting disease outbreaks in free-living birds and confirm detailed responses using leukocyte composition and large-scale transcriptomics. We simulated natural infections by viral and bacterial pathogens in captive mallards (Anas platyrhynchos), an important natural vector for avian influenza virus. We show that body temperature, heart rate and leukocyte composition change reliably during an acute phase immune response. Using genome-wide gene expression profiling of whole blood across time points we confirm that immunostimulants activate pathogen-specific gene regulatory networks. By reporting immune response related changes in physiological and behavioural traits that can be studied in free-ranging populations, we provide baseline information with importance to the global monitoring of zoonotic diseases.

Metatranscriptomic virome assessment of Rhipicephalus microplus from Colombia.

Ticks (Ixodida) are hematophagous ectoparasites that harbor and transmit diverse species of viruses, some of which cause serious diseases with worldwide veterinary and human health consequences. Rhipicephalus microplus is an important cattle tick in Colombia, where it causes significant economic losses. Despite the importance of this tick, its viral profile is unknown. RNA sequencing was used in this study as a surveillance method for virus detection in R. microplus. Most of the viral origin contigs were assigned to two putative viruses: one chuvirus (Wuhan tick virus 2) and one phlebovirus-like (Lihan tick virus). In addition, viral contigs corresponding to two jingmenviruses previously reported in R. microplus from China and Brazil were detected, as well as a novel putative tymovirus, named here as Antioquia tymovirus-like 1 (ATV-like 1). The presence of some of these viruses across numerous regions in the world could have several explanations, including i) a long-term association between those viruses and R. microplus and ii) a consequence of livestock historical trade. Our results shed new light on the virus diversity of this tick species and provide a basis for further studies on the evolutionary history and pathogenic potential of these interesting viruses.

Expression and function assessment of two serpin-type serine protease inhibitors from Haemaphysalis doenitzi.

Serine protease inhibitors (serpins) in ticks are implicated in the modulation of the vertebrate host response to the tick bite. Experimentally, it has been demonstrated that serpins interfere with tick-borne pathogen transmission. However, knowledge on serpins in the tick Haemaphysalis doenitzi is lacking. In this study, the expression of two serpin genes, named HDS1 and HDS2, were assessed in H. doenitzi, and their roles in immune regulation were further investigated. The expression of HDS1 and HDS2 showed no tissue specificity, with maximum expression levels detected in the hemolymph and salivary gland, respectively. Among the developmental stages, the highest expression of HDS1 and HDS2 were detected in larvae and adults, respectively. The recombinant protein rHDS1 displayed obvious inhibitory effects on trypsin and thrombin, whereas rHDS2 clearly inhibited thrombin only. In addition, rHDS1 and rHDS2 showed certain inhibitory activities against bacteria and fungi. The female engorgement body weight, female engorgement rate, and egg hatchability were significantly decreased after injection of double-stranded RNA (dsRNA) of HDS1 gene, whereas no significant effects were observed concerning the feeding period or attachment rate at 24 h after introduction via rabbit ears. When injected with dsRNA of HDS2 gene, no significant effect was observed on the attachment rate at 24 h after introduction into the rabbit ears, but the engorgement body weight and engorgement rate of female ticks were significantly decreased, and no egg hatchment occurred. The above results contribute to better understanding the function of serpins in the development and innate immunity of H. doenitzi.

Taxonomic and functional assessment using metatranscriptomics reveals the effect of Angus cattle on rumen microbial signatures.

A greater understanding of the rumen microbiota and its function may help find new strategies to improve feed efficiency in cattle. This study aimed to investigate whether the cattle breed affects specific ruminal taxonomic microbial groups and functions associated with feed conversion ratio (FCR), using two genetically related Angus breeds as a model. Total RNA was extracted from 24 rumen content samples collected from purebred Black and Red Angus bulls fed the same forage diet and then subjected to metatranscriptomic analysis. Multivariate discriminant analysis (sparse partial least square discriminant analysis (sPLS-DA)) and analysis of composition of microbiomes were conducted to identify microbial signatures characterizing Black and Red Angus cattle. Our analyses revealed relationships among bacterial signatures, host breeds and FCR. Although Black and Red Angus are genetically similar, sPLS-DA detected 25 bacterial species and 10 functions that differentiated the rumen microbial signatures between those two breeds. In Black Angus, we identified bacterial taxa Chitinophaga pinensis, Clostridium stercorarium and microbial functions with large and small subunits ribosomal proteins L16 and S7 exhibiting a higher abundance in the rumen microbiome. In Red Angus, nonetheless, we identified the poorly characterized bacterial taxon Oscillibacter valericigenes with a higher abundance and pathways related to carbohydrate metabolism. Analysis of composition of microbiomes revealed that C. pinensis and C. stercorarium exhibited a higher abundance in Black Angus compared to Red Angus associated with FCR, suggesting that these bacterial species may play a key role in the feed conversion efficiency of forage-fed bulls. This study highlights how the discovery of signatures of bacterial taxa and their functions can be used to harness the full potential of the rumen microbiome in Angus cattle.

Assessment of RNA Stability in Postmortem Tissue from New-Born Lambs.

The recovery of high quality RNA from postmortem tissue is crucial to gene expression analyses. The acquisition of postmortem tissue has inherent time delays and, hence, understanding the temporal variation in the stability of total RNA is imperative. This experiment aimed: ( 1 ) to qualitatively and quantitatively assess the integrity of total RNA derived from a range of new-born ovine tissues (liver, spleen, thyroid, skeletal muscle, ileum, and perirenal adipose tissue) which were stored at ambient temperature until extraction at 0, 3, 6, and 9 h postmortem; and ( 2 ) to analyze the stability of the reference gene(s) and expression of specific target genes in these tissues. Postmortem sampling time resulted in variable reductions in the relative integrity number (RIN) values across the tissues, ranging from 0.9 to 1.8% in liver, spleen, skeletal muscle, and ileum to 5.7-11.1% in the thyroid and perirenal adipose tissues, respectively (P < 0.05). In conclusion, tissues with small reductions in RIN value can exhibit disproportionately large differences in the normalization factor used to calculate the target gene expression. Hence, changes in transcript abundance due to RNA degradation are not always sufficiently buffered through normalization with reference genes. The normalization factor should be presented alongside the RIN value in postmortem tissue studies.

Assessment of reference genes for reliable analysis of gene transcription by RT-qPCR in ovine leukocytes.

With the availability of genetic sequencing data, quantitative reverse transcription PCR (RT-qPCR) is increasingly being used for the quantification of gene transcription across species. Too often there is little regard to the selection of reference genes and the impact that a poor choice has on data interpretation. Indeed, RT-qPCR provides a snapshot of relative gene transcription at a given time-point, and hence is highly dependent on the stability of the transcription of the reference gene(s). Using ovine efferent lymph cells and peripheral blood mono-nuclear cells (PBMCs), the two most frequently used leukocytes in immunological studies, we have compared the stability of transcription of the most commonly used ovine reference genes: YWHAZ, RPL-13A, PGK1, B2M, GAPDH, HPRT, SDHA and ACTB. Using established algorithms for reference gene normalization "geNorm" and "Norm Finder", PGK1, GAPDH and YWHAZ were deemed the most stably transcribed genes for efferent leukocytes and PGK1, YWHAZ and SDHA were optimal in PBMCs. These genes should therefore be considered for accurate and reproducible RT-qPCR data analysis of gene transcription in sheep.

Comprehensive cross production system assessment of the impact of in vitro microenvironment on the expression of messengers and long non-coding RNAs in the bovine blastocyst.

In vitro production (IVP) of cattle embryos over the past two decades has revealed several negative impacts that have been attributed to the artificial microenvironment. Studies on embryos produced in vitro clearly point to aberrant gene expression levels. So far, the causal association between phenotype and measured gene expression has not led to substantial improvement of IVP systems. The aim of this study was to generate a unique dataset composed of microarray-derived relative transcript abundance values for blastocysts produced in ten in vitro systems differing primarily in culture medium formulation. Between-group comparisons determine the level of overall similarity among systems relative to in vivo reference embryos. The use of the dataset to contrast all in vitro treatments with the in vivo blastocysts pointed to a single common gene network. The 'boutique' array contained a panel of novel uncharacterized transcripts that were variably expressed depending on the medium in which the blastocysts were produced. These novel transcripts were differentially expressed in blastocysts even as carryover from conditions encountered 7 days earlier during oocyte maturation. All of the selected novel candidates thus expressed were from intergenic regions. The function of this long non-coding RNA remains unknown but clearly points to an additional level of complexity in early embryo development.