Genetic variation analysis of Guanling cattle based on whole-genome resequencing

Article information

Anim Biosci. 2024;37(12):2044-2053

Publication date (electronic) : 2024 June 27

doi : https://doi.org/10.5713/ab.24.0181

Longxin Xu¹, Xin Wang¹, Junda Wu¹, Hua Wang¹, Wenzhang Zhou¹, Jing Liu¹, Mengmeng Ni^2,3, Kaikai Zhang¹, Bo Yu^1,, Ruiyi Lin^4,

¹Institute of Animal Husbandry and Veterinary Science, Guizhou Academy of Agricultural Sciences, Guiyang 550005, China

²College of Animal Sciences, Guizhou University, Guiyang 550000, China

³Guizhou Yellow Cattle Industry Group Co., Ltd, Guiyang 550001, China

⁴College of Animal Sciences, Fujian Agriculture and Forestry University, Fuzhou 350002, China

^*Corresponding Authors: Bo Yu, Tel: +86-15085916661, E-mail: yuboky@163.com, Ruiyi Lin, Tel: +86-14759103866, E-mail: linruiyi@fafu.edu.cn

Received 2024 March 26; Revised 2024 May 12; Accepted 2024 June 18.

Abstract

INTRODUCTION

In the extensive history of cattle breeding in Guizhou, a diverse range of local cattle breeds has developed, with the Guanling cattle being a prominent example. These cattle primarily inhabit the vast mountainous regions within the Panjiang River basin, spanning across Guizhou, Yunnan, and northern Guangxi, with Guanling County being their most renowned habitat [1]. Guanling cattle represent a superior local breed in Guizhou Province and are among the 78 nationally protected local livestock breeds. They are highly regarded for their culinary, economic, and developmental potential. Previous research has employed biochemical and molecular techniques to analyze the genetic characteristics of the Guanling breed [2]. These studies have revealed significant genetic variation and growth potential. However, the advancement and utilization of these breeds have faced limitations attributed to the introduction of foreign breeds, crossbreeding enhancements, environmental degradation, and various natural factors, which have led to a decline in their population.

Advances in sequencing technologies have greatly facilitated comprehensive and in-depth genome analysis. And, sequencing of the bovine genome and HapMap projects have revealed a substantial number of genetic variations [3], with single nucleotide polymorphisms (SNPs) being the most extensively studied variant. SNPs serve as valuable tools for identifying genomic regions through association analysis, as they exhibit linkage disequilibrium with quantitative trait loci influencing target traits, a phenomenon observed in various animal species. For instance, Eck et al [4] identified 2.4 million SNPs in Holstein cattle using the Illumina HiSeq platform. Stothard et al [5] employed SOLiD technology to successfully map genomic variations between Black Angus bulls and Holstein bulls, identifying approximately 7 million SNPs and 790 copy number variations. The utilization of SNPs as selection criteria for meat traits in marker-assisted selection can significantly enhance cattle selection and breeding programs.

Recent studies on bovine genomic variation have been extensive, however, whole-genome studies on Guanling cattle have yet to be reported. Here, we resequenced the Guanling cattle and compared them with the Limousin, Simmental, and Angus breeds to reveal their genomic characteristics and variations. This study’s primary objective was to unravel the genetic traits of Guanling cattle, pinpoint genes advantageous for muscle growth, and provide essential genomic data to support further analysis of genetic mechanisms tied to economic traits and the preservation of cattle breed genetic diversity.

MATERIALS AND METHODS

Quality control and data filtering

In order to ensure the quality of the data, the sequencing raw data should go through quality control and data filtering. The quality of the raw data (three Guanling cattle, one of each Simmental, Angus, and Limousin cattle) is controlled by analyzing the composition and quality values of the bases (Table 1). The raw reads with some joint or low quality reads was filtered by fastp [6] to obtain high quality clean data, and subsequent analysis was based on clean data. Data filtering is mainly about removing paired reads with joints; removing paired reads with N base (N indicates uncertain base information) greater than 10%; removing paired reads with low quality (mass value Q 7) bases exceeding 30% of the total number of reads.

Table 1

Quality evaluation of sequencing data

RESULTS

Sequencing and alignment

A quality assessment was performed on the initial data generated from whole-genome sequencing (WGS), which included the removal of potential PCR duplicates and realignment around insertions and deletions. After duplicate removal, an average data output of 26 Gb was obtained. The alignment rate of all the samples to the reference genome exceeded 99% (Table 2), with an average coverage of more than 97% (Table 3).

Table 2

The result of sample mapping

Table 3

The result of sequencing depth and coverage

Functional annotation and potential functional exploration of non-synonymous single nucleotide polymorphisms

