Characterization of miRNA profiles in the mammary tissue of dairy cattle in response to heat stress

Circ_002033 Regulates Proliferation, Apoptosis, and Oxidative Damage of Bovine Mammary Epithelial Cells via the miR-199a-5p-MAP3K11 Axis in Heat Stress

Heat stress is becoming the major factor regarding dairy cow health and milk quality because of global warming. Circular RNAs (circRNAs) represent a special type of noncoding RNAs, which are related to regulating many biological processes. Nonetheless, little is known concerning their effects on heat-stressed bovine mammary epithelial cells (BMECs). Here, this study found a novel circRNA, circ_002033, using RNA sequencing (RNA-seq) and explored the role and underlying regulatory mechanism in proliferation, apoptosis, and oxidative damage in a heat-stressed bovine mammary epithelial cell line (MAC-T). According to the previous RNA-seq analysis, the abundance of circ_002033 in mammary gland tissue of heat-stressed cows increased relative to nonheat-stressed counterparts. This study found that the knockdown of circ_002033 promoted proliferation and alleviated apoptosis and oxidative damage in heat-stressed MAC-T. Mechanistically, circ_002033 localizes to miR-199a-5p in the cytoplasm of MAC-T to regulate mitogen-activated protein kinase kinase 11 (MAP3K11) expression. Meanwhile, miR-199a-5p and MAP3K11 are also involved in regulating the proliferation and apoptosis of heat-stressed MAC-T. Importantly, circ_002033 knockdown promoted the expression of miR-199a-5p while decreasing that of MAP3K11, thereby enhancing proliferation while alleviating apoptosis and oxidative damage in heat-stressed MAC-T. In summary, we found that circ_002033 regulates the proliferation, apoptosis, and oxidative damage of heat-stressed BMECs through the miR-199a-5p/MAP3K11 axis, providing the theoretical molecular foundation for mitigating heat stress of dairy cows.

Identification of Hub Genes and Target miRNAs Crucial for Milk Production in Holstein Friesian Dairy Cattle

Dairy milk production is a quantitative trait that is controlled by many biological and environmental factors. This study employs a network-driven systems approach and clustering algorithm to uncover deeper insights into its genetic associations. We analyzed the GSE33680 dataset from the GEO database to understand the biological importance of milk production through gene expression and modules. In this study, we employed CytoNCA and ClusterONE plugins within Cytoscape for network analysis. Moreover, miRWalk software was utilized to detect miRNAs, and DAVID was employed to identify gene ontology and pathways. The results revealed 140 up-regulated genes and 312 down-regulated genes. In addition, we have identified 91 influential genes and 47 miRNAs that are closely associated with milk production. Through our examination of the network connecting these genes, we have found significant involvement in important biological processes such as calcium ion transit across cell membranes, the BMP signaling pathway, and the regulation of MAPK cascade. The conclusive network analysis further reveals that GAPDH, KDR, CSF1, PYGM, RET, PPP2CA, GUSB, and PRKCA are closely linked to key pathways essential for governing milk production. Various mechanisms can control these genes, making them valuable for breeding programs aiming to enhance selection indexes.

Potential role of miR-8159-x in heat stress response in rainbow trout (Oncorhynchus mykiss)

Rainbow trout (Oncorhynchus mykiss) is a representative species of cold-water fish. Elevated temperatures during summer often result in significant high mortality rates. MicroRNAs (miRNAs) are class of small non-coding RNAs that play a crucial role as post-transcriptional regulators in various biological processes. Emerging evidence suggests that miRNAs are important regulators role during heat stress. Analyzing previously obtained miRNA-sequencing data, we observed substantial down regulation of miR-8159-x in the liver tissue of heat stressed rainbow trout. In this study, we conducted a dual luciferase reporter assay to validate that miR-8159-x target, a key gene involved in heat stress in rainbow trout. By examining the expression patterns of miR-8159-x and hsp90a1 in the liver tissue at 18 °C (CG) and 24 °C (HS) groups, we propose that miR-8159-x may negatively regulate hsp90a1. Furthermore, in vitro hepatocyte assay, transfection with miR-8159-x mimics significantly reduced the expression level of hsp90a1, whereas transfection with a miR-8159-x inhibitor yielded the opposite effect. Additionally, overexpression of miR-8159-x inhibited cell proliferation and induced apoptosis in normal rainbow trout hepatocytes. We further investigated the effects of miR-8159-x overexpression or inhibition on the mRNA and protein levels of the target gene hsp90a1 under heat stress conditions. In conclusion, our findings suggest that miR-8159-x participates in the biological response to heat stress by targeting hsp90a1. These results contribute to a better understanding of the molecular mechanisms underlying heat stress in rainbow trout and provide valuable insights for future research.

Impact of heat stress on dairy cattle and selection strategies for thermotolerance: a review

Climate change is a problem that causes many environmental issues that impact the productivity of livestock species. One of the major issues associated with climate change is an increase of the frequency of hot days and heat waves, which increases the risk of heat stress for livestock species. Dairy cattle have been identified as being susceptible to heat stress due to their high metabolic heat load. Studies have shown heat stress impacts several biological processes that can result in large economic consequences. When heat stress occurs, dairy cattle employ several physiological and cellular mechanisms in order to dissipate heat and protect cells from damage. These mechanisms require an increase and diversion in energy toward protection and away from other biological processes. Therefore, in turn heat stress in dairy cattle can lead numerous issues including reductions in milk production and reproduction as well as increased risk for disease and mortality. This indicates a need to select dairy cattle that would be thermotolerant. Various selection strategies to confer thermotolerance have been discussed in the literature, including selecting for reduced milk production, crossbreeding with thermotolerant breeds, selecting based on physiological traits and most recently selecting for enhanced immune response. This review discusses the various issues associated with heat stress in dairy cattle and the pros and cons to the various selection strategies that have been proposed to select for thermotolerance in dairy cattle.

Extracellular vesicle-microRNAs mediated response of bovine ovaries to seasonal environmental changes

BACKGROUND

Among the various seasonal environmental changes, elevated ambient temperature during the summer season is a main cause of stress in dairy and beef cows, leading to impaired reproductive function and fertility. Follicular fluid extracellular vesicles (FF-EVs) play an important role in intrafollicular cellular communication by, in part, mediating the deleterious effects of heat stress (HS). Here we aimed to investigate the changes in FF-EV miRNA cargoes in beef cows in response to seasonal changes: summer (SUM) compared to the winter (WIN) season using high throughput sequencing of FF-EV-coupled miRNAs. In addition to their biological relevance, the potential mechanisms involved in the packaging and release of those miRNAs as a response to environmental HS were elucidated.

RESULTS

Sequencing analysis revealed that an average of 6.6% of the EV-RNA mapped reads were annotated to bovine miRNAs. Interestingly, miR-148a, miR-99a-5p, miR-10b, and miR-143 were the top four miRNAs in both groups accounting for approximately 52 and 62% of the total miRNA sequence reads in the SUM and WIN groups, respectively. A group of 16 miRNAs was up-regulated and 8 miRNAs were down-regulated in the SUM compared to the WIN group. Five DE-miRNAs (miR-10a, miR-10b, miR-26a, let-7f, and miR-1246) were among the top 20 expressed miRNA lists. Sequence motif analysis revealed the appearance of two specific motifs in 13 out of the 16 upregulated miRNAs under HS conditions. Both motifs were found to be potentially bonded by specific RNA binding proteins including Y-box binding proteins (YBX1 and YBX2) and RBM42.

