Characterization and comparison of the leukocyte transcriptomes of three cattle breeds

Differential gene expression in peripheral leukocytes of pre-weaned Holstein heifer calves with respiratory disease

Bovine respiratory disease (BRD) is a leading cause of calf morbidity and mortality, and prevalence remains high despite current management practices. Differential gene expression (DGE) provides detailed insight into individual immune responses and can illuminate enriched pathways and biomarkers that contribute to disease susceptibility and outcomes. The aims of this study were to investigate differences in peripheral leukocyte gene expression in Holstein preweaned heifer calves 1) with and without BRD, and 2) across weeks of age. Calves were enrolled for this short-term longitudinal study on two commercial dairies in Washington State. Calves were assessed every two weeks throughout the pre-weaning period using clinical respiratory scoring (CRS) and thoracic ultrasonography (TUS), and blood samples were collected. Calves were selected that were either healthy (n = 10) or had BRD diagnosed by CRS (n = 7), TUS (n = 6), or both (n = 6) in weeks 5 or 7 of life). Three consecutive time point samples were analyzed for each BRD calf consisting of PRE, ONSET, and POST samples. Nineteen genes of interest were selected based on previous gene expression studies in cattle: ALOX15, BPI, CATHL6, CXCL8, DHX58, GZMB, HPGD, IFNG, IL17D, IL1R2, ISG15, LCN2, LIF, MX1, OAS2, PGLYRP1, S100A8, SELP, and TNF. Comparisons were made between age and disease time point matched BRD and healthy calves as well as between calf weeks of age. No DGE was observed between diseased and healthy calves; however, DGE was observed between calf weeks of age regardless of disease state. Developmental differences in leukocyte gene expression, phenotype, and functionality make pre-weaned calves immunologically distinct from mature cattle, and early life shifts in calf leukocyte populations likely contribute to the age-related gene expression differences we observed. Age overshadows disease impacts to influence gene expression in young calves, and immune development progresses upon a common trajectory regardless of disease during the preweaning period.

Circular RNA Expression and Regulation Profiling in Testicular Tissues of Immature and Mature Wandong Cattle (Bos taurus)

Wandong cattle are an autochthonous Chinese breed used extensively for beef production. The breed tolerates extreme weather conditions and raw feed and is resistant to tick-borne diseases. However, the genetic basis of testis development and sperm production as well as breeding management is not well established in local cattle. Therefore, improving the reproductive efficiency of bulls via genetic selection is crucial as a single bull can breed thousands of cows through artificial insemination (AI). Testis development and spermatogenesis are regulated by hundreds of genes and transcriptomes. However, circular RNAs (circRNAs) are the key players in many biological developmental processes that have not been methodically described and compared between immature and mature stages in Bovine testes. In this study, we performed total RNA-seq and comprehensively analyzed the circRNA expression profiling of the testis samples of six bulls at 3 years and 3 months of developmental age. In total, 17,013 circRNAs were identified, of which 681 circRNAs (p-adjust < 0.05) were differentially expressed (DE). Among these DE circRNAs, 579 were upregulated and 103 were downregulated in calf and bull testes. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses revealed that the identified target genes were classified into three broad functional categories, including biological process, cellular component, and molecular function, and were enriched in the lysine degradation, cell cycle, and cell adhesion molecule pathways. The binding interactions between DE circRNAs and microRNAs (miRNAs) were subsequently constructed using bioinformatics approaches. The source genes ATM, CCNA1, GSK3B, KMT2C, KMT2E, NSD2, SUCLG2, QKI, HOMER1, and SNAP91 were found to be actively associated with bull sexual maturity and spermatogenesis. In addition, a real-time quantitative polymerase chain reaction (RT-qPCR) analysis showed a strong correlation with the sequencing data. Moreover, the developed model of Bovine testes in the current study provides a suitable framework for understanding the mechanism of circRNAs in the development of testes and spermatogenesis.

Allele Specific Expression (ASE) analysis between Bos Taurus and Bos Indicus cows using RNA-Seq data at SNP level and gene level

In the current study, allele specific expression analysis was performed in two subspecies cows (Bos taurus and Bos indicus) at SNP and gene levels. RNA-Seq data of 21,078,477 and 20940063 paired end reads from pooling of whole blood samples (Leukocyte) from 40 US Holstein (Bos Taurus) and 45 Cholistani cows (Bos indicus) obtained from SRA database in NCBI. Quality control and trimming of row RNA-Seq data were processed by FASTQC and Trimmomatic softwares. The transcriptome was assembled by TopHat2 software in two cow's population by aligning and mapping the RNA-Seq reads on bovine reference genome. The SNPs were discovered by Samtools software and ASE analysis was performed by Chi-square test. Results showed that 50183 and 137954 SNPs were discovered on the assembled transcriptome of Holstein and Cholistani cow samples, respectively, and 15308 SNPs were common in both breeds. 10158 SNPs from 50183 (20%) in Holstein and 31523 SNPs from 137954 (23%) in Cholistani cows were identified as ASE-SNPs. Reference allele and alternative allele count in Holstein and Cholistani cows were 3041 and 7155, respectively. Among 131 discovered SNPs in 41 genes with different expression in Holstein and Cholistani cows, 31 ASE-SNPs (5 in Holstein; 26 in Cholistani cows) were discovered.

