RumimiR: a detailed microRNA database focused on ruminant species

Genetic parameters, genomic prediction, and identification of regulatory regions located on chromosome 14 for weight traits in Nellore cattle

This study aimed to investigate functional variants in chromosome 14 (BTA14) and its impact in genomic selection for birth weight (BW), weaning weight (WW), and yearling weight (YW) in Nellore cattle. Genetic parameter estimation and the weighted single-step genomic best linear unbiased prediction (WssGBLUP) analyses were performed. Direct additive heritability estimates were high for WW and YW, and moderate for BW. Trait-associated variants distributed across multiple regions on BTA14 were observed in the weighted single-step genome-wide association studies (WssGWAS) results, implying a polygenic genetic architecture for weight in different ages. Several genes have been found in association with the weight traits, including the CUB And Sushi multiple domains 3 (CSMD3), thyroglobulin (TG), and diacylglycerol O-acyltransferase 1 (DGAT1) genes. The variance explained per SNP was higher in six functional classes of gene regulatory regions (5UTR, CpG islands, downstream, upstream, long non-coding RNA, and transcription factor binding sites (TFBS)), highlighting their importance for weight traits in Nellore cattle. A marginal increase in accuracy was observed when the selected functional variants (SV) information was considered in the WssGBLUP method, probably because of the small number of SV available on BTA14. The identified genes, pathways, and functions contribute to a better understanding of the genetic and physiological mechanisms regulating weight traits in the Nellore breed.

Unraveling the microRNAs Involved in Fasciolosis: Master Regulators of the Host-Parasite Crosstalk

Fasciolosis is a neglected tropical disease caused by helminth parasites of the genus Fasciola spp., including Fasciola hepatica (F. hepatica) and Fasciola gigantica (F. gigantica), being a major zoonotic problem of human and animal health. Its control with antihelminthics is becoming ineffective due to the increase in parasite resistance. Developing new therapeutic protocols is crucial to a deeper knowledge of the molecular bases in the host-parasite interactions. The high-throughput omics technologies have dramatically provided unprecedented insights into the complexity of the molecular host-parasite crosstalk. MicroRNAs (miRNAs) are key players as critical regulators in numerous biological processes, modifying the gene expression of cells by degradation of messenger RNA (mRNA), regulating transcription and translation functions, protein positioning, cell cycle integrity, differentiation and apoptosis. The large-scale exploration of miRNAs, including the miRNome, has offered great scientific knowledge of steps in fasciolosis, further scrutinizing the pathogenesis, the growth and development of their strains and their interaction with the host for the survival of the different parasite stages. This review compiles the updated knowledge related to miRNAs involved in fasciolosis and the generated miRNome, highlighting the importance of these key molecules in the host-parasite interactions and the pathogenesis of Fasciola spp. directing towards the development of new biotherapeutic protocols for the control of fasciolosis.

Association Between Dairy Production System and Milk Functionality Based on Analysis of miRNAs in Exosomes from Milk

Dairy farming practices significantly affect the nutritional and functional properties of milk. This study compares miRNAs in milk exosomes from extensive and intensive dairy systems and explores their potential implications for human consumers. Extensive systems are believed to produce milk of higher quality with better animal welfare compared to intensive systems. Milk samples from eight extensive and nine intensive dairy farms were analysed. Milk-derived extracellular vesicles were isolated through sequential ultracentrifugation and characterised through Dynamic Light Scattering and Nanosight to determine the size and the concentration of the extracellular vesicles, in addition to immunoblotting to ensure the presence of exosome-specific proteins in their membrane. miRNA levels were quantified using RT-qPCR, and metabolic pathways associated with miRNAs showing significant differences between farm groups were analysed. EVs from extensive farms had higher concentrations. Notably, bta-miR-451 levels were significantly higher in milk from extensive farms (p = 0.021). Like human miRNA hsa-miR-451, it is linked to pathways related to Parkinson's disease and cancer. Our research suggests that milk production in extensive systems not only provides socioeconomic and environmental benefits but may also have positive effects on human health. Further research is warranted to explore the bioactive potential of these miRNAs and their implications for human health.

Profiling several key milk miRNAs and analysing their signalling pathways in dairy sheep breeds during peak and late lactation

BACKGROUND

The comprehensive understanding of microRNAs (miRNAs) in sheep milk during various lactation periods and their impact on milk yield and composition remains limited.

OBJECTIVES

This study aimed to investigate the expression patterns of four highly expressed miRNAs in sheep milk and their association with milk composition and yield parameters during peak and late lactation stages.

METHODS

