Selection signatures of Fuzhong Buffalo based on whole-genome sequences

A comprehensive review of livestock development: insights into domestication, phylogenetics, diversity, and genomic advances

Livestock plays an essential role in sustaining human livelihoods, offering a diverse range of species integral to food security, economic stability, and cultural traditions. The domestication of livestock, which began over 10,000 years ago, has driven significant genetic changes in species such as cattle, buffaloes, sheep, goats, and pigs. Recent advancements in genomic technologies, including next-generation sequencing (NGS), genome-wide association studies (GWAS), and genomic selection, have dramatically enhanced our understanding of these genetic developments. This review brings together key research on the domestication process, phylogenetics, genetic diversity, and selection signatures within major livestock species. It emphasizes the importance of admixture studies and evolutionary forces like natural selection, genetic drift, and gene flow in shaping livestock populations. Additionally, the integration of machine learning with genomic data offers new perspectives on the functional roles of genes in adaptation and evolution. By exploring these genomic advancements, this review provides insights into genetic variation and evolutionary processes that could inform future approaches to improving livestock management and adaptation to environmental challenges, including climate change.

Identification of novel genetic loci related to dromedary camel (Camelus dromedarius) morphometrics, biomechanics, and behavior by genome-wide association studies

In the realm of animal breeding for sustainability, domestic camels have traditionally been valued for their milk and meat production. However, key aspects such as zoometrics, biomechanics, and behavior have often been overlooked in terms of their genetic foundations. Recognizing this gap, the present study perfomed genome-wide association analyses to identify genetic markers associated with zoometrics-, biomechanics-, and behavior-related traits in dromedary camels (Camelus dromedarius). 16 and 108 genetic markers were significantly associated (q < 0.05) at genome and chromosome-wide levels of significance, respectively, with zoometrics- (width, length, and perimeter/girth), biomechanics- (acceleration, displacement, spatial position, and velocity), and behavior-related traits (general cognition, intelligence, and Intelligence Quotient (IQ)) in dromedaries. In most association loci, the nearest protein-coding genes are linkedto neurodevelopmental and sensory disorders. This suggests that genetic variations related to neural development and sensory perception play crucial roles in shaping a dromedary camel's physical characteristics and behavior. In summary, this research advances our understanding of the genomic basis of essential traits in dromedary camels. Identifying specific genetic markers associated with zoometrics, biomechanics, and behavior provides valuable insights into camel domestication. Moreover, the links between these traits and genes related to neurodevelopmental and sensory disorders highlight the broader implications of domestication and modern selection on the health and welfare of dromedary camels. This knowledge could guide future breeding strategies, fostering a more holistic approach to camel husbandry and ensuring the sustainability of these animals in diverse agricultural contexts.

Whole-Genome Resequencing Analysis of the Camelus bactrianus (Bactrian Camel) Genome Identifies Mutations and Genes Affecting Milk Production Traits

Milk production is an important trait that influences the economic value of camels. However, the genetic regulatory mechanisms underlying milk production in camels have not yet been elucidated. We aimed to identify candidate molecular markers that affect camel milk production. We classified Junggar Bactrian camels (9-10-year-old) as low-yield (<1.96 kg/d) or high-yield (>2.75 kg/d) based on milk production performance. Milk fat (5.16 ± 0.51 g/100 g) and milk protein (3.59 ± 0.22 g/100 g) concentrations were significantly lower in high-yielding camels than those in low-yielding camels (6.21 ± 0.59 g/100 g, and 3.93 ± 0.27 g/100 g, respectively) (p < 0.01). There were no apparent differences in gland tissue morphology between the low- and high-production groups. Whole-genome resequencing of 12 low- and 12 high-yield camels was performed. The results of selection mapping methods, performed using two methods (FST and θπ), showed that 264 single nucleotide polymorphism sites (SNPs) overlapped between the two methods, identifying 181 genes. These genes were mainly associated with the regulation of oxytocin, estrogen, ErbB, Wnt, mTOR, PI3K-Akt, growth hormone synthesis/secretion/action, and MAPK signaling pathways. A total of 123 SNPs were selected, based on significantly associated genomic regions and important pathways for SNP genotyping, for verification in 521 additional Bactrian camels. This analysis showed that 13 SNPs were significantly associated with camel milk production yield and 18 SNPs were significantly associated with camel milk composition percentages. Most of these SNPs were located in coding regions of the genome. However, five and two important mutation sites were found in the introns of CSN2 (β-casein) and CSN3 (κ-casein), respectively. Among the candidate genes, NR4A1, ADCY8, PPARG, CSN2, and CSN3 have previously been well studied in dairy livestock. These observations provide a basis for understanding the molecular mechanisms underlying milk production in camels as well as genetic markers for breeding programs aimed at improving milk production.

