Genome-wide detection of genetic structure and runs of homozygosity analysis in Anhui indigenous and Western commercial pig breeds using PorcineSNP80k data.

BACKGROUND: Runs of homozygosity (ROH) are continuous homozygous regions typically located in the DNA sequence of diploid organisms. Identifications of ROH that lead to reduced performance can provide valuable insight into the genetic architecture of complex traits. Here, we systematically investigated the population genetic structure of five Anhui indigenous pig breeds (AHIPs), and compared them to those of five Western commercial pig breeds (WECPs). Furthermore, we examined the occurrence and distribution of ROHs in the five AHIPs and estimated the inbreeding coefficients based on the ROHs (FROH) and homozygosity (FHOM). Finally, we identified genomic regions with high frequencies of ROHs and annotated candidate genes contained therein. RESULTS: The WECPs and AHIPs were clearly differentiated into two separate clades consistent with their geographical origins, as revealed by the population structure and principal component analysis. We identified 13,530 ROHs across all individuals, of which 4,555 and 8,975 ROHs were unique to AHIPs and WECPs, respectively. Most ROHs identified in our study were short (< 10 Mb) or medium (10-20 Mb) in length. WECPs had significantly higher numbers of short ROHs, and AHIPs generally had longer ROHs. FROH values were significantly lower in AHIPs than in WECPs, indicating that breed improvement and conservation programmes were successful in AHIPs. On average, FROH and FHOM values were highly correlated (0.952-0.991) in AHIPs and WECPs. A total of 27 regions had a high frequency of ROHs and contained 17 key candidate genes associated with economically important traits in pigs. Among these, nine candidate genes (CCNT2, EGR2, MYL3, CDH13, PROX1, FLVCR1, SETD2, FGF18, and FGF20) found in WECPs were related to muscular and skeletal development, whereas eight candidate genes (CSN1S1, SULT1E1, TJP1, ZNF366, LIPC, MCEE, STAP1, and DUSP) found in AHIPs were associated with health, reproduction, and fatness traits. CONCLUSION: Our findings provide a useful reference for the selection and assortative mating of pig breeds, laying the groundwork for future research on the population genetic structures of AHIPs, ultimately helping protect these local varieties.

Enzymatic Synthesis of the Unnatural Nucleotide 2'-Deoxyisoguanosine 5'-Monophosphate.

Naturally occurring DNA contains four canonical bases, forming two Watson-Crick base pairs (adenine-thymine, guanine-cytosine). Efforts over the past decades have led to the development of several unnatural base pairs, enabling the synthesis of unnatural DNA with an expanded genetic alphabet. The engineering of organisms capable of de novo biosynthesis of unnatural DNA would have significant technological and philosophical implications, but remains a challenge. Here we report the enzymatic conversion of 2'-deoxyxanthosine 5'-monophosphate (dXMP) into deoxyisoguanosine monophosphate (dBMP), a precursor of the unnatural isoguanine-isocytosine base pair. The reaction is catalyzed by the bacteriophage enzyme PurZ and bacterial PurB, and is a key addition to the toolbox for de novo biosynthesis of unnatural DNA.

Copy number variation association studies for sheep tail-relevant traits in Hulunbuir sheep.

Tail type is an important trait that influences meat quality and consumer purchasing attitudes. As a novel genetic marker, the study of genomic copy number variations (CNVs) provides a new research method to study the genetic mechanisms underlying trait formation. In the present paper, we conducted CNV-based association studies for sheep tail type and growth traits in Hulunbuir sheep. A total of 167 CNVs and 288 individuals were analyzed at the marker and individual level, respectively. Association studies of 10 traits were carried out using two mixed linear models. We found that three, one and one CNV loci were significantly associated with tail type, tail length and tail fat weight, respectively, which together contained a total of 52 candidate genes. Out of these genes, five have been documented to be associated with fat metabolism in sheep. Our findings provide a reference for further studies on fat deposition in sheep as well as a theoretical basis for breeding Hulunbuir sheep with the required tail type.

