RNA-seq profiling of skin in temperate and tropical cattle

Skin transcriptome analysis in Brangus cattle under heat stress

Heat stress is a major factor that negatively affects animal welfare and production systems. Livestock should adapt to tropical and subtropical areas and to meet this, composite breeds have been developed. This work aimed to evaluate gene expression profiles in the skin of Brangus cattle under heat stress using a case-control design, and to correlate this with skin histological characteristics. Two groups of bulls were set using rectal temperature as a criterion to define stress conditions: stressed (N = 5) and non-stressed (N = 5) groups. Skin transcriptomics was performed and correlations between breed composition, phenotypic and skin histological traits were evaluated. Results showed 4309 differentially expressed genes (P < 0.01), 2113 downregulated and 2196 upregulated. Enrichment and ontology analyses revealed 132 GO terms and 67 pathways (P < 0.01), including thermogenesis, glycolysis, gluconeogenesis, mitochondrial activity, antioxidant and immune response, and apoptosis. The identity of the terms and pathways indicated the diversity of mechanisms directed to relieve the animals' suffering, acting from simple passive mechanisms (conduction, convection and radiation) to more complex active ones (behavioural changes, evaporation, vasodilation and wheezing). Furthermore, significant differences between phenotypic and skin histological traits and correlations between pairs of traits suggested a direction towards heat dissipation processes. In this sense, number of vessels was positively correlated with number of sweat glands (P < 0.001) and both were positively correlated with zebuine genetic content (P < 0.05 and P < 0.01, respectively), gland size was positively correlated with epidermal thickness and negatively with hair length (P < 0.05), and epidermal thickness was negatively correlated with gland-epidermis distance (P < 0.0005). These results support the notion that response to heat stress is physiologically complex, producing significant changes in the expression of genes involved in several biological pathways, while the animal's ability to face it depends greatly on their skin features.

Multi-Omics Approaches to Improve Meat Quality and Taste Characteristics

With rapid advances in meat science in recent decades, changes in meat quality during the pre-slaughter phase of muscle growth and the post-slaughter process from muscle to meat have been investigated. Commonly used techniques have evolved from early physicochemical indicators such as meat color, tenderness, water holding capacity, flavor, and pH to various omic tools such as genomics, transcriptomics, proteomics, and metabolomics to explore fundamental molecular mechanisms and screen biomarkers related to meat quality and taste characteristics. This review highlights the application of omics and integrated multi-omics in meat quality and taste characteristics studies. It also discusses challenges and future perspectives of multi-omics technology to improve meat quality and taste. Consequently, multi-omics techniques can elucidate the molecular mechanisms responsible for changes of meat quality at transcriptome, proteome, and metabolome levels. In addition, the application of multi-omics technology has great potential for exploring and identifying biomarkers for meat quality and quality control that can make it easier to optimize production processes in the meat industry.

Heat stress effects on milk yield traits and metabolites and mitigation strategies for dairy cattle breeds reared in tropical and sub-tropical countries

Heat stress is an important problem for dairy industry in many parts of the world owing to its adverse effects on productivity and profitability. Heat stress in dairy cattle is caused by an increase in core body temperature, which affects the fat production in the mammary gland. It reduces milk yield, dry matter intake, and alters the milk composition, such as fat, protein, lactose, and solids-not-fats percentages among others. Understanding the biological mechanisms of climatic adaptation, identifying and exploring signatures of selection, genomic diversity and identification of candidate genes for heat tolerance within indicine and taurine dairy breeds is an important progression toward breeding better dairy cattle adapted to changing climatic conditions of the tropics. Identifying breeds that are heat tolerant and their use in genetic improvement programs is crucial for improving dairy cattle productivity and profitability in the tropics. Genetic improvement for heat tolerance requires availability of genetic parameters, but these genetic parameters are currently missing in many tropical countries. In this article, we reviewed the HS effects on dairy cattle with regard to (1) physiological parameters; (2) milk yield and composition traits; and (3) milk and blood metabolites for dairy cattle reared in tropical countries. In addition, mitigation strategies such as physical modification of environment, nutritional, and genetic development of heat tolerant dairy cattle to prevent the adverse effects of HS on dairy cattle are discussed. In tropical climates, a more and cost-effective strategy to overcome HS effects is to genetically select more adaptable and heat tolerant breeds, use of crossbred animals for milk production, i.e., crosses between indicine breeds such as Gir, white fulani, N'Dama, Sahiwal or Boran to taurine breeds such as Holstein-Friesian, Jersey or Brown Swiss. The results of this review will contribute to policy formulations with regard to strategies for mitigating the effects of HS on dairy cattle in tropical countries.

