Identification and characterization of novel and differentially expressed microRNAs in peripheral blood from healthy and mastitis Holstein cattle by deep sequencing

Review of exosomes and their potential for veterinary medicine

Small extracellular vesicles called exosomes are released by almost all cell types and play a crucial role in both healthy and pathological circumstances. Exosomes, found in biological fluids (including plasma, urine, milk, semen, saliva, abdominal fluid and cervical vaginal fluid) and ranging in size from 50 to 150 nm, are critical for intercellular communication. Analysis of exosomal cargos, including micro RNAs (miRNAs), proteins and lipids, has been proposed as valuable diagnostic and prognostic biomarkers of disease. Exosomes can also be used as novel, cell-free, treatment strategies. In this review, we discuss the role, significance and application of exosomes and their cargos in diseases of animals.

Emerging roles of noncoding micro RNAs and circular RNAs in bovine mastitis: Regulation, breeding, diagnosis, and therapy

Bovine mastitis is one of the most troublesome and costly problems in the modern dairy industry, which is not only difficult to monitor, but can also cause economic losses while having significant implications on public health. However, efficacious preventative methods and therapy are still lacking. Moreover, new drugs and therapeutic targets are in increasing demand due to antibiotic restrictions. In recent years, noncoding RNAs have gained popularity as a topic in pathological and genetic studies. Meanwhile, there is growing evidence that they play a role in regulating various biological processes and developing novel treatment platforms. In light of this, this review focuses on two types of noncoding RNAs, micro RNAs and circular RNAs, and summarizes their characterizations, relationships, potential applications as selection markers, diagnostic or treatment targets and potential applications in RNA-based therapy, in order to shed new light on further research.

Vetinformatics from functional genomics to drug discovery: Insights into decoding complex molecular mechanisms of livestock systems in veterinary science

Having played important roles in human growth and development, livestock animals are regarded as integral parts of society. However, industrialization has depleted natural resources and exacerbated climate change worldwide, spurring the emergence of various diseases that reduce livestock productivity. Meanwhile, a growing human population demands sufficient food to meet their needs, necessitating innovations in veterinary sciences that increase productivity both quantitatively and qualitatively. We have been able to address various challenges facing veterinary and farm systems with new scientific and technological advances, which might open new opportunities for research. Recent breakthroughs in multi-omics platforms have produced a wealth of genetic and genomic data for livestock that must be converted into knowledge for breeding, disease prevention and management, productivity, and sustainability. Vetinformatics is regarded as a new bioinformatics research concept or approach that is revolutionizing the field of veterinary science. It employs an interdisciplinary approach to understand the complex molecular mechanisms of animal systems in order to expedite veterinary research, ensuring food and nutritional security. This review article highlights the background, recent advances, challenges, opportunities, and application of vetinformatics for quality veterinary services.

Non-coding transcriptomic profiles in the sheep mammary gland during different lactation periods

Sheep milk production is a dynamic and multifactorial trait regulated by diverse biological mechanisms. To improve the quality and production of sheep milk, it is necessary to understand the underlying non-coding transcriptomic mechanisms. In this study, ribonucleic acid-sequencing (RNA-seq) was used to profile the expression of microRNAs (miRNAs) and circular RNAs (circRNAs) in the sheep mammary gland at three key lactation time points (perinatal period, PP; early lactation, EL; and peak lactation, PL). A total of 2,369 novel circRNAs and 272 miRNAs were profiled, of which 348, 373, and 36 differentially expressed (DE) circRNAs and 30, 34, and 7 DE miRNAs were detected in the comparison of EL vs. PP, PL vs. PP, and PL vs. EL, respectively. A series of bioinformatics analyses including functional enrichment, machine learning prediction, and competing endogenous RNA (ceRNA) network analyses were conducted to identify subsets of the potential candidate miRNAs (e.g., oar_miR_148a, oar_miR_362, and oar_miR_432) and circRNAs (e.g., novel_circ_0011066, novel_circ_0010460, and novel_circ_0006589) involved in sheep mammary gland development. Taken together, this study offers a window into the dynamics of non-coding transcriptomes that occur during sheep lactation and may provide further insights into miRNA and circRNA that influence sheep mammary gland development.

Emerging Roles of Noncoding RNAs in Bovine Mastitis Diseases

Non-coding RNAs (ncRNAs) are an abundant class of RNA with varying nucleotide lengths. They have been shown to have great potential in eutherians/human disease diagnosis and treatments and are now gaining more importance for the improvement of diseases in livestock. To date, thousands of ncRNAs have been discovered in the bovine genome and the continuous advancement in deep sequencing technologies and various bioinformatics tools has enabled the elucidation of their roles in bovine health. Among farm animals' diseases, mastitis, a common inflammatory disease in cattle, has caused devastating economic losses to dairy farmers over the last few decades. Here, we summarize the biology of bovine mastitis and comprehensively discuss the roles of ncRNAs in different types of mastitis infection. Based on our findings and relevant literature, we highlighted various evidence of ncRNA roles in mastitis. Different approaches (in vivo versus in vitro) for exploring ncRNA roles in mastitis are emphasized. More particularly, the potential applications of emerging genome editing technologies, as well as integrated omics platforms for ncRNA studies and implications for mastitis are presented.

