Identification of the acclimation genes in transcriptomic responses to heat stress of White Pekin duck

Using a metabolomics approach to investigate the sensitivity of a potential Arctic-invader and its Arctic sister-species to marine heatwaves and traditional harvesting disturbances

Coastal species are threatened by fishing practices and changing environmental conditions, such as marine heatwaves (MHW). The mechanisms that confer tolerance to such stressors in marine invertebrates are poorly understood. However, differences in tolerance among different species may be attributed to their geographical distribution. To test the tolerance of species occupying different thermal ranges, we used two closely related bivalves the softshell clam Mya arenaria (Linnaeus, 1758), a cold-temperate invader with demonstrated potential for establishment in the Arctic, and the blunt gaper Mya truncata (Linnaeus, 1758), a native polar species. Clams were subjected to a thermal stress, mimicking a MHW, and harvesting stress in a controlled environment. Seven acute temperature changes (2, 7, 12, 17, 22, 27, and 32 °C) were tested at two harvesting disturbance intensities (with, without). Survival was measured after 12 days and three tissues (gills, mantle, and posterior adductor muscle) collected from surviving individuals for targeted metabolomic profiling. MHW tolerance differed significantly between species: 26.9 °C for M. arenaria and 17.8 °C for M. truncata, with a negligeable effect of harvesting. At the upper thermal limit, M. arenaria displayed a more profound metabolomic remodelling when compared to M. truncata, and this varied greatly between tissue types. Network analysis revealed differences in pathway utilization at the upper MHW limit, with M. arenaria displaying a greater reliance on multiple DNA repair and expression and cell signalling pathways, while M. truncata was limited to fewer pathways. This suggests that M. truncata is ill equipped to cope with warming environments. MHW patterning in the Northwest Atlantic may be a strong predictor of population survival and future range shifts in these two clam species. As polar environments undergo faster rates of warming compared to the global average, M. truncata may be outcompeted by M. arenaria expanding into its native range.

Longitudinal genomic analyses of automatically-recorded vaginal temperature in lactating sows under heat stress conditions based on random regression models

BACKGROUND

Automatic and continuous recording of vaginal temperature (TV) using wearable sensors causes minimal disruptions to animal behavior and can generate data that enable the evaluation of temporal body temperature variation under heat stress (HS) conditions. However, the genetic basis of TV in lactating sows from a longitudinal perspective is still unknown. The objectives of this study were to define statistical models and estimate genetic parameters for TV in lactating sows using random regression models, and identify genomic regions and candidate genes associated with HS indicators derived from automatically-recorded TV.

RESULTS

Heritability estimates for TV ranged from 0.14 to 0.20 over time (throughout the day and measurement period) and from 0.09 to 0.18 along environmental gradients (EG, - 3.5 to 2.2, which correspond to dew point values from 14.87 to 28.19 ˚C). Repeatability estimates of TV over time and along EG ranged from 0.57 to 0.66 and from 0.54 to 0.77, respectively. TV measured from 12h00 to 16h00 had moderately high estimates of heritability (0.20) and repeatability (0.64), indicating that this period might be the most suitable for recording TV for genetic selection purposes. Significant genotype-by-environment interactions (GxE) were observed and the moderately high estimates of genetic correlations between pairs of extreme EG indicate potential re-ranking of selection candidates across EG. Two important genomic regions on chromosomes 10 (59.370-59.998 Mb) and16 (21.548-21.966 Mb) were identified. These regions harbor the genes CDC123, CAMK1d, SEC61A2, and NUDT5 that are associated with immunity, protein transport, and energy metabolism. Across the four time-periods, respectively 12, 13, 16, and 10 associated genomic regions across 14 chromosomes were identified for TV. For the three EG classes, respectively 18, 15, and 14 associated genomic windows were identified for TV, respectively. Each time-period and EG class had uniquely enriched genes with identified specific biological functions, including regulation of the nervous system, metabolism and hormone production.

CONCLUSIONS

TV is a heritable trait with substantial additive genetic variation and represents a promising indicator trait to select pigs for improved heat tolerance. Moderate GxE for TV exist, indicating potential re-ranking of selection candidates across EG. TV is a highly polygenic trait regulated by a complex interplay of physiological, cellular and behavioral mechanisms.

Genetic diversity, population structure, and selective signature of sheep in the northeastern Tarim Basin

