Genetic architectures and selection signatures of body height in Chinese indigenous donkeys revealed by next-generation sequencing

Comprehensive genomic analysis of genetic diversity, body size, and origins of the Hetian Gray donkey

BACKGROUND

The Hetian Gray donkey, characterized by its large stature and gray coat color, is a relatively recently identified breed. This breed plays an important role in local agricultural practices and cultural heritage. However, with the widespread adoption of mechanization, its use as a draft animal has gradually diminished. Additionally, the insufficient development of its meat and milk products has contributed to a continuous decline in its population, making its conservation increasingly urgent. Despite its significance, research on the origin and ancestral composition of the Hetian Gray donkey remains limited, which is essential for developing effective conservation strategies. To address this, whole-genome resequencing data of 208 donkeys from seven Chinese indigenous donkey breeds, Equus africanus somaliensis, and Equus kiang were analyzed to investigate the relationships between the Hetian Gray donkey and other breeds in the present study. The findings of this study will provide valuable scientific evidence for preserving the genetic diversity and unique characteristics of the Hetian Gray donkey.

RESULTS

The analysis revealed that wild donkeys had no significant impact on the formation of the Hetian Gray donkey. Among the breeds studied, the Xinjiang donkey and the Guanzhong donkey exhibit the closest genetic affinity with the Hetian Gray donkeys. Additionally, genome-wide selection signature analysis between the Hetian Gray donkey and the Xinjiang donkey identified positive selection for genes associated with body size traits.

CONCLUSION

In summary, the results confirm that geographical distance plays a crucial role in determining the genetic relationships among the Hetian Gray donkey and other populations. The Xinjiang donkey and the Guanzhong donkey played an important role in the formation of the Hetian Gray donkey, and genes related to body size were positively selected during development of the breed.

Whole-genome sequencing reveals patterns of runs of homozygosity underlying genetic diversity and selection in domestic rabbits

BACKGROUND

Runs of homozygosity (ROH) are continuous segments of homozygous genotypes inherited from both parental lineages. These segments arise due to the transmission of identical haplotypes. The genome-wide patterns and hotspot regions of ROH provide valuable insights into genetic diversity, demographic history, and selection trends. In this study, we analyzed whole-genome resequencing data from 117 rabbits to identify ROH patterns and inbreeding level across eleven rabbit breeds, including seven Chinese indigenous breeds and four exotic breeds, and to uncover selective signatures based on ROH islands.

RESULTS

We detected a total of 31,429 ROHs across the autosomes of all breeds, with the number of ROHs (NROH) per breed ranging from 1316 to 7476. The mean sum of ROHs length (SROH) per individual was 493.84 Mb, covering approximately 22.79% of the rabbit autosomal genome. The majority of the detected ROHs ranged from 1 to 2 Mb in length, with an average ROH length (LROH) of 1.84 Mb. ROHs longer than 6 Mb constituted only 0.83% of the detected ROHs. The average inbreeding coefficient derived from ROHs (FROH) was 0.23, with FROH values ranging from 0.14 to 0.38 across breeds. Among Chinese indigenous breeds, the Jiuyishan rabbit exhibited the highest values of NROH, SROH, LROH, and FROH, whereas the Fujian Yellow rabbit had the lowest FROH values. In exotic rabbit breeds, the Japanese White rabbit displayed the highest values for NROH, SROH, LROH, and FROH, while the Flemish Giant rabbit had the lowest values for these metrics. Additionally, we identified 17 ROH islands in Chinese indigenous breeds and 22 ROH islands in exotic rabbit breeds, encompassing 124 and 186 genes, respectively. In Chinese indigenous breeds, we identified prominent genes associated with reproduction, including CFAP206, RNF133, CPNE4, ASTE1, and ATP2C1, as well as genes related to adaptation, namely CADPS2, FEZF1, and EPHA7. In contrast, the exotic breeds exhibited a prevalence of genes associated with fat deposition, such as ELOVL3 and NPM3, as well as growth and body weight related genes, including FAM184B, NSMCE2, and TWNK.

CONCLUSIONS

This study enhances our understanding of genetic diversity and selection pressures in domestic rabbits, offering valuable implications for breeding management and conservation strategies.

Revealing the Transcriptional and Metabolic Characteristics of Sebocytes Based on the Donkey Cell Transcriptome Atlas

Worldwide, donkeys (Equus asinus) are valued for their meat and milk, and in China also for the medical value of their skin. Physiological characteristics are key to the donkey's adaptability, including their digestive, respiratory, and reproductive systems, which enable them to survive and work in a variety of environments. However, the understanding of donkey physiological characteristics at the cellular level remains poor. Thus, single-cell transcriptome sequencing is used to construct a detailed transcriptional atlas based on 20 tissues from the Dezhou donkey (in total 84 cell types and 275 050 high quality cells) to perform an in-depth investigation of molecular physiology. Cross-species and cross-tissue comparative analyses reveal SOX10 to be an evolutionally conserved regulon in oligodendrocytes and illuminate the distinctive transcriptional patterns of donkey sebocytes. Moreover, through multispecies skin metabolomics, highly abundant, species-specific metabolites in donkey skin are identified, such as arachidonic acid and gamma-glutamylcysteine, and the pivotal role of sebocytes in donkey skin metabolism is highlighted. In summary, this work offers new insights into the unique metabolic patterns of donkey skin and provides a valuable resource for the conservation of donkey germplasm and the advancement of selective breeding programs.