Identification of imprinted genes in the skeletal muscle of newborn piglets by high-throughput sequencing.

Imprinted genes - exhibiting parent-specific transcription - play essential roles in the process of mammalian development and growth. Skeletal muscle growth is crucial for meat production. To further understand the role of imprinted genes during the porcine skeletal muscle growth, DNA-seq and RNA-seq were used to explore the characteristics of imprinted genes from porcine reciprocal crosses. A total of 584 545 single-nucleotide variations were discovered in the DNA-seq data of F0 parents, heterozygous in two pig breeds (Yorkshire and Min pigs) but homozygous in each breed. These single-nucleotide variations were used to determine the allelic-specific expression in F1 individuals. Finally, eight paternal expression sites and three maternal expression sites were detected, whereas two paternally expressed imprinted genes (NDN and IGF2) and one maternally expressed imprinted gene (H1-3) were validated by Sanger sequencing. DNA methylation regulates the expression of imprinted genes, and all of the identified imprinted genes in this study were predicted to possess CpG islands. PBX1 and YY1 binding motifs were discovered in the promoter regions of all three imprinted genes, which were candidate elements regulating the transcription of imprinted genes. For these identified imprinted genes, IGF2 and NDN promoted muscle growth whereas H1-3 inhibited cell proliferation, corroborating the 'parental conflict' theory that paternally expressed imprinted genes assisted descendants' growth whereas maternally expressed imprinted genes inhibited it. This study discovered porcine imprinted genes in skeletal muscle and was the first to reveal that H1-3 was expressed by the maternal allele to our knowledge. Our findings provided valuable resources for the potential utilization of imprinted genes in pig breeding.

Genome-wide association study revealed ABCD4 on SSC7 and GREB1L and MIB1 on SSC6 as crucial candidate genes for rib number in Beijing Black pigs.

As one of the few animals with variation in the number of rib pairs (RIB), the pig is a good model to study the mechanism of RIB regulation. Quantitative trait loci (QTL) for porcine RIB are present on Sus scrofa chromosome 7 (SSC7). Although several candidate genes exist in this QTL region on SSC7, the causal gene has yet to be verified. Beijing Black pig with 14-17 RIB is a good population for candidate gene mining and 1104 individuals were genotyped using the Illumina Porcine 50K BeadChip. A total of 14 SNPs from 95.49 to 97.78 Mb on SSC7 showed genome-wide significant association with RIB. On SSC7, a locuszoom plot using pairwise linkage disequilibrium displayed the narrowest linkage region encompassing only two genes, ABCD4 and VRTN. In mice, a transcriptome expression profile was obtained using three embryos at E9.5 (the critical period for rib formation). ABCD4 was highly expressed, but no expression of VRTN was detected. On SSC6, there were four genome-wide significant SNPs from 106.42 to 106.92 Mb associated with RIB. GREB1L and MIB1, in this region, were regarded as novel candidate genes. These results revealed a crucial candidate causal gene on SSC7 and novel genes on SSC6 for rib number and provided interesting new insights into its genetic basis.

GWAS of Reproductive Traits in Large White Pigs on Chip and Imputed Whole-Genome Sequencing Data.

Total number born (TNB), number of stillborn (NSB), and gestation length (GL) are economically important traits in pig production, and disentangling the molecular mechanisms associated with traits can provide valuable insights into their genetic structure. Genotype imputation can be used as a practical tool to improve the marker density of single-nucleotide polymorphism (SNP) chips based on sequence data, thereby dramatically improving the power of genome-wide association studies (GWAS). In this study, we applied Beagle software to impute the 50 K chip data to the whole-genome sequencing (WGS) data with average imputation accuracy (R2) of 0.876. The target pigs, 2655 Large White pigs introduced from Canadian and French lines, were genotyped by a GeneSeek Porcine 50K chip. The 30 Large White reference pigs were the key ancestral individuals sequenced by whole-genome resequencing. To avoid population stratification, we identified genetic variants associated with reproductive traits by performing within-population GWAS and cross-population meta-analyses with data before and after imputation. Finally, several genes were detected and regarded as potential candidate genes for each of the traits: for the TNB trait: NOTCH2, KLF3, PLXDC2, NDUFV1, TLR10, CDC14A, EPC2, ORC4, ACVR2A, and GSC; for the NSB trait: NUB1, TGFBR3, ZDHHC14, FGF14, BAIAP2L1, EVI5, TAF1B, and BCAR3; for the GL trait: PPP2R2B, AMBP, MALRD1, HOXA11, and BICC1. In conclusion, expanding the size of the reference population and finding an optimal imputation strategy to ensure that more loci are obtained for GWAS under high imputation accuracy will contribute to the identification of causal mutations in pig breeding.