A total of 166,411 non-synonymous SNPs (nsSNPs) were identified, and their mutations potentially play a pivotal role in altering economic traits in Guanling cattle. All of these mutations were annotated and subjected to GO functional enrichment and Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analyses. GO functional enrichment analysis revealed 251 terms, mainly including adenyl ribonucleotide binding, ATP binding, adenyl nucleotide binding and so on (Figure 1A).

Figure 1

Enrichment analysis of nsSNPs. (A) Gene ontology enrichment of common SNPs between the two groups; (B) Kyoto encyclopedia of genes and genomes enrichment of common SNPs between the two groups. nsSNPs, non-synonymous single nucleotide polymorphisms; SNP, single nucleotide polymorphism.

While KEGG pathway enrichment analysis revealed 60 pathways associated with genes whose expression was significantly enriched (p<0.05), including Extracellular matrix organization, Focal adhesion, extracellular matrix (ECM) proteoglycans, ECM-receptor interaction, Non-integrin membrane-ECM interactions, Transmembrane transport of small molecules and so on. Notably, the “Transmembrane transport of small molecules” pathway had the highest number of genes involved, with the highest enrichment found in the “Extracellular matrix organization” (Figure 1B).

Indel detection

Typically, deletions and insertions over 50 base pairs (bp) are considered structural variations, while those under 50 bp are collectively referred to as insertions and deletions (indels) [14]. In Guanling cattle, a total of 42,423 indels were identified, comprising 22,101 deletions and 20,322 insertions, accounting for 52.09% and 47.91%, respectively. GO functional enrichment analysis yielded 105 terms, including autophagosome maturation, kinesin complex, vacuole fusion, microtubule associated complex, intracellular ligand-gated ion channel activity and so on (Figure 2A). After significant enrichment analysis of the KEGG pathways, 48 pathways were found to be notably enriched (p<0.05), including Intra-Golgi and retrograde Golgi-to-ER traffic, COPI-dependent Golgi-to-ER retrograde traffic, Kinesins and so on (Figure 2B).

Figure 2

Enrichment analysis of indels. (A) Represents the gene ontology enrichment of common indels between the two groups; and (B) represents the Kyoto encyclopedia of genes and genomes enrichment of common indels between the two groups.

Candidate gene selection

Unique 1,520 homozygous SNP loci between the two groups were screened, excluding those on sex chromosomes and mitochondria. This led to the identification of 91 loci related to meat quality traits. We further selected loci with missense mutations and utilized the Ensemble database for annotating non-synonymous SNPs, resulting in the identification of 7 loci related to meat quality traits involving 6 different genes (Table 4).

Table 4

Single nucleotide polymorphism associated with meat quality traits

Single nucleotide polymorphism detection

To determine the distribution of the aforementioned 7 SNPs within the Guanling cattle population, high-multiplex PCR was employed for locus scanning in 55 cattle. The results (Table 5) showed that all 7 loci were in Hardy–Weinberg equilibrium. The identified mutation loci in the LEPR, SIX4, SPIDR, and FASN genes were dominant alleles in the Guanling cattle population.

Table 5

Detection of single nucleotide polymorphism locus

DISCUSSION

The number of Guanling cattle is small, so it is particularly important to select individuals which can represent varieties for sequencing. In order to avoid individual differences, more individuals are gathered at a low cost to reflect the population genetic diversity of Guanling cattle varieties. Therefore, three Guanling cattle samples were used for resequencing, and 55 Guanling cattle samples were used for further SNP detection. Finally, we get an average of 26 Gb raw data, 99.5% reads aligned to the reference genome, with high single base correctness, similar to the previous sequencing results of ordinary cattle [4,5]. High sequencing depth, and the detected variants are fully credible [15].

Here, we found 598,688 SNPs and 42,423 indels in 29 autosomes and X chromosomes. Of the total SNP, heterozygous SNPs 574,082 (95.89%), homozygous SNPs 24,606 (4.11%) and the ratio of heterozygous/homozygous SNPs was 23.33, significantly higher than that of Japanese cattle (1.24) [16] and Korean cattle (1.63) [17]. From a sequencing perspective, pure sequences and SNPs exhibit distinct characteristics. In mixed samples, the base types at a particular site are consistent with the reference genome. However, heterozygous SNPs indicate multiple base types at the same site across all mixed samples. Guanling cattle display a low homozygosity ratio for SNPs, suggesting high heterozygosity. This observation may reflect sequencing variations among individuals and could also be attributed to low levels of cattle breeding. Additionally, there appears to be increased gene communication between Guanling cattle and other breeds, potentially resulting in the loss of specific functional genes. Protecting the genetic diversity of Guanling cattle varieties is therefore of paramount importance. By comparing SNPs and indels, 13,769 new SNPs and 2,206 new indels were found, accounting for 2.3% of the total SNPs and 5.2% of the total indels, respectively. Due to the development of genome sequencing in recent years, more new SNPs and indels have been found, and the database is more and more perfect. So, the proportion of comparison is significantly increased, and the number of new discoveries gradually less. Most indels are short in length, with deletions ranging from 1 to 29 kb, insertions from 1 to 44 kb, the number of deletions and insertions concentrated in 1 to 10 bp, with 1 to 3 bp being the most, and similar phenomena are observed from human genome data [18]. In Guanling cattle data, nearly 85.6% insertions and 77.9% deletions were less than 5 kb. SNPs and indels detected on 29 autosomes are proportional to chromosome length and the expected results with the lowest X chromosome mutation rate of 4.61%, and on small population studies, the X chromosome has a lower mutation rate compared to the autosomes [19].