A total of 40 healthy 4-year-old Akkaraman (n = 20) and Awassi (n = 20) ewes registered with the Ministry of Agriculture and Forestry of the Republic of Türkiye were used in the present study. For miRNA isolation from milk, the Qiagen miRNeasy Serum/Plasma Advanced Kit was utilised following the manufacturer's instructions. The expression levels of miRNAs were assessed using Qiagen miRNA PCR Assays.

RESULTS

The significant fold changes in the expression levels of oar-miR-30a-5p, oar-miR-148a and oar-miR-181a were observed between peak and late lactation periods in the Awassi sheep breed. Conversely, only oar-miR-30a-5p and oar-miR-148a exhibited statistically significant changes in the Akkaraman sheep breed during the same lactation periods. Furthermore, oar-miR-21-5p demonstrated a significant fold change exclusively in peak lactation compared to Akkaraman and Awassi ewes.

CONCLUSIONS

The findings suggest that the expression of the analysed miRNAs is influenced by both the lactation stage and different sheep breeds. This study offers valuable insights into the relationship between key miRNA expressions in sheep milk and milk composition and yield parameters during peak and late lactation, contributing to the existing knowledge in this field.

The Potential Roles of Host Cell miRNAs in Fine-Tuning Bovine Coronavirus (BCoV) Molecular Pathogenesis, Tissue Tropism, and Immune Regulation

Bovine coronavirus (BCoV) infection causes significant economic loss to the dairy and beef industries worldwide. BCoV exhibits dual tropism, infecting the respiratory and enteric tracts of cattle. The enteric BCoV isolates could also induce respiratory manifestations under certain circumstances. However, the mechanism of this dual tropism of BCoV infection has not yet been studied well. MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression and play a dual role in virus infection, mediating virus or modulating host immune regulatory genes through complex virus-host cell interactions. However, their role in BCoV infection remains unclear. This study aims to identify bovine miRNAs crucial for regulating virus-host interaction, influencing tissue tropism, and explore their potential as biomarkers and therapeutic agents against BCoV. We downloaded 18 full-length BCoV genomes (10 enteric and eight respiratory) from GenBank. We applied several bioinformatic tools to study the host miRNAs targeting various regions in the viral genome. We used the criteria of differential targeting between the enteric/respiratory isolates to identify some critical miRNAs as biological markers for BCoV infection. Using various online bioinformatic tools, we also searched for host miRNA target genes involved in BCoV infection, immune evasion, and regulation. Our results show that four bovine miRNAs (miR-2375, miR-193a-3p, miR-12059, and miR-494) potentially target the BCoV spike protein at multiple sites. These miRNAs also regulate the host immune suppressor pathways, which negatively impacts BCoV replication. Furthermore, we found that bta-(miR-2338, miR-6535, miR-2392, and miR-12054) also target the BCoV genome at certain regions but are involved in regulating host immune signal transduction pathways, i.e., type I interferon (IFN) and retinoic acid-inducible gene I (RIG-I) pathways. Moreover, both miR-2338 and miR-2392 also target host transcriptional factors RORA, YY1, and HLF, which are potential diagnostic markers for BCoV infection. Therefore, miR-2338, miR-6535, miR-2392, and miR-12054 have the potential to fine-tune BCoV tropism and immune evasion and enhance viral pathogenesis. Our results indicate that host miRNAs play essential roles in the BCoV tissue tropism, pathogenesis, and immune regulation. Four bovine miRNAs (miR-2375, bta-miR-193a-3p, bta-miR-12059, and bta-miR-494) target BCoV-S glycoprotein and are potentially involved in several immune suppression pathways during the viral infection. These miRNA candidates could serve as good genetic markers for BCoV infection. However, further studies are urgently needed to validate these identified miRNAs and their target genes in the context of BCoV infection and dual tropism and as genetic markers.

Milk-derived exosome nanovesicles: recent progress and daunting hurdles

Raw milk is the foundation of quality and safety in the dairy industry, and improving milk source management is the fundamental guarantee. Milk-derived exosomes (MDEs) are nanoscale information transfer molecules secreted by mammary cells with unique content and high stability, which can be used not only as potential markers to analyze key traits of lactation, reproduction, nutrition and health of animals, but also help farm managers to take timely interventions to improve animal welfare, milk quality, and functional traits. Our review first outlines the latest advances in MDEs isolation and purification, compositional analysis and characterization tools. We then provide a comprehensive summary of recent applications of MDEs liquid biopsy in breed selection, disease prevention and control, and feeding management. Finally, we evaluate the impact of processing on the stability of MDEs to offer guidance for dairy production and storage. The limitations and challenges in the development and use of MDEs markers are also discussed. As a noninvasive marker with high sensitivity and specificity, the MDEs-mediated assay technology is expected to be a powerful tool for measuring cow health and raw milk quality, enabling dynamic and precise regulation of dairy cows and full traceability of raw milk.