Assessment of selective breeding effects and selection signatures in Qingyuan partridge chicken and its strains

Qingyuan partridge chicken (QYM) is a highly regarded native breed in China, highly esteemed for its exceptional breeding characteristics. However, the investigation into the selection signatures and its strains remains largely unexplored. In this study, blood sampling, DNA extracting, and high-depth resequencing were performed in 27 QYMs. Integrating the genomic data of 14 chicken (70 individuals) breeds from other researches, to analyze the genetic structure, selection signatures, and effects of selective breeding within QYM and its 3 strains (QYMA, QYMB, and QYMC). Population structure analysis revealed an independent QYM cluster, which exhibited distinct from other breeds, with each of its 3 strains displaying distinct clustering patterns. Linkage disequilibrium analysis highlighted QYMB's notably slower decay rate, potentially influenced by selection pressure from various production indicators. Examination of selection signatures uncovered genes and genetic mechanisms associated with genomic changes resulting from extensive selective breeding within the QYM and its strains. Intriguingly, diacylglycerol kinase beta (DGKB) and catenin alpha 2 (CTNNA2) were identified as commonly selected genes across the 3 QYM strains, linked to energy metabolism, muscle development, and fat metabolism. Our research validates the substantial impact of selective breeding on QYM and its strains, concurrently identifying genomic regions and signaling pathways associated with their distinctive characters. This research also establishes a fundamental framework for advancing yellow-feathered broiler breeding strategies.

Identification of genomic diversity and selection signatures in Luxi cattle using whole-genome sequencing data

OBJECTIVE

The objective of this study was to investigate the genetic diversity, population structure and whole-genome selection signatures of Luxi cattle to reveal its genomic characteristics in terms of meat and carcass traits, skeletal muscle development, body size, and other traits.

METHODS

To further analyze the genomic characteristics of Luxi cattle, this study sequenced the whole-genome of 16 individuals from the core conservation farm in Shandong region, and collected 174 published genomes of cattle for conjoint analysis. Furthermore, three different statistics (pi, Fst, and XP-EHH) were used to detect potential positive selection signatures related to selection in Luxi cattle. Moreover, gene ontology and Kyoto encyclopedia of genes and genomes pathway enrichment analyses were performed to reveal the potential biological function of candidate genes harbored in selected regions.

RESULTS

The results showed that Luxi cattle had high genomic diversity and low inbreeding levels. Using three complementary methods (pi, Fst, and XP-EHH) to detect the signatures of selection in the Luxi cattle genome, there were 2,941, 2,221 and 1,304 potentially selected genes identified, respectively. Furthermore, there were 45 genes annotated in common overlapping genomic regions covered 0.723 Mb, including PLAG1 zinc finger (PLAG1), dedicator of cytokinesis 3 (DOCK3), ephrin A2 (EFNA2), DAZ associated protein 1 (DAZAP1), Ral GTPase activating protein catalytic subunit alpha 1 (RALGAPA1), mediator complex subunit 13 (MED13), and decaprenyl diphosphate synthase subunit 2 (PDSS2), most of which were enriched in pathways related to muscle growth and differentiation and immunity.