Local sheep in the northeastern Tarim Basin can adapt to dry and low-rainfall regional environments. In this study, three local sheep breeds in the northeastern Tarim Basin, LOP (LOP) sheep, Bayinbuluke (BYK) sheep, and Kunlun (KUN, also known as the Qiemo sheep) sheep, and three introduced sheep breeds, Suffolk (SUF) sheep, Dorset (APD) sheep, and Texel (TEX) sheep, were analyzed for genetic diversity, population structure, and selective signature using the Illumina OvineSNP50K BeadChip. We found that LOP, BYK, and KUN had lower observed heterozygosity and expected heterozygosity than TEX, SUF, and ADP, which were differentiated based on geographic distribution. We performed fixation index (FST) analysis on three local sheep breeds in the northeastern Tarim Basin (LOP, BYK, and KUN) and introduced sheep breeds (TEX, SUF, and ADP) to measure genetic differentiation. Nucleotide diversity (PI) analysis was performed on single-nucleotide polymorphism (SNP) data of LOP, BYK, and KUN. A total of 493 candidate genes were obtained by taking the intersection at a threshold of 5%. Among them, SMAD2, ESR2, and HAS2 were related to reproductive traits. PCDH15, TLE4, and TFAP2B were related to growth traits. SOD1, TSHR, and DNAJB5 were related to desert environmental adaptation. Analyzing the genetic patterns of local sheep in the northeastern Tarim Basin can protect the germplasm resources of local sheep and promote the development and utilization of sheep genetic resources.

Identification of Runs of Homozygosity Islands and Functional Variants in Wenchang Chicken

Wenchang chickens, a native breed in the Hainan province of China, are famous for their meat quality and adaptability to tropical conditions. For effective management and conservation, in the present study, we systematically investigated the characteristics of genetic variations and runs of homozygosity (ROH) along the genome using re-sequenced whole-genome sequencing data from 235 Wenchang chickens. A total of 16,511,769 single nucleotide polymorphisms (SNPs) and 53,506 ROH segments were identified in all individuals, and the ROH of Wenchang chicken were mainly composed of short segments (0-1 megabases (Mb)). On average, 5.664% of the genome was located in ROH segments across the Wenchang chicken samples. According to several parameters, the genetic diversity of the Wenchang chicken was relatively high. The average inbreeding coefficient of Wenchang chickens based on F_HOM, F_GRM, and F_ROH was 0.060 ± 0.014, 0.561 ± 0.020, and 0.0566 ± 0.01, respectively. A total of 19 ROH islands containing 393 genes were detected on 9 different autosomes. Some of these genes were putatively associated with growth performance (AMY1a), stress resistance (THEMIS2, PIK3C2B), meat traits (MBTPS1, DLK1, and EPS8L2), and fat deposition (LANCL2, PPARγ). These findings provide a better understanding of the degree of inbreeding in Wenchang chickens and the hereditary basis of the characteristics shaped under selection. These results are valuable for the future breeding, conservation, and utilization of Wenchang and other chicken breeds.

Improving Genomic Selection for Heat Tolerance in Dairy Cattle: Current Opportunities and Future Directions

Heat tolerance is the ability of an animal to maintain production and reproduction levels under hot and humid conditions and is now a trait of economic relevance in dairy systems worldwide because of an escalating warming climate. The Australian dairy population is one of the excellent study models for enhancing our understanding of the biology of heat tolerance because they are predominantly kept outdoors on pastures where they experience direct effects of weather elements (e.g., solar radiation). In this article, we focus on evidence from recent studies in Australia that leveraged large a dataset [∼40,000 animals with phenotypes and 15 million whole-genome sequence variants] to elucidate the genetic basis of thermal stress as a critical part of the strategy to breed cattle adapted to warmer environments. Genotype-by-environment interaction (i.e., G × E) due to temperature and humidity variation is increasing, meaning animals are becoming less adapted (i.e., more sensitive) to changing environments. There are opportunities to reverse this trend and accelerate adaptation to warming climate by 1) selecting robust or heat-resilient animals and 2) including resilience indicators in breeding goals. Candidate causal variants related to the nervous system and metabolic functions are relevant for heat tolerance and, therefore, key for improving this trait. This could include adding these variants in the custom SNP panels used for routine genomic evaluations or as the basis to design specific agonist or antagonist compounds for lowering core body temperature under heat stress conditions. Indeed, it was encouraging to see that adding prioritized functionally relevant variants into the 50k SNP panel (i.e., the industry panel used for genomic evaluation in Australia) increased the prediction accuracy of heat tolerance by up to 10% units. This gain in accuracy is critical because genetic improvement has a linear relationship with prediction accuracy. Overall, while this article used data mainly from Australia, this could benefit other countries that aim to develop breeding values for heat tolerance, considering that the warming climate is becoming a topical issue worldwide.

