Shulin Yang

Genetic variation and relationships of eighteen Chinese indigenous pig breeds

Chinese indigenous pig breeds are recognized as an invaluable component of the worlds pig genetic resources and are divided traditionally into six types. Twenty-six microsatellite markers recommended by the FAO (Food and Agriculture Organization) and ISAG (International Society of Animal Genetics) were employed to analyze the genetic diversity of 18 Chinese indigenous pig breeds with 1001 individuals representing five types, and three commercial breeds with 184 individuals. The observed heterozygosity, unbiased expected heterozygosity and the observed and effective number of alleles were used to estimate the genetic variation of each indigenous breed. The unbiased expected heterozygosity ranged between 0.700 (Mashen) and 0.876 (Guanling), which implies that there is an abundant genetic variation stored in Chinese indigenous pig breeds. Breed differentiation was shown by fixation indices (FIT, FIS, and FST). The FSTper locus varied from 0.019 (S0090) to 0.170 (SW951), and the average FSTof all loci was 0.077, which means that most of the genetic variation was kept within breeds and only a little of the genetic variation exists between populations. The Neighbor-Joining tree was constructed based on the Nei DA(1978) distances and one large cluster with all local breeds but the Mashen breed, was obtained. Four smaller sub-clusters were also found, which included two to four breeds each. These results, however, did not completely agree with the traditional type of classification. A Neighbor-Joining dendrogram of individuals was established from the distance of – ln(proportions of shared alleles); 92.14% of the individuals were clustered with their own breeds, which implies that this method is useful for breed demarcation. This extensive research on pig genetic diversity in China indicates that these 18 Chinese indigenous breeds may have one common ancestor, helps us to better understand the relative distinctiveness of pig genetic resources, and will assist in developing a national plan for the conservation and utilization of Chinese indigenous pig breeds.

Highly efficient CRISPR/Cas9-mediated transgene knockin at the H11 locus in pigs

Transgenic pigs play an important role in producing higher quality food in agriculture and improving human health when used as animal models for various human diseases in biomedicine. Production of transgenic pigs, however, is a lengthy and inefficient process that hinders research using pig models. Recent applications of the CRISPR/Cas9 system for generating site-specific gene knockout/knockin models, including a knockout pig model, have significantly accelerated the animal model field. However, a knockin pig model containing a site-specific transgene insertion that can be passed on to its offspring remains lacking. Here, we describe for the first time the generation of a site-specific knockin pig model using a combination of CRISPR/Cas9 and somatic cell nuclear transfer. We also report a new genomic “safe harbor” locus, named pH11, which enables stable and robust transgene expression. Our results indicate that our CRISPR/Cas9 knockin system allows highly efficient gene insertion at the pH11 locus of up to 54% using drug selection and 6% without drug selection. We successfully inserted a gene fragment larger than 9u2009kb at the pH11 locus using the CRISPR/Cas9 system. Our data also confirm that the gene inserted into the pH11 locus is highly expressed in cells, embryos and animals.

Systematic identification and characterization of long intergenic non-coding RNAs in fetal porcine skeletal muscle development.

Long intergenic non-coding RNAs (lincRNAs) play important roles in many cellular processes. Here, we present the first systematic identification and characterization of lincRNAs in fetal porcine skeletal muscle. We obtained a total of 55.02 million 90-bp paired-end reads and assembled 54,550 transcripts using cufflinks. We developed a pipeline to identify 570 multi-exon lincRNAs by integrating a set of previous approaches. These putative porcine lincRNAs share many characteristics with mammalian lincRNAs, such as a relatively short length, small number of exons and low level of sequence conservation. We found that the porcine lincRNAs were preferentially located near genes mediating transcriptional regulation rather than those with developmental functions. We further experimentally analyzed the features of a conserved mouse lincRNA gene and found that isoforms 1 and 4 of this lincRNA were enriched in the cell nucleus and were associated with polycomb repressive complex 2 (PRC2). Our results provide a catalog of fetal porcine lincRNAs for further experimental investigation of the functions of these genes in the skeletal muscle developmental process.

Characterization, expression profiles, intracellular distribution and association analysis of porcine PNAS-4 gene with production traits.

