Long Noncoding RNA Lnc-SEMT Modulates IGF2 Expression by Sponging miR-125b to Promote Sheep Muscle Development and Growth.

BACKGROUND/AIMS: Long noncoding RNAs (lncRNAs) are RNA transcripts that are more than 200 nt long but have little protein-coding potential. Within the last few years, thousands of lncRNAs have been identified and their functions in biological processes have begun to be understood. Although many studies havebegun to examine the functions of many noncoding RNAs, very little is known about the functions of long noncoding (lncRNA) function of livestock production and molecular mechanisms of their functions are still lackingrelated to livestock production. METHODS: Expression of sheep enhanced muscularityTranscript lncRNA (lnc-SEMT) and miR-125b were examined in sheep using quantitative reverse-transcription polymerase chain reaction. Expression of Myod (myogenic determination factor), Myog (myoglobin) and Insulin-like growth factor 2 (IGF2)were examined by Western Blot.Luciferase reporter assays were performedto confirm the relationship between lnc-SEMT and miR-125b. RESULTS: Here, we identified a novel lnc-SEMT that promote sheep myoblast differentiation in vitro and enhanced sheep muscularity in vivo. Functional analyses showed that lnc-SEMT accelerates sheep myoblast differentiation in vitro. lnc-SEMT transgenic sheep exhibit a muscle hypertrophy phenotype characterized by increased body weight, and increased the number of muscle fibers indicating that lnc-SEMT play an important role in the regulation of skeletal muscle differentiation in vivo. Our results show that lnc-SEMT acts as a molecular sponge by antagonizing miR-125b to control IGF2 protein labundance in vitro and in vivo. CONCLUSION: In brief, lnc-SEMT is the first example of a lncRNA could be a useful candidate for improving biological growth traits such as skeletal muscle production in sheep.

Identification and Characterization of Long Noncoding RNAs in Ovine Skeletal Muscle.

Long noncoding RNAs (lncRNAs) are increasingly being recognized as key regulators in many cellular processes. However, few reports of them in livestock have been published. Here, we describe the identification and characterization of lncRNAs in ovine skeletal muscle. Eight libraries were constructed from the gastrocnemius muscle of fetal (days 85 and 120), newborn and adult Texel and Ujumqin sheep. The 2002 identified transcripts shared some characteristics, such as their number of exons, length and distribution. We also identified some coding genes near these lncRNA transcripts, which are particularly associated with transcriptional regulation- and development-related processes, suggesting that the lncRNAs are associated with muscle development. In addition, in pairwise comparisons between the libraries of the same stage in different breeds, a total of 967 transcripts were differentially expressed but just 15 differentially expressed lncRNAs were common to all stages. Among them, we found that TCONS_00013201 exhibited higher expression in Ujumqin samples, while TCONS_00006187 and TCONS_00083104 were higher in Texel samples. Moreover, TCONS_00044801, TCONS_00008482 and TCONS_00102859 were almost completely absent from Ujumqin samples. Our results suggest that differences in the expression of these lncRNAs may be associated with the muscular differences observed between Texel and Ujumqin sheep breeds.

MiR-378 and BMP-Smad can influence the proliferation of sheep myoblast.

MicroRNA (miRNA) is a sort of endogenous ~20-25 nt non-coding RNAs, and it can regulate a variety of biological events. We found the miR-378 may involve in regulating the muscle development of sheep during our previous research. However, the molecular mechanism of miR-378 regulating myoblast proliferation is still unclear. In this research, we predicted that BMP2 (Bone morphogenetic protein 2) was the target gene of miR-378 and the BMP-Smad signal pathway that BMP2 participated in playing an important role in the muscle development. Therefore, we tried to determine whether miR-378 influence myoblast proliferation of sheep through the BMP-Smad signal pathway. The results indicated that inhibit BMP-Smad signal pathway by interfering Smad4 to promote proliferation of sheep myoblasts; promote BMP-Smad signal pathway by interfering Smad7 to inhibit proliferation of sheep myoblasts; over-expression miR-378 promotes BMP-Smad signal pathway and myoblast proliferation in sheep; interfering miR-378 inhibits BMP-Smad signal pathway and myoblast proliferation in sheep. However, when both of which functioned at the myoblast, miR-378 could not fully depend on BMP-Smad signal pathway to regulate myoblast proliferation. In sum, both miR-378 and BMP-Smad can influence the proliferation of myoblast, but miR-378 does not target the 3' UTR of sheep BMP2.

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.

