Screening of Genes Related to Growth, Development and Meat Quality of Sahan Crossbred F1 Sheep Based on RNA-Seq Technology

Estimation of Genetic Parameters of Early Growth Traits in Dumeng Sheep

This study aimed to estimate the genetic and non-genetic factors that affect the nine early growth traits of Dumeng sheep, as well as to estimate the variance components and genetic parameters associated with these traits. A dataset containing detailed information on 17,896 preweaning trait records of 4474 lambs was collected. In addition, 5015 postweaning trait records of 1003 lambs were documented. The effects of recipient dam age, sex, year, season, and herd on the early growth traits were assessed using the general linear model procedure of the statistical analysis system, revealing different levels of significance across different traits. To determine the most suitable model for estimating the genetic parameters, the likelihood ratio (LR) test was employed, fitting six animal models that either excluded or included maternal genetic and maternal permanent environmental effects within the average information restricted maximum likelihood (AIREML) framework using WOMBAT software (Version: 23/11/23). The model incorporating direct additive genetic effects, maternal genetic effects, and maternal permanent environment effects as random effects (model 6) provided the best fit for birth weight (BW) estimation. In contrast, the model combining direct additive genetic effects and maternal permanent environment effects as random effects (model 2) demonstrated a superior fit for estimating the genetic parameters of weaning weight (WW), average daily gain weight from birth to weaning (ADG1), and Kleiber ratio from birth to weaning (KR1). With regard to the genetic parameters of body weight at 6 months of age (6MW), average daily gain weight from weaning to 6 months (ADG2), average daily gain weight from birth to 6 months (ADG3), Kleiber ratio from weaning to 6 months (KR2), and Kleiber ratio from birth to 6 months (KR3), model 1, which incorporates only direct additive genetic effects, was identified as the optimal choice. With the optimal model, the heritability estimates ranged from 0.010 ± 0.033 for 6MW to 0.1837 ± 0.096 for KR3. The bivariate analysis method was employed to estimate the correlation between various traits using the most suitable model. The absolute values of genetic correlation coefficients among the traits spanned a range from 0.1460 to 0.9998, highlighting both weak and strong relationships among the studied traits. Specifically, the estimated genetic correlations between WW and ADG1, ADG3, KR1, and KR3 were 0.9859, 0.9953, 0.9911, and 0.9951, respectively, while the corresponding phenotypic correlations were 0.9752, 0.7836, 0.8262, and 0.5767. These findings identified that WW could serve as an effective selection criterion for enhancing early growth traits.

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.

Genomic structure analysis and construction of DNA fingerprint for four sheep populations

The Yongdeng Qishan sheep (QS) is a sheep population found locally in China. To gain in-depth knowledge of its population characteristics, three control groups were chosen, comprising the Lanzhou fat-tailed sheep (LFT), TAN sheep (TAN), and Minxian black fur sheep (MBF), inhabiting the nearby environments. This study genotyped a total of 120 individuals from four sheep populations: QS, LFT, TAN, and MBF. Using Specific-Locus Amplified Fragment Sequencing, we conducted genetic diversity, population structure, and selective sweep analysis, and constructed the fingerprint of each population. In total, there were 782 535 single nucleotide polymorphism (SNP) variations identified, with most being situated within regions that are intergenic or intronic. The genetic diversity analysis revealed that the QS population exhibited lower genetic diversity compared to the other three populations. Consistent results were obtained from the principal component, phylogenetic tree, and population structure analysis, indicating significant genetic differences between QS and the other three populations. However, a certain degree of differentiation was observed within the QS population. The linkage disequilibrium (LD) patterns among the four populations showed clear distinctions, with the QS group demonstrating the most rapid LD decline. Kinship analysis supported the findings of population structure, dividing the 90 QS individuals into two subgroups consisting of 23 and 67 individuals. Selective sweep analysis identified a range of genes associated with reproduction, immunity, and adaptation to high-altitude hypoxia. These genes hold potential as candidate genes for marker-assisted selection breeding. Additionally, a total of 86 523 runs of homozygosity (ROHs) were detected, showing non-uniform distribution across chromosomes, with chromosome 1 having the highest coverage percentage and chromosome 26 the lowest. In the high-frequency ROH islands, 79 candidate genes were associated with biological processes such as reproduction and fat digestion and absorption. Furthermore, a DNA fingerprint was constructed for the four populations using 349 highly polymorphic SNPs. In summary, our research delves into the genetic diversity and population structure of QS population. The construction of DNA fingerprint profiles for each population can provide valuable references for the identification of sheep breeds both domestically and internationally.