BackgroundIn a previous screen to identify differentially expressed genes associated with embryonic development, the porcine PNAS-4 gene had been found. Considering differentially expressed genes in early stages of muscle development are potential candidate genes to improve meat quality and production efficiency, we determined how porcine PNAS-4 gene regulates meat production. Therefore, this gene has been sequenced, expression analyzed and associated with meat production traits.ResultsWe cloned the full-length cDNA of porcine PNAS-4 gene encoding a protein of 194 amino acids which was expressed in the Golgi complex. This gene was mapped to chromosome 10, q11–16, in a region of conserved synteny with human chromosome 1 where the human homologous gene was localized. Real-time PCR revealed that PNAS-4 mRNA was widely expressed with highest expression levels in skeletal muscle followed by lymph, liver and other tissues, and showed a down-regulated expression pattern during prenatal development while a up-regulated expression pattern after weaning. Association analysis revealed that allele C of SNP A1813C was prevalent in Chinese indigenous breeds whereas A was dominant allele in Landrace and Large White, and the pigs with homozygous CC had a higher fat content than those of the pigs with other genotypes (P < 0.05).ConclusionPorcine PNAS-4 protein tagged with green fluorescent protein accumulated in the Golgi complex, and its mRNA showed a widespread expression across many tissues and organs in pigs. It may be an important factor affecting the meat production efficiency, because its down-regulated expression pattern during early embryogenesis suggests involvement in increase of muscle fiber number. In addition, the SNP A1813C associated with fat traits might be a genetic marker for molecular-assisted selection in animal breeding.

Characterization of porcine MMP-2 and its association with immune traits

Matrix metalloproteinase-2 (MMP-2) plays important roles in inflammation and immunity besides its basic role in degrading and remodelling extracellular matrix (ECM). The expression of MMP-2 is up-regulated in many human as well as animal models of inflammatory and immune diseases. In this study, we cloned the 5-upstream sequence, 3-downstream sequence as well as other missed genomic sequences of porcine MMP-2, the genomic structure and the promotor sequence were analyzed and found to share high similarity with those of human MMP-2. Porcine MMP-2 was assigned to SSC6p14-p15, and closely linked to microsatellite SW1108 (53cR, LOD score 7.59) by IMpRH panel. Real-time PCR analysis revealed that the expression of porcine MMP-2 was remarkably different in diverse tissues, a high level expression was observed in the testis and uterus, relatively low expression in other tissues. Allele frequencies determination in different pig breeds and association study were performed on the selected SNP and indel. The results showed that the SNP AcyI in exon 12 was significantly associated with white blood cell count (WBC) of neonate piglets at 0 day (P=0.0079), and classical swine fever virus antibody level (CSFV-AB) of pigs at 17 days (P=0.0461), the indel MspI in intron 4 had remarkable correlation with mean corpuscular hemoglobin (MCH) of pigs at 17 days (P<0.0001).

Hyperinsulinemia shifted energy supply from glucose to ketone bodies in early nonalcoholic steatohepatitis from high-fat high-sucrose diet induced Bama minipigs.

The minipig can serve as a good pharmacological model for human subjects. However, the long-term pathogenesis of high-calorie diet-induced metabolic syndromes, including NASH, has not been well described in minipigs. We examined the development of metabolic syndromes in Bama minipigs that were fed a high-fat, high-sucrose diet (HFHSD) for 23 months, by using histology and serum biochemistry and by profiling the gene expression patterns in the livers of HFHSD pigs compared to controls. The pathology findings revealed microvesicular steatosis, iron overload, arachidonic acid synthesis, lipid peroxidation, reduced antioxidant capacity, increased cellular damage, and inflammation in the liver. RNA-seq analysis revealed that 164 genes were differentially expressed between the livers of the HFHSD and control groups. The pathogenesis of early-stage NASH was characterized by hyperinsulinemia and by de novo synthesis of fatty acids and nascent triglycerides, which were deposited as lipid droplets in hepatocytes. Hyperinsulinemia shifted the energy supply from glucose to ketone bodies, and the high ketone body concentration induced the overexpression of cytochrome P450 2E1 (CYP2E1). The iron overload, CYP2E1 and alcohol dehydrogenase 4 overexpression promoted reactive oxygen species (ROS) production, which resulted in arachidonic and linoleic acid peroxidation and, in turn, led to malondialdehyde production and a cellular response to ROS-mediated DNA damage.

Identification and analysis of pig chimeric mRNAs using RNA sequencing data

