Selection signatures of Fuzhong Buffalo based on whole-genome sequences

A Multifaceted Giant Protein Microtubule-Actin Cross-Linking Factor 1

Microtubule-actin cross-linking factor 1 (MACF1), also known as actin cross-linking family protein 7 (ACF7), is a giant cytolinker protein with multiple conserved domains that can orchestrate cytoskeletal networks of actin and microtubules. MACF1 is involved in various biological processes, including cell polarity, cell-cell connection, cell proliferation, migration, vesicle transport, signal transduction, and neuronal development. In this review, we updated the physiological and pathological roles of MACF1, highlighting the components and signaling pathways involved. Novel evidence showed that MACF1 is involved in diverse human diseases, including multiple neuronal diseases, congenital myasthenic syndrome, premature ovarian insufficiency, spectraplakinopathy, osteoporosis, proliferative diabetic retinopathy, and various types of cancer. We also reviewed the physiological roles of MACF1, including its involvement in adhesome formation, bone formation, neuronal aging, and tooth development. In addition, MACF1 plays other roles, functioning as a biomarker for the prediction of infections in patients with burns and as a marker for genome selection breeding. These studies reinforce the idea that MACF1 is a bona fide versatile, multifaceted giant protein. Identifying additional MACF1 functions would finally help with the treatment of diseases caused by MACF1 defects.

BuffExDb: web-based tissue-specific gene expression resource for breeding and conservation programmes in Bubalus bubalis

Amidst the global challenge of extreme poverty, the livestock sector can significantly contribute to global sustainable development goals by enhancing resilience, smallholder productivity, and market participation. The Indian livestock sector is one of the largest in the world with a total livestock population of 535.82 million, ∼10.7% of the world's livestock population. Buffalo (Bubalus bubalis) holds significant importance in India and other Asian countries, notably contributing to their economies by surpassing cattle in milk production and providing various valuable products. The limited availability of genomic and transcriptomic resources for buffaloes hinders the efforts to enhance their traits for increased milk and meat production. To address this gap, this study adopted the state-of-the-art bioinformatics tools to analyse 2429 transcriptomes representing 438 BioSamples from 23 BioProjects obtained from a public domain database, representing 76 different types of tissues and cell types from all major organ systems in buffalo species (river and swamp). The outcome of this exhaustive genomic data led to the development of a relational buffalo expression database based on a three-tier architecture named as BuffExDb (http://46.202.167.198/buffex/). The user-friendliness and flexibilities in retrieval of tissue-specific genes (TSGs) and their functional annotation are the major characteristics of BuffExDb. This is the first of its kind that offers an effortlessly navigable and filterable database, enabling users to examine and visualize the expression levels of each tissue across multiple samples, simultaneously. It also provides the Tau score parameter for the identification of TSGs along with their essential roles in tissue development, maintenance, and function as observed through the enrichment test for gene ontologies. The exhaustive outcome of this work would pave the way for the biological, functional, and evolutionary studies for easy access. This prior information based on tissue-specific mechanisms can be used for future genomic research, especially in association studies in endeavour of enhanced buffalo breeding and conservation programmes. Database URL: http://46.202.167.198/buffex/.

Effects of moderate intensity exercise on liver metabolism in mice based on multi-omics analysis

Physical exercise is beneficial to keep physical and mental health. The molecular mechanisms underlying exercise are still worth exploring. The healthy adult mice after six weeks of moderate-intensity exercise (experimental group) and sedentary mice (control group) were used to perform transcriptomic, proteomic, lactylation modification, and metabolomics analysis. In addition, gene sets related to hypoxia, glycolysis, and fatty acid metabolism were used to aid in the screening of hub genes. The mMCP-counter was employed to evaluate infiltration of immune cells in murine liver tissues. Transcriptomics analysis revealed 82 intersection genes related to hypoxia, glycolysis, and fatty acid metabolism. Proteomics and lactylation modification analysis identified 577 proteins and 141 differentially lactylation modification proteins. By overlapping 82 intersection genes with 577 differentially expressed proteins and 141 differentially lactylation modification proteins, three hub genes (Aldoa, Acsl1, and Hadhb) were obtained. The immune infiltration analysis revealed a decreased score for monocytes/macrophages and an increased score for endothelial cells in the experimental group. Then, 459 metabolites in positive mode and 181 metabolites in negative mode were identified. The "Metabolic pathways" (mmu01100) was a common pathway between intersection genes-enriched pathways and metabolites-enriched pathways. These findings highlight the pivotal roles of hub genes in the glycolysis and fatty acid metabolism under the context of chronic exercise.