Identification of MEF2B and TRHDE Gene Polymorphisms Related to Growth Traits in a New Ujumqin Sheep Population.

2 SNPs were discovered in our previous genome-wide association study (GWAS): s58995.1 (rs420767326 A>G) in MEF2B gene and OAR3_115712045.1 (rs401775061 A>C) in TRHDE gene, which were significantly associated with post-weaning gain in sheep. Herein, we performed a replication experiment to investigate single nucleotide polymorphisms (SNPs) within the MEF2B and TRHDE gene exons, the 5'untranslated regions (within 1000bp), the 3' untranslated regions (within 1000bp) and their associations with Ujumqin sheep growth traits in 4-month age and 6-month age, respectively. Finally,3 SNPs were selected to be investigated including 1 SNP in 3'untranslated regions in MEF2B gene (rs417014745 A>G) and 2 SNPs in TRHDE gene (rs426980328 T>C and rs430810656 G>A).The χ2 test showed all the 3 variations were in Hardy-Weinberg equilibrium (P>0.05) status. Association analysis suggested that rs426980328 T>C was significantly associated with body weight and chest girth in 4-month age (P<0.05). rs430810656 G>A exhibited extremely significant association with body weight and chest girth in 4-month age (P<0.01). rs417014745 A>G was extremely significantly associated with body weight and chest girth in 4-month age and chest girth in 6-month age (P<0.01), and it was also significantly associated with body weight in 6-month age (P<0.05). Combined effect analysis indicated significant associations between the combinations of rs426980328-rs417014745, rs430810656-rs417014745 and several growth traits (P<0.05). These results suggested MEF2B and TRHDE genes affected growth traits in Ujumqin sheep and the combination effect of the two genes also played a significant effective role. These SNPs might have potential value as genetic markers for growth traits and it could be used in Ujumqin sheep breeding in future. Further studies are necessary to confirm our findings.

Genome-wide analysis reveals adaptation to high altitudes in Tibetan sheep.

Tibetan sheep have lived on the Tibetan Plateau for thousands of years; however, the process and consequences of adaptation to this extreme environment have not been elucidated for important livestock such as sheep. Here, seven sheep breeds, representing both highland and lowland breeds from different areas of China, were genotyped for a genome-wide collection of single-nucleotide polymorphisms (SNPs). The FST and XP-EHH approaches were used to identify regions harbouring local positive selection between these highland and lowland breeds, and 236 genes were identified. We detected selection events spanning genes involved in angiogenesis, energy production and erythropoiesis. In particular, several candidate genes were associated with high-altitude hypoxia, including EPAS1, CRYAA, LONP1, NF1, DPP4, SOD1, PPARG and SOCS2. EPAS1 plays a crucial role in hypoxia adaption; therefore, we investigated the exon sequences of EPAS1 and identified 12 mutations. Analysis of the relationship between blood-related phenotypes and EPAS1 genotypes in additional highland sheep revealed that a homozygous mutation at a relatively conserved site in the EPAS1 3' untranslated region was associated with increased mean corpuscular haemoglobin concentration and mean corpuscular volume. Taken together, our results provide evidence of the genetic diversity of highland sheep and indicate potential high-altitude hypoxia adaptation mechanisms, including the role of EPAS1 in adaptation.

Rosa26-targeted sheep gene knock-in via CRISPR-Cas9 system.

Recent advances in our ability to design DNA binding factors with specificity for desired sequences have resulted in a revolution in genetic engineering, enabling directed changes to the genome to be made relatively easily. Technologies that facilitate specific and precise genome editing, such as knock-in, are critical for determining the functions of genes and for understanding fundamental biological processes. The CRISPR/Cas9 system has recently emerged as a powerful tool for functional genomic studies in mammals. Rosa26 gene can encode a non-essential nuclear RNA in almost all organizations, and become a hot point of exogenous gene insertion. Here, we describe efficient, precise CRISPR/Cas9-mediated Integration using a donor vector with tGFP sequence targeted in the sheep genomic Rosa26 locus. We succeeded in integrating with high efficiency an exogenous tGFP (turboGFP) gene into targeted genes in frame. Due to its simplicity, design flexibility, and high efficiency, we propose that CRISPR/Cas9-mediated knock-in will become a standard method for the generation transgenic sheep.