CONCLUSION

Our findings indicate that FF EV-coupled miRNA profile varies under seasonal changes. These miRNAs could be a good indicator of the cellular mechanism in mediating HS response and the potential interplay between miRNA motifs and RNA binding proteins can be one of the mechanisms governing the packaging and release of miRNAs via EVs to facilitate cellular survival.

The alteration of N6-methyladenosine (m6A) modification at the transcriptome-wide level in response of heat stress in bovine mammary epithelial cells

Background

Heat stress has a substantial negative economic impact on the dairy industry. N6-methyladenosine (m6A) is the most common internal RNA modification in eukaryotes and plays a key role in regulating heat stress response in animals. In dairy cows, however, this modification remains largely unexplored. Therefore, we examined the effects of heat stress on the m6A modification and gene expression in bovine mammary epithelial cells to elucidate the mechanism of heat stress response. In this study, Mammary alveolar cells-large T antigen (MAC-T) cells were incubated at 37 °C (non-heat stress group, NH) and 40 °C (heat stress group, H) for 2 hours, respectively. HSP70, HSF1, BAX and CASP3 were up regulated in H group compared with those in the NH group.

Results

Methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) were conducted to identify m6A peaks and to produce gene expression data of MAC-T cells in the H and NH groups. In total, we identified 17,927 m6A peaks within 9355 genes in the H group, and 18,974 peaks within 9660 genes in the NH groups using MeRIP-seq. Compared with the NH group, 3005 significantly differentially enriched m6A peaks were identified, among which 1131 were up-regulated and 1874 were down-regulated. In addition, 1502 significantly differentially expressed genes were identified using RNA-seq, among which 796 were up-regulated and 706 were down-regulated in the H group compared to the NH group. Furthermore, 199 differentially expressed and synchronously differentially methylated genes were identified by conjoint analysis of the MeRIP-seq and RNA-seq data, which were subsequently divided into four groups: 47 hyper-up, 53 hyper-down, 59 hypo-up and 40 hypo-down genes. In addition, GO enrichment and KEGG analyses were used to analyzed the potential functions of the genes in each section.

Conclusion

The comparisons of m6A modification patterns and conjoint analyses of m6A modification and gene expression profiles suggest that m6A modification plays a critical role in the heat stress response by regulating gene expression.

Programming effects of late gestation heat stress in dairy cattle

The final weeks of gestation represent a critical period for dairy cows that can determine the success of the subsequent lactation. Many physiological changes take place and additional exogenous stressors can alter the success of the transition into lactation. Moreover, this phase is pivotal for the final stage of intrauterine development of the fetus, which can have negative long-lasting postnatal effects. Heat stress is widely recognised as a threat to dairy cattle welfare, health, and productivity. Specifically, late gestation heat stress impairs the dam's productivity by undermining mammary gland remodelling during the dry period and altering metabolic and immune responses in early lactation. Heat stress also affects placental development and function, with relevant consequences on fetal development and programming. In utero heat stressed newborns have reduced birth weight, growth, and compromised passive immune transfer. Moreover, the liver and mammary DNA of in utero heat stressed calves show a clear divergence in the pattern of methylation relative to that of in utero cooled calves. These alterations in gene regulation might result in depressed immune function, as well as altered thermoregulation, hepatic metabolism, and mammary development jeopardising their survival in the herd and productivity. Furthermore, late gestation heat stress appears to exert multigenerational effects, influencing milk yield and survival up to the third generation.

Negative effects of heat stress on ovarian tissue in female rabbit

Numerous studies have highlighted the role of miRNA in the deformation and necrosis of cells of ovarian tissue caused by heat stress (HS), which ultimately affects ovarian function. Although the role of small RNAs has been investigated in alterations in ovarian tissue functioning in response to HS, the expression profile of ovarian miRNA has been explored to a lesser extent. In this study, female rabbits were subject to HS treatment by using electrical heater. The current work demonstrated that HS could significantly change physiological performance of female rabbits including body weight, rectal temperature and relative ovary weight, and significantly reduce serum IL-2, IL-8, CAT, and GSH-Px concentrations by enzyme-linked immunosorbent assay (ELISA) technique. As a result, an increase in apoptosis in ovarian cells, as well as unhealthy follicles, were observed by Hematoxylin-eosin (HE) and TUNEL staining. Additionally, small RNA-seq revealed changes in the miRNA expression profile of rabbit ovaries under HS. Five hundred fourteen miRNAs were obtained including known miRNAs 442 and novel miRNAs 72. Among these miRNAs, 23 miRNAs were significantly expressed under HS. Eleven differentially expressed miRNAs (DE miRNAs) and 9 their predicted targets were confirmed by qPCR, which were expected miRNA-mRNA negative regulation pattern. Among the DE miRNAs and targets, miR-141-39 may target COQ6, miR-449a-5p and miR-34c-5p may control RFC5 and RTN2 together, miR-449a-5p may target ACADVL, miR-34c-5p potentially targets Bcl-2 and miR-196b-5p potentially regulates CASK and HOXB6. Thus, the current work suggested the negative effects of HS on the ovarian tissue of female rabbits, and in conclusion these changes could be caused by decreased serum IL-2, IL-8, CAT and GSH-Px levels, increased ovarian apoptosis, and changed the expression of miRNAs.

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.

Assessment of microRNA profiles in small extracellular vesicles isolated from bovine colostrum with different immunoglobulin G concentrations

•

The RNA concentration in sEV isolated from bovine colostrum was highest using the combination of the miRCURY Cell/Urine/CSF and miRNeasy Mini kits.

•

The top 50 miRNA were the same using miRDeep2 and sRNAbench, predominated by let-7b, let-7a-5p, miR-30a-5p, and miR-148a.

•

Predicted target genes of the top 50 miRNA regulate PI3K-Akt and MAPK signaling pathways, axon guidance, and focal adhesion.

•

The abundance of miR-27a-3p was higher in colostrum with high IgG concentrations.

The consumption of bovine colostrum by newborn calves during the first days of life is essential to ensure the transfer of passive immunity. In addition to critical IgG, colostrum also contains non-IgG biomolecules, including microRNA (miRNA). The present study investigated the profiles of miRNA in small extracellular vesicles (sEV) isolated from bovine colostrum with high (256.5 ± 5.7 mg/mL, mean ± standard deviation, n = 4) and low (62.8 ± 3.6 mg/mL, n = 4) concentrations of IgG. Different combination of sEV extraction methods and bioinformatic pipelines (miRDeep2 and sRNAbench) for miRNA analysis were evaluated. Results showed that miRCURY exosome Cell/Urine/CSF and miRNeasy Mini kits yielded the highest RNA concentration. The miRNA-seq data analysis showed miRDeep2 yielded more comprehensive miRNAome compared with sRNAbench (527 versus 392 unique miRNA), whereas 389 shared miRNA were identified using both approaches. The profiles of top 50 miRNA were the same using both approaches, and their abundance contributed to 91.7% and 94.3% of total abundance of miRNA using miRDeep2 and sRNAbennch, respectively. These core miRNA were predicted to target 2,655 genes, which regulate 78 KEGG (Kyoto Encyclopedia of Genes and Genomes) level-3 pathways including PI3K-Akt and MAPK signaling pathway, axon guidance, and focal adhesion. The expression profiles of sEV-associated miRNA were similar between high- and low-IgG colostrum samples, despite the fact that the abundance of miR-27a-3p was higher in colostrum with high concentrations of IgG. In conclusion, a core miRNAome in bovine colostrum may play a role in regulating health and developmental stages in neonatal calves, independent of IgG concentration.

