Selection Signatures in South African Nguni and Bonsmara Cattle Populations Reveal Genes Relating to Environmental Adaptation

Climate change is a major influencing factor in beef production. The greenhouse gases produced from livestock production systems contribute to the overall greenhouse gas emissions. The aim of this study was to identify selection signatures within and between Nguni and Bonsmara cattle in relation to production and adaptation. For this purpose, genomic 150 K single nucleotide polymorphism data from Nguni (n = 231) and Bonsmara (n = 252) cattle in South Africa were used. Extended haplotype homozygosity (EHH) based analysis was executed within each population using integrated haplotype score (iHS). The R package rehh was used for detecting selection signatures across the two populations with cross population EHH (XP-EHH). Total of 121 regions of selection signatures were detected (p < 0.0001) in the Bonsmara and Nguni populations. Several genes relating to DNA methylation, heat stress, feed efficiency and nitrogen metabolism were detected within and between each population. These regions also included QTLs associated with residual feed intake, residual gain, carcass weight, stature and body weight in the Bonsmara, while QTLs associated with conception rate, shear force, tenderness score, juiciness, temperament, heat tolerance, feed efficiency and age at puberty were identified in Nguni. Based on the results of the study it is recommended that the Nguni and Bonsmara be utilized in crossbreeding programs as they have beneficial traits that may allow them to perform better in the presence of climate change. Results of this study coincide with Nguni and Bonsmara breed characteristics and performance, and furthermore support informative crossbreeding programs to enhance livestock productivity in South Africa.

Introduction

Among different factors influencing production in the agricultural sector, the impact of climate change could be one of the most challenging (Calzadilla et al., 2014). The Intergovernmental Panel on Climate Change (IPCC) is predicting a 1.5°C increase above pre-industrial levels between 2030 and 2052 if warming continues to increase at the current rate (Intergovernmental Panel on Climate Change, 2018). Increasing occurrences of heat stress, drought and flooding events are likely, and these will have adverse effects on crop and livestock productivity over and above the impacts due to fluctuations in mean variable temperature (Intergovernmental Panel on Climate Change, 2019). For example, during the 2015 drought season 40,000 head of cattle died due to extreme heat waves experienced in the Kwa-Zulu Natal province (Willemse and Mare, 2016). The impact of climate change may be limited by using smaller framed animals, as well as using locally adapted breeds in commercial production systems (Scholtz et al., 2014). Locally adapted breeds are known for their ability to withstand changing climate conditions as well as resistance to parasites and diseases (Chan et al., 2010). The Nguni, a distinct Sanga type, is an indigenous cattle breed that displays adaptation to changing environments and is efficient in coping with various stressors such as parasites, heat spells, and lack of water availability (Mwai et al., 2015; Nyamushamba et al., 2017). Furthermore, Nguni cattle have lower metabolic energy requirements in heat challenged environments, as they regulate internal metabolic heat production with lowered feed intake (Sejian et al., 2018; Jiyana, 2019). Nguni cattle have smaller frames, with various horn shapes and hide patterns and are used for beef production. The Bonsmara is a local composite breed derived from crossing 5/8 Afrikaner (indigenous Sanga or Bos taurus africanus) and 3/16 Shorthorn and 3/16 Hereford (Bos taurus) (Strydom et al., 2001). Nguni cattle have smaller frames, with various horn shapes and hide patterns and are used for beef production. The Bonsmara is a local composite breed derived from crossing 5/8 Afrikaner (indigenous Sanga or Bos taurus africanus) and 3/16 Shorthorn and 3/16 Hereford (Bos taurus) (Strydom et al., 2001; Corbet et al., 2006). Bonsmara cattle are medium to large framed, smooth coated with tolerance to heat and ticks (Muchenje et al., 2008). The Bonsmara is known for beef production, high fertility rates and efficient performance in feedlot systems (Esterhuizen et al., 2008; Asizua et al., 2017).

Natural and artificial selection leads to changes in certain regions of the genome resulting in selection signatures that can reveal genes associated with the selected traits (Xu et al., 2015). Signatures of selection are the changes in allele frequencies, linkage disequilibrium and genetic variation over time and generate unique patterns of specific regions of the genome (Aliloo et al., 2020). Long-term selection pressures can be investigated through the detection of signatures of selection. With the utilization of advanced genetic technology, a more in depth view of the genetic status of each breed is possible. There are two main categories for identifying selection signatures. Firstly, there are signatures that can be detected within a population and secondly, signatures than can be identified across two populations (Voight et al., 2006; Tang et al., 2007). To detect selection signatures, an estimator called extended haplotype homozygosity (EHH) to search for genetic footprints of positive selection was developed by Sabeti et al. (2007). Several studies have employed this approach in composite cattle (Yurchenko et al., 2018; Singh et al., 2020) and indigenous breeds (Bahbahani et al., 2018; Mastrangelo et al., 2020). For instance, Saravanan et al. (2021) performed a scan of selection signatures in Indian taurine and indicine breeds and identified several candidate genes relating to adaptation such as heat stress (HSPA1B) and disease susceptibility (DEFB4). While, Yurchenko et al. (2018) identified signatures of selection in nine Russian cattle breeds that related to several genes such as feed intake (TMEM68, PLAG1), reproduction (HFM1, CSF2) and acclimation (RETREG1, AQP5). Interestingly, the study by Ben-Jemaa et al. (2020) identified several genes to be under positive selection in North African cattle populations, such as adaptive immune response (MAP3K3, IRF8), blood pressure regulation and heart contraction related genes were theorised to be involved with adaptation to extreme temperatures, which may explain the presence of hypoxia (BCL2, CBFA2T3) and heat (HSPH1, MVD) related genes in these cattle. Collectively, these studies show that regions under selection may harbour genes relating to environmental adaption in cattle.

Makina et al. (2015) studied the signatures of selection in six South African Sanga cattle breeds, which includes the Bonsmara and Nguni breeds, using the Illumina BovineSNP50 Bead chip. With the availability of high density SNP data, there is an opportunity to evaluate the genetic structure and selection pressures present in South African cattle, particularly in the Nguni and Bonsmara breeds, as this will serve as a baseline for crossbreeding studies of the reciprocal crosses of these populations. An observation made during hot and dry year (2015/16), when temperatures increased above normal, the Nguni-sired calves performed better than the Bonsmara-sired calves (Pyoos, 2018). Additionally, the identification of genes targeted by selection can be helpful in finding genes related to environmental traits (e.g., genes relating to adaptation) that are difficult to find experimentally. Thus, the aim of this study was to detect signatures of selection in Nguni and Bonsmara cattle using the GGP Bovine 150 K SNP data.