Integrated Analysis of Transcriptome mRNA and miRNA Profiles Reveals Self-Protective Mechanism of Bovine MECs Induced by LPS

Many studies have investigated the molecular crosstalk between mastitis-pathogens and cows by either miRNA or mRNA profiles. Here, we employed both miRNA and mRNA profiles to understand the mechanisms of the response of bovine mammary epithelial cells (bMECs) to lipopolysaccharide (LPS) by RNA-Seq. The total expression level of miRNAs increased while mRNAs reduced after LPS treatment. About 41 differentially expressed mRNAs and 45 differentially expressed miRNAs involved in inflammation were screened out. We found the NFκB-dependent chemokine, CXCL1, CXCL3, CXCL6, IL8, and CX3CL1 to be strongly induced. The anti-apoptosis was active because BCL2A1 and BIRC3 significantly increased with a higher expression. The effects of anti-microbe and inflammation were weakly activated because TNF, IL1, CCL20, CFB, S100A, MMP9, and NOS2A significantly increased but with a low expression, IL6 and β-defensin decreased. These activities were supervised by the NFKBIA to avoid excessive damage to bMECs. The bta-let-7a-5p, bta-miR-30a-5p, bta-miR-125b, and bta-miR-100 were essential to regulate infection process in bMECs after LPS induction. Moreover, the lactation potential of bMECs was undermined due to significantly downregulated SOSTDC1, WNT7B, MSX1, and bta-miR-2425-5p. In summary, bMECs may not be good at going head-to-head with the pathogens; they seem to be mainly charged with sending out signals for help and anti-apoptosis for maintaining lives after LPS induction.

The Role of microRNAs in the Mammary Gland Development, Health, and Function of Cattle, Goats, and Sheep

Milk is an integral and therefore complex structural element of mammalian nutrition. Therefore, it is simple to conclude that lactation, the process of producing milk, is as complex as the mammary gland, the organ responsible for this biochemical activity. Nutrition, genetics, epigenetics, disease pathogens, climatic conditions, and other environmental variables all impact breast productivity. In the last decade, the number of studies devoted to epigenetics has increased dramatically. Reports are increasingly describing the direct participation of microRNAs (miRNAs), small noncoding RNAs that regulate gene expression post-transcriptionally, in the regulation of mammary gland development and function. This paper presents a summary of the current state of knowledge about the roles of miRNAs in mammary gland development, health, and functions, particularly during lactation. The significance of miRNAs in signaling pathways, cellular proliferation, and the lipid metabolism in agricultural ruminants, which are crucial in light of their role in the nutrition of humans as consumers of dairy products, is discussed.

Genome-wide association study on milk production and somatic cell score for Thai dairy cattle using weighted single-step approach with random regression test-day model

Genome-wide association studies are a powerful tool to identify genomic regions and variants associated with phenotypes. However, only limited mutual confirmation from different studies is available. The objectives of this study were to identify genomic regions as well as genes and pathways associated with the first-lactation milk, fat, protein, and total solid yields; fat, protein, and total solid percentage; and somatic cell score (SCS) in a Thai dairy cattle population. Effects of SNPs were estimated by a weighted single-step GWAS, which back-solved the genomic breeding values predicted using single-step genomic BLUP (ssGBLUP) fitting a single-trait random regression test-day model. Genomic regions that explained at least 0.5% of the total genetic variance were selected for further analyses of candidate genes. Despite the small number of genotyped animals, genomic predictions led to an improvement in the accuracy over the traditional BLUP. Genomic predictions using weighted ssGBLUP were slightly better than the ssGBLUP. The genomic regions associated with milk production traits contained 210 candidate genes on 19 chromosomes [Bos taurus autosome (BTA) 1 to 7, 9, 11 to 16, 20 to 21, 26 to 27 and 29], whereas 21 candidate genes on 3 chromosomes (BTA 11, 16, and 21) were associated with SCS. Many genomic regions explained a small fraction of the genetic variance, indicating polygenic inheritance of the studied traits. Several candidate genes coincided with previous reports for milk production traits in Holstein cattle, especially a large region of genes on BTA14. We identified 141 and 5 novel genes related to milk production and SCS, respectively. These novel genes were also found to be functionally related to heat tolerance (e.g., SLC45A2, IRAG1, and LOC101902172), longevity (e.g., SYT10 and LOC101903327), and fertility (e.g., PAG1). These findings may be attributed to indirect selection in our population. Identified biological networks including intracellular cell transportation and protein catabolism implicate milk production, whereas the immunological pathways such as lymphocyte activation are closely related to SCS. Further studies are required to validate our findings before exploiting them in genomic selection.

Characterization of microRNA expression in B cells derived from Japanese black cattle naturally infected with bovine leukemia virus by deep sequencing