Investigating circulating miRNA in transition dairy cows: What miRNAomics tells about metabolic adaptation

In the current study, we investigated dairy cows’ circulating microRNA (miRNA) expression signature during several key time points around calving, to get insights into different aspects of metabolic adaptation. In a trial with 32 dairy cows, plasma samples were collected on days −21, 1, 28, and 63 relative to calving. Individually extracted total RNA was subjected to RNA sequencing using NovaSeq 6,000 (Illumina, CA) on the respective platform of IGA Technology Services, Udine, Italy. MiRDeep2 was used to identify known and novel miRNA according to the miRbase collection. Differentially expressed miRNA (DEM) were assessed at a threshold of fold-change > 1.5 and false discovery rate < 0.05 using the edgeR package. The MiRWalk database was used to predict DEM targets and their associated KEGG pathways. Among a total of 1,692 identified miRNA, 445 known miRNA were included for statistical analysis, of which 84, 59, and 61 DEM were found between days −21 to 1, 1 to 28, and 28 to 63, respectively. These miRNA were annotated to KEGG pathways targeting the insulin, MAPK, Ras, Wnt, Hippo, sphingolipid, T cell receptor, and mTOR signaling pathways. MiRNA-mRNA network analysis identified miRNA as master regulators of the biological process including miR-138, miR-149-5p, miR-2466-3p, miR-214, miR-504, and miR-6523a. This study provided new insights into the miRNA signatures of transition to the lactation period. Calving emerged as a critical time point when miRNA were most affected, while the following period appeared to be recovering from massive parturition changes. The primarily affected pathways were key signaling pathways related to establishing metabolic and immune adaptations.

Weighted single-step GWAS and RNA sequencing reveals key candidate genes associated with physiological indicators of heat stress in Holstein cattle

Background

The study of molecular processes regulating heat stress response in dairy cattle is paramount for developing mitigation strategies to improve heat tolerance and animal welfare. Therefore, we aimed to identify quantitative trait loci (QTL) regions associated with three physiological indicators of heat stress response in Holstein cattle, including rectal temperature (RT), respiration rate score (RS), and drooling score (DS). We estimated genetic parameters for all three traits. Subsequently, a weighted single-step genome-wide association study (WssGWAS) was performed based on 3200 genotypes, 151,486 phenotypic records, and 38,101 animals in the pedigree file. The candidate genes located within the identified QTL regions were further investigated through RNA sequencing (RNA-seq) analyses of blood samples for four cows collected in April (non-heat stress group) and four cows collected in July (heat stress group).

Results

The heritability estimates for RT, RS, and DS were 0.06, 0.04, and 0.03, respectively. Fourteen, 19, and 20 genomic regions explained 2.94%, 3.74%, and 4.01% of the total additive genetic variance of RT, RS, and DS, respectively. Most of these genomic regions are located in the Bos taurus autosome (BTA) BTA3, BTA6, BTA8, BTA12, BTA14, BTA21, and BTA24. No genomic regions overlapped between the three indicators of heat stress, indicating the polygenic nature of heat tolerance and the complementary mechanisms involved in heat stress response. For the RNA-seq analyses, 2627 genes were significantly upregulated and 369 downregulated in the heat stress group in comparison to the control group. When integrating the WssGWAS, RNA-seq results, and existing literature, the key candidate genes associated with physiological indicators of heat stress in Holstein cattle are: PMAIP1, SBK1, TMEM33, GATB, CHORDC1, RTN4IP1, and BTBD7.

Conclusions

Physiological indicators of heat stress are heritable and can be improved through direct selection. Fifty-three QTL regions associated with heat stress indicators confirm the polygenic nature and complex genetic determinism of heat tolerance in dairy cattle. The identified candidate genes will contribute for optimizing genomic evaluation models by assigning higher weights to genetic markers located in these regions as well as to the design of SNP panels containing polymorphisms located within these candidate genes.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00748-6.

A Survey of Genome-Wide Genetic Characterizations of Crossbred Dairy Cattle in Local Farms in Cambodia

Simple Summary

The Cambodian dairy sector collapsed due to the civil war. Thereby, the Cambodia dairy sector has not obtained basic information, such as cattle breed and milk production. To encourage the recovery of the Cambodian dairy sector, the present study aimed to reveal the genetic variation and the milk production in Cambodian crossbred dairy cattle. Initially, we conducted interviews to understand the breeding background and milk production of two local farms (Farm R and Farm M) in Cambodia. The percentage (%) of milk fat content in Farm R was higher than that in Farm M. A genome-wide analysis of the genetic characterization for 75 cows in the two dairy farms implies that some cows in Farm R are genetically far from other crossbred cattle and retain a higher proportion of genetic background derived from Cambodian-native cattle. The present study indicated genetic characteristics and milk composition in Cambodian crossbred cattle. Gene character in Cambodian local crossbred cattle could contribute to milk production in the Cambodian dairy system. Thus, our basic genetic study could provide a new breeding strategy by using local Cambodian crossbred dairy cattle to establish an adequate dairy strain.

Abstract

To improve the dairy sector in Cambodia in the future, we aimed to reveal the genetic variation and the milk production in Cambodian crossbred dairy cattle. We calculated the percent (%) milk fat content and the average milk yield per cow (L/day) for two farms (Farm R and M) based on the farmers’ records and interviews. The crossbred cows originated from Cambodian local farmers and Thailand breeders in Farm R, whereas the crossbred cows originated in Thailand breeders in Farm M. Then, we performed genetic characterization for 75 individuals from the two farms and an individual Japanese pure Holstein-Friesian cow based on 133,705 single nucleotide polymorphisms (SNPs) obtained by the GRAS-Di method. The milk fat contents in the bulk milk in the dry season and the average milk yield per cow on Farm R were 3.77 ± 0.98% and 7.81 ± 2.66 L/day, respectively, and were higher than those on Farm M (3.35 ± 0.54% and 6.5–7.5 L/day). Cattle originating in Cambodia in Farm R possessed a unique genetic character different from cattle from Thailand in Farm M. The present study suggests that the differences in milk fat content between the two farms might be explained by the genetic differences in crossbred cows.

Fecal microbiota and their association with heat stress in Bos taurus

Background