Materials and Methods

Ethical Approval

Ethical approval was obtained from the Agricultural Research Council, Animal Production Ethics Committee (APIEC21/06). In addition, approval for the use of the data was obtained from the Bonsmara SA and Nguni Cattle Breeders’ Societies respectively, and the Beef Genomics Programme.

SNP Data and Quality Control

SNP data generated using the GeneSeek Genomic Profiler 150 K (GGP Bovine 150 K) bead chip of two South African cattle breeds, Nguni (n = 258) and Bonsmara (n = 248) was used in this study. Individual datasets were merged using PLINK v1.9 (Purcell et al., 2007). Additionally, only autosomal SNPs that were common across all datasets were retained while SNPs present within mitochondrial and sex chromosomes were excluded to avoid any bias. Quality control was performed on the merged data with the parameters involving removal of SNPs with less than 95% call rate; SNPs with less than 2% minor allele frequencies (MAF); samples with more than 10% missing genotypes and SNPs that failed the Hardy Weinberg Equilibrium test (p < 0.001). Identity-by-descent (IBD) relatedness between individuals within each population was tested using the “genome” parameter in PLINK v 1.9 (Purcell et al., 2007). For this test, PI_HAT (Proportion of IBD), i.e. P(IBD = 2) + 0.5*P(IBD = 1), where P denotes the probability)—a variable that reveals the extended haplotypes shared between distantly related individuals was determined. Closely related individuals with a PI_HAT value greater than 0.25 were removed.

Genetic Diversity

The expected heterozygosity (He) and observed heterozygosity (Ho) were computed using the “hardy” parameter in PLINK v1.9 (Purcell et al., 2007) to assess the level of genetic diversity. Briefly, the Ho estimates for each population was determined from observed genotype frequencies obtained as follows: (N–O)/N, where N is the number of non-missing genotypes and O is the number of observed homozygous genotypes for a given individual.

Population Structure

Principal Component Analysis (PCA) was performed on the Nguni and Bonsmara cattle populations to validate breed separation in the merged (two breeds) datasets. More specifically, the relationship between the first two principal components was examined to show relationship between the animals. The Eigen values were generated using PLINK v1.9, thereafter a PCA plot was generated in R v 4.0.3 (R Core Team, 2019) using the package Tidyverse (Wickham et al., 2019). The estimation of population genetic structure was performed on both populations with the R package LEA (Frichot and François, 2015). Binary files were generated in PLINK v1.9 (Purcell et al., 2007). Thereafter, “.ped” and “.map” files were transformed into the “.geno” format required as input by sNMF using the “ped2geno” function. Individual ancestry coefficients were determined based on the calculation of cross-entropy criterion using the sparse non-negative matrix factorization (snmf) function from the LEA package with K (i.e., number of hypothetical ancestors) values ranging from 1 to 10 with five iterations each. The optimal number of ancestors had the lowest cross-entropy criterion. Finally, the individual ancestry coefficients were plotted using the plot function in R v 4.0.2 (R Core Team, 2019).

Signatures of Selection

In this study, signatures of selection were identified within and between the two populations. The within population selection signatures was used to show the regions under selection in each of the cattle populations. While the across population selection signatures was used to show the regions under positive selection in both populations. To prepare input files for signatures of selection testing, binary files were created for each breed per chromosome using PLINK v1.9 (Purcell et al., 2007). The pedigree and map files generated were converted to fastPHASE formats using the fcGENE software (Roshyara and Scholz, 2014). The binary files were converted using the “oformat./fastPhase” function in fcGENE, producing the output files in the “.inp” format. Haplotypes for the iHS and XP-EHH analyses were derived with fastPHASE (Scheet and Stephens, 2006) using the output files generated by fcGENE with 10 starts (T10) and 25 iterations (C25) of the expectation-maximization (EM) algorithm. The fastPhase output files were used to detect signatures of selection using the R v 4.0.3 (R Core Team, 2019) and package rehh (Gautier et al., 2017).

The within population signatures of selection in the Nguni and Bonsmara population were determined. The iHS scores were calculated for both the Nguni and Bonsmara using the “ihs2ihs” function in rehh package with the minimum MAF at 0.02. This test detects evidence of recent positive selection at a locus based on the differential levels of LD surrounding a positively selected (derived) allele compared to the background (ancestral) allele at the same position within a population. The iHS candidate regions were identified at p < 0.0001, with a window size of 1 Mb allowing for a 10 Kb overlap and the region must contain a minimum of 2 markers. The iHS scores for both the Nguni and Bonsmara were plotted using the “manhattanplot” function from package rehh (Gautier et al., 2017). For visualisation and comparison of selection signals, | iHS| scores were transformed into - log 10 [2Φ - |iHS|] where Φ represents the Gaussian cumulative distribution functionand a value of 4 is equivalent to p < 0.0001 (Gautier and Naves, 2011). Cross-population extended haplotype heterozygosity (XP-EHH), statistics was determined to compare EHH profiles between the Nguni and Bonsmara populations. This statistic tests whether the genome site is homozygous in the one population (e.g., Nguni) but polymorphic in the other population (e.g., Bonsmara) through the evaluation of the two populations on one core SNP (Maiorano et al., 2018). The XP-EHH was determined between the Nguni and Bonsmara population using the “ies2xpehh” function in rehh (Sabeti et al., 2007). The candidate regions and plots were identified and visualized with the same parameters as for the iHS test.