Signature selection analysis reveals candidate genes associated with production traits in Iranian sheep breeds

BACKGROUND

Sheep were among the first animals to be domesticated. They are raised all over the world and produce a major scale of animal-based protein for human consumption and play an important role in agricultural economy. Iran is one of the important locations for sheep genetic resources in the world. Here, we compared the Illumina Ovine SNP50 BeadChip data of three Iranian local breeds (Moghani, Afshari and Gezel), as a population that does not undergone artificial breeding programs as yet, and five other sheep breeds namely East Friesian white, East Friesian brown, Lacaune, DorsetHorn and Texel to detect genetic mechanisms underlying economical traits and daptation to harsh environments in sheep.

RESULTS

To identify genomic regions that have been targeted by positive selection, we used fixation index (Fst) and nucleotide diversity (Pi) statistics. Further analysis indicated candidate genes involved in different important traits such as; wool production included crimp of wool (PTPN3, NBEA and KRTAP20-2 genes), fiber diameter (PIK3R4 gene), hair follicle development (LHX2 gene), the growth and development of fiber (COL17A1 gene)), adaptation to hot arid environments (CORIN gene), adaptive in deficit water status (CPQ gene), heat stress (PLCB4, FAM107B, NBEA, PIK3C2B and USP43 genes) in sheep.

CONCLUSIONS

We detected several candidate genes related to wool production traits and adaptation to hot arid environments in sheep that can be applicable for inbreeding goals. Our findings not only include the results of previous researches, but also identify a number of novel candidate genes related to studied traits. However, more works will be essential to acknowledge phenotype- genotype relationships of the identified genes in our study.

Landscape genomics of the streamside salamander: Implications for species management in the face of environmental change

Understanding spatial patterns of genetic differentiation and local adaptation is critical in a period of rapid environmental change. Climate change and anthropogenic development have led to population declines and shifting geographic distributions in numerous species. The streamside salamander, Ambystoma barbouri, is an endemic amphibian with a small geographic range that predominantly inhabits small, ephemeral streams. As A. barbouri is listed as near‐threatened by the IUCN, we describe range‐wide patterns of genetic differentiation and adaptation to assess the species’ potential to respond to environmental change. We use outlier scans and genetic‐environment association analyses to identify genomic variation putatively underlying local adaptation across the species’ geographic range. We find evidence for adaptation with a polygenic architecture and a set of candidate SNPs that identify genes putatively contributing to local adaptation. Our results build on earlier work that suggests that some A. barbouri populations are locally adapted despite evidence for asymmetric gene flow between the range core and periphery. Taken together, the body of work describing the evolutionary genetics of range limits in A. barbouri suggests that the species may be unlikely to respond naturally to environmental challenges through a range shift or in situ adaptation. We suggest that management efforts such as assisted migration may be necessary in future.

New loci and neuronal pathways for resilience to heat stress in cattle

While understanding the genetic basis of heat tolerance is crucial in the context of global warming's effect on humans, livestock, and wildlife, the specific genetic variants and biological features that confer thermotolerance in animals are still not well characterized. We used dairy cows as a model to study heat tolerance because they are lactating, and therefore often prone to thermal stress. The data comprised almost 0.5 million milk records (milk, fat, and proteins) of 29,107 Australian Holsteins, each having around 15 million imputed sequence variants. Dairy animals often reduce their milk production when temperature and humidity rise; thus, the phenotypes used to measure an individual's heat tolerance were defined as the rate of milk production decline (slope traits) with a rising temperature-humidity index. With these slope traits, we performed a genome-wide association study (GWAS) using different approaches, including conditional analyses, to correct for the relationship between heat tolerance and level of milk production. The results revealed multiple novel loci for heat tolerance, including 61 potential functional variants at sites highly conserved across 100 vertebrate species. Moreover, it was interesting that specific candidate variants and genes are related to the neuronal system (ITPR1, ITPR2, and GRIA4) and neuroactive ligand-receptor interaction functions for heat tolerance (NPFFR2, CALCR, and GHR), providing a novel insight that can help to develop genetic and management approaches to combat heat stress.

