This study aimed to investigate the genes and pathways that respond to heat stress in Holstein bull calves exposed to severe ranges of temperature and humidity. A total of ten animals from 4 to 6 months of age were subjected to heat stress at 37°C and 90 % humidity for 12 h. Skin and rectal temperatures were measured before and after heat stress; while no correlation was found between them before heat stress, a moderate correlation was detected after heat stress, confirming rectal temperature to be a better barometer for monitoring heat stress. RNAseq analysis identified 8567 genes to be differentially regulated, out of which 465 genes were significantly upregulated (≥2-fold, P < 0.05) and 49 genes were significantly downregulated (≤2-fold, P < 0.05) in response to heat stress. Significant terms and pathways enriched in response to heat stress included chaperones, cochaperones, cellular response to heat stress, phosphorylation, kinase activation, immune response, apoptosis, Toll-like receptor signaling pathway, Pi3K/AKT activation, protein processing in endoplasmic reticulum, interferon signaling, pathways in cancer, estrogen signaling pathway, and MAPK signaling pathway. The differentially expressed genes were validated by quantitative realtime PCR analysis, which confirmed the tendency of the expression. The genes and pathways identified in this analysis extend our understanding of transcriptional response to heat stress and their likely functioning in adapting the animal to hyperthermic stress. The identified genes could be used as candidate genes for association studies to select and breed animals with improved heat tolerance.

Non-synonymous SNPs and protein coding SNPs within the promoter region of genes (regulatory SNPs) might have a significant effect on carcass traits. Imputed sequence level data of 10,215 Hanwoo bulls, annotated and filtered to include only regulatory SNPs (450,062 SNPs), were used in a genome-wide association study (GWAS) to identify loci associated with backfat thickness (BFT), carcass weight (CWT), eye muscle area (EMA), and marbling score (MS). A total of 15, 176, and 1 SNPs were found to be significantly associated (p < 1.11 × 10 −7 ) with BFT, CWT, and EMA, respectively. The significant loci were BTA4 (CWT), BTA6 (CWT), BTA14 (CWT and EMA), and BTA19 (BFT). BayesR estimated that 1.1%~1.9% of the SNPs contributed to more than 0.01% of the phenotypic variance. So, the GWAS was complemented by a gene-set enrichment (GSEA) and protein-protein interaction network (PPIN) analysis in identifying the pathways affecting carcass traits. At p < 0.005 (~2,261 SNPs), 25 GO and 18 KEGG categories, including calcium signaling, cell proliferation, and folate biosynthesis, were found to be enriched through GSEA. The PPIN analysis showed enrichment for 81 candidate genes involved in various pathways, including the PI3K-AKT, calcium, and FoxO signaling pathways. Our finding provides insight into the effects of regulatory SNPs on carcass traits.Genes 2020, 11, 316 2 of 22 and eye muscle area (EMA) [4]. Though substantial improvement in carcass and meat quality have been achieved, due to market requirement for higher quality, and for improving the economic value of Hanwoo, continuous improvement of economically important trait is required [5,6]. A genome-wide association study (GWAS) is an affordable and powerful tool to discover candidate genes and loci associated with quantitative traits [7]. GWASs in livestock, including in Hanwoo [8][9][10][11], have resulted in remarkable insights into the genetic architecture of carcass traits. Genetic variation in complex traits such as carcass and meat quality traits are, however, due to the contribution of many mutations with small effects [1,12] (polygenic effect). Though some of these mutations have been successfully identified through GWASs, the high significance thresholds required to correct for the multiple testing problem results in the identification of only SNPs with a large effect size [7,12]. Further, a GWAS does not make use of the fact that genes work together in a network, and multi-allelic QTL might not be captured due to the bi-allelic nature of SNPs [13]. Moreover, epistasis is an important genetic component underlying phenotypic variation that also accounts for missing heritability [14]. Therefore, a GWAS alone might result in only limited understanding of the nature of complex traits [13]. Suggested solutions to overcome this limitation and understand the genetic complexities regulating complex traits are to complement a GWAS with gene-set enrichment, a protein-protein interaction network (PPIN), and pathway analyses [15][16][17][18]. In GSEA and pathway analysis, a g...

This study shows that microperimetry can be a useful tool for objective evaluation of macular function and progression of the disease.

No abstract