Humans have been influencing climate changes by burning fossil fuels, farming livestock, and cutting down rainforests, which has led to global temperature rise. This problem of global warming affects animals by causing heat stress, which negatively affects their health, biological functions, and reproduction. On the molecular level, it has been proved that heat stress changes the expression level of genes and therefore causes changes in proteome and metabolome. The importance of a microbiome in many studies showed that it is considered as individuals’ “second genome”. Physiological changes caused by heat stress may impact the microbiome composition.

Results

In this study, we identified fecal microbiota associated with heat stress that was quantified by three metrics – rectal temperature, drooling, and respiratory scores represented by their Estimated Breeding Values. We analyzed the microbiota from 136 fecal samples of Chinese Holstein cows through a 16S rRNA gene sequencing approach. Statistical modeling was performed using a negative binomial regression. The analysis revealed the total number of 24 genera and 12 phyla associated with heat stress metrics. Rhizobium and Pseudobutyrivibrio turned out to be the most significant genera, while Acidobacteria and Gemmatimonadetes were the most significant phyla. Phylogenetic analysis revealed that three heat stress indicators quantify different metabolic ways of animals’ reaction to heat stress. Other studies already identified that those genera had significantly increased abundance in mice exposed to stressor-induced changes.

Conclusions

This study provides insights into the analysis of microbiome composition in cattle using heat stress measured as a continuous variable. The bacteria highly associated with heat stress were highlighted and can be used as biomarkers in further microbiological studies.

Transcriptome and iTRAQ-Based Proteome Reveal the Molecular Mechanism of Intestinal Injury Induced by Weaning Ewe's Milk in Lambs

Early feeding regime has a substantial lifelong effect on lambs and weaning ewe's milk can lead to the intestinal injury of lambs. To explore the molecular regulatory mechanism of intestinal injury of lambs under weaning stress, the jejunum was conducted transcriptome and then integrated analyzed with our previous proteome data. A total of 255 upregulated genes and 285 downregulated genes were significantly identified. These genes showed low overlapping with differentially expressed proteins identified by isobaric tags for relative and absolute quantification (iTRAQ). However, according to their functions, the differentially expressed genes (DEGs) and proteins with the same expression trend were enriched for the similar Gene Ontology (GO) terms and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, such as intestinal lipid absorption, urea cycle, peroxisome proliferator-activated receptor (PPAR) signaling pathway, and ferroptosis. Furthermore, the DEGs, including FABP2, ACSL3, APOA2, APOC3, and PCK1, might play essential roles in intestinal lipid absorption and immune response through the PPAR signaling pathway and ferroptosis. This study could provide new insights into early lamb breeding at the molecular level.

Differential Expression of MicroRNAs in Dark-Cutting Meat from Beef Carcasses

“Dark-cutting” meat in beef carcasses can result from conditions such as long-term stress and depleted glycogen stores, but some aspects of the physiological mechanisms that cause dark-cutting phenotypes remain poorly understood. Certain responses to stress factors in fully developed tissues are known to be regulated by specific microRNAs. We investigated microRNA expression in Longissimus lumborum biopsies from carcasses derived from a contemporary group of 78 steers from which a high incidence of dark-cutting meat occurred. Our objective was to identify any potential microRNA signatures that reflect the impact of environmental factors and stresses on genetic signaling networks and result in dark-cutting beef (also known as dark, firm, and dry, or DFD) in some animals. MicroRNA expression was quantified by Illumina NextSeq small RNA sequencing. When RNA extracts from DFD muscle biopsy samples were compared with normal, non-DFD (NON) samples, 29 differentially expressed microRNAs were identified in which expression was at least 20% different in the DFD samples (DFD/NON fold ratio ≤0.8 or ≥1.2). When correction for multiple testing was applied, a single microRNA bta-miR-2422 was identified at a false discovery probability (FDR) of 5.4%. If FDR was relaxed to 30%, additional microRNAs were differentially expressed (bta-miR-10174-5p, bta-miR-1260b, bta-miR-144, bta-miR-142-5p, bta-miR-2285at, bta-miR-2285e, bta-miR-3613a). These microRNAs may play a role in regulating aspects of stress responses that ultimately result in dark-cutting beef carcasses.

Analysis of Genomic Alternative Splicing Patterns in Rat under Heat Stress Based on RNA-Seq Data

Heat stress is one of the most severe challenges faced in livestock production in summer. Alternative splicing as an important post-transcriptional regulation is rarely studied in heat-stressed animals. Here, we performed and analyzed RNA-sequencing assays on the liver of Sprague-Dawley rats in control (22 °C, n = 5) and heat stress (4 °C for 120 min, H120; n = 5) groups, resulting in the identification of 636 differentially expressed genes. Identification analysis of the alternative splicing events revealed that heat stress-induced alternative splicing events increased by 20.18%. Compared with other types of alternative splicing events, the alternative start increased the most (43.40%) after heat stress. Twenty-eight genes were differentially alternatively spliced (DAS) between the control and H120 groups, among which Acly, Hnrnpd and mir3064 were also differentially expressed. For DAS genes, Srebf1, Shc1, Srsf5 and Ensa were associated with insulin, while Cast, Srebf1, Tmem33, Tor1aip2, Slc39a7 and Sqstm1 were enriched in the composition of the endoplasmic reticulum. In summary, our study conducts a comprehensive profile of alternative splicing in heat-stressed rats, indicating that alternative splicing is one of the molecular mechanisms of heat stress response in mammals and providing reference data for research on heat tolerance in mammalian livestock.

Effects of nutrient restriction during early or mid-gestation in bovine on placental development and miRNA expression in the cotyledon

The objective of this study was to determine effects of maternal nutrient restriction (NR) during early or mid-gestation on uterine composition and miRNA expression in cotyledons. Primiparous Angus-cross cows (n = 38) were synchronized and inseminated using male sexed semen, blocked by body condition score and body weight (BW), and assigned to treatments. Animals were fed either: control (CON; gain 1 kg/week) or NR (55% maintenance energy and crude protein requirements) based on BW. An initial set of animals were fed either NR (n = 8) or CON (n = 8) from day 30-110 of gestation. A second set of animals were fed CON (n = 8) d 30-190 (CON/CON); NR (n = 7) day 30-110 followed by CON day 110-190 (NR/CON); or CON (n = 7) day 30-110 followed by NR day 110-190 (CON/NR). Cows were harvested on day 110 or 190 of gestation to collect placental tissues. RNA was isolated from cotyledon samples (3 animals/group) prior to microarray analysis using known Bos taurus microRNA sequences. Relative microRNA abundance was analyzed via ANOVA. Maternal NR increased (P < 0.05) cotyledon weight and total placentome surface area irrespective of gestational day. At day 110 of gestation, 51 microRNAs were reduced while 91 microRNAs observed greater abundance (P < 0.05) in NR verses CON cotyledons. At day 190 of gestation, 40 microRNAs were reduced and 26 microRNAs were increased (P < 0.05) in both NR/CON and CON/NR verses CON cotyledons. Top KEGG pathway analysis included: axon guidance, endocytosis, neuroactive ligand receptor interaction, and MAPK signaling pathway. Early-gestation maternal NR altered microRNA abundance to a greater extent than mid-gestation NR.

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.

