Expression patterns of candidate genes with important functions in animal metabolism can help to identify potential molecular markers for cattle production traits. Reverse transcription followed by polymerase chain reaction is a method for rapid and accurate mRNA quantification. However, for exact comparison of mRNA quantity in various samples or tissues, it is important to choose appropriate reference genes. In cattle, little information is available on the expression stability of housekeeping genes (HKGs). The aim of the present study is to develop a set of reference genes that can be used for normalization of concentrations of mRNAs of genes expressed in the bovine liver, kidney, pituitary and thyroid. The study was performed on 6-, 9-, and 12-month-old bulls of dairy and meat cattle breeds. Six HKGs were investigated: ACTB, GAPDH, HPRTI, SDHA, TBP, and YWHAZ. The most stably expressed potential reference HKGs differed among tissues/organs examined: ACTB, TBP, YWHAZ, GAPDH, HPRTI, and SDHA in the liver; GAPDH and YWHAZ in the kidney; GAPDH and SDHA in the pituitary; and TBP and HPRTI in the thyroid. The results showed that the use of a single gene for normalization may lead to relatively large errors, so it is important to use multiple control genes based on a survey of potential reference genes applied to representative samples from specific experimental conditions.

Examination of bovine pituitary gland transcriptome by strand-specific RNA-seq allows detection of putative single nucleotide polymorphisms (SNPs) within potential candidate genes (CGs) or QTLs regions as well as to understand the genomics variations that contribute to economic trait. Here we report a breed-specific model to successfully perform the detection of SNPs in the pituitary gland of young growing bulls representing Polish Holstein-Friesian (HF), Polish Red, and Hereford breeds at three developmental ages viz., six months, nine months, and twelve months. A total of 18 bovine pituitary gland polyA transcriptome libraries were prepared and sequenced using the Illumina NextSeq 500 platform. Sequenced FastQ databases of all 18 young bulls were submitted to NCBI-SRA database with NCBI-SRA accession numbers SRS1296732. For the investigated young bulls, a total of 113,882,3098 raw paired-end reads with a length of 156 bases were obtained, resulting in an approximately 63 million paired-end reads per library. Breed-wise, a total of 515.38, 215.39, and 408.04 million paired-end reads were obtained for Polish HF, Polish Red, and Hereford breeds, respectively. Burrows-Wheeler Aligner (BWA) read alignments showed 93.04%, 94.39%, and 83.46% of the mapped sequencing reads were properly paired to the Polish HF, Polish Red, and Hereford breeds, respectively. Constructed breed-specific SNP-db of three cattle breeds yielded at 13,775,885 SNPs. On an average 765,326 breed-specific SNPs per young bull were identified. Using two stringent filtering parameters, i.e., a minimum 10 SNP reads per base with an accuracy ≥ 90% and a minimum 10 SNP reads per base with an accuracy = 100%, SNP-db records were trimmed to construct a highly reliable SNP-db. This resulted in a reduction of 95,7% and 96,4% cut-off mark of constructed raw SNP-db. Finally, SNP discoveries using RNA-Seq data were validated by KASP™ SNP genotyping assay. The comprehensive QTLs/CGs analysis of 76 QTLs/CGs with RNA-seq data identified KCNIP4, CCSER1, DPP6, MAP3K5 and GHR CGs with highest SNPs hit loci in all three breeds and developmental ages. However, CAST CG with more than 100 SNPs hits were observed only in Polish HF and Hereford breeds.These findings are important for identification and construction of novel tissue specific SNP-db and breed specific SNP-db dataset by screening of putative SNPs according to QTL db and candidate genes for bovine growth and reproduction traits, one can develop genomic selection strategies for growth and reproductive traits.

MYOG and MYF6 belong to the MyoD gene family. They code for the bHLH transcription factors playing a key role in later stages of myogenesis: differentiation and maturation of myotubes. Three SNPs in porcine MYF6 and two in porcine MYOG were analysed in order to establish associations with chosen carcass quality and growth rate traits in Polish Landrace, Polish Large White and line 990 sows. No statistically significant effect of SNP in the promoter region of the MYF6 gene on its expression measured on mRNA level was found. Associations between the genotype at the MYF6 locus and carcass quality traits appeared to be breed-dependent. The C allele in the case of SNP in the promoter region and GC haplotype in exon 1 were advantageous for right carcass side weight in Polish Landrace sows and disadvantageous for this trait in Polish Large White sows. These gene variants were also the most advantageous for loin and ham weight in sows of line 990. The mutation in exon 1 of the MYOG gene had no statistically significant association with carcass quality traits and the mutation in the 3'-flanking region had the breed-dependent effect as well. These results suggest that SNPs analysed in this study are not causative mutations, but can be considered as markers of some other, still unrevealed genetic polymorphism that influences the physiological processes and phenotypic traits considered in this study.