Genes contained within the identified candidate genomic regions with signatures of selection have been annotated using the Bovine UMD3.1 reference genome (https://bovinegenome.elsiklab.missouri.edu/). The Bovine QTL database (Hu et al., 2013; http://www.animalgenome.org/cgi-bin/QTLdb/BT/search) was used to identify the overlap with previously published bovine QTLs within the candidate regions. BioMart (Zhang et al., 2011), a program in Ensembl was used to identify genes found within the candidate regions (https://www.ensembl.org/biomart/martview/717be2236eefa78b5aecbc78a6aee7ea).

Results

SNP Data and Quality Control

The relatedness test showed no individuals with PI_HAT value greater than 0.25, hence no individuals were removed based on this test. Due to low call rates, 27 Nguni and 32 Bonsmara animals were removed from the dataset. Furthermore, 40,224 and 28,409 SNPs were excluded based on failure to meet the requirements of the missing genotypes, MAF<0.02, the significance of deviation from Hardy Weinberg Equilibrium (p < 0.001) in the Nguni and Bonsmara, respectively. Finally, 231 Nguni with 100,952 SNPs and 252 Bonsmara animals with 110,556 SNPs were used for further analysis. A total of 94,147 SNPs were overlapping in the two populations while 6,805 SNPs were unique to the Nguni and 16,419 SNPs were unique to the Bonsmara population.

Genetic Diversity

The mean He and Ho were higher in the Nguni population compared to the Bonsmara population. The mean He coefficient in the Nguni was 0.363 ± 0.126 while the mean Ho was 0.362 ± 0.130, which was slightly higher compared to the Bonsmara. However, in the Bonsmara population, a difference of 0.003 was observed between the mean He = 0.361 ± 0.132 and the mean Ho = 0.358 ± 0.133.

Population Structure

The PCA analysis divided the study population into two distinct clusters (Figure 1). The first principal component (PC1) accounted for 31.4% of the total variation and separated the Nguni population from the Bonsmara population. The second principal component (PC2) accounted for 7.59% of the total variation. The clustering was further confirmed by the ancestry matrix of the two populations (Figure 2). For the optimal number of ancestors, the lowest cross-entropy value of 0.851 was chosen with the number of ancestral populations (K) of 2. The ancestry matrix shows that the Nguni population has a different genetic structure when compared to the Bonsmara population.

FIGURE 1

FIGURE 1. Principal component analysis (PCA) plot showing clustering of Nguni and Bonsmara cattle. Each point represents an individual animal coloured as per breed.

FIGURE 2

FIGURE 2. Ancestry matrix of Nguni and Bonsmara populations. Shown in the figure is the ancestry bar plot using sparse nonnegative matrix factorization (sNMF).

Genomic Regions of Signatures of Selection

Positive signatures of selection were identified within and across the Nguni and Bonsmara populations. The iHS analysis for the Nguni revealed 118 regions of selection across all chromosomes except on BTA 8, 20, 21, 22, and 28 (Figure 3). The highest iHS score of 6.047 was located on BTA 12 (p < 0.0001). In contrast to the Nguni, only three significant regions were found for the Bonsmara. The iHS analysis for Bonsmara only identified regions on BTA 5, 12 and 14 (p < 0.0001, Figure 4). For the across population analysis, three regions of selection were identified across the two populations located on BTA 3, 12 and 14 (p < 0.0001, Figure 5). The highest XP-EHH score of 7.243 was detected on BTA14. The genes identified within the candidate regions are summarised in Table 1 and Supplementary Table S1. Additionally, all the putative genes identified across all candidate are summarized in Supplementary Table S2.

FIGURE 3

FIGURE 3. Genome-wide distribution of the integrated haplotype score (iHS) in the Nguni Population. Plot depicting the candidate regions of Nguni population with p < 0.0001 indicated with a dotted line The Φ represents the Gaussian cumulative distribution function for iHS Nguni scores.

FIGURE 4

FIGURE 4. Genome-wide distribution of the integrated haplotype score (iHS) in the Bonsmara Population. Plot depicting the candidate regions of Bonsmara population with p < 0.0001 indicated with a dotted line. The Φ represents the Gaussian cumulative distribution function for iHS Bonsmara scores.

FIGURE 5

FIGURE 5. Manhattan plots showing the candidate regions of cross-population extended haplotype homozygosity (XP-EHH) analysis between Bonsmara and Nguni cattle populations with the p < 0.0001 threshold indicated with a dotted line. The Φ represents the Gaussian cumulative distribution function for XP-EHH scores.

TABLE 1

TABLE 1. Genomic candidate regions detected by integrated haplotype score (iHS) and cross-population extended haplotype homozygosity (XP-EHH) analysis in the Nguni and Bonsmara population, p < 0.0001.

Overlapping regions were further interrogated through the comparison of all the regions of selection identified in this study. Based on the criteria of Voight et al. (2006), the overlapping regions of selection are defined as those regions located above the cut-off threshold and in the same chromosomal location. Two significant selection regions were detected by iHS Bonsmara, iHS Nguni and XP-EHH tests. The two regions were detected on BTA12 (32.2–34.1 Mb) and BTA14 (23.5–25.4 Mb) (Table 2). Furthermore, one overlapping region was detected by the iHS analysis between the Nguni and Bonsmara (Table 3).

TABLE 2

TABLE 2. Genes and QTLs found in the overlapping regions detected by integrated haplotype score (iHS) and cross-population extended haplotype homozygosity (XP-EHH) in the analysis Nguni and Bonsmara cattle (p < 0.0001).

TABLE 3

TABLE 3. One overlapping region detected by comparing integrated haplotype score (iHS) analysis in Nguni and iHS analysis Bonsmara populations (p < 0.0001).

Discussion

Genetic exploration of beef cattle breeds can provide vital information that may contribute to sustainable meat production systems such as crossbreeding. The main goal was to unravel the genetic variability and determine the signatures of selection in Nguni, an indigenous cattle breed and the Bonsmara, a composite breed using bovine 150 K SNP data. The results of this study showed that both Nguni and Bonsmara populations have a high level of genetic diversity. More specifically, the mean expected heterozygosity values were comparable in Nguni (He = 0.362) and Bonsmara (He = 0.361) populations (Table 1). In addition, the level of genetic diversity in Bonsmara populations was similar to those reported by Bosman et al. (2017) with He = 0.365. However, the Nguni population had higher level of genetic diversity than previously reported by Makina et al. (2014) with He = 0.28, although this may be due to the low number of samples in the previous study. Furthermore, PCA and ancestry matrix analyses agreed with previous findings of Makina et al. (2016), where the Nguni cluster is distinct from the Bonsmara cluster (Figures 1, 2). These PCA results are also similar to those found by Madilindi et al. (2020), where the South African Nguni population clustered away from the Bonsmara reference population when using microsatellite data. Overall, these findings display that the Nguni population is significantly different to the Bonsmara, and each population is genetically variable.

Candidate regions that show preferential selection in the genome of Nguni and Bonsmara cattle were discovered using the two statistical methods (iHS and XP-EHH), successfully applied in multiple studies to identify signatures of selection in cattle (Bahbahani et al., 2018; Cheruiyot et al., 2018; Alshawi et al., 2019; Aliloo et al., 2020). While the iHS analysis detected 118 regions of selection within the Nguni population, only three regions of selection signatures were found within the Bonsmara population. This may be due to Nguni experiencing selection pressures over longer period of time than the Bonsmara cattle (Makina et al., 2015). The XP-EHH analysis across the Nguni and Bonsmara population identified 3 regions under selection (BTA 3, 12 and 14) that were previously found to be under selection. For example, the study by Ben-Jemaa et al. (2020) detected the same region (BTA 12:32.2-34.1 Mb) using the XP-EHH tests in Indicine (Gir and Nelore) vs. North African cattle analysis. Another study by Bahbahani et al. (2015) identified the same region on BTA 12: 32.2-34.1 Mb comparing East African Shorthorn Zebu (EASZ) vs. European Taurine (Holstein-Friesian and Jersey) and comparing EASZ vs. African Taurine (N’Dama). Cheruiyot et al. (2018) also identified the region on BTA 14: 23.49-23.86 Mb through XP-EHH analysis in Tanzanian crossbred cattle vs. EASZ. The genes detected in this region include, PLAG1, CHCHD&, TOX, and TMEM68 involved in pleiotropic traits associated with milk and growth in cattle. Reported studies displayed regions of signatures of selection present on BTA 12 and 14 and coincide with the findings of this study.

Climate change may cause increased thermal stress experienced by livestock species, resulting in reduced meat production and decreased fertility that will have a negative impact on food security (Sejian et al., 2018). In this study, several genes related to heat stress were identified in the genomic regions displaying positive signatures of selection in both populations. This is important, as heat stress is the physiological imbalance between the heat produced within the body and its expenditure, overall negatively impacting growth, milk production and reproduction (Belhadj Slimen et al., 2016). Genes relating to heat shock included HSP90 which is one of the members of the heat shock protein family that plays a key role in influencing heat stress in animals. Charoensook et al. (2012) identified polymorphisms in the bovine HSP90AB1 associated with heat tolerance in Thai indigenous cattle, which may support the presence of HSP90AB1 (BTA23: 17.75-17.75 Mb) in the Nguni cattle, and not being detected in the Bonsmara cattle. Another DnaJ heat shock protein family (Hsp40) member B13 (DNAJB13) gene located on BTA 15: 53.39-53.40 Mb was found in the Nguni population. However, this gene was shown to be expressed in both Nguni and Bonsmara populations during low energy diet treatment (Linde, 2018). The heat shock protein family H (Hsp110) member 1 (HSPH1) detected on BTA12: 29.79- 29.81 Mb has shown to be involved in Nguni and Bonsmara body temperature regulation in the presence of climatic stress (Howard et al., 2014). Heat shock protein family B (small) member 9 (HSPB9) located on BTA 19 was previously detected in Nguni cattle and showed to be associated to tropical adaptation in Zebu cattle (Makina et al., 2015). This region (BTA 19: 42.25-42.26 Mb) is under selection in the current Nguni population. The expression of additional heat shock proteins in the Nguni could play a role in its ability to tolerate high temperatures better than the Bonsmara cattle.

DNA methylation is one of the necessary processes involved in epigenetic gene regulation (Tian et al., 2017). Hence, genes relating to DNA methylation were also located in the regions displaying signatures of selection. The de novo and maintenance methylation reactions carried out by DNA methyl transferases have been studied extensively (Chen and Li, 2006; Yakovlev, 2018). Despite the stable and heritable features of DNA methylation patterns, genome-wide DNA demethylation occurs both in developing germ cells and in fertilized oocytes. Mbd3l1 expressed predominantly in testis and found specifically in round spermatids was identified in Nguni cattle (BTA17: 44.25-44.26 Mb). Expression of Mbd3l1 is specific to haploid male germ cells. This protein could be incorporated into sperm chromatin, and after fertilization could participate in demethylation of the paternal genome (Jin et al., 2008). Methyltransferase like 21A (METTL21A) located on BTA2, is responsible for the methylation of HSPA1-K561, a heat shock protein that was identified in the Nguni population. This interaction displays the vital role of DNA methylation in response to heat stress (Jakobsson et al., 2016). Protein arginine methyl transferase 8 (PRMT8) was shown to play a role in lactation persistency in Canadian Holstein cattle and was found in both Nguni and Bonsmara cattle, located on BTA 5 (Do et al., 2017). This further supports the detection of Lactation persistency (QTL: 125,220) located on BTA 5:105.6-107.5 Mb. Methyltransferase like 21C (METTL21C), located on BTA12 was shown to be differently expressed in Bonsmara and Nguni breeds by Linde (2018). However, in this study, METTL21C was detected only in the Nguni cattle. Additionally, METTL21C showed to be an important modulator of protein degradation in skeletal muscle in mice by Wiederstein et al. (2018). The XP-EHH revealed BTA12: 35.79-35.81 Mb to harbour EEF1AKMT1 to play a role in DNA methylation in both populations. The presence of DNA methylation regulating genes in the Nguni and Bonsmara are important for its adaptation to thermally challenged environments.

Maternal effects, which result from variation in the environment provided by dams to their progeny, are recognized as being important for a variety of traits expressed at, and after birth in beef cattle (Campêlo et al., 2004; Zhang et al., 2021). In addition, the study by Linde (2018) showed that the impact of a low energy diet on the expression of DNA methyltransferase 3B (DNMT3B) increased expression in the Bonsmara compared to the Nguni. However, in this study, DNMT3B located on BTA 13: 62.14-62.17 Mb was only identified in the Nguni population. This could possibly imply that the DNMT3B is under strong selection in the Nguni population and not in the Bonsmara. An interesting finding by Pyoos (2018) showed the maternal effect of the Nguni on birth weight and pre-weaning growth was lower compared to the Bonsmara, possibly due to reduced feed intake on the part of the Nguni cow and this could be another reason for the detection of DNMT3B only in the Nguni population.

Animal growth is dependent on feed intake and metabolism of feed into useful energy for the animal (Dong et al., 2019). Furthermore, the need for animals to utilize feed efficiently is a necessity in order to reduce secondary by-products such as methane emissions that are harmful to the environment (Dong et al., 2019). In this study candidate regions that display signatures of selection, include QTLs associated with residual feed intake, conception rate, body weight gain, carcass weight, residual feed intake and heat tolerance. The two overlapping regions that were detected by all EHH analysis located on BTA 12 and 14, include QTLs associated with traits such as immunity, meat quality, disease, growth, and reproduction. For instance, the QTL for feed conversion ratio (FCR) (QTL: 213,190) located on BTA14 was common across both populations, implying that both breeds have the ability to convert fibrous diet into reasonable weight gain (Shike, 2013). Jiyana (2019) found that the Nguni need less feed than the Bonsmara to gain the same weight. This is in line with the ability of the Nguni breed to survive on pasture grazing for a long time while maintaining its weight (Nyamushamba et al., 2017). The impact of using Nguni cattle with lower FCR in crossbreeding programs may reduce the amount of feed needed to sustain them for longer periods. This observation is supported by the feed efficiency genes such as serpin family H member 1 (SERPINH1, BTA15: 54.73-54.74 Mb), plasminogen activator, urokinase receptor (PLAUR, BTA18:51.86-51.87 Mb), RELB proto-oncogene, NF-kB subunit (RELB, BTA18:52.69-52.72 Mb) detected to be under selection in the Nguni population which have previously been detected. SERPINH1, PLAUR and RELB have shown to be associated to feed conversion ratio, dry matter intake, average daily gain and residual feed intake in beef cattle (Sun et al., 2019). Additionally, the XP-EHH analysis revealed two genes, the family with sequence similarity 110 member B (FAM110B) and UBX domain protein 2B (UBXN2B) relating to feed efficiency to be found on BTA 14:24.36-24.62 Mb in both Nguni and Bonsmara. For sufficient feed utilization, nitrogen metabolism plays a key role in the animal’s utilization of nitrogen from feed. In this study several genes under selection relating to nitrogen metabolism were identified only in the Nguni population. Keratin 17(KRT17) located on BTA19:41.83-41.84 Mb plays a role in the regulation of urine urea nitrogen yield and urine volume, while during milk protein synthesis Apolipoprotein E (APOE) located on BTA18:52.63-52.64 Mb improved ribosome biosynthesis resulting in less nitrogen being excreted as urine urea as previously detected by Wang et al. (2021). Other genes include Alpha-2-macroglobulin like 1 (A2ML1 located on BTA5:101.04-101.09 Mb) and 3-hydroxy-3-methylglutaryl-CoA synthase 2 (HMGCS2 located on BTA3: 23.53-23.56 Mb) were found to increase urine urea nitrogen and also previously detected by Wang et al. (2021). Previous studies observed that Nguni cattle were more capable of maintaining their body weight during winter than other breeds. Moreover they had higher blood urea (N) and ruminal ammonia levels (Linington, 1990) than other breeds. Ultimately, based on the presence of these genes in the Nguni and Bonsmara cattle, it clearly shows that they may be involved with effective feed regulation while also reducing the impact of nitrogen released into the environment.

Conclusion

With the use of high-density genome-wide SNP genotype data it was possible to understand the genetic diversity, population structure and selection signatures present in the South African Nguni and Bonsmara cattle populations. Broadly, the findings from this study revealed that these two populations are genetically distinct and have different selection signatures influencing their population structure. Candidate regions under positive selection were associated with heat tolerance, feed efficiency, nitrogen metabolism and DNA methylation. The QTLs relating to meat quality, disease resistance, growth, reproduction and immunity were located in the candidate regions under selection in the Nguni and Bonsmara populations. Furthermore, the results also indicate that different ancestral backgrounds (indigenous vs. composite breed genotypes) are advantageous in different regions of the genome, thus, they may assist in designing breeding programs to enhance performance of both the Nguni and Bonsmara cattle, as well as using the two breeds in crossbreeding programs. The genes identified may be used in improving the current Nguni and Bonsmara cattle population through controlled breeding strategies and selection of specific traits. The findings in this study present an opportunity in identifying causative variants that confer adaptation of indigenous and composite cattle populations to the South African environment.

Data Availability Statement

The datasets presented in this article are not readily available because of a binding memorandum of agreement between the South African Stud Book (SA Studbook) and the Agricultural Research Council (ARC). Requests to access the datasets should be directed to SA Studbook.

Ethics Statement

Ethical approval was obtained from the Agricultural Research Council, Animal Production Ethics Committee (APIEC21/06). In addition, approval for the use of the data was obtained from the Bonsmara SA and Nguni Cattle Breeders’ Societies respectively, and the Beef Genomics Programme.

Author Contributions

BK, PS, MN, MS, and FN conceived the idea and developed the theory. BK, PS, and MN performed the computations and verified the analytical methods. All authors discussed the results and contributed to the final manuscript.

Funding

This work is based on research that was partly supported by the National Research Foundation of South Africa (NRF), under grant UID 135438. The Grant-holder acknowledges that opinions, findings and conclusions or recommendations expressed in any publication generated by the NRF supported research are that of the authors and that the NRF accepts no liability whatsoever in this regard.

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.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors would like to thank the Bonsmara SA Cattle Breeders’ Society, Nguni Cattle Breeders’ Society and the Beef Genomics Programme for the use of SNP genotype data.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.909012/full#supplementary-material

References

Aliloo, H., Mrode, R., Okeyo, A. M., and Gibson, J. P. (2020). Ancestral Haplotype Mapping for GWAS and Detection of Signatures of Selection in Admixed Dairy Cattle of Kenya. Front. Genet. 11, 544. doi:10.3389/fgene.2020.00544

PubMed Abstract | CrossRef Full Text | Google Scholar

Alshawi, A., Essa, A., Al-Bayatti, S., and Hanotte, O. (2019). Genome Analysis Reveals Genetic Admixture and Signature of Selection for Productivity and Environmental Traits in Iraqi Cattle. Front. Genet. 10, 609. doi:10.3389/fgene.2019.00609

PubMed Abstract | CrossRef Full Text | Google Scholar

Asizua, D., Mpairwe, D., Kabi, F., Mutetikka, D., Hvelplund, T., Weisbjerg, M. R., et al. (2017). Effects of Grazing and Feedlot Finishing Duration on the Performance of Three Beef Cattle Genotypes in Uganda. Livest. Sci. 199, 25–30. doi:10.1016/j.livsci.2017.03.006

CrossRef Full Text | Google Scholar

Bahbahani, H., Clifford, H., Wragg, D., Mbole-Kariuki, M. N., Van Tassell, C., Sonstegard, T., et al. (2015). Signatures of Positive Selection in East African Shorthorn Zebu: A Genome-wide Single Nucleotide Polymorphism Analysis. Sci. Rep. 5 (1), 1–13. doi:10.1038/srep11729

CrossRef Full Text | Google Scholar

Bahbahani, H., Afana, A., and Wragg, D. (2018). Genomic Signatures of Adaptive Introgression and Environmental Adaptation in the Sheko Cattle of Southwest Ethiopia. PloS One 13 (8), e0202479. doi:10.1371/journal.pone.0202479

PubMed Abstract | CrossRef Full Text | Google Scholar

Belhadj Slimen, I., Najar, T., Ghram, A., and Abdrrabba, M. (2016). Heat Stress Effects on Livestock: Molecular, Cellular and Metabolic Aspects, a Review. J. Anim. Physiol. Anim. Nutr. 100 (3), 401–412. doi:10.1111/jpn.12379

CrossRef Full Text | Google Scholar

Ben-Jemaa, S., Mastrangelo, S., Lee, S.-H., Lee, J. H., and Boussaha, M. (2020). Genome-wide Scan for Selection Signatures Reveals Novel Insights into the Adaptive Capacity in Local North African Cattle. Sci. Rep. 10 (1), 1–14. doi:10.1038/s41598-020-76576-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Bhuiyan, M. S., Lim, D., Park, M., Lee, S., Kim, Y., Gondro, C., et al. (2018). Functional partitioning of genomic variance and genome-wide association study for carcass traits in Korean Hanwoo cattle using imputed sequence level SNP data. Front. Genet. 9, 217.

PubMed Abstract | CrossRef Full Text | Google Scholar

Bosman, L., van Marle-Köster, E., Van der Westhuizen, R. R., Visser, C., and Berry, D. P. (2017). Short Communication: Population Structure of the South African Bonsmara Beef Breed Using High Density Single Nucleotide Polymorphism Genotypes. Livest. Sci. 197, 102–105. doi:10.1016/j.livsci.2017.01.012

CrossRef Full Text | Google Scholar

Calzadilla, A., Zhu, T., Rehdanz, K., Tol, R. S. J., and Ringler, C. (2014). Climate Change and Agriculture: Impacts and Adaptation Options in South Africa. Water Resour. Econ. 5, 24–48. doi:10.1016/j.wre.2014.03.001

CrossRef Full Text | Google Scholar

Campêlo, J. E. G., Lopes, P. S., Torres, R. D. A., Silva, L. O. C. D., Euclydes, R. F., Araújo, C. V. D., et al. (2004). Maternal Effects on the Genetic Evaluation of Tabapuã Beef Cattle. Genet. Mol. Biol. 27, 517–521. doi:10.1590/S1415-47572004000400009

CrossRef Full Text | Google Scholar

Chan, E. K. F., Nagaraj, S. H., and Reverter, A. (2010). The Evolution of Tropical Adaptation: Comparing Taurine and Zebu Cattle. Anim. Genet. 41 (5), 467–477. doi:10.1111/j.1365-2052.2010.02053.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Charoensook, R., Gatphayak, K., Sharifi, A. R., Chaisongkram, C., Brenig, B., and Knorr, C. (2012). Polymorphisms in the Bovine HSP90AB1 Gene are Associated with Heat Tolerance in Thai Indigenous Cattle. Trop. Anim. Health Prod. 44 (4), 921–928. doi:10.1007/s11250-011-9989-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, T., and Li, E. (2006). “Establishment and Maintenance of DNA Methylation Patterns in Mammals,” in DNA Methylation: Basic Mechanisms. Berlin, Heidelberg: Springer, 179–201. doi:10.1007/3-540-31390-7_6

CrossRef Full Text | Google Scholar

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. Front. Genet. 9, 607. doi:10.3389/fgene.2018.00607

PubMed Abstract | CrossRef Full Text | Google Scholar

Corbet, N. J., Shepherd, R. K., Burrow, H. M., Van der Westhuizen, J., Strydom, P. E., and Bosman, D. J. (2006). Evaluation of Bonsmara and Belmont Red Cattle Breeds in South Africa. 1. Productive Performance. Aust. J. Exp. Agric. 46 (2), 199–212. doi:10.1071/EA05223

CrossRef Full Text | Google Scholar

Do, D. N., Bissonnette, N., Lacasse, P., Miglior, F., Sargolzaei, M., Zhao, X., et al. (2017). Genome-wide Association Analysis and Pathways Enrichment for Lactation Persistency in Canadian Holstein Cattle. J. Dairy Sci. 100 (3), 1955–1970. doi:10.3168/jds.2016-11910

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, L., Li, B., and Diao, Q. (2019). Effects of Dietary Forage Proportion on Feed Intake, Growth Performance, Nutrient Digestibility, and Enteric Methane Emissions of Holstein Heifers at Various Growth Stages. Animals 9 (10), 725. doi:10.3390/ani9100725

PubMed Abstract | CrossRef Full Text | Google Scholar

Esterhuizen, J., Groenewald, I. B., Strydom, P. E., and Hugo, A. (2008). The Performance and Meat Quality of Bonsmara Steers Raised in a Feedlot, on Conventional Pastures or on Organic Pastures. South Afr. J. Animal Sci. 38 (4), 303–314. Available at: https://hdl.handle.net/10520/EJC94590.

Google Scholar

Frichot, E., and François, O. (2015). LEA: An R Package for Landscape and Ecological Association Studies. Methods Ecol. Evol. 6 (8), 925–929. doi:10.1111/2041-210X.12382

CrossRef Full Text | Google Scholar

Gautier, M., and Naves, M. (2011). Footprints of Selection in the Ancestral Admixture of a New World Creole Cattle Breed. Mol. Ecol. 20 (15), 3128–3143. doi:10.1111/j.1365-294x.2011.05163.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Gautier, M., Klassmann, A., and Vitalis, R. (2017). Rehh 2.0: A Reimplementation of the R Package Rehh to Detect Positive Selection from Haplotype Structure. Mol. Ecol. Resour. 17 (1), 78–90. doi:10.1111/1755-0998.12634

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, X., Liu, X., Xu, X., Wu, M., Zhang, X., Li, Q., et al. (2012). The Expression Levels of DNMT3a/3b and Their Relationship with Meat Quality in Beef Cattle. Mol. Biol. Rep. 39 (5), 5473–5479. doi:10.1007/s11033-011-1349-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamey, J. J., and Wilkins, M. R. (2018). Methylation of Elongation Factor 1A: Where, Who, and Why? Trends Biochem. Sci. 43 (3), 211–223. doi:10.1016/j.tibs.2018.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Howard, J. T., Kachman, S. D., Snelling, W. M., Pollak, E. J., Ciobanu, D. C., Kuehn, L. A., et al. (2014). Beef Cattle Body Temperature during Climatic Stress: A Genome-wide Association Study. Int. J. Biometeorol. 58 (7), 1665–1672. doi:10.1007/s00484-013-0773-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, Z.-L., Park, C. A., Wu, X.-L., and Reecy, J. M. (2013). Animal QTLdb: An Improved Database Tool for Livestock Animal QTL/association Data Dissemination in the Post-genome Era. Nucleic acids Res. 41 (D1), D871–D879. doi:10.1093/nar/gks1150

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, J., Hsu, Y. H., Mo, C., Abreu, E., Kiel, D. P., Bonewald, L. F., et al. (2014). METTL21C is a Potential Pleiotropic Gene for Osteoporosis and Sarcopenia Acting Through the Modulation of the NF‐κB Signaling Pathway. J. Bone Miner. Res. 29 (7), 1531–1540.

PubMed Abstract | CrossRef Full Text | Google Scholar

Intergovernmental Panel on Climate Change (IPCC) (2018). Global Warming of 1.5°C: An IPCC Special Report. IPCC. Available at: https://www.ipcc.ch/site/assets/uploads/sites/2/2019/06/SR15_Chapter3_Low_Res.pdf. (Accessed January 18, 2022).

Google Scholar

Intergovernmental Panel on Climate Change (IPCC) (2019). Climate Change and Land. IPCC. Available at: https://www.ipcc.ch/site/assets/uploads/2019/11/SRCCL-Full-Report-Compiled-191128.pdf (Accessed January 20, 2022).

Google Scholar

Jakobsson, M. E., Moen, A., and Falnes, P. Ø. (2016). Correspondence: on the Enzymology and Significance of HSPA1 Lysine Methylation. Nat. Commun. 7 (1), 1–3. doi:10.1038/ncomms11464

CrossRef Full Text | Google Scholar

Jin, S.-G., Tsark, W., Szabó, P. E., and Pfeifer, G. P. (2008). Haploid Male Germ Cell‐and Oocyte‐specific Mbd3l1 and Mbd3l2 Genes are Dispensable for Early Development, Fertility, and Zygotic DNA Demethylation in the Mouse. Dev. Dyn. 237 (11), 3435–3443. doi:10.1002/dvdy.21767

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiyana, S. T. (2019). The Effects of Different Dietary Fibre Levels on Methane Production and Growth Performance of Bonsmara and Nguni Steers [Master’s thesis]. South Africa: University of South Africa.

Google Scholar

Linde, D. A. (2018). Gene Expression Profiles in Longissimus Dorsi Muscle of Bonsmara and Nguni Cattle Fed Two Different Diets. [Master’s thesis]. Pretoria: University of Pretoria. Available at: http://hdl.handle.net/2263/65906.

Google Scholar

Linington, M. J. (1990). The Use of Sanga Cattle in Beef Production. Tech. Commun. Depart. Agric. Dev. 223, 31–37.

Google Scholar

Madilindi, M. A., Banga, C. B., Bhebhe, E., Sanarana, Y. P., Nxumalo, K. S., Taela, M. G., et al. (2020). Genetic Diversity and Relationships Among Three Southern African Nguni Cattle Populations. Trop. Anim. Health Prod. 52 (2), 753–762. doi:10.1007/s11250-019-02066-y

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (8), e0200694. doi:10.1371/journal.pone.0200694

PubMed Abstract | CrossRef Full Text | Google Scholar

Makina, S. O., Muchadeyi, F. C., van Marle-Köster, E., MacNeil, M. D., and Maiwashe, A. (2014). Genetic Diversity and Population Structure Among Six Cattle Breeds in South Africa Using a Whole Genome SNP Panel. Front. Genet. 5, 333. doi:10.3389/fgene.2014.00333

PubMed Abstract | CrossRef Full Text | Google Scholar

Makina, S. O., Muchadeyi, F. C., van Marle-Köster, E., Taylor, J. F., Makgahlela, M. L., and Maiwashe, A. (2015). Genome-wide Scan for Selection Signatures in Six Cattle Breeds in South Africa. Genet. Sel. Evol. 47 (1), 1–14. doi:10.1186/s12711-015-0173-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Makina, S. O., Whitacre, L. K., Decker, J. E., Taylor, J. F., MacNeil, M. D., Scholtz, M. M., et al. (2016). Insight into the Genetic Composition of South African Sanga Cattle Using SNP Data from Cattle Breeds Worldwide. Genet. Sel. Evol. 48 (1), 1–7. doi:10.1186/s12711-016-0266-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Mastrangelo, S., Ben Jemaa, S., Ciani, E., Sottile, G., Moscarelli, A., Boussaha, M., et al. (2020). Genome‐wide Detection of Signatures of Selection in Three Valdostana Cattle Populations. J. Anim. Breed. Genet. 137 (6), 609–621. doi:10.1111/jbg.12476

PubMed Abstract | CrossRef Full Text | Google Scholar

Muchenje, V., DzamaChimonyo, K. M., Chimonyo, M., Raats, J. G., and Strydom, P. E. (2008). Meat Quality of Nguni, Bonsmara and Aberdeen Angus Steers Raised on Natural Pasture in the Eastern Cape, South Africa. Meat Sci. 79 (1), 20–28. doi:10.1016/j.meatsci.2007.07.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Mwai, O., Hanotte, O., Kwon, Y.-J., and Cho, S. (2015). African Indigenous Cattle: Unique Genetic Resources in a Rapidly Changing World. Asian Australas. J. Anim. Sci. 28 (7), 911–921. doi:10.5713/ajas.15.0002r

PubMed Abstract | CrossRef Full Text | Google Scholar

Nyamushamba, G. B., Mapiye, C., Tada, O., Halimani, T. E., and Muchenje, V. (2017). Conservation of Indigenous Cattle Genetic Resources in Southern Africa's Smallholder Areas: Turning Threats into Opportunities - A Review. Asian Australas J. Anim. Sci. 30 (5), 603–621. doi:10.5713/ajas.16.0024

PubMed Abstract | CrossRef Full Text | Google Scholar

Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A. R., Bender, D., et al. (2007). PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 81 (3), 559–575. doi:10.1086/519795

PubMed Abstract | CrossRef Full Text | Google Scholar

Pyoos, G. M. (2018). Crossbreeding Effects on Cow Efficiency and Component Traits. [Master’s thesis]. Bloemfontein: University of the Free State.

Google Scholar

R Core Team (2019). R Foundation for Statistical Computing. New York: Chapman and Hall/CRC.: R: A language and environment for statistical computing.

Google Scholar

Roshyara, N. R., and Scholz, M. (2014). fcGENE: A Versatile Tool for Processing and Transforming SNP Datasets. PloS one 9 (7), e97589. doi:10.1371/journal.pone.0097589

PubMed Abstract | CrossRef Full Text | Google Scholar

Sabeti, P. C., Varilly, P., Varilly, P., Fry, B., Lohmueller, J., Hostetter, E., et al. (2007). Genome-wide Detection and Characterization of Positive Selection in Human Populations. Nature 449 (7164), 913–918. doi:10.1038/nature06250

PubMed Abstract | CrossRef Full Text | Google Scholar

Saravanan, K. A., Panigrahi, M., Kumar, H., Parida, S., Bhushan, B., Gaur, G. K., et al. (2021). Genomic Scans for Selection Signatures Revealed Candidate Genes for Adaptation and Production Traits in a Variety of Cattle Breeds. Genomics 113 (3), 955–963. doi:10.1016/j.ygeno.2021.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (4), 629–644. doi:10.1086/502802

PubMed Abstract | CrossRef Full Text | Google Scholar

Scholtz, M. M., Schönfeldt, H. C., Neser, F. W. C., and Schutte, G. M. (2014). Research and Development on Climate Change and Greenhouse Gases in Support of Climate-Smart Livestock Production and a Vibrant Industry. South Afr. J. Animal Sci. 44 (5), 1–7.

Google Scholar

Sejian, V., Bhatta, R., Gaughan, J. B., Dunshea, F. R., and Lacetera, N. (2018). Adaptation of Animals to Heat Stress. Animal 12 (s2), s431–s444. doi:10.1017/S1751731118001945

PubMed Abstract | CrossRef Full Text | Google Scholar

Shike, D. W. (2013). Beef Cattle Feed Efficiency. Available at: https://dr.lib.iastate.edu/handle/20.500.12876/20785 (Accessed November 17, 2021).

Google Scholar

Singh, A., Mehrotra, A., Gondro, C., Romero, A. R. d. S., Pandey, A. K., Karthikeyan, A., et al. (2020). Signatures of Selection in Composite Vrindavani Cattle of India. Front. Genet. 11, 589496. doi:10.3389/fgene.2020.589496

PubMed Abstract | CrossRef Full Text | Google Scholar

Strydom, P. E., Naude, R. T., Smith, M. F., Kotze, A., Scholtz, M. M., and van Wyk, J. B. (2001). Relationships between Production and Product Traits in Subpopulations of Bonsmara and Nguni Cattle. South Afr. J. Animal Sci. 31 (3), 181–194. doi:10.4314/sajas.v31i3.3801

CrossRef Full Text | Google Scholar

Sun, H.-Z., Zhao, K., Zhou, M., Chen, Y., and Guan, L. L. (2019). Landscape of Multi-Tissue Global Gene Expression Reveals the Regulatory Signatures of Feed Efficiency in Beef Cattle. Bioinformatics 35 (10), 1712–1719. doi:10.1093/bioinformatics/bty883

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, K., Thornton, K. R., and Stoneking, M. (2007). A New Approach for Using Genome Scans to Detect Recent Positive Selection in the Human Genome. PLoS Biol. 5 (7), e171. doi:10.1371/journal.pbio.0050171

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, S., Slager, S. L., Gaonkar, K. S., and Yan, H. (2017). “Analyses of Genome-wide Histone Modifications in the Mammalian Genome,” in Handbook of Epigenetics. (Rochester: Academic Press), 135–152. doi:10.1016/B978-0-12-805388-1.00010-9

CrossRef Full Text | Google Scholar

Voight, B. F., Kudaravalli, S., Wen, X., and Pritchard, J. K. (2006). A Map of Recent Positive Selection in the Human Genome. PLoS Biol. 4 (3), e72. doi:10.1371/journal.pbio.0040072

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, D., Cai, J., Wang, B., Ding, S., Guan, L. L., and Liu, J. (2021). Integrative Network Analysis Revealed Molecular Mechanisms of Urine Urea Output in Lactating Dairy Cows: Potential Solutions to Reduce Environmental Nitrate Contamination. Genomics 113 (3), 1522–1533. doi:10.1016/j.ygeno.2021.03.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Wickham, H., Averick, M., Bryan, J., Chang, W., McGowan, L., François, R., et al. (2019). Welcome to the Tidyverse. J. Open Source Softw. 4 (43), 1686. doi:10.21105/joss.01686

CrossRef Full Text | Google Scholar

Wiederstein, J. L., Nolte, H., Günther, S., Piller, T., Baraldo, M., Kostin, S., et al. (2018). Skeletal Muscle-specific Methyltransferase METTL21C Trimethylates P97 and Regulates Autophagy-Associated Protein Breakdown. Cell Rep. 23 (5), 1342–1356. doi:10.1016/j.celrep.2018.03.136

PubMed Abstract | CrossRef Full Text | Google Scholar

Willemse, J., and Mare, F. (2016). The 2015 Drought and the Red Meat Industry. FarmBiz 2 (3), 8–13. Available at: https://hdl.handle.net/10520/EJC185523.

Google Scholar

Xu, L., Bickhart, D. M., Cole, J. B., Schroeder, S. G., Song, J., Tassell, C. P. V., et al. (2015). Genomic Signatures Reveal New Evidences for Selection of Important Traits in Domestic Cattle. Mol. Biol. Evol. 32 (3), 711–725. doi:10.1093/molbev/msu333

PubMed Abstract | CrossRef Full Text | Google Scholar

Yakovlev, A. F. (2018). Epigenetic Effects in Livestock Breeding. Russ. J. Genet. 54 (8), 897–909. doi:10.1134/S1022795418080148

CrossRef Full Text | Google Scholar

Yurchenko, A. A., Daetwyler, H. D., Yudin, N., Schnabel, R. D., Vander Jagt, C. J., Soloshenko, V., et al. (2018). Scans for Signatures of Selection in Russian Cattle Breed Genomes Reveal New Candidate Genes for Environmental Adaptation and Acclimation. Sci. Rep. 8 (1), 1–16. doi:10.1038/s41598-018-31304-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Haider, S., Baran, J., Cros, A., Guberman, J. M., Hsu, J., et al. (2011). BioMart: A Data Federation Framework for Large Collaborative Projects. Database 2011, bar038. doi:10.1093/database/bar038

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, M., Chaney, H. L., Current, J. Z., and Yao, J. (2021). Identification of the Core Promoter of ZNFO, an Oocyte-specific Maternal Effect Gene in Cattle. Gene 791, 145717. doi:10.1016/j.gene.2021.145717

PubMed Abstract | CrossRef Full Text | Google Scholar