Differential DNA methylation analysis reveals key genes in Chinese Qingyu and Landrace pigs

The Chinese Qingyu pig is a typical domestic fatty pig breed and an invaluable indigenous genetic resource in China. Compared with the Landrace pig, the Qingyu pig has unique meat characteristics, including muscle development, intramuscular fat, and other meat quality traits. At present, few studies have explored epigenetic differences due to DNA methylation between the Qingyu pig and the Landrace pig. In this study, 30 Qingyu pigs and 31 Landrace pigs were subjected to reduced representation bisulfite sequencing (RRBS). Genome-wide differential DNA methylation analysis was conducted. Six genomic regions, including regions on Sus scrofa chromosome (SSC) 1: 266.09-274.23 Mb, SSC5: 0.88-10.68 Mb, SSC8: 41.23-48.51 Mb, SSC12: 45.43-54.38 Mb, SSC13: 202.15-207.95 Mb, and SSC14: 126.43-139.85 Mb, were regarded as key regions that may be associated with phenotypic differences between the Qingyu pig and the Landrace pig. Furthermore, according to further analysis, five differentially methylated genes (ADCY1, FUBP3, GRIN2B, KIT, and PIK3R6) were identified as key candidate genes that might be associated with meat characteristics. Our findings provide new insights into the differences in DNA methylation between the Qingyu pig and the Landrace pig. These results enrich the epigenetic research of the Chinese Qingyu pig.

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

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

Liver Transcriptome Changes of Hyla Rabbit in Response to Chronic Heat Stress

Simple Summary

It has been widely acknowledged in farm animals that environmental heat stress would have comprehensive influences on many kinds of physiological aspects, including the metabolic characteristics, production performances, welfare concerns, etc. The rabbit is a small herbivore and needs to regulate the body temperature in a fine mechanism. Little is known, however, about the genes and pathways that are involved in the regulatory responses under chronic heat stress conditions. In the present study, we investigated the liver transcriptome changes in response to chronic heat stress for Hyla rabbit, that is a commercial meat breed recently introduced into China. We successfully revealed the differentially expressed genes that were significantly enriched in heat stress related biological processes. The results would help us for better understanding the molecular mechanisms underlying physiological responses against heat stress in rabbits.

Abstract

Rabbit is an economically important farm animal in China and also is a widely used animal model in biological researches. Rabbits are very sensitive to the environmental conditions, therefore we investigated the liver transcriptome changes in response to chronic heat stress in the present study. Six Hyla rabbits were randomly divided into two groups: chronic heat stress (HS) and controls without heat stress (CN). Six RNA-Seq libraries totally yielded 380 million clean reads after the quality filtering. Approximately 85.07% of reads were mapped to the reference genome. After assembling transcripts and quantifying gene expression levels, we detected 51 differentially expressed genes (DEGs) between HS and CN groups with thresholds of the adjusted p-value < 0.05 and |log2(FoldChange)| > 1. Among them, 33 and 18 genes were upregulated and downregulated, respectively. Gene ontology analyses further revealed that these DEGs were mainly associated with metabolism of lipids, thyroid hormone metabolic process, and cellular modified amino acid catabolic process. The upregulated ACACB, ACLY, LSS, and CYP7A1 genes were found to be inter-related through biological processes of thioester biosynthetic process, acyl-CoA biosynthetic process, acetyl-CoA metabolic process, and others. Six DEGs were further validated by quantitative real-time PCR analysis. The results revealed the candidate genes and biological processes that will potentially be considered as important regulatory factors involved in the heat stress response in rabbits.

Transcriptome analysis reveals corresponding genes and key pathways involved in heat stress in Hu sheep