Effect of selenium (Se) supplementation on the selenoprotein and lipid metabolism gene expression patterns in ruminants, especially in lambs is not yet fully understood. The aim of study was to evaluate the effect of Se supplementation on the messenger RNA (mRNA) expression patterns of selected selenoproteins and genes related to lipid metabolism in growing lambs. The experiment was conducted on 48 Polish Merino lambs divided into two groups (n = 24): control (C)—lambs fed with a basal diet (BD) with no Se supplementation, and supplemented (S)—lambs fed with a BD, supplemented with 0.5 mg Se/kg as sodium selenate for 8 weeks. Expression of 12 selenoproteins and six genes related to lipid metabolism was analyzed in the liver and longissimus dorsi (LD) muscle of growing lambs by qPCR. Significant differences were found in the expression of GPX1, GPX2, SEPM, SEPW1, SEP15, SEPGS2, and TXNRD1 in the liver, and GPX1, SEPP1, SEPN1, SEPW1, SEP15, and MSRB1 in the LD muscle between S and C lambs. Se supplementation mainly upregulated SEPW1, SEP15 (P < 0.001; P < 0.01) mRNA expression in the liver, and GPX1, SEPP1, SEPN1, SEPW1 (P < 0.001; P < 0.01) in the muscle of S group. On the other hand, significant decrease in GPX2 (P < 0.01), SEPM (P < 0.001), and SEPHS2 (P < 0.01) mRNA expression levels were observed in the liver of S group of lambs. Se supplementation did not affect PON1, LXRα, and PPARα mRNA expression levels, but a significant increase in mRNA levels of APOE and LPL in the LD muscle (P < 0.05) as well as LPL (P < 0.05) in the liver were noticed in the group of Se supplemented lambs. Our study confirmed that, in lambs, similarly to other species, mRNA expression patterns of several selenoproteins highly depend on dietary Se levels, and their expression is ruled by hierarchical principles and tissue-specific mechanisms. Moreover, the study showed that changes Se intake leads to different levels of genes expression related with lipid metabolism.

The optimal ratio of omega-6 to omega-3 polyunsaturated fatty acids (PUFAs) is important for keeping the homeostasis of biological processes and metabolism, yet the underlying biological mechanism is poorly understood. The objective of this study was to identify changes in the pig liver transcriptome induced by a diet enriched with omega-6 and omega-3 fatty acids and to characterize the biological mechanisms related to PUFA metabolism.Polish Landrace pigs (n = 12) were fed diet enriched with linoleic acid (LA, omega-6) and α-linolenic acid (ALA, omega-3) or standard diet as a control. The fatty acid profiling was assayed in order to verify how feeding influenced the fatty acid content in the liver, and subsequently next-generation sequencing (NGS) was used to identify differentially expressed genes (DEG) between transcriptomes between dietary groups. The biological mechanisms and pathway interaction networks were identified using DAVID and Cytoscape tools. Fatty acid profile analysis indicated a higher contribution of PUFAs in the liver for LA- and ALA-enriched diet group, particularly for the omega-3 fatty acid family, but not omega-6. Next-generation sequencing identified 3565 DEG, 1484 of which were induced and 2081 were suppressed by PUFA supplementation. A low ratio of omega-6/omega-3 fatty acids resulted in the modulation of fatty acid metabolism pathways and over-representation of genes involved in energy metabolism, signal transduction, and immune response pathways.In conclusion, a diet enriched with omega-6 and omega-3 fatty acids altered the transcriptomic profile of the pig liver and would influence animal health status.Electronic supplementary materialThe online version of this article (doi:10.1186/s12263-016-0517-4) contains supplementary material, which is available to authorized users.

Background: The level of omega-6 and omega-3 polyunsaturated fatty acids can affect many cellular systems and function via nuclear receptors or the bioactive lipid regulation of gene expression. The objective of this study was to investigate changes in the muscle transcriptome and the biological functions regulated by increased consumption of omega-3 and omega-6 fatty acids in the pig gluteus medius muscle. Results: The transcriptome of the gluteus medius muscle was studied for pigs subjected to either a control diet or a diet supplemented with linseed and rapeseed oil to increase polyunsaturated fatty acid content. Next-generation sequencing (NGS) was used to generate the muscle tissue transcriptome database pointing differentially expressed genes (DEG). Comparative expression analyses identified 749 genes significantly differing at least in the twofold of change between two groups of animals fed with divergent level of omega-3 and omega-6 fatty acids. The expression of 219 genes was upregulated, and the expression of 530 genes was downregulated in the group of pigs supplemented with omega-3 and omega-6 fatty acids in relation to control group pigs. Results of RNA-seq indicated a role of fatty acid in the regulation of the expression of genes which are essential for muscle tissue development and functioning. Functional analysis revealed that the identified genes were important for a number of biological processes including inflammatory response, signaling, lipid metabolism, and homeostasis. Conclusions: Summarizing, obtained results provide strong evidence that omega-6 and omega-3 fatty acids regulate fundamental metabolic processes in muscle tissue development and functioning.