Bovine leukemia virus (BLV) is the causative agent of enzootic bovine leukosis (EBL), a malignant B cell lymphoma. However, the mechanisms of BLV-associated lymphomagenesis remain poorly understood. Here, after deep sequencing, we performed comparative analyses of B cell microRNAs (miRNAs) in cattle infected with BLV and those without BLV. In BLV-infected cattle, BLV-derived miRNAs (blv-miRNAs) accounted for 38% of all miRNAs in B cells. Four of these blv-miRNAs (blv-miR-B1-5p, blv-miR-B2-5p, blv-miR-B4-3p, and blv-miR-B5-5p) had highly significant positive correlations with BLV proviral load (PVL). The read counts of 90 host-derived miRNAs (bta-miRNAs) were significantly down-regulated in BLV-infected cattle compared to those in uninfected cattle. Only bta-miR-375 had a positive correlation with PVL in BLV-infected cattle and was highly expressed in the B cell lymphoma tissue of EBL cattle. There were a few bta-miRNAs that correlated with BLV tax/rex gene expression; however, BLV AS1 expression had a significant negative correlation with many of the down-regulated bta-miRNAs that are important for tumor development and/or tumor suppression. These results suggest that BLV promotes lymphomagenesis via AS1 and blv-miRNAs, rather than tax/rex, by down-regulating the expression of bta-miRNAs that have a tumor-suppressing function, and this downregulation is linked to increased PVL.

Expression Profiling of microRNA From Peripheral Blood of Dairy Cows in Response to Staphylococcus aureus-Infected Mastitis

As the main pathogen causing dairy cow mastitis, Staphylococcus aureus can cause subclinical mastitis, which is difficult to be diagnosed. It seriously affects milk quality and the economic benefits of the dairy industry. Therefore, it is very necessary to find biomarkers for early diagnosis of S. aureus-infected mastitis in peripheral blood of dairy cows. In this study, S. aureus was used to infect the mammary gland tissues of dairy cows, and a mastitis model was successfully constructed. The RNAseq technology was used to determine the expression profiles of microRNA (miRNA) from peripheral blood of dairy cows infected with S. aureus at 0, 1, 3, 5, and 7 days. A total of 288 differentially expressed miRNAs (DIE-miRNAs) were found, of which 108 were known miRNAs and 180 were novel predicted miRNAs. Bioinformatics analysis results showed that the above DIE-miRNAs might be involved in 10 immune system-related signaling pathways (i.e., chemokine signaling pathway, leukocyte transendothelial migration, natural killer cell-mediated cytotoxicity, toll-like receptor signaling pathway, Jak-STAT signaling pathway, MAPK signaling pathway, Wnt signaling pathway, cell adhesion molecules, cytokine-cytokine receptor interaction, and ECM-receptor interaction), thus regulating the process of S. aureus mastitis. It was also found that the expression variation of up-regulated expression of miR-320a, miR-19a, and miR-19b as well as down-regulated expression of miR-143, miR-205, and miR-24 reached a significant level on the 5th and 7th day of infection, suggesting that they might play an important biological role in mastitis and provide a direction for the research and development of molecular therapy technology for mastitis. However, at different times after S. aureus infection, miR-1301 was significantly up-regulated in peripheral blood. miR-2284r was significantly down-regulated, suggesting that these two miRNAs might be the new blood biomarkers for S. aureus-infected dairy cow mastitis. The above results laid a new foundation for the research and development of molecular diagnosis and biological therapy technology for S. aureus-infected mastitis in dairy cow.

The Regulation of Staphylococcus aureus-Induced Inflammatory Responses in Bovine Mammary Epithelial Cells

Mastitis, an inflammatory disease, causes severe economic loss in the dairy industry, which is mainly infected by bacteria. Staphylococcus aureus (S. aureus), the major pathogenic microorganism, derived from lipoteichoic acid (LTA) has been identified to activate inflammatory responses, but the cellular or intercellular regulatory mechanism is unclear. This study mainly focused on the effects of LTA in bovine mammary epithelial cells (Mac-T) and elaborated the regulation of microRNAs (miRNAs). The results showed that LTA enhanced the messenger RNA (mRNA) expression and production of tumor necrosis factor α (TNF-α) and interleukin (IL)-6. Furthermore, LTA could activate Toll-like receptor (TLR)2/MyD88-mediated phosphoinositide 3-kinase (PI3K)/AKT pathway, and TLR2 plays a pivotal role in LTA-induced inflammatory responses. The results of qRT-PCR showed that miRNA levels increased and reached the highest at 3 h and then gradually decreased over time in Mac-T cells. In exosomes, the levels of 11 and three miRNAs were upregulated and downregulated at 24 h, respectively. In addition, miR-23a showed the highest increase in Mac-T cells treated with LTA and targeted PI3K to regulate inflammatory responses. Furthermore, Mac-T cell-derived exosomes were identified to play a cell–cell communication by promoting M1 polarization of bovine macrophages. In summary, our study demonstrated that LTA could activate inflammatory responses via TLR2/MyD88/PI3K/AKT signaling pathway, and miR-23a inhibited it by targeting PI3K. Furthermore, we found that Mac-T cell-derived exosomes might be associated with inflammatory responses by promoting M1 polarization of bovine macrophages.

miRNA Regulatory Functions in Farm Animal Diseases, and Biomarker Potentials for Effective Therapies