Novel Insights to Assess Climate Resilience in Goats Using a Holistic Approach of Skin-Based Advanced NGS Technologies

A novel study was conducted to elucidate heat-stress responses on a number of hair- and skin-based traits in two indigenous goat breeds using a holistic approach that considered a number of phenotypic and genomic variables. The two goat breeds, Kanni Aadu and Kodi Aadu, were subjected to a simulated heat-stress study using the climate chambers. Four groups consisting of six goats each (KAC, Kaani Aadu control; KAH, Kanni Aadu heat stress; KOC, Kodi Aadu control; and KOH, Kodi Aadu heat stress) were considered for the study. The impact of heat stress on caprine skin tissue along with a comparative assessment of the thermal resilience of the two goat breeds was assessed. The variables considered were hair characteristics, hair cortisol, hair follicle quantitative PCR (qPCR), sweating (sweating rate and active sweat gland measurement), skin histometry, skin-surface infrared thermography (IRT), skin 16S rRNA V3-V4 metagenomics, skin transcriptomics, and skin bisulfite sequencing. Heat stress significantly influenced the hair fiber characteristics (fiber length) and hair follicle qPCR profile (Heat-shock protein 70 (HSP70), HSP90, and HSP110). Significantly higher sweating rate, activated sweat gland number, skin epithelium, and sweat gland number (histometry) were observed in heat stressed goats. The skin microbiota was also observed to be significantly altered due to heat stress, with a relatively higher alteration being noticed in Kanni Aadu goats than in Kodi Aadi goats. Furthermore, the transcriptomics and epigenetics analysis also pointed towards the significant impact of heat stress at the cellular and molecular levels in caprine skin tissue. The higher proportion of differentially expressed genes (DEGs) along with higher differentially methylated regions (DMRs) in Kanni Aadu goats due to heat stress when compared to Kodi Aadu goats pointed towards the better resilience of the latter breed. A number of established skin, adaptation, and immune-response genes were also observed to be significantly expressed/methylated. Additionally, the influence of heat stress at the genomic level was also predicted to result in significant functional alterations. This novel study thereby highlights the impact of heat stress on the caprine skin tissue and also the difference in thermal resilience exhibited by the two indigenous goat breeds, with Kodi Aadu goats being more resilient.

Effects of the bovine SLICK1 mutation in PRLR on sweat gland area, FOXA1 abundance, and global gene expression in skin

The SLICK1 mutation in the prolactin receptor (PRLR) results in a short-hair coat and increased ability to regulate body temperature during heat stress. It is unclear whether the mutation affects capacity for sweating. The objective of this observational study was to evaluate whether the SLICK1 mutation in PRLR alters characteristics of skin related to sweat gland abundance or function. Skin biopsies from 31 Holstein heifers, including 14 wild-type (SL^-/-) and 17 heterozygous slick (SL^+/-), were subjected to histological analysis to determine the percent of the surface area of skin sections that are occupied by sweat glands. We detected no effect of genotype on this variable. Immunohistochemical analysis of the forkhead transcription factor A1 (FOXA1), a protein essential for sweating in mice, from 6 SL^-/- and 6 SL^+/- heifers indicated twice as much FOXA1 in sweat glandular epithelia of SL^+/- heifers as in SL^-/- heifers. Results from RNA sequencing of skin biopsies from 5 SL^-/- and 7 SL^+/- heifers revealed few genes that were differentially expressed and none that have been associated with sweat gland development or function. In conclusion, results do not support the idea that the SLICK1 mutation changes the abundance of sweat glands in skin, but do show that functional properties of sweat glands, as indicated by increased abundance of immunoreactive FOXA1, are modified by inheritance of the mutation in PRLR.

Transcriptome-wide analysis reveals gluten-induced suppression of small intestine development in young chickens