Common and Differential Dynamics of the Function of Peripheral Blood Mononuclear Cells between Holstein and Jersey Cows in Heat-Stress Environment

Simple Summary

Seasonal change, particularly changing to hot and humid season, has a negative effect on dairy cows in various ways, including productivity, reproduction, metabolism, and immunity. In high-temperature and humid weather, dairy cows are vulnerable to diseases by weakened immune system. However, the cause of this has not been fully described. Therefore, this study aims to understand changes of specific gene expression and immune pathways based on transcriptome analysis from peripheral blood mononuclear cells of Holstein and Jersey dairy cows between normal and heat-stress environmental conditions. We observed that the two breeds of dairy cow have common and different immune shifts according to the changes of temperature and humidity condition. Overall, the findings of this study improve the understanding of the underlying mechanisms by which seasonal changes affect dairy cow immunity.

Abstract

Heat stress has been reported to affect the immunity of dairy cows. However, the mechanisms through which this occurs are not fully understood. Two breeds of dairy cow, Holstein and Jersey, have distinct characteristics, including productivity, heat resistance, and disease in high-temperature environments. The objective of this study is to understand the dynamics of the immune response of two breeds of dairy cow to environmental change. Ribonucleic acid sequencing (RNA-seq) results were analyzed to characterize the gene expression change of peripheral blood mononuclear cells (PBMCs) in Holstein and Jersey cows between moderate temperature-humidity index (THI) and high THI environmental conditions. Many of the differentially expressed genes (DEGs) identified are associated with critical immunological functions, particularly phagocytosis, chemokines, and cytokine response. Among the DEGs, CXCL3 and IL1A were the top down-regulated genes in both breeds of dairy cow, and many DEGs were related to antimicrobial immunity. Functional analysis revealed that cytokine and chemokine response-associated pathways in both Holstein and Jersey PBMCs were the most important pathways affected by the THI environmental condition. However, there were also breed-specific genes and pathways that altered according to THI environmental condition. Collectively, there were both common and breed-specific altered genes and pathways in Holstein and Jersey cows. The findings of this study expand our understanding of the dynamics of immunity in different breeds of dairy cow between moderate THI and high THI environmental conditions.

Identification of candidate genes for milk production traits by RNA sequencing on bovine liver at different lactation stages

Background

RNA-sequencing was performed to explore the bovine liver transcriptomes of Holstein cows to detect potential functional genes related to lactation and milk composition traits in dairy cattle. The bovine transcriptomes of the nine liver samples from three Holstein cows during dry period (50-d prepartum), early lactation (10-d postpartum), and peak of lactation (60-d postpartum) were sequenced using the Illumina HiSeq 2500 platform.

Results

A total of 204, 147 and 81 differentially expressed genes (DEGs, p < 0.05, false discovery rate q < 0.05) were detected in early lactation vs. dry period, peak of lactation vs. dry period, and peak of lactation vs. early lactation comparison groups, respectively. Gene ontology and KEGG pathway analysis showed that these DEGs were significantly enriched in specific biological processes related to metabolic and biosynthetic and signaling pathways of PPAR, AMPK and p53 (p < 0.05). Ten genes were identified as promising candidates affecting milk yield, milk protein and fat traits in dairy cattle by using an integrated analysis of differential gene expression, previously reported quantitative trait loci (QTL), data from genome-wide association studies (GWAS), and biological function information. These genes were APOC2, PPP1R3B, PKLR, ODC1, DUSP1, LMNA, GALE, ANGPTL4, LPIN1 and CDKN1A.

Conclusion

This study explored the complexity of the liver transcriptome across three lactation periods in dairy cattle by performing RNA sequencing. Integrated analysis of DEGs and reported QTL and GWAS data allowed us to find ten key candidate genes influencing milk production traits.

Systematic identification of long intergenic non-coding RNAs expressed in bovine oocytes

BACKGROUND

Long non-coding RNAs (lncRNAs) are key regulators of diverse cellular processes. Although a number of studies have reported the identification of bovine lncRNAs across many tissues, very little is known about the identity and characteristics of lncRNAs in bovine oocytes.

