A sheep pangenome reveals the spectrum of structural variations and their effects on tail phenotypes.

Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here we generated high-quality genome assemblies for 15 individuals from genetically diverse sheep breeds using Pacific Biosciences (PacBio) high-fidelity sequencing, discovering 130.3 Mb nonreference sequences, from which 588 genes were annotated. A total of 149,158 biallelic insertions/deletions, 6531 divergent alleles, and 14,707 multiallelic variations with precise breakpoints were discovered. The SV spectrum is characterized by an excess of derived insertions compared to deletions (94,422 vs. 33,571), suggesting recent active LINE expansions in sheep. Nearly half of the SVs display low to moderate linkage disequilibrium with surrounding single-nucleotide polymorphisms (SNPs) and most SVs cannot be tagged by SNP probes from the widely used ovine 50K SNP chip. We identified 865 population-stratified SVs including 122 SVs possibly derived in the domestication process among 690 individuals from sheep breeds worldwide. A novel 168-bp insertion in the 5' untranslated region (5' UTR) of HOXB13 is found at high frequency in long-tailed sheep. Further genome-wide association study and gene expression analyses suggest that this mutation is causative for the long-tail trait. In summary, we have developed a panel of high-quality de novo assemblies and present a catalog of structural variations in sheep. Our data capture abundant candidate functional variations that were previously unexplored and provide a fundamental resource for understanding trait biology in sheep.

Whole genome sequencing of simmental cattle for SNP and CNV discovery.

BACKGROUD: The single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) are two major genomic variants, which play crucial roles in evolutionary and phenotypic diversity. RESULTS: In this study, we performed a comprehensive analysis to explore the genetic variations (SNPs and CNVs) of high sperm motility (HSM) and poor sperm motility (PSM) Simmental bulls using the high-coverage (25×) short-read next generation sequencing and single-molecule long reads sequencing data. A total of ~ 15 million SNPs and 2,944 CNV regions (CNVRs) were detected in Simmental bulls, and a set of positive selected genes (PSGs) and CNVRs were found to be overlapped with quantitative trait loci (QTLs) involving immunity, muscle development, reproduction, etc. In addition, we detected two new variants in LEPR, which may be related to the artificial breeding to improve important economic traits. Moreover, a set of genes and pathways functionally related to male fertility were identified. Remarkably, a CNV on SPAG16 (chr2:101,427,468 - 101,429,883) was completely deleted in all poor sperm motility (PSM) bulls and half of the bulls in high sperm motility (HSM), which may play a crucial role in the bull-fertility. CONCLUSIONS: In conclusion, this study provides a valuable genetic variation resource for the cattle breeding and selection programs.

Association of SNPs within PTPN3 gene with wool production and growth traits in a dual-purpose sheep population.

Protein tyrosine phosphatase non-receptor type 3 (PTPN3), a member of the membrane-associated non-receptor protein tyrosine phosphatase (PTP) family, plays significant roles in the cytoplasm and affects the development and growth of skin and hair. A recent study identified the PTPN3 as the potential gene related to sheep wool quality. To detect single-nucleotide polymorphisms (SNPs) of PTPN3 and elucidate its association with wool production and growth traits in fine wool sheep a total of 644 healthy SG (South African mutton merino♂ × Gansu alpine fine-wool sheep♀, SG) and SSG (South African mutton merino♂ × SG♀, SSG) hybrid sheep were selected. Pooled-DNA sequencing and SNPscan methods were used to scan and genotype SNPs within PTPN3. Association analyses between SNPs and wool production and growth traits were implemented. Consequently, the results revealed that PTPN3 has six SNPs (two missense mutations, one synonymous mutation, and three intron mutations), of which four loci (SNP2, SNP3, SNP4, and SNP5) were significantly positively correlated with growth and wool traits (p < 0.05). SNP4 was significantly (p < 0.05) linked with thigh wool length, and SNP6 was significantly (p < 0.05) associated with abdomen wool length. Moreover, one strongly linked SNP block was identified to be correlated with wool production and growth traits (body weight and body size). The significant SNPs founded by this study could serve as useful genetic markers for breeding fine-wool sheep.®.