Liver transcriptome profiling and functional analysis of intrauterine growth restriction (IUGR) piglets reveals a genetic correction and sexual-dimorphic gene expression during postnatal development.

BACKGROUND: Intrauterine growth restriction (IUGR) remains a major problem associated with swine production. Thus, understanding the physiological changes of postnatal IUGR piglets would aid in improving growth performance. Moreover, liver metabolism plays an important role in the growth and survival of neonatal piglets. RESULTS: By profiling the transcriptome of liver samples on postnatal Days 1, 7, and 28, our study focused on characterizing the growth, function, and metabolism in the liver of IUGR neonatal piglets. Our study demonstrates that the livers of IUGR piglets were associated with a series of complications, including inflammatory stress and immune dysregulation; cytoskeleton and membrane structure disorganization; dysregulated transcription events; and abnormal glucocorticoid metabolism. In addition, the abnormal liver function index in the serum [alanine aminotransferase (ALT), aspartate aminotransferase (AST), and total protein (TP)], coupled with hepatic pathological and ultrastructural morphological changes are indicative of liver damage and dysfunction in IUGR piglets. Moreover, these results reveal the sex-biased developmental dynamics between male and female IUGR piglets, and that male IUGR piglets may be more sensitive to disrupted metabolic homeostasis. CONCLUSIONS: These observations provide a detailed reference for understanding the mechanisms and characterizations of IUGR liver functions, and suggest that the potential strategies for improving the survival and growth performance of IUGR offspring should consider the balance between postnatal catch-up growth and adverse metabolic consequences. In particular, sex-specific intervention strategies should be considered for both female and male IUGR piglets.

Detection of copy number variation and selection signatures on the X chromosome in Chinese indigenous sheep with different types of tail.

OBJECTIVE: Chinese indigenous sheep breeds can be classified into the following three categories by their tail morphology: fat-tailed, fat-rumped and thin-tailed sheep. The typical sheep breeds corresponding to fat-tailed, fat-rumped, and thin-tailed sheep are large-tailed Han, Altay, and Tibetan sheep, respectively. Detection of copy number variation (CNV) and selection signatures provides information on the genetic mechanisms underlying the phenotypic differences of the different sheep types. METHODS: In this study, PennCNV software and F-statistics (FST) were implemented to detect CNV and selection signatures, respectively, on the X chromosome in three Chinese indigenous sheep breeds using ovine high-density 600K single nucleotide polymorphism arrays. RESULTS: In large-tailed Han, Altay, and Tibetan sheep, respectively, a total of six, four and 22 CNV regions (CNVRs) with lengths of 1.23, 0.93, and 7.02 Mb were identified on the X chromosome. In addition, 49, 34, and 55 candidate selection regions with respective lengths of 27.49, 16.47, and 25.42 Mb were identified in large-tailed Han, Altay, and Tibetan sheep, respectively. The bioinformatics analysis results indicated several genes in these regions were associated with fat, including dehydrogenase/reductase X-linked, calcium voltage-gated channel subunit alpha1 F, and patatin like phospholipase domain containing 4. In addition, three other genes were identified from this analysis: the family with sequence similarity 58 member A gene was associated with energy metabolism, the serine/arginine-rich protein specific kinase 3 gene was associated with skeletal muscle development, and the interleukin 2 receptor subunit gamma gene was associated with the immune system. CONCLUSION: The results of this study indicated CNVRs and selection regions on the X chromosome of Chinese indigenous sheep contained several genes associated with various heritable traits.

