- 1Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
- 2Key Laboratory of Animal Genetics, Breeding and Reproduction, MARA – National Engineering Laboratory for Animal Breeding – College of Animal Science and Technology, China Agricultural University, Beijing, China
- 3Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, China
- 4Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada
- 5College of Life Sciences and Bioengineering, School of Science, Beijing Jiaotong University, Beijing, China
Understanding the biological mechanisms of climatic adaptation is of paramount importance for the optimization of breeding programs and conservation of genetic resources. The aim of this study was to investigate genetic diversity and unravel genomic regions potentially under selection for heat and/or cold tolerance in thirty-two worldwide cattle breeds, with a focus on Chinese local cattle breeds adapted to divergent climatic conditions, Datong yak (Bos grunniens; YAK), and Bali (Bos javanicus) based on dense SNP data. In general, moderate genetic diversity levels were observed in most cattle populations. The proportion of polymorphic SNP ranged from 0.197 (YAK) to 0.992 (Mongolian cattle). Observed and expected heterozygosity ranged from 0.023 (YAK) to 0.366 (Sanhe cattle; SH), and from 0.021 (YAK) to 0.358 (SH), respectively. The overall average inbreeding (±SD) was: 0.118 ± 0.028, 0.228 ± 0.059, 0.194 ± 0.041, and 0.021 ± 0.004 based on the observed versus expected number of homozygous genotypes, excess of homozygosity, correlation between uniting gametes, and runs of homozygosity (ROH), respectively. Signatures of selection based on multiple scenarios and methods (FST, HapFLK, and ROH) revealed important genomic regions and candidate genes. The candidate genes identified are related to various biological processes and pathways such as heat-shock proteins, oxygen transport, anatomical traits, mitochondrial DNA maintenance, metabolic activity, feed intake, carcass conformation, fertility, and reproduction. This highlights the large number of biological processes involved in thermal tolerance and thus, the polygenic nature of climatic resilience. A comprehensive description of genetic diversity measures in Chinese cattle and YAK was carried out and compared to 24 worldwide cattle breeds to avoid potential biases. Numerous genomic regions under positive selection were detected using three signature of selection methods and candidate genes potentially under positive selection were identified. Enriched function analyses pinpointed important biological pathways, molecular function and cellular components, which contribute to a better understanding of the biological mechanisms underlying thermal tolerance in cattle. Based on the large number of genomic regions identified, thermal tolerance has a complex polygenic inheritance nature, which was expected considering the various mechanisms involved in thermal stress response.
Introduction
Climate change is a major concern around the world as it negatively impacts the welfare and sustainability of livestock production and, consequently, food security in a rapidly growing population (Polsky and Keyserlingk, 2017; Thamo et al., 2017; Henry et al., 2018; Bernabucci, 2019). Understanding the adaptation of different cattle (Bos taurus) breeds to local environments is of paramount importance for the design and optimization of breeding programs and conservation of genetic resources (Hoffmann, 2010; Mwai et al., 2015; Naskar et al., 2015). Domestication, followed by spatial dispersion, has resulted in the formation of multiple cattle breeds that are raised in a large range of environmental conditions and production systems (Ajmone-Marsan et al., 2010; Yurchenko et al., 2018). Consequently, long-term (natural and artificial) selection pressure has operated in the genomic regions that control traits for adaptive fitness and animal resilience (Kristensen et al., 2015; Berry, 2017).
Commercial cattle breeds are mainly used in intensive farming systems, with a limited number of traditional cattle breeds raised throughout the world (Kukučková et al., 2017; Mastrangelo et al., 2018; Upadhyay et al., 2019). In contrast, traditional cattle breeds have a long history of adaptation to their local environments (Mwai et al., 2015; Zhang et al., 2018). Whereas commercial cattle far exceed the traditional (or native) cattle breeds in terms of milk or meat production, the latter breeds have great cultural value, are usually adapted to local environmental and climatic conditions, and are sometimes genetically superior for specific production or functional traits (e.g., fertility, disease resistance, workability) compared to commercial cattle breeds (Medugorac et al., 2009; Mastrangelo et al., 2018; Upadhyay et al., 2019).
Cattle were domesticated in Asia and there is a wide variety of breeds adapted to different environmental conditions, which can be clustered into groups based on their geographical origin (Lei et al., 2006). For instance, Chinese cattle breeds are mainly composed by Bos taurus taurus (especially in Northern China) and Bos taurus indicus (mainly in the South of China) or hybrids between both subspecies (Li et al., 2013; Xu et al., 2019b). In addition to the Bos taurus species, domestic yak (Bos grunniens) and Bali cattle (BLI; Bos javanicus) have a substantial economic and social role in their respective habitat. In the mountainous regions of Asia, Datong yak (YAK) provides milk, meat, fiber and hide, and transportation (Wiener et al., 2003; Wiener, 2013). YAKs have the ability to survive in extremely harsh environments at high altitudes, overcoming hypoxia, severely cold winters, and cool moist summers, with limited feeding sources (Wiener, 2013; Chen et al., 2018). However, little is known about the specific adaptive mechanisms of YAKs. BLI cattle are found in Indonesia and Malaysia, and carries unique characteristics such as ability to survive dry seasons and higher heat tolerance, high conception rate, and greater resistance to diseases when compared to buffalo and cattle breeds (Noor et al., 2001; Syed-Shabthar et al., 2013). Other breeds around the world have also adapted to survive in extreme conditions, such as African cattle breeds that have adapted to high temperatures and food shortages (Mwai et al., 2015), and others that have adapted to extreme cold temperatures (Weldenegodguad et al., 2019).
Understanding how selection and evolution act on livestock genomes can benefit the optimization of breeding programs in order to improve economically important traits (Hiemstra et al., 2010; Mastrangelo et al., 2018), and especially animal resilience. For instance, a genomic evaluation for improved heat tolerance has been recently developed for Australian dairy cattle (Nguyen et al., 2016). There is an increased interest in identifying and exploring signatures of selection (also known as selective sweeps) and the genomic diversity resulting from adaptation to environment and artificial selective pressure (Igoshin et al., 2019; Osei-Amponsah et al., 2019; Upadhyay et al., 2019). Different approaches can be used to scan the genome for regions of homozygosity, as well as to estimate differences in the frequency of alleles or haplotypes between divergent populations or even generations within a population (Boschiero et al., 2018; Ceballos et al., 2018).
Several statistical methods have been proposed to identify signatures of selection, such as the fixation index (FST) statistics (Wright, 1922), runs of homozygosity (ROH; Ku et al., 2010), and hapFLK (Fariello et al., 2013). ROH are genomic regions containing contiguous homozygous genotypes (Ceballos et al., 2018; Dixit et al., 2020). The FST is based on the measure of population differentiation due to locus-specific allele frequencies between populations (Wright, 1922; Maiorano et al., 2018) and is an alternative to analyze genetic differences across individuals (populations, breeds, or groups of animals) that may have been caused by opposing selection pressure (Beaumont, 2005; Willing et al., 2012). High FST values indicate local positive adaptation while low FST values suggest negative or neutral selection (Kullo and Ding, 2007). The hapFLK approach is based on haplotype differentiation between populations and it incorporates hierarchical structure of populations extended to account for haplotype structure in the samples (Fariello et al., 2013; Yurchenko et al., 2018). These methods have been applied in various livestock species [e.g., cattle (Randhawa et al., 2016; Cheruiyot et al., 2018), sheep (Moradi et al., 2012; Mcrae et al., 2014), goats (Brito et al., 2017; Onzima et al., 2018), pigs (Ai et al., 2013, 2014), chickens (Stainton et al., 2014; Walugembe et al., 2019)] and revealed multiple genomic regions associated with environmental adaptation (as reviewed by Randhawa et al., 2016). Therefore, the main objectives of this study were to investigate genetic diversity levels based on multiple metrics in cattle, with a focus on Chinese local breeds, and unravel genomic regions and candidate genes potentially under selection for thermal tolerance in cattle breeds adapted to divergent climatic conditions worldwide, domestic yak (Datong breed; Bos grunniens), and Bali (Bos javanicus) populations. To our best knowledge, no signatures of selection have been reported in the literature for some important Chinese cattle breeds [Yunnan Humped (DH), Dengchuan (DC), Wannan (WN), Wenling Humped (WL), Mongolian (MG), Sanhe (SH), Hazake (KZK), Xinjiang Brown (XIN)], which are unique genetic resources adapted to divergent climatic conditions. Furthermore, adding outgroups (closely related species) to the analyses can contribute to a better understanding of evolutive processes related to climatic adaptation that shaped cattle speciation, which has not been well investigated due to previous lack of genomic datasets.
Materials and Methods
No Animal Care Committee approval was necessary for the purposes of this study, as all information required was obtained from pre-existing databases.
Datasets
For the different analysis incorporated in this study different datasets were considered (specification of the populations used for the calculation of each diversity metric are described latter on). In summary, a total of 986 animals from 32 different cattle (Bos taurus) breeds, Datong yak (Bos grunniens) and Bali (Bos javanicus) were included in this study (Supplementary Table 1). Three sources of genotypes were used: (1) a set of 17 populations [DH, DC, WL, WN, Arabic Zebu (ZAR), BLI, Bangladesh zebu (BAD), Boran (BAN), Canchin (CAN), Creole from Guadalupe (CGU), East African Shorthorn zebu (EAZ), Gir (GIR), N’Dama (NDA), Senepol (SEN), Zebu Bororo (ZBO), Zebu from Madagascar (ZMA), and Zebu Fulani (ZFU)], considered as heat tolerant; (2) a set of 13 populations [SH, MG, KZK, XIN, Angus (ANG), Hereford (HFD), Holsteins (HOL), Jersey (JER), Brown Swiss (BSW), Kalmykian (KAL), Negra Andaluza (NAN), Santa Gertrudis (SGT), and Shorthorn (SHO)] considered as cold tolerant; and, (3) YAK, Busha (BUS), Yanbian (YAN), and Yakutian (YAKU) that are considered as resilient to cold and harsh environmental conditions. All animals were genotyped with the Illumina BovineSNP50K BeadChip (50K; Illumina, San Diego, CA, United States), except XIN and YAK animals that were genotyped using the GeneSeek GGP-HDV3 SNP panel (139K; GeneSeek, Lincoln, NE, United States). The 50K and 139K single nucleotide polymorphisms (SNP) panels contained 53,714 and 139,376 SNP, respectively. From the populations used in this study, the genomic information of 25 (see Supplementary Table 1) were downloaded from the Web-Interfaced next generation Database for genetic Diversity Exploration (Sempéré et al., 2015). A limitation of the other two Bos genus (YAK and BLI) datasets are the genomic datasets, in which genotyping platforms developed for cattle (a closely related Bovidae species) were used. This might have resulted in ascertainment bias (Lachance and Tishkoff, 2013; Malomane et al., 2018; Geibel et al., 2019) in the results of YAK and BLI, as discussed later. Genotypic quality control (QC) was implemented for all breeds together. SNP with minor allele frequency (MAF) lower than 0.01, animal and SNP call rate lower than 90%, and SNP located on the sex chromosomes or with unknown position in the genome were excluded from the analyses. Genotypic QC was performed using the software PLINK 1.9 (Purcell et al., 2007). The number of common overlapping SNP between SNP panels that remained after QC was 35,505. Several signatures of selection studies have successfully used similar SNP panel densities (e.g., Flori et al., 2012; Rothammer et al., 2013; Randhawa et al., 2014, 2016; Makina et al., 2015).
Genetic Diversity Metrics
Heterozygosity
Observed heterozygosity (HO) per animal and within population was calculated based on the number of heterozygotes divided by the total number of genotypes. HO values were compared to the expected heterozygosity under Hardy Weinberg Equilibrium (HE). Both metrics were calculated using the – hardy flag (Purcell et al., 2007; Chang et al., 2015).
Proportion of Polymorphic SNP (PSNP)
The PSNP was calculated as the proportion of SNP with MAF greater than 1% within each breed, indicating the proportion of SNPs segregating in each population. Both PSNP and MAF calculations were performed before the genotyping QC.
Average Pairwise Genetic Distance (D)
The average proportion of alleles shared between two individuals was calculated as , where IBS1 and IBS2 are the number of loci that share 1 or 2 identical-by-state (IBS) alleles, respectively, and m is the total number of loci. The genetic distance between all pair-wise combinations of individuals was calculated as: D = 1 − DST.
Runs of Homozygosity
Runs of homozygosity were identified using the – homozyg flag (Purcell et al., 2007; Chang et al., 2015). The following criteria were considered to define a ROH segment: (i) scanning window of 50 SNP across the genome; (ii) minimum density of 50 kb/SNP; (iii) minimum ROH length of 500 kb (Xu et al., 2019a), and, (iv) maximum gap between consecutive homozygous SNP of 100 kb. To investigate genomic regions under possible selection, non-overlapping ROH segments between the heat and cold adaptation groups, and that were shared at least between 50% of the animals in the same breed were selected and compared with the other signature of selection metrics.
Inbreeding Coefficients
Four different measurements of genomic inbreeding were calculated for all cattle populations and YAK: (1) based on the observed versus expected number of homozygous genotypes (FEH), (2) genomic inbreeding based on excess of homozygosity (FHOM), (3) correlation between uniting gametes (FU), and (4) ROH-based inbreeding (FROH). FEH, FHOM, and FU were calculated using all genotyped animals and SNP that remained after the QC performed individually for each population. FROH was calculated as the genome length covered by ROH divided by the total genome length covered by SNPs across all 29 autosomes. Different inbreeding metrics are expected to differ as they capture complementary populational diversity characteristics (Meuwissen et al., 2020).
Linkage Disequilibrium (LD)
The extent of LD (measured as r2) was calculated for each breed group using the –r2 flag (Purcell et al., 2007), as following (Sved, 1971):
where f(AB), f(A), f(a), f(B), and f(b) are observed frequencies of AB, A, a, B, and b, respectively, and D2 = [f(AB) − f(A)∗f(B)]2. Within-population quality control was used to calculate LD (or r2) for each breed. The LD decay was then analyzed for the following different genetic distance classes between SNP pairs (0, <10, <20, <40, <60, <100, <200, <500, <1,000 kb).
Consistency of Gametic Phase
Consistency of gametic phase was determined by calculating the square root of the LD values (r2) and assigning the appropriate sign based on the calculated disequilibrium (D) metric, as used in the calculation of LD. The D values were calculated using the –dprime-signed option (Purcell et al., 2007; Chang et al., 2015). Thereafter, the consistency of gametic phase was assumed as the Pearson correlation coefficient between each two breed-group pair, using the signed-squared-root values. The breakdown in the consistency of gametic phase across distances was determined based on the same bins described above. Only SNP in common (after within-population QC) among the cattle populations were used to calculate consistency of gametic phase.
Ancestral and Recent Effective Population Size (Ne)
The LD combined with distance between markers can be used to estimate the approximate Ne at a given time point in the past. To estimate Ne through time, the formula as described by Sved (1971) was used:
where, Ne is the effective population size, and r2 is the average r2 value at a given distance c. Each genetic distance (c) corresponds to a value of t generations in the past. This value was calculated as (Hayes et al., 2003). The Ne in generation 5 was assumed as the recent Ne and ancient Ne was investigated up to 500 generations in the past.
For all the above-mentioned metrics 35,505 SNPs from 312 individuals from eight Chinese cattle breeds (DH, DC, WL, WN, SH, MG, KZK, and XIN) and YAK population were used. These metrics were only calculated to these populations as all the other datasets were extracted from public databases and similar studies have been published elsewhere (e.g., Chen et al., 2016; Boushaba et al., 2018; Edea et al., 2018).
Principal Component Analysis (PCA)
For the PCA and genomic population tree (described below) all individuals were considered (n = 986) from the 34 populations (Supplementary Table 1). In order to better visualize the genetic similarity among populations, PCA was assessed using the –pca flag (Purcell et al., 2007; Chang et al., 2015). Two scenarios were analyzed: (1) including all the 34 populations; and (2) considering only the 32 cattle populations (i.e., excluding YAK and BLI).
Genomic Population Tree
The neighbor-joining algorithm was used to plot population trees using pair-wise Reynolds distance between populations (Fariello et al., 2013). The genomic population tree was created using all genome-wide SNP that passed the QC. The software PLINK 1.9 (Purcell et al., 2007; Chang et al., 2015) was used for calculating observed and expected heterozygosity, proportion of polymorphic SNP, genetic distance, ROH, genomic inbreeding coefficients, LD, gametic phase, and PCA.
Signatures of Selection Analyses
Fixation Index
To identify population-specific loci under positive selection, FST was calculated for each of the 35,505 informative SNP along the genome using different contrasting groups to estimate the differences in allelic frequencies. In brief, FST was calculated as the squared deviation of the average frequency in a specific population from the average frequency across all populations divided by the allele frequency variance (p∗q, in which p and q are the frequencies of each alternate SNP alleles). In order to identify putative genomic regions under selection, the analyses were performed by grouping the cattle breeds in divergent clusters defined as heat, mild cold, and extreme cold tolerant populations. Preliminary analyses with balanced number of individuals per group were performed, but similar results were observed (results not shown). Therefore, the results with the complete datasets are presented here. Genomic regions under a mean plus four SD threshold were considered as candidate regions. To identify putative genome regions under selection, the analyses were performed under four scenarios of contrasting models:
(1) Groups based on climate adaptation (FST1): The groups (n = 3) included all populations and were created based on climate adaptation: heat adaptation (DH, DC, WL, WN, ZAR, BLI, BAD, BAN, CAN, CGU, EAZ, GIR, NDA, SEN, ZBO, ZMA, and ZFU), mild cold/thermoneutral adaptation (SH, MG, KZK, XIN, ANG, HFD, HOL, JER, BSW, KAL, NAN, SGT and SHO) and harsh cold adaptation (YAK, BUS, YAN, and YAKU).
(2) Groups based on climate adaptation and breeding purpose excluding YAK and BLI populations (FST2): The groups (n = 3) included all populations mentioned for FST1, but YAK and BLI were excluded.
(3) Most climatic-divergent individual breeds (FST3): The three most climatic-divergent populations (one for each group described above) based on FST values were compared against each other. The populations were DH (heat group), HOL (mild cold/thermoneutral group), and YAKU (harsh cold group).
HapFLK
The hapFLK method was used for calculation of distance matrices and the estimation of hapFLK as described by Fariello et al. (2013). Three scenarios were considered and included all populations (Supplementary Table 1). For the first scenario, YAK was considered as an outgroup; for the second BLI was considered as an outgroup; and for the third scenario there was no outgroup and only cattle breeds (32 breeds; Supplementary Table 1) were included. A Reynolds distance matrix was calculated to estimate the hierarchical population structure within each population set. The fastPHASE (Scheet and Stephens, 2006) cross-validation procedure was used to determine the number of underlying latent states. A total of 21 haplotype clusters were needed to capture the total haplotype diversity.
Candidate Genes and Enriched Function Analyses
The positional candidate genes located in genomic regions containing relevant SNP based on FST, ROH, and HapFLK analyses were mapped using the Biomart tool (Kinsella et al., 2011), embedded in the Ensembl Genes database version 991. Based on the start and end chromosomal positions, important genomic regions were further investigated to identify biological processes related to climatic adaptation and to define prospective functional genes. Complete gene functions were obtained from the National Center for Biotechnology Information Database2 based on the latest Bos taurus reference genome ARS-UCD1.2 (Rosen et al., 2020) and from the Animal QTL database (Hu et al., 2018). Biological functions and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (Kanehisa, 2000; Kanehisa et al., 2016, 2018) in which the identified candidate genes are involved were assessed using the Database for Annotation, Visualization and Integrated Discover tool (DAVID v6.8; Huang et al., 2008a, b). Four different scenarios were analyzed during the enriched function analyses considering high classification stringency taking into consideration: (1) all candidate genes found in all methods together (heat and cold adaptation; SCEN1), (2) candidate genes identified in methods based on heat adaptation (SCEN2), (3) candidate genes identified in methods based on cold adaptation (SCEN3), and (4) candidate genes identified for YAK (SCEN4).
Results
Genetic Diversity Metrics
The following genetic diversity metrics were estimated from eight Chinese cattle breeds (DH, DC, WL, WN, SH, MG, KZK, and XIN) and YAK population. As previously mentioned, these metrics were only calculated to these populations as all the other datasets were extracted from public databases and similar studies have been published elsewhere (e.g., Chen et al., 2016; Boushaba et al., 2018; Edea et al., 2018).
Heterozygosity
HO ranged from 0.023 (YAK) to 0.366 (SH). The average HO was lower than the HE for DH, DC, Wannam, and MG populations, indicating slightly lower genetic diversity levels than expected.
Proportion of Polymorphic SNP
The PSNP was high in all analyzed populations (ranged from 0.659 to 0.992), except in YAK (0.197). The distribution of SNP percentage by MAF range is shown in Figure 1. In general, the distribution of SNP percentage was variable across the different populations, especially for YAK, which showed the highest value for MAF less than 0.01 and the lowest values for the other MAF ranges.
Figure 1. Distribution of single nucleotide polymorphisms (SNP percentage) by minor allele frequency ranges in eight Chinese cattle breeds and Datong yak.
Average Pairwise Genetic Distance
D ranged from 0.016 (YAK) to 0.290 (MG) (Table 1). Higher values indicate greater within-population genetic variation (heterogeneity).
Table 1. Genetic diversity metrics estimated within each one of the eight Chinese cattle breeds and Datong yak.
Inbreeding Coefficients
Four different approaches to estimate genomic inbreeding coefficients were used. The average genomic inbreeding coefficients estimated using different approaches are shown in Table 1. The overall average genomic inbreeding (±SD) coefficients were: 0.254 ± 0.028 (FHOM), 0.250 ± 0.051 (FEH), 0.194 ± 0.041(FU), and 0.021 ± 0.004 (FROH). The average inbreeding coefficients differed among cattle breeds and YAK had the highest values for FHOM, FU, and FROH metrics. Table 2 presents the Pearson correlations among the different measures of genomic inbreeding coefficients. In general moderate to high correlations were observed among all methods (0.57–0.89), with exception of FEH and FHOM, which surprisingly had a negative correlation (−0.71).
Table 2. Pearson correlation among inbreeding coefficients in eight Chinese cattle breeds and Datong yak.
Runs of Homozygosity
The proportion of ROH segments in each length category for the nine populations are shown in Figure 2 and Supplementary Table 1. The majority of ROH segments was of short length (i.e., segments were shorter than 2,500 kb or between 2,500 and 5,000 kb), indicating ancient inbreeding. DH, DC, WN, and XIN had the highest proportion of short segments compared to the other populations. Only a small proportion of ROH segments were longer than 15,000 kb, except for YAK, which had the highest proportion of longer ROH segments, and MG cattle population. The KZK population had the highest proportion of short ROH segments and, in contrast, no ROH segment higher than 15,000 kb.
Figure 2. Proportion of runs of homozygosity segments in each length category for eight Chinese cattle breeds and Datong yak.
A descriptive summary of ROH is presented in Table 3. The average number of ROH segments for each animal within breed ranged from 2.57 ± 1.4 (KZK) to 163.5 ± 16.39 (YAK), with a maximum of 247 segments in a YAK. The average length of genome contained within ROH segments ranged from 3,505 kb ± 2,409.51 (KZK) to 14,669 kb ± 16,059.5 (YAK). The animal with the longest proportion of its genome characterized as ROH was also observed in the YAK population (123,041 kb). The average number of SNP in ROH run per population ranged from 81.28 ± 100.79 (WN) to 664.9 ± 705.51 (YAK), with minimum and maximum values of 50 and 5,398 SNP, respectively. The average SNP density (SNP per kb) was similar for all populations (∼43 SNP/kb), except for XIN and YAK which were ∼21 SNP/kb. For all populations, the proportion of homozygous sites was higher than 97%.
Table 3. Descriptive analysis of the runs of homozygosity per population in eight Chinese cattle breeds and Datong yak.
Linkage Disequilibrium
The LD decay pattern for all populations is shown in Figure 3. In general, the highest LD was observed for YAK [ranged from 0.642 (at 0.02 Mb) to 0.290 (at 1 Mb)] and DH [ranged from 0.453 (at 0.01 Mb) to 0.074 (at 1 Mb)]. The LD decay pattern for WN was similar to DC cattle. WN [ranged from 0.346 (at 0.01 Mb) to 0.069 (at 1 Mb)] and DC [ranged from 0.347 (at 0.01 Mb) to 0.045 (at 1 Mb)] had the lowest LD values.
Figure 3. Average linkage disequilibrium (LD) at given distances for eight Chinese cattle breeds and Datong yak.
Consistency of Gametic Phase
As an example, the consistency of gametic phase between DH and SH with the other Chinese cattle breeds are shown in Figure 4. Among all populations DC, WN, and Weiling Humped had the highest consistency of gametic phase with DH (Figure 4A) at all analyzed distances, ranging from 0.594 (at 10 Kb) to 0.047 (at 1,000 Kb), from 0.618 (at 10 Kb) to 0.061 (at 1,000 Kb), and from 0.628 (at 10 Kb) to 0.067 (at 1,000 Kb), respectively. On the other hand, MG and XIN populations had the lowest consistency of gametic phase with DH, ranging from 0.319 (at 10 Kb) to 0.004 (at 1,000 Kb) and from 0.298 (at 10 Kb) to 0.018 (at 1,000 Kb), respectively. Similar trends to DH have been observed for DC and WN with the other populations. It was observed that the consistency of gametic phase values between SH and the other cattle breeds (Figure 4B) have similar values after 100 Kb distance. The taurine breeds have similar consistency of gametic phase values with each other and the same was observed with the hybrid populations. The consistency of gametic phase between all population pairs are shown in Supplementary Table 2.
Figure 4. Consistency of gametic phase at given distances between (A) Yunnan Humped and seven Chinese cattle breeds, and (B) Sanhe and seven Chinese cattle breeds.
Ancestral Effective Population Size
The historical Ne ranged from 3,032 (MG breed) at 500 generations ago to 47 (KZK) at five generations ago. The Ne five generations ago (considered as recent Ne) were quite variable across populations, ranging from 147 (DC) to 47 (KZK). YAK had the smallest Ne 500 generations ago (Ne = 195) and 66 five generation ago. The Ne values from 5 to 500 generations ago are shown in Table 4.
Table 4. Ancestor effective population size based on linkage disequilibrium (r2), calculated from generation 5 to 500 for eight Chinese cattle breeds and Datong yak.
Principal Component Analysis
The genetic relationship among 32 cattle breeds, YAK and BLI populations revealed by PCA are shown in Figure 5. The first, second, and third principal components explained 43.80%, 24.68%, and 20.86% of the total variability among cattle breeds, YAK and BLI, respectively. In general, all three components show a clear differentiation of groups adapted to heat and cold climatic conditions. Populations adapted to heat were grouped together and the same occurred with those adapted to cold, and YAK and BLI were clustered as two separate group from the other populations, as expected. Breeds with similar development were clustered together such as XIN and BSW. Two clear clusters were observed for the Zebu breeds (Bos indicus; ZBO, ZFU, ZMA, and EAZ) and another group of Chinese cattle breeds adapted to heat (DH, DC, WL, and WN), which have small contributions from Bos taurus indicus from crossing over centuries ago. The PCA plots considering only the 32 cattle populations are shown in Supplementary Figure 1. The first, second, and third principal components explained 58.61%, 24.82%, and 19.94% of the total variability among cattle breeds, respectively. It was observed that the populations remained grouped into common clusters based on adaptation to heat or cold stress, despite of their sub-species. As in the first PCA scenario, there were clusters for grouping XIN and BSW, Zebu breeds and Chinese heat adapted breeds. In the two PCA scenarios, specially between first and second, and second and third PC, it is also possible to see that breeds adapted to the cold/thermoneutral temperatures are located together (e.g., HOL, JER, and NAN).
Figure 5. Principal component decomposition of the genomic relationship matrix colored by population*. Letters in the figure represent the decomposition of the: (A) first and second; (B) first and third; and (C) second and third principal components, respectively (*see Supplementary Table 1 for abbreviations).
Genomic Population Tree
The genomic population tree (constructed based on whole-genome SNPs) for 32 cattle breeds is presented in Figure 6. In summary, the populations were grouped according to their specific sub-species, i.e., hybrid (Bos taurus indicus × Bos taurus taurus), taurine (Bos taurus taurus) or indicine (Bos taurus indicus) breeds were clustered separately. Also, populations with similar development background were clustered together (e.g., XIN and BSW).
Figure 6. Genomic population tree comparing the genomic distance between 32 cattle breeds*. (*see Supplementary Table 1 for abbreviations).
Signatures of Selection Analyses
Fixation Index
High FST values indicate potential positive selection while low FST values suggest negative or neutral selection. Considering 32 cattle breeds, YAK and BLI populations there were several genomic regions potentially under selection in at least one of the populations studied. These selective sweeps were distributed on all 29 autosomal chromosomes. In scenario 1 (FST1), where we contrasted climate adaptation groups (n = 3), we observed 303 significant SNP and 3 shared SNP were found between the groups (one between mild cold/thermoneutral and heat, and two between heat and harsh cold). The highest number of SNP identified was on BTA1 (15 SNP) in the harsh cold adapted group and the lowest was on BTA1 (1 SNP) in the cold adapted group. In FST2, which contrasted climate adaptation groups (n = 3), but excluding YAK and BLI from the comparison, 858 significant SNP were identified. Thirty-one SNP were shared between harsh cold and the heat adapted group, and a total of 827 unique SNP were observed. The highest number of SNP identified was on BTA1 (40 SNP) in the harsh cold adapted group and one SNP were found to be significant in BTA 3 in the cold/thermoneutral group. For FST3, where the most climatic-divergent breeds were compared, 987 significant SNP were observed. For the mild cold/thermoneutral and heat groups there were seven overlapping SNP, 37 between heat and harsh cold groups and 11 between mild cold/thermoneutral and harsh cold groups. The highest number of SNP identified was on BTA2 (62 SNP) in DH and the lowest was on BTA28 (7 SNP) in the MG population. As an example of the FST process, Figure 7 shows scans for signatures of selection based on FST1, considering climate adaptation and including YAK and BLI genotypes. Supplementary Figures 2, 3 show the genome scans for FST2 and FST3, respectively.
Figure 7. SCENARIO 1: Whole genome scans for signatures of selection using the FST approach for 32 cattle breeds and yak, considering climate adaptation. Letters in the figure represent: (A) Heat adapted populations, (B) Mild cold/thermoneutral adapted populations, and (C) Harsh cold adapted populations.
HapFLK
The hapFLK method enables the identification of genomic regions of outlying differentiation between populations while accounting for their hierarchical structure (Fariello et al., 2013). Significant peaks were identified using the hapFLK metric for assessing haplotype differentiation between populations. Peaks with p-values < 0.005 were considered significant. These peaks were located in BTA3, BTA4, BTA10, and BTA14, for the first scenario, using YAK as an outgroup. In the second scenario, which had BLI population as outgroup, peaks were also found on BTA3, BTA4, BTA10, and BTA14, as well as on BTA15. In the third scenario, which considered only the 32 cattle populations and no outgroup, there were also significant peaks in BTA3, BTA4, BTA10, and BTA14. The highest peaks for the three scenarios were located on BTA10. The haplotype based hapFLK metric and −log (p-values) for the whole genome are shown in Figure 8 for the first scenario and Supplementary Figures 4, 5 for the second and third scenario, respectively.
Figure 8. Whole genome scans for signatures of selection using the haplotype based hapFLK metric considering 32 cattle breeds and Yak as outgroup (first scenario). Odd and even numbered chromosomes are shown in yellow and black, respectively. SNP number is given on the x-axis, and the genome-wide threshold corresponding to P < 0.001, P < 0.005, and P < 0.01 is shown as horizontal blue, green and red lines, respectively.
Runs of Homozygosity
The highest number of ROH segments in the cattle populations were found for XIN and Weiling Humped (1,887 and 1,804, respectively) and the lowest (36) was found in KZK. The highest number of ROH segments were observed in BTA6 in XIN (179 ROH segments) and BTA1 in the Weiling Humped (153 ROH segments). The lowest number of ROH segments was observed in BTA28 in the DC, SH, and MG breeds (2 ROH segments). No ROH segments were found for KZK in BTA7, BTA15, BTA16, BTA18, BTA19, BTA26, BTA27, BTA28 and BTA29. The number of ROH segments observed in YAK was substantially greater (7,358 ROH segments) in comparison to the cattle populations. The number of ROH segments per chromosome in YAK ranged from 145 (BTA25) to 487 (BTA5).
Positively Selected Regions
The genomic regions with FST values greater than the mean plus four SD for FST1, FST2, and FST3; hapFLK p-values smaller than 0.005, and regions in ROH were further investigated to identify potential genes under positive selection for thermal tolerance. A summary of the number of regions identified and the number of genes for each method and scenario is presented in Table 5. The largest number of candidate genes were found based on ROH and FST1, 8,702 and 1,206 candidate genes, respectively. On the other hand, the lowest number (30) was found for hapFLK when considering YAK as an outgroup. Some of the genes located in at least two different metrics were considered as relevant genes for further investigation. For instance, some of these genes are related to heat shock proteins (e.g., DNAJB13), oxygen transport (e.g., HBB and HBG; Neupane et al., 2017), mitochondrial DNA maintenance (e.g., PIF1; Foury and Kolodynski, 1983), metabolic activity (e.g., HSPE1; Kumar et al., 2015), feed intake (e.g., SNORA70; Pitt et al., 2018), carcass conformation (e.g., SNORA70; Pitt et al., 2018; PLCL1; Shin et al., 2015), fertility and reproduction traits (e.g., RBPMS2, OAZ2; Fernández et al., 2019; IL6ST; Song et al., 2009; IGSF8; Glazar and Evans, 2009), and olfactory receptors (e.g., OR51A7; Olivieri et al., 2016).
Table 5. Summary of the number of regions identified and the number of genes (inside brackets) under positive selection for individual (diagonal) method/scenario and between methods/scenarios (off diagonal).
Several novel candidate genes were found when combining different methods. Two candidate genes were identified in YAK [RBPMS2 (BTA10: 45.31–45.33 Mb) and PIF1 (BTA10: 45.27–45.27 Mb)] and one was present only in MG cattle [ATG16L2 (BTA15: 52.35–52.36 Mb)]. These same two candidate genes present in YAK were identified between FST2 and hapFLK. Two novel candidate genes were also identified based on hapFLK and FST3, including OAZ2 (BTA10: 45.36–45.37 Mb) and CSNK1G1 (BTA10: 45.66–45.81 Mb). Larger number of overlapping regions were identified when comparing scenarios based either on FST or hapFLK to the ROH-based approach. For instance, the gene SNORA70 (BTA11: 99.006–99.007 Mb) was identified in DC cattle based on ROH and FST1. All genomic regions and candidate genes are reported in Supplementary Tables 4, 5 for individual method and overlapping regions among methods, respectively, and in Supplementary Table 6 are the candidate genes identified between the three methods (or scenarios).
The functionally enriched categories of candidate genes are shown in Supplementary Tables 7–10. The main functionally enriched categories with the respective genes within each cluster are shown in Table 6. Considering the first scenario (SCEN1), the candidate genes associated with thermal stress were significantly enriched into 116 functional categories (P < 0.05) including hemoglobin complex (GOTERM_CC), oxygen binding (GOTERM_MF), oxygen transport (GOTERM_BP), epidermal growth factor-like domain (INTERPRO), carbohydrate biosynthetic process (GOTERM_BP), steroid hormone receptor and steroid hormone receptor activity (INTERPRO and GOTERM_MF), folate receptor and folate receptor-like (INTERPRO and GOTERM_MF), folic acid transporter activity and folic acid binding (INTERPRO and GOTERM_MF), biosynthesis of unsaturated fatty acids (KEGG_PATHWAY), sperm capacitation and fertilization (GOTERM_BP and UP_KEYWORDS), immunoglobulin-like domain and immunoglobulin subtype (INTERPRO), and immunoglobulin I-set and immunoglobulin subtype 2 (INTERPRO). In the second scenario (SCEN2), where candidate genes identified in the signatures of selection analyses in the heat adaptation groups were analyzed, 65 significantly enriched functional categories were identified including hemostasis and blood coagulation (UP_KEYWORDS), platelet activation (GOTERM_BP), fertilization (UP_KEYWORDS), sperm capacitation and single fertilization (GOTERM_BP), immunoglobulin-like domain and immunoglobulin subtype (INTERPRO), sugar/inositol transporter, general substrate transporter and sugar transporter (INTERPRO), myosin head [motor domain (INTERPRO)], myosin complex (GOTERM_CC), myosin (UP_KEYWORDS), steroid hormone receptor and nuclear hormone receptor [ligand-binding-core (INTERPRO)], steroid hormone receptor activity (GOTERM_MF), collagen triple helix repeat and collagen trimer (INTERPRO and GOTERM_CC), collagen (UP_KEYWORDS) and respiratory chain (UP_KEYWORDS).
Table 6. Main enriched function categories and candidate genes for thermal tolerance in eight local Chinese cattle and Datong yak populations.
Considering the third scenario (SCEN3), where candidate genes identified in the signatures of selection analyses in the cold adaptation groups were analyzed, 36 significantly enriched functional categories were identified including myosin head [motor domain (INTERPRO)], myosin complex (GOTERM_CC), innate immunity and innate immune response (UP_KEYWORDS and GOTERM_BP), immunity (UP_KEYWORDS), cell–cell adhesion and cell–cell adherens junction (GOTERM_BP and GOTERM_CC), and immunoglobulin I-set and immunoglobulin V-set (INTERPRO). Lastly, we analyzed candidate genes identified in the signatures of selection analyses for YAK (SCEN4) and 35 significantly enriched functional categories were identified including salivary secretion and gastric acid secretion (KEGG_PATHWAY), oxygen binding (GOTERM_MF), hemoglobin complex (GOTERM_CC), oxygen transport (UP_KEYWORDS), oxygen transport activity (GOTERM_MF), immunoglobulin subtype 2 and immunoglobulin-like domain (INTERPRO), myosin head [motor domain (INTERPRO)], myosin complex (GOTERM_CC), secondary metabolites biosynthesis, transport, and catabolism (COG_ONTOLOGY), steroid hormone receptor and steroid hormone receptor activity (INTERPRO and GOTERM_MF), and prolactin signaling pathway (KEGG_PATHWAY).
Discussion
Genetic Diversity Metrics
Population Characterization and PSNP
In general, the distribution of SNP percentage was approximately constant by MAF bins (Figure 1), and the PSNP was high in the majority of the eight Chinese cattle breeds (DH, DC, WL, WN, SH, MG, KZK, and XIN) and YAK population studied in this section (Table 1). These results indicate the appropriateness of the SNP panels used for genomic studies in these populations as well as substantial genetic diversity, which could be due to the fact that these local breeds have not undergone intensive artificial selection, as experienced in commercial and worldwide spread populations. The higher PSNP value found for taurine compared to hybrid breeds is in agreement with the literature (Lin et al., 2010; Edea et al., 2013). For instance, Edea et al. (2013) studying Korean Hanwoo cattle (Bos taurus taurus) also found high PSNP (0.952). YAK was an exception, both for the distribution of SNP percentage and for the proportion of polymorphic SNP. These differences between YAK and the cattle breeds may be due to the fact that YAK was not included in the development of the SNP chip panel and the fact that YAK and cattle are different species (Bos gruniens and Bos taurus), even though closely related in the evolutionary tree. To minimize and avoid this potential issue, a SNP chip panel could be developed specifically for YAKs. Alternatively, genotyping-by-sequencing (GBS) could be employed (Qiu et al., 2015), which can be performed at a low cost per sample and tend to avoid ascertainment bias. The candidate genes obtained based on the scenarios including YAKs are biologically reasonable based on its known adaptation processes. Therefore, the results obtained are useful and worth reporting in this study.
Heterozygosity and Genetic Distance
Heterozygosity measures the level of genetic variation within a population. In this context, the levels of HO and HE were moderate to high (>0.166; Table 1), except for YAK where the values were 0.023 and 0.021 for HO and HE, respectively. Similar levels of HO and HE were also reported by Zhang et al. (2018), studying XIN (HO = 0.313 and HE = 0.288) and MG (HO = 0.331 and HE = 0.332) cattle breeds. Similar heterozygosity values (HO = 0.340 and HE = 0.341) compared to those found for taurine breeds in this study have been reported in Australian HOL cattle (Zenger et al., 2007) as well as in African cattle breeds (HO = 0.28–0.33 and HE = 0.29–0.34; Boushaba et al., 2018). HO was lower than HE in the DH, DC, WN, and MG breeds indicating increased inbreeding and reduced genetic diversity levels (Sharma et al., 2016). The differences in heterozygosity levels observed among populations, especially for YAK, can be partially explained by the SNP array design, which did not include all the cattle populations evaluated in this study and therefore, ascertainment bias might have influenced the estimates. It is worth noting that cattle SNP chip panels have been used for genotyping animals from other Bovidae species (e.g., bison, water buffalo, YAK), which generated very relevant results (e.g., Michelizzi et al., 2011; Wu et al., 2012; Oleński et al., 2015).
The D value ranged from 0.149 to 0.290 in the eight Chinese cattle populations (Table 1), suggesting variation on the levels of genetic diversity within each population. For YAK, this value was lower, indicating lower genetic distance between individuals or it may be due to the fact of the SNP chip panel used. The results for genetic distance corroborate with the PCA results, where low genetic variation between individuals was found within each cattle breed and practically no variation was found between YAK individuals. The PCA showed that some of the populations were more clearly separated while others were clustered more closely. The animals, in this study, classified as adapted to heat and belonging to the same subspecies were grouped together. In all PC, genetic uniformity was observed for hybrid (taurine × indicine) animals adapted to heat conditions. For taurine animals, there was greater heterogeneity between and within populations. In all PC, greater heterogeneity was observed in the XIN population.
Genomic Inbreeding Coefficients
Inbreeding depression has been reported to affect reproductive and productive traits, longevity, health, and the ability of the individuals to cope with environmental challenges (Parland et al., 2009; Leroy, 2014; Pereira et al., 2016; Reverter et al., 2017; Doekes et al., 2019). Therefore, monitoring and controlling inbreeding is important to avoid inbreeding depression. The levels of inbreeding varied among populations and differed substantially among methods. Overall, the levels of genomic inbreeding were lower based on FEH and FROH, and higher based on FHOM. Inbreeding based on ROH is highly dependent on the ROH length, which can change according to the population analyzed (Rodríguez-Ramilo et al., 2019), sample size, and parameters used for the detection of ROH. MG, XIN, and YAK had the highest level of FROH, a fact that can be attributed to the greater number of longer ROH regions.
Moderate correlation was observed between FROH and FEH, which can be explained by the fact that both FROH and FEH directly reflect homozygosity levels in the genome (Zhang et al., 2015). Zhang et al. (2015) reported correlation between FROH and FHOM (0.61) and between FROH and FU (0.15), which are similar to the findings in this study. A similar correlation to our study (−0.07 versus −0.10) were found between FHOM and FU in Danish Red cattle and a negative correlation, but higher than found here (−0.07 versus −0.44), was observed between FHOM and FU in Danish HOL cattle bulls (Zhang et al., 2015). As discussed by Zhang et al. (2015), FEH tends to be less accurate for populations with a high level of heterozygosity, as observed for SH, MG, KZK, and XIN. The authors indicated that a large sample size would be required to obtain a better estimate of genomic inbreeding in these populations. Solé et al. (2017) studying genomic inbreeding levels in Belgian Blue cattle found slightly higher (0.97) correlation between FHOM and FROH than that found in this study (0.73). High correlation between FHOM and FROH (0.85) has also been reported in other Chinese cattle breeds (Xu et al., 2019a). The differences in inbreeding coefficients can be related to sampling, sample size, accuracy of observed allele frequencies, statistical method used, and the SNP array (or genotyping platform) used.
The degree of relatedness between individuals within a population can be associated with ROH levels. Thus, population history events can also be inferred based on ROH analyses. In this context, short ROH are most likely correlated to ancestral inheritance or potential ancient bottleneck, whereas long ROH segments are more likely associated with relatively recent inbreeding (Purfield et al., 2012). The number of longer ROH in the MG breed may be attributed to more recent inbreeding. The same fact was observed for YAK, with the highest proportion of long ROH, but in this case it may also be due to the SNP chip used, as previously discussed.
LD and Consistency of Gametic Phase
The high average LD estimation for YAK are similar to those found in the literature. Wang et al. (2014) studying the genome variation between wild and domestic YAK using GBS found higher values of LD and lower LD decay for domesticated YAK than for wild YAK. Higher LD estimates have also been reported for Jinchuan yak (Wang et al., 2020). The higher LD values and lower LD decay pattern observed for YAK may be due to a potential bottleneck event during the domestication process (Wang et al., 2014).
The consistency of gametic phase evaluates the association between SNPs and QTL alleles across breeds (Brito et al., 2017). Thus, the linkage phase will not be consistent across populations over long distances in the genome if the genetic distance between populations is large. The consistency of gametic phase between breeds indicates whether or not different breeds could potentially be pooled into one common training population to better estimate SNP effects for genomic prediction of breeding values. The highest consistency of gametic phase was found between the hybrid cattle breeds, suggesting a greater level of relatedness between these populations. In dairy cattle, de Roos et al. (2009) reported that the benefits of combining populations in a training set were higher when the populations had diverged for only a few generations ago, when the marker density was high, and when the trait heritability was low.
In this study, the extent of LD at different distances between markers was computed and the ancestral effective population sizes was estimated. The use of LD pattern can be explored to comprehensively understand the population evolutionary history (Xu et al., 2019b). Knowledge of the effective size of a population can facilitate the design of selection schemes in animal breeding and the management of populations for endangered breeds or genetic lines (Mokry et al., 2014). Evaluating Ne in Chinese cattle populations, Xu et al. (2019b) reported recent Ne for South Chinese cattle populations (Ne = 129) similar to the ones found in this study for DC, DH, Weiling Humped, and WN cattle. There was an increase in Ne for YAK in the most recent generations (∼50 generations ago), which is probably due to the introduction of outbred animals or a high level of crossbreeding with other populations. To ensure a long-term viable genetic diversity, a threshold of Ne equals to 100 has been suggested (Meuwissen, 2009). In our study, only three populations are above the threshold (Table 4), indicating that care should be taken in this regard to ensure that the Ne and, consequently, reasonable genetic diversity levels are maintained.
PCA and Genomic Tree
The PCA showed a discrimination between heat adapted, cold adapted, YAK and BLI populations. As expected, the populations adapted to heat conditions clustered together as well as the cold ones, and YAK and BLI clustered separately (Figure 3). However, the first and third principal components showed an overlap among individuals from all different populations. These findings suggest that there is moderate genetic variability between the three groups and genetic similarity between populations clustered in the same group (i.e., heat and cold resilient and YAK). Especially on the first and third, and third and fourth components, it is possible to observe that XIN overlap with the BSW population which is also adapted to cold, and one of the breeds that XIN derived from Figures 5B,C. Another interesting cluster is between breeds from China (DH, DC, WN, WL) adapted to heat environments and are formed by crosses between taurines and indicine breeds. The Bos taurus indicus breeds formed another cluster, clearly visualized in both PCA scenarios. When considering just cattle population in the analysis (Supplementary Figure 1) the cluster from the previous scenario remains.
The genomic population tree constructed based on the genomic distance estimated between the cattle populations showed that there is differentiation between heat and cold adapted cattle populations. In this context, the breeds adapted to heat appear to be more related among themselves as well as those adapted to cold conditions. Our findings will be useful for policy makers, geneticists, and cattle breeders when making decisions toward the improvement of genetic resources, while conserving local breeds around the world.
Signatures of Selection and Functional Analyses
Identifying signatures of selection contribute to better elucidate the mechanisms of selection to adaptation and climatic resilience and to identify candidate genes of interest to breeding programs and molecular studies. Using different methods, various genomic regions that are potentially under selection in at least one population or breed group were identified. Due to the fact that climate adaptation is a complex trait, influenced by several genes and metabolic pathways, the detection of only few loci under positive selection can be challenging (Kemper et al., 2014). The use of multiple methods reduces the chance of false positive results. In this study, there were many overlapping regions identified between at least two approaches (Supplementary Table 5), which provides stronger evidence of selection for specific traits (Hohenlohe et al., 2010; Oleksyk et al., 2010; Abied et al., 2020). It is worth mentioning that the use of different methods can result in the identification of different regions (Fariello et al., 2013), as methods might be based on different source of information (e.g., single-SNP allele frequency or haplotype structure). This grouping strategy facilitates the interpretation of the signatures of selection results and minimizes sources of variations due to selection for other group of traits.
Another important point is the use or not of outgroups when performing hapFLK analyses. The different genomic regions identified in both hapFLK scenarios might be due to the fact that an outgroup can provide the connection between groups of populations. As discussed by Brito et al. (2017), the use of single-SNP tests may fail to identify signals of selection when a single SNP is not in high LD with the causal mutation for the trait under selection. Another problem that may arise is that hapFLK and ROH may not capture ancient signatures of selection, for which the mutation-carrying haplotype is small and do not include many SNP on the SNP chip panel.
Despite the limitation of the ascertainment bias and reduced genomic coverage of the SNP dataset used, an interesting list of potential genes under selection in multiple cattle breeds, YAK, and BLI, which will be the foundation for future investigations could be still provided. For instance, the RBPMS2 gene, found in YAK is related to age at first calving in Romosinuano and Costeño-con-Cuernos cattle breeds (Fernández et al., 2019) and the gene ATG16L2 found in MG cattle is related to mastitis recovery in Danish HOL cows (Welderufael et al., 2017). Angrecka and Herbut (2016) observed that cows in excessive cold temperatures may show fertility impairments and the potential candidate gene RBPMS2 found in YAK might be an indication of adaptation for the cold environments. Shathele (2009) observed a clear pattern of higher incidence of mastitis in dairy cows when the ambient temperature was lower than 21°C. The candidate gene ATG16L2 found in the MG breed can be an indicator of cold adaptation. Also, the candidate gene IL6ST identified in KZK cattle, may be another indicator of fertility adaptation in cold environment, since this gene is implicated in uterine receptivity to conceptus implantation in sheep (Song et al., 2009). Two genes (WARS2 and TBX15) were identified in ROH and were demonstrated by previous studies (Fumagalli et al., 2015; Hallmark et al., 2018) as having an important role in cold adaptation in Siberia Inuit (human) populations. Polymorphisms in or near WARS2 and TBX15 have been shown to be associated with numerous human phenotypes among individuals of European ancestry, including waist-hip ratio and fat distribution in Europeans (Heid et al., 2010). TBX15 gene plays a role in the differentiation of brown (subcutaneous) and brite (typically inguinal) adipocytes, that in the case of exposure to cold can differentiate into cells capable of expressing uncoupling protein 1, producing heat by lipid oxidation (Gburcik et al., 2012). Therefore, WARS2 and TBX15 are candidate genes that might be associated with cold adaptation in both cattle and YAK.
The candidate gene PLA2G2A identified in the heat group in FST1 is involved in lipid metabolism and was previously reported by Hallmark et al. (2018) in human inhabitants of Siberia, in which expression patterns was consistent with selection for cold climate and/or diet. The difference in adaptation environments found in this study and in theirs suggests that this candidate gene (PLA2G2A) may also be positively selected for lipid metabolism in warm environments. The gene TMC6, identified for mild cold/thermoneutral group based on FST3 and ROH, has been reported to influence milk quality in cold environment in Northern China breeds (Yang et al., 2017). Thus, this gene may play an important role in adaptation to cold climate. Yang et al. (2017) also identified the gene COL27A1 (found here in FST2 and ROH), which has been proposed to be involved in cartilage calcification, which is another important aspect in cold climates. Studying signatures of selection for candidate genes for body temperature maintenance under the cold stress in Siberian cattle populations, Igoshin et al. (2019) found that the gene GRIA4 contribute to cold-stress resistance due to indirect involvement in body thermoregulation. In our study a gene from the same family (GRIA1) was found for YAK in the FST3 and ROH scenarios. This indicates that GRIA1 may also be involved in cattle adaptation to cold climates. Howard et al. (2013) performed a GWAS for body temperature regulation in cattle during a five-day period of heat and cold stresses and identified candidate genes involved in vasculogenesis (RASA1) and ion regulation (CACNG3, SLC22A23, and TRPC4), which were also identified in our study for scenarios FST2 and ROH (RASA1), and FST3 and ROH (CACNG3), FST3 and ROH (SLC22A23), and FST1 and ROH (TRPC4).
High altitudes, where YAKs are raised, might result in hematopoiesia, where the reduced oxygen tension leads to an increase of erythrocytes as an adaptive mechanism to low oxygen levels (Jain, 1993; Ding et al., 2014). The candidate genes HBM and HBG, related with oxygen transportation were identified for YAK in FST1, FST2, and ROH (HBM) and in hapFLK and ROH (HBG). These genes were previous identified by Mehla et al. (2014) in a GWAS of heat stress response in Zebu (Sahiwal) cattle. In their study, HBM and HBG were suppressed in the presence of heat stress. The fact that these genes were found here as under positive selection in YAK highlights the importance of these genes in oxygen transportation in harsh and high-altitude environments.
The genes SLC38A4, AGXT2, STEAP3, CSMD3, ABCB4, ZDBF2, and LAMA2 identified in FST1 and ROH, were previously reported in a study of functional characterization of mammalian genes involved in adaptation to Arctic or Antarctic environments (Yudin et al., 2017) and had traces of positive selection in human, mammoth, and polar bear (SLC38A4), mammoth and polar bear (AGXT2, STEAP3, CSMD3 and ABCB4), mammoth, whale, and polar bear (ZDBF2), and human and mammoth (LAMA2). The gene SLC38A4 encodes a transporter protein that is expressed predominantly in the liver and participates in the transfer of amino acids (Bröer, 2013) and is potentially involved in gluconeogenesis (González-Renteria et al., 2013). Regarding the genes AGXT2 and STEAP3, Yudin et al. (2017) reported that these genes had two enriched terms related to regulation of reactive oxygen species metabolic process and transition metal ion transport. Similar to the AGXT2 and STEAP3 genes, Yudin et al. (2017) reported ABCB4 to be associated with the enriched term bile acid secretion, and likely related to fat-rich diets of mammals in the Arctic and Antarctic regions. The gene ZDBF2 was paternally expressed in lymphocytes, but bi-allelically expressed in the placenta in humans (Kobayashi et al., 2009) and related to various genetic abnormalities (Petit et al., 2004). The positive selection in the ZDBF2 gene might be related with physiological or immunological adaptations of animals to low temperatures. These findings represent a strong indication of the role of these genes (SLC38A4, AGXT2, STEAP3, ABCB4, ZDBF2, and LAMA2) in the metabolism and resistance of animals exposed to harsh environments (e.g., low temperatures and high altitudes).
The stress caused by high temperatures has a great impact on the animal’s life and can induce behavioral and metabolic changes in cattle that are intended to maintain homeothermy, often at the expense of decreased productivity, welfare, health, and profitability (West, 1994; Gao S. T. et al., 2017). In Weiling Humped cattle, which is adapted to hot environments, the gene HSPE1, which regulates maintenance of metabolic activity and survival (Kumar et al., 2015) was identified. Pitt et al. (2018), studying the local adaptation of Creole cattle genome diversity in the tropics, reported an association between SNORA70 with feed intake, body conformation, and live weight. Here this gene was identified in populations adapted either to cold or heat conditions, indicating its importance for metabolic adaptation in both thermal conditions. In DH and DC cattle, both adapted to hot environments, the gene EIF4E was identified, which is associated with conceptus development and/or return to cyclicity in cattle (Forde et al., 2011). The gene SNORD116 identified in DC cattle, has been reported to regulate food intake and body weight in mice (Qi et al., 2016). Three genes (LOX, HSPB9, and CSF3) identified based on FST1, FST3, and ROH in DH were previous reported to harbor selective sweeps for environmental adaptation in EAZ (Bahbahani et al., 2017). Another identified gene was the DNAJC8 which is a member of the DnaJ family which is known to act as cofactors for other heat shock proteins (i.e., HSP70) to maintain protein folding under heat stress (Coleman et al., 1995; Kampinga and Craig, 2010). The LOX gene, on the other hand, is important to maintain optimum growth and development and it is associated with lung and blood vessels (Mäki et al., 2005). The CSF3 gene was related with innate and adaptive immune response in African cattle, as this gene acts as a positive regulator for macrophages to induce their antimicrobial effects (Bronsvoort et al., 2013; Bahbahani et al., 2017). The presence of this gene in African cattle was expected since they are exposed to a diversity of pathogens and physiological stressors in their surrounding environment, e.g., endoparasites, haemoparasites, and bacteria (Bronsvoort et al., 2013; Murray et al., 2013; Thumbi et al., 2014).
Two genes identified between FST1 and ROH (NRXN3 and CACNB2) were previously reported to be related to fertility traits in cattle (Reverter et al., 2017). NRXN3 harbors a large number of mutations associated with heifer puberty in cattle (Dias et al., 2017). CACNB2 is related to the secretion of the follicle stimulating hormone, and variants of this gene have been found to be associated with conception rate in cows (Sugimoto et al., 2013). Reverter et al. (2017) also identified XKR4 as an important candidate gene for climate adaptation in Australian beef cattle. This gene (XKR4) was previously reported to be associated with fat deposition (Porto-Neto et al., 2014; Silva et al., 2017), feed intake and growth (Lindholm-Perry et al., 2011) phenotypes. The genes SLC9A4, PLCB1, FTO, ITPR2 identified in several scenarios for the heat adaptation group were previously reported by Taye et al. (2017) as under positive selection for thermotolerance in African cattle. SLC9A4 and ITPR2 are involved in cellular process in thermal sweating in African cattle. ITPR2 promotes the release of calcium from extracellular interstitial fluid and release of intracellular Ca2+ stores needed for normal sweat production (Cui and Schlessinger, 2015) and SLC9A4 is Na+/H+ exchangers that are expressed in the duct and secretary sweat glands in humans (Cui and Schlessinger, 2015). Another gene identified by Taye et al. (2017) is the PLCB1 gene, which hydrolyzes phospholipids into fatty acids and other lipophilic molecules. PLCB1 was previously reported as positively selected in sheep and goats in response to adaptation to dry-arid environments (Kim et al., 2015). Taye et al. (2017) also found the FTO gene as related to metabolism and feed intake, being FTO a nuclear protein of the AlkB related nonhaem iron and 2-oxoglutarate-dependent oxygenase superfamily, which is known to contribute to the regulation of the global metabolic rate, energy expenditure, energy homeostasis and feed intake (Zhang et al., 2011; Zielke et al., 2013). The above-mentioned genes (SLC9A4, PLCB1, FTO, ITPR2) might contribute to improved production and reproduction under heat stress conditions.
Immune and anti−infectious response related genes were identified in FST1 and ROH (LEF1 and MLST8), and FST2, FST3, and ROH (SMYD3). Flori et al. (2019) indicated that LEF1, MLST8, and SMYD3 are genes related with immune and anti−infectious response in Mediterranean cattle breeds. LEF1 could be a strong candidate gene to heat adaptation since it is also related to thermotolerance and ultra-violet protection, as suggested by its central role in hair pigmentation (Gautier et al., 2009). Studying environmental adaptation of indigenous Chinese cattle, Gao Y. et al. (2017) also identified the gene LEF1 near several known color genes. TICAM1 and CXCR4 genes (identified in FST1 and ROH) are involved in antigen recognition and related to the development of immune response in West African cattle (Gautier et al., 2009). TICAM1 possesses a high interferon type I inducing activity during viral infection and thus plays a role in innate immunity (Seya et al., 2005; Gautier et al., 2009) and CXCR4 is related to antigen presentation to T-cells in the context of major histocompatibility complex, a dense region of genes in the mammalian genome involved in the immune system, autoimmunity, and reproductive success (Contento et al., 2008; Gautier et al., 2009). The presence of these genes related to the immune response in populations adapted to heat highlights the importance of these genes to the immune response of animals during heat stress events.
A slick coat is associated to higher thermotolerance and higher milk production on crossbred animals under tropical environment and can be considered as an indicator or indirect phenotype to some highly important production traits (Dikmen et al., 2014; Casas and Kehrli, 2016; Carabaño et al., 2017; Porto-Neto et al., 2018). A deletion within the last exon of PRLR gene that truncates a portion of the cytoplasmic domain of the protein is responsible for this very smooth coat and were first identified in the SEN breed (Flori et al., 2012; Littlejohn et al., 2014). The PRLR gene was also identified in our study (based on FST1 and ROH) in WN and SEN, breeds with short coat hair. This finding in WN might be linked to thermotolerance in this cattle breed, an adaptation to high temperature environment. Another important trait in thermal tolerance adaptation in livestock is coat color (Mcmanus et al., 2008; Abdurehman, 2019). Cows with light coat colors in tropical regions reflect solar radiation; thereby protecting the animal from the adverse effects of solar radiation (Fanta, 2017; Gaughan et al., 2018). The gene PMEL, identified in FST1 and ROH in DC cattle (which is characterized by a yellow coat; Table 1 and Figure 9), is essential for melanosome development and is responsible for lightening or diluting the base color defined by the Extension locus (MC1R) in some cattle breeds (Schmutz and Dreger, 2012; Flori et al., 2014). The positive selection identified in DC cattle might be due to the association of this gene (PMEL) with thermotolerance.
Figure 9. Local Chinese cattle populations and yak: (a) Yunnan Humped; (b) Dengchuan; (c) Wenling Humped; (d) Wannan; (e) Sanhe; (f) Mongolian; (g) Hazake; (h) Xinjiang Brown; and (i) Datong Yak. Photo credits: Animal genetic resources in China: Bovines.
Heat shock factors have the function to act as chaperones assisting in protein folding, thereby avoiding protein aggregation and, as result, protein homeostasis during cellular response to heat stress (Moseley, 1997; Hightower, 2004). In our study, gene members of the heat shock family were identified (DNAJA1, DNAJB1, HSPH1, and HSPA4). Srikanth et al. (2016), studying genes and pathways linked to heat stress in HOL calves, identified 20 upregulated genes encoding molecular chaperones in response to heat stress. Among them, the members of HSP40 (DNAJA1 and DNAJB1) that regulates the ATPase activity of HSP70 (HSPH1 and HSPA4) by interacting with the J domain of the HSP70 proteins and also conferring thermo-tolerance to cells on exposure to heat stress (Maio, 1999). O’Brien et al. (2014), assessing signatures of selection through variation in linkage disequilibrium between taurine and indicine cattle, found the gene DNAJA1 in BSW and ANG cattle. A previous study (Bernard et al., 2007) has reported DNAJA1 to be linked to meat tenderness, since this gene explained up to 63% of the phenotypic variability of tenderness in Charolais cattle. In our study, DNAJA1 was also found for SH cattle, which are animals with good meat quality, and thus, DNAJA1 might be a good indicator of selection for meat tenderness in this breed.
The high number of genes identified in our study (Tables 5, 6 and Supplementary Material) makes it challenging to prioritize the most important candidate genes. The majority of genes were linked to resilience, adaptation response, health, fertility, and production. Important biological pathways and functions have been identified and have important functions related to adaptation to thermal stress and can help to identify possible candidate genes. Four genes (ALB, HBM, HBZ, and HBQ1) were clustered in the “oxygen binding” molecular function that was identified for YAK in scenarios SCEN1 and SCEN4. At high altitudes, where YAKs are raised, immediate means of acclimatization such as hyperventilation, hemoconcentration, and erythropoiesis are initiated to increase total oxygen carrying capacity (Burtscher et al., 2018). In such environments, it is observed increased oxygen binding by hemoglobin even at low arterial partial oxygen pressure (Mairbäurl and Weber, 2012; Burtscher et al., 2018). In addition, the candidate genes HBM, HBZ, and HBQ1 were grouped in other two clusters identified for YAK, which are related to oxygen transport and hemoglobin complex, indicating that these three genes may be of interest when selecting animals adapted to high altitudes where the oxygen level is lower and biological adaptations are required. Three candidate genes (FOLR1, FOLR2, and FOLR3) were grouped into two groups related to folate receptor and folic acid transporter activity. Koyama et al. (2012), studying the effects of folic acid on the development and oxidative stress of mouse embryos exposed to heat stress found that the inhibitory effects of heat stress on the development of mouse pre-implantation embryos are ameliorated by folic acid.
Heat stress can be considered as cytotoxic, as it alters biological molecules, disturbs cell functions, modulates metabolic reactions, induces oxidative cell damage, and activates both apoptosis and necrosis pathways (Pandey et al., 2012; Slimen et al., 2015). Four genes (NDUFS3, NDUFB6, CYCS, and CYCT) were clustered together and are related to respiratory chain. Ganaie et al. (2013) and Nisar et al. (2013) suggested that heat stress is responsible for inducing oxidative stress during summer in livestock animals, which enhances reactive oxygen species production and induces oxidative stress, what can lead to cytotoxicity (Slimen et al., 2015). The identification of these four genes related to respiratory chain in the scenarios of animals adapted to heat stress might indicate that these genes play an important role in preventing oxidative stress and cytotoxicity caused by high temperatures. High temperature is also considered to be one of the most important factors for subfertility in cattle (Takahashi, 2011; Das et al., 2016). Three genes (BSP1, BSP3, and BSP5) related to fertilization and sperm capacitation were identified in scenario SCEN2, based on the heat adaptation candidate genes group. Manjunath and Thérien (2002) reported the binder of sperm protein (BSP) family as a key step to sperm capacitation in ejaculation, as BSP bind to sperm and induce efflux of membrane cholesterol and phospholipids (Therrien et al., 2013). The BSP family genes identified here might indicate that under heat stress these genes play an important role on sperm capacitation and consequently, bull fertility.
The large number of genomic regions and candidate genes identified in this study was not unexpected. As there are multiple mechanisms involved in cold and heat stress response, natural and artificial selection are expected to cause weak selection pressure at many sites across the genome. The weak selection pressure is due to the small effect of each of the many variants or QTL underlying the phenotype, which results in small changes to allele frequencies at the many loci, but which result in the adaptation response observed across breeds (Kemper et al., 2014).
The information reported here will be useful for the implementation of genomic selection and other genomic studies in local Chinese cattle populations and YAK. Genomic selection is expected to increase the genetic gain along with the decrease in generation interval (Farah et al., 2015; Kariuki et al., 2017; Mehrban et al., 2017; Weller et al., 2017; Doublet et al., 2019). Due to the decrease in generation interval, a decrease in genetic diversity could be expected (Doublet et al., 2019). Eynard et al. (2018) indicated that optimal contribution strategy is an alternative to avoid the loss of genetic diversity and obtain genetic gain, as it takes into account both the relationships and genetic merit of the individuals. Further investigation using additional breeds, more complete datasets (e.g., larger number of breeds and phenotypes), and larger numbers of SNP are recommended to validate the role and the specific function of the highlighted candidate genes associated with climate adaptation in this study. Alternative “-omic” technologies (e.g., transcriptomics) and functional studies should be conducted to validate the function of the candidate genes identified.
Lastly, some of the potential limitations of this study should be highlighted. Having additional breeds from the same subspecies, but adapted to different climatic conditions, would be beneficial. Ascertainment bias might have influenced the estimates due the fact the SNP array design did not include all the cattle populations evaluated and especially any population of YAK or BLI. Moreover, due to the large number of candidate genes identified, additional studies should prioritize specific genes for more in-depth studies for validation purposes.
Conclusion
A comprehensive description of genetic diversity measures in thirty-two worldwide cattle populations (with a focus on Chinese local breeds), Datong yak, and Bali populations is presented. Moderate genetic diversity was observed within each population, but optimal contribution strategies should be adopted to avoid further reduction in the genetic diversity of these populations. The breeds adapted to heat seems to be moderately related to each other and the same was observed for those adapted to cold environments. This study provides a high-confidence map of positive selection signatures for thermal tolerance in cattle, and Datong yak using three alternative methods (FST, hapFLK, and ROH). Various candidate genes were identified as potentially under selection for thermal tolerance based on multiple methods and scenarios. Enrichment analyses identified important biological pathways, molecular function and cellular components, which greatly contribute to our understanding of the genetic background of thermal tolerance in cattle and Datong yak. These results highlight the diversity and complexity of mechanisms that are involved in response to harsh environmental conditions. Further studies are needed to validate these results by using larger and independent populations, additional methods for specific thermal stress indicators (e.g., GWAS), and other tools such as whole-genome sequencing. These novel findings are useful for the conservation of genetic resources, especially under global climate changes, and contribute to the potential discovery of functional genes and genetic mechanisms related to thermal tolerance in cattle, Datong yak, and Bali populations.
Data Availability Statement
All relevant information supporting the results of this article are included within the article and its additional files. Public genotypic data used for our analysis can be found at WIDDE database: http://widde.toulouse.inra.fr/widde/widde/main.do;jsessionid=73509E5BFFDDE224635B91CA70E9FC82?module=cattle.
Author Contributions
LB, YW, QX, and FS participated in the design of the study. PF performed all the analyses and results interpretation, was involved in the discussions, and wrote the first draft of the manuscript. LB provided training to the first author, contributed to the analyses and results interpretation, was involved in all the discussions, and edited all the versions of the manuscript. YW and YZ provided the cattle genotypes, while PY provided the YAK genotypes used for the study. HO participated in the results interpretation, was involved in the discussions and gave editorial assistance. FS, YW, QX, and PY edited the final version of the manuscript and provided valuable suggestions. All authors have approved the final manuscript.
Funding
This work was supported by the Research Fund for International Young Scientists by the National Natural Science Foundation of China (Grant Number: 31750110459), China Agriculture Research System (CARS-36), the Program for Changjiang Scholar, and Innovation Research Team in University (IRT_15R62).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors are grateful to Prof. Xixia Huang for sharing the genotype data of XIN cattle. YZ was supported by National Natural Scientific Foundation of China (31561143010). FS thanks the support from the research program Food from Thought: Agricultural Systems for a Healthy Planet led by the University of Guelph and funded in part by the Canada First Research Excellence Fund.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.604823/full#supplementary-material
Footnotes
References
Abdurehman, A. (2019). Physiological and Anatomical Adaptation Characteristics of Borana Cattle to Pastoralist Lowland Environments. Asian J. Biol. Sci. 12, 364–372. doi: 10.3923/ajbs.2019.364.372
Abied, A., Bagadi, A., Bordbar, F., Pu, Y., Augustino, S. M., Xue, X., et al. (2020). Genomic Diversity, Population Structure, and Signature of Selection in Five Chinese Native Sheep Breeds Adapted to Extreme Environments. Genes 11:494. doi: 10.3390/genes11050494
Ai, H., Huang, L., and Ren, J. (2013). Genetic Diversity, Linkage Disequilibrium and Selection Signatures in Chinese and Western Pigs Revealed by Genome-Wide SNP Markers. PLoS One 8:e56001. doi: 10.1371/journal.pone.0056001
Ai, H., Yang, B., Li, J., Xie, X., Chen, H., and Ren, J. (2014). Population history and genomic signatures for high-altitude adaptation in Tibetan pigs. BMC Genom. 15:834. doi: 10.1186/1471-2164-15-834
Ajmone-Marsan, P., Garcia, J. F., and Lenstra, J. A. (2010). On the origin of cattle: How aurochs became cattle and colonized the world. Evol. Anthropol. 19, 148–157. doi: 10.1002/evan.20267
Angrecka, S., and Herbut, P. (2016). Conditions for cold stress development in dairy cattle kept in free stall barn during severe frosts. Czech J. Anim. Sci. 60, 81–87. doi: 10.17221/7978-cjas
Bahbahani, H., Tijjani, A., Mukasa, C., Wragg, D., Almathen, F., Nash, O., et al. (2017). Signatures of Selection for Environmental Adaptation and Zebu × Taurine Hybrid Fitness in East African Shorthorn Zebu. Front. Genet. 8:68. doi: 10.3389/fgene.2017.00068
Beaumont, M. A. (2005). Adaptation and speciation: what can Fst tell us? Trends Ecol. Evolut. 20, 435–440. doi: 10.1016/j.tree.2005.05.017
Bernabucci, U. (2019). Climate change: impact on livestock and how can we adapt. Anim. Front. 9, 3–5. doi: 10.1093/af/vfy039
Bernard, C., Cassar-Malek, I., Cunff, M. L., Dubroeucq, H., Renand, G., and Hocquette, J.-F. (2007). New Indicators of Beef Sensory Quality Revealed by Expression of Specific Genes. J. Agricul. Food Chem. 55, 5229–5237. doi: 10.1021/jf063372l
Berry, D. P. (2017). Symposium review: Breeding a better cow-Will she be adaptable? J. Dairy Sci. 10, 3665–3685.
Boschiero, C., Moreira, G. C. M., Gheyas, A. A., Godoy, T. F., Gasparin, G., Mariani, P. D. S. C., et al. (2018). Genome-wide characterization of genetic variants and putative regions under selection in meat and egg-type chicken lines. BMC Genom. 19:83. doi: 10.1186/s12864-018-4444-0
Boushaba, N., Boujenane, I., Moazami-Goudarzi, K., Flori, L., Saïdi-Mehtar, N., Tabet-Aoul, N., et al. (2018). Genetic diversity and relationships among six local cattle populations in semi-arid areas assessed by a bovine medium-density single nucleotide polymorphism data. Animal 13, 8–14. doi: 10.1017/s1751731118001179
Brito, L. F., Kijas, J. W., Ventura, R. V., Sargolzaei, M., Porto-Neto, L. R., Cánovas, A., et al. (2017). Genetic diversity and signatures of selection in various goat breeds revealed by genome-wide SNP markers. BMC Genom. 18:229. doi: 10.1186/s12864-017-3610-0
Bröer, S. (2013). The SLC38 family of sodium–amino acid co-transporters. Eur. J. Physiol. 466, 155–172. doi: 10.1007/s00424-013-1393-y
Bronsvoort, B. M. D. C., Thumbi, S. M., Poole, E. J., Kiara, H., Auguet, O. T., Handel, I. G., et al. (2013). Design and descriptive epidemiology of the Infectious Diseases of East African Livestock (IDEAL) project, a longitudinal calf cohort study in western Kenya. BMC Vet. Res. 9:171. doi: 10.1186/1746-6148-9-171
Burtscher, M., Gatterer, H., Burtscher, J., and Mairbäurl, H. (2018). Extreme Terrestrial Environments: Life in Thermal Stress and Hypoxia. A Narrative Review. Front. Physiol. 9:572. doi: 10.3389/fphys.2018.00572
Carabaño, M. J., Ramón, M., Díaz, C., Molina, A., Pérez-Guzmán, M. D., and Serradilla, J. M. (2017). BREEDING AND GENETICS SYMPOSIUM: Breeding for resilience to heat stress effects in dairy ruminants. A comprehensive review. J. Anim. Sci. 95:1813. doi: 10.2527/jas2016.1114
Casas, E., and Kehrli, M. E. (2016). A Review of Selected Genes with Known Effects on Performance and Health of Cattle. Front. Vet. Sci. 3:113. doi: 10.3389/fvets.2016.00113
Ceballos, F. C., Joshi, P. K., Clark, D. W., Ramsay, M., and Wilson, J. F. (2018). Runs of homozygosity: windows into population history and trait architecture. Nat. News 4:19. Available online at: https://www.nature.com/articles/nrg.2017.109
Chang, C. C., Chow, C. C., Tellier, L. C., Vattikuti, S., Purcell, S. M., and Lee, J. J. (2015). Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience 4:7. doi: 10.1186/s13742-015-0047-8
Chen, M., Pan, D., Ren, H., Fu, J., Li, J., Su, G., et al. (2016). Identification of selective sweeps reveals divergent selection between Chinese Holstein and Simmental cattle populations. Genet. Select. Evol. 48:76. doi: 10.1186/s12711-016-0254-5
Chen, M., Sun, Y., Yang, C., Zeng, G., Li, Z., and Zhang, J. (2018). The road to wild yak protection in China. Science 360:866. doi: 10.1126/science.aat6749
Cheruiyot, E. K., Bett, R. C., Amimo, J. O., Zhang, Y., Mrode, R., and Mujibi, F. D. N. (2018). Signatures of Selection in Admixed Dairy Cattle in Tanzania. Frontiers 9:607. doi: 10.3389/fgene.2018.00607/full
Coleman, J. S., Heckathorn, S. A., and Hallberg, R. L. (1995). Heat-shock proteins and thermotolerance: linking molecular and ecological perspectives. Trends Ecol. Evolut. 10, 305–306.
Contento, R. L., Molon, B., Boularan, C., Pozzan, T., Manes, S., Marullo, S., et al. (2008). CXCR4–CCR5: A couple modulating T cell functions. PNAS 105, 10101–10106. doi: 10.1073/pnas.0804286105
Cui, C. Y., and Schlessinger, D. (2015). Eccrine sweat gland development and sweat secretion. Exp. Dermatol. 24, 644–50. doi: 10.1111/exd.12773
Das, R., Sailo, L., Verma, N., Bharti, P., Saikia, J., Imtiwati, et al. (2016). Impact of heat stress on health and performance of dairy animals: A review. Vet. World 9, 260–268. doi: 10.14202/vetworld.2016.260-268
Dias, M. M., Cánovas, A., Mantilla-Rojas, C., Riley, D. G., Luna-Nevarez, P., Coleman, S. J., et al. (2017). SNP detection using RNA-sequences of candidate genes associated with puberty in cattle. Genet. Mole. Res. 1:16.
Dikmen, S., Khan, F. A., Huson, H. J., Sonstegard, T. S., Moss, J. I., Dahl, G. E., et al. (2014). The SLICK hair locus derived from Senepol cattle confers thermotolerance to intensively managed lactating Holstein cows. J. Dairy Sci. 97, 5508–5520.
Ding, X. Z., Liang, C. N., Guo, X., Wu, X. Y., Wang, H. B., Johnson, K. A., et al. (2014). Physiological insight into the high-altitude adaptations in domesticated yaks (Bos grunniens) along the Qinghai-Tibetan Plateau altitudinal gradient. Livestock Sci. 162:2014.
Dixit, S. P., Singh, S., Ganguly, I., Bhatia, A. K., Sharma, A., Kumar, N. A., et al. (2020). Genome-Wide Runs of Homozygosity Revealed Selection Signatures in Bos indicus. Frontiers 11:92. doi: 10.3389/fgene.2020.00092/full
Doekes, H. P., Veerkamp, R. F., Bijma, P., de Jong, G., Hiemstra, S. J., and Windig, J. J. (2019). Inbreeding depression due to recent and ancient inbreeding in Dutch Holstein-Friesian dairy cattle. GSE. 51:54.
Doublet, A.-C., Croiseau, P., Fritz, S., Michenet, A., Hozé, C., Danchin-Burge, C., et al. (2019). The impact of genomic selection on genetic diversity and genetic gain in three French dairy cattle breeds. GSE 51:52. doi: 10.1186/s12711-019-0495-1
Edea, Z., Dadi, H., Dessie, T., Uzzaman, M. R., Rothschild, M. F., Kim, E.-S., et al. (2018). Genome-wide scan reveals divergent selection among taurine and zebu cattle populations from different regions. Anim. Genet. 49, 550–563. doi: 10.1111/age.12724
Edea, Z., Dadi, H., Kim, S.-W., Dessie, T., Lee, T., Kim, H., et al. (2013). Genetic diversity, population structure and relationships in indigenous cattle populations of Ethiopia and Korean Hanwoo breeds using SNP markers. Front. Genet. 4:35. doi: 10.3389/fgene.2013.00035
Eynard, S. E., Croiseau, P., Laloë, D., Fritz, S., Calus, M. P. L., and Restoux, G. (2018). Which Individuals To Choose To Update the Reference Population? Minimizing the Loss of Genetic Diversity in Animal Genomic Selection Programs. G3, 113–121. Available online at: https://www.g3journal.org/content/8/1/113
Fanta, M. (2017). Physiological adaptation of Holstein Frisian dairy cattle in Ethiopia: review article. J. Biol. Agric. Health. 7, 67–78.
Farah, M. M., Swan, A. A., Fortes, M. R. S., Fonseca, R., Moore, S. S., and Kelly, M. J. (2015). Accuracy of genomic selection for age at puberty in a multi−breed population of tropically adapted beef cattle. Anim. Genet. 47, 3–11. doi: 10.1111/age.12362
Fariello, M. I., Boitard, S., Naya, H., SanCristobal, M., and Servin, B. (2013). Detecting signatures of selection through haplotype differentiation among hierarchically structured populations. Genetics 193, 929–941.
Fernández, J., Pérez, J., Herrera, N., Martínez, R., Bejarano, D., and Rocha, J. (2019). Research Article Genomic association study for age at first calving and calving interval in Romosinuano and Costeño con Cuernos cattle. Genet. Mole. Res. 18:2019. doi: 10.4238/gmr18258
Flori, L., Gonzatti, M. I., Thevenon, S., Chantal, I., Pinto, J., Berthier, D., et al. (2012). A quasi-exclusive European ancestry in the Senepol tropical cattle breed highlights the importance of the slick locus in tropical adaptation. PLoS One 7:e36133. doi: 10.1371/journal.pone.0036133
Flori, L., Moazami−Goudarzi, K., Alary, V., Araba, A., Boujenane, I., Boushaba, N., et al. (2019). A genomic map of climate adaptation in Mediterranean cattle breeds. Mol. Ecol. 28, 1009–1029. doi: 10.1111/mec.15004
Flori, L., Thevenon, S., Dayo, G. K., Senou, M., Sylla, S., Berthier, D., et al. (2014). Adaptive admixture in the West African bovine hybrid zone: insight from the Borgou population. Mol Ecol. 23, 3241–3257. doi: 10.1111/mec.12816
Forde, N., Carter, F., Spencer, T. E., Bazer, F. W., Sandra, O., Mansouri-Attia, N., et al. (2011). Conceptus-Induced Changes in the Endometrial Transcriptome: How Soon Does the Cow Know She Is Pregnant? 1. OUP Acad. 85, 144–156.
Foury, F., and Kolodynski, J. (1983). pif mutation blocks recombination between mitochondrial rho and rho- genomes having tandemly arrayed repeat units in Saccharomyces cerevisiae. Proc. Natl. Acad. Sci. U. S. A. 80, 5345–5349.
Fumagalli, M., Moltke, I., Grarup, N., Racimo, F., Bjerregaard, P., Jørgensen, M. E., et al. (2015). Greenlandic Inuit show genetic signatures of diet and climate adaptation. Science 349, 1343–1347. doi: 10.1126/science.aab2319
Ganaie, A. H., Ghasura, R. S., Mir, N. A., Bumla, N. A., Sankar, G., and Wani, S. A. (2013). Biochemical and Physiological Changes during Thermal Stress in Bovines: A Review. Iranian J. Appl. Anim. Sci. 3, 423–430.
Gao, S. T., Guo, J., Quan, S. Y., Nan, X. M., Fernandez, M. V. S., Baumgard, L. H., et al. (2017). The effects of heat stress on protein metabolism in lactating Holstein cows. J. Dairy Sci. 100, 5040–5049.
Gao, Y., Gautier, M., Ding, X., Zhang, H., Wang, Y., Wang, X., et al. (2017). Species composition and environmental adaptation of indigenous Chinese cattle. Nat. News. 7:16196. doi: 10.1038/s41598-017-16438-7
Gaughan, J. B., Sejian, V., Mader, T. L., and Dunshea, F. R. (2018). Adaptation strategies: ruminants. OUP Acad. 9, 47–53.
Gautier, M., Flori, L., Riebler, A., Jaffrézic, F., Laloé, D., Gut, I., et al. (2009). A whole genome Bayesian scan for adaptive genetic divergence in West African cattle. BMC Genom. 10:550. doi: 10.1186/1471-2164-10-550
Gburcik, V., Cawthorn, W. P., Nedergaard, J., Timmons, J. A., and Cannon, B. (2012). An essential role for Tbx15 in the differentiation of brown and “brite” but not white adipocytes. Am. J. Physiol. Endocrinol. Metabol. 303, E1053–E1060.
Geibel, J., Reimer, C., Weigend, S., Weigend, A., Pook, T., and Simianer, H. (2019). How Array Design Affects SNP Ascertainment Bias. bioRxiv. doi: 10.1101/833541v1
Glazar, A. I., and Evans, J. P. (2009). Immunoglobulin superfamily member IgSF8 (EWI-2) and CD9 in fertilisation: evidence of distinct functions for CD9 and a CD9-associated protein in mammalian sperm-egg interaction. Reprod. Fert. Devel. 21, 293–303. doi: 10.1071/rd08158
González-Renteria, S. M., Loera-Castañeda, V., Chairez-Hernández, I., Sosa-Macias, M., Paniagua-Castro, N., Lares-Aseff, I., et al. (2013). Association of the polymorphisms 292C>T and 1304G>A in the SLC38A4 gene with hyperglycaemia. Diab. Metabol. Res. Rev. 29, 39–43. doi: 10.1002/dmrr.2344
Hallmark, B., Karafet, T. M., Hsieh, P. H., Osipova, L. P., Watkins, J. C., and Hammer, M. F. (2018). Genomic Evidence of Local Adaptation to Climate and Diet in Indigenous Siberians. OUP Acad. 36, 315–327. doi: 10.1093/molbev/msy211
Hayes, B. J., Visscher, P. M., McPartlan, H. C., and Goddard, M. E. (2003). Novel multilocus measure of linkage disequilibrium to estimate past effective population size. Genome Res. 13, 635–643.
Heid, I. M., Jackson, A. U., Randall, J. C., Winkler, T. W., Qi, L., Steinthorsdottir, V., et al. (2010). Meta-analysis identifies 13 new loci associated with waist-hip ratio and reveals sexual dimorphism in the genetic basis of fat distribution. Nat. News 42, 949–960.
Henry, B. K., Eckard, R. J., and Beauchemin, K. A. (2018). Review: Adaptation of ruminant livestock production systems to climate changes. Animal 12, s445–s456. doi: 10.1017/S1751731118001301
Hiemstra, S. J., Has, Y., Mäki-Tanila, A., and Gandini, G. (2010). Local cattle breeds in Europe: development of policies and strategies for self-sustaining breeds.
Hightower, L. E. (2004). Heat shock, stress proteins, chaperones, and proteotoxicity. Cell 66, 191–197. doi: 10.1016/0092-8674(91)90611-2
Hoffmann, I. (2010). Climate change and the characterization, breeding and conservation of animal genetic resources. Anim. Genet. 41, 32–46. doi: 10.1111/j.1365-2052.2010.02043.x
Hohenlohe, P. A., Phillips, P. C., and Cresko, W. A. (2010). Using Population Genomics to Detect Selection in Natural Populations: Key Concepts and Methodological Considerations. Int. J. Plant Sci. 171, 1059–1071. doi: 10.1086/656306
Howard, J. T., Kachman, S. D., Snelling, W. M., Pollak, E. J., Ciobanu, D. C., Kuehn, L. A., et al. (2013). Beef cattle body temperature during climatic stress: a genome-wide association study. Int. J. Biometeorol. 58, 1665–1672. doi: 10.1007/s00484-013-0773-5
Hu, Z.-L., Park, C. A., and Reecy, J. M. (2018). Building a livestock genetic and genomic information knowledgebase through integrative developments of Animal QTLdb and CorrDB. Nucleic Acids Res. 47:1084. doi: 10.1093/nar/gky1084
Huang, D. W., Sherman, B. T., and Lempicki, R. A. (2008a). Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 37, 1–13. doi: 10.1093/nar/gkn923
Huang, D. W., Sherman, B. T., and Lempicki, R. A. (2008b). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57. doi: 10.1038/nprot.2008.211
Igoshin, A. V., Yurchenko, A. A., Belonogova, N. M., Petrovsky, D. V., Aitnazarov, R. B., Soloshenko, V. A., et al. (2019). Genome-wide association study and scan for signatures of selection point to candidate genes for body temperature maintenance under the cold stress in Siberian cattle populations. BMC Genet. 20:26. doi: 10.1186/s12863-019-0725-0
Kampinga, H. H., and Craig, E. A. (2010). Erratum: The HSP70 chaperone machinery: J proteins as drivers of functional specificity. Nat. Rev. Mole. Cell Biol. 11, 750–750. doi: 10.1038/nrm2982
Kanehisa, M. (2000). KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 28, 27–30. doi: 10.1093/nar/28.1.27
Kanehisa, M., Furumichi, M., Tanabe, M., Sato, Y., and Morishima, K. (2016). KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 45, D353–D361. doi: 10.1093/nar/gkw1092
Kanehisa, M., Sato, Y., Furumichi, M., Morishima, K., and Tanabe, M. (2018). New approach for understanding genome variations in KEGG. Nucleic Acids Res. 47, D590–D595. doi: 10.1093/nar/gky962
Kariuki, C., Brascamp, E., Komen, H., Kahi, A., and Arendonk, J. V. (2017). Economic evaluation of progeny-testing and genomic selection schemes for small-sized nucleus dairy cattle breeding programs in developing countries. J. Dairy Sci. 100, 2258–2268. doi: 10.3168/jds.2016-11816
Kemper, K. E., Saxton, S. J., Bolormaa, S., Hayes, B. J., and Goddard, M. E. (2014). Selection for complex traits leaves little or no classic signatures of selection. BMC Genom. 15:246. doi: 10.1186/1471-2164-15-246
Kim, E.-S., Elbeltagy, A. R., Aboul-Naga, A. M., Rischkowsky, B., Sayre, B., Mwacharo, J. M., et al. (2015). Multiple genomic signatures of selection in goats and sheep indigenous to a hot arid environment. Heredity 116, 255–264. doi: 10.1038/hdy.2015.94
Kinsella, R. J., Kahari, A., Haider, S., Zamora, J., Proctor, G., Spudich, G., et al. (2011). Ensembl BioMarts: a hub for data retrieval across taxonomic space. Database 2011:bar030. doi: 10.1093/database/bar030
Kobayashi, H., Yamada, K., Morita, S., Hiura, H., Fukuda, A., Kagami, M., et al. (2009). Identification of the mouse paternally expressed imprinted gene Zdbf2 on chromosome 1 and its imprinted human homolog ZDBF2 on chromosome 2. Genomics 93, 461–472. doi: 10.1016/j.ygeno.2008.12.012
Koyama, H., Ikeda, S., Sugimoto, M., and Kume, S. (2012). Effects of Folic Acid on the Development and Oxidative Stress of Mouse Embryos Exposed to Heat Stress. Reproduct. Domest. Anim. 47, 921–927. doi: 10.1111/j.1439-0531.2012.01992.x
Kristensen, T. N., Hoffmann, A. A., Pertoldi, C., and Stronen, A. V. (2015). What can livestock breeders learn from conservation genetics and vice versa? Front. Genet. 6:38. doi: 10.3389/fgene.2015.00038
Ku, C. S., Naidoo, N., Teo, S. M., and Pawitan, Y. (2010). Regions of homozygosity and their impact on complex diseases and traits. Hum. Genet. 129, 1–15. doi: 10.1007/s00439-010-0920-6
Kukučková, V., Moravčíková, N., Ferenčaković, M., Simčič, M., Mčszáros, G., Sölkner, J., et al. (2017). Genomic characterization of Pinzgau cattle: genetic conservation and breeding perspectives. Conservat. Genet. 18, 893–910. doi: 10.1007/s10592-017-0935-9
Kullo, I. J., and Ding, K. (2007). Patterns of population differentiation of candidate genes for cardiovascular disease. BMC Genet. 8:48. doi: 10.1186/1471-2156-8-48
Kumar, A., Ashraf, S., Goud, T. S., Grewal, A., Singh, S., Yadav, B., et al. (2015). Expression profiling of major heat shock protein genes during different seasons in cattle (Bos indicus) and buffalo (Bubalus bubalis) under tropical climatic condition. J. Therm. Biol. 51, 55–64. doi: 10.1016/j.jtherbio.2015.03.006
Lachance, J., and Tishkoff, S. A. (2013). SNP ascertainment bias in population genetic analyses: Why it is important, and how to correct it. BioEssays 35, 780–786. doi: 10.1002/bies.201300014
Lei, C. Z., Chen, H., Zhang, H. C., Cai, X., Liu, R. Y., Luo, L. Y., et al. (2006). Origin and phylogeographical structure of Chinese cattle. Anim. Genet. 37, 579–582. doi: 10.1111/j.1365-2052.2006.01524.x
Leroy, G. (2014). Inbreeding depression in livestock species: review and meta-analysis. Anim. Genet. 45, 618–628. doi: 10.1111/age.12178
Li, R., Zhang, X. M., Campana, M. G., Huang, J. P., Chang, Z. H., Qi, X. B., et al. (2013). Paternal origins of Chinese cattle. Anim. Genet. 44, 446–449. doi: 10.1111/age.12022
Lin, B. Z., Sasazaki, S., and Mannen, H. (2010). Genetic diversity and structure in Bos taurus and Bos indicus populations analyzed by SNP markers. Anim. Sci. J. 81, 281–289. doi: 10.1111/j.1740-0929.2010.00744.x
Lindholm-Perry, A. K., Kuehn, L. A., Smith, T. P. L., Ferrell, C. L., Jenkins, T. G., Freetly, H. C., et al. (2011). A region on BTA14 that includes the positional candidate genes LYPLA1, XKR4 and TMEM68 is associated with feed intake and growth phenotypes in cattle1. Anim. Genet. 43, 216–219. doi: 10.1111/j.1365-2052.2011.02232.x
Littlejohn, M. D., Henty, K. M., Tiplady, K., Johnson, T., Harland, C., Lopdell, T., et al. (2014). Functionally reciprocal mutations of the prolactin signalling pathway define hairy and slick cattle. Nat. Commun. 5:5861. doi: 10.1038/ncomms6861
Maio, A. D. (1999). Heat Shock Proteins: Facts, Thoughts, And Dreams. Shock 11, 1–12. doi: 10.1097/00024382-199901000-00001
Maiorano, A. M., Lourenco, D. L., Tsuruta, S., Ospina, A. M. T., Stafuzza, N. B., Masuda, Y., et al. (2018). Assessing genetic architecture and signatures of selection of dual purpose Gir cattle populations using genomic information. PLoS One 13:e0200694. doi: 10.1371/journal.pone.0200694
Mairbäurl, H., and Weber, R. E. (2012). Oxygen Transport by Hemoglobin. Comprehen. Physiol. 2, 1463–1489. doi: 10.1002/cphy.c080113
Mäki, J. M., Sormunen, R., Lippo, S., Kaarteenaho-Wiik, R., Soininen, R., and Myllyharju, J. (2005). Lysyl Oxidase Is Essential for Normal Development and Function of the Respiratory System and for the Integrity of Elastic and Collagen Fibers in Various Tissues. Am. J. Pathol. 167, 927–936. doi: 10.1016/s0002-9440(10)61183-2
Makina, S. O., Muchadeyi, F. C., Marle-Köster, E. V., Taylor, J. F., Makgahlela, M. L., and Maiwashe, A. (2015). Genome-wide scan for selection signatures in six cattle breeds in South Africa. Genet. Select. Evol. 47:92. doi: 10.1186/s12711-015-0173-x
Malomane, D. K., Reimer, C., Weigend, S., Weigend, A., Sharifi, A. R., and Simianer, H. (2018). Efficiency of different strategies to mitigate ascertainment bias when using SNP panels in diversity studies. BMC Genom. 19:22. doi: 10.1186/s12864-017-4416-9
Manjunath, P., and Thérien, I. (2002). Role of seminal plasma phospholipid-binding proteins in sperm membrane lipid modification that occurs during capacitation. J. Reprod. Immunol. 53, 109–119. doi: 10.1016/s0165-0378(01)00098-5
Mastrangelo, S., Ciani, E., Marsan, P. A., Bagnato, A., Battaglini, L., Bozzi, R., et al. (2018). Conservation status and historical relatedness of Italian cattle breeds. Genet. Select. Evol. 50:2018. doi: 10.1186/s12711-018-0406-x
Mcmanus, C., Paludo, G. R., Louvandini, H., Gugel, R., Sasaki, L. C. B., and Paiva, S. R. (2008). Heat tolerance in Brazilian sheep: Physiological and blood parameters. Trop. Anim. Health Product. 41, 95–101. doi: 10.1007/s11250-008-9162-1
Mcrae, K. M., Mcewan, J. C., Dodds, K. G., and Gemmell, N. J. (2014). Signatures of selection in sheep bred for resistance or susceptibility to gastrointestinal nematodes. BMC Genom. 15:637. doi: 10.1186/1471-2164-15-637
Medugorac, I., Medugorac, A., Russ, I., Veit-Kensch, C. E., Taberlet, P., Luntz, B., et al. (2009). Genetic diversity of European cattle breeds highlights the conservation value of traditional unselected breeds with high effective population size. Mole. Ecol. 18, 3394–3410. doi: 10.1111/j.1365-294x.2009.04286.x
Mehla, K., Magotra, A., Choudhary, J., Singh, A., Mohanty, A., Upadhyay, R., et al. (2014). Genome-wide analysis of the heat stress response in Zebu (Sahiwal) cattle. Gene 533, 500–507. doi: 10.1016/j.gene.2013.09.051
Mehrban, H., Lee, D. H., Moradi, M. H., Ilcho, C., Naserkheil, M., and Ibáñez-Escriche, N. (2017). Predictive performance of genomic selection methods for carcass traits in Hanwoo beef cattle: impacts of the genetic architecture. Genet. Select. Evol. 49:1. doi: 10.1186/s12711-016-0283-0
Meuwissen, T. (2009). Genetic management of small populations: A review. Anim. Sci. 59, 71–79. doi: 10.1080/09064700903118148
Meuwissen, T. H. E., Sonesson, A. K., Gebregiwergis, G., and Woolliams, J. A. (2020). Management of Genetic Diversity in the Era of Genomics. Front. Genet. 11:880. doi: 10.3389/fgene.2020.00880
Michelizzi, V. N., Wu, X., Dodson, M. V., Michal, J. J., Zambrano-Varon, J., Mclean, D. J., et al. (2011). A Global View of 54,001 Single Nucleotide Polymorphisms (SNPs) on the Illumina BovineSNP50 BeadChip and Their Transferability to Water Buffalo. Int. J. Biol. Sci. 7, 18–27. doi: 10.7150/ijbs.7.18
Mokry, F., Buzanskas, M., Mudadu, M. D. A., Grossi, D. D. A., Higa, R., Ventura, R., et al. (2014). Linkage disequilibrium and haplotype block structure in a composite beef cattle breed. BMC Genom. 15:S6. doi: 10.1186/1471-2164-15-s7-s6
Moradi, M. H., Nejati-Javaremi, A., Moradi-Shahrbabak, M., Dodds, K. G., and Mcewan, J. C. (2012). Genomic scan of selective sweeps in thin and fat tail sheep breeds for identifying of candidate regions associated with fat deposition. BMC Genet. 13:10. doi: 10.1186/1471-2156-13-10
Moseley, P. L. (1997). Heat shock proteins and heat adaptation of the whole organism. J. Appl. Physiol. 83, 1413–1417. doi: 10.1152/jappl.1997.83.5.1413
Murray, G. G., Woolhouse, M. E., Tapio, M., Mbole-Kariuki, M. N., Sonstegard, T. S., Thumbi, S. M., et al. (2013). Genetic susceptibility to infectious disease in East African Shorthorn Zebu: a genome-wide analysis of the effect of heterozygosity and exotic introgression. BMC Evolut. Biol. 13:246. doi: 10.1186/1471-2148-13-246
Mwai, O., Hanotte, O., Kwon, Y.-J., and Cho, S. (2015). - Invited Review - African Indigenous Cattle: Unique Genetic Resources in a Rapidly Changing World. Asian Austral. J. Anim. Sci. 28, 911–921. doi: 10.5713/ajas.15.0002r
Naskar, S., Gowane, G. R., and Chopra, A. (2015). Strategies to Improve Livestock Genetic Resources to Counter Climate Change Impact. Adapt. Mitigat. 332, 441–475. doi: 10.1007/978-81-322-2265-1_25
Neupane, M., Geary, T. W., Kiser, J. N., Burns, G. W., Hansen, P. J., Spencer, T. E., et al. (2017). Loci and pathways associated with uterine capacity for pregnancy and fertility in beef cattle. PLoS One 12:e0188997. doi: 10.1371/journal.pone.0188997
Nguyen, T. T., Bowman, P. J., Haile-Mariam, M., Pryce, J. E., and Hayes, B. J. (2016). Genomic selection for tolerance to heat stress in Australian dairy cattle. J. Dairy Sci. 99, 2849–2862. doi: 10.3168/jds.2015-9685
Nisar, N. A., Sultana, M., Waiz, H. A., Para, P. A., and Dar, S. A. (2013). Oxidative stress threat to animal health and production. Int. J. Livest. Res. 3, 76–83.
Noor, R. R., Farajallaip, A., and Karmitai, M. (2001). Pengujian kemurnian sapi Bali dengan analisis hemoglobin dengan metode lsoelectric focusing. Hayati 8, 107–111.
O’Brien, A. M. P., Utsunomiya, Y. T., Mészáros, G., Bickhart, D. M., Liu, G. E., Tassell, C. P. V., et al. (2014). Assessing signatures of selection through variation in linkage disequilibrium between taurine and indicine cattle. Genet. Select. Evol. 46:19. doi: 10.1186/1297-9686-46-19
Oleksyk, T. K., Smith, M. W., and O’brien, S. J. (2010). Genome-wide scans for footprints of natural selection. Philosophical. Trans. Roy. Soc. B 365, 185–205. doi: 10.1098/rstb.2009.0219
Oleński, K., Tokarska, M., Hering, D., Puckowska, P., Ruść, A., Pertoldi, C., et al. (2015). Genome-wide association study for posthitis in the free-living population of European bison (Bison bonasus). Biol. Direct. 10:2. doi: 10.1186/s13062-014-0033-6
Olivieri, B. F., Mercadante, M. E. Z., Joslaine Noely, Dos Santos, Gonçalves Cyrillo, Branco, R. H., et al. (2016). Genomic Regions Associated with Feed Efficiency Indicator Traits in an Experimental Nellore Cattle Population. PLoS One 11:e0164390. doi: 10.1371/journal.pone.0164390
Onzima, R. B., Upadhyay, M. R., Doekes, H. P., Brito, L. F., Bosse, M., Kanis, E., et al. (2018). Genome-Wide Characterization of Selection Signatures and Runs of Homozygosity in Ugandan Goat Breeds. Front. Genet. 9:318. doi: 10.3389/fgene.2018.00318
Osei-Amponsah, R., Chauhan, S. S., Leury, B. J., Cheng, L., Cullen, B., Clarke, I. J., et al. (2019). Genetic Selection for Thermotolerance in Ruminants. Animals 9:948. doi: 10.3390/ani9110948
Pandey, N., Kataria, N., Kumar, K. A., Joshi, A., Narayan, S. L., Asopa, S., et al. (2012). Extreme ambiances vis−à−vis endogenous antioxidants of Marwari goat from arid tracts in India. ELBA Bioflux 4, 29–33.
Parland, S. M., Kearney, F., and Berry, D. P. (2009). Purging of inbreeding depression within the Irish Holstein-Friesian population. Genet. Select. Evolut. 41:16. doi: 10.1186/1297-9686-41-16
Pereira, R., Santana, M., Ayres, D., Bignardi, A., Menezes, G., Silva, L., et al. (2016). Inbreeding depression in Zebu cattle traits. J. Anim. Breed. Genet. 133, 523–533. doi: 10.1111/jbg.12219
Petit, F. M., Gajdos, V., Parisot, F., Capel, L., Aboura, A., Lachaux, A., et al. (2004). Paternal isodisomy for chromosome 2 as the cause of Crigler–Najjar type I syndrome. Eur. J. Hum. Genet. 13, 278–282. doi: 10.1038/sj.ejhg.5201342
Pitt, D., Bruford, M. W., Barbato, M., Orozco−Terwengel, P., Martínez, R., and Sevane, N. (2018). Demography and rapid local adaptation shape Creole cattle genome diversity in the tropics. Evol. Appl. 12, 105–122. doi: 10.1111/eva.12641
Polsky, L., and Keyserlingk, M. A. V. (2017). Invited review: Effects of heat stress on dairy cattle welfare. J. Dairy Sci. 100, 8645–8657. doi: 10.3168/jds.2017-12651
Porto-Neto, L. R., Bickhart, D. M., Landaeta-Hernandez, A. J., Utsunomiya, Y. T., Pagan, M., Jimenez, E., et al. (2018). Convergent Evolution of Slick Coat in Cattle through Truncation Mutations in the Prolactin Receptor. Front. Genet. 9:57. doi: 10.3389/fgene.2018.00057
Porto-Neto, L. R., Reverter, A., Prayaga, K. C., Chan, E. K. F., Johnston, D. J., Hawken, R. J., et al. (2014). The Genetic Architecture of Climatic Adaptation of Tropical Cattle. PLoS One 9:e113284. doi: 10.1371/journal.pone.0113284
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., et al. (2007). PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 81, 559–575. doi: 10.1086/519795
Purfield, D. C., Berry, D. P., Mcparland, S., and Bradley, D. G. (2012). Runs of homozygosity and population history in cattle. BMC Genet. 13:70. doi: 10.1186/1471-2156-13-70
Qi, Y., Purtell, L., Fu, M., Lee, N. J., Aepler, J., Zhang, L., et al. (2016). Snord116 is critical in the regulation of food intake and body weight. Scient. Rep. 6:18614. doi: 10.1038/srep18614
Qiu, Q., Wang, L., Wang, K., Yang, Y., Ma, T., Wang, Z., et al. (2015). Yak whole-genome resequencing reveals domestication signatures and prehistoric population expansions. Nat. Commun. 6:10283. doi: 10.1038/ncomms10283
Randhawa, I. A. S., Khatkar, M. S., Thomson, P. C., and Raadsma, H. W. (2016). A Meta-Assembly of Selection Signatures in Cattle. PLoS One 11:e0153013. doi: 10.1371/journal.pone.0153013
Randhawa, I. A., Khatkar, M., Thomson, P., and Raadsma, H. (2014). Composite selection signals can localize the trait specific genomic regions in multi-breed populations of cattle and sheep. BMC Genet. 15:34. doi: 10.1186/1471-2156-15-34
Reverter, A., Porto-Neto, L. R., Fortes, M. R. S., Kasarapu, P., Cara, M. A. R. D., Burrow, H. M., et al. (2017). Genomic inbreeding depression for climatic adaptation of tropical beef cattle. J. Anim. Sci. 95:3809. doi: 10.2527/jas2017.1643
Rodríguez-Ramilo, S., Elsen, J., and Legarra, A. (2019). Inbreeding and effective population size in French dairy sheep: Comparison between genomic and pedigree estimates. J. Dairy Sci. 102, 4227–4237. doi: 10.3168/jds.2018-15405
Roos, A. P. W. D., Hayes, B. J., and Goddard, M. E. (2009). Reliability of Genomic Predictions Across Multiple Populations. Genetics 183, 1545–1553. doi: 10.1534/genetics.109.104935
Rosen, B. D., Bickhart, D. M., Schnabel, R. D., et al. (2020). De novo assembly of the cattle reference genome with single-molecule sequencing. GigaScience 9:giaa021. doi: 10.1093/gigascience/giaa021
Rothammer, S., Seichter, D., Förster, M., and Medugorac, I. (2013). A genome-wide scan for signatures of differential artificial selection in ten cattle breeds. BMC Genom. 14:908. doi: 10.1186/1471-2164-14-908
Scheet, P., and Stephens, M. (2006). A Fast and Flexible Statistical Model for Large-Scale Population Genotype Data: Applications to Inferring Missing Genotypes and Haplotypic Phase. Am. J. Hum. Genet. 78, 629–644. doi: 10.1086/502802
Schmutz, S. M., and Dreger, D. L. (2012). Interaction ofMC1RandPMELalleles on solid coat colors in Highland cattle. Anim. Genet. 44, 9–13. doi: 10.1111/j.1365-2052.2012.02361.x
Sempéré, G., Moazami-Goudarzi, K., Eggen, A., Laloë, D., Gautier, M., and Flori, L. (2015). WIDDE: a Web-Interfaced next generation Database for genetic Diversity Exploration, with a first application in cattle. BMC Genom 16:940.
Seya, T., Oshiumi, H., Sasai, M., Akazawa, T., and Matsumoto, M. (2005). TICAM-1 and TICAM-2: toll-like receptor adapters that participate in induction of type 1 interferons. Int. J. Biochem. Cell Biol. 37, 524–529. doi: 10.1016/j.biocel.2004.07.018
Sharma, R., Kumar, B., Arora, R., Ahlawat, S., Mishra, A., and Tantia, M. (2016). Genetic diversity estimates point to immediate efforts for conserving the endangered Tibetan sheep of India. Meta Gene 8, 14–20. doi: 10.1016/j.mgene.2016.01.002
Shathele, M. (2009). Weather Effect on Bacterial Mastitis in Dairy Cows. Int. J. Dairy Sci. 4, 57–66. doi: 10.3923/ijds.2009.57.66
Shin, D., Lee, C., Park, K.-D., Kim, H., and Cho, K.-H. (2015). Genome-association analysis of Korean Holstein milk traits using genomic estimated breeding value. Asian Austral. J. Anim. Sci. 30, 309–319. doi: 10.5713/ajas.15.0608
Silva, R. M. D. O., Stafuzza, N. B., Fragomeni, B. D. O., Camargo, G. M. F. D., Ceacero, T. M., Joslaine Noely, et al. (2017). Genome-Wide Association Study for Carcass Traits in an Experimental Nelore Cattle Population. PLoS One 12:24. doi: 10.1371/journal.pone.0169860
Slimen, I. B., Najar, T., Ghram, A., and Abdrrabba, M. (2015). Heat stress effects on livestock: molecular, cellular and metabolic aspects, a review. J. Anim. Physiol. Anim. Nutr. 100, 401–412. doi: 10.1111/jpn.12379
Solé, M., Gori, A.-S., Faux, P., Bertrand, A., Farnir, F., Gautier, M., et al. (2017). Age-based partitioning of individual genomic inbreeding levels in Belgian Blue cattle. Genet. Select. Evolut. 49:92. doi: 10.1186/s12711-017-0370-x
Song, G., Satterfield, M. C., Kim, J., Bazer, F. W., and Spencer, T. E. (2009). Progesterone and interferon tau regulate leukemia inhibitory factor receptor and IL6ST in the ovine uterus during early pregnancy. Reproduction 137, 553–565. doi: 10.1530/rep-08-0437
Srikanth, K., Kwon, A., Lee, E., and Chung, H. (2016). Characterization of genes and pathways that respond to heat stress in Holstein calves through transcriptome analysis. Cell Stress Chaper. 22, 29–42. doi: 10.1007/s12192-016-0739-8
Stainton, J. J., Haley, C. S., Charlesworth, B., Kranis, A., Watson, K., and Wiener, P. (2014). Detecting signatures of selection in nine distinct lines of broiler chickens. Anim. Genet. 46, 37–49. doi: 10.1111/age.12252
Sugimoto, M., Sasaki, S., Gotoh, Y., Nakamura, Y., Aoyagi, Y., Kawahara, T., et al. (2013). Genetic variants related to gap junctions and hormone secretion influence conception rates in cows. Proc. Natl. Acad. Sci. 110, 19495–19500. doi: 10.1073/pnas.1309307110
Sved, J. (1971). Linkage disequilibrium and homozygosity of chromosome segments in finite populations. Theoret. Popul. Biol. 2, 125–141. doi: 10.1016/0040-5809(71)90011-6
Syed-Shabthar, S. M. F., Rosli, M. K. A., Mohd-Zin, N. A. A., Romaino, S. M. N., Fazly-Ann, Z. A., Mahani, M. C., et al. (2013). The molecular phylogenetic signature of Bali cattle revealed by maternal and paternal markers. Mole. Biol. Rep. 40, 5165–5176.
Takahashi, M. (2011). Heat stress on reproductive function and fertility in mammals. Reproduct. Med. Biol. 11, 37–47. doi: 10.1007/s12522-011-0105-6
Taye, M., Lee, W., Caetano-Anolles, K., Dessie, T., Hanotte, O., Mwai, O. A., et al. (2017). Whole genome detection of signature of positive selection in African cattle reveals selection for thermotolerance. Anim. Sci. J. 88, 1889–1901. doi: 10.1111/asj.12851
Thamo, T., Addai, D., Pannell, D. J., Robertson, M. J., Thomas, D. T., and Young, J. M. (2017). Climate change impacts and farm-level adaptation: Economic analysis of a mixed cropping–livestock system. Agricult. Syst. 150, 99–108. doi: 10.1016/j.agsy.2016.10.013
Therrien, A., Manjunath, P., and Lafleur, M. (2013). Chemical and physical requirements for lipid extraction by bovine binder of sperm BSP1. Biochim. Biophys. Acta 1828, 543–551. doi: 10.1016/j.bbamem.2012.08.020
Thumbi, S. M., Bronsvoort, B. M. D. C., Poole, E. J., Kiara, H., Toye, P. G., Mbole-Kariuki, M. N., et al. (2014). Parasite Co-Infections and Their Impact on Survival of Indigenous Cattle. PLoS One 9:e76324. doi: 10.1371/journal.pone.0076324
Upadhyay, M., Eriksson, S., Mikko, S., Strandberg, E., Stålhammar, H., Groenen, M. A. M., et al. (2019). Genomic relatedness and diversity of Swedish native cattle breeds. Genet. Select. Evolut. 51:96. doi: 10.1186/s12711-019-0496-0
Walugembe, M., Bertolini, F., Dematawewa, C. M. B., Reis, M. P., Elbeltagy, A. R., Schmidt, C. J., et al. (2019). Detection of Selection Signatures Among Brazilian, Sri Lankan, and Egyptian Chicken Populations Under Different Environmental Conditions. Front. Genet. 9:737. doi: 10.3389/fgene.2018.00737
Wang, K., Hu, Q., Ma, H., Wang, L., Yang, Y., Luo, W., et al. (2014). Genome-wide variation within and between wild and domestic yak. Mole. Ecol. Resour. 14, 794–801. doi: 10.1111/1755-0998.12226
Wang, Y., Cai, H., Luo, X., Ai, Y., Jiang, M., and Wen, Y. (2020). Insight into unique somitogenesis of yak (Bos grunniens) with one additional thoracic vertebra. BMC Genom. 21:201. doi: 10.1186/s12864-020-6598-9
Weldenegodguad, M., Popov, R., Pokharel, K., Ammosov, I., Ming, Y., Ivanova, Z., et al. (2019). Whole-Genome Sequencing of Three Native Cattle Breeds Originating From the Northernmost Cattle Farming Regions. Front. Genet. 9:728. doi: 10.3389/fgene.2018.00728
Welderufael, B., Janss, L., Koning, D. D., Sørensen, L., Løvendahl, P., and Fikse, W. (2017). Bivariate threshold models for genetic evaluation of susceptibility to and ability to recover from mastitis in Danish Holstein cows. J. Dairy Sci. 100, 4706–4720. doi: 10.3168/jds.2016-11894
Weller, J., Ezra, E., and Ron, M. (2017). Invited review: A perspective on the future of genomic selection in dairy cattle. J. Dairy Sci. 100, 8633–8644. doi: 10.3168/jds.2017-12879
West, J. (1994). Interactions of Energy and Bovine Somatotropin with Heat Stress. J. Dairy Sci. 77, 2091–2102. doi: 10.3168/jds.s0022-0302(94)77152-6
Wiener, G. (2013). The yak, an essential element of the high altitude regions of Central Asia. Études mongoles et sibériennes. Centrasiatiques et Tibétaines 43–44, 2–10. doi: 10.4000/emscat.2038
Wiener, G., Jianlin, H., Long, R., and Li, C. (2003). The yak. Bangkok, Thailand: Regional Office for Asia and the Pacific. Thailand: Food and Agricultural Organization of the United Nations.
Willing, E.-M., Dreyer, C., and Oosterhout, C. V. (2012). Estimates of Genetic Differentiation Measured by FST Do Not Necessarily Require Large Sample Sizes When Using Many SNP Markers. PLoS One 7:e42649. doi: 10.1371/journal.pone.0042649
Wright, S. (1922). Coefficients of Inbreeding and Relationship. Am. Natur. 56, 330–338. doi: 10.1086/279872
Wu, J. J., Song, L. J., Wu, F. J., Liang, X. W., Yang, B. Z., Wathes, D. C., et al. (2012). Investigation of transferability of BovineSNP50 BeadChip from cattle to water buffalo for genome wide association study. Mole. Biol. Rep. 40, 743–750. doi: 10.1007/s11033-012-1932-1
Xu, L., Zhao, G., Yang, L., Zhu, B., Chen, Y., Zhang, L., et al. (2019a). Genomic Patterns of Homozygosity in Chinese Local Cattle. Scient. Rep. 9:16977. doi: 10.1038/s41598-019-53274-3
Xu, L., Zhu, B., Wang, Z., Xu, L., Liu, Y., Chen, Y., et al. (2019b). Evaluation of Linkage Disequilibrium, Effective Population Size and Haplotype Block Structure in Chinese Cattle. Animals 9:83. doi: 10.3390/ani9030083
Yang, L., Xu, L., Zhu, B., Niu, H., Zhang, W., Miao, J., et al. (2017). Genome-wide analysis reveals differential selection involved with copy number variation in diverse Chinese Cattle. Scient. Rep. 7:1. doi: 10.1038/s41598-017-14768-0
Yudin, N. S., Larkin, D. M., and Ignatieva, E. V. (2017). A compendium and functional characterization of mammalian genes involved in adaptation to Arctic or Antarctic environments. BMC Genet. 18:111. doi: 10.1186/s12863-017-0580-9
Yurchenko, A. A., Daetwyler, H. D., Yudin, N., Schnabel, R. D., Jagt, C. J. V., Soloshenko, V., et al. (2018). Scans for signatures of selection in Russian cattle breed genomes reveal new candidate genes for environmental adaptation and acclimation. Scient. Rep. 8:12984. doi: 10.1038/s41598-018-31304-w
Zenger, K. R., Khatkar, M. S., Cavanagh, J. A. L., Hawken, R. J., and Raadsma, H. W. (2007). Genome-wide genetic diversity of Holstein Friesian cattle reveals new insights into Australian and global population variability, including impact of selection. Anim. Genet. 38, 7–14. doi: 10.1111/j.1365-2052.2006.01543.x
Zhang, B., Zhang, Y., Zhang, L., Wang, J., Li, Z., and Chen, H. (2011). Allelic Polymorphism Detected in the Bovine FTO Gene. Mole. Biotechnol. 49, 257–262. doi: 10.1007/s12033-011-9400-z
Zhang, Q., Calus, M. P., Guldbrandtsen, B., Lund, M. S., and Sahana, G. (2015). Estimation of inbreeding using pedigree, 50k SNP chip genotypes and full sequence data in three cattle breeds. BMC Genet. 16:1. doi: 10.1186/s12863-015-0227-7
Zhang, W., Gao, X., Zhang, Y., Zhao, Y., Zhang, J., Jia, Y., et al. (2018). Genome-wide assessment of genetic diversity and population structure insights into admixture and introgression in Chinese indigenous cattle. BMC Genet. 19:114. doi: 10.1186/s12863-018-0705-9
Keywords: cold stress, climate resilience, genetic resources, heat stress, heat tolerance, selective sweep
Citation: Freitas PHF, Wang Y, Yan P, Oliveira HR, Schenkel FS, Zhang Y, Xu Q and Brito LF (2021) Genetic Diversity and Signatures of Selection for Thermal Stress in Cattle and Other Two Bos Species Adapted to Divergent Climatic Conditions. Front. Genet. 12:604823. doi: 10.3389/fgene.2021.604823
Received: 10 September 2020; Accepted: 15 January 2021;
Published: 03 February 2021.
Edited by:
Vinicius De França Carvalho Fonsêca, Federal University of Paraíba, BrazilReviewed by:
Lingyang Xu, Chinese Academy of Agricultural Sciences, ChinaMarkus Neuditschko, Agroscope, Switzerland
Copyright © 2021 Freitas, Wang, Yan, Oliveira, Schenkel, Zhang, Xu and Brito. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yi Zhang, zhangyi@cau.edu.cn; Qing Xu, qingxu@bjtu.edu.cn; Luiz F. Brito, britol@purdue.edu