Search

Reproductive traits of dairy cattle are bound to the actual efficiency of farm operation, which therefore show great economic importance. Among them, some traits were deemed to be simultaneously affected by service sire and mating cow. Service sires are proved to play an important role in reproduction process of cows. However, limited study explored the genetic effect of service sire (GESS), let alone the genomic prediction of this effect. In the present study, 2244 genotyped bulls together with phenotypic records were used to predict the GESS on conception rate, 56-day non-return rate, calving ease, stillbirth and gestation length. The feasibilities of multi-step genomic best linear unbiased predictor (msGBLUP) and single-step genomic best linear unbiased predictor (ssGBLUP) were investigated under different scenarios, that is, different marker densities and validation population. The predictive accuracies and unbiasedness for GESS ranged from 0.159 to 0.647 and from 0.202 to 2.018, respectively, when validated on young bulls, while the accuracies and unbiasedness ranged from 0.409 to 0.802 and 0.333 to 1.146 when validated on random split data sets. It is feasible to predict GESS on reproductive traits by using a linear mixed model and genomic data, and high-density marker panel had limited contribution to the prediction. This research investigated the potential factors that influence the genomic prediction of GESS on reproductive traits and indicated the possibility of genomic selection on GESS, both in ideal and practical circumstances.

The main objectives of this study were to estimate genetic parameters for milk urea nitrogen (MUN) in Holstein cattle and to conduct a single-step (ss)GWAS to identify candidate genes associated with MUN. Phenotypic measurements from 24,435 Holstein cows were collected from March 2013 to July 2019 in 9 dairy farms located in the Beijing area, China. A total of 2,029 cows were genotyped using the Illumina 150K Bovine Bead Chip, containing 121,188 SNP. A single-trait repeatability model was used to evaluate the genetic background of MUN. We found that MUN is a trait with low heritability (0.06 ± 0.004) and repeatability (0.12). Considering similar milk production levels, a lower MUN concentration indicates higher nitrogen digestibility. The genetic correlations between MUN and milk yield, net energy concentration, fat percentage, protein percentage, and lactose percentage were positive and ranged from 0.02 to 0.26. The genetic correlation between MUN and somatic cell score (SCS) was negative (-0.18), indicating that animals with higher MUN levels tend to have lower SCS. Both ssGWAS and pathway enrichment analyses were used to explore the genetic mechanisms underlying MUN. A total of 18 SNP (located on BTA11, BTA12, BTA14, BTA17, and BTA18) were found to be significantly associated with MUN. The genes CFAP77, CAMSAP1, CACNA1B, ADGRB1, FARP1, and INTU are considered to be candidate genes for MUN. These candidate genes are associated with important biological processes such as protein and lipid metabolism and binding to specific proteins. This set of candidate genes, metabolic pathways, and their functions provide a better understanding of the genomic architecture and physiological mechanisms underlying MUN in Holstein cattle.

Enhancing the immune response through breeding is regarded as an effective strategy for improving animal health, as dairy cattle identified as high immune responders are reported to have a decreased prevalence of economically significant diseases. The identification of differentially expressed genes (DEGs) associated with immune responses might be an effective tool for breeding healthy dairy cattle. In this study, antibody-mediated immune responses (AMIRs) were induced by the immunization of hen egg white lysozyme (HEWL) in six Chinese Holstein dairy bulls divided into high- and low-AMIR groups based on their HEWL antibody level. Then, RNA-seq was applied to explore the transcriptome of peripheral whole blood between the two comparison groups. As a result, several major upregulated and downregulated genes were identified and attributed to the regulation of locomotion, tissue development, immune response, and detoxification. In addition, the result of the KEGG pathway analysis revealed that most DEGs were enriched in pathways related to disease, inflammation, and immune response, including antigen processing and presentation, Staphylococcus aureus infection, intestinal immune network for IgA production, cytokine–cytokine receptor interaction, and complement and coagulation cascades. Moreover, six genes (BOLA-DQA5, C5, CXCL2, HBA, LTF, and COL1A1) were validated using RT-qPCR, which may provide information for genomic selection in breeding programs. These results broaden the knowledge of the immune response mechanism in dairy bulls, which has strong implications for breeding cattle with an enhanced AMIR.

Heat stress has been a big challenge for animal survival and health due to global warming. However, the molecular processes driving heat stress response were unclear. In this study, we exposed the control group rats (n = 5) at 22 °C and the other three heat stress groups (five rats in each group) at 42 °C lasting 30, 60, and 120 min, separately. We performed RNA sequencing in the adrenal glands and liver and detected the levels of hormones related to heat stress in the adrenal gland, liver, and blood tissues. Weighted gene co-expression network analysis (WGCNA) was also performed. Results showed that rectal temperature and adrenal corticosterone levels were significantly negatively related to genes in the black module, which was significantly enriched in thermogenesis and RNA metabolism. The genes in the green-yellow module were strongly positively associated with rectal temperature and dopamine, norepinephrine, epinephrine, and corticosterone levels in the adrenal glands and were enriched in transcriptional regulatory activities under stress. Finally, 17 and 13 key genes in the black and green-yellow modules were identified, respectively, and shared common patterns of changes. Methyltransferase 3 (Mettl3), poly(ADP-ribose) polymerase 2 (Parp2), and zinc finger protein 36-like 1 (Zfp36l1) occupied pivotal positions in the protein–protein interaction network and were involved in a number of heat stress-related processes. Therefore, Parp2, Mettl3, and Zfp36l1 could be considered candidate genes for heat stress regulation. Our findings shed new light on the molecular processes underpinning heat stress.

Udder performance and health are important traits in dairy cattle worldwide. The present study aimed to explore the association of single nucleotide polymorphisms (SNPs) in the STAT5B gene in dairy cattle with milk performance and mastitis related traits. The study included 201 cows of three pure breeds (i.e. Holstein Friesian, Jersey and Achai) and two crossbred cattle at four established dairy farms in Khyber Pakhtunkhwa, Pakistan. The milk samples were analysed for somatic cell count (SCC) and milk composition (i.e. fat , protein and lactose percentages). The generalized linear model was deployed for association analysis using SAS. The pool DNA sequencing showed four (three synonymous and a 3/ UTR) SNPs in STAT5B. These SNPs were further validated in all DNA samples using SNaPshot assay. The breed-wise analysis showed that most of the SNPs were consistent with Hardy-Weinberg equilibrium (P>0.05). The association analysis revealed a significantly higher protein percentage in TT genotype and lower SCC in CC genotypes of SNP 1 (exon 2, C>T), whereas in SNP 2 (exon 16, T>C) the TT genotypes revealed significantly lower SCC and SCS compared with other genotypes (PT) showed significantly lower SCC and frequency of clinical mastitis in the heterozygous (CT) genotype compared to the homozygous genotypes. The in silico predictions revealed changes in the RNA secondary structure for SNP 2 and SNP 3. The study suggests that STAT5B should be considered as a candidate gene, and the variants identified as useful genetic markers for improved milk composition and udder health.