Genome-wide association studies revealed candidate genes for tail fat deposition and body size in the Hulun Buir sheep.

Fat-tailed sheep have a unique characteristic of depositing fat in their tails. In the present study, we conducted genome-wide association studies (GWAS) on traits related to tail fat deposition and body size in the Hulun Buir sheep. A total number of 300 individuals belonging to two fat-tailed lines of the Hulun Buir sheep breed genotyped with the Ovine Infinium HD SNP BeadChip were included in the current study. Two mixed models, one for continuous and one for binary phenotypic traits, were employed to analyse ten traits, that is, body length (BL), body height (BH), chest girth (CG), tail length (TL), tail width (TW), tail circumference (TC), carcass weight (CW), tail fat weight (TF), ratio of CW to TF (RCT) and tail type (TT). We identified 7, 6, 7, 2, 10 and 1 SNPs significantly associated with traits TF, CW, RCT, TW, TT and CG, respectively. Their associated genomic regions harboured 42 positional candidate genes. Out of them, 13 candidate genes including SMURF2, FBF1, DTNBP1, SETD7 and RBM11 have been associated with fat metabolism in sheep. The RBM11 gene has already been identified in a previous study on signatures of selection in this specific sheep population. Two more genes, that is, SMARCA5 and GAB1 were associated with body size in sheep. The present study has identified candidate genes that might be implicated in tail fat deposition and body size in sheep.

Advances in genome-wide association studies for important traits in sheep and goats.

Genome-wide association study (GWAS), an effective strategy to identify genetic variants associated with complex traits, has been used to study candidate genes of economical traits in animals. With the recent completion of sheep and goat genomes, single nucleotide polymorphism (SNP) chips of different densities are developed and commercialized. All these advances have enlarged the collection of molecular markers and also shed new light on the genetics of traits of interest in sheep and goats. In this review, we focus on the adoption of GWAS for important traits in sheep and goats, such as horn types, wool, dairy, growth and meat, reproduction and disease types, etc., and summarize the populations, major statistical methods and results of the GWAS analysis. Moreover, we also discuss the current state of GWAS, aiming to provide a reference for further studies on the genetic background of the important traits of sheep and goats by GWAS.

Screening of reproduction-related single-nucleotide variations from MeDIP-seq data in sheep.

Extensive variation in reproduction has arisen in Chinese Mongolian sheep during recent domestication. Hu and Small-tailed Han sheep, for example, have become non-seasonal breeders and exhibit higher fecundity than Tan and Ujumqin breeds. We therefore scanned reproduction-related single-nucleotide variations from methylated DNA-immunoprecipitation sequencing data generated from each of those four breeds to uncover potential mechanisms underlying this breed variation. We generated a high-quality map of single nucleotide variations (SNVs) in DNA methylation enriched regions, and found that the majority of variants are located within non-coding regions. We identified 359 SNVs within the Sheep Quantitative Trait Locus (QTL) database. Nineteen of these SNVs associated with the Aseasonal Reproduction QTL, and 10 out of the 19 reside close to genes with known reproduction functions. We also identified the well-known FecB mutation in high-fecundity sheep (Hu and Small-tailed Han sheep). When we applied these FecB finding to our breeding system, we improved lambing rate by 175%. In summary, this study provided strong candidate SNVs associated with sheep fecundity that can serve as targets for functional testing and to enhance selective breeding strategies. Mol. Reprod. Dev. 83: 958-967, 2016 © 2016 Wiley Periodicals, Inc.

Effect of polymorphisms in the CAMKMT gene on growth traits in Ujumqin sheep.