MiRNome variations in milk fractions during feed restrictions of different intensities in dairy cows

BACKGROUND

In dairy cows, diet is one factor that can affect their milk production and composition. However, the effect of feed restriction on milk miRNome has not yet been described. Indeed, milk is the body fluid with the highest RNA concentration, which includes numerous microRNA. Its presence in the four different milk fractions, whole milk, fat globules, mammary epithelial cells and extracellular vesicles, is still poorly documented. This study aimed to describe the effects of different feed restrictions on the miRNome composition of different milk fractions.

RESULTS

Two feed restrictions were applied to lactating dairy cows, one of high intensity and one of moderate intensity. 2,896 mature microRNA were identified in the different milk fractions studied, including 1,493 that were already known in the bovine species. Among the 1,096 microRNA that were sufficiently abundant to be informative, the abundance of 1,027 of them varied between fractions: 36 of those were exclusive to one milk fraction. Feed restriction affected the abundance of 155 microRNA, with whole milk and milk extracellular vesicles being the most affected, whereas milk fat globules and exfoliated mammary epithelial cells were little or not affected at all. The high intensity feed restriction led to more microRNA variations in milk than moderate restriction. The target prediction of known microRNA that varied under feed restriction suggested the modification of some key pathways for lactation related to milk fat and protein metabolisms, cell cycle, and stress responses.

CONCLUSIONS

This study highlighted that the miRNome of each milk fraction is specific, with mostly the same microRNA composition but with variations in abundance between fractions. These specific miRNomes were affected differently by feed restrictions, the intensity of which appeared to be a major factor modulating milk miRNomes. These findings offer opportunities for future research on the use of milk miRNA as biomarkers of energy status in dairy cows, which is affected by feed restrictions.

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.

FishmiRNA: An evolutionarily supported microRNA annotation and expression database for ray-finned fishes

microRNAs are important post-transcriptional regulators of gene expression involved in countless biological processes and are widely studied across metazoans. While miRNA research continues to grow, the large community of fish miRNA researchers lacks exhaustive resources consistent among species. To fill this gap, we developed FishmiRNA, an evolutionarily supported microRNA annotation and expression database for ray-finned fishes: www.fishmirna.org. The self-explanatory database contains detailed, manually-curated miRNA annotations with orthology relationships rigorously established by sequence similarity and conserved syntenies, and expression data provided for each detected mature miRNA. In just few clicks, users can download the annotation and expression database in several convenient formats either in its entirety or a subset. Simple filters and BLAST search options also permit the simultaneous exploration and visual comparison of expression data for up to any ten mature miRNAs across species and organs. FishmiRNA was specifically designed for ease of use to reach a wide audience.

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.

Establishment of Repertoire of Placentome-Associated MicroRNAs and Their Appearance in Blood Plasma Could Identify Early Establishment of Pregnancy in Buffalo (Bubalus bubalis)

Precise early pregnancy diagnosis in dairy animals is of utmost importance for an efficient dairy production system. Not detecting a dairy animal pregnant sufficiently early after the breeding results to extending the unproductive time of their milk production cycle and causes substantial economic loss for a dairy producer. At present, the most conventional and authentic pregnancy confirmation practice in cows and buffaloes is rectal palpation of the reproductive organs at Days 35–40 after insemination, which sometime leads to considering an animal as false pregnant. Other alternative methods available for early pregnancy diagnosis lack either accuracy or reproducibility or require elaborate instrumentation and laboratory setup not feasible to practice at farmers’ doorstep. The present study was aimed at establishment of the microRNA (miRNA) repertoire of the placentome in buffaloes, which could capture the event of the cross talk between a growing embryo and a dam, through fetal cotyledons and maternal caruncles, and thus could hint at the early pregnancy establishment event in ruminants. Total RNA was isolated from buffalo placentome tissues during early stages of pregnancy (at Day < 25 and Days 30–35), and global small RNA analysis was performed by using Illumina single-end read chemistry and Bubalus bubalis genome. A total of 2,199 miRNAs comprising 1,620 conserved and 579 non-conserved miRNAs were identified. Stringent functional miRNA selection criteria could predict 20 miRNAs worth evaluating for their abundance in the plasma of pregnant, non-pregnant, cyclic non-bred, and non-cyclic prepubertal animals. Eight of them (viz., miR-195-5p, miR-708-3p, miR-379-5p, miR-XX1, miR-XX2, miR-130a-3p, miR-200a-3p, and miR-27) displayed typical abundance patterns in the plasma samples of the animals on Day 19 as well as Day 25 post-insemination, thus making them ambiguous candidates for early pregnancy detection. Similarly, higher abundance of miR-200a-3p and miR130a-3p in non-pregnant animals was indicative of their utility for detecting the animals as not pregnant. Most interestingly, miR-XX1 and miR-XX2 were very characteristically abundant only in pregnant animals. In silico target prediction analysis confirmed that these two miRNAs are important regulators of cyclooxygenase-2 (COX-2) and cell adhesion molecule-2 (CADM-2), both of which play a significant role in the implantation process during feto-maternal cross talk. We interpret that circulatory miR-XX1 and miR-XX2 in blood plasma could be the potential biomarkers for early pregnancy detection in buffaloes.