Non-coding RNA in raw and commercially processed milk and putative targets related to growth and immune-response

BACKGROUND

Bovine milk contains extracellular vesicles (EVs) that play a role in cellular communication, acting in either an autocrine, paracrine, or an exocrine manner. The unique properties of the EVs protect the cargo against degradation. We profiled the ncRNAs (non-coding RNA) present in the EVs from seven dairy products - raw whole milk, heat-treated skim milk, homogenized heat-treated skim milk, pasteurized homogenized skim milk, pasteurized heavy whipping cream, sweet cream buttermilk and cultured buttermilk with four replicates each, obtained at different processing steps from a commercial dairy plant. EVs and their cargo were extracted by using a validated commercial kit that has been shown to be efficient and specific for EVs. Further, to find the annotation of ncRNAs, we probed bovine and other organism repositories(such as miRBase, miRTarBase, Ensemble) to find homolog ncRNA annotation in case the annotations of ncRNA are not available in Bos Taurus database.

RESULTS

Specifically, 30 microRNAs (miRNAs), were isolated throughout all the seven milk samples, which later when annotated with their corresponding 1546 putative gene targets have functions associated with immune response and growth and development. This indicates the potential for these ncRNAs to beneficially support mammary health and growth for the cow as well as neonatal gut maturation. The most abundant miRNAs were bta-miR-125a and human homolog miR-718 based on the abundance values of read count obtained from the milk samples.bta-miR-125a is involved in host bacterial and viral immune response, and human homolog miR-718 is involved in the regulation of p53, VEGF, and IGF signaling pathways, respectively. Sixty-two miRNAs were up-regulated and 121 miRNAs were down-regulated throughout all the milk samples when compared to raw whole milk. In addition, our study explored the putative roles of other ncRNAs which included 88 piRNAs (piwi-interacting RNA), 64 antisense RNAs, and 105 lincRNAs (long-intergenic ncRNAs) contained in the bovine exosomes.

CONCLUSION

Together, the results indicate that bovine milk contains significant numbers of ncRNAs with putative regulatory targets associated with immune- and developmental-functions important for neonatal bovine health, and that processing significantly affects the ncRNA expression values; but statistical testing of overall abundance(read counts) of all miRNA samples suggests abundance values aren't much affected. This can be attributed to the breakage of exosomal vesicles during the processing stages. It is worth noting, however, that these gene regulatory targets are putative, and further evidence could be generated through experimental validation.

Acute heat stress reduces viability but increases lactate secretion of porcine immature Sertoli cells through transcriptome reprogramming

Sertoli cells, important constituents of the somatic niche, supports the growth and development of spermatogonia. Heat stress (HS), among multiple intrinsic and external factors, can induce physiological and biochemical changes in Sertoli cells. However, the underlying molecular mechanism remains largely unclear. Here, we showed that acute heat stress (43 °C, 0.5 h) could reduce cell viability, promote apoptosis, and increase the lactate production of porcine immature Sertoli cells (iSCs) cultured in vitro. Then, transcriptome sequencing identified 126 immediately and 3372 prolonged responded differentially expressed genes (DEGs) after acute heat stress (43 °C, 0.5 h) (HS0.5), and 36 h recovery culture following heat stress (HS0.5-R36), respectively. Enrichment analyses found different signaling pathways: immediate changes including cell response to heat, regulation of cellular response to stress, heat shock protein binding, chaperon-mediated protein folding, and sterol biosynthetic process, but prolonged changes mainly involving cell cycle, regulation of apoptotic process/cell proliferation, reproductive process, P53, PI3K-Akt and Glycolysis/Gluconeogenesis. Furthermore, transcriptional patterns of 9 DEGs (Dnajb1, Traf6, Insig1, Gadd45g, Hdac6, Fkbp4, Serpine1, Pfkp and Galm), and 6 heat shock proteins (HSPs) (Hspa6, Hspb1, Hspd1, HSP90aa1, HSP90ab1 and Hsph1) were validated, as well as the protein pattern of HSP90AA1 via immunostaining and western blot. Taken together, heat stress could initiate immediate changes of heat shock-related genes, and reprogram transcriptome and signaling pathways affecting the viability, apoptosis and metabolite production of pig iSCs.

Can miRNAs Serve as Potential Markers in Thermal Burn Injury: An In Silico Approach

Burn injury has been a major cause of morbidity at global levels. They can occur by multiple agents, such as thermal radiation and chemicals. Among all, thermal burn is predominant and may require specialized treatment in some patients. Although various biomarkers are reportedly used in thermal burn for understanding the pathophysiology of the injury, their limitations prompt for the search of suitable markers that can address the depth and severity of the burn. MicroRNAs (miRNAs) are conserved noncoding molecules that seem to be the promising marker due to their role in multiple pathways and participation in different physiological processes of the body. The present review highlights the role of miRNAs in the repair of the wound and their interaction with specific genes in response to burn stress. Key miR candidates include miR-21, miR-29a, miR-378a-5p, miR-100, miR-27b, miR-200c, miR-150, miR-499-5p, miR-92a, miR-194, and miR-146b, which are identified for their respective targets involved in wound repair. Furthermore, bioinformatics and computational tools were used to confirm the miRNAs and their specific targets. Gene and miRNA expression data sets were downloaded from Research Collaboratory for Structural Bioinformatics Protein Data Bank Database and RNAComposer, respectively, and docked by PatchDock. The possible implications of the identified miRNAs could be in understanding the mechanism of burn injury. These can also be studied with the available drugs being used for burn injury. Apart from that, new intended molecules may also be tested for their effect on these miRNAs.

Genome-wide detection of copy number variations in Tharparkar cattle

Copy number variations (CNVs) are major forms of genetic variation with an increasing importance in animal genomics. This study used the Illumina BovineSNP 50 K BeadChip to detect the genome-wide CNVs in the Tharparkar cattle. With the aid of PennCNV software, we noticed a total of 447 copy number variation regions (CNVRs) across the autosomal genome, occupying nearly 2.17% of the bovine genome. The average size of detected CNVRs was found to be 122.2 kb, the smallest CNVR being 50.02 kb in size, to the largest being 1,232.87 Kb. Enrichment analyses of the genes in these CNVRs gave significant associations with molecular adaptation-related Gene Ontology (GO) terms. Most CNVR genes were significantly enriched for specific biological functions; signaling pathways, sensory responses to stimuli, and various cellular processes. In addition, QTL analysis of CNVRs described them to be linked with economically essential traits in cattle. The findings here provide crucial information for constructing a more comprehensive CNVR map for the indigenous cattle genome.

Heat shock during in vitro maturation of bovine oocytes disturbs bta-miR-19b and DROSHA transcripts abundance after in vitro fertilization

