A genome-wide association study revealed five SNPs affecting 8-month weight in sheep

Genome-wide association studies of body size traits in Tibetan sheep

BACKGROUND

Elucidating the genetic variation underlying phenotypic diversity will facilitate improving production performance in livestock species. The Tibetan sheep breed in China holds significant historical importance, serving as a fundamental pillar of Qinghai's animal husbandry sector. The Plateau-type Tibetan sheep, comprising 90% of the province's population, are characterized by their tall stature and serve as the primary breed among Tibetan sheep. In contrast, Zhashijia sheep exhibit larger size and superior meat quality. These two species provide an excellent model for elucidating the genetic basis of body size variation. Therefore, this study aims to conduct a comprehensive genome-wide association study on these two Tibetan sheep breeds to identify single nucleotide polymorphism loci and regulatory genes that influence body size traits in Tibetan sheep.

RESULT

In this study, the phenotypic traits of body weight, body length, body height, chest circumference, chest depth, chest width, waist angle width, and pipe circumference were evaluated in two Tibetan sheep breeds: Plateau-type sheep and Zhashijia Tibetan sheep. Whole genome sequencing generated 48,215,130 high-quality SNPs for genome-wide association study. Four methods were applied and identified 623 SNPs significantly associated with body size traits. The significantly associated single nucleotide polymorphisms identified in this study are located near or within 111 candidate genes. These genes exhibit enrichment in the cAMP and Rap1 signaling pathways, significantly affecting animal growth, and body size. Specifically, the following genes were associated: ASAP1, CDK6, FRYL, NAV2, PTPRM, GPC6, PTPRG, KANK1, NTRK2 and ADCY8.

CONCLUSION

By genome-wide association study, we identified 16 SNPs and 10 candidate genes associated with body size traits in Tibetan sheep, which hold potential for application in genomic selection breeding programs in sheep. Identifying these candidate genes will establish a solid foundation for applying molecular marker-assisted selection in sheep breeding and improve our understanding of body size control in farmed animals.

NCAPG Regulates Myogenesis in Sheep, and SNPs Located in Its Putative Promoter Region Are Associated with Growth and Development Traits

Previously, NCAPG was identified as a candidate gene associated with sheep growth traits. This study aimed to investigate the direct role of NCAPG in regulating myogenesis in embryonic myoblast cells and to investigate the association between single-nucleotide polymorphisms (SNPs) in its promoter region and sheep growth traits. The function of NCAPG in myoblast proliferation and differentiation was detected after small interfering RNAs (siRNAs) knocked down the expression of NCAPG. Cell proliferation was detected using CCK-8 assay, EdU proliferation assay, and flow cytometry cell cycle analysis. Cell differentiation was detected via cell immunofluorescence and the quantification of myogenic regulatory factors (MRFs). SNPs in the promoter region were detected using Sanger sequencing and genotyped using the improved multiplex ligation detection reaction (iMLDR®) technique. As a result, a notable decrease (p < 0.01) in the percentage of EdU-positive cells in the siRNA-694-treated group was observed. A significant decrease (p < 0.01) in cell viability after treatment with siRNA-694 for 48 h and 72 h was detected using the CCK-8 method. The quantity of S-phase cells in the siRNA-694 treatment group was significantly decreased (p < 0.01). After interfering with NCAPG in myoblasts during induced differentiation, the relative expression levels of MRFs were markedly (p < 0.05 or p < 0.01) reduced compared with the control group on days 5-7. The myoblast differentiation in the siRNA-694 treatment group was obviously suppressed compared with the control group. SNP1, SNP2, SNP3, and SNP4 were significantly (p < 0.05) associated with all traits except body weight measured at birth and one month of age. SNP5 was significantly (p < 0.05) associated with body weight, body height, and body length in six-month-old sheep. In conclusion, interfering with NCAPG can inhibit the proliferation and differentiation of ovine embryonic myoblasts. SNPs in its promoter region can serve as potential useful markers for selecting sheep growth traits.

Uncovering the candidate genes related to sheep body weight using multi-trait genome-wide association analysis