Whole-Genome Scanning for Selection Signatures Reveals Candidate Genes Associated with Growth and Tail Length in Sheep

Compared to Chinese indigenous sheep, Western sheep have rapid growth rate, larger physique, and higher meat yield. These excellent Western sheep were introduced into China for crossbreeding to expedite the enhancement of production performance and mutton quality in local breeds. Here, we investigated population genetic structure and genome-wide selection signatures among the Chinese indigenous sheep and the introduced sheep based on whole-genome resequencing data. The PCA, N-J tree and ADMIXTURE results showed significant genetic difference between Chinese indigenous sheep and introduced sheep. The nucleotide diversity (π) and linkage disequilibrium (LD) decay results indicated that the genomic diversity of introduced breeds were lower. Then, Fst & π ratio, XP-EHH, and de-correlated composite of multiple signals (DCMS) methods were used to detect the selection signals. The results showed that we identified important candidate genes related to growth rate and body size in the introduced breeds. Selected genes with stronger selection signatures are associated with growth rate (CRADD), embryonic development (BVES, LIN28B, and WNT11), body size (HMGA2, MSRB3, and PTCH1), muscle development and fat metabolism (MSTN, PDE3A, LGALS12, GGPS1, and SAR1B), wool color (ASIP), and hair development (KRT71, KRT74, and IRF2BP2). Thus, these genes have the potential to serve as candidate genes for enhancing the growth traits of Chinese indigenous sheep. We also identified tail-length trait-related candidate genes (HOXB13, LIN28A, PAX3, and VEGFA) in Chinese long-tailed breeds. Among these genes, HOXB13 is the main candidate gene for sheep tail length phenotype. LIN28A, PAX3, and VEGFA are related to embryonic development and angiogenesis, so these genes may be candidate genes for sheep tail type traits. This study will serve as a foundation for further genetic improvement of Chinese indigenous sheep and as a reference for studies related to growth and development of sheep.

Genome-Wide Identification and Characterization of Long Non-Coding RNAs in Longissimus dorsi Skeletal Muscle of Shandong Black Cattle and Luxi Cattle

There is an increasing understanding of the possible regulatory role of long non-coding RNAs (LncRNA). Studies on livestock have mainly focused on the regulation of cell differentiation, fat synthesis, and embryonic development. However, there has been little study of skeletal muscle of domestic animals and the potential role of lncRNA. In this study, the transcriptome numbers of longissimus muscle of different beef cattle (Shandong black catle and Luxi catlle) were used to construct muscle related lncRNAs-miRNA-mRNA interaction network through bioinformatics analysis. This is helpful to clarify the molecular mechanism of bovine muscle development, and can be used to promote animal husbandry and improve animal husbandry production. According to the screening criteria of |FC|≧2 and q < 0.05, a total of 1,415 transcripts (of which 480 were LncRNAs) were differentially expressed (q < 0.05) in the different breeds. Further, we found that the most differentially expressed LncRNAs were found on chromosome 9, in which the differentially expressed LncRNAs targeted 1,164 protein coding genes (MYORG, Wnt4, PAK1, ADCY7,etc) (upstream and downstream<50 Kb). In addition, Pearson’s correlation coefficients of co-expression levels indicated a potential trans regulatory relationship between the differentially expressed LncRNAs and 43844 mRNAs (r > 0.9). The identified co-expressed mRNAs (MYORG, Dll1, EFNB2, SOX6, MYOCD, and MYLK3) are related to the formation of muscle structure, and enriched in muscle system process, strained muscle cell differentiation, muscle cell development, striated muscle tissue development, calcium signaling, and AMPK signaling. Additionally, we also found that some LncRNAs (LOC112444238, LOC101903367, LOC104975788, LOC112441863, LOC112449549, and LOC101907194) may interact with miRNAs related to cattle muscle growth and development. Based on this, we constructed a LncRNAs-miRNA-mRNA interaction network as the putative basis for biological regulation in cattle skeletal muscle. Interestingly, a candidate differential LncRNA (LOC104975788) and a protein-coding gene (Pax7) contain miR-133a binding sites and binding was confirmed by luciferase reporter assay. LOC104975788 may combined miR-133a competitively with Pax7, thus relieving the inhibitory effect of miR-133a on Pax7 to regulate skeletal muscle development. These results will provide the theoretical basis for further study of LncRNA regulation and activity in different cattle breeds.