Heat stress (HS) seriously affects animal performance. In view of global warming, it is essential to understand the regulatory mechanisms by which animals adapt to heat stress. In this study, our aim was to explore the genes and pathways involved in heat stress in sheep. To this end, we used transcriptome analysis to understand the molecular responses to heat stress and thereby identify means to protect sheep from heat shock. To obtain an overview of the effects of heat stress on sheep, we used the hypothalamus for transcriptome sequencing and identified differentially expressed genes (DEGs; false discovery rate (FDR) < 0.01; fold change > 2) during heat stress. A total of 1423 DEGs (1122 upregulated and 301 downregulated) were identified and classified into Gene Ontology (GO) categories and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Heat stress triggered dramatic and complex alterations in gene expression in the hypothalamus. We hypothesized that heat stress induced apoptosis and dysfunction in cells and vital organs and affected growth, development, reproduction, and circadian entrainment via the calcium signaling pathway, which influences ribosome assembly and function. Real-time PCR was used to evaluate the expression of the genes regulating important biological functions or whose expression profiles were significantly changed after acute heat stress (FDR < 0.01; fold change > 4), and the results showed that the expression patterns of these genes were consistent with the results of transcriptome sequencing, indicating that the credibility of the sequencing results. Our data indicated that heat stress induced calcium dyshomeostasis, blocked biogenesis, caused ROS accumulation, impaired the antioxidant system and innate defense, and induced apoptosis through the P53 signaling pathway activated by PEG3, decreased growth and development, and enhanced organ damage. These data is very important and helpful to elucidate the molecular mechanism of heat stress and finally to find ways to deal with heat stress damage in sheep.

A Discovery of Relevant Hepatoprotective Effects and Underlying Mechanisms of Dietary Clostridium butyricum Against Corticosterone-Induced Liver Injury in Pekin Ducks

Clostridium butyricum (C. butyricum) can attenuate oxidative stress, inflammation, and hepatic fatty deposition in poultry, however, the underlying mechanisms for this in Pekin ducks remain unclear. This study evaluated these hepatoprotective effects and the underlying mechanisms in a corticosterone (CORT)-induced liver injury model in Pekin ducks fed a C. butyricum intervention diet. A total of 500 Pekin ducks were randomly divided into five groups: one group (CON group) was only provided with a basal diet, three groups were provided a basal diet with 200 mg/kg (LCB group), 400 mg/kg (MCB group), or 600 mg/kg (HCB group) C. butyricum, respectively, and one group was provided a basal diet with 150 mg/kg aureomycin (ANT group) for 42 d. At 37 days-old, all ducks received daily intraperitoneal injections of CORT for five days to establish a liver injury model. C. butyricum intervention alleviated liver injury by decreasing the liver organ indices, hepatic steatosis and hepatocyte necrosis, and improving liver function, antioxidant capacity, and inflammatory factors. Hepatic RNA-seq revealed 365 differentially expressed genes (DEGs) between the MCB and CON groups, with 229 up- and 136 down-regulated DEGs in the MCB group. Between the MCB and ANT groups, 407 DEGs were identified, including 299 up- and 108 down-regulated genes in MCB group. Some DEGs in the MCB group related to oxidative stress and inflammatory responses such as Sod3, Tlr2a/b, and Il10, which were up-regulated, while Apoa1, Cyp7a1, Acsl1/5, Fasn, Ppar-γ, and Scd, which are involved in lipid metabolism, were down-regulated, indicating that these genes were responsive to dietary C. butyricum for the alleviation of corticosterone-induced hepatic injury. Toll-like receptor signaling, PI3K-Akt signaling pathway, cytokine-cytokine receptor interaction, peroxisome proliferator-activated receptor (PPAR) signaling pathway, adipocytokine and glycerophospholipid metabolism signaling pathway were significantly enriched in the MCB group. These findings indicate that C. butyricum intervention can protect Pekin ducks from corticosterone-induced liver injury by the modulation of immunoregulatory- and lipid metabolism-related genes and pathways.

Transcriptomic profiles of muscle, heart, and spleen in reaction to circadian heat stress in Ethiopian highland and lowland male chicken

Temperature stress impacts both welfare and productivity of livestock. Global warming is expected to increase the impact, especially in tropical areas. We investigated the biological mechanisms regulated by temperature stress due to the circadian temperature cycle in temperature adapted and non-adapted chicken under tropical conditions. We studied transcriptome profiles of heart, breast muscle, and spleen tissues of Ethiopian lowland chicken adapted to high circadian temperatures and non-adapted Ethiopian highland chicken under lowland conditions at three points during the day: morning, noon, and evening. Functional annotations and network analyses of genes differentially expressed among the time points of the day indicate major differences in the reactions of the tissues to increasing and decreasing temperatures, and also the two chickens lines differ. However, epigenetic changes of chromatin methylation and histone (de)acetylation seemed to be central regulatory mechanisms in all tissues in both chicken lines. Finally, all tissues showed differentially expressed genes between morning and evening times indicating biological mechanisms that need to change during the night to reach morning levels again the next day.