Application of Omics in Donkey Meat Research: A Review

This review comprehensively examines the molecular basis of donkey meat characteristics and growth-associated genes, integrating findings from multiple omics approaches. This study examines the nutritional profile of donkey meat, which is rich in protein, essential amino acids, and unsaturated fatty acids. Through a systematic literature search across Web of Science, Google Scholar, PubMed, and Scopus databases (2000-2024), we collected and analyzed data from 400 research articles using predefined inclusion criteria focused on nutritional composition, omics approaches, and meat quality parameters in donkey populations. The study also evaluates various factors affecting meat quality, including breed differences, age, feeding management, and storage conditions. Advanced genomic and transcriptomic analyses have revealed numerous candidate genes, such as ACTN3, BMP7, NR6A1, Wnt7a, HOXC8, LCORL, TPM2, and TPM3, associated with growth traits and meat quality characteristics, providing valuable insights for genetic improvement programs. Furthermore, the review discusses various authentication methods for ensuring donkey meat quality and preventing adulteration, highlighting the integration of traditional and modern analytical approaches.

A review of genetic resources and trends of omics applications in donkey research: focus on China

Omics methodologies, such as genomics, transcriptomics, proteomics, metabolomics, lipidomics and microbiomics, have revolutionized biological research by allowing comprehensive molecular analysis in livestock animals. However, despite being widely used in various animal species, research on donkeys has been notably scarce. China, renowned for its rich history in donkey husbandry, plays a pivotal role in their conservation and utilization. China boasts 24 distinct donkey breeds, necessitating conservation efforts, especially for smaller breeds facing extinction threats. So far, omics approaches have been employed in studies of donkey milk and meat, shedding light on their composition and quality. Similarly, omics methods have been utilized to explore the molecular basis associated with donkey growth, meat production, and quality traits. Omics analysis has also unraveled the critical role of donkey microbiota in health and nutrition, with gut microbiome studies revealing associations with factors such as pregnancy, age, transportation stress, and altitude. Furthermore, omics applications have addressed donkey health issues, including infectious diseases and reproductive problems. In addition, these applications have also provided insights into the improvement of donkey reproductive efficiency research. In conclusion, omics methodologies are essential for advancing knowledge about donkeys, their genetic diversity, and their applications across various domains. However, omics research in donkeys is still in its infancy, and there is a need for continued research to enhance donkey breeding, production, and welfare in China and beyond.

Genome-wide association analysis identify candidate genes for feed efficiency and growth traits in Wenchang chickens

BACKGROUND

Wenchang chickens are one of the most popular local chicken breeds in the Chinese chicken industry. However, the low feed efficiency is the main shortcoming of this breed. Therefore, there is a need to find a more precise breeding method to improve the feed efficiency of Wenchang chickens. In this study, we explored important candidate genes and variants for feed efficiency and growth traits through genome-wide association study (GWAS) analysis.

RESULTS

Estimates of genomic heritability for growth and feed efficiency traits, including residual feed intake (RFI) of 0.05, average daily food intake (ADFI) of 0.21, average daily weight gain (ADG) of 0.24, body weight (BW) at 87, 95, 104, 113 days of age (BW87, BW95, BW104 and BW113) ranged from 0.30 to 0.44. Important candidate genes related to feed efficiency and growth traits were identified, such as PLCE1, LAP3, MED28, QDPR, LDB2 and SEL1L3 genes.

CONCLUSION

The results identified important candidate genes for feed efficiency and growth traits in Wenchang chickens and provide a theoretical basis for the development of new molecular breeding technology.

Genomic analyses of withers height and linear conformation traits in German Warmblood horses using imputed sequence-level genotypes

BACKGROUND

Body conformation, including withers height, is a major selection criterion in horse breeding and is associated with other important traits, such as health and performance. However, little is known about the genomic background of equine conformation. Therefore, the aim of this study was to use imputed sequence-level genotypes from up to 4891 German Warmblood horses to identify genomic regions associated with withers height and linear conformation traits. Furthermore, the traits were genetically characterised and putative causal variants for withers height were detected.

RESULTS

A genome-wide association study (GWAS) for withers height confirmed the presence of a previously known quantitative trait locus (QTL) on Equus caballus (ECA) chromosome 3 close to the LCORL/NCAPG locus, which explained 16% of the phenotypic variance for withers height. An additional significant association signal was detected on ECA1. Further investigations of the region on ECA3 identified a few promising candidate causal variants for withers height, including a nonsense mutation in the coding sequence of the LCORL gene. The estimated heritability for withers height was 0.53 and ranged from 0 to 0.34 for the conformation traits. GWAS identified significantly associated variants for more than half of the investigated conformation traits, among which 13 showed a peak on ECA3 in the same region as withers height. Genetic parameter estimation revealed high genetic correlations between these traits and withers height for the QTL on ECA3.