METHODS

A bovine oocyte cDNA library was constructed and sequenced using the Illumina HiSeq 2000 sequencing system. The oocyte transcriptome was constructed using the ab initio assembly software Scripture and Cufflinks. The assembled transcripts were categorized to identify the novel intergenic transcripts, and the coding potential of these novel transcripts was assessed using CPAT and PhyloCSF. The resulting candidate long intergenic non-coding RNAs (lincRNAs) transcripts were further evaluated to determine if any of them contain any known protein coding domains in the Pfam database. RT-PCR was used to analyze the expression of oocyte-expressed lincRNAs in various bovine tissues.

RESULTS

A total of 85 million raw reads were generated from sequencing of the bovine oocyte library. Transcriptome reconstruction resulted in the assembly of a total of 42,396 transcripts from 37,678 genomic loci. Analysis of the assembled transcripts using the step-wide pipeline resulted in the identification of 1535 oocyte lincRNAs corresponding to 1183 putative non-coding genes. A comparison of the oocyte lincRNAs with the lncRNAs reported in other bovine tissues indicated that 970 of the 1535 oocyte lincRNAs appear to be unique to bovine oocytes. RT-PCR analysis of 5 selected lincRNAs showed either specific or predominant expression of 4 lincRNAs in the fetal ovary. Functional prediction of the oocyte-expressed lincRNAs suggested their involvement in oogenesis through regulating their neighboring protein-coding genes.

CONCLUSIONS

This study provides a starting point for future research aimed at understanding the roles of lncRNAs in controlling oocyte development and early embryogenesis in cattle.

Fasciola hepatica Infection in Cattle: Analyzing Responses of Peripheral Blood Mononuclear Cells (PBMC) Using a Transcriptomics Approach

The parasitic helminth Fasciola hepatica (liver fluke) causes economic loss to the livestock industry globally and also causes zoonotic disease. New control strategies such as vaccines are urgently needed, due to the rise of drug resistance in parasite populations. Vaccine development requires a comprehensive understanding of the immunological events during infection. Previous in vivo studies by our group have investigated global differentially expressed genes (DEGs) in ovine peripheral blood mononuclear cells (PBMC) in response to both acute and chronic F. hepatica infection. This work demonstrated that pathways involved in the pathogenesis of ovine fasciolosis included fibrosis, inhibition of macrophage nitric oxide production, and antibody isotype switching, among others. Transcriptomic changes in PBMC populations following F. hepatica infection in cattle, in which the disease phenotype is quite different, have not yet been examined. Using RNA sequencing we investigated gene expression changes in PBMC isolated from 9 non-infected and 11 F. hepatica-experimentally-infected calves immediately before infection, at 1 and at 14 weeks post-infection. Longitudinal time-course comparisons between groups revealed 21 and 1,624 DEGs driven exclusively by F. hepatica infection in cattle at acute and chronic stages, respectively. These results show that fewer DEGs at the acute stage of infection can be identified in cattle, as compared with sheep. In addition, the log₂ fold-changes of these DEGs were relatively low (−1 to 3) reflecting the different clinical presentation of F. hepatica infection in cattle. Gene pathways for hepatic fibrosis and hepatic cholestasis along with apoptosis of antigen-presenting cells were enriched at chronic stages. Our results reflect the major differences in the disease phenotype between cattle and sheep and may indicate pathways to target in vaccine development.

Comparative Analysis of the Liver Transcriptome among Cattle Breeds Using RNA-seq