Supernumerary teats (SNT) are a common epidermal abnormality of udders in mammals. The SNT negatively affect machine milking ability, udder health, and animal welfare and sometimes act as reservoirs for undesirable bacteria, resulting in economic losses on calves and lactating cows due to the cost of SNT removal surgery, early culling, and low milk yield. This study aimed to analyze the incidence and genetic parameter of SNT and detect SNT-related genes in Chinese Holstein cattle. In this study, the incidence of SNT was recorded in 4,670 Chinese Holstein cattle (born between 2008 and 2017) from 2 farms, including 734 genotyped cows with 114,485 SNPs. The SNT had a total frequency of 9.8% and estimated heritability of 0.22 (SE = 0.07), which were obtained using a threshold model in the studied Chinese Holstein population. Furthermore, we calculated approximate genetic corre-lations between SNT and the following indicator traits: 12 milk production, 28 body conformation, 5 fertility and reproduction, 5 health, and 9 longevity. Generally, the estimated correlations, such as 305-d milk yield for third parity (-0.55; SE = 0.02) and age at first calving in heifer (0.19; SE = 0.03), were low to moderate. A single-step GWAS was implemented, and 10 genes asso-ciated with SNT located in BTA4 were identified. The region (112.70-112.90 Mb) on BTA4 showed the high -est genetic variance for SNT. The quantitative trait loci on BTA4 was mapped into the RARRES2 gene, which was previously shown to affect adipogenesis and hormone secretion. The WIF1 gene, which was located in BTA5, was also considered as a candidate gene for SNT. Overall, these findings provide useful information for breeders who are interested in reducing SNT.

Climate changes affects animal physiology. In particular, rising ambient temperatures reduce animal vitality due to heat stress and this can be observed at various levels which included genome, transcriptome, and microbiome. In a previous study, microbiota highly associated with changes in cattle physiology were identified under heat stress conditions which included rectal temperature, drooling score, and respiratory score. In the present study, genes differentially expressed between individuals were selected representing different additive genetic effects towards the heat stress response in cattle in their production condition. Moreover, a correlation network analysis was performed to identify interactions between the transcriptome and microbiome for 71 Chinese Holstein cows sequenced for mRNA from blood samples and for 16S rRNA genes from fecal samples. Bioinformatics analysis was performed comprising: i) clustering and classification of 16S rRNA sequence reads, ii) mapping cows’ transcripts to the reference genome and their expression quantification, and iii) statistical analysis of both data types – including differential gene expression analysis and gene set enrichment analysis. A weighted co-expression network analysis was carried out to assess changes in the association between gene expression and microbiota abundance as well as to find hub genes/microbiota responsible for the regulation of gene expression under heat stress. Results showed 1,851 differentially expressed genes were found that were shared by three heat stress phenotypes. These genes were predominantly associated with the cytokine-cytokine receptor interaction pathway. The interaction analysis revealed three modules of genes and microbiota associated with rectal temperature with which two hubs of those modules were bacterial species, demonstrating the importance of the microbiome in the regulation of gene expression during heat stress. Genes and microbiota from the significant modules can be used as biomarkers of heat stress in cattle.

In our previous GWAS of Chinese and Nordic dairy cattle, genes CACNB2, SLC39A12, and ZEB1 locating on BTA 13 were suggested as candidate genes for reproduction. In this study, validation of these associations was performed in an independent population with records of nine reproductive traits. More importantly, functions of these genes in the reproductive process were verified by employing the expression data of ovarian follicles. The potential variants within the three genes were firstly detected in 68 Chinese Holstein bulls, and then screened in 1,588 Chinese Holstein cows using the KASP (Kompetitive allele-specific PCR) method. There were nine variants with polymorphisms in CACNB2, five in SLC39A12, and four in ZEB1, respectively, of which one SNP was in the upstream regulatory region, two in exon region, four in downstream regulatory region, and 11 SNPs in intronic regions. Amongst the 18 variants, g.33267056T/G in CACNB2 explained the largest phenotypic variance for age at first calving (0.011%), interval from first to last insemination (0.004%), and calving ease (0.002%), while g.32751518G/A in SLC39A12 contributed the most to stillbirth in heifers (0.038%). Two haplotype blocks were constructed for CACNB2 while one each for SLC39A12 and ZEB1, which were significantly associated with five reproductive traits, including age at the first service, age at the first calving, calving ease in heifers and cows, and the interval from calving to the first insemination. We then studied the profile of gene expression in granulosa cells isolated from four developmental stages of ovarian follicles from eight dairy cows. All three genes were differentially expressed between ovarian follicles with different sizes (p < 0.05), indicating their potential roles in the reproductive process of dairy cows. This study successfully demonstrated the associations of three BTA 13 genes CACNB2, SLC39A12, and ZEB1 with reproduction and further examined their expression levels in ovarian follicles directly. These findings can be beneficial for the ongoing genomic selection program for reproductive traits which have long been considered as traits that are difficult to achieve genetic improvement due to the lack of efficient genetic markers.

Heat stress is a major cause of welfare issues and economic losses to the worldwide dairy cattle industry. Genetic selection for heat tolerance has a great potential to positively affect the dairy industry, as the gains are permanent and cumulative over generations. Rectal temperature (RT) is hypothesized to be a good indicator trait of heat tolerance. Therefore, this study investigated the genetic architecture of RT by estimating genetic parameters, performing genome-wide association studies, and biologically validating potential candidate genes identified to be related to RT in Holstein cattle. A total of 33,013 RT records from 7,598 cows were used in this study. In addition, 1,114 cows were genotyped using the Illumina 150K Bovine BeadChip (Illumina, San Diego, CA). Rectal temperature measurements taken in the morning (AMRT) and in the afternoon (PMRT) are moderately heritable traits, with estimates of 0.09 ± 0.02 and 0.04 ± 0.01, respectively. These 2 traits are also highly genetically correlated (r = 0.90 ± 0.08). A total of 10 SNPs (located on BTA3, BTA4, BTA8, BTA13, BTA14, and BTA29) were found to be significantly associated with AMRT and PMRT. Subsequently, gene expression analyses were performed to validate the key functional genes identified (SPAG17, FAM107B, TSNARE1, RALYL, and PHRF1). This was done through in vitro exposure of peripheral blood mononuclear cells (PBMC) to different temperatures (37°C, 39°C, and 42°C). The relative mRNA expression of 2 genes, FAM107B and PHRF1, significantly changed between the control and heat stressed PBMC. In summary, RT is heritable, and enough genetic variability exists to enable genetic improvement of heat tolerance in Holstein cattle. Important genomic regions were identified and biologically validated; FAM107B and PHRF1 are the main candidate genes identified to influence heat stress response in dairy cattle.