Wheat gluten is an increasingly common ingredient in poultry diets but its impact on the small intestine in chicken is not fully understood. This study aimed to identify effects of high-gluten diets on chicken small intestines and the variation of their associated transcriptional responses by age. A total of 120 broilers (Ross Strain) were used to perform two animal experiments consisting of two gluten inclusion levels (0% or 25%) by bird’s age (1 week or 4 weeks). Transcriptomics and histochemical techniques were employed to study the effect of gluten on their duodenal mucosa using randomly selected 12 broilers (3 chicks per group). A reduction in feed intake and body weight gain was found in the broilers fed a high-gluten containing diet at both ages. Histochemical photomicrographs showed a reduced villus height to crypt depth ratio in the duodenum of gluten-fed broilers at 1 week. We found mainly a significant effect on the gene expression of duodenal mucosa in gluten-fed broilers at 1 week (289 differentially expressed genes [DEGs]). Pathway analyses revealed that the significant DEGs were mainly involved in ribosome, oxidative phosphorylation, and peroxisome proliferator-activated receptor (PPAR) signaling pathways. These pathways are involved in ribosome protein biogenesis, oxidative phosphorylation and fatty acid metabolism, respectively. Our results suggest a pattern of differential gene expression in these pathways that can be linked to chronic inflammation, suppression of cell proliferation, cell cycle arrest and apoptosis. And via such a mode of action, high-gluten inclusion levels in poultry diets could lead to the observed retardation of villi development in the duodenal mucosa of young broiler chicken.

An analysis of skin thickness in the Dezhou donkey population and identification of candidate genes by RNA-seq

The aim of the present study was to analyze the main factors that have a significant impact on skin thickness, and to further identify the genes and signaling pathways regulating skin growth by RNA-seq in Dezhou donkeys. Skin samples from different body regions of 15 slaughtered donkeys were obtained to study variations in skin thickness over the bodies. Skin thickness data for another 514 donkeys was obtained by minimally invasive skin sampling from the back, and measurements of the donkeys' body size traits and pedigree data were also collected. These data were used to analyze changes in skin thickness and estimate genetic parameters. In addition, transcriptomic analysis was conducted on the skin tissues of individuals from two groups with significant differences in skin thickness. Our results showed that skin thickness over the bodies ranged from 1.08 to 4.36 mm. The skin from the back was the thickest and had the highest correlation with that of other regions of the body. The skin thickness decreased from the back to the side of the ventral abdomen, and the skin thickness on the limbs increased from the proximal end to the distal end. The results also showed that the skin from the same body regions of jacks was thicker than that of jennies in the same age group. The skin thickness of jennies increased from birth to the age of 2 and then clearly decreased after 2 years of age. The estimated heritability of skin thickness was 0.15, and the genetic correlations between skin thickness and body size traits were negligible. Transcriptome analysis showed that the thick-skin group had 65 up-regulated genes and 38 down-regulated genes compared with the thin-skin group. The differentially expressed genes were highly enriched in epidermal development and cell adhesion molecule signaling pathways. We identified the candidate genes responsible for variations in skin thickness in the Dezhou donkey, including KRT10, KRT1, CLDN9, MHCII and MMP28. These results contribute to a better understanding of the growth and development of donkey skin, reveal the molecular mechanism responsible for donkey skin thickness and suggest directions for genetic selection in the Dezhou donkey population.

Historical Introgression from Wild Relatives Enhanced Climatic Adaptation and Resistance to Pneumonia in Sheep

How animals, particularly livestock, adapt to various climates and environments over short evolutionary time is of fundamental biological interest. Further, understanding the genetic mechanisms of adaptation in indigenous livestock populations is important for designing appropriate breeding programs to cope with the impacts of changing climate. Here, we conducted a comprehensive genomic analysis of diversity, interspecies introgression, and climate-mediated selective signatures in a global sample of sheep and their wild relatives. By examining 600K and 50K genome-wide single nucleotide polymorphism data from 3,447 samples representing 111 domestic sheep populations and 403 samples from all their seven wild relatives (argali, Asiatic mouflon, European mouflon, urial, snow sheep, bighorn, and thinhorn sheep), coupled with 88 whole-genome sequences, we detected clear signals of common introgression from wild relatives into sympatric domestic populations, thereby increasing their genomic diversities. The introgressions provided beneficial genetic variants in native populations, which were significantly associated with local climatic adaptation. We observed common introgression signals of alleles in olfactory-related genes (e.g., ADCY3 and TRPV1) and the PADI gene family including in particular PADI2, which is associated with antibacterial innate immunity. Further analyses of whole-genome sequences showed that the introgressed alleles in a specific region of PADI2 (chr2: 248,302,667-248,306,614) correlate with resistance to pneumonia. We conclude that wild introgression enhanced climatic adaptation and resistance to pneumonia in sheep. This has enabled them to adapt to varying climatic and environmental conditions after domestication.