Lixin Du

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.

Expression patterns of peroxisome proliferator-activated receptor gamma 1 versus gamma 2, and their association with intramuscular fat in goat tissues

Intramuscular fat (IMF) shortage causes the lack of juiciness and tenderness of goat meat, while peroxisome proliferator-activated receptor gamma 1 (PPARγ1) and gamma 2 (PPARγ2) play key roles in lipid metabolism. Nevertheless, their expression patterns and the relationship with IMF have been poorly exposed. Using quantitative polymerase chain reaction (qPCR), classical Soxhlet extraction, and in situ hybridization, we demonstrated that among 13 goat tissues, expression of PPARγ1 was dramatically higher than that of PPARγ2 except for lung. We further demonstrated the expression patterns of PPARγ1 and PPARγ2 and their negative association with intramuscular fat content in three goat muscles with kids growing. Meanwhile, PPARγ expression was located in the connective tissues. These results suggest that PPARγ1 is rather active for most tissues of goat, and closely related with the muscular fat metabolism during early postnatal life, but a more direct proof remains to be provided.

Genome-wide detection of CNVs in Chinese indigenous sheep with different types of tails using ovine high-density 600K SNP arrays.

Chinese indigenous sheep can be classified into three types based on tail morphology: fat-tailed, fat-rumped, and thin-tailed sheep, of which the typical breeds are large-tailed Han sheep, Altay sheep, and Tibetan sheep, respectively. To unravel the genetic mechanisms underlying the phenotypic differences among Chinese indigenous sheep with tails of three different types, we used ovine high-density 600K SNP arrays to detect genome-wide copy number variation (CNV). In large-tailed Han sheep, Altay sheep, and Tibetan sheep, 371, 301, and 66 CNV regions (CNVRs) with lengths of 71.35 Mb, 51.65 Mb, and 10.56 Mb, respectively, were identified on autosomal chromosomes. Ten CNVRs were randomly chosen for confirmation, of which eight were successfully validated. The detected CNVRs harboured 3130 genes, including genes associated with fat deposition, such as PPARA, RXRA, KLF11, ADD1, FASN, PPP1CA, PDGFA, and PEX6. Moreover, multilevel bioinformatics analyses of the detected candidate genes were significantly enriched for involvement in fat deposition, GTPase regulator, and peptide receptor activities. This is the first high-resolution sheep CNV map for Chinese indigenous sheep breeds with three types of tails. Our results provide valuable information that will support investigations of genomic structural variation underlying traits of interest in sheep.

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.

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.

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

Optimal Design for Marker-assisted Gene Pyramiding in Cross Population.

Marker-assisted gene pyramiding aims to produce individuals with superior economic traits according to the optimal breeding scheme which involves selecting a series of favorite target alleles after cross of base populations and pyramiding them into a single genotype. Inspired by the science of evolutionary computation, we used the metaphor of hill-climbing to model the dynamic behavior of gene pyramiding. In consideration of the traditional cross program of animals along with the features of animal segregating populations, four types of cross programs and two types of selection strategies for gene pyramiding are performed from a practical perspective. Two population cross for pyramiding two genes (denoted II), three population cascading cross for pyramiding three genes(denoted III), four population symmetry (denoted IIII-S) and cascading cross for pyramiding four genes (denoted IIII-C), and various schemes (denoted cross program-A–E) are designed for each cross program given different levels of initial favorite allele frequencies, base population sizes and trait heritabilities. The process of gene pyramiding breeding for various schemes are simulated and compared based on the population hamming distance, average superior genotype frequencies and average phenotypic values. By simulation, the results show that the larger base population size and the higher the initial favorite allele frequency the higher the efficiency of gene pyramiding. Parents cross order is shown to be the most important factor in a cascading cross, but has no significant influence on the symmetric cross. The results also show that genotypic selection strategy is superior to phenotypic selection in accelerating gene pyramiding. Moreover, the method and corresponding software was used to compare different cross schemes and selection strategies.

Double-Muscled Phenotype in Mutant Sheep Directed by the CRISPRCas9 System

Myostatin (MSTN) is a well-known negative regulator of muscle growth. The double-muscled sheep caused by natural loss-of-function mutations of MSTN have very strong skeletal muscle. In this study, our results demonstrate the successful generation of MSTN mutant sheep via specific targeting of an exon 1 site using Cas9 technology. The MSTN-knockout sheep in our study had increased muscle significantly just like double-muscled phenotype. Our study suggests that the direct injection of Cas9: sgRNA into zygotes could be widely used to create gene knockouts in large domestic animals. Notably, on the basis of our findings, sheep can be added to the growing list of species for which genome editing is now practical. The generation of MSTN mutant sheep has implications for the genetic improvement of local sheep varieties, and also for the usage of sheep as a model for large animal medical research.

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.