While microRNAs (miRNAs) are a class of non-coding RNAs important for embryo development, the relationship between them and heat stress during oocyte maturation has not yet been established. This study investigated the effect of heat shock during in vitro maturation (IVM) on the abundance of bta-miR-20a, -27b, -103, -21-5p, -19b, -1246 miRNAs and DROSHA and DICER1 mRNAs, previously reported for being involved in oocyte maturation, response to heat stress and miRNA biogenesis. Oocytes were exposed for 12h to heat shock during IVM, fertilized in vitro and the presumptive zygotes cultured for eight days. The relative quantification of miRNAs and mRNAs was performed by real-time PCR in vitro-matured oocytes and 8-cell stage embryos. Progression of meiosis, embryonic development and apoptotic indices was also evaluated. Heat shock compromised (p < .05) oocyte nuclear maturation, cleavage and embryo development, with a higher (p < .05) embryonic apoptotic index than the control group. The abundance of bta-miR-19b increased (p < .05) whereas the abundance of DROSHA transcripts decreased (p < .05) in embryos derived from heat-shocked oocytes. In conclusion, heat shock during IVM influences the abundance of bta-miR-19b and DROSHA in pre-implantation embryos, indicating a persistent effect of heat shock that can be associated with impaired embryo development.

Novel and known miRNAs in zebu (Tharparkar) and crossbred (Karan-Fries) cattle under heat stress

MicroRNAs (miRNAs) are small single-stranded non-coding RNAs that act as the master regulator of animal growth and development. RNA-RNA interaction is an important mechanism of gene regulation during biotic and abiotic stress. Heat stress (HS) is one of the most important abiotic factors which affect the growth, milk yield and reproductive health of the dairy animals. In the present investigation, we identified 387 known and 77 novel miRNAs from Tharparkar (TH) and Karan Fries (KF) cattle under HS condition. Family distribution analysis showed the identified miRNAs belong to more than 15 different families in which miR-2284 was the most abundant. We identified 42,350 targets for the known miRNAs reported in cattle. Pathway analysis of the identified targets showed most of the target genes were involved in cancer, mitogen-activated protein kinase (MAPK) signaling, calcium signaling, Ras signaling, and cAMP signaling pathways. Differential gene expression showed more than 344 miRNAs changed their expression significantly between control and HS condition. Heat map was generated for the top 20 most up and downregulated miRNAs. Ten miRNAs were validated using qRT-PCR to be heat responsive, based on read count value and differential gene expression. These novel miRNAs are new addition to the miRNA database of cattle. This study provides an overview of miRNA profile and their interaction with the target genes which leads to further understanding in deciphering the thermotolerance mechanism in cattle.

Elucidation of the Host Bronchial Lymph Node miRNA Transcriptome Response to Bovine Respiratory Syncytial Virus

Bovine respiratory disease (BRD) causes substantial morbidity and mortality, affecting cattle of all ages. One of the main causes of BRD is an initial inflammatory response to bovine respiratory syncytial virus (BRSV). MicroRNAs are novel and emerging non-coding small RNAs that regulate many biological processes and are implicated in various inflammatory diseases. The objective of the present study was to elucidate the changes in the bovine bronchial lymph node miRNA transcriptome in response to BRSV following an experimental viral challenge. Holstein-Friesian calves were either administered a challenge dose of BRSV (103.5 TCID50/ml × 15 ml) (n = 12) or were mock inoculated with sterile phosphate buffered saline (n = 6). Daily scoring of clinical signs was performed and calves were euthanized at day 7 post-challenge. Bronchial lymph nodes were collected for subsequent RNA extraction and sequencing (75 bp). Read counts for known miRNAs were generated using the miRDeep2 package using the UMD3.1 reference genome and the bovine mature miRNA sequences from the miRBase database (release 22). EdgeR was used for differential expression analysis and Targetscan was used to identify target genes for the differentially expressed (DE) miRNAs. Target genes were examined for enriched pathways and gene ontologies using Ingenuity Pathway Analysis (Qiagen). Multi-dimensional scaling (MDS) based on miRNA gene expression changes, revealed a clearly defined separation between the BRSV challenged and control calves, although the clinical manifestation of disease was only mild. One hundred and nineteen DE miRNAs (P < 0.05, FDR < 0.1, fold change > 1.5) were detected between the BRSV challenged and control calves. The DE miRNAs were predicted to target 465 genes which were previously found to be DE in bronchial lymph node tissue, between these BRSV challenged and control calves. Of the DE predicted target genes, 455 had fold changes that were inverse to the corresponding DE miRNAs. There were eight enriched pathways among the DE predicted target genes with inverse fold changes to their corresponding DE miRNA including: granulocyte and agranulocyte adhesion and diapedesis, interferon signalling and role of pathogen recognition receptors in recognition of bacteria and viruses. Functions predicted to be increased included: T cell response, apoptosis of leukocytes, immune response of cells and stimulation of cells. Pathogen recognition and proliferation of cytotoxic T cells are vital for the recognition of the virus and its subsequent elimination.

Characterization of Short-Term Heat Stress in Holstein Dairy Cows Using Altered Indicators of Metabolomics, Blood Parameters, Milk MicroRNA-216 and Characteristics

Simple Summary

In this study, we characterize the influence of short-term (4 days) heat stress on Holstein cows during early lactation. The use of indicators, such as production performance, physiological variables, blood parameters, micro-RNA expression, and metabolomes, in heat-stressed cows during early lactation—which is a high-stress phase—may provide insights into how to deal with the level of damage to dairy cows, through appropriate nutritional and management strategies. We identify that short-term heat stress has a negative effect, to some extent, on feed and water intake, rectal temperature, heart rate, blood hematology and metabolites, milk characteristics, miRNA expression in milk, and metabolomics in blood.

Abstract

This study aims to characterize the influence of short-term heat stress (HS; 4 day) in early lactating Holstein dairy cows, in terms of triggering blood metabolomics and parameters, milk yield and composition, and milk microRNA expression. Eight cows (milk yield = 30 ± 1.5 kg/day, parity = 1.09 ± 0.05) were homogeneously housed in environmentally controlled chambers, assigned into two groups with respect to the temperature humidity index (THI) at two distinct levels: approximately ~71 (low-temperature, low-humidity; LTLH) and ~86 (high-temperature, high-humidity; HTHH). Average feed intake (FI) dropped about 10 kg in the HTHH group, compared with the LTLH group (p = 0.001), whereas water intake was only numerically higher (p = 0.183) in the HTHH group than in the LTLH group. Physiological parameters, including rectal temperature (p = 0.001) and heart rate (p = 0.038), were significantly higher in the HTHH group than in the LTLH group. Plasma cortisol and haptoglobin were higher (p < 0.05) in the HTHH group, compared to the LTLH group. Milk yield, milk fat yield, 3.5% fat-corrected milk (FCM), and energy-corrected milk (ECM) were lower (p < 0.05) in the HTHH group than in the LTLH group. Higher relative expression of milk miRNA-216 was observed in the HTHH group (p < 0.05). Valine, isoleucine, methionine, phenylalanine, tyrosine, tryptophan, lactic acid, 3-phenylpropionic acid, 1,5-anhydro-D-sorbitol, myo-inositol, and urea were decreased (p < 0.05). These results suggest that early lactating cows are more vulnerable to short-term (4 day) high THI levels—that is, HTHH conditions—compared with LTLH, considering the enormous negative effects observed in measured blood metabolomics and parameters, milk yield and compositions, and milk miRNA-216 expression.

Genome-wide identification and functional prediction of circular RNAs in response to heat stress in Chinese Holstein cows