Proteomic profile of pre-implantational ovine embryos produced in vivo

The present study was conducted to decipher the proteome of in vivo-produced pre-implantation ovine embryos. Ten locally adapted Morana Nova ewes received hormonal treatment and were inseminated 12 hr after ovulation. Six days later, 54 embryos (morula and blastocyst developmental state) were recovered from eight ewes and pooled to obtain sufficient protein for proteomic analysis. Extracted embryo proteins were analysed by LC-MS/MS, followed by identification based on four database searches (PEAKS, Proteome Discoverer software, SearchGUI software, PepExplorer). Identified proteins were analysed for gene ontology terms, protein clusters and interactions. Genes associated with the ovine embryo proteome were screened for miRNA targets using data sets of TargetScan (http://www.targetscan.org) and mIRBase (http://www.mirbase.org) servers. There were 667 proteins identified in the ovine embryos. Biological processes of such proteins were mainly related to cellular process and regulation, and molecular functions, to binding and catalytic activity. Analysis of the embryo proteins revealed 49 enriched functional clusters, linked to energy metabolism (TCA cycle, pyruvate and glycolysis metabolism), zona pellucida (ZP), MAPK signalling pathway, tight junction, binding of sperm to ZP, translation, proteasome, cell cycle and calcium/phospholipid binding. Sixteen miRNAs were related to 25 pre-implantation ovine embryo genes, all conserved in human, bovine and ovine species. The interaction network generated by miRNet showed four key miRNAs (hsa-mir-106b-5p; hsa-mir-30-5p; hsa-mir-103a-5p and hsa-mir-106a-5p) with potential interactions with embryo-expressed genes. Functional analysis of the network indicated that miRNAs modulate genes related to cell cycle, regulation of stem cell and embryonic cell differentiation, among others. Retrieved miRNAs also modulate the expression of genes involved in cell signalling pathways, such as MAPK, Wnt, TGF-beta, p53 and Toll-like receptor. The current study describes the first major proteomic profile of 6-day-old ovine embryos produced in vivo, setting a comprehensive foundation for our understanding of embryo physiology in the ovine species.

Species-Specificity of Transcriptional Regulation and the Response to Lipopolysaccharide in Mammalian Macrophages

Mammalian macrophages differ in their basal gene expression profiles and response to the toll-like receptor 4 (TLR4) agonist, lipopolysaccharide (LPS). In human macrophages, LPS elicits a temporal cascade of transient gene expression including feed forward activators and feedback regulators that limit the response. Here we present a transcriptional network analysis of the response of sheep bone marrow-derived macrophages (BMDM) to LPS based upon RNA-seq at 0, 2, 4, 7, and 24 h post-stimulation. The analysis reveals a conserved transcription factor network with humans, and rapid induction of feedback regulators that constrain the response at every level. The gene expression profiles of sheep BMDM at 0 and 7 h post LPS addition were compared to similar data obtained from goat, cow, water buffalo, horse, pig, mouse and rat BMDM. This comparison was based upon identification of 8,200 genes annotated in all species and detected at >10TPM in at least one sample. Analysis of expression of transcription factors revealed a conserved transcriptional millieu associated with macrophage differentiation and LPS response. The largest co-expression clusters, including genes encoding cell surface receptors, endosome-lysosome components and secretory activity, were also expressed in all species and the combined dataset defines a macrophage functional transcriptome. All of the large animals differed from rodents in lacking inducible expression of genes involved in arginine metabolism and nitric oxide production. Instead, they expressed inducible transporters and enzymes of tryptophan and kynurenine metabolism. BMDM from all species expressed high levels of transcripts encoding transporters and enzymes involved in glutamine metabolism suggesting that glutamine is a major metabolic fuel. We identify and discuss transcripts that were uniquely expressed or regulated in rodents compared to large animals including ACOD1, CXC and CC chemokines, CD163, CLEC4E, CPM, CSF1, CSF2, CTSK, MARCO, MMP9, SLC2A3, SLC7A7, and SUCNR1. Conversely, the data confirm the conserved regulation of multiple transcripts for which there is limited functional data from mouse models and knockouts. The data provide a resource for functional annotation and interpretation of loci involved in susceptibility to infectious and inflammatory disease in humans and large animal species.

MicroRNA Milk Exosomes: From Cellular Regulator to Genomic Marker

Recent advances in ruminants' milk-derived exosomes (EXO) have indicated a role of microRNAs (miRNAs) in cell-to-cell communication in dairy ruminants. The miRNAs EXO retain peculiar mechanisms of uptake from recipient cells, which enables the selective delivery of cargos, with a specific regulation of target genes. Although many studies have been published on the miRNAs contained in milk, less information is available on the role of miRNAs EXO, which are considered stable over time and resistant to digestion and milk processing. Several miRNAs EXO have been implicated in the cellular signaling pathway, as in the regulation of immune response. Moreover, they exert epigenetic control, as extenuating the expression of DNA methyltransferase 1. However, the study of miRNAs EXO is still challenging due to the difficulty of isolating EXO. In fact, there are not agreed protocols, and different methods, often time-consuming, are used, making it difficult to routinely process a large number of samples. The regulation of cell functions in mammary glands by miRNAs EXO, and their applications as genomic markers in livestock, is presented.

Genome wide association analysis on semen volume and milk yield using different strategies of imputation to whole genome sequence in French dairy goats

Background

Goats were domesticated 10,500 years ago to supply humans with useful resources. Since then, specialized breeds that are adapted to their local environment have been developed and display specific genetic profiles. The VarGoats project is a 1000 genomes resequencing program designed to cover the genetic diversity of the Capra genus. In this study, our main objective was to assess the use of sequence data to detect genomic regions associated with traits of interest in French Alpine and Saanen breeds.

Results

Direct imputation from the GoatSNP50 BeadChip genotypes to sequence level was investigated in these breeds using FImpute and different reference panels: within-breed, all Capra hircus sequenced individuals, European goats and French mainland goats. The best results were obtained with the French goat panel with allele and genotype concordance rates reaching 0.86 and 0.75 in the Alpine and 0.86 and 0.73 in the Saanen breed respectively. Mean correlations tended to be low in both breeds due to the high proportion of variants with low frequencies.

For association analysis, imputation was performed using FImpute for 1129 French Alpine and Saanen males using within-breed and French panels on 23,338,436 filtered variants. The association results of both imputation scenarios were then compared. In Saanen goats, a large region on chromosome 19 was significantly linked to semen volume and milk yield in both scenarios. Significant variants for milk yield were annotated for 91 genes on chromosome 19 in Saanen goats. For semen volume, the annotated genes include YBOX2 which is related to azoospermia or oligospermia in other species. New signals for milk yield were detected on chromosome 2 in Alpine goats and on chromosome 5 in Saanen goats when using a multi-breed panel.

Conclusion

Even with very small reference populations, an acceptable imputation quality can be achieved in French dairy goats. GWAS on imputed sequences confirmed the existence of QTLs and identified new regions of interest in dairy goats. Adding identified candidates to a genotyping array and sequencing more individuals might corroborate the involvement of identified regions while removing potential imputation errors.

Characterization of Holstein and Normande whole milk miRNomes highlights breed specificities

The concept of milk as a healthy food has opened the way for studies on milk components, from nutrients to microRNAs, molecules with broad regulatory properties present in large quantities in milk. Characterization of these components has been performed in several species, such as humans and bovine, depending on the stages of lactation. Here, we have studied the variation in milk microRNA composition according to genetic background. Using high throughput sequencing, we have characterized and compared the milk miRNomes of Holstein and Normande cattle, dairy breeds with distinct milk production features, in order to highlight microRNAs that are essential for regulation of the lactation process. In Holstein and Normande milk, 2,038 and 2,030 microRNAs were identified, respectively, with 1,771 common microRNAs, of which 1,049 were annotated and 722 were predicted. The comparison of the milk miRNomes of two breeds allowed to highlight 182 microRNAs displaying significant differences in the abundance. They are involved in the regulation of lipid metabolism and mammary morphogenesis and development, which affects lactation. Our results provide new insights into the regulation of molecular mechanisms involved in milk production.