Variants of ADPGK gene and its effect on the male reproductive organ parameters and sperm count in Hu sheep.

ADP-dependent glucokinase (ADPGK) plays an important role instead of hexokinase in regulating energy metabolism via the Embden-Meyerhof-Parnas Pathway. And energy provided via glycolysis promotes testis development and spermatogenesis. In this study, 466 Hu sheep were screened for mutations in the ADPGK gene to examine the association of the ADPGK gene polymorphisms with the testis traits and spermatogenesis. The NC_056060.1: g.31295 C > T SNP was found in the 3'-UTR region, resulting in two genotypes CC and TC type with genotypic frequencies of 0.66 and 0.34, respectively. This mutation was significantly associated with testis weight, testis long circumference, testis short girth, epididymis weight, and sperm concentration (p < 0.05). Moreover, TC genotype individuals had an increased tendency in the expression of the ADPGK gene and had significant reproductive performance advantages compared with CC genotype individuals in the study. And compared with the small testes (<50 g), the ADPGK gene expression of big testes (>160 g) increased significantly. This indicates an association between the ADPGK gene and reproductive organ parameters and sperm count in selected Hu sheep breed, and this SNP may serve as an effective DNA molecular marker for marker-assisted selection in Hu sheep breeding programs.

Genomic analysis reveals a KIT-related chromosomal translocation associated with the white coat phenotype in yak.

White coat pigmentation is a striking phenotype of many domesticated species and has various genetic controls. The Tianzhu White yak, an indigenous breed with a complete white coat, has fascinated Tibetans for centuries. However, the genetic basis of this trait remains unknown. Here, we conducted population genomics analysis and genome-wide association study based on the whole-genome sequencing data of 38 white and 59 non-white-coated yak. The results revealed the presence of KIT-linked Cs alleles characterized by the translocations between chromosomes 6 and 29 in all-white yak. Furthermore, structural variations showed additional duplications of the Cs alleles in white yak compared with colour-sidedness cattle. Interestingly, the Cs alleles associated with the white coat phenotype in yak were found to have introgressed from taurine cattle. Our findings unveil the shared genetic control of the white coat phenotype and its evolution in closely related bovine species.

A mutation modulating DDX3Y gene expression cosegregates with the major Y-chromosomal haplogroups and with testis size in Hu sheep.

Variations in the Y-chromosome are usually correlated with male-specific traits. However, this condition has been described only sporadically, even in human genetics. The present study was conducted to clone the full-length gene sequence of ovine DEAD-box helicase 3, Y-linked (DDX3Y), and investigate the effect of the expression and variation of DDX3Y on the reproductive traits of Hu sheep. Consequently, we identified the full coding sequence and genomic sequence of ovine DDX3Y. Quantitative PCR (qPCR) analysis showed that ovine DDX3Y was highly expressed in testis, and the expression level increased during testicular development. Furthermore, individuals with larger testis at 6 months expressed significantly more DDX3Y mRNA in the testis than individuals with smaller testis. Notably, a novel SNP (g. 12657 C>A) in the 3' untranslated region was identified in Hu sheep and Tan sheep according to the investigation of the full DDX3Y genomic sequence of 1069 individuals from nine sheep breeds. Association analysis revealed that the SNP was significantly related to testis size in Hu sheep. Meanwhile, Hu rams with the derived C allele showed significantly higher expression levels of DDX3Y in testis than those with the ancestral A allele. In addition, data mining in a previous study showed that the C allele cosegregated with the globally major Y-chromosomal haplogroups y-HA and y-HC, and the A allele is found in all rams with haplogroups y-HB1, y-HB2 and y-HD. This study suggests that the association of the Y-chromosomal haplogroups with testis size in Hu sheep can be extrapolated to the sheep population worldwide.

α-Linolenic acid induced TM4 Sertoli cells proliferation and enhanced total antioxidant capacity.