MicroRNAs (miRNAs) are small endogenous RNAs that regulate gene expression post-transcriptionally by targeting either the 3′ untranslated or coding regions of genes. They have been reported to play key roles in a wide range of biological processes. The recent remarkable developments of transcriptomics technologies, especially next-generation sequencing technologies and advanced bioinformatics tools, allow more in-depth exploration of messenger RNAs (mRNAs) and non-coding RNAs (ncRNAs), including miRNAs. These technologies have offered great opportunities for a deeper exploration of miRNA involvement in farm animal diseases, as well as livestock productivity and welfare. In this review, we provide an overview of the current knowledge of miRNA roles in major farm animal diseases with a particular focus on diseases of economic importance. In addition, we discuss the steps and future perspectives of using miRNAs as biomarkers and molecular therapy for livestock disease management as well as the challenges and opportunities for understanding the regulatory mechanisms of miRNAs related to disease pathogenesis.

MicroRNAs as Biomarkers for Animal Health and Welfare in Livestock

MicroRNAs (miRNAs) are small and highly conserved non-coding RNA molecules that orchestrate a wide range of biological processes through the post-transcriptional regulation of gene expression. An intriguing aspect in identifying these molecules as biomarkers is derived from their role in cell-to-cell communication, their active secretion from cells into the extracellular environment, their high stability in body fluids, and their ease of collection. All these features confer on miRNAs the potential to become a non-invasive tool to score animal welfare. There is growing interest in the importance of miRNAs as biomarkers for assessing the welfare of livestock during metabolic, environmental, and management stress, particularly in ruminants, pigs, and poultry. This review provides an overview of the current knowledge regarding the potential use of tissue and/or circulating miRNAs as biomarkers for the assessment of the health and welfare status in these livestock species.

Anti-mastitis SNV identification of NFκB1 in Chinese Holstein cows and the possible anti-inflammation role of NFκB1/p105 in bovine MECs

NFκB1/p105 is the critical member of the NFκB family which can suppress inflammation, ageing, and cancer when p50/p50 homodimer is formed. Currently, the research about the role of NFκB1/p105 during cow mastitis is limited. Here, we analyzed the correlation of six single-nucleotide variants of the NFκB1 gene with somatic cell count, milk yield, milk fat content, and milk protein content in 547 Chinese Holstein cows, and explored the mRNA expression profiles of the NFκB family and ubiquitin ligases (βTrCP1, βTrCP2, KPC1, KPC2) in LPS-induced bovine mammary epithelial cells (MECs) by transcriptome-Seq. The association analysis showed that cows with SNV2-TT and SNV6-CC in the NFκB1 gene had significantly higher milk protein content (P < 0.05), while cows with SNV5-TT in the NFκB1 gene had significantly lower somatic cell score (SCS), but CC genotype at SNV5 locus was not detected in our Holstein cows. The transcriptome-Seq results demonstrated the mRNA expression of NFκB1 was increased and peaked at 4 h post-induction, while the mRNA expressions of both KPC1 and BCL3 that promote the anti-inflammation function of NFκB1/p105 were decreased in LPS-induced bovine MECs. TNFAIP3, an inhibitor of both degradation and processing of p105 precursor, was markedly increased by more than 3 folds. Furthermore, bta-miR-125b which targets at the 3'UTR of TNFAIP3 was reduced by 50%. These results indicated that SNV5-TT of the NFκB1 gene with lower SCS may be an anti-mastitis genotype that could cope with infection more efficiently in Chinese Holstein cows. In addition, the anti-inflammation role of NFκB1/p105 seemed to be inhibited in LPS-induced-bovine MECs because the formation of the p50/p50 homodimer was arrested. This study provides a new perspective to understand the inflammatory mechanism in dairy mastitis.

Bovine serum miR-21 expression affected by mastitis

The expression levels of circulating microRNAs (miRNAs) can be affected by disease. The miRNA released from cells within exosomes can act as a remote communication tool and can participate in inflammatory response regulation. Therefore, circulating miRNA has the potential to be an indicator of local disease. The objective of this study was to investigate the serum level of bovine mastitis-related miRNAs. We found that miR-16 expression in serum was affected by hemolysis. The expression levels of miR-21 in serum were increased significantly in cows with mastitis compared with unaffected controls; however, the expression levels of miR-146a, miR-155, miR-222 and miR-383 in cows with mastitis were unchanged. We further verified the upregulation of miR-21 in the serum of cows with mastitis using a digital PCR system. Although the sensitivity and specificity of miR-21 in the serum to detect bovine mastitis was inferior to miRNA biomarkers in the milk, the significant increase of miR-21 in serum may reflect the impact of local inflammation on the systemic reaction.

Bta-miR-223 Targeting CBLB Contributes to Resistance to Staphylococcus aureus Mastitis Through the PI3K/AKT/NF-κB Pathway