CONCLUSION

In this study, we provided a series of genes associated with important economic traits were found in positive selection regions, and a scientific basis for the scientific conservation and genetic improvement of Luxi cattle.

Candidate genes for litter size in Xinjiang sheep identified by Specific Locus Amplified Fragment (SLAF) sequencing

The aim of this study was to investigate the selection signatures at a genome-wide level in 'Pishan' sheep using Specific Locus Amplified Fragment (SLAF)-seq. Blood samples from 126 ewes were sequenced using SLAF tags, and the ovarian tissues from 8 ewes (Bashbay sheep, a single litter size group (SG group); 'Pishan' sheep, double litter size group (DG group)) were collected to detect expression levels by quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR). Selection signature analysis was performed using global fixation index (Fst) and nucleotide diversity (π) ratio. A total of 1,192,168 high-quality SLAFs were identified. Notably, 2380 candidate regions under selection using two approaches were identified. A total of 2069 genes were identified, which were involved in dopaminergic synapses, thyroid hormone synthesis, ovarian steroidogenesis and thyroid hormone signalling pathways. Furthermore, Growth Differentiation Factor 9 (GDF9), Period Circadian Regulator 2 (PER2), Thyroid Stimulating Hormone Receptor (TSHR), and Nuclear Receptor Coactivator 1 (NCOA1) reside within these regions and pathways. The expression levels of GDF9 and PER2 genes in sheep tissue of the DG group were significantly higher than those in the SG group. These genes are interesting candidates for litter size and provide a starting point for further identification of conservation strategies for 'Pishan' sheep.

Effects of different load weights on the work performance and physiological and hematobiochemical responses in working water buffalo

Background and Aim

Working animals are important in agriculture because they play a role in various agricultural activities, including milk and meat production. Thus, they contribute to the development of rural communities. In this study, we aimed to evaluate the effects of different load weights on the work performance and the physiological and hematological responses of working water buffalo (Bubalus bubalis).

Materials and Methods

The work performances of 12 buffaloes (average weight 782.16 ± 21.62 kg) transporting 200, 350, and 500 kg of African palm fruits in metal baskets placed on their backs were evaluated. Work performance variables evaluated immediately after work were as follows: total number of trips (n), total weight (kg), distance traveled (km/day), working period (h), stopped time (h), and average speed (km/h). In addition, we evaluated physiological, biochemical, and hematological variables at three different times: before the start of work, immediately after work, and on the rest day.

Results

Among the load weights, 500 kg (total load carried = 4,138.88 kg) improved work efficiency compared to 200 kg loads (total load = 3,322.22 kg) (p = 0.0281). However, 500 kg loads resulted in slower average speed (2.4 km, p = 0.0164), shorter working period (2.39 h, p < 0.0001) and distance traveled (7.29 km, p < 0.0001), and less total number of trips (8.27 trips, p < 0.0001) compared to 350 and 200 kg load weights (3.45 and 3.52 km/h, 2.55 and 2.79 h, 8.71 and 9.75 km, 10.94 and 16.61 trips, respectively); and the heaviest loads resulted in significantly higher (p < 0.005) respiratory rate, pulse, heart rate, rectal temperature, glucose, lactate dehydrogenase, creatine kinase, total protein, white blood cell count, neutrophils, lymphocytes, monocytes, eosinophils, and basophils. In contrast, was associated with lower levels of red blood cells, hemoglobin, and hematocrit compared to lower loads. All differences were more pronounced in 500 kg (p < 0.005) compared to 200 kg loads.

Conclusion

Working buffaloes responded to work related to the transport of African palm fruits through various physiological, biochemical, and hematological changes. However, some variables remained close to the reference values reported in the literature for water buffaloes, and in general, all variables were reestablished during the rest day, indicating that these animals have adapted to working conditions.