The present study was undertaken to elucidate the direct ALA effects on mice TM4 Sertoli cells proliferation in vitro. Our results showed that TM4 cells viability was significantly stimulated by ALA (p < 0.05). The 50 µM ALA increased the concentration of total antioxidant capacity, induced the mitochondrial membrane hyperpolarized, and markedly decreased the number of apoptosis cells (p < 0.05). ALA also up-regulated G2/Mitotic-specific cyclin-B1 gene and apoptosis suppressive gene Bcl2 expression (p < 0.05). In conclusion, those results indicated that ALA could increase TM4 Sertoli cells antioxidant capacity, induced the mitochondrial membrane hyperpolarized, inhibited cells apoptosis and stimulated TM4 Sertoli cells proliferation in vitro.

Expression profile of FASN gene and association of its polymorphisms with intramuscular fat content in Hu sheep.

The content of intramuscular fat (IMF) is one of the most important factors that has a large impact on meat quality, and it is an effective way to improve IMF according to marker-assisted selection (MAS). Fatty-acid synthase (FASN) is a key gene in meat lipid deposition and fatty acid composition. Thus, this study was conducted to investigate the expression profile of FASN in mRNA and protein levels using real-time quantitative PCR (RT-qPCR) and western-blot methods. In addition, single nucleotide polymorphisms (SNPs) within FASN in 921 Hu rams with IMF content records were investigated using DNA-pooling sequencing and improved multiple ligase detection reaction (iMLDR) methods. Consequently, the highest mRNA expression level of FASN was observed in the perinephric fat, and the lowest in the liver among the 11 tissues analyzed, while no significant difference was found in mRNA and protein expression levels in longissimus dorsi among individuals with different IMF contents. A total of 10 putative SNPs were identified within FASN, and 9 of them can be genotyped by iMLDR method. Notably, two SNPs were significantly associated with IMF content, including NC_040262.1: g.5157 A > G in intron 5 (p = 0.046) and NC_040262.1: g.9413 T > C in intron 16 (p = 0.041), which supply molecular markers for improving meat quality in sheep breeding.

Expression quantitative trait loci in sheep liver and muscle contribute to variations in meat traits.

BACKGROUND: Variants that regulate transcription, such as expression quantitative trait loci (eQTL), have shown enrichment in genome-wide association studies (GWAS) for mammalian complex traits. However, no study has reported eQTL in sheep, although it is an important agricultural species for which many GWAS of complex meat traits have been conducted. Using RNA sequence data produced from liver and muscle from 149 sheep and imputed whole-genome single nucleotide polymorphisms (SNPs), our aim was to dissect the genetic architecture of the transcriptome by associating sheep genotypes with three major molecular phenotypes including gene expression (geQTL), exon expression (eeQTL) and RNA splicing (sQTL). We also examined these three types of eQTL for their enrichment in GWAS of multi-meat traits and fatty acid profiles. RESULTS: Whereas a relatively small number of molecular phenotypes were significantly heritable (h2 > 0, P < 0.05), their mean heritability ranged from 0.67 to 0.73 for liver and from 0.71 to 0.77 for muscle. Association analysis between molecular phenotypes and SNPs within ± 1 Mb identified many significant cis-eQTL (false discovery rate, FDR < 0.01). The median distance between the eQTL and transcription start sites (TSS) ranged from 68 to 153 kb across the three eQTL types. The number of common variants between geQTL, eeQTL and sQTL within each tissue, and the number of common variants between liver and muscle within each eQTL type were all significantly (P < 0.05) larger than expected by chance. The identified eQTL were significantly (P < 0.05) enriched in GWAS hits associated with 56 carcass traits and fatty acid profiles. For example, several geQTL in muscle mapped to the FAM184B gene, hundreds of sQTL in liver and muscle mapped to the CAST gene, and hundreds of sQTL in liver mapped to the C6 gene. These three genes are associated with body composition or fatty acid profiles. CONCLUSIONS: We detected a large number of significant eQTL and found that the overlap of variants between eQTL types and tissues was prevalent. Many eQTL were also QTL for meat traits. Our study fills a gap in the knowledge on the regulatory variants and their role in complex traits for the sheep model.