Our previous genome-wide association study in sheep revealed that OAR3-84073899.1 (SNP31) in intron 8 of the CAMKMT gene was significantly associated with post-weaning gain at the genomic level. Herein, we performed a replication study to investigate single nucleotide polymorphisms (SNPs) within the CAMKMT gene exons, and 1000 bp of the 5'- and 3'-intranslated regions (UTRs) and their associations with growth traits in Ujumqin sheep. Five SNPs were identified through DNA pool sequencing technology: SNP26 in the 5'-UTR, SNP06 in exon 5, SNP07 in exon 8 and SNP27 and SNP28 in the 3'-UTR. Six SNPs, including SNP31 in intron 8, were genotyped in the validation group of 343 Ujumqin sheep, and each SNP was classified into three genotypes. The chi-square test suggested that all the variations were in Hardy-Weinberg equilibrium (P > 0.05) except for SNP28 and SNP31. Linkage disequilibrium analysis showed that SNP07 and SNP31 were strongly linked. An association analysis suggested that SNP06 was significantly associated with chest girth at 6 months of age (P < 0.05). SNP07 exhibited significant correlation with body weight and chest girth at 4 months of age and with body weight, chest girth and chest width at 6 months of age (P < 0.05). SNP27 was highly associated with body weight and chest girth at 4 months of age (P < 0.05), and SNP28 was extremely significantly associated with body weight and chest girth at 4 months of age and with chest girth at 6 months of age (P < 0.01). SNP31 was significantly associated with body weight and shin circumference at 4 months of age and with post-weaning gain (P < 0.05). Association analysis of the combined effect of SNP07 and SNP31 showed significant correlation with body weight and chest girth at four of months of age (P < 0.05) and with body weight and chest girth at 6 months of age (P < 0.05). These results indicate that the SNPs could be used as meritorious and available genetic markers in growth traits breeding and that the CAMKMT gene may be one of the key candidate genes that affect Ujumqin economic traits.

Complete Mitochondrial Genome Sequences of Chinese Indigenous Sheep with Different Tail Types and an Analysis of Phylogenetic Evolution in Domestic Sheep.

China has a long history of sheep (Ovis aries [O. aries]) breeding and an abundance of sheep genetic resources. Knowledge of the complete O. aries mitogenome should facilitate the study of the evolutionary history of the species. Therefore, the complete mitogenome of O. aries was sequenced and annotated. In order to characterize the mitogenomes of 3 Chinese sheep breeds (Altay sheep [AL], Shandong large-tailed sheep [SD], and small-tailed Hulun Buir sheep [sHL]), 19 sets of primers were employed to amplify contiguous, overlapping segments of the complete mitochondrial DNA (mtDNA) sequence of each breed. The sizes of the complete mitochondrial genomes of the sHL, AL, and SD breeds were 16,617 bp, 16,613 bp, and 16,613 bp, respectively. The mitochondrial genomes were deposited in the GenBank database with accession numbers KP702285 (AL sheep), KP981378 (SD sheep), and KP981380 (sHL sheep) respectively. The organization of the 3 analyzed sheep mitochondrial genomes was similar, with each consisting of 22 tRNA genes, 2 rRNA genes (12S rRNA and 16S rRNA), 13 protein-coding genes, and 1 control region (D-loop). The NADH dehydrogenase subunit 6 (ND6) and 8 tRNA genes were encoded on the light strand, whereas the rest of the mitochondrial genes were encoded on the heavy strand. The nucleotide skewness of the coding strands of the 3 analyzed mitogenomes was biased toward A and T. We constructed a phylogenetic tree using the complete mitogenomes of each type of sheep to allow us to understand the genetic relationships between Chinese breeds of O. aries and those developed and utilized in other countries. Our findings provide important information regarding the O. aries mitogenome and the evolutionary history of O. aries inside and outside China. In addition, our results provide a foundation for further exploration of the taxonomic status of O. aries.

Genome-wide analysis reveals population structure and selection in Chinese indigenous sheep breeds.