BackgroundGene fusion is ubiquitous over the course of evolution. It is expected to increase the diversity and complexity of transcriptomes and proteomes through chimeric sequence segments or altered regulation. However, chimeric mRNAs in pigs remain unclear. Here we identified some chimeric mRNAs in pigs and analyzed the expression of them across individuals and breeds using RNA-sequencing data.ResultsThe present study identified 669 putative chimeric mRNAs in pigs, of which 251 chimeric candidates were detected in a set of RNA-sequencing data. The 618 candidates had clear trans-splicing sites, 537 of which obeyed the canonical GU-AG splice rule. Only two putative pig chimera variants whose fusion junction was overlapped with that of a known human chimeric mRNA were found. A set of unique chimeric events were considered middle variances in the expression across individuals and breeds, and revealed non-significant variance between sexes. Furthermore, the genomic region of the 5′ partner gene shares a similar DNA sequence with that of the 3′ partner gene for 458 putative chimeric mRNAs. The 81 of those shared DNA sequences significantly matched the known DNA-binding motifs in the JASPAR CORE database. Four DNA motifs shared in parental genomic regions had significant similarity with known human CTCF binding sites.ConclusionsThe present study provided detailed information on some pig chimeric mRNAs. We proposed a model that trans-acting factors, such as CTCF, induced the spatial organisation of parental genes to the same transcriptional factory so that parental genes were coordinatively transcribed to give birth to chimeric mRNAs.

Identification, bioinformatic analysis and expression profiling of candidate mRNA-like non-coding RNAs in Sus scrofa

Messenger RNA-like non-coding RNAs (mlncRNAs) are a newly identified group of non-coding RNAs (ncRNAs) that may be involved in a number of critical cellular events. In this study, 93 candidate porcine mlncRNAs were obtained by computational prediction and screening, among which 72 were mapped to the porcine genome. Further analysis of 8 representative candidates revealed that these mlncRNA candidates are not highly conserved among species. Remarkably, one of the candidates, sTF35495, was found to be precursor of a putative porcine microRNA. By RACE PCR, we determined that the full length of sTF35495 was 3 kb. The protein-coding potential of this RNA was tested in silico with no significant finding. Semi-quantitative RT-PCR analysis of the subgroup of 8 candidates revealed two distinct expression profiles and two molecules were further validated by real-time PCR. The predicted pre-microRNA sequence in this study provides a potentially interesting insight into the in vivo function of porcine mlncRNAs and our findings suggest that they play key biological roles in Sus scrofa.

Molecular Cloning, Mapping, and Expression Analysis of the EIF4A2 Gene in Pig

A full-length cDNA clone encoding the eukaryotic translation initiation factor 4A, isoform 2 (EIF4A2), was cloned from the fetal skeletal cDNA library from the pig (Sus scrofa). EIF4A2 is a highly conserved gene for one of the protein-synthesis initiation factors involved in the binding of mRNA to the ribosome. Based on this cDNA sequence, the deduced protein of 407 amino acids contains the characteristic motifs shared by the DEAD-box supergene family. The genomic nucleotide sequence of this gene was determined and a single nucleotide polymorphism located in the 5′ untranslated region was genotyped. The porcine EIF4A2 was expressed in all tissues examined but in variable amounts. The EIF4A2 expression level in muscle was upregulated through embryonic and neonatal development until adult, suggesting that porcine EIF4A2 was possibly involved in translation regulation of other muscle-related genes in muscle formation and development. In addition, we mapped porcine EIF4A2 to q4.1 of SSC13, in agreement with comparative mapping data.

Dynamic Expression of MicroRNA-127 During Porcine Prenatal and Postnatal Skeletal Muscle Development

MicroRNAs (miRNAs), evolutionarily conserved non-coding RNAs in length 21-24 bp, play a critical role in skeletal muscle development. In this study, to explore the function of mircoRNA-127 in porcine skeletal muscle development, eight tissue samples from adult pigs and longissimus muscle samples at 26 developmental stages were collected from Tongcheng and Landrace pigs. The spatial-temporal expression profiles of miRNA-127 were carried out using step-loop quantitative real-time PCR (stem-loop RT-PCR). To explore the molecular functions of miRNA-127, we predicted its target genes and performed functional annotation using bioinformatics methods. Results suggested that miRNA-127 was abundantly expressed in heart, ovary, uterus and spleen tissues and was weakly expressed in liver, lung, kidney and small intestine in both Tongcheng and Landrace pigs. And miRNA-127 showed significant expression differences in heart, ovary, spleen and uterus tissues between these two breeds. miRNA-127 basically kept at a relatively stable high level in middle and later embryonic stages and a low expression level in early embryonic stages and postnatal stages, but the expression levels of miRNA-127 were higher in Tongcheng pigs than in Landrace at most developmental stages. miRNA-127 potentially regulated 240 candidate genes. Results of Gene Ontology and KEGG pathway analysis indicated that these genes could be involved in many molecular functions and mechanisms, such as regulation of the force of heart contraction, regulation of transcription, regulation of T cell differentiation, MAPK signaling pathway and GnRH signaling pathway. Many significantly enriched GO terms and KEGG pathways were related to skeletal muscle development. This study will be helpful to understand the biological function for miRNA-127 and identify candidate gene associated with meat production traits in pigs.