Global gene expression in liver transcriptome varies among cattle breeds. The present investigation was aimed to identify the differentially expressed genes (DEGs), metabolic gene networks and metabolic pathways in bovine liver transcriptome of young bulls. In this study, we comparatively analyzed the bovine liver transcriptome of dairy (Polish Holstein Friesian (HF); n = 6), beef (Hereford; n = 6), and dual purpose (Polish-Red; n = 6) cattle breeds. This study identified 895, 338, and 571 significant (p < 0.01) differentially expressed (DE) gene-transcripts represented as 745, 265, and 498 hepatic DE genes through the Polish-Red versus Hereford, Polish-HF versus Hereford, and Polish-HF versus Polish-Red breeds comparisons, respectively. By combining all breeds comparisons, 75 hepatic DE genes (p < 0.01) were identified as commonly shared among all the three breed comparisons; 70, 160, and 38 hepatic DE genes were commonly shared between the following comparisons: (i) Polish-Red versus Hereford and Polish-HF versus Hereford; (ii) Polish-Red versus Hereford and Polish-HF versus Polish-Red; and (iii) Polish-HF versus Hereford and Polish-HF versus Polish-Red, respectively. A total of 440, 82, and 225 hepatic DE genes were uniquely observed for the Polish-Red versus Hereford, Polish-HF versus Hereford, and Polish-Red versus Polish-HF comparisons, respectively. Gene ontology (GO) analysis identified top-ranked enriched GO terms (p < 0.01) including 17, 16, and 31 functional groups and 151, 61, and 140 gene functions that were DE in all three breed liver transcriptome comparisons. Gene network analysis identified several potential metabolic pathways involved in glutamine family amino-acid, triglyceride synthesis, gluconeogenesis, p38MAPK cascade regulation, cholesterol biosynthesis (Polish-Red versus Hereford); IGF-receptor signaling, catecholamine transport, lipoprotein lipase, tyrosine kinase binding receptor (Polish-HF versus Hereford), and PGF-receptor binding, (Polish-HF versus Polish-Red). Validation results showed that the relative expression values were consistent to those obtained by RNA-seq, and significantly correlated between the quantitative reverse transcription PCR (RT-qPCR) and RNA-seq (Pearson’s r  > 0.90). Our results provide new insights on bovine liver gene expressions among dairy versus dual versus beef breeds by identifying the large numbers of DEGs markers submitted to NCBI gene expression omnibus (GEO) accession number GSE114233, which can serve as useful genetic tools to develop the gene assays for trait-associated studies as well as, to effectively implement in genomics selection (GS) cattle breeding programs in Poland.

Transcriptome analysis of adipose tissues from two fat-tailed sheep breeds reveals key genes involved in fat deposition

BACKGROUND

The level of fat deposition in carcass is a crucial factor influencing meat quality. Guangling Large-Tailed (GLT) and Small-Tailed Han (STH) sheep are important local Chinese fat-tailed breeds that show distinct patterns of fat depots. To gain a better understanding of fat deposition, transcriptome profiles were determined by RNA-sequencing of perirenal, subcutaneous, and tail fat tissues from both the sheep breeds. The common highly expressed genes (co-genes) in all the six tissues, and the genes that were differentially expressed (DE genes) between these two breeds in the corresponding tissues were analyzed.

RESULTS

Approximately 47 million clean reads were obtained for each sample, and a total of 17,267 genes were annotated. Of the 47 highly expressed co-genes, FABP4, ADIPOQ, FABP5, and CD36 were the four most highly transcribed genes among all the known genes related to adipose deposition. FHC, FHC-pseudogene, and ZC3H10 were also highly expressed genes and could, thus, have roles in fat deposition. A total of 2091, 4233, and 4131 DE genes were identified in the perirenal, subcutaneous, and tail fat tissues between the GLT and STH breeds, respectively. Gene Ontology (GO) analysis showed that some DE genes were associated with adipose metabolism. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that PPAR signaling pathway and ECM-receptor interaction were specifically enriched. Four genes, namely LOC101102230, PLTP, C1QTNF7, and OLR1 were up-regulated and two genes, SCD and UCP-1, were down-regulated in all the tested tissues of STH. Among the genes involved in ECM-receptor interaction, the genes encoding collagens, laminins, and integrins were quite different depending on the depots or the breeds. In STH, genes such as LAMB3, RELN, TNXB, and ITGA8, were identified to be up regulated and LAMB4 was observed to be down regulated.

CONCLUSIONS

This study unravels the complex transcriptome profiles in sheep fat tissues, highlighting the candidate genes involved in fat deposition. Further studies are needed to investigate the roles of the candidate genes in fat deposition and in determining the meat quality of sheep.

Performance Measurement and Comparative Transcriptome Analysis Revealed the Efforts on Hybrid Improvement of Qinchuan Cattle

Crossbreeding can provide productive gains through heterosis, however, surveys about the effects of crossbreeding through global transcriptomic sequencing are few. This study revealed that Angus × Qinchuan cattle (AQF) have improved performance characteristics compared to Qinchuan cattle (QCF). We performed RNA-seq on the subcutaneous fat tissue of QCF and AQF. More than 42.2 million clean reads were obtained in each sample. We detected 40 and 21 breed-specific highly expressed genes (FPKM > 500) in QCF and AQF, respectively. Furthermore, a total of 353 differentially expressed genes (DEGs, |log2 ratio| ≥ 1 and Probability ≥ 0.8) were found between these two groups, of which 227 genes were upregulated in AQF and 126 genes were upregulated in QCF. Functional enrichment analyses showed that breed-specific highly expressed genes and DEGs were closely related to terms such as development in AQF, and adaption or immune in QCF. In addition, we also identified the novel transcript units, alternative splicing events, single-nucleotide polymorphisms and Indels. Our results revealed differences in inherent characteristics and genetic differences when comparing QCF with AQF.

