The Genome Analysis Toolkit (GATK) is a popular set of programs for discovering and genotyping variants from next-generation sequencing data. The current GATK recommendation for RNA sequencing (RNA-seq) is to perform variant calling from individual samples, with the drawback that only variable positions are reported. Versions 3.0 and above of GATK offer the possibility of calling DNA variants on cohorts of samples using the HaplotypeCaller algorithm in Genomic Variant Call Format (GVCF) mode. Using this approach, variants are called individually on each sample, generating one GVCF file per sample that lists genotype likelihoods and their genome annotations. In a second step, variants are called from the GVCF files through a joint genotyping analysis. This strategy is more flexible and reduces computational challenges in comparison to the traditional joint discovery workflow. Using a GVCF workflow for mining SNP in RNA-seq data provides substantial advantages, including reporting homozygous genotypes for the reference allele as well as missing data. Taking advantage of RNA-seq data derived from primary macrophages isolated from 50 cows, the GATK joint genotyping method for calling variants on RNA-seq data was validated by comparing this approach to a so-called “per-sample” method. In addition, pair-wise comparisons of the two methods were performed to evaluate their respective sensitivity, precision and accuracy using DNA genotypes from a companion study including the same 50 cows genotyped using either genotyping-by-sequencing or with the Bovine SNP50 Beadchip (imputed to the Bovine high density). Results indicate that both approaches are very close in their capacity of detecting reference variants and that the joint genotyping method is more sensitive than the per-sample method. Given that the joint genotyping method is more flexible and technically easier, we recommend this approach for variant calling in RNA-seq experiments.

Spermiogenesis represents the transition from haploid spermatids to spermatozoa. This process entails an extreme condensation of the nucleus and a loss of nearly all cytoplasmic content. The presence of messenger RNAs in the spermatozoa has previously been shown. Generally, these transcripts are considered to be remnants of spermiogenesis. However, it has recently been proposed that there may exist a function for these sperm-associated RNAs. To address the possibility of a functional role for these transcripts, we sought to investigate and characterize the RNA pool found in bovine spermatozoa. The main goals of this study were to examine RNA integrity and survey the mRNA found in spermatids and spermatozoa. Assessment of mRNAs integrity was performed by three approaches: microelectrophoresis, comparative smearing after global amplification, and PCR amplification of target sequences located either in the 5′ or the 3′ ends, while mRNAs survey was performed by microarray hybridizations. RNA integrity studies in the spermatozoa showed a majority of low molecular size fragments indicating a natural segmentation of the mRNA population. The mRNA survey indicated that the sperm transcriptome harbors a complex mixture of messengers implicated in a wide array of cell functions and representing a large subset of transcripts found in spermatids. Subsequently, such sperm RNA profiling could allow the molecular diagnosis of male gamete quality.

Regarding bull fertility, establishing an association between in vitro findings and field fertility requires a multi-parametric approach that measures the integrity of various structures and dynamic functions, such as motion characteristics, among others. The heterogeneous RNA pattern of spermatozoa could be used in genomic analysis for evaluating both spermatogenesis and fertility potential of semen, mainly because of the static status of the transcriptome of this particular differentiated cell. In a previous study, we determined that some spermatozoal transcripts identified by PCR-based cDNA subtraction are associated with non-return rate, a field fertility index. In the present study, the microarray technology was used in conjunction with differential RNA transcript extraction. We have shown that among these genes, some transcripts are also associated with the motility status of a population of sperm cells fractionated from the same ejaculate. We highlighted a systematic data analysis and validation scheme important for the identification of significant transcripts in this context. With such an approach, we found that transcripts encoding a serine/threonine testis-specific protein kinase (TSSK6) and a metalloproteinase non coding RNA (ADAM5P) are associated with high-motility status (P<0.001), also confirmed by quantitative PCR (P=0.0075). This association was found only when transcripts were extracted using the hot-TRIzol protocol, whereas the cold-TRIzol RNA extract comprised mitochondrial transcripts. These results demonstrate that some transcripts previously identified in association with field fertility are also found associated with in vitro motility provided that a stringent RNA extraction protocol is used.