Gallic acid caused cultured mice TM4 Sertoli cells apoptosis and necrosis.

OBJECTIVE: The study was designed to determine the cytotoxic effect of gallic acid (GA), obtained by the hydrolysis of tannins, on mice TM4 Sertoli cells apoptosis. METHODS: In the present study, non-tumorigenic mice TM4 Sertoli cells were treated with different concentrations of GA for 24 h. After treatment, cell viability was evaluated using WST-1, mitochondrial dysfunction, cells apoptosis and necrosis was detected using JC-1, Hoechst 33342 and propidium iodide staining. The expression levels of Cyclin B1, proliferating cell nuclear antigen (PCNA), Bcl-2-associated X protein (BAX), and Caspase-3 were also detected by quantitative real-time polymerase chain reaction and Western-blotting. RESULTS: The results showed that 20 to 400 µM GA inhibited viability of TM4 Sertoli cells in a dose-dependent manner. Treatment with 400 µM GA significantly inhibited PCNA and Cyclin B1 expression, however up-regulated BAX and Caspase-3 expression, caused mitochondrial membrane depolarization, activated Caspase-3, and induced DNA damage, thus, markedly increased the numbers of dead cells. CONCLUSION: Our findings showed that GA could disrupt mitochondrial function and caused TM4 cells to undergo apoptosis and necrosis.

Comprehensive investigation of nucleotide diverdity in yaks.

To understand the maternal genetic diversity of Tianzhu white yak better, we analyzed mtDNA D-loop sequences of 209 Tianzhu white yaks, which are from the central region of Tianzhu white yak habitat. Accordingly, a total of 45 haplotypes were identified in Tianzhu white yaks in this study, and 18 of them were unique. The nucleotide diversity and haplotype diversity of population studied were 0.024 ± 0.003 and 0.946 ± 0.007 respectively, revealing that Tianzhu white yak possess a relatively high genetic diversity. The phylogenetic analysis, combining D-loop sequences in this study with 533 previous published D-loop sequences of 13 yak breeds, indicated that Tianzhu white yaks fell mainly into haplogroup A and that a small portion belonged to haplogroups B, C, D and E. Moreover, six haplotypes of 20 individuals identified in Tianzhu white yak were in the taurine haplogroup, indicating hybridization between Bos taurus and Tianzhu white yaks. In summary, this study supplies a comprehensive maternal genetic pattern for Tianzhu white yak and provides a basic reference for future breeding programs to conserve the purebred Tianzhu white yak.

A limited number of Y chromosome lineages is present in North American Holsteins.

Holsteins are the most numerous dairy cattle breed in North America and the breed has undergone intensive selection for improving milk production and conformation. Theoretically, this intensive selection could lead to a reduction of the effective population size and reduced genetic diversity. The objective of this study was to investigate the effective population size of the Holstein Y chromosome and the effects of limited Y chromosome lineages on male reproduction and the future of the breed. Paternal pedigree information of 62,897 Holstein bulls born between 1950 and 2013 in North America and 220,872 bulls evaluated by multiple-trait across-country genetic evaluations of Interbull (Uppsala, Sweden) were collected and analyzed. The results indicated that the number of Y chromosome lineages in Holsteins has undergone a dramatic decrease during the past 50 years because of artificial selection and the application of artificial insemination (AI) technology. All current Holstein AI bulls in North America are the descendants of only 2 ancestors (Hulleman and Neptune H) born in 1880. These 2 ancestral Y-lineages are continued through 3 dominant pedigrees from the 1960s; namely, Pawnee Farm Arlinda Chief, Round Oak Rag Apple Elevation, and Penstate Ivanhoe Star, with a contribution of 48.78, 51.06, and 0.16% to the Holstein bull population in the 2010s, respectively. The Y-lineage of Penstate Ivanhoe Star is almost eliminated from the breed. The genetic variations in the 2 ancestral Y-lineages were evaluated among 257 bulls by determining the copy number variations (CNV) of 3 Y-linked gene families: PRAMEY, HSFY, and ZNF280BY, which are spread along the majority (95%) of the bovine Y chromosome male-specific region (MSY). No significant difference was found between the 2 ancestral Y-lineages, although large CNV were observed within each lineage. This study suggests minimal genetic diversity on the Y chromosome in Holsteins and provides a starting point for investigating the effect of the extremely limited number of Y-lineages on male reproduction and other traits important for the future of the Holstein breed.