BACKGROUND: Traditionally, Chinese indigenous sheep were classified geographically and morphologically into three groups: Mongolian, Kazakh and Tibetan. Herein, we aimed to evaluate the population structure and genome selection among 140 individuals from ten representative Chinese indigenous sheep breeds: Ujimqin, Hu, Tong, Large-Tailed Han and Lop breed (Mongolian group); Duolang and Kazakh (Kazakh group); and Diqing, Plateau-type Tibetan, and Valley-type Tibetan breed (Tibetan group). RESULTS: We analyzed the population using principal component analysis (PCA), STRUCTURE and a Neighbor-Joining (NJ)-tree. In PCA plot, the Tibetan and Mongolian groups were clustered as expected; however, Duolang and Kazakh (Kazakh group) were segregated. STRUCTURE analyses suggested two subpopulations: one from North China (Kazakh and Mongolian groups) and the other from the Southwest (Tibetan group). In the NJ-tree, the Tibetan group formed an independent branch and the Kazakh and Mongolian groups were mixed. We then used the d i statistic approach to reveal selection in Chinese indigenous sheep breeds. Among the 599 genome sequence windows analyzed, sixteen (2.7%) exhibited signatures of selection in four or more breeds. We detected three strong selection windows involving three functional genes: RXFP2, PPP1CC and PDGFD. PDGFD, one of the four subfamilies of PDGF, which promotes proliferation and inhibits differentiation of preadipocytes, was significantly selected in fat type breeds by the Rsb (across pairs of populations) approach. Two consecutive selection regions in Duolang sheep were obviously different to other breeds. One region was in OAR2 including three genes (NPR2, SPAG8 and HINT2) the influence growth traits. The other region was in OAR 6 including four genes (PKD2, SPP1, MEPE, and IBSP) associated with a milk production quantitative trait locus. We also identified known candidate genes such as BMPR1B, MSRB3, and three genes (KIT, MC1R, and FRY) that influence lambing percentage, ear size and coat phenotypes, respectively. CONCLUSIONS: Based on the results presented here, we propose that Chinese native sheep can be divided into two genetic groups: the thin type (Tibetan group), and the fat type (Mongolian and Kazakh group). We also identified important genes that drive valuable phenotypes in Chinese indigenous sheep, especially PDGFD, which may influence fat deposition in fat type sheep.

Effect of thyroglobulin gene polymorphisms on growth, carcass composition and meat quality traits in Chinese beef cattle.

The thyroglobulin (TG) gene has been studied as an important gene related to fat deposition, since not only does TG gene locate in a quantitative trait locus with an effect on fat deposition, but also it encodes the precursor of thyroid hormones which have crucial biological functions in energy metabolism. In the present study, we identified four novel SNPs at the 5' flanking region of the bovine TG gene. Association analysis indicated that the G275A, G277C, G280A and C281G SNPs were significantly associated with average daily gain (ADG, P < 0.01 for G275A and G277C, P < 0.05 for G280A and C281G). Five haplotypes for the 4 SNPs were constructed and their effects on growth, carcass composition and meat quality traits were evaluated. The results showed no significant effect of haplotype on ADG. Meanwhile, no significant association was found between 4 SNPs and other growth, carcass composition and meat quality traits including intramuscular fat. Bioinformatics analysis showed that 4 SNPs may results in potential transcription factor binding site changes. Results of this study suggest that TG gene-specific SNPs may be a useful marker for growth traits in marker assisted selection programs in beef cattle.

Detection of selection signatures in dairy and beef cattle using high-density genomic information.