Previous studies reported the physical, transcriptomics and metabolomics changes in in-vitro acute heat stressed bovine granulosa cells. Granulosa cells exhibited transient proliferation senescence, oxidative stress, increased rate of apoptosis, and decline in steroidogenic activity. This study performs joint integration and network analysis of metabolomics and transcriptomics data to further narrow down and elucidate the role of differentially expressed genes, important metab-olites and relevant cellular and metabolic pathways in acute heat-stressed granulosa cells. Among significant (Raw P-value

Heat stress affects the granulosa cells (GCs) and ovarian follicular microenvironment, causing poor oocyte developmental competence and fertility. This study aimed to investigate the physical responses and global transcriptomic changes in bovine GCs to acute heat stress (43 ℃ for 2 h) in-vitro and gave essential insights into the general interaction at cell–stress nexus. Heat-stressed GCs exhibited transient proliferation senescence, resumed proliferation at 48 h post-stress. While post-stress immediate culture-media change had a relatively positive effect on proliferation resumption. Increased accumulation of reactive oxygen species and apoptosis was observed in heat stress group. In spite of the upregulation of pro-apoptotic and caspase executioner genes, antioxidants and anti-apoptotic genes were also upregulated in heat-stressed GCs. Progesterone and Estrogen hormones along with steroidogenic genes expression, declined significantly, in spite of the upregulation of genes involved in cholesterol synthesis. Out of 12385 differentially expressed genes (DEGs), 330 significant DEGs (75 upregulated, 225 downregulated) were subjected to KEGG functional pathway annotation, gene ontology enrichment, and STRING network analyses. Based on the manual query of DEGs, pathway and enrichment analyses, a vast interplay observed among all major signaling pathways strongly evidence the repression of cellular transcriptional and proliferation activity, averting the effects of heat stress through remodeling of cellular structural proteins and energetic-homeostasis. This study presents detailed responses of acute heat-stressed GCs at physical, transcriptional, and pathway levels and presents interesting insights into future studies regarding GCs adaptation and their interaction with oocyte and reproductive system at ovarian level.

In our previous GWAS of Chinese and Nordic dairy cattle, genes

Heat stress affects granulosa cells and the ovarian follicular microenvironment, ultimately resulting in poor oocyte developmental competence. This study aims to investigate the metabo-lomics response of bovine granulosa cells (bGCs) to in vitro acute heat stress of 43 °C. Heat stress triggers oxidative stress-mediated apoptosis in cultured bGCs. Heat-stressed bGCs exhibited a time-dependent recovery of proliferation potential by 48 h. A total of 119 metabolites were identified through LC–MS/MS-based metabolomics of the spent culture media, out of which, 37 metabolites were determined as differentially involved in metabolic pathways related to bioenergetics support mechanisms and the physical adaptations of bGCs. Multiple analyses of metabolome data identified choline, citric acid, 3-hydroxy-3-methylglutaric acid, glutamine, and glycocyamine as being upregulated, while galactosamine, AICAR, ciliatine, 16-hydroxyhexadecanoic acid, lysine, succinic acid, uridine, xanthine, and uraconic acid were the important downregulated metabolites in acute heat stress. These differential metabolites were implicated in various important metabolic pathways directed towards bioenergetics support mechanisms including glycerophospholipid metabolism, the citrate cycle (TCA cycle), glyoxylate and dicarboxylate metabolism, and serine, threonine, and tyrosine metabolism. Our study presents important metabolites and metabolic pathways involved in the adaptation of bGCs to acute heat stress in vitro.

Simple Summary Milk production is an important trait in the breeding and genetic improvement of Xinjiang Brown cattle. To obtain the best strategy for improving the reliability of the breeding value estimation for each trait, we used single-trait and multitrait models based on the A-array pedigree-based best linear unbiased prediction (PBLUP) and H-array single-step genomic best linear unbiased prediction (ssGBLUP) to perform the genetic evaluation of different strategies using the restricted maximum likelihood (REML) and Bayesian methods. Upon comparison, the ssGBLUP calculation results of the multitrait models obtained using the REML and Bayesian methods were better than those of other strategies. Considering the calculation time, the multitrait model REML method is recommended for ssGBLUP calculation to accurately predict the breeding value of young animals; thus, this strategy should be used for the early breeding selection of Xinjiang Brown cattle. Abstract One-step genomic selection is a method for improving the reliability of the breeding value estimation. This study aimed to compare the reliability of pedigree-based best linear unbiased prediction (PBLUP) and single-step genomic best linear unbiased prediction (ssGBLUP), single-trait and multitrait models, and the restricted maximum likelihood (REML) and Bayesian methods. Data were collected from the production performance records of 2207 Xinjiang Brown cattle in Xinjiang from 1983 to 2018. A cross test was designed to calculate the genetic parameters and reliability of the breeding value of 305 daily milk yield (305 dMY), milk fat yield (MFY), milk protein yield (MPY), and somatic cell score (SCS) of Xinjiang Brown cattle. The heritability of 305 dMY, MFY, MPY, and SCS estimated using the REML and Bayesian multitrait models was approximately 0.39 (0.02), 0.40 (0.03), 0.49 (0.02), and 0.07 (0.02), respectively. The heritability and estimated breeding value (EBV) and the reliability of milk production traits of these cattle calculated based on PBLUP and ssGBLUP using the multitrait model REML and Bayesian methods were higher than those of the single-trait model REML method; the ssGBLUP method was significantly better than the PBLUP method. The reliability of the estimated breeding value can be improved from 0.9% to 3.6%, and the reliability of the genomic estimated breeding value (GEBV) for the genotyped population can reach 83%. Therefore, the genetic evaluation of the multitrait model is better than that of the single-trait model. Thus, genomic selection can be applied to small population varieties such as Xinjiang Brown cattle, in improving the reliability of the genomic estimated breeding value.

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.

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

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. In this study, we identified fecal microbiota associated with heat stress that was quantified by three metrics – rectal temperature, drooling, and respiratory scores and represented by their Estimated Breeding Values. For this purpose, the 16S rRNA sequencing technique was used. We analyzed the microbiota from 136 fecal samples of Chinese Holstein cows through a 16S rRNA gene sequencing approach. Sequence data were processed using a pipeline involving QIIME2 software together with SILVA database. 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 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. Moreover, identified microbiota significantly associated with heat stress measures were mostly mesophilic, so their association seems to be due to heat stress-induced secondary, metabolic changes, and not directly by temperature. Moreover, high fold changes of many genera suggest that they may be used as biomarkers for monitoring the level of heat stress in cattle.

BackgroundCharacterization of the molecular mechanisms underlying hair follicle development is of paramount importance in the genetic improvement of wool-related traits in sheep and skin-related traits in humans. The Merino is the most important breed of fine-wooled sheep in the world. In this study, we systematically investigated the complexity of sheep hair follicle development by integrating transcriptome and methylome datasets from Merino sheep skin.ResultsWe analysed 72 sequence datasets, including DNA methylome and the whole transcriptome of four gene types, i.e. protein-coding genes (PCGs), lncRNAs, circRNAs, and miRNAs, across four embryonic days (E65, E85, E105, and E135) and two postnatal days (P7 and P30) from the skin tissue of 18 Merino sheep. We revealed distinct expression profiles of these four gene types across six hair follicle developmental stages, and demonstrated their complex interactions with DNA methylation. PCGs with stage-specific expression or regulated by stage-specific lncRNAs, circRNAs, and miRNAs were significantly enriched in epithelial differentiation and hair follicle morphogenesis. Regulatory network and gene co-expression analyses identified key transcripts controlling hair follicle development. We further predicted transcriptional factors (e.g. KLF4, LEF1, HOXC13, RBPJ, VDR, RARA, and STAT3) with stage-specific involvement in hair follicle morphogenesis. Through integrating these stage-specific genomic features with results from genome-wide association studies (GWAS) of five wool-related traits in 7135 Merino sheep, we detected developmental stages and genes that were relevant with wool-related traits in sheep. For instance, genes that were specifically upregulated at E105 were significantly associated with most of wool-related traits. A phenome-wide association study (PheWAS) demonstrated that candidate genes of wool-related traits (e.g.SPHK1,GHR,PPP1R27,CSRP2,EEF1A2, andPTPN1) in sheep were also significantly associated with dermatological, metabolic, and immune traits in humans.ConclusionsOur study provides novel insights into the molecular basis of hair follicle morphogenesis and will serve as a foundation to improve breeding for wool traits in sheep. It also indicates the importance of studying gene expression in the normal development of organs in understanding the genetic architecture of economically important traits in livestock. The datasets generated here are useful resources for functionally annotating the sheep genome, and for elucidating early skin development in mammals, including humans.

In this study, we estimated the genetic parameters for 6 composite traits and 27 body conformation traits of 1016 dual-purpose Simmental cattle reared in northwestern China from 2010 to 2019 using a linear animal mixed model. To integrate these traits, a variety of methods were used as follows: (1) genetic parameters estimates for composite and individual body conformation traits based on the pedigree relationship matrix (A) and combined genomic-pedigree relationship matrix (H); (2) factor analysis to explore the relationships among body conformation traits; and (3) genetic parameters of factor scores estimated using A and H, and the correlations of EBVs of the factor scores and EBVs of the composite traits. Heritability estimates of the composite traits using A and H were low to medium (0.07–0.47). The 24 common latent factors explained 96.13% of the total variance. Among factors with eigenvalues ≥ 1, F1 was mainly related to body frame, muscularity, and rump; F2 was related to feet and legs; F3, F4, F5, and F6 were related to teat placement, teat size, udder size, and udder conformation; and F7 was related to body frame. Single-trait analysis of factor scores yielded heritability estimates that were low to moderate (0.008–0.43 based on A and 0.04–0.43 based on H). Spearman and Pearson correlations, derived from the best linear unbiased prediction analysis of composite traits and factor scores, showed a similar pattern. Thus, incorporating factor analysis into the morphological evaluation to simplify the assessment of body conformation traits may improve the genetics of dual-purpose Simmental cattle.

Molecular responses to heat stress are multifaceted and under a complex cellular post-transcriptional control. This study explores the epigenetic and transcriptional alterations induced by heat stress (42 °C for 120 min) in the liver of rats, by integrating ATAC-seq, RNA-Seq, and WGBS information. Out of 2586 differential ATAC-seq peaks induced by heat stress, 36 up-regulated and 22 down-regulated transcript factors (TFs) are predicted, such as Cebpα, Foxa2, Foxo4, Nfya and Sp3. Furthermore, 150,189 differentially methylated regions represent 2571 differentially expressed genes (DEGs). By integrating all data, 45 DEGs are concluded as potential heat stress response markers in rats. To comprehensively annotate and narrow down predicted markers, they are integrated with GWAS results of heat stress parameters in cows, and PheWAS data in humans. Besides better understanding of heat stress responses in mammals, INSR, MAPK8, RHPN2 and BTBD7 are proposed as candidate markers for heat stress in mammals.

Background Genetic improvement of wool and growth traits is a major goal in the sheep industry, but their underlying genetic architecture remains elusive. To improve our understanding of these mechanisms, we conducted a weighted single-step genome-wide association study (WssGWAS) and then integrated the results with large-scale transcriptome data for five wool traits and one growth trait in Merino sheep: mean fibre diameter (MFD), coefficient of variation of the fibre diameter (CVFD), crimp number (CN), mean staple length (MSL), greasy fleece weight (GFW), and live weight (LW). Results Our dataset comprised 7135 individuals with phenotype data, among which 1217 had high-density (HD) genotype data (n = 372,534). The genotypes of 707 of these animals were imputed from the Illumina Ovine single nucleotide polymorphism (SNP) 54 BeadChip to the HD Array. The heritability of these traits ranged from 0.05 (CVFD) to 0.36 (MFD), and between-trait genetic correlations ranged from − 0.44 (CN vs. LW) to 0.77 (GFW vs. LW). By integrating the GWAS signals with RNA-seq data from 500 samples (representing 87 tissue types from 16 animals), we detected tissues that were relevant to each of the six traits, e.g. liver, muscle and the gastrointestinal (GI) tract were the most relevant tissues for LW, and leukocytes and macrophages were the most relevant cells for CN. For the six traits, 54 quantitative trait loci (QTL) were identified covering 81 candidate genes on 21 ovine autosomes. Multiple candidate genes showed strong tissue-specific expression, e.g. BNC1 (associated with MFD) and CHRNB1 (LW) were specifically expressed in skin and muscle, respectively. By conducting phenome-wide association studies (PheWAS) in humans, we found that orthologues of several of these candidate genes were significantly (FDR BNC1 was significantly associated with MFD in sheep and with hair colour in humans, and CHRNB1 was significantly associated with LW in sheep and with body mass index in humans. Conclusions Our findings provide novel insights into the biological and genetic mechanisms underlying wool and growth traits, and thus will contribute to the genetic improvement and gene mapping of complex traits in sheep.

Genetic selection for improved climatic resilience is paramount to increase the long-term sustainability of high-producing dairy cattle, especially in face of climate change. Various physiological indicators, such as rectal temperature (RT), respiration rate score (RR), and drooling score (DS), can be used to genetically identify animals with more effective coping mechanisms in response to heat stress events. In this study, we investigated genetic parameters for RT, RR (score from 1–3), and DS (score from 1–3). Furthermore, we assessed the genetic relationship among these indicators and other economically important traits for the dairy cattle industry. After data editing, 59,265 (RT), 30,290 (RR), and 30,421 (DS) records from 13,592 lactating Holstein cows were used for the analyses. Variance components were estimated based on a multiple-trait repeatability animal model. The heritability ± standard error estimate for RT, RR, and DS was 0.06 ± 0.01, 0.04 ± 0.01, and 0.02 ± 0.01, respectively, whereas their repeatability was 0.19, 0.14, and 0.14, respectively. Moderate genetic correlations of RR with RT and DS (0.26 ± 0.11 and 0.25 ± 0.16) and nonsignificant correlation between RT and DS (−0.11 ± 0.14) were observed. Furthermore, the approximate genetic correlations between RT, RR, and DS with 12 production, 29 conformation, 5 fertility and reproduction, 5 health, and 9 longevity-indicator traits were assessed. In general, the approximate genetic correlations calculated were low to moderate. In summary, 3 physiological indicators of heat stress response were measured in a large number of animals and shown to be lowly heritable. There is a value in developing a selection index including all the 3 indicators to improve heat tolerance in dairy cattle. All the unfavorable genetic relationships observed between heat tolerance and other economically important traits can be accounted for in a selection index to enable improved climatic resilience while also maintaining or increasing productivity in Holstein cattle.

Due to the complexity of longevity trait in dairy cattle, two groups of trait definitions are widely used to measure longevity, either covering the full lifespan or representing only a part of it to achieve an early selection. Usually, only one group of longevity definition is used in breeding program for one population, and genetic studies on the comparisons of two groups of trait definitions are scarce. Based on the data of eight traits well representing the both groups of trait definitions, the current study investigated genetic parameters and genetic architectures of longevity in Holsteins. Heritabilities and correlations of eight longevity traits were estimated using single-trait and multi-trait animal models, with the data from 103,479 cows. Among the cows with phenotypes, 2,630 cows were genotyped with the 150K-SNP panel. A single-trait fixed and random Circuitous Probability Unification model was performed to detect candidate genes for eight longevity traits. Generally, all eight longevity traits had low heritabilities, ranging from 0.038 for total productive life and herd life to 0.090 for days from the first calving to the end of first lactation or culling. High genetic correlations were observed among the traits within the same definition group: from 0.946 to 0.997 for three traits reflecting full lifespan and from 0.666 to 0.997 for five traits reflecting partial productive life. Genetic correlations between two groups of traits ranged from 0.648 to 0.963, and increased gradually with the extension of lactations number regarding the partial productive life traits. A total of 55 SNPs located on 25 chromosomes were found genome-wide significantly associated with longevity, in which 12 SNPs were associated with more than one trait, even across traits of different definition groups. This is the first study to investigate the genetic architecture of longevity representing both full and the partial lifespan simultaneously, which will assist the selection of an appropriate trait definition for genetic improvement of longevity. Because of high genetic correlations with the full lifespan traits and higher heritability, the partial productive life trait measured as the days from the first calving to the end of the third lactation or culling could be a good alternative for early selection on longevity. The candidate genes identified by this study, such as RPRM, GRIA3, GTF2H5, CA5A, CACNA2D1, FGF10, and DNAJA3, could be used to pinpoint causative mutations for longevity and further benefit the genomic improvement of longevity in dairy cattle.

Chronic stress is a risk factor for a variety of physiological disorders because of its increased activation of the hypothalamic-pituitary-adrenal (HPA) axis; however, it is difficult to reveal environmental and genetic effects contributing to long-term HPA activity because of the complexity of chronic stress. The hair cortisol concentration (HCC) can be used to reflect the accumulation of HPA axis activity over time. Some studies suggest that the HCC might be associated with the protein concentration (PC) in the hair shaft; however, no studies have revealed a dynamic relationship between them. In the present study, 1,086 hair samples from 418 Holstein cows were collected, and the effects of environmental factors on HCC, PC, and ratio of HCC to PC (HCCP) were studied. Subsequently, regression analysis and curve fitting were used to identify for better-performing indicators of chronic stress. Additionally, univariate and bivariate genetic evaluation were used to estimate the genetic components of cortisol traits and genotype by environment interactions (G × E) under different environmental and physiological states. The results showed that HCC and PC are significantly affected by hair color, sampling year, and season, whereas HCCP is not influenced by hair color. Adjusted PC and HCCP, where confounding effects are excluded, were moderately related with chronic stress indicators. Moderate to high heritabilities were obtained for HCC (0.347 and 0.390 for winter and summer, respectively), PC (0.402 and 0.495 for winter and summer, respectively) and HCCP (0.289 and 0.460 for winter and summer, respectively) when animals in the same season were evaluated. A moderate G × E interaction was detected in this study, as indicated by the low or negative genetic correlation for the same cortisol trait in different environments (e.g. heat stress condition and thermoneutral condition). In conclusion, HCCP is not affected by hair color compared with the other 2 traits; thus, it has potential as an indicator of chronic stress. Hair cortisol traits could monitor stress response process in cattle, as well as provide a better understanding of genetic mechanism for long-term HPA activity.

Background The effect of heat stress on livestock production is a worldwide issue. Animal performance is influenced by exposure to harsh environmental conditions potentially causing genotype-by-environment interactions (G × E), especially in highproducing animals. In this context, the main objectives of this study were to (1) detect the time periods in which heifer fertility traits are more sensitive to the exposure to high environmental temperature and/or humidity, (2) investigate G × E due to heat stress in heifer fertility traits, and, (3) identify genomic regions associated with heifer fertility and heat tolerance in Holstein cattle. Results Phenotypic records for three heifer fertility traits (i.e., age at first calving, interval from first to last service, and conception rate at the first service) were collected, from 2005 to 2018, for 56,998 Holstein heifers raised in 15 herds in the Beijing area (China). By integrating environmental data, including hourly air temperature and relative humidity, the critical periods in which the heifers are more sensitive to heat stress were located in more than 30 days before the first service for age at first calving and interval from first to last service, or 10 days before and less than 60 days after the first service for conception rate. Using reaction norm models, significant G × E was detected for all three traits regarding both environmental gradients, proportion of days exceeding heat threshold, and minimum temperature-humidity index. Through single-step genome-wide association studies, PLAG1, AMHR2, SP1, KRT8, KRT18, MLH1, and EOMES were suggested as candidate genes for heifer fertility. The genes HCRTR1, AGRP, PC, and GUCY1B1 are strong candidates for association with heat tolerance. Conclusions The critical periods in which the reproductive performance of heifers is more sensitive to heat stress are trait-dependent. Thus, detailed analysis should be conducted to determine this particular period for other fertility traits. The considerable magnitude of G × E and sire re-ranking indicates the necessity to consider G × E in dairy cattle breeding schemes. This will enable selection of more heat-tolerant animals with high reproductive efficiency under harsh climatic conditions. Lastly, the candidate genes identified to be linked with response to heat stress provide a better understanding of the underlying biological mechanisms of heat tolerance in dairy cattle.

Previous studies reported the physical, transcriptome, and metabolome changes in in vitro acute heat-stressed (38 °C versus 43 °C for 2 h) bovine granulosa cells. Granulosa cells exhibited transient proliferation senescence, oxidative stress, an increased rate of apoptosis, and a decline in steroidogenic activity. In this study, we performed a joint integration and network analysis of metabolomic and transcriptomic data to further narrow down and elucidate the role of differentially expressed genes, important metabolites, and relevant cellular and metabolic pathways in acute heat-stressed granulosa cells. Among the significant (raw p-value < 0.05) metabolic pathways where metabolites and genes converged, this study found vitamin B6 metabolism, glycine, serine and threonine metabolism, phenylalanine metabolism, arginine biosynthesis, tryptophan metabolism, arginine and proline metabolism, histidine metabolism, and glyoxylate and dicarboxylate metabolism. Important significant convergent biological pathways included ABC transporters and protein digestion and absorption, while functional signaling pathways included cAMP, mTOR, and AMPK signaling pathways together with the ovarian steroidogenesis pathway. Among the cancer pathways, the most important pathway was the central carbon metabolism in cancer. Through multiple analysis queries, progesterone, serotonin, citric acid, pyridoxal, L-lysine, succinic acid, L-glutamine, L-leucine, L-threonine, L-tyrosine, vitamin B6, choline, and CYP1B1, MAOB, VEGFA, WNT11, AOX1, ADCY2, ICAM1, PYGM, SLC2A4, SLC16A3, HSD11B2, and NOS2 appeared to be important enriched metabolites and genes, respectively. These genes, metabolites, and metabolic, cellular, and cell signaling pathways comprehensively elucidate the mechanisms underlying the intricate fight between death and survival in acute heat-stressed bovine granulosa cells and essentially help further our understanding (and will help the future quest) of research in this direction.

Heat stress affects granulosa cells (GCs) and the ovarian follicular microenvironment, causing poor oocyte developmental competence and fertility. This study aimed to investigate the physical responses and global transcriptomic changes in bovine GCs to acute heat stress (43 °C for 2 h) in vitro. Heat-stressed GCs exhibited transient proliferation senescence and resumed proliferation at 48 h post-stress, while post-stress immediate culture-media change had a relatively positive effect on proliferation resumption. Increased accumulation of reactive oxygen species and apoptosis was observed in the heat-stress group. In spite of the upregulation of inflammatory (CYCS, TLR2, TLR4, IL6, etc.), pro-apoptotic (BAD, BAX, TNFSF9, MAP3K7, TNFRSF6B, FADD, TRADD, RIPK3, etc.) and caspase executioner genes (CASP3, CASP8, CASP9), antioxidants and anti-apoptotic genes (HMOX1, NOS2, CAT, SOD, BCL2L1, GPX4, etc.) were also upregulated in heat-stressed GCs. Progesterone and estrogen hormones, along with steroidogenic gene expression, declined significantly, in spite of the upregulation of genes involved in cholesterol synthesis. Out of 12,385 differentially expressed genes (DEGs), 330 significant DEGs (75 upregulated, 225 downregulated) were subjected to KEGG functional pathway annotation, gene ontology enrichment, STRING network analyses and manual querying of DEGs for meaningful molecular mechanisms. High inflammatory response was found to be responsible for oxidative-stress-mediated apoptosis of GCs and nodes towards the involvement of the NF-κB pathway and repression of the Nrf2 pathway. Downregulation of MDM4, TP53, PIDD1, PARP3, MAPK14 and MYC, and upregulation of STK26, STK33, TGFB2, CDKN1A and CDKN2A, at the interface of the MAPK and p53 signaling pathway, can be attributed to transient cellular senescence and apoptosis in GCs. The background working of the AMPK pathway through upregulation of AKT1, AMPK,SIRT1, PYGM, SLC2A4 and SERBP1 genes, and downregulation of PPARGCIA, IGF2,PPARA, SLC27A3, SLC16A3, TSC1/2, KCNJ2, KCNJ16, etc., evidence the repression of cellular transcriptional activity and energetic homeostasis modifications in response to heat stress. This study presents detailed responses of acute-heat-stressed GCs at physical, transcriptional and pathway levels and presents interesting insights into future studies regarding GC adaptation and their interaction with oocytes and the reproductive system at the ovarian level.

The gastrointestinal microbiota greatly affects the health status and production performance of bovines. Presently, many studies have used high-throughput sequencing methods to investigate the gastrointestinal microbiome in bovines. However, the microbiome profile of crossbred cattle across the whole gastrointestinal tract (GIT) has not been thoroughly reported. In this study, the digesta at ten regions (including the rumen, reticulum, omasum, abomasum, duodenum, jejunum, ileum, cecum, colon, and rectum) of the GIT were collected in three Simmental × Holstein crossbred heifers aged 17 months, and microbial DNA was extracted and amplified for sequencing of the V3–V4 regions of the 16S rRNA gene. Functional orthologs of the microbiota genome were predicted and analyzed. We found that samples were categorized into three groups (the stomach, small intestine, and large intestine) by principal coordinate analysis (PCoA) based on Bray–Curtis dissimilarity in both the bacterial composition and functional profile. Samples from small intestine had the lowest alpha diversity of bacteria composition and highest alpha diversity of the functional composition. Three groups of GIT regions were characterized by several microbiome features. The stomach was characterized by Bacteroidetes and Fibrobacteres at the phylum level, and KEGG pathways related to the metabolism of cofactors and vitamins, glycan biosynthesis, and metabolism were enriched in the stomach. The small intestine was characterized by Actinobacteria and Patescibacteria at the phylum level, and KEGG pathways related to xenobiotics biodegradation and metabolism were enriched in the small intestine. The large intestine featured Ruminococcaceae, Rikenellaceae, and Bacteroidacea at the family level, and KEGG pathways, including steroid hormone biosynthesis, linoleic acid metabolism, and cysteine and methionine metabolism were enriched in the large intestine. The results of the current study revealed the spatial heterogeneity of microbiota across the GIT in Simmental × Holstein crossbreeds and identified microbial biomarkers of different regions. The results can provide useful information for the study of the gastrointestinal microbiome in bovines.

Genomic evaluations are a method for improving the accuracy of breeding value estimation. This study aimed to compare estimates of genetic parameters and the accuracy of breeding values for wool traits in Merino sheep between pedigree-based best linear unbiased prediction (PBLUP) and single-step genomic best linear unbiased prediction (ssGBLUP) using Bayesian inference. Data were collected from 28,391 yearlings of Chinese Merino sheep (classified in 1992&ndash, 2018) at the Xinjiang Gonaisi Fine Wool Sheep-Breeding Farm, China. Subjectively-assessed wool traits, namely, spinning count (SC), crimp definition (CRIM), oil (OIL), and body size (BS), and objectively-measured traits, namely, fleece length (FL), greasy fleece weight (GFW), mean fiber diameter (MFD), crimp number (CN), and body weight pre-shearing (BWPS), were analyzed. The estimates of heritability for wool traits were low to moderate. The largest h2 values were observed for FL (0.277) and MFD (0.290) with ssGBLUP. The heritabilities estimated for wool traits with ssGBLUP were slightly higher than those obtained with PBLUP. The accuracies of breeding values were low to moderate, ranging from 0.362 to 0.573 for the whole population and from 0.318 to 0.676 for the genotyped subpopulation. The correlation between the estimated breeding values (EBVs) and genomic EBVs (GEBVs) ranged from 0.717 to 0.862 for the whole population, and the relative increase in accuracy when comparing EBVs with GEBVs ranged from 0.372% to 7.486% for these traits. However, in the genotyped population, the rank correlation between the estimates obtained with PBLUP and ssGBLUP was reduced to 0.525 to 0.769, with increases in average accuracy of 3.016% to 11.736% for the GEBVs in relation to the EBVs. Thus, genomic information could allow us to more accurately estimate the relationships between animals and improve estimates of heritability and the accuracy of breeding values by ssGBLUP.

Background Heat stress is known to affect follicular dynamics, oocyte maturation, and fertilization by impairing steroidogenic ability and viability of bovine granulosa cell (bGCs). The present study explored the physiological and molecular response of bGCs to different heat stress intensities in-vitro. We exposed the primary bGCs to heat stress (HS) at 39 °C, 40 °C and 41 °C along with control samples (38 °C) for 2 h. To evaluate the impact of heat stress on bGCs, several in vitro cellular parameters including cell apoptosis, intracellular reactive oxygen species (ROS) accumulation and HSP70 kinetics were assessed by flow cytometry, florescence microscopy and western blot, respectively. Furthermore, the ELISA was performed to confirm the 17β-estradiol (E2) and progesterone (P4) levels. In addition, the RNA sequencing (RNA-Seq) method was used to get the molecular based response of bGCs to different heat treatments. Results Our findings revealed that the HS significantly decreased the cell viability, E2 and P4 levels in bGCs, whereas, increased the cellular apoptosis and ROS. Moreover, the RNA-Seq experiments showed that all the treatments (39 °C, 40 °C and 41 °C) significantly regulated many differentially expressed genes (DEGs) i.e. BCL2L1, STAR, CYP11A1, CASP3, SOD2, HSPA13, and MAPK8IP1 and pathways associated with heat stress, apoptosis, steroidogenesis, and oxidative stress. Conclusively, our data demonstrated that the impact of 40 °C treatment was comparatively detrimental for cell viability, apoptosis and ROS accumulation. Notably, a similar trend of gene expression was reported by RT-qPCR for RNA-seq data. Conclusions Our study presented a worthy strategy for the first time to characterize the cellular and transcriptomic adaptation of bGCs to heat stress (39, 40 and 41 °C) in-vitro. The results infer that these genes and pathways reported in present study could be useful candidates/indicators for heat stress research in dairy cattle. Moreover, the established model of bGCs to heat stress in the current study provides an appropriate platform to understand the mechanism of how heat-stressed bGCs can affect the quality of oocytes and developing embryo.

Additional file 7: Covering letter.

Additional file 7: Covering letter.

Additional file 6: Figures.

Additional file 5: Table of gene primers for RT-qPCR.

The objective of this study was to investigate the role of liquorice extract (LE) as a dietary antioxidant supplement in the regulation of muscle antioxidant enzyme in lambs. Fifty Tan male lambs were randomly allocated among five groups with LE supplementation at levels of 0, 1000, 2000, 3000 and 4000 mg/kg feed, respectively. After 120 days of feeding, the lambs were slaughtered and the longissimus dorsi muscle was sampled for measuring the antioxidant enzymes activities and genes expression levels. In addition, superoxide anion (O 2 − ), hydrogen peroxide (H 2 O 2 ) and thiobarbituric acid reactive substance (TBARS) levels were also tested in our study. The results showed that supplementation with LE increased ( P 0.05) glutathione peroxidase (GPx) and catalase (CAT) activities, but LE had no effect on Cu/Zn superoxide dismutase (Cu/Zn SOD) and Mn superoxide dismutase (Mn SOD) activities. For gene expression, no significant differences were observed in mRNA levels of Cu/Zn SOD, Mn SOD, GPx and CAT genes among the five groups. In addition, LE supplementation decreased ( P 0.05) O 2 − , H 2 O 2 and TBARS levels in muscle. The results demonstrated that LE as a supplement could improve GPx and CAT activities, then prevent muscle against lipid oxidation.

Inflammaging is well understood in the study of humans; however, it is rarely reported for dairy cows. To understand the changing pattern of the gut microbiota, inflammatory status and milk production performance during the aging process in cows, we grouped 180 cows according to their lactation period: L1 (n = 60, 1st lactation), L3 (n = 60, 3rd lactation), and L5+ (n = 60, at least 5th lactation) and analyzed their milk components and daily milk yields to evaluate the changing pattern of milk production. The microbiota was analyzed using high-throughput sequencing of amplicons of 16S rRNA, which also allowed us to predict the functions of microbes and then study the changing pattern of the ruminal and fecal microbiota. Serum cytokines, including TNF-α, IL-6, IL-10, and TGF-β were measured to study the progress of inflammaging in the cows. We found that old cows (L5+) suffered from a long-term and low-level chronic inflammation, as indicated by significantly higher levels of inflammatory cytokines IL-10, TNF-α, and TGF-β in the L5+ group (p < 0.001). We also observed a significant decrease in daily milk yield and milk lactose, as well as a significant increase in somatic cell score, among the cows in the L5+ group. For the gut microbiota, most of the genera belonging to Prevotellaceae and Lachnospiraceae, which had a higher abundance among cows of both the L1 and L3 groups (LEfSe, LDA > 2), showed a similar change pattern during the aging process, both in the rumen and in feces, and across the six farms. Beneficial bacteria, like Bacteroidaceae, Eubacterium, and Bifidobacterium, displayed lower abundance in the feces of the L5+ group (LEfSe, LDA > 2). Reconstruction of the fecal bacteria community indicated transformation of the fermenting pattern of older cows’ (L5+) feces microbiota, with increased functions related the protein metabolism and fewer functions related to carbohydrate and lipid metabolism compared with those in L1 (p < 0.05). Finally, the connections among these changing patterns were revealed using redundancy analysis and network analysis. The results support the hypothesis of prolonging a cows’ productive life and improve dairy cow milk productive performances by manipulating the gut microbiota.

The skin has many important roles in dairy cattle, and skinfold thickness could be used as an indicator of body fat deposition. The objectives of this study were to estimate genetic parameters of skinfold thickness and to explore its association with body condition score (BCS) and milk production traits in a Chinese Holstein population. Skinfold thicknesses over the neck (STN) and the last rib (STR), BCS, and test-day records of milk production traits were available for 6,416 lactating Holstein cows in the summers of 2015 and 2016 in Beijing, China. Multi-trait animal models were used to estimate variance and covariance components using the DMU software. The average STN was 7.15 ± 1.28 mm, and the average STR was 11.76 ± 1.95 mm (mean ± standard deviation). Estimated heritability was 0.13 ± 0.03 for STN and 0.26 ± 0.04 for STR. We detected a high genetic correlation (0.79 ± 0.08; heritability ± standard error) between STN and STR. Genetic correlations between skinfold thickness and BCS were low to moderate: 0.18 between STR and BCS, and 0.33 between STN and BCS. Genetic correlations between skinfold thickness and milk yield, milk fat percentage, and milk protein percentage were negligible, ranging from −0.02 to 0.15. Collectively, skinfold thickness is characterized as a trait with moderate heritability. Skinfold thickness is sensitive to changes in body condition or fat deposition across parities and lactation stages in milking cows, and we confirmed the complementary nature of skinfold thickness and BCS genetically as well as phenotypically by comparing their changing trends throughout lactation and across lactations. The use of skinfold thickness, together with BCS, can assist in the monitoring of changes in body fat deposition to achieve higher management precision.

The rates of mortality and involuntary culling of dairy calves and replacement heifers have great economic implications on the dairy cattle industry around the world. The main objectives of this study were: (1) to obtain population parameters of mortality and involuntary culling rates of dairy calves and replacement heifers, and, (2) to investigate the factors affecting mortality and involuntary culling rates in Chinese Holstein cattle. Two datasets containing records of birth, calving, and culling events from 142,833 Holstein cattle born between 1991 and 2018 were used in this study. The population parameters were obtained using dataset 1, which consisted of dairy calves and replacement heifers that died or were involuntarily culled. Three survival traits were defined in dataset 2, which consisted of females born from 1999 to 2018. A binomial logistic regression was used to analyze the risk factors on the survival traits. The mortality rate of dairy calves and replacement heifers from day 3 to 60, 61 to 365, and 366 to first calving was 5.5%, 7.4%, and 8.7%, and an unfavorable increasing trend was observed. Health events associated with digestive and respiratory or circulatory systems were the main death reasons. Herd-birth year, birth season, and dam parity had significant effects on survival traits. The results from this study will help farmers to better manage calves and replacement heifers and highlight the need to include survival traits in dairy calves and replacement heifers as part of national genetic evaluation schemes.

To investigate the potential use of n-alkanes (alkanes), long-chain alcohols (alcohols) and long-chain fatty acids (acids) for estimating the diet composition of sheep, in a feeding trial. A total of 18 sheep were assigned randomly to three different diets (diet A, diet B and diet C) containing up to eight herbage species (Leymus chinensis, Leymus dasystachys, Elymus sibiricum, Chenopodium album, Puccinellia chinampoensis, Medicago sativa, Saussurea sinuata and Bromus inermis). Faecal recoveries of alkanes, alcohols and acids were determined, and diet compositions were estimated using different combinations of alkanes, alcohols and acids. The faecal concentrations of individual alkanes, alcohols and acids were corrected using the mean recovery of the dietary treatment that the respective animal belonged to (diet recovery), or the mean recovery across all dietary treatments (general recovery). In general, diets did not affect the faecal recovery values for alkanes, alcohols and acids, and no difference in accuracy was found between diet composition estimates based on dietary recovery and general recovery. The accuracy of diet composition estimates declined as the number of dietary components increased from four to eight herbage species (P < 0.001). Better (P < 0.05) estimates of diet composition were obtained with the combinations of two or three marker types instead of alkanes alone. Moreover, results showed that excluding minor diet components from the calculations decreased (P < 0.05) the accuracy of diet composition estimates, whereas including extra non-grazed herbage species did not reduce (P > 0.05) the quality of diet composition estimates. These results confirmed the usefulness of alkanes, alcohols and acids as markers for determining complex diet composition of sheep. However, a negative impact on the accuracy of diet composition estimates, caused by missing minor diet components from the calculation of diet composition, could happen when plant wax markers are used to estimate the diet composition of free-ranging animals.

The limitations of the alkane technique in estimating the diet components of herbivores call for the introduction of new diet composition markers. Recently, long-chain alcohols (alcohols) and long-chain fatty acids (acids) have received the most attention and show great potential, when combined with alkanes, to estimate composition of complex diets. In the current study, faecal recoveries of alcohols and acids were determined in sheep in four different live weight groups fed three herbage species, either Leymus chinensis, L. dasystachys or Elymus sibiricum. Analysis of variance (ANOVA) was used to examine the effects of herbage species and live weight of sheep on faecal recoveries of individual alcohols and acids. Further, an indoor experiment with six sheep fed a diet of equal proportions, on dry matter (DM) basis, of three herbages was performed, allowing to assess the accuracy of alcohols and/or acids in combination with alkanes, to estimate diet composition. A one-sample t-test was carried out to test the accuracy of these estimates. Results of the first experiment indicated that the faecal recoveries of alcohols and acids were significantly affected by herbage species (P < 0.05). While the effects were significant or near significant for the faecal recoveries of some alcohols (C24-ol, C30-ol and C26-ol) (P ⩽ 0.05), no effect of live weight on faecal recoveries of acids was observed (P > 0.05). Therefore, adjustments based on diet-specific faecal recoveries might improve diet composition estimates. This was illustrated by the results of the second experiment. The diet composition estimated from alcohols or all combinations of alcohols with other marker types, after diet-specific correction of faecal recoveries, did not significantly differ from the actual composition (P > 0.05). However, using acids as additional markers resulted in poorer diet composition estimates. This study confirmed the utility of alcohols, combined with alkanes, as markers to estimate composition of complex diets. Although corrections based on mean faecal recoveries, average over animals and diets, resulted in some accuracy loss, results were still satisfactory and better than without recovery correction.