Y-Single Nucleotide Polymorphisms Diversity in Chinese Indigenous Horse.

In contrast to high genetic diversity of mitochondrial DNA (mtDNA), equine Y chromosome shows extremely low variability, implying limited patrilines in the domesticated horse. In this study, we applied direct sequencing and restriction fragment length polymorphism (RFLP) methods to investigate the polymorphisms of 33 Y chromosome specific loci in 304 Chinese indigenous horses from 13 breeds. Consequently, two Y-single nucleotide polymorphisms (SNPs) (Y-45701/997 and Y-50869) and one Y-indel (Y-45288) were identified. Of those, the Y-50869 (T>A) revealed the highest variation frequency (24.67%), whereas it was only 3.29% and 1.97% in Y-45288 (T/-) and Y-45701/997 (G>T) locus, respectively. These three mutations accounted for 27.96% of the total samples and identified five Y-SNP haplotypes, demonstrating genetic diversity of Y chromosome in Chinese horses. In addition, all the five Y-SNP haplotypes were shared by different breeds. Among 13 horse breeds analyzed, Balikun horse displayed the highest nucleotide diversity (π = 5.6×10(-4)) and haplotype diversity (h = 0.527), while Ningqiang horse showed the lowest nucleotide diversity (π = 0.00000) and haplotype diversity (h = 0.000). The results also revealed that Chinese horses had a different polymorphic pattern of Y chromosome from European and American horses. In conclusion, Chinese horses revealed genetic diversity of Y chromosome, however more efforts should be made to better understand the domestication and paternal origin of Chinese indigenous horses.

Copy number variations of the extensively amplified Y-linked genes, HSFY and ZNF280BY, in cattle and their association with male reproductive traits in Holstein bulls.

BACKGROUND: Recent transcriptomic analysis of the bovine Y chromosome revealed at least six multi-copy protein coding gene families, including TSPY, HSFY and ZNF280BY, on the male-specific region (MSY). Previous studies indicated that the copy number variations (CNVs) of the human and bovine TSPY were associated with male fertility in men and cattle. However, the relationship between CNVs of the bovine Y-linked HSFY and ZNF280BY gene families and bull fertility has not been investigated. RESULTS: We investigated the copy number (CN) of the bovine HSFY and ZNF280BY in a total of 460 bulls from 15 breeds using a quantitative PCR approach. We observed CNVs for both gene families within and between cattle breeds. The median copy number (MCN) of HSFY among all bulls was 197, ranging from 21 to 308. The MCN of ZNF280BY was 236, varying from 28 to 380. Furthermore, bulls in the Bos taurus (BTA) lineage had a significantly higher MCN (202) of HSFY than bulls in the Bos indicus (BIN) lineage (178), while taurine bulls had a significantly lower MCN (231) of ZNF280BY than indicine bulls (284). In addition, the CN of ZNF280BY was positively correlated to that of HSFY on the BTAY. Association analysis revealed that the CNVs of both HSFY and ZNF280BY were correlated negatively with testis size, while positively with sire conception rate. CONCLUSION: The bovine HSFY and ZNF280BY gene families have extensively expanded on the Y chromosome during evolution. The CN of both gene families varies significantly among individuals and cattle breeds. These variations were associated with testis size and bull fertility in Holstein, suggesting that the CNVs of HSFY and ZNF280BY may serve as valuable makers for male fertility selection in cattle.

Molecular characterization of SPATA6 and association of its SNPs with testicular size in sheep.