Whole-genome sequence (WGS) data are increasingly being applied into genomic predictions, offering a higher predictive ability by including causal mutations or single-nucleotide polymorphisms (SNPs) putatively in strong linkage disequilibrium with causal mutations affecting the trait. This study aimed to improve the predictive performance of the customized Hanwoo 50 k SNP panel for four carcass traits in commercial Hanwoo population by adding highly predictive variants from sequence data. A total of 16,892 Hanwoo cattle with phenotypes (i.e., backfat thickness, carcass weight, longissimus muscle area, and marbling score), 50 k genotypes, and WGS imputed genotypes were used. We partitioned imputed WGS data according to functional annotation [intergenic (IGR), intron (ITR), regulatory (REG), synonymous (SYN), and non-synonymous (NSY)] to characterize the genomic regions that will deliver higher predictive power for the traits investigated. Animals were assigned into two groups, the discovery set (7324 animals) used for predictive variant detection and the cross-validation set for genomic prediction. Genome-wide association studies were performed by trait to every genomic region and entire WGS data for the pre-selection of variants. Each set of pre-selected SNPs with different density (1000, 3000, 5000, or 10,000) were added to the 50 k genotypes separately and the predictive performance of each set of genotypes was assessed using the genomic best linear unbiased prediction (GBLUP). Results showed that the predictive performance of the customized Hanwoo 50 k SNP panel can be improved by the addition of pre-selected variants from the WGS data, particularly 3000 variants from each trait, which is then sufficient to improve the prediction accuracy for all traits. When 12,000 pre-selected variants (3000 variants from each trait) were added to the 50 k genotypes, the prediction accuracies increased by 9.9, 9.2, 6.4, and 4.7% for backfat thickness, carcass weight, longissimus muscle area, and marbling score compared to the regular 50 k SNP panel, respectively. In terms of prediction bias, regression coefficients for all sets of genotypes in all traits were close to 1, indicating an unbiased prediction. The strategy used to select variants based on functional annotation did not show a clear advantage compared to using whole-genome. Nonetheless, such pre-selected SNPs from the IGR region gave the highest improvement in prediction accuracy among genomic regions and the values were close to those obtained using the WGS data for all traits. We concluded that additional gain in prediction accuracy when using pre-selected variants appears to be trait-dependent, and using WGS data remained more accurate compared to using a specific genomic region.

Heat stress (HS) negatively affects chicken performance. Agricultural expansion will happen in regions that experience high ambient temperatures, where fast-growing commercial chickens are vulnerable. Indigenous chickens of such regions, due to generations of exposure to environmental challenges, might have higher thermal tolerance. In this study, two indigenous chicken ecotypes, from the hot and humid Mombasa (lowland) and the colder Naivasha (highland) regions, were used to investigate the effects of acute (5 h, 35°C) and chronic (3 days of 35°C for 8 h/day) HS on the cardiac and skeletal muscle, through RNA sequencing. The rectal temperature gain and the number of differentially expressed genes (DEGs) [False Discovery Rate (FDR) < 0.05] were two times higher in the acute stage than in the chronic stage in both ecotypes, suggesting that cyclic exposure to HS can lead to adaptation. A tissue- and stage-specific difference in response to HS was observed, with peroxisome proliferator-activated-receptor (PPAR) signaling and mitogen-activate protein kinase (MAPK) signaling pathways, enriched in heart and skeletal muscle, respectively, and the p53 pathway enriched only in the acute stage in both tissues. The acute and chronic stage DEGs were integrated by a region-specific gene coexpression network (GCN), and genes with the highest number of connections (hub genes) were identified. The hub genes in the lowland network were CCNB2, Crb2, CHST9, SESN1, and NR4A3, while COMMD4, TTC32, H1F0, ACYP1, and RPS28 were the hub genes in the highland network. Pathway analysis of genes in the GCN showed that p53 and PPAR signaling pathways were enriched in both low and highland networks, while MAPK signaling and protein processing in endoplasmic reticulum were enriched only in the gene network of highland chickens. This shows that to dissipate the accumulated heat, to reduce heat induced apoptosis, and to promote DNA damage repair, the ecotypes activated or suppressed different genes, indicating the differences in thermal tolerance and HS response mechanisms between the ecotypes. This study provides information on the HS response of chickens, adapted to two different agro climatic environments, extending our understanding of the mechanisms of HS response and the effect of adaptation in counteracting HS.

No abstract

RNA-Seq analysis was used to characterize transcriptome response of Holstein calves to thermal stress. A total of eight animals aged between 2 and 3 months were randomly selected and subjected to thermal stress corresponding to a temperature humidity index of 95 in an environmentally controlled house for 12 h consecutively for 3 days. A set of 15,787 unigenes were found to be expressed and after a threshold of threefold change, and a Q value <0.05; 502, 394, and 376 genes were found to be differentially expressed on days 1, 2, and 3 out of which 343, 261 and 256 genes were upregulated and 159, 133, and 120 genes were downregulated. Only 356 genes out of these were expressed on all 3 days, and only they were considered as significantly differentially expressed. KEGG pathway analysis revealed that ten pathways were significantly enriched; the top two among them were protein processing in endoplasmic reticulum and MAPK signaling pathways. These results suggest that thermal stress triggered a complex response in Holstein calves and the animals adjusted their physiological and metabolic processes to survive. Many of the genes identified in this study have not been previously reported to be involved in thermal stress response. The results of this study extend our understanding of the animal's response to thermal stress and some of the identified genes may prove useful in the efforts to breed Holstein cattle with superior thermotolerance, which might help in minimizing production loss due to thermal stress.