Background: Characterization 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.Results: We 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.Conclusions: Our 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.

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

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.

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.

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.