The testis is an important organ for maintaining fertility in males, and testis size is positively correlated with ejaculate volume, sperm motility, thus fertility. Spermatogenesis-associated 6 (SPATA6) is an evolutionarily conserved testis-specific gene reported in many species. However, the effect of SPATA6 expression levels on testicular development and the effect of single nucleotide polymorphisms (SNPs) on testis and epididymis phenotype in sheep have not been studied. The purpose of the research was to investigate the expression profile of SPATA6 and its effect on testicular development and to confirm the effect of SNPs on the testis and epididymis phenotype. In this study, we detected a 1245bp coding sequence (CDS) of SPATA6 and encoded 414 amino acids. The expression levels of SPATA6 were significantly higher in the testis than in other tissues and gradually increased with testis development. Moreover, the expression level in the large testis was significantly higher than that in the small testis at six months. A total of 11 SNPs were detected in the coding region of SPATA6 by cDNA-pooling sequencing and improved multiplex ligation detection reaction (iMLDR) methods. Correlation analysis showed that SNP2 (c. 3631C > G) significantly affected left epididymis weight (LEW) and right epididymis weight (REW), and SNP10 (c. 937 A > G) significantly affected REW. And the combined genotype of SNP1 (c. 4245 G > A) and SNP2 significantly affected REW. The current study concluded that SPATA6 plays an important role in testicular development and the SNPs significantly associated with the epididymis phenotype can provide molecular markers for the early selection of high-fertility Hu sheep.

SNPS within the SLC27A6 gene are highly associated with Hu sheep fatty acid content.

Fatty acids (FA) are an important factor affecting meat quality and human health, and the important role of the solute carrier family 27 member 6 (SLC27A6) in FA metabolism has been demonstrated in several species. However, the expression profile of the SLC27A6 in different tissues and the effect of its polymorphism on FA in sheep are currently unknown. This study aimed to explore the differences in FAs in the longissimus dorsi (LD) of 1,085 Hu sheep, the expression profile of SLC27A6, and confirm the effect of single nucleotide polymorphisms (SNPs) on FA phenotypes. We found that many FA phenotypes differ significantly across different seasons, and winter promoted the deposition of polyunsaturated fatty acids (PUFA). The mRNA expression level of SLC27A6 in the lung was significantly higher than that in the heart, testis, and LD. A total of 16 SNPs were detected in the SLC27A6, and 14 SNPs were successfully genotyped by improved multiplex ligase detection reaction (iMLDR) technology. Correlation analysis showed that 7 SNPs significantly affected at least one FA phenotype. Among them, SNP14 contributes to the selection of lamb with low saturated fatty acid content and high PUFA content. Combined genotypes also significantly affected a variety of beneficial FAs such as C18:3n3, C20:4n6, C22:6n3, and monounsaturated fatty acids. This study suggests that SLC27A6 plays an important role in FA metabolism and SNPs that are significantly associated with FA phenotype could be used as potential molecular markers for later targeted regulation of FA profiles in sheep.

Transcriptomic analysis of the HPG axis-related tissues reveals potential candidate genes and regulatory pathways associated with testicular size in Hu sheep.

Testicular size is an excellent proxy for selecting high-fertility rams. The hypothalamus-pituitary-gonadal (HPG) axis plays an important role in regulating reproductive capacity in vertebrates, while key genes and regulatory pathways within the HPG axis associated with testicular size remain largely unknown in sheep. This study comprehensively compared the transcriptomic profiles in the hypothalamus, pituitary and testis of rams after sexual maturity between the large-testis group (LTG, testicular weight = 454.29 ± 54.24 g) and the small-testis group (STG, testicular weight = 77.29 ± 10.76 g). In total, 914, 795 and 10518 differentially expressed genes (DEGs) were identified in the hypothalamus, pituitary and testis between LTG and STG, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that these DEGs were mainly involved in the biological processes of reproduction, biological regulation, and development process. Notably, the neuroactive ligand-receptor interaction and cAMP signaling pathways, commonly enriched by the DEGs in the hypothalamus and pituitary between two groups, were considered as two key signal pathways regulating testicular development through the HPGs axis. Weighted gene co-expression network analysis (WGCNA) identified two modules that were significantly associated with testicular size, and 97 key genes were selected with high module membership (MM) and gene significance (GS) in these two modules. Finally, a protein-protein interaction (PPI) network was constructed, and ten genes with the highest degree were represented as hub genes, including FOS, NPY, SST, F2, AGT, NTS, OXT, EDN1, VIP and TAC1. Taken together, these results provide new insights into the molecular mechanism underlying the HPG axis regulating testicular size of Hu sheep.