Transcriptome analysis of the pectoral muscles of local chickens and commercial broilers using Ribo-Zero ribonucleic acid sequencing

BACKGROUND

The molecular mechanisms underlying meat quality and muscle growth are not clear. The meat quality and growth rates of local chickens and commercial broilers are very different. The Ribo-Zero RNA-Seq technology is an effective means of analyzing transcript groups to clarify molecular mechanisms. The aim of this study was to provide a reference for studies of the differences in the meat quality and growth of different breeds of chickens.

RESULTS

Ribo-Zero RNA-Seq technology was used to analyze the pectoral muscle transcriptomes of Gushi chickens and AA broilers. Compared with AA broilers, 1649 genes with annotated information were significantly differentially expressed (736 upregulated and 913 downregulated) in Gushi chickens with Q≤0.05 (Q is the P-value corrected by multiple assumptions test) at a fold change ≥2 or ≤0.5. In addition, 2540 novel significantly differentially expressed (SDE) genes (1405 upregulated and 1135 downregulated) were discovered. The results showed that the main signal transduction pathways that differed between Gushi chickens and AA broilers were related to amino acid metabolism. Amino acids are important for protein synthesis, and they regulate key metabolic pathways to improve the growth, development and reproduction of organisms.

CONCLUSION

This study showed that differentially expressed genes in the pectoral tissues of Gushi chickens and AA broilers were related to fat metabolism, which affects meat. Additionally, a large number of novel genes were found that may be involved in fat metabolism and thus may affect the formation of meat, which requires further study. The results of this study provide a reference for further studies of the molecular mechanisms of meat formation.

Single nucleotide polymorphism discovery in bovine liver using RNA-seq technology

BACKGROUND

RNA-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.

RESULTS

The 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.

CONCLUSIONS

This 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.

Evaluation of SCD and FASN Gene Expression in Baluchi, Iran-Black, and Arman Sheep

BACKGROUND

With the increasing concern for health and nutrition, dietary fat has attracted considerable attention. The composition of fatty acids in the diet is important because they are associated with major diseases including cancers, diabetes, and cardiovascular disease. The fatty acid synthase (FASN) and stearoyl-CoA desaturase (delta-9-desaturase) (SCD) genes affect fatty acid composition (1). The expression of SCD and FASN genes is related to an increase in conjugated linoleic acid (CLA) in dairy products, which benefits human health.The aim of current study was to investigate expression changes of SCD and FASN genes that resulted from crossbreeding the local Baluchi sheep with alien breeds.

METHODS

We collected tissue samples from the mammary glands of 24 single-born ewes from local Baluchi and synthetic Iran-Black and Arman sheep breeds in the Abbas Abad breeding center. After RNA extraction and cDNA synthesis, real-time PCR was performed with all samples in triplicate.

RESULTS

The maximum and minimum expression of SCD and FASN genes was in the local Baluchi sheep and the crossbred Arman sheep, respectively.

CONCLUSION

With the highest SCD and FASN gene expression in local Baluchi sheep and relatively less expression of these genes in synthetic Iran-Black and Arman Sheep breeds, it may be necessary to consider the consequences of crossbreeding local sheep and the fatty acid composition of their dairy products.

Single Nucleotide Polymorphism Discovery in Bovine Pituitary Gland Using RNA-Seq Technology

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.

The developmental transcriptome landscape of bovine skeletal muscle defined by Ribo-Zero ribonucleic acid sequencing

Ribonucleic acid sequencing (RNA-Seq) libraries are normally prepared with oligo(dT) selection of poly(A)+ mRNA, but it depends on intact total RNA samples. Recent studies have described Ribo-Zero technology, a novel method that can capture both poly(A)+ and poly(A)- transcripts from intact or fragmented RNA samples. We report here the first application of Ribo-Zero RNA-Seq for the analysis of the bovine embryonic, neonatal, and adult skeletal muscle whole transcriptome at an unprecedented depth. Overall, 19,893 genes were found to be expressed, with a high correlation of expression levels between the calf and the adult. Hundreds of genes were found to be highly expressed in the embryo and decreased at least 10-fold after birth, indicating their potential roles in embryonic muscle development. In addition, we present for the first time the analysis of global transcript isoform discovery in bovine skeletal muscle and identified 36,694 transcript isoforms. Transcriptomic data were also analyzed to unravel sequence variations; 185,036 putative SNP and 12,428 putative short insertions-deletions (InDel) were detected. Specifically, many stop-gain, stop-loss, and frameshift mutations were identified that probably change the relative protein production and sequentially affect the gene function. Notably, the numbers of stage-specific transcripts, alternative splicing events, SNP, and InDel were greater in the embryo than in the calf and the adult, suggesting that gene expression is most active in the embryo. The resulting view of the transcriptome at a single-base resolution greatly enhances the comprehensive transcript catalog and uncovers the global trends in gene expression during bovine skeletal muscle development.