Bovine mastitis is an inflammatory condition of the mammary gland often caused by (Staphylococcus aureus) S. aureus infection. The aim of this study was to identify mastitis-related miRNAs and their downstream target genes, and therefore elucidate the regulatory mechanisms involved in disease progression and resistance. Three healthy and three mastitic cows were identified on the basis of the somatic cell count and bacterial culture of their milk, and the histological examination of udder tissues. High-throughput RNA sequencing and bioinformatic analyses revealed that 48 differentially expressed miRNAs (DEMs) in the mastitic udder tissues relative to the healthy tissues. Among 48 DEMs, the expression level of bta-miR-223 was the most up-regulated. Overexpression of the bta-miR-223 in Mac-T cells mitigated the inflammatory pathways induced by S. aureus-derived lipoteichoic acid (LTA). The Cbl proto-oncogene B (CBLB) was identified as the target gene of bta-miR-223, and the direct binding of the miRNA to the CBLB promoter was confirmed by dual luciferase reporter assay using wild-type and mutant 3'-UTR constructs. Furthermore, overexpression of CBLB in the LTA-stimulated Mac-T cells significantly upregulated PI3K, AKT, and phosphorylated NF-κB p65, whereas CBLB knockdown had the opposite effect. Consistent with the in vitro findings, the mammary glands of mice infected with 10⁸CFU/100 μL S. aureus showed high levels of CBLB, PI3K, AKT, and p-NF-κB p65 48 h after infection. Taken together, bta-miR-223 is a predominant miRNA involved in mastitis, and bta-miR-223 likely mitigates the inflammatory progression by targeting CBLB and inhibiting the downstream PI3K/AKT/NF-κB pathway.

Expression profiling of microRNAs in the Mycoplasma bovis infected mammary gland tissue in Holstein Friesian cattle

The immune response associated with mastitis caused by Mycoplasma bovis is a very complicated biological process in several type of cells, including immune cells, mammary epithelial cells, and endothelial cells. Thus, revealing of the microRNAs in the Mycoplasma bovis infected mammary gland tissues is particularly important for the immune response mechanism to Mycoplasma bovis. Firstly, 20 mammary gland tissue samples were collected from Holstein Friesian cattle that was located in Erzurum province at 2018 and screened for Mycoplasma bovis. Then, total RNA was isolated from Mycoplasma bovis infected tissues and high-throughput sequencing was performed. After bioinformatics analysis, GO and KEGG analysis of target genes of identified microRNAs were conducted. In this study, a total of 616 microRNAs were found. Our results revealed that 24 of the known microRNAs were expressed differently and 13 of the novel microRNAs were expressed differently in Mycoplasma bovis positive tissues. The target genes of these microRNAs were found to be associated with especially inflammation pathways, including B cell and T cell receptor signaling, Fc gamma R-mediated, phagocytosis/chemokine signaling, and MAPK signaling. In conclusion, this study demonstrated that identified miRNAs may be involved in the signaling pathways during mastitis caused by Mycoplasma bovis.

Discovery and annotation of novel microRNAs in the porcine genome by using a semi-supervised transductive learning approach

Despite the broad variety of available microRNA (miRNA) prediction tools, their application to the discovery and annotation of novel miRNA genes in domestic species is still limited. In this study we designed a comprehensive pipeline (eMIRNA) for miRNA identification in the yet poorly annotated porcine genome and demonstrated the usefulness of implementing a motif search positional refinement strategy for the accurate determination of precursor miRNA boundaries. The small RNA fraction from gluteus medius skeletal muscle of 48 Duroc gilts was sequenced and used for the prediction of novel miRNA loci. Additionally, we selected the human miRNA annotation for a homology-based search of porcine miRNAs with orthologous genes in the human genome. A total of 20 novel expressed miRNAs were identified in the porcine muscle transcriptome and 27 additional novel porcine miRNAs were also detected by homology-based search using the human miRNA annotation. The existence of three selected novel miRNAs (ssc-miR-483, ssc-miR484 and ssc-miR-200a) was further confirmed by reverse transcription quantitative real-time PCR analyses in the muscle and liver tissues of Göttingen minipigs. In summary, the eMIRNA pipeline presented in the current work allowed us to expand the catalogue of porcine miRNAs and showed better performance than other commonly used miRNA prediction approaches. More importantly, the flexibility of our pipeline makes possible its application in other yet poorly annotated non-model species.

Bovine milk transcriptome analysis reveals microRNAs and RNU2 involved in mastitis

Mastitis is a common inflammatory infectious disease in dairy cows. To understand the microRNA (miRNA) expression profile changes during bovine mastitis, we undertook a genome-wide miRNA study of normal milk and milk that tested positive on the California mastitis test for bovine mastitis (CMT+). Twenty-five miRNAs were differentially expressed (23 miRNAs upregulated and two downregulated) during bovine mastitis relative to their expression in normal milk. Upregulated mature miR-1246 probably derived from a U2 small nuclear RNA rather than an miR-1246 precursor. The significantly upregulated miRNA precursors and RNU2 were significantly enriched on bovine chromosome 19, which is homologous to human chromosome 17. A gene ontology analysis of the putative mRNA targets of the significantly upregulated miRNAs showed that these miRNAs were involved in binding target mRNA transcripts and regulating target gene expression, and a Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that the upregulated miRNAs were predominantly related to cancer and immune system pathways. Three novel miRNAs were associated with bovine mastitis and were relatively highly expressed in milk. We confirmed that one of the novel mastitis-related miRNAs was significantly upregulated using a digital PCR system. The differentially expressed miRNAs were involved in human cancers, infections, and immune-related diseases. The genome-wide analysis of miRNA profiles in this study provides insight into bovine mastitis and inflammatory diseases. DATABASES: The miRNAseq generated for this study can be found in the Sequence Read Archive (SRA) under BioProject Number PRJNA421075 and SRA Study Number SRP126134 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA421075).