In functional genomics, transcript measurement is of fundamental importance. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) assays are the most popular technology and depend on the initial molecular step, the reverse transcription (RT). This study provides a complex overview of the influence of elements such as RT systems, amount of background RNA, and transcript abundance on the efficiency of qRT-PCR. Using qRT-PCR, we compared the efficiency of some commonly used RT systems and measured the production of PCR-amplifiable products and the influence of PCR inhibitor contents.

Bovine milk fat composition is responsive to dietary manipulation providing an avenue to modify the content of fatty acids and especially some specific unsaturated fatty acid (USFA) isomers of benefit to human health. MicroRNAs (miRNAs) regulate gene expression but their specific roles in bovine mammary gland lipogenesis are unclear. The objective of this study was to determine the expression pattern of miRNAs following mammary gland adaptation to dietary supplementation with 5 % linseed or safflower oil using next generation RNA-sequencing. Twenty-four Canadian Holstein dairy cows (twelve per treatment) in mid lactation were fed a control diet (total mixed ration of corn:grass silages) for 28 days followed by a treatment period (control diet supplemented with 5 % linseed or safflower oil) of 28 days. Milk samples were collected weekly for fat and individual fatty acid determination. RNA from mammary gland biopsies harvested on day-14 (control period) and on days +7 and +28 (treatment period) from six randomly selected cows per treatment was subjected to small RNA sequencing. Milk fat percentage decreased significantly (P < 0.001) during treatment with the two diets as compared to the control period. The individual saturated fatty acids C4:0, C6:0, C8:0, C14:0 and C16:0 decreased significantly (P < 0.05) while five USFAs (C14:1, C18:1n11t, C20:3n3, C20:5n3 and CLA:t10c12) increased remarkably (P < 0.05) in response to both treatments. Analysis of 361 million sequence reads generated 321 known bovine miRNAs and 176 novel miRNAs. The expression of fourteen and twenty-two miRNAs was affected (P < 0.05) by linseed and safflower oil treatments, respectively. Seven miRNAs including six up-regulated (bta-miR-199c, miR-199a-3p, miR-98, miR-378, miR-148b and miR-21-5p) and one down-regulated (bta-miR-200a) were found to be regulated (P < 0.05) by both treatments, and thus considered core differentially expressed (DE) miRNAs. The gene targets of core DE miRNAs have functions related to gene expression and general cellular metabolism (P < 0.05) and are enriched in four pathways of lipid metabolism (3-phosphoinositide biosynthesis, 3-phosphoinositide degradation, D-myo-inisitol-5-phosphate metabolism and the superpathway of inositol phosphate compounds). Our results suggest that DE miRNAs in this study might be important regulators of bovine mammary lipogenesis and metabolism. The novel miRNAs identified in this study will further enrich the bovine miRNome repertoire and contribute to understanding mammary gland biology.

This study evaluated the effect of the probiotics Pediococcus acidilactici and Saccharomyces cerevisiae boulardii on the intestinal colonization of O149 enterotoxigenic Escherichia coli harbouring the F4 (K88) fimbriae (ETEC F4) and on the expression of ileal cytokines in weaned pigs. At birth, different litters of pigs were randomly assigned to one of the following treatments: 1) control without antibiotics or probiotics (CTRL); 2) reference group in which chlortetracycline and tiamulin were added to weanling feed (ATB); 3) P. acidilactici; 4) S. cerevisiae boulardii; or 5) P. acidilactici + S. cerevisiae boulardii. Probiotics were administered daily (1 × 10⁹ CFU per pig) during the lactation period and after weaning (day 21). At 28 days of age, all pigs were orally challenged with an ETEC F4 strain, and a necropsy was performed 24 h later. Intestinal segments were collected to evaluate bacterial colonization in the small intestine and ileal cytokine expressions. Attachment of ETEC F4 to the intestinal mucosa was significantly reduced in pigs treated with P. acidilactici or S. cerevisiae boulardii in comparison with the ATB group (P = 0.01 and P = 0.03, respectively). In addition, proinflammatory cytokines, such as IL-6, were upregulated in ETEC F4 challenged pigs treated with P. acidilactici alone or in combination with S. cerevisiae boulardii compared with the CTRL group. In conclusion, the administration of P. acidilactici or S. cerevisiae boulardii was effective in reducing ETEC F4 attachment to the ileal mucosa, whereas the presence of P. acidilactici was required to modulate the expression of intestinal inflammatory cytokines in pigs challenged with ETEC F4.