Heat stress (HS) leads to substantial economic loss of dairy industry each year. The negative effect of HS in dairy cows is becoming one of the more urgent issue due to accelerating side-effects of global warming. Various genes are involved in HS response but the information about the role of noncoding RNAs, especially circular RNAs (circRNAs) is largely unknown. In our study, we aimed to investigate the different expression profile of circRNAs between HS and Non-heat-stressed condition (NC) of Chinese Holstein cow's mammary gland. CircRNAs were identified using RNA sequencing and bioinformatics analysis. In total, 37405 circRNAs were detected and 95 were differentially expressed (DE), including 15 downregulated and 80 upregulated circRNAs in HS group compared to NC. Eight circRNAs were randomly selected to verify the RNA sequencing result. Further, Sanger sequencing validated the backsplicing site of the eight circRNAs. Moreover, results obtained from the Quantitative real time PCR (qRT-PCR) showed consistent expression trend with that of RNA sequencing. GO annotation and KEGG analysis suggested that these DE circRNAs probably involved in the energy metabolic regulation. Furthermore, we constructed ceRNA network and the result indicated that these DE circRNAs could regulate lactation through IGF1 and PRL signaling pathway.

Effects of seasonal ambient heat stress on expression of microRNAs in the mammary gland of Holstein cows

This study was conducted to assess the link of miRNA expressions in cow's mammary gland undergoing heat stress. Twelve Holstein cows were allocated either to undergo heat stress (HS) or remain in a thermoneutral environment (non-heat stress, NS), respectively. The experiment with HS cows was carried out in August, and the experiment with NS cows was done in November. After a month, three cows from each group were slaughtered, and mammary gland samples were obtained, and then miRNA were extracted from the samples for later sequencing. From the miRNA-seq, we obtained a total of 124 differentially expressed miRNAs in HS and NS cows' mammary gland. The differentially expressed miRNA could be predicted to influence multiple target genes. The target interleukin-1 (IL-1), which play a role in regulating the function of mammary gland in dairy cows, could be affected by bta-let-7c, bta-let-7e, bta-miR-181d, bta-miR-452, and bta-miR-31. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that mitogen-activated protein kinase (MAPK) pathway plays an important role in the mammary glands of dairy cows and bta-miR-25 and bta-miR-382 may influence MAPK pathway through c-Jun N-terminal kinase (JNK) gene to affect the function of mammary gland in HS cows. In conclusion, this study characterized expression profile of miRNAs in the Holstein cows' mammary gland under summer heat stress or not. We observed miRNA expression during heat stress, which was significantly different from non-heat stress states. A comprehensive analysis of the miRNA's expression will be helpful to further study the link of miRNAs with mechanisms regulating heat stress in the cow mammary gland.

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.

The Mediation of miR-34a/miR-449c for Immune Cytokines in Acute Cold/Heat-Stressed Broiler Chicken

Simple Summary

In the intensive and scale poultry industry, the level of heat stress (HS) directly affects the growth, development, and production performance of poultry. To alleviate the adverse effects of stress in broilers, microRNA (miRNA) was regarded as a potential regulator of immune cytokines. In this study, through the sequencing analysis of spleens after cold/heat stress, we found that 33 and 37 miRNA were differentially expressed in the heat stress group compared with the normal (NS) group and cold stress (CS) group, respectively. The differential miRNA were mainly involved in biological processes such as the cytokine–cytokine receptor interaction. To further understand the miRNA-mediated effect of heat stress on the immune level of chickens, we selected miR-34a and miR-449c as the research objects, predicted and verified that interleukin 2 (IL-2) and interleukin 12α (IL-12α) were the target genes of miR-34a and miR-449c. Coupled with the analysis of the expression of other cytokines, we found that miRNA could change the expression of immune cytokines directly or indirectly. This discovery provides a new insight into the mediation of miRNA for immune cytokines in acute cold/heat stressed broiler chicken.

Abstract

An increasing amount of evidence has revealed that microRNAs (miRNAs) participated in immune regulation and reaction to acute cold and heat stresses. As a new type of post-transcriptional regulatory factor, miRNA has received widespread attention; However, the specific mechanism used for this regulation still needs to be determined. In this study, thirty broilers at the same growth period were divided into three groups and treated with different temperature and humidity of CS (10–15 °C and 90% Relative Humidity (RH)), HS (39 °C and 90% RH), and NS (26 °C and 50–60% RH) respectively. After 6 h, splenic tissues were collected from all study groups. miRNA sequencing was performed to identify the differentially expressed miRNAs (DEMs) between HS, CS, and NS. We found 33, 37, and 7 DEMs in the HS-NS, HS-CS, CS-NS group. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that DEMs were significantly enriched in cytokine–cytokine receptor interaction and functioned as the cellular responders to stress. We chose two miRNA, miR-34a and miR-449c, from the same family and differential expressed in HS-CS and HS-NS group, as the research objects to predict and verify the target genes. The dual-luciferase reporter assay and quantitative real-time PCR (qRT-PCR) confirmed that two cytokines, IL-2 and IL-12α, were the direct target genes of miR-34a and miR-449c. To further understand the mediation mechanism of miRNAs in acute cold/heat-stressed broiler chicken, a splenic cytokines profile was constructed. The results showed that IL-1β was strongly related to acute heat stress in broiler chicken, and from this we predicted that the increased expression of IL-1β might promote the expression of miR-34a, inducing the upregulation of interferon-γ (INF-γ) and IL-17. Our finds have laid a theoretical foundation for the breeding of poultry resistance and alleviation of the adverse effects of stress.

Significance of molecular chaperones and micro RNAs in acquisition of thermo-tolerance in dairy cattle

Ambient temperature is considered as the major abiotic factor which regulates body physiological mechanisms of all living creatures across the globe. Variation in ambient temperature which emulates thermoneutral zone culminates in heat stress. Heat stress has been emerged as major ultimatum to livestock's growth, development, production and reproduction across the world. Livestock's responds to the heat stress via different mechanisms such as behavioral, physiological, biochemical, endocrine and molecular mechanisms. Amongst the aforementioned mechanisms, molecular mechanism plays crucial role to achieve thermo-tolerance via expression of highly conserved family of proteins known as heat shock proteins (HSPs) across livestock species. HSPs serve as molecular chaperones to ameliorate the menace of heat stress in domestic species. In addition, microRNAs are small non-coding RNA which down regulates post-transcriptional gene expression by targeting various HSPs to regulate the thermoregulatory responses in livestock species. Despite of thermal adaptation mechanisms, heat stress breaches animal body homeostasis thereby depresses their production and productivity. Therefore, veterinary researches have been targeting to explore different repertoire of HSPs and microRNAs expression to counteract the rigors of heat stress thereby confer thermo-tolerance in livestock species. The present review highlights the significance of molecular chaperones and microRNAs in the acquisition of thermo-tolerance in dairy cattle.

Identification of Ovine Serum miRNAs Following Bacterial Lipopolysaccharide Challenge