In sheep, body weight is an economically important trait. This study sought to map genetic loci related to weaning weight and yearling weight. To this end, a single-trait and multi-trait genome-wide association study (GWAS) was performed using a high-density 600 K single nucleotide polymorphism (SNP) chip. The results showed that 43 and 56 SNPs were significantly associated with weaning weight and yearling weight, respectively. A region associated with both weaning and yearling traits (OARX: 6.74-7.04 Mb) was identified, suggesting that the same genes could play a role in regulating both these traits. This region was found to contain three genes (TBL1X, SHROOM2 and GPR143). The most significant SNP was Affx-281066395, located at 6.94 Mb (p = 1.70 × 10-17), corresponding to the SHROOM2 gene. We also identified 93 novel SNPs elated to sheep weight using multi-trait GWAS analysis. A new genomic region (OAR10: 76.04-77.23 Mb) with 22 significant SNPs were discovered. Combining transcriptomic data from multiple tissues and genomic data in sheep, we found the HINT1, ASB11 and GPR143 genes may involve in sheep body weight. So, multi-omic anlaysis is a valuable strategy identifying candidate genes related to body weight.

Genome-Wide Association Analysis of Muscle pH in Texel Sheep × Altay Sheep F2 Resource Population

pH was one of the important meat quality traits, which was an important factor affecting the storage/shelf life and quality of meat in meat production. In order to find a way to extend the storage/shelf life, the pH values (pH_45min, pH_24h, pH_48h and pH_72h) of the longissimus dorsi muscles in F₂ individuals of 462 Texel sheep × Altay sheep were determined, genotyping was performed using Illumina Ovine SNP 600 K BeadChip and whole genome resequencing technology, a genome-wide association analysis (GWAS) was used to screen the candidate genes and molecular markers for pH values related to the quality traits of mutton, and the effects of population stratification were detected by Q-Q plots. The results showed that the pH population stratification analysis did not find significant systemic bias, and there was no obvious population stratification effect. The results of the association analysis showed that 28 SNPs significantly associated with pH reached the level of genomic significance. The candidate gene associated with pH_45min was identified as the CCDC92 gene by gene annotation and a search of the literature. Candidate genes related to pH_24h were KDM4C, TGFB2 and GOT2 genes. The candidate genes related to pH_48h were MMP12 and MMP13 genes. The candidate genes related to pH_72h were HILPDA and FAT1 genes. Further bioinformatics analyses showed 24 gene ontology terms and five signaling pathways that were significantly enriched (p ≤ 0.05). Many terms and pathways were related to cellular components, processes of protein modification, the activity of protein dimerization and hydrolase activity. These identified SNPs and genes could provide useful information about meat and the storage/shelf life of meat, thereby extending the storage/shelf life and quality of meat.

Assessing Genomic Diversity and Signatures of Selection in Chinese Red Steppe Cattle Using High-Density SNP Array

Chinese Red Steppe Cattle (CRS), a composite cattle breed, is well known for its milk production, high slaughter rate, carcass traits, and meat quality. Nowadays, it is widely bred in Jilin and Hebei Province and the Inner Mongolia Autonomous region. However, the population structure and the genetic basis of prominent characteristics of CRS are still unknown. In this study, we systematically describe their population structure, genetic diversity, and selection signature based on genotyping data from 61 CRS individuals with GGP Bovine 100 K chip. The results showed that CRS cattle had low inbreeding levels and had formed a unique genetic structure feature. Using two complementary methods (including comprehensive haplotype score and complex likelihood ratio), we identified 1291 and 1285 potentially selected genes, respectively. There were 141 genes annotated in common 106 overlapping genomic regions covered 5.62 Mb, including PLAG1, PRKG2, DGAT1, PARP10, TONSL, ADCK5, and BMP3, most of which were enriched in pathways related to muscle growth and differentiation, milk production, and lipid metabolism. This study will contribute to understanding the genetic mechanism behind artificial selection and give an extensive reference for subsequent breeding.

Genome-wide identification of SNPs associated with body weight in yak

BACKGROUND

The yak is the most important livestock in the Qinghai-Tibet Plateau, and body weight directly affects the economic values of yak. Up to date, the genome-wide profiling of single-nucleotide polymorphisms (SNPs) associating with body weight has not been reported in yak. In the present study, the SNPs in 480 yaks from three breeds were analyzed using the commercial high-density (600 K) yak SNP chips.

RESULTS