Characterisation of the Whole Blood mRNA Transcriptome in Holstein-Friesian and Jersey Calves in Response to Gradual Weaning

Weaning of dairy calves is an early life husbandry management practice which involves the changeover from a liquid to a solid feed based diet. The objectives of the study were to use RNA-seq technology to examine the effect of (i) breed and (ii) gradual weaning, on the whole blood mRNA transcriptome of artificially reared Holstein-Friesian and Jersey calves. The calves were gradually weaned over 14 days (day (d) -13 to d 0) and mRNA transcription was examined one day before gradual weaning was initiated (d -14), one day after weaning (d 1), and 8 days after weaning (d 8). On d -14, 550 genes were differentially expressed between Holstein-Friesian and Jersey calves, while there were 490 differentially expressed genes (DEG) identified on d 1, and 411 DEG detected eight days after weaning (P < 0.05; FDR < 0.1). No genes were differentially expressed within breed, in response to gradual weaning (P > 0.05). The pathways, gene ontology terms, and biological functions consistently over-represented among the DEG between Holstein-Friesian and Jersey were associated with the immune response and immune cell signalling, specifically chemotaxis. Decreased transcription of several cytokines, chemokines, immunoglobulin-like genes, phagocytosis-promoting receptors and g-protein coupled receptors suggests decreased monocyte, natural killer cell, and T lymphocyte, chemotaxis and activation in Jersey compared to Holstein-Friesian calves. Knowledge of breed-specific immune responses could facilitate health management practices better tailored towards specific disease sensitivities of Holstein-Friesian and Jersey calves. Gradual weaning did not compromise the welfare of artificially-reared dairy calves, evidenced by the lack of alterations in the expression of any genes in response to gradual weaning.

Microarray analysis of tick-infested skin in resistant and susceptible cattle confirms the role of inflammatory pathways in immune activation and larval rejection

Tick bites promote activation of an inflammatory process that is influenced by bovine genetic composition and its history of previous exposure. Taurine and indicine breeds are known to differ on its immune response development against Rhipicephalus microplus. Nevertheless, further investigation about the complex molecular pathways involved in the development of immune response to tick infestation in cattle presenting the same genetic background is mandatory. The aim of this work was to access the early immune response triggered by R. microplus larvae attachment in previously selected resistant and susceptible animals in a bovine F2 population derived from Gyr (Bos indicus)×Holstein (Bos taurus) crosses. Microarray data analysis of RNA samples from tick infested skin was used to evaluate the gene expression at 0, 24 and 48h after R. microplus larvae attachment. Our experimental design allowed us to deeply explore the immune response related to R. microplus infestation avoiding the innate differences between these breeds. The differentially expressed genes found reveal networks and pathways that suggest a key role of lipid metabolism in inflammation control and impairment of tick infestation in resistant animals. Acute phase response also seems to be impaired in susceptible animals. These results provide new insights about early immune response against ticks and raise the possibility of using immunomodulation processes to improve and develop novel tools for tick control.

Identification of large intergenic non-coding RNAs in bovine muscle using next-generation transcriptomic sequencing

Background

The advent of large-scale gene expression technologies has helped to reveal in eukaryotic cells, the existence of thousands of non-coding transcripts, whose function and significance remain mostly poorly understood. Among these non-coding transcripts, long non-coding RNAs (lncRNAs) are the least well-studied but are emerging as key regulators of diverse cellular processes. In the present study, we performed a survey in bovine Longissimus thoraci of lincRNAs (long intergenic non-coding RNAs not overlapping protein-coding transcripts). To our knowledge, this represents the first such study in bovine muscle.

Results

To identify lincRNAs, we used paired-end RNA sequencing (RNA-Seq) to explore the transcriptomes of Longissimus thoraci from nine Limousin bull calves. Approximately 14–45 million paired-end reads were obtained per library. A total of 30,548 different transcripts were identified. Using a computational pipeline, we defined a stringent set of 584 different lincRNAs with 418 lincRNAs found in all nine muscle samples. Bovine lincRNAs share characteristics seen in their mammalian counterparts: relatively short transcript and gene lengths, low exon number and significantly lower expression, compared to protein-encoding genes. As for the first time, our study identified lincRNAs from nine different samples from the same tissue, it is possible to analyse the inter-individual variability of the gene expression level of the identified lincRNAs. Interestingly, there was a significant difference when we compared the expression variation of the 418 lincRNAs with the 10,775 known selected protein-encoding genes found in all muscle samples. In addition, we found 2,083 pairs of lincRNA/protein-encoding genes showing a highly significant correlated expression. Fourteen lincRNAs were selected and 13 were validated by RT-PCR. Some of the lincRNAs expressed in muscle are located within quantitative trait loci for meat quality traits.