BACKGROUND: Artificial selection for economically important traits in cattle is expected to have left distinctive selection signatures on the genome. Access to high-density genotypes facilitates the accurate identification of genomic regions that have undergone positive selection. These findings help to better elucidate the mechanisms of selection and to identify candidate genes of interest to breeding programs. RESULTS: Information on 705 243 autosomal single nucleotide polymorphisms (SNPs) in 3122 dairy and beef male animals from seven cattle breeds (Angus, Belgian Blue, Charolais, Hereford, Holstein-Friesian, Limousin and Simmental) were used to detect selection signatures by applying two complementary methods, integrated haplotype score (iHS) and global fixation index (FST). To control for false positive results, we used false discovery rate (FDR) adjustment to calculate adjusted iHS within each breed and the genome-wide significance level was about 0.003. Using the iHS method, 83, 92, 91, 101, 85, 101 and 86 significant genomic regions were detected for Angus, Belgian Blue, Charolais, Hereford, Holstein-Friesian, Limousin and Simmental cattle, respectively. None of these regions was common to all seven breeds. Using the FST approach, 704 individual SNPs were detected across breeds. Annotation of the regions of the genome that showed selection signatures revealed several interesting candidate genes i.e. DGAT1, ABCG2, MSTN, CAPN3, FABP3, CHCHD7, PLAG1, JAZF1, PRKG2, ACTC1, TBC1D1, GHR, BMP2, TSG1, LYN, KIT and MC1R that play a role in milk production, reproduction, body size, muscle formation or coat color. Fifty-seven common candidate genes were found by both the iHS and global FST methods across the seven breeds. Moreover, many novel genomic regions and genes were detected within the regions that showed selection signatures; for some candidate genes, signatures of positive selection exist in the human genome. Multilevel bioinformatic analyses of the detected candidate genes suggested that the PPAR pathway may have been subjected to positive selection. CONCLUSIONS: This study provides a high-resolution bovine genomic map of positive selection signatures that are either specific to one breed or common to a subset of the seven breeds analyzed. Our results will contribute to the detection of functional candidate genes that have undergone positive selection in future studies.

Detection of Selection Signatures on the X Chromosome in Three Sheep Breeds.

Artificial selection has played a critical role in animal breeding. Detection of artificial selection footprints in genomic regions can provide insights for understanding the function of specific phenotypic traits and better guide animal breeding. To more fully understand the relationship between genomic composition and phenotypic diversity arising from breed development, a genome-wide scan was conducted using an OvineSNP50 BeadChip and integrated haplotype score and fixation index analyses to detect selection signatures on the X chromosome in three sheep breeds. We identified 49, 34, and 55 candidate selection regions with lengths of 27.49, 16.47, and 25.42 Mb in German Mutton, Dorper, and Sunit sheep, respectively. Bioinformatics analysis showed that some of the genes in these regions with selection signatures, such as BMP15, were relevant to reproduction. We also identified some selection regions harboring genes that had human orthologs, including BKT, CENPI, GUCY2F, MSN, PCDH11X, PLP1, VSIG4, PAK3, WAS, PCDH19, PDHA1, and SRPX2. The VSIG4 and PCDH11X genes are associated with the immune system and disease, PDHA1 is associated with biosynthetic related pathways, and PCDH19 is expressed in the nervous system and skin. These genes may be useful as candidate genes for molecular breeding.

AAV2-mediated follistatin overexpression induces ovine primary myoblasts proliferation.

BACKGROUND: Follistatin (FST) has been shown to bind to some TGF-ß family members and can function as a potent myostatin (MSTN) antagonist. Recent studies have revealed that over-expression of FST by adeno-associated viruses increases muscle growth in mice, humans and nonhuman primates. In the present study, to determine the effect of FST on ovine primary myoblast (OPM) proliferation, FST was over-expressed using an adeno-associated virus serotype 2 (AAV 2) vector. RESULTS: Western blot results showed that AAV induced the expression of FST protein in transduced OPM cells. Real-time quantitative PCR results indicated that over-expression of FST resulted in a dramatic increase in Akt I and CDK2 expression and a decrease in p21 expression. Moreover, cell cycle analysis confirmed that FST down-regulated p21, a CDK inhibitor, and increased the level of CDK2 expression in OPM cells. Hence, follistatin positively regulated the G1 to S progression. Our results showed that FST induced proliferation through a down-regulation of p21, as only the p21 expression level was down-regulated as a result of FST over-expression in myoblasts, whereas no change was observed in the level of p57 expression. CONCLUSIONS: These results expanded our understanding of the regulation mechanism of FST in ovine primary myoblasts. Our results provide the first evidence that the AAV viral system can be used for gene transfer in ovine myoblast cells. Moreover, the results showed that an AAV vector can successfully induce the expression of FST in OPM cells in vitro. These findings demonstrated that FST over-expression induces proliferation through a down-regulation of the p21 gene under proliferating conditions.