CONCLUSIONS

The use of imputed sequence-level genotypes from a large study cohort led to the discovery of novel QTL associated with conformation traits in German Warmblood horses. The results indicate the high relevance of the QTL on ECA3 for various conformation traits, including withers height, and contribute to deciphering causal mutations for body size in horses.

Genome-wide association study of growth traits and validation of key mutations (MSTN c.C861T) associated with the muscle mass of meat pigeons

Selective breeding of meat pigeons is primarily based on growth traits, especially muscle mass (MM). Identification of functional genes and molecular markers of growth and slaughter traits through a genome-wide association study (GWAS) will help to elucidate the underlying molecular mechanisms and provide a theoretical basis for the selective breeding of meat pigeons. The phenotypic data of body weight (BW) and body size (BS) of 556 meat pigeons at 52 and 80 weeks of age were collected. In total, 160 434 high-quality single nucleotide polymorphism sites were obtained by restriction site-associated DNA sequencing. The GWAS analysis revealed that MSTN, IGF2BP3 and NCAPG/LCORL were important candidate genes affecting the growth traits of meat pigeons. IGF2BP3 and NCAPG/LCORL were highly correlated to BW and BS, which are related to overall growth and development, while MSTN was associated with pectoral thickness and BW. Phenotypic association validation with the use of two meat pigeon populations found that the MSTN mutation c.C861T determines the MM. These results provide new insights into the genetic mechanisms underlying phenotypic variations of growth traits and MM in commercial meat pigeons. The identified markers and genes provide a theoretical basis for the selective breeding of meat pigeons.

A mutation in POLR2A gene associated with body size traits in Dezhou donkeys revealed with GWAS

The Dezhou donkey is a famous local donkey breed in China. The aim of the present study was to identify the genes associated with the body size traits of the Dezhou donkey and facilitate the breeding activities of the donkeys. A total of 349 donkeys from 2 generations (113 individuals in F0 and 236 in F1) were analyzed with restriction-site-associated DNA sequencing. A genome-wide association study revealed that the region between 13.7 and 15.6 Mb of chromosome 13 is significantly associated with body sizes. Candidate genes related to body size development, including POLR2A, CHRNB1, FGF11, and ZBTB4, were identified. The results of GO and KEGG analysis indicated that the genes involved in many GO terms were related to metabolic processes and developmental processes. Additionally, a T>C mutation (Chr13:14312485) was found at intron 10 of the POLR2A gene. The association analysis showed significant differences among genotypes for the size traits. The body size of the individuals with the TT genotype was significantly higher than that with the CC genotype. The results showed that the polymorphism of POLR2A has the potential to be used as a marker in the breeding programs of the Dezhou donkeys.

Genomic Selection for Live Weight in the 14th Month in Alpine Merino Sheep Combining GWAS Information

Alpine Merino Sheep is a novel breed reared from Australian Merino Sheep as the father and Gansu Alpine Fine-Wool Sheep as the mother, living all year in cold and arid alpine areas with exceptional wool quality and meat performance. Body weight is an important economic trait of the Alpine Merino Sheep, but there is limited research on identifying the genes associated with live weight in the 14th month for improving the accuracy of the genomic prediction of this trait. Therefore, this study's sample comprised 1310 Alpine Merino Sheep ewes, and the Fine Wool Sheep 50K Panel was used for genome-wide association study (GWAS) analysis to identify candidate genes. Moreover, the trial population (1310 ewes) in this study was randomly divided into two groups. One group was used as the population for GWAS analysis and screened for the most significant top 5%, top 10%, top 15%, and top 20% SNPs to obtain prior marker information. The other group was used to estimate the genetic parameters based on the weight assigned by heritability combined with different prior marker information. The aim of this study was to compare the accuracy of genomic breeding value estimation when combined with prior marker information from GWAS analysis with the optimal linear unbiased prediction method for genome selection (GBLUP) for the breeding value of target traits. Finally, the accuracy was evaluated using the five-fold cross-validation method. This research provides theoretical and technical support to improve the accuracy of sheep genome selection and better guide breeding. The results demonstrated that eight candidate genes were associated with GWAS analysis, and the gene function query and literature search results suggested that FAM184B, NCAPG, MACF1, ANKRD44, DCAF16, FUK, LCORL, and SYN3 were candidate genes affecting live weight in the 14th month (WT), which regulated the growth of muscle and bone in sheep. In genome selection analysis, the heritability of GBLUP to calculate the WT was 0.335-0.374, the accuracy after five-fold cross-verification was 0.154-0.190, and after assigning different weights to the top 5%, top 10%, top 15%, and top 20% of the GWAS results in accordance with previous information to construct the G matrix, the accuracy of the WT in the GBLUP model was improved by 2.59-7.79%.