The results identified 12 and 4 SNPs potentially associated with body weight in male and female yaks, respectively. Among them, 9 and 2 SNPs showed significant difference in yak body weight between different genotypes at each locus in male and female yaks, respectively. Further exploration found 33 coding genes within the 100 kbp upstream or downstream to the SNP loci, which might be potentially affected by the variation of SNPs. Among them, G protein-coupled receptor kinase 4 (GRK4) might be potentially affected by the SNP AX-174555047, which has been reported to affect the functioning of two body-weight associated hormones (parathyroid hormone, PTH, and adrenomedullin, ADM). Determination of PTH and ADM levels in yak revealed positive relationship between PTH level and body weight, negative relationship between ADM level and body weight along with the variation of AX-174555047 mutation.

CONCLUSIONS

These results suggested that the SNP AX-174555047 might potentially affect body weight through mediating GRK4 expression and then PTH and ADM functioning. Thus, the SNP AX-174555047 might be used as a biomarker for molecular breeding of yak. More investigations are required to validate this point.

Sheep fecal transplantation affects growth performance in mouse models by altering gut microbiota

Animal growth traits are important and complex traits that determine the productivity of animal husbandry. There are many factors that affect growth traits, among which diet digestion is the key factor. In the process of animal digestion and absorption, the role of gastrointestinal microbes is essential. In this study, we transplanted two groups of sheep intestinal microorganisms with different body weights into the intestines of mice of the same age to observe the effect of fecal bacteria transplantation on the growth characteristics of the mouse model. The results showed that receiving fecal microbiota transplantation (FMT) had an effect on the growth traits of recipient mice (P < 0.05). Interestingly, only mice receiving high-weight donor microorganisms showed differences. Use 16S rDNA sequencing technology to analyze the stool microorganisms of sheep and mice. The microbial analysis of mouse feces showed that receiving FMT could improve the diversity and richness of microorganisms (P < 0.05), and the microbial composition of mouse feces receiving low-weight donor microorganisms was similar to that of the control group, which was consistent with the change trend of growth traits. The feces of high-weight sheep may have higher colonization ability. The same five biomarkers were identified in the donor and recipient, all belonging to Firmicutes, and were positively correlated with the body weight of mice at each stage. These results suggest that FMT affects the growth traits of receptors by remodeling their gut microflora.

Identification of Genomic Regions and Candidate Genes Associated with Body Weight and Body Conformation Traits in Karachai Goats

The objective of this study was to identify the SNPs and candidate genes related to body weight and seven body conformation traits at the age of 8 months in the Russian aboriginal Karachai goats (n = 269) by conducting genome-wide association studies (GWAS), using genotypes generated by Goat SNP BeadChip (Illumina Inc., USA). We identified 241 SNPs, which were significantly associated with the studied traits, including 47 genome-wide SNPs (p < 10−5) and 194 suggestive SNPs (p < 10−4), distributed among all goat autosomes except for autosome 23. Fifty-six SNPs were common for two and more traits (1 SNP for six traits, 2 SNPs for five traits, 12 SNPs for four traits, 20 SNPs for three traits, and 21 SNPs for two traits), while 185 SNPs were associated with single traits. Structural annotation within a window of 0.4 Mb (±0.2 Mb from causal SNPs) revealed 238 candidate genes. The largest number of candidate genes was identified at Chr13 (33 candidate genes for the five traits). The genes identified in our study were previously reported to be associated with growth-related traits in different livestock species. The most significant genes for body weight were CRADD, HMGA2, MSRB3, MAX, HACL1 and RAB15, which regulate growth processes, body sizes, fat deposition, and average daily gains. Among them, the HMGA2 gene is a well-known candidate for prenatal and early postnatal development, and the MSRB3 gene is proposed as a candidate gene affecting the growth performance. APOB, PTPRK, BCAR1, AOAH and ASAH1 genes associated with withers height, rump height and body length, are involved in various metabolic processes, including fatty acid metabolism and lipopolysaccharide catabolism. In addition, WDR70, ZBTB24, ADIPOQ, and SORCS3 genes were linked to chest width. KCNG4 was associated with rump height, body length and chest perimeter. The identified candidate genes can be proposed as molecular markers for growth trait selection for genetic improvement in Karachai goats.

New Insight into the Role of the Leucine Aminopeptidase 3 (LAP3) in Cell Proliferation and Myogenic Differentiation in Sheep Embryonic Myoblasts