Host-pathogen interactions are complex and influenced by host genetic and epigenetic modifications. Recently, the significance of microRNAs (miRNAs) in pathogenic infection and the regulation of immune response has been highlighted. However, information on miRNAs' role in the course of inflammation is still very limited in small ruminants. The present study was intended to identify changes in the expression of circulatory miRNAs post-lipopolysaccharide (LPS)-challenge. In this study, young ewes (n = 18) were challenged with Escherichia coli LPS (400 ng/kg i.v.) and blood samples were collected for serum miRNA isolation at two-time points; prior to challenge (T0), and 4 h (T4) post-challenge, reflecting the peak cortisol response. A total of 91 miRNAs were profiled, including 84 miRNAs on a commercial ovine miRNA-PCR array, and seven individual miRNAs. Forty five miRNAs were differentially expressed (DE) with 35 being up-regulated (Fold regulation, FR > 2) and 10 being down-regulated (FR < 1, p < 0.05) at T4. Among the up-regulated miRNAs, 14 were significantly (p < 0.05) induced, including oar-miRs: 369-3p, 495-3p, 376a-3p, 543-3p, 668-3p, 329a-3p, 655-3p, 411a-5p, and 154a-3p, which were located on ovine chromosome 18 forming four miRNA clusters within 10 kb. The elevated miRNAs belonged to different functional classes, playing roles in activating the hypothalamic-pituitary-adrenal axis; increasing cell survival and differentiation; and inducing inflammatory responses and targeted PI3K-Akt and MAPK signaling and chemokine signaling pathways. In summary, these results reveal the dynamic nature of ovine serum miRNAs during LPS-induced stress and highlight the potential role of identified miRNA-clusters on chromosome 18 to understand the regulation of the acute-phase response. Some of these identified circulating miRNAs may also serve as stress biomarkers for livestock in the future.

Effect of heat stress on udder health of dairy cows

This Research Reflection short review presents an overview of the effects of heat stress on dairy cattle udder health and discusses existing heat stress mitigation strategies for a better understanding and identification of appropriate abatement plans for future stress management. Due to high ambient temperatures with high relative humidity in summer, dairy cows respond by changes of physical, biochemical and biological pathways to neutralize heat stress resulting in decreased production performance and poorer immunity resulting in an increased incidence of intramammary infections (IMI) and a higher somatic cell count (SCC). In vitro studies on bovine polymorphonuclear cells (PMN) suggested that heat stress reduces the phagocytosis capacity and oxidative burst of PMN and alters the expression of apoptotic genes and miRNA which, together with having a negative effect on the immune system, may explain the increased susceptibility to IMI. Although there are limited data regarding the incidence rate of clinical mastitis in many countries or regions, knowledge of SCC at the cow or bulk tank level helps encourage farmers to improve herd health and to develop strategies for infection prevention and cure. Therefore, more research into bulk tank SCC and clinical mastitis rates is needed to explain the effect of heat stress on dairy cow udder health and functions that could be influenced by abatement plans.

Analysis of circulating-microRNA expression in lactating Holstein cows under summer heat stress

Korean peninsula weather is rapidly becoming subtropical due to global warming. In summer 2018, South Korea experienced the highest temperatures since the meteorological observations recorded in 1907. Heat stress has a negative effect on Holstein cows, the most popular breed of dairy cattle in South Korea, which is susceptible to heat. To examine physiological changes in dairy cows under heat stress conditions, we analyzed the profiles circulating microRNAs isolated from whole blood samples collected under heat stress and non-heat stress conditions using small RNA sequencing. We compared the expression profiles in lactating cows under heat stress and non-heat stress conditions to understand the regulation of biological processes in heat-stressed cows. Moreover, we measured several heat stress indicators, such as rectal temperature, milk yield, and average daily gain. All these assessments showed that pregnant cows were more susceptible to heat stress than non-pregnant cows. In addition, we found the differential expression of 11 miRNAs (bta-miR-19a, bta-miR-19b, bta-miR-30a-5p, and several from the bta-miR-2284 family) in both pregnant and non-pregnant cows under heat stress conditions. In target gene prediction and gene set enrichment analysis, these miRNAs were found to be associated with the cytoskeleton, cell junction, vasculogenesis, cell proliferation, ATP synthesis, oxidative stress, and immune responses involved in heat response. These miRNAs can be used as potential biomarkers for heat stress.

Effects of chronic heat stress on mRNA and miRNA expressions in dairy cows

Heat stress has a negative impact on dairy cow productivity. In order to reveal the mechanisms of heat-stress response, the mRNA and miRNA expression profiles of five cows under chronic heat-stress and thermoneutral conditions were assayed in blood by high-throughput sequencing technology. A total of 540 mRNAs and 9 miRNAs were expressed differently under heat-stress and thermoneutral conditions. Functional analyses revealed that MAPK signaling pathway, cellular senescence, circadian entrainment, aldosterone synthesis and secretion, and pathways in cancer were enriched for differently expressed mRNAs; meanwhile cGMP-PKG signaling pathway, thermogenesis, and protein digestion and absorption were enriched for differently expressed miRNAs. In addition, GADD45G, TGFB2, and GNG11 may play roles in the heat stress, and bta-miR-423-5p might be one of the regulators of heat-stress response in cows as potential mediators of chronic heat-stress response. In conclusion, the present study described the mRNA and miRNA expression patterns in blood extracted from cows during the transition from heat-stress to thermoneutral conditions. The results provide new data that could help in identifying mechanisms that mediate cows' response to chronic heat stress.

Short communication: Characterizing ovine serum stress biomarkers during endotoxemia

Breeding stress-resilient livestock is a potential strategy to help mitigate the negative effect of environmental and pathogenic stressors. The hypothalamic-pituitary-adrenal axis and immune system are activated during stress events and release mediators into the circulation that help restore physiological homeostasis. The purpose of this study was to assess a comprehensive set of circulatory mediators released in response to an acute immune stress challenge to identify candidate biomarkers that can be used for the selection of stress-resilient animals. Fifteen female lambs were stress challenged with an intravenous bolus of lipopolysaccharide (LPS; 400 ng/kg), and blood was collected from the jugular vein at 0, 2, 4, and 6 h after LPS challenge to identify and monitor candidate stress biomarkers; temperature was also recorded over time. Biomarker responses were evaluated with a repeated-measures model to compare time points with baseline values. As expected, all sheep had a monophasic febrile response to LPS challenge, and cortisol increased and returned to baseline by 6 h. The cytokines tumor necrosis factor-α, IL-6, IFN-γ (proinflammatory), and IL-10 (anti-inflammatory) increased, but only tumor necrosis factor-α returned to baseline during the monitoring period. The cytokines IL-1α, IL-1β, IL-17α (proinflammatory), and IL-4 (anti-inflammatory) did not respond to LPS challenge. All chemokines (CCL2, CCL3, CCL4, CXCL10, and IL-8) responded to LPS challenge; however, only CCL2, CCL3, CCL4, and CXCL10 increased over time, and only CCL3, CCL4, and CXCL10 returned to baseline during the monitoring period. MicroRNA (miR-145, miR-233, and miR-1246) also increased and remained elevated during the study. In summary, the LPS challenge induced a strong stress response in Rideau-Dorset sheep that could be monitored with a distinct profile of circulatory biomarkers.