Characterizing binding sites of heat responsive microRNAs and their expression pattern in heat stressed PBMCs of native cattle, exotic cattle and riverine buffaloes

It is generally believed that due to evolutionary differences and adaptation to tropical conditions, Indian native cattle has superior heat tolerant ability than Bos taurus cattle. In the present study, 3'-UTR of two most important heat responsive genes i.e., heat shock protein 70.1 (HSP70.1) and heat shock factor- 1 (HSF-1) were sequence characterized in different breeds of Indian native cattle to identify the variations and miRNA binding sites. In addition, the impact of heat stress was assessed in a total of 57 PBMCs samples of native Sahiwal cows (Bos indicus), exotic Holstein cows (Bos taurus) and Murrah buffaloes (Bubalus bubalis) using various cellular parameters like cell viability, cytotoxicity and apoptosis. Further, expression profile of 12 heat responsive miRNAs were also evaluated in unstressed and stressed PBMCs to understand post transcriptional changes in native cows, exotic cows and Murrah buffaloes. The sequence data showed 3'-UTR of HSP70.1 gene of Indian cattle to be exactly similar to Bos taurus with no miRNA binding site. Whereas, sequencing of 3'-UTR of HSF-1 gene revealed 3 SNPs at positions G1762T; C1811T and C1983T with 7 well conserved miRNA binding sites. The impact of heat stress on various cellular parameters in terms of cell viability, cytotoxicity and apoptosis was highest in PBMCs of Holstein cows followed by Murrah buffaloes and Sahiwal cows. Further, in contrast to Holstein Frisian cows and Murrah buffaloes, the expression pattern of 12 heat responsive miRNAs, in heat stressed PBMCs of Sahiwal cows were quite distinct. There was a significant (p < 0.05) induction in expression of most of the miRNAs after heat stress in PBMCs of Sahiwal cows followed by a rapid decline. The distinct cellular response and pattern of miRNA expression across cattle types and buffaloes might be influencing their PBMCs tolerance level to heat stress.

In silico identification of conserved miRNAs and their selective target gene prediction in indicine (Bos indicus) cattle

The modern cattle was domesticated from aurochs, sharing its physiological traits into two subspecies Bos taurus and Bos indicus. MicroRNAs (miRNAs) are a class of non-coding short RNAs of ~22nt which have a key role in the regulation of many cellular and physiological processes in the animal. The current study was aimed to predict and annotate the potential mutations in indicine miRNAs throughout the genome using de novo and homology-based in silico approaches. Genome-wide mapping was performed in available indicine assembly by the homology-based approach and 768 miRNAs were recovered out of 808 reported taurine miRNAs belonging to 521 unique mature miRNA families. While 42 precursors were dropped due to lack of secondary miRNA structure, increasing stringency or decreasing similarity between the two genomes' miRNA. Increasing tendency of miRNAs incidence was observed on chr5, chr7, chr8, chr12 and chr21 with 19 polycistronic miRNA within 1-kilobase distance throughout the indicine genome. Notably, 12 miRNAs showed copy number variation. Eighteen miRNAs showed a mutation in their mature sequences in which eight were found in their seed region. Whilst in de novo based approach, 12 novel potential miRNAs on Y chromosome in indicine cattle along with a new miRNA (bind-miR-1264) on chrX were found. The final data set is annotated and explains the impending target genes that are responsible for enhanced immunity, heat tolerance and disease tolerance regulation in indicine. The study conforms to better understanding and perceptive approach towards indicine genome.

Expression profiling of peripheral blood miRNA using RNAseq technology in dairy cows with Escherichia coli-induced mastitis

E. coli is the main causative agent of mastitis in dairy cows, but the mechanism of molecular regulation underlying the occurrence and development of mastitis has not yet been fully elucidated. In this study, an E. coli-induced mastitis model was created and RNASeq technology was used to measure the miRNA expression profiles at different times post-infection (0, 1, 3, 5, 7 dpi), as well as to screen for differentially expressed miRNA. The results show detection of 2416 miRNAs, including 628 known miRNAs and 1788 newly discovered miRNAs. A total of 200 differentially expressed miRNAs were found at different time points. Bioinformatics analysis showed that these differentially expressed miRNAs may regulate the occurrence and development of mastitis in dairy cows through seven signal transduction pathways, namely cytokine-cytokine receptor interaction, MAPK signaling pathway, chemokine signaling pathway, leukocyte transendothelial migration, T cell receptor signaling pathway, Toll-like receptor signaling pathway, and cell adhesion molecules. In addition, bta-miR-200a, bta-miR-205, bta-miR-122, bta-miR-182 and the newly discovered conservative_15_7229 might be involved in immune process in late stage of E. coli-induced mastitis. The results of this study lay the foundation for molecular network analysis of mastitis and molecular breeding of dairy cows.