DNA methylation Landscape of body size variation in sheep.

Sub-populations of Chinese Mongolian sheep exhibit significant variance in body mass. In the present study, we sequenced the whole genome DNA methylation in these breeds to detect whether DNA methylation plays a role in determining the body mass of sheep by Methylated DNA immunoprecipitation - sequencing method. A high quality methylation map of Chinese Mongolian sheep was obtained in this study. We identified 399 different methylated regions located in 93 human orthologs, which were previously reported as body size related genes in human genome-wide association studies. We tested three regions in LTBP1, and DNA methylation of two CpG sites showed significant correlation with its RNA expression. Additionally, a particular set of differentially methylated windows enriched in the "development process" (GO: 0032502) was identified as potential candidates for association with body mass variation. Next, we validated small part of these windows in 5 genes; DNA methylation of SMAD1, TSC1 and AKT1 showed significant difference across breeds, and six CpG were significantly correlated with RNA expression. Interestingly, two CpG sites showed significant correlation with TSC1 protein expression. This study provides a thorough understanding of body size variation in sheep from an epigenetic perspective.

Genome-Wide Specific Selection in Three Domestic Sheep Breeds.

BACKGROUND: Commercial sheep raised for mutton grow faster than traditional Chinese sheep breeds. Here, we aimed to evaluate genetic selection among three different types of sheep breed: two well-known commercial mutton breeds and one indigenous Chinese breed. RESULTS: We first combined locus-specific branch lengths and di statistical methods to detect candidate regions targeted by selection in the three different populations. The results showed that the genetic distances reached at least medium divergence for each pairwise combination. We found these two methods were highly correlated, and identified many growth-related candidate genes undergoing artificial selection. For production traits, APOBR and FTO are associated with body mass index. For meat traits, ALDOA, STK32B and FAM190A are related to marbling. For reproduction traits, CCNB2 and SLC8A3 affect oocyte development. We also found two well-known genes, GHR (which affects meat production and quality) and EDAR (associated with hair thickness) were associated with German mutton merino sheep. Furthermore, four genes (POL, RPL7, MSL1 and SHISA9) were associated with pre-weaning gain in our previous genome-wide association study. CONCLUSIONS: Our results indicated that combine locus-specific branch lengths and di statistical approaches can reduce the searching ranges for specific selection. And we got many credible candidate genes which not only confirm the results of previous reports, but also provide a suite of novel candidate genes in defined breeds to guide hybridization breeding.

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.

Signals of Ezh2, Src, and Akt Involve in myostatin-Pax7 pathways regulating the myogenic fate determination during the sheep myoblast proliferation and differentiation.

Myostatin and Pax7 have been well documented individually, however, the mechanism by which Myostatin regulates Pax7 is seldom reported. Here, based on muscle transcriptome analysis in Texel (Myostatin mutant) and Ujumqin (wild type) sheep across the five fetal stages, we constructed and examined the Myostatin-Pax7 pathways in muscle. Then we validated the signals by RNAi in the proliferating and differentiating sheep myoblasts in vitro at mRNA, protein, and cell morphological levels. We reveal that Myostatin signals to Pax7 at least through Ezh2, Src, and Akt during the sheep myoblast proliferation and differentiation. Other signals such as p38MAPK, mTOR, Erk1/2, Wnt, Bmp2, Smad, Tgfb1, and p21 are most probably involved in the Myostatin-affected myogenic events. Myostatin knockdown significantly reduces the counts of nucleus and myotube, but not the fusion index of myoblasts during cell differentiation. In addition, findings also indicate that Myostatin is required for normal myogenic differentiation of the sheep myoblasts, which is different from the C2C12 myoblasts. We expand the regulatory network of Myostatin-Pax7 pathways and first illustrate that Myostatin as a global regulator participates in the epigenetic events involved in myogenesis, which contributes to understand the molecular mechanism of Myostatin in regulation of myogenesis.