A multi-omics database of buffaloes from Yangtze valley reveals diversity of water buffalo (Bubalus bubalis)

Asian water buffalo (Bubalus bubalis) is the fundamental livestock resource for local rural populations and holds a promising prospect of their milk and meat. Xuyi mountain (XYM) and Haizi (HZ) buffaloes from Yangtze valley comprises of species diversity of Asian water buffaloes. Current multi-omics enables identification of causal genes and elucidation of genetic regulatory mechanisms underlying complex traits in buffaloes. Here, we conducted the integrated analysis of metabolome and metagenome of rumen fluid, transcriptome and metabolome of blood, and whole genome sequence data from XYM (n = 7) and HZ (n = 10) male buffaloes. Our results revealed the apparent diversity of multi-layer omics profiles between two buffalo species. The built-up multi-omics database supports the discoveries of diversity in Asian water buffalo and potentially serves valuable resources for studying causal regulatory variants and their mechanisms.

Selection signature analysis reveals genes associated with tail phenotype in sheep

Tail type of sheep, which may be affected by many genes with a complex mechanism, is an important economic trait concerned by both raiser and consumers. Here, we employed two sheep breeds with extreme phenotypes - Mongolian sheep (short-fat-tailed) and Bamei Mutton sheep (long-thin-tailed) to analyze the genetic differences at the genomic level and find candidate genes associated with tail phenotype. The results of population structure analysis showed that the LD decay rate of Mongolian sheep was greater than that of Bamai Mutton sheep. When K = 2, the two populations were obviously separated with a certain degree of mixing. From 49 sheep individuals, 20,270,930 and 2,479,474 SNPs and Indels were identified, respectively. Selection signals were detected based on F ST , π-Ratio, and XP-EHH. These three methods identified 85 candidate genes, of which PDGFD, GLIS1, AR, and FGF9 were reported to be associated with tail fat deposition, while VRTN associated with tail length in sheep tail phenotype; the others were novel genes that may play important roles in sheep tail phenotype formation. Gene annotation revealed that these candidate genes mainly participate in pathways associated with fat deposition or lipid metabolism. This study provided insight into sheep tail type development and a guide for molecular breeding.

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.

The genome landscape of the Xinglong buffalo

BACKGROUND

Xinglong buffalo, as an indigenous breed in Hainan province of China, possesses characteristics such as high humidity tolerance, disease resistance and high reproductive capacity. Combined with whole genome sequencing technology, comprehensive investigation can be undertaken to elucidate the genomic characteristics, functions and genetic variation of Xinglong buffalo population.

RESULTS

Xinglong buffalo has the highest genetic diversity, lowest runs of homozygosity average length, and fasted decay of linkage disequilibrium in our study population. Phylogenetic tree results revealed that Xinglong buffalo was gathered together with Fuzhong buffalo firstly. The population genetic structure analysis indicates that at K = 3, the Xinglong buffalo for the first time showed a distinct ancestral origin from other water buffalo. Furthermore, compared to different populations, candidate genes displaying significantly distinct patterns of single nucleotide polymorphisms (SNPs) (e.g., RYR2, COX15, PCDH9, DTWD2, FCRL5) distribution have been identified in the Xinglong buffalo.

CONCLUSIONS

Based on the whole genome sequencing data, this study identified a substantial number of SNPs and assessed the genetic diversity and selection signatures within the Xinglong buffalo population. These results contribute to understanding the genomic characteristics of Xinglong buffalo and their genetic evolutionary status. However, the practical significance of these signatures for genetic enhancement still requires confirmation through additional samples and further experimental validation.

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.