Solexa sequencing and custom microRNA chip reveal repertoire of microRNAs in mammary gland of bovine suffering from natural infectious mastitis

Identification of microRNAs (miRNAs), target genes and regulatory networks associated with innate immune and inflammatory responses and tissue damage is essential to elucidate the molecular and genetic mechanisms for resistance to mastitis. In this study, a combination of Solexa sequencing and custom miRNA chip approaches was used to profile the expression of miRNAs in bovine mammary gland at the late stage of natural infection with Staphylococcus aureus, a widespread mastitis pathogen. We found 383 loci corresponding to 277 known and 49 putative novel miRNAs, two potential mitrons and 266 differentially expressed miRNAs in the healthy and mastitic cows' mammary glands. Several interaction networks and regulators involved in mastitis susceptibility, such as ALCAM, COL1A1, APOP4, ITIH4, CRP and fibrinogen alpha (FGA), were highlighted. Significant down-regulation and location of bta-miR-26a, which targets FGA in the mastitic mammary glands, were validated using quantitative real-time PCR, in situ hybridization and dual-luciferase reporter assays. We propose that the observed miRNA variations in mammary glands of mastitic cows are related to the maintenance of immune and defense responses, cell proliferation and apoptosis, and tissue injury and healing during the late stage of infection. Furthermore, the effect of bta-miR-26a in mastitis, mediated at least in part by enhancing FGA expression, involves host defense, inflammation and tissue damage.

Using RNA sequencing to identify putative competing endogenous RNAs (ceRNAs) potentially regulating fat metabolism in bovine liver

RNA sequencing has been extensively used to study specific gene expression patterns to discover potential key genes related to complex traits of interest in animals. Of note, a new regulatory mechanism builds a large-scale regulatory network among transcriptome, where lncRNAs act as competing endogenous RNAs (ceRNAs) to sponge miRNAs to regulate the expression of miRNA target genes post-transcriptionally. In this study, we sequenced the cDNA and sRNA libraries of nine liver samples from three Holstein cows during dry period, early lactation, and peak of lactation with HiSeq platform. As a result, we identified 665 genes, 57 miRNAs and 33 lncRNAs that displayed differential expression patterns across periods. Subsequently, a total of 41ceRNA pairs (lncRNA-mRNA) sharing 11 miRNAs were constructed including 30 differentially expressed genes. Importantly, 12 among them were presented in our large metabolic networks, and predicted to influence the lipid metabolism through insulin, PI3K-Akt, MAPK, AMPK, mTOR, and PPAR signaling pathways, thus, these genes were considered as the most promising candidates for milk fat formation. To our knowledge, this is first investigation to profile the ceRNA regulatory networks of liver transcriptome that could affect milk fat synthesis in bovine, providing a new view of the regulatory mechanism of RNAs.

Identification and expression analysis of <i>miR-144-5p</i> and <i>miR-130b-5p</i> in dairy cattle

Abstract. MicroRNAs (miRNAs) can coordinate the main pathways involved in innate and adaptive immune responses by regulating gene expression. To explore the resistance to mastitis in cows, miR-144-5p and miR-130b-5p were identified in bovine mammary gland tissue and 14 potential target genes belonging to the chemokine signaling pathway, the arginine and proline metabolism pathway and the mRNA surveillance pathway were predicted. Subsequently, we estimated the relative expression of miR-144-5p and miR-130b-5p in cow mammary tissues by using stem-loop quantitative real-time polymerase chain reaction. The results showed that the relative expression of miR-144-5p and miR-130b-5p in the mastitis-infected mammary tissues (n = 5) was significantly downregulated 0.14-fold (p < 0. 01) and upregulated 3.34-fold (p < 0. 01), respectively, compared to healthy tissues (n = 5). Our findings reveal that miR-144-5p and miR-130b-5p may have important roles in resistance to mastitis in dairy cattle.

MicroRNA-214 regulates immunity-related genes in bovine mammary epithelial cells by targeting NFATc3 and TRAF3

In human, microRNA-214 (miR-214) plays crucial roles in mechanisms of immunity. However, the potential importance of miR-214 in immune mechanisms in dairy cows has not been investigated. In this study, we assessed potential immunity-related functions of miR-214 in human 293A cells and in bovine mammary epithelial cells (BMECs). We found that NFATc3 and TRAF3 could be targeted by miR-214 in both 293A cells and BMECs. We also found that miR-214 indirectly inhibited the expression of MAP3K14, TBK1 and inflammatory cytokines IL-6 and IL-1β. Taken together, our data revealed miR-214 regulated immunity-related genes by targeting NFATc3 and TRAF3, which provides insight into the molecular basis of immunity.

Identification and characterization of known and novel microRNAs in three tissues of Chinese giant salamander base on deep sequencing approach