Preliminary genetic parameter estimates of meat quality traits in Hu sheep.

Because substantial numbers of Chinese consumers are prepared to pay for tender and quality lamb, meat quality traits are becoming more relevant for breeding programs for Chinese sheep breeds. The current study estimated heritabilities and genetic correlations for 13 meat quality traits recorded on lamb loins from Hu sheep. Heritability estimates ranged from 0.04 ± 0.06 for meat redness at 45 min to 0.57 ± 0.10 for drip loss, with most of the meat quality traits having moderate heritabilities. Positive genetic correlations were observed among meat color traits. Intramuscular fat (IMF) was genetically correlated with most meat quality traits, indicating that increasing IMF can favor meat pH, color, and tenderness, but would lead to increased cooking loss. Direct selection to increase IMF of loins is recommended to be included in breeding programs for Hu sheep, as it was more efficient than indirect selection on the other meat quality traits. The genetic parameters presented in this preliminary study provide valuable genetic information needed to design a breeding program aimed at improving the quality of lamb meat from Hu sheep.

Localization and expression of phospholipase A2 and polyunsaturated fatty acid profile in the testis tissues of Hu sheep.

The fatty acid content and the localization and expression of phospholipase A2 (PLA2) in the testis of Hu sheep were investigated. A total of 18 six-month-old Hu sheep were divided into small group (S, with left testis weight < 50 g), medium group (M, with left testis weight among 90-110 g), and large group (L, with left testis weight >160 g), which had six individuals each. The expression of PLA2 in testicular tissues of different sizes was analyzed by immunohistochemistry, RT-qPCR, and Western blot. The fatty acid profile was detected by gas chromatography. Immunohistochemical labeling determined that PLA2 protein was expressed in the Leydig and Sertoli cells of testis, and the immunohistochemical average optional density in the S group was significantly greater than the L group (P < 0.05). RT-qPCR and Western blot analysis showed that PLA2 in the S group was greater than that in the L group (P < 0.05). Docosahexaenoic acid, ω-3 polyunsaturated fatty acid (PUFA), and total PUFA content in the testis of the L group were significantly less than those of the S and M groups (P < 0.01). This study showed that PLA2 content in the S group was greater than that in the L group.

Identification of Y-SNPs within ovine MSY region and their association with testicular size.

The screening of genomic variations within the male-specific region of the mammalian Y chromosome (MSY) is one of the most effective ways to investigate paternal evolutionary history and identify molecular markers related to male fertility. The current study was to identify single nucleotide polymorphisms (SNPs) within single-copy genes of the ovine MSY, and confirm whether they are associated with testicular size. A total of 21 Y-specific gene fragments were successfully amplified to screen Y-SNPs in 956 rams across nine sheep breeds. Three Y-SNPs, including SRY16: g.88 A > G in South African Mutton Merino sheep, ZFY16: g.146 C > T in Suffolk and South African Mutton Merino sheep, and EIF2S3Y2: g.77 C > G in Hu and Tan sheep, were identified using DNA-pooled sequencing and PCR restriction fragment length polymorphism (PCR-RFLP) methods. The investigation of the global distribution for three Y-SNPs showed that the C allele of ZFY16: g.146 C > T co-segregated with haplogroup y-HC, and the C/G allele of EIF2S3Y2: g.77 C > G co-segregated with haplogroup y-HA/y-HB1 in Hu sheep according to data mining from a previous study. In addition, association analysis revealed that ZFY16: g.146 C > T had a significant effect on yearling scrotal circumference in Suffolk sheep, and EIF2S3Y2: g.77 C > G was significantly associated with testicular and epididymis weight in Hu sheep (P ≤ 0.05). The current study concluded that Y-SNPs were associated with testicular size in specific sheep, which provides valuable candidate makers for selecting elite rams at an early age.