Whole genome sequencing reveals signals of adaptive admixture in Creole cattle

The Creole cattle from Guadeloupe (GUA) are well adapted to the tropical environment. Its admixed genome likely played an important role in such adaptation. Here, we sought to detect genomic signatures of selection in the GUA genome. For this purpose, we sequenced 23 GUA individuals and combined our data with sequenced genomes of 99 animals representative of European, African and indicine groups. We detect 17,228,983 single nucleotide polymorphisms (SNPs) in the GUA genome, providing the most detailed exploration, to date, of patterns of genetic variation in this breed. We confirm the higher level of African and indicine ancestries, compared to the European ancestry and we highlight the African origin of indicine ancestry in the GUA genome. We identify five strong candidate regions showing an excess of indicine ancestry and consistently supported across the different detection methods. These regions encompass genes with adaptive roles in relation to immunity, thermotolerance and physical activity. We confirmed a previously identified horn-related gene, RXFP2, as a gene under strong selective pressure in the GUA population likely owing to human-driven (socio-cultural) pressure. Findings from this study provide insight into the genetic mechanisms associated with resilience traits in livestock.

Analysis of genomic copy number variations through whole-genome scan in Chinese Qaidam cattle

Qaidam cattle (CDM) are indigenous breed inhabiting Northwest China. In the present study, we newly sequenced 20 Qaidam cattle to investigate the copy number variants (CNVs) based on the ARS-UMD1.2 reference genome. We generated the CNV region (CNVR) datasets to explore the genomic CNV diversity and population stratification. The other four cattle breeds (Xizang cattle, XZ; Kazakh cattle, HSK; Mongolian cattle, MG; and Yanbian cattle, YB) from the regions of North China embracing 43 genomic sequences were collected and are distinguished from each of the other diverse populations by deletions and duplications. We also observed that the number of duplications was significantly more than deletions in the genome, which may be less harmful to gene formation and function. At the same time, only 1.15% of CNVRs overlapped with the exon region. Population differential CNVRs and functional annotations between the Qaidam cattle population and other cattle breeds revealed the functional genes related to immunity (MUC6), growth (ADAMTSL3), and adaptability (EBF2). Our analysis has provided numerous genomic characteristics of some Chinese cattle breeds, which are valuable as customized biological molecular markers in cattle breeding and production.

Genome-Wide Association Study for Udder Conformation Traits in Chinese Holstein Cattle

Udder conformation traits are one of the most economic traits in dairy cows, greatly affecting animal health, milk production, and producer profitability in the dairy industry. Genetic analysis of udder structure and scores have been developed in Holstein cattle. In our research, we conducted a genome-wide association study for five udder traits, including anterior udder attachment (AUA), central suspensory ligament (CSL), posterior udder attachment height (PUAH), posterior udder attachment width (PUAW), and udder depth (UD), in which the fixed and random model circulating probability unification (FarmCPU) model was applied for the association analysis. The heritability and the standard errors of these five udder traits ranged from 0.04 ± 0.00 to 0.49 ± 0.03. Phenotype data were measured from 1000 Holstein cows, and the GeneSeek Genomic Profiler (GGP) Bovine 100 K SNP chip was used to analyze genotypic data in Holstein cattle. For GWAS analysis, 984 individual cows and 84,407 single-nucleotide polymorphisms (SNPs) remained after quality control; a total of 18 SNPs were found at the GW significant threshold (p < 5.90 × 10−7). Many candidate genes were identified within 200kb upstream or downstream of the significant SNPs, which include MGST1, MGST2, MTUS1, PRKN, STXBP6, GRID2, E2F8, CDH11, FOXP1, SLF1, TMEM117, SBF2, GC, ADGRB3, and GCLC. Pathway analysis revealed that 58 Gene Ontology (GO) terms and 18 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched with adjusted p values, and these GO terms and the KEGG pathway analysis were associated with biological information, metabolism, hormonal growth, and development processes. These results could give valuable biological information for the genetic architecture of udder conformation traits in dairy Holstein cattle.