Conclusions

Our study provides a glimpse into the lincRNA content of bovine muscle and will facilitate future experimental studies to unravel the function of these molecules. It may prove useful to elucidate their effect on mechanisms underlying the genetic variability of meat quality traits. This catalog will complement the list of lincRNAs already discovered in cattle and therefore will help to better annotate the bovine genome.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-499) contains supplementary material, which is available to authorized users.

Comparison of transcriptional profiles of Clostridium thermocellum grown on cellobiose and pretreated yellow poplar using RNA-Seq

The anaerobic, thermophilic bacterium, Clostridium thermocellum, secretes multi-protein enzyme complexes, termed cellulosomes, which synergistically interact with the microbial cell surface and efficiently disassemble plant cell wall biomass. C. thermocellum has also been considered a potential consolidated bioprocessing (CBP) organism due to its ability to produce the biofuel products, hydrogen, and ethanol. We found that C. thermocellum fermentation of pretreated yellow poplar (PYP) produced 30 and 39% of ethanol and hydrogen product concentrations, respectively, compared to fermentation of cellobiose. RNA-seq was used to analyze the transcriptional profiles of these cells. The PYP-grown cells taken for analysis at the late stationary phase showed 1211 genes up-regulated and 314 down-regulated by more than two-fold compared to the cellobiose-grown cells. These affected genes cover a broad spectrum of specific functional categories. The transcriptional analysis was further validated by sub-proteomics data taken from the literature; as well as by quantitative reverse transcription-PCR (qRT-PCR) analyses of selected genes. Specifically, 47 cellulosomal protein-encoding genes, genes for 4 pairs of SigI-RsgI for polysaccharide sensing, 7 cellodextrin ABC transporter genes, and a set of NAD(P)H hydogenase and alcohol dehydrogenase genes were up-regulated for cells growing on PYP compared to cellobiose. These genes could be potential candidates for future studies aimed at gaining insight into the regulatory mechanism of this organism as well as for improvement of C. thermocellum in its role as a CBP organism.

Identification of novel transcripts and noncoding RNAs in bovine skin by deep next generation sequencing

BACKGROUND

Deep RNA sequencing (RNAseq) has opened a new horizon for understanding global gene expression. The functional annotation of non-model mammalian genomes including bovines is still poor compared to that of human and mouse. This particularly applies to tissues without direct significance for milk and meat production, like skin, in spite of its multifunctional relevance for the individual. Thus, applying an RNAseq approach, we performed a whole transcriptome analysis of pigmented and nonpigmented bovine skin to describe the comprehensive transcript catalogue of this tissue.

RESULTS

A total of 39,577 unique primary skin transcripts were mapped to the bovine reference genome assembly. The majority of the transcripts were mapped to known transcriptional units (65%). In addition to the reannotation of known genes, a substantial number (10,884) of unknown transcripts (UTs) were discovered, which had not previously been annotated. The classification of UTs was based on the prediction of their coding potential and comparative sequence analysis, subsequently followed by meticulous manual curation. The classification analysis and experimental validation of selected UTs confirmed that RNAseq data can be used to amend the annotation of known genes by providing evidence for additional exons, untranslated regions or splice variants, by approving genes predicted in silico and by identifying novel bovine loci. A large group of UTs (4,848) was predicted to potentially represent long noncoding RNA (lncRNA). Predominantly, potential lncRNAs mapped in intergenic chromosome regions (4,365) and therefore, were classified as potential intergenic lncRNA. Our analysis revealed that only about 6% of all UTs displayed interspecies conservation and discovered a variety of unknown transcripts without interspecies homology but specific expression in bovine skin.

CONCLUSIONS

The results of our study demonstrate a complex transcript pattern for bovine skin and suggest a possible functional relevance of novel transcripts, including lncRNA, in the modulation of pigmentation processes. The results also indicate that the comprehensive identification and annotation of unknown transcripts from whole transcriptome analysis using RNAseq data remains a tremendous future challenge.

Gene-based single nucleotide polymorphism discovery in bovine muscle using next-generation transcriptomic sequencing

Background