The objective of this study was to determine hepatic expression levels of GHR, IGF1R, IGF1 and IGF2 genes in young growing gilts at different developmental ages (60-210 days) in five pig breeds: Polish Large White (PLW), Polish Landrace (PL), Pulawska (Pul), Duroc (Dur) and Pietrain (Pie). We studied the differences among pig breeds as well as within each breed for pigs in different developmental ages. Obtained results revealed major differences among breeds in hepatic gene expression of porcine GHR, IGF1R, IGF1 and IGF2 genes in different developmental ages. The differences among breeds of GHR expression were significantly higher in PLW, PL at the age of 60, 90, 120 days as compared to Pul, Dur and Pie. In turn, the highest level of IGF1R expression was observed in PL at age of 150, 180 and 210 days, whereas in case of IGF1 the highest level was recorded in Pie gilts at the age of 60 and 90 days. Moreover trait associated study revealed highly significant correlations between hepatic expressions of IGF1R and IGF2 genes and carcass composition traits (P < 0.01) The results of study suggest that porcine GHR, IGF1R, IGF1 and IGF2 genes may be potential candidate genes for postnatal growth and carcass composition traits. Therefore, the implementation of the hepatic expression of GH/IGF genes into the pig breeding and gene assisted selection program in different pig breeds should be considered. However, further population wide study is needed to clarify the hepatic expression association with economic traits, such as body growth, meat quality and carcass composition traits.

BackgroundRNA-seq is a useful next-generation sequencing (NGS) technology that has been widely used to understand mammalian transcriptome architecture and function. In this study, a breed-specific RNA-seq experiment was utilized to detect putative single nucleotide polymorphisms (SNPs) in liver tissue of young bulls of the Polish Red, Polish Holstein-Friesian (HF) and Hereford breeds, and to understand the genomic variation in the three cattle breeds that may reflect differences in production traits.ResultsThe RNA-seq experiment on bovine liver produced 107,114,4072 raw paired-end reads, with an average of approximately 60 million paired-end reads per library. Breed-wise, a total of 345.06, 290.04 and 436.03 million paired-end reads were obtained from the Polish Red, Polish HF, and Hereford breeds, respectively. Burrows-Wheeler Aligner (BWA) read alignments showed that 81.35%, 82.81% and 84.21% of the mapped sequencing reads were properly paired to the Polish Red, Polish HF, and Hereford breeds, respectively. This study identified 5,641,401 SNPs and insertion and deletion (indel) positions expressed in the bovine liver with an average of 313,411 SNPs and indel per young bull. Following the removal of the indel mutations, a total of 195,3804, 152,7120 and 205,3184 raw SNPs expressed in bovine liver were identified for the Polish Red, Polish HF, and Hereford breeds, respectively. Breed-wise, three highly reliable breed-specific SNP-databases (SNP-dbs) with 31,562, 24,945 and 28,194 SNP records were constructed for the Polish Red, Polish HF, and Hereford breeds, respectively. Using a combination of stringent parameters of a minimum depth of ≥10 mapping reads that support the polymorphic nucleotide base and 100% SNP ratio, 4,368, 3,780 and 3,800 SNP records were detected in the Polish Red, Polish HF, and Hereford breeds, respectively. The SNP detections using RNA-seq data were successfully validated by kompetitive allele-specific PCR (KASPTM) SNP genotyping assay. The comprehensive QTL/CG analysis of 110 QTL/CG with RNA-seq data identified 20 monomorphic SNP hit loci (CARTPT, GAD1, GDF5, GHRH, GHRL, GRB10, IGFBPL1, IGFL1, LEP, LHX4, MC4R, MSTN, NKAIN1, PLAG1, POU1F1, SDR16C5, SH2B2, TOX, UCP3 and WNT10B) in all three cattle breeds. However, six SNP loci (CCSER1, GHR, KCNIP4, MTSS1, EGFR and NSMCE2) were identified as highly polymorphic among the cattle breeds.ConclusionsThis study identified breed-specific SNPs with greater SNP ratio and excellent mapping coverage, as well as monomorphic and highly polymorphic putative SNP loci within QTL/CGs of bovine liver tissue. A breed-specific SNP-db constructed for bovine liver yielded nearly six million SNPs. In addition, a KASPTM SNP genotyping assay, as a reliable cost-effective method, successfully validated the breed-specific putative SNPs originating from the RNA-seq experiments.