Guanling cattle are mainly used for cultivated land, and have gradually developed towards the direction of serving meat. The LEPR gene is located on bovine chromosome 3 and encodes a protein belonging to the cytokine receptor gp130 family. It plays a crucial role in regulating fat metabolism and is a novel hematopoietic pathway essential for normal lymphocyte production [20]. Previous research by Raza et al [21] investigated that three SNPs of LEPR gene were associated with backfat thickness, and intramuscular fat content in Chinese beef cattle breeds. Moreover, the expression level of LEPR was significant difference in the longissimus dorsi muscle of Yunling and Simmental cattle [22]. On bovine chromosome 4, the AKAP9 gene, a member of the AKAP family, encodes a kinase anchoring protein. Raza et al [10] findings suggest that bovine AKAP9 may be involved in regulating fat formation, growth traits, differentiation of adipose tissue, regeneration of skeletal muscle, and metabolism.

Like the bovine chromosome 10 gene, the SIX4 gene, which is part of the Sine Oculis/Six gene family, plays a crucial role in skeletal muscle development. Genetic variations or deletions in SIX4 can have implications for pituitary function [23]. Wang et al [24] research revealed significant correlations between three SNPs in the SIX4 gene of Qinchuan cattle and body measurement traits, suggesting that SIX4 is a candidate gene influencing cattle body size traits. Huang et al [25] demonstrated that SPIDR regulates the assembly or stability of RAD51/DMC1 on single-stranded DNA, a vital recombination factor in meiotic recombination in mammals. Zhang et al [26] demonstrated the significant effects of the SPIDR gene on growth parameters and carcass traits in 1,173 Chinese Simmental beef cattle.

Along bovine chromosome 16, the PRG4 gene encodes a proteoglycan-like glycoprotein synthesized by various tissues, including joint cartilage, the meniscus, the synovial lining, and tendon cells [27]. Research by Abubacker et al [28] emphasized the importance of the intermolecular disulfide bond polymers of PRG4 for its adsorption onto cartilage surfaces and function as a boundary lubricant. These findings suggested that PRG4 may play a role in reducing friction in tissues and spaces of bovine knee joints, contributing to effective weight-bearing systems, particularly during weight gain in Guanling cattle.

Fatty acid synthase (FASN) encodes a versatile protein primarily responsible for catalyzing the synthesis of palmitic acid from acetyl-CoA and malonyl-CoA, facilitated by nicotinamide adenine dinucleotide phosphate (NADPH), resulting in the formation of long-chain saturated fatty acids [29]. Previous studies have identified FASN as a pivotal candidate gene influencing the composition of fat in both milk and meat [30,31]. Chu et al [32] observed significant variations in intramuscular fat content among Datong yaks based on different FASN gene genotypes. Individuals with the HH and HG genotypes exhibited notably greater intramuscular fat content than did those with the GG genotype [32]. FASN mRNA expression levels in subcutaneous fat and abdominal fat in Yan yellow cattle were significantly higher than that in Yanbian yellow cattle [33]. Genome-wide association analysis revealed that g.841G>C SNP of FASN gene showed significant associations with the percentages of C14:0, C14:1, C16:1, and C18:1 at 5% genome-wide significance level in Japanese Black cattle [34]. In our study, we screened seven SNP sites spanning six genes within the Guanling cattle population. The dominant alleles in this population were found to be A at rs43347904 within the LEPR gene, A at rs133120166, T at rs109170670 within the SIX4 gene, C at rs208094969 within the SPIDR gene, and A at rs715140536 within the FASN gene. However, the biological functions of these mutation sites require further investigation.

CONCLUSION

This study conducted genome-wide genotyping of Guanling cattle, which yielded an abundance of genetic markers. These markers, which are crucial for selection purposes, serve as powerful tools for enhancing the breeding and management of Guanling cattle, offering valuable insights for breeding decisions. The integration of resequencing with advanced breeding techniques to identify gene polymorphisms and characteristics holds great promise in elevating the breeding standards of economic animals, with a specific focus on improving growth and meat quality traits. Our research successfully pinpointed several genes which may influence meat quality traits, thus contributing essential knowledge about the variety identification and optimal breeding strategies for the Guanling cattle breed.

Article information Continued