The mutation rate used in the previous analyses of pig evolution and demographics was cursory and hence invited potential bias in inferring evolutionary history. Herein, we estimated the de novo mutation rate of pigs as 3.6 × 10−9 per base per generation using high-quality whole-genome sequencing data from nine individuals in a three-generation pedigree through stringent filtering and validation. Using this mutation rate, we re-investigated the evolutionary history of pigs. The estimated divergence time of ∼ 10 kiloyears ago (KYA) between European wild and domesticated pigs was consistent with the domestication time of European pigs based on archaeological evidence. However, other divergence events inferred here were not as ancient as previously described. Our estimates suggest that Sus speciation occurred ∼ 1.36 million years ago (MYA); European wild pigs split from Asian wild pigs only ∼ 219 KYA; and south and north Chinese wild pigs split ∼ 25 KYA. Meanwhile, our results showed that the most recent divergence event between Chinese wild and domesticated pigs occurred in the Hetao Plain, northern China, approximately 20 KYA, supporting the possibly independent domestication in northern China along the middle Yellow River. We also found that the maximum effective population size of pigs was ∼ 6 times larger than estimated before. An archaic migration from other Sus species originating ∼ 2 MYA to European pigs was detected during western colonization of pigs, which may affect the accuracy of previous demographic inference. Our de novo mutation rate estimation and its consequences for demographic history inference reasonably provide a new vision regarding the evolutionary history of pigs.

Mucopolysaccharidosis type IIIB (MPS IIIB) is a rare genetic disorder caused by loss-of-function mutations in the NAGLU gene. Pigs are an ideal large animal model for human diseases; however, a porcine model of MPS IIIB has not been reported. We have previously generated a heterozygous NAGLU-deficient (NAGLU^+/-) Large White boar via a transgenic approach. Here we characterized phenotypes of the F₁ offspring of this founder to establish a pig model for MPS IIIB. qRT-PCR revealed that the NAGLU expression level was significantly decreased in a variety of tissues in NAGLU^+/− pigs. ELISA assays showed obvious deficiency of NAGLU and higher (P<0.05) glycosaminoglycan levels in multiple tissues from NAGLU^+/− pigs. NAGLU^+/− pigs grew at a significantly (P<0.05) slower rate than control animals (NAGLU^+/+). Death, mostly sudden death, occurred at all ages in NAGLU^+/− pigs, most of which died within two years. Necropsy findings included pleural adhesions, lung shrinkage and abnormalities in the pericardium and mild hepatomegaly in NAGLU^+/− pigs. Notable pathological changes were observed in the sections of brain, liver, spleen and kidney from NAGLU^+/− pigs. Brain atrophy, ventriculomegaly, cerebellar atrophy and abnormalities in the intracerebral capsule, parietal lobes and the thalamus were also evident in NAGLU^+/− pigs. Together, NAGLU^+/− pigs show typical symptoms of human MPS IIIB patients and thus represent a novel large animal model for the disease. This article has an associated First Person interview with the first author of the paper.

In an attempt to generate g.A746G substitution in the BMPR-IB gene, we unexpectedly obtained BMPR-IB homozygous knockout piglets (BMPR-IB^-/-) and heterogeneous knockout piglets with one copy of the A746G mutation (BMPR-IB^-/746G) via CRISPR/Cas9 editing. Polymerase chain reaction (PCR) and sequencing revealed complex genomic rearrangements in the target region. All BMPR-IB-disrupted piglets showed an inability to stand and walk normally. Both BMPR-IB^-/- and BMPR-IB^-/746G piglets exhibited severe skeletal dysplasia characterized by distorted and truncated forearms (ulna, radius) and disordered carpal, metacarpal, and phalangeal bones in the forelimbs. The piglets displayed more severe deformities in the hindlimbs by visual inspection, including fibular hemimelia, enlarged tarsal bone, and disordered toe joint bones. Limb deformities were more profound in BMPR-IB^-/- piglets than in the BMPR-IB^-/746G piglets. Proteomic analysis identified 139 differentially expressed proteins (DEPs) in the hindlimb fibula of BMPR-IB^-/746G piglets compared to the wild-type (WT) controls. Most DEPs are involved in skeletal or embryonic development and/or the TGF-β pathway and tumor progression. Gene Ontology (GO) and protein domain enrichment analysis suggested alterations in these processes. Of the top 50 DEPs, a large proportion, e.g., C1QA, MYO1H, SRSF1, P3H1, GJA1, TCOF1, RBM10, SPP2, MMP13, and PHAX, were significantly associated with skeletal development. Our study provides novel findings on the role of BMPR-IB in mammalian limb development.

P53 is the most frequently mutated gene in human cancers, with over half of all tumors harboring mutation at this locus. R248 and R249 (corresponding to porcine R241 and R242), are among the hotspot mutations frequently mutated in liver, lung, breast, and some other cancers. In this study, p53 gene was knocked out or point-edited (R241 and R242 were converted to 241W and 242S) in porcine fetal fibroblast (PFF) cells via CRISPR-Cas9 technique. High throughput sequencing of miRNA and mRNA uncovered a total of 225 differentially expressed miRNAs (DEMs) and 738 differentially expressed genes (DEGs) in the p53 knockout (p53-KO) cells, and a total of 211 DEMs and 722 DEGs in the point-modified (p53-241W242S) cells. Totally 28 annotated DEMs were found to overlap between p53-KO/p53-WT and p53-241W242S/p53-WT miRNAs datasets, of which miR-34 c, miR-218, miR-205, miR-105-1, miR-105-2, miR-206, miR-224 and miR-429 play important roles in p53 regulatory network. Among the top 10 DEGs in p53-KO and p53-241W242S cells, most genes were reported to be involved in tumors, cell proliferation or cell migration. p53-KO and p53-241W242S cells showed a significantly higher (P < 0.01) proliferation rate compared with p53-WT cells. In conclusion, genetic modifications of p53 gene significantly affect the expression levels of a large number of genes and miRNAs in the PFF cells. The p53-edited PFF cells could be used as non-tumor cell models for investigating the p53 signaling network, and as donor cells for somatic nuclear transfer, with the aim to develop porcine models with the corresponding p53 mutations. Abbreviations: CRISPR-Cas9: Clustered regularly interspaced short palindromic repeats-associated protein 9; PFF: porcine fetal fibroblasts; SCNT: somatic cell nuclear transfer; RNA sequencing: small RNA sequencing and mRNA sequencing; DEGs: differentially expressed mRNAs; DEMs: differentially expressed miRNAs.

Cytosine base editors (CBEs) and adenine base editors (ABEs) are recently developed CRISPR-mediated genome-editing tools that do not introduce double-strand breaks. In this study, five ABEs, ABE7.10, ABEmax, NG-ABEmax, ABE8e and NG-ABE8e, were used to generate A-to-G (T-to-C) conversions in five genome loci in porcine fetal fibroblasts (PFFs). Variable yet appreciable editing efficiencies and variable activity windows were observed in these targeting regions via these five editors. The strategy of two sgRNAs in one vector exhibited superior editing efficiency to that of using two separate sgRNA expression vectors. ABE-mediated start-codon mutation in APOE silenced its expression of protein and, unexpectedly, eliminated the vast majority of its mRNA. No off-target DNA site was detected for these editors. Substantial off-target RNA events were present in the ABE-edited cells, but no KEGG pathway was found to be significantly enriched. Our study supports that ABEs are powerful tools for A-to-G (T-to-C) point-mutation modification in porcine cells.