Genetic information based on molecular markers has increasingly being used in cattle breeding improvement programmes, as a mean to improve conventionally phenotypic selection. Advances in molecular genetics have led to the identification of several genetic markers associated with genes affecting economic traits. Until recently, the identification of the causative genetic variants involved in the phenotypes of interest has remained a difficult task. The advent of novel sequencing technologies now offers a new opportunity for the identification of such variants. Despite sequencing costs plummeting, sequencing whole-genomes or large targeted regions is still too expensive for most laboratories. A transcriptomic-based sequencing approach offers a cheaper alternative to identify a large number of polymorphisms and possibly to discover causative variants. In the present study, we performed a gene-based single nucleotide polymorphism (SNP) discovery analysis in bovine Longissimus thoraci, using RNA-Seq. To our knowledge, this represents the first study done in bovine muscle.

Results

Messenger RNAs from Longissimus thoraci from three Limousin bull calves were subjected to high-throughput sequencing. Approximately 36–46 million paired-end reads were obtained per library. A total of 19,752 transcripts were identified and 34,376 different SNPs were detected. Fifty-five percent of the SNPs were found in coding regions and ~22% resulted in an amino acid change. Applying a very stringent SNP quality threshold, we detected 8,407 different high-confidence SNPs, 18% of which are non synonymous coding SNPs. To analyse the accuracy of RNA-Seq technology for SNP detection, 48 SNPs were selected for validation by genotyping. No discrepancies were observed when using the highest SNP probability threshold. To test the usefulness of the identified SNPs, the 48 selected SNPs were assessed by genotyping 93 bovine samples, representing mostly the nine major breeds used in France. Principal component analysis indicates a clear separation between the nine populations.

Conclusions

The RNA-Seq data and the collection of newly discovered coding SNPs improve the genomic resources available for cattle, especially for beef breeds. The large amount of variation present in genes expressed in Limousin Longissimus thoracis, especially the large number of non synonymous coding SNPs, may prove useful to study the mechanisms underlying the genetic variability of meat quality traits.

RNA-seq reveals activation of both common and cytokine-specific pathways following neutrophil priming

Neutrophils are central to the pathology of inflammatory diseases, where they can damage host tissue through release of reactive oxygen metabolites and proteases, and drive inflammation via secretion of cytokines and chemokines. Many cytokines, such as those generated during inflammation, can induce a similar “primed” phenotype in neutrophils, but it is unknown if different cytokines utilise common or cytokine-specific pathways to induce these functional changes. Here, we describe the transcriptomic changes induced in control human neutrophils during priming in vitro with pro-inflammatory cytokines (TNF-α and GM-CSF) using RNA-seq. Priming led to the rapid expression of a common set of transcripts for cytokines, chemokines and cell surface receptors (CXCL1, CXCL2, IL1A, IL1B, IL1RA, ICAM1). However, 580 genes were differentially regulated by TNF-α and GM-CSF treatment, and of these 58 were directly implicated in the control of apoptosis. While these two cytokines both delayed apoptosis, they induced changes in expression of different pro- and anti-apoptotic genes. Bioinformatics analysis predicted that these genes were regulated via differential activation of transcription factors by TNF-α and GM-CSF and these predictions were confirmed using functional assays: inhibition of NF-κB signalling abrogated the protective effect of TNF-α (but not that of GM-CSF) on neutrophil apoptosis, whereas inhibition of JAK/STAT signalling abrogated the anti-apoptotic effect of GM-CSF, but not that of TNF-α (p<0.05). These data provide the first characterisation of the human neutrophil transcriptome following GM-CSF and TNF-α priming, and demonstrate the utility of this approach to define functional changes in neutrophils following cytokine exposure. This may provide an important, new approach to define the molecular properties of neutrophils after in vivo activation during inflammation.

Application of Top-Down and Bottom-up Systems Approaches in Ruminant Physiology and Metabolism

Systems biology is a computational field that has been used for several years across different scientific areas of biological research to uncover the complex interactions occurring in living organisms. Applications of systems concepts at the mammalian genome level are quite challenging, and new complimentary computational/experimental techniques are being introduced. Most recent work applying modern systems biology techniques has been conducted on bacteria, yeast, mouse, and human genomes. However, these concepts and tools are equally applicable to other species including ruminants (e.g., livestock). In systems biology, both bottom-up and top-down approaches are central to assemble information from all levels of biological pathways that must coordinate physiological processes. A bottom-up approach encompasses draft reconstruction, manual curation, network reconstruction through mathematical methods, and validation of these models through literature analysis (i.e., bibliomics). Whereas top-down approach encompasses metabolic network reconstructions using ‘omics’ data (e.g., transcriptomics, proteomics) generated through DNA microarrays, RNA-Seq or other modern high-throughput genomic techniques using appropriate statistical and bioinformatics methodologies. In this review we focus on top-down approach as a means to improve our knowledge of underlying metabolic processes in ruminants in the context of nutrition. We also explore the usefulness of tissue specific reconstructions (e.g., liver and adipose tissue) in cattle as a means to enhance productive efficiency.