Previous genome-wide association studies (GWAS) have found that LAP3 may have the potential function to impact sheep muscle development. In order to further explore whether LAP3 expression has an important role in the development of sheep embryonic myoblasts, we conducted the spatiotemporal expression profile analysis of LAP3 at the tissue and cellular level. Then we used small interfering RNA and eukaryotic recombinant vectors to perform gain/loss-of-function analysis of LAP3. CCK-8 detection, EdU staining, and flow cytometry were used to investigate the impact of LAP3 knockdown or overexpression on the proliferation of embryonic myoblasts. In addition, cell phenotype observation, MyHC indirect immunofluorescence, and quantitative detection of the expression changes of myogenic regulatory factors (MRFs) were used to explore the effect of LAP3 on myogenic differentiation. The results showed that the LAP3 expression level in muscle tissue of fetuses was significantly higher than that in newborn lambs and adult sheep, and its expression level on day 3 of differentiation was also significantly higher than that in the proliferation phase and other differentiation time points. LAP3 silencing could significantly increase cell viability and EdU-positive cells, as well as prolonging the length of S phase of myoblasts to promote proliferation, while the results were reversed when LAP3 was overexpressed. Moreover, LAP3 silencing significantly hindered myotube formation and down-regulated the expression levels of MRFs from day 5 to day 7 of terminal differentiation, while the results were reversed when LAP3 was highly expressed. Overall, our results suggested that the expression of LAP3 impacts on the development of sheep embryonic myoblasts which provides an important theoretical basis for molecular breeding of meat production in sheep.

Whole genome sequencing of Luxi Black Head sheep for screening selection signatures associated with important traits

Objective

Luxi Black Head sheep (LBH) is the first crossbreed specialized for meat production and was developed by crossbreeding Black Head Dorper sheep (DP) and Small Tailed Han sheep (STH) in the farming areas of northern China. Research on the genomic variations and selection signatures of LBH caused by continuous artificial selection is of great significance for identifying the genetic mechanisms of important traits of sheep and for the continuous breeding of LBH.

Methods

We explored the genetic relationships of LBH, DP, and several Mongolian sheep breeds by constructing phylogenetic tree, principal component analysis and linkage disequilibrium analysis. In addition, we analysed 29 whole genomes of sheep. The genome-wide selection signatures have been scanned with four methods: heterozygosity (H_P), fixation index (F_ST), cross-population extended haplotype homozygosity (XP-EHH) and the nucleotide diversity (θ_π) ratio.

Results

The genetic relationships analysis showed that LBH appeared to be an independent cluster closer to DP. The candidate signatures of positive selection in sheep genome revealed candidate genes for developmental process (HoxA gene cluster, BCL2L11, TSHR), immunity (CXCL6, CXCL1, SKAP2, PTK6, MST1R), growth (PDGFD, FGF18, SRF, SOCS2), and reproduction (BCAS3, TRIM24, ASTL, FNDC3A). Moreover, two signalling pathways closely related to reproduction, the thyroid hormone signalling pathway and the oxytocin signalling pathway, were detected.

Conclusion

The selective sweep analysis of LBH genome revealed candidate genes and signalling pathways associated with developmental process, immunity, growth, and reproduction. Our findings provide a valuable resource for sheep breeding and insight into the mechanisms of artificial selection.

Genome-wide association study and inbreeding depression on body size traits in Qira black sheep (Ovis aries)

Qira black sheep is a famous indigenous sheep breed in China. The objectives of this study are to identify candidate genes related to body size, and to estimate the level of inbreeding depression on body size based on runs of homozygosity in Qira black sheep. Here, 188 adult Qira black sheep were genotyped with a high density (630 K) SNP chip and genome-wide association study for body weight and body size traits (including withers height, body slanting length, tail length, chest girth, chest width, and chest depth) were performed using an additive linear model. In consequence, 12 genome- and chromosome-wide significant SNPs and, accordingly, six candidate genes involved in muscle differentiation, metabolism and cell processes were identified. Of them, ZNF704 (zinc finger protein 704) was identified for body weight; AK2 (adenylate kinase 2) and PARK2 (parkin RBR E3 ubiquitin protein ligase) for tail length; MOCOS (molybdenum cofactor sulfurase) and ELP2 (elongator acetyltransferase complex subunit 2) for chest width; and MFAP1 (microfibril associated protein 1) for chest girth. Additionally, inbreeding depressions on body size were observed in the current herd. These results will provide insightful understandings into the genetic mechanisms of adult body size, and into the conservation and utilization of Qira black sheep.