Comparison of single-trait and multiple-trait genomic prediction models.

BACKGROUND: In this study, a single-trait genomic model (STGM) is compared with a multiple-trait genomic model (MTGM) for genomic prediction using conventional estimated breeding values (EBVs) calculated using a conventional single-trait and multiple-trait linear mixed models as the response variables. Three scenarios with and without missing data were simulated; no missing data, 90% missing data in a trait with high heritability, and 90% missing data in a trait with low heritability. The simulated genome had a length of 500 cM with 5000 equally spaced single nucleotide polymorphism markers and 300 randomly distributed quantitative trait loci (QTL). The true breeding values of each trait were determined using 200 of the QTLs, and the remaining 100 QTLs were assumed to affect both the high (trait I with heritability of 0.3) and the low (trait II with heritability of 0.05) heritability traits. The genetic correlation between traits I and II was 0.5, and the residual correlation was zero. RESULTS: The results showed that when there were no missing records, MTGM and STGM gave the same reliability for the genomic predictions for trait I while, for trait II, MTGM performed better that STGM. When there were missing records for one of the two traits, MTGM performed much better than STGM. In general, the difference in reliability of genomic EBVs predicted using the EBV response variables estimated from either the multiple-trait or single-trait models was relatively small for the trait without missing data. However, for the trait with missing data, the EBV response variable obtained from the multiple-trait model gave a more reliable genomic prediction than the EBV response variable from the single-trait model. CONCLUSIONS: These results indicate that MTGM performed better than STGM for the trait with low heritability and for the trait with a limited number of records. Even when the EBV response variable was obtained using the multiple-trait model, the genomic prediction using MTGM was more reliable than the prediction using the STGM.

Analysis of 206 whole-genome resequencing reveals selection signatures associated with breed-specific traits in Hu sheep.

As an invaluable Chinese sheep germplasm resource, Hu sheep are renowned for their high fertility and beautiful wavy lambskins. Their distinctive characteristics have evolved over time through a combination of artificial and natural selection. Identifying selection signatures in Hu sheep can provide a straightforward insight into the mechanism of selection and further uncover the candidate genes associated with breed-specific traits subject to selection. Here, we conducted whole-genome resequencing on 206 Hu sheep individuals, each with an approximate 6-fold depth of coverage. And then we employed three complementary approaches, including composite likelihood ratio, integrated haplotype homozygosity score and the detection of runs of homozygosity, to detect selection signatures. In total, 10 candidate genomic regions displaying selection signatures were simultaneously identified by multiple methods, spanning 88.54 Mb. After annotating, these genomic regions harbored collectively 92 unique genes. Interestingly, 32 candidate genes associated with reproduction were distributed in nine genomic regions detected. Out of them, two stood out as star candidates: BMPR1B and GNRH2, both of which have documented associations with fertility, and a HOXA gene cluster (HOXA1-5, HOXA9, HOXA10, HOXA11 and HOXA13) had also been linked to fertility. Additionally, we identified other genes that are related to hair follicle development (LAMTOR3, EEF1A2), ear size (HOXA1, KCNQ2), fat tail formation (HOXA10, HOXA11), growth and development (FAF1, CCNDBP1, GJB2, GJA3), fat deposition (ACOXL, JAZF1, HOXA3, HOXA4, HOXA5, EBF4), immune (UBR1, FASTKD5) and feed intake (DAPP1, RNF17, NPBWR2). Our results offer novel insights into the genetic mechanisms underlying the selection of breed-specific traits in Hu sheep and provide a reference for sheep genetic improvement programs.