Comparative analyses of copy number variations between swamp and river buffalo

Domestic buffalo is an important livestock in the tropical and sub-tropical region, including two types: swamp and river buffalo. The swamp buffalo is mainly used as draft animal, while the river buffalo is raised for milk production. In this study, based on the new high-quality buffalo reference genome UOA_WB_1, we firstly investigated the copy number variants in buffalo using whole-genome Illumina sequencing. A total of 3,734 CNV regions (CNVRs) were detected in 106 buffalo population with a total length of 23,429,066 bp, corresponding to ∼ 0.88% of the water buffalo genome (UOA_WB_1). Our results revealed a clear population differentiation in CNV between swamp and river buffalo. In addition, a total of 667 highly differentiated CNVRs (covering 886 genes) were detected between river and swamp buffalo population. We detected a set of CNVR-overlapping genes associated with exercise, immunity, nerve, and milk trait which exhibited different copy numbers between swamp and river buffalo population. This study provides valuable genome variation resources for buffalo and would contribute to understanding the genetic differences between swamp and river buffalo.

Whole-Genome Resequencing of Xiangxi Cattle Identifies Genomic Diversity and Selection Signatures

Understanding the genetic diversity in Xiangxi cattle may facilitate our efforts toward further breeding programs. Here we compared 23 Xiangxi cattle with 78 published genomes of 6 worldwide representative breeds to characterize the genomic variations of Xiangxi cattle. Based on clustering models in population structure analysis, we displayed that Xiangxi cattle had a mutual genome ancestor with Chinese indicine, Indian indicine, and East Asian taurine. Population genetic diversity was analyzed by four methods (nucleotide diversity, inbreeding coefficient, linkage disequilibrium decay and runs of homozygosity), and we found that Xiangxi cattle had higher genomic diversity and weaker artificial selection than commercial breed cattle. Using four testing methods (θπ, CLR, F_ST, and XP-EHH), we explored positive selection regions harboring genes in Xiangxi cattle, which were related to reproduction, growth, meat quality, heat tolerance, and immune response. Our findings revealed the extent of sequence variation in Xiangxi cattle at the genome-wide level. All of our fruitful results can bring about a valuable genomic resource for genetic studies and breed protection in the future.

Assessing Genomic Diversity and Productivity Signatures in Dianzhong Cattle by Whole-Genome Scanning

Dianzhong cattle is a classic Chinese indigenous cattle breed with historical records dating back to 200 BC. But with its genomic differences having not been clearly elucidated, the quest for genomic characterization will be an essential step towards understanding the genomic basis of productivity and adaptation to survival under Chinese farming systems. Here we compared 10 Dianzhong cattle (four newly sequenced and six downloaded) with 29 published genomes of three underlying ancestral populations (Chinese zebu, Indian zebu, and Yanbian cattle) to characterize the genomic variations of Dianzhong cattle. Dianzhong cattle has a high nucleotide diversity (0.0034), second only to Chinese zebu. Together with analyses of linkage disequilibrium decay and runs of homozygosity, Dianzhong cattle displayed higher genomic diversity and weaker artificial selection compared with Yanbian cattle. From a selective sweep analysis by four methods (Fst, π-ratio, XP-CLR, and XP-EHH), the positive selective signals were mainly manifested in candidate genes and pathways related to heat resistance, growth and development, fat deposition, and male reproduction. Missense mutations were detected in candidate genes, SDS (c.944C > A and p.Ala315Glu), PDGFD (c.473A > G and p.Lys158Arg), and DDX4 (rs460251486, rs722912933, and rs517668236), which related to heat resistance, fat deposition, and spermatogenesis, respectively. Our findings unravel, at the genome-wide level, the unique diversity of Dianzhong cattle while emphasizing the opportunities for improvement of livestock productivity in further breeding programs.