MicroRNAs (miRNAs) are a class of non-coding small RNAs that play important roles in the regulation of gene expression. Although animal miRNAs have been extensively studied in model systems, less is known in other animal with limited genome sequence data, including Chinese giant salamander (Andrias davidianus). The identification of the full repertoire of miRNAs expressed in the liver, spleen and muscle of A. davidianus would significantly increase our understanding for physiological function of A. davidianus, in this ancient and endangered urodele amphibian. In this study, three independent small RNA libraries were constructed from the liver, spleen and muscle of A. davidianus. The libraries were subjected to high-throughput sequencing by using the Illumina deep sequencing. As a result, a total of 12,831,239, 13,592,195 and 9,887,531 raw reads representing 2,240,771, 1,363,266 and 1,964,252 clean reads per library were obtained separately. Through bioinformatics analysis, we identified total of 553 known miRNAs and 44 putative novel miRNAs in our small RNA dataset from liver, spleen and muscle tissues. Five known miRNAs (gga-miR-10a-5p, pma-miR-29d-5p, aca-miR-338-3p, hsa-miR-455-3p and ssa-miR-2184-5p_R-1) and three novel miRNAs (PC-5p-891_1763, PC-5p-32538_50 and PC-3p-33645_48) showed different expression in eight different tissues as revealed by stem-loop qPCR analysis. This study characterized the miRNA of A. davidianus for the first time, which provides an opportunity for further understanding of miRNA regulation function in A. davidianus ranked as living fossils.

MicroRNA Transcriptome of Poly I:C-Stimulated Peripheral Blood Mononuclear Cells Reveals Evidence for MicroRNAs in Regulating Host Response to RNA Viruses in Pigs

MicroRNAs (miRNAs) are one family of small noncoding RNAs that function to modulate the activity of specific mRNA targets in animals. To understand the role of miRNAs in regulating genes involved in the host immune response to RNA viruses, we profiled and characterized the miRNAs of swine peripheral blood mononuclear cells (PBMC) stimulated with poly I:C, a synthetic dsRNA analog, by miRNA-sequencing (miRNA-seq). We identified a total of 905 miRNAs, of which 503 miRNAs were firstly exploited herein with no annotation in the latest miRBase 21.0. Expression analysis demonstrated that poly I:C stimulation can elicit significantly differentially expressed (DE) miRNAs in Dapulian (n = 20), one Chinese indigenous breed, as well as Landrace (n = 23). By integrating the mRNA expression profiles of the same sample with miRNA profiles, we carried out function analyses of the target genes of these DE miRNAs, with the results indicating that target genes were most enriched in some immune-related pathways and gene ontology (GO) terms, suggesting that DE miRNAs play an important role in the regulation of host to poly I:C stimulation. Furthermore, we also detected 43 and 61 significantly DE miRNAs between the two breeds in the control sample groups and poly I:C stimulation groups, respectively, which may be involved in regulation of the different characteristics of the two breeds. This study describes for the first time the PBMC miRNA transcriptomic response to poly I:C stimulation in pigs, which not only contributes to a broad view of the pig miRNAome but improves our understanding of miRNA function in regulating host immune response to RNA viruses.

Expression patterns of miR-146a and miR-146b in mastitis infected dairy cattle

This study reports a significant up-regulation of bta-miR-146a and bta-miR-146b expression levels in bovine mammary tissues infected with subclinical, clinical and experimental mastitis. Potential target genes are involved in multiple immunological pathways. These results suggest a regulatory function of both miRNAs for the bovine inflammatory response in mammary tissue.

Identification and bioinformatics analysis of microRNAs associated with stress and immune response in serum of heat-stressed and normal Holstein cows

MicroRNAs (miRNAs) are small single-stranded non-coding RNAs that have an important regulatory function in animal growth and developmental processes. However, the differential expression of miRNA and the role of these miRNAs in heat-stressed Holstein cows are still unknown. In this study, the profile of differentially expressed miRNAs and the target genes analysis in the serum of heat-stressed and normal Holstein cows were investigated by a Solexa deep-sequencing approach and bioinformatics. The data identified 52 differentially expressed miRNAs in 486 known miRNAs which were changed significantly between heat-stressed and normal Holstein cows (fold change >2, P < 0.001). Target genes analysis showed that at least 7 miRNAs (miR-19a, miR-19b, miR-146a, miR-30a-5p, miR-345-3p, miR-199a-3p, and miR-1246) were involved in the response to stress, oxidative stress, development of the immune system, and immune response among the identified 52 differentially expressed miRNAs. Five miRNAs (miR-27b, miR-181a, miR-181b, miR-26a, and miR-146b) were involved in stress and immune responses and the expression of five miRNAs was striking (P < 0.001). In addition, RT-qPCR and deep-sequencing methods showed that 8 miRNAs among the 12 selected miRNAs (miR-19a, miR-19b, miR-27b, miR-30a-5p, miR-181a, miR-181b, miR-345-3p, and miR-1246) were highly expressed in the serum of heat-stressed Holstein cows. GO and KEGG pathway analysis showed that these differentially expressed miRNAs were involved in a pathway that may differentially regulate the expression of stress response and immune response genes. Our study provides an overview of miRNAs expression profile and the interaction between miRNAs and their target genes, which will lead to further understanding of the important roles of miRNAs in heat-stressed Holstein cows.