Exploring the genetic signature of body size in Yucatan miniature pig

Parental Phasing Study Identified Lineage-Specific Variants Associated with Gene Expression and Epigenetic Modifications in European-Chinese Hybrid Pigs

Understanding how hybrids integrate lineage-specific regulatory variants at the haplotype level is crucial for elucidating the genetic basis of heterosis in livestock. In this study, we established three crossbred pig families derived from distant genetic lineages and systematically identified variants from different lineages, including single nucleotide polymorphisms (SNPs) and structural variations (SVs). At the phase level, we quantitatively analyzed gene expression, four histone modifications (H3K4me3, H3K27ac, H3K4me1, and H3K27me3), and the binding strength of transcription factor (CTCF) in backfat (BF) and longissimus dorsi (LD) muscle. By colocalization analysis of phased genetic variants with phased gene expression levels and with phased epigenetic modifications, we identified 18,670 expression quantitative trait loci (eQTL) (FDR < 0.05) and 8,652 epigenetic modification quantitative trait loci (epiQTL) (FDR < 0.05). The integration of eQTL and epiQTL allowed us to explore the potential regulatory mechanisms by which lineage-specific genetic variants simultaneously influence gene expression and epigenetic modifications. For example, we identified a Large White lineage-specific duplication (DUP) encompassing the KIT gene that was significantly associated with its promoter activity (FDR = 7.83 × 10-4) and expression levels (FDR = 9.03 × 10-4). Additionally, we found that a Duroc lineage-specific SNP located upstream of AMIGO2 was significantly associated with a Duroc-specific H3K27ac peak (FDR = 0.035) and also showed a significant association with AMIGO2 expression levels (FDR = 5.12 × 10-4). These findings underscore the importance of phased regulatory variants in shaping lineage-specific transcriptional programs and highlight how the haplotype-resolved integration of eQTL and epigenetic signals can reveal the mechanistic underpinnings of hybrid regulatory architecture. Our results offer insights for molecular marker development in precision pig breeding.

Development of a whole-exome sequencing kit to facilitate porcine biomedical research

BACKGROUND

It is important for porcine models to replicate gene mutations present in human diseases to improve the translatability of animal studies. In this study, the high efficacy of a whole exome sequencing kit was demonstrated for the improved pig reference genome (Sus scrofa 11.1) to profile biomedically relevant swine breeds and enable high-depth sequencing required for intratumor heterogeneity profiling.

RESULTS

We identify a total of 751,624 single nucleotide variants (SNVs) and 113,597 insertions and deletions (INDELs) across 93 samples from 12 porcine breeds. The identified mutations and affected pathways are correlated to muscle-to-fat ratios between different porcine breeds and further inform their utility as models of obesity and cardiovascular disease. Finally, 7935 SNVs and 358 INDELs are present in an Oncopig hepatocellular carcinoma (HCC) cell line and samples from a single Oncopig HCC tumor, with pathways related to hepatic fibrosis, WNT/B-catenin, ATM signaling, and p53 signaling enriched.

CONCLUSIONS

These results demonstrate the kit's high efficacy and utility for identifying mutations in the context of obesity, cardiovascular disease, and cancer across a range of pig models used in biomedical research.

Genomic insights into the mechanisms of body size evolution in Serpentes

BACKGROUND

Body size is a critical trait that influences an animal's physiology, behavior, and ecology. However, the molecular mechanisms underlying its evolution remain poorly understood, particularly in snakes. Snakes exhibit an extremely wide range of body sizes and strong ecological adaptability. Among snake species, the maximum body mass exceeds the minimum by over 200,000-fold, while the maximum body length surpasses the minimum by more than 110-fold.

RESULTS

Through phylogenomic and comparative genomic analyses of 26 snake genomes, we identified 77 body size-associated genes (BSAGs) related to body length or body mass, highlighting key genetic drivers of body size evolution. Functional enrichment analyses revealed that metabolic pathways, particularly fatty acid metabolism and oxidoreductase activity, underwent significant expansion and positive selection, suggesting metabolic adaptations crucial for meeting the energetic demands of increased body size. Immune system-related genes, including those involved in antigen processing and presentation, similarly showed signatures of expansion and adaptive evolution, highlighting strengthened immune defenses in large-bodied snakes. Additionally, key candidate genes, such as YAP1, PLAG1, MGAT1 and SPRY1, exhibited both strong selection signals and correlation signals, and are functionally involved in developmental pathways critical for growth regulation.

CONCLUSIONS

Our findings reveal a complex interplay of sensory, immune, metabolic, and growth-related genetic adaptations driving large body size evolution in snakes. This study provides novel insights into the molecular underpinnings of snake body size diversification and advances our understanding of their evolutionary history.

Multi-tool copy number detection highlights common body size-associated variants in miniature pig breeds from different geographical regions

BACKGROUND

Copy number variations (CNVs) represent a common and highly specific type of variation in the genome, potentially influencing genetic diversity and mammalian phenotypic development. Structural variants, such as deletions, duplications, and insertions, have frequently been highlighted as key factors influencing traits in high-production pigs. However, comprehensive CNV analyses in miniature pig breeds are limited despite their value in biomedical research.

RESULTS

This study performed whole-genome sequencing in 36 miniature pigs from nine breeds from America, Asia and Oceania, and Europe. By employing a multi-tool approach (CNVpytor, Delly, GATK gCNV, Smoove), the accuracy of CNV identification was improved. In total, 34 homozygous CNVs overlapped with exonic regions in all samples, suggesting a role in expressing specific phenotypes such as uniform growth patterns, fertility, or metabolic function. In addition, 386 copy number variation regions (CNVRs) shared by all breeds were detected, covering 33.6 Mb (1.48% of the autosomal genome). Further, 132 exclusive CNVRs were identified for American breeds, 47 for Asian and Oceanian breeds, and 114 for European breeds. Functional enrichment analysis revealed genes within the common CNVRs involved in body height determination and other growth-related parameters. Exclusive CNVRs were located in the region of genes enriched for lipid metabolism in American minipigs, reproductive traits in Asian and Oceanian breeds, and cardiovascular features and body height in European breeds. In the selected groups, quantitative trait loci associated with body size, meat quality, reproduction, and disease susceptibility were highlighted.

CONCLUSION

This investigation of the CNV landscape of minipigs underlines the impact of selective breeding on structural variants and its role in the development of specific breed phenotypes across geographical areas. The multi-tool approach provides a valuable resource for future studies on the effects of artificial selection on livestock genomes.

Genomic consequences of intensive inbreeding in miniature inbred pigs

BACKGROUND

Inbreeding, a central theme in evolutionary and conservation biology, is a crucial practice in breeding to stabilize and enhance the specific traits or to establish inbred lines. It also carries the risk of inbreeding depression, reduced fitness, and increased potential for extinction. Nevertheless, inbreeding has been extensively studied in small and endangered populations but its effects in large domesticated animals are poorly understood. Here, we aim to investigate the genomic consequences of inbreeding in the Banna miniature inbred pig (BN), a breed that has been inbred for over 40 years.

RESULTS

We have sequenced 41 genomes of BN and Diannan miniature pig (DN) at high-coverage (> 31×) and combined them with published whole-genomes of swine to comprehensively investigate the genetic consequences of inbreeding. We find that BN is genetically closely related to DN, which is consistent with breeding records. All families of BN have undergone an extreme bottleneck due to intensive inbreeding, resulting in higher genomic inbreeding coefficients, reduced genetic diversity, and a lower effective population size (Ne) compare to non-inbred pigs. Furthermore, BN and DN exhibit an increased genetic load relative to Asian wild boars. Prolonged inbreeding and bottlenecks have led to some purging of deleterious mutations in BN compared to DN, and a conversion from masked load to realized load.

CONCLUSIONS

We present a comprehensive analysis to understand and assess the consequences of inbreeding in miniature inbred pigs from a perspective of population genomics. Utilizing genomic measurements proves effective in estimating the consequences of inbreeding, especially when a detailed and accurate historical record of pedigree are lacking. Our results provide valuable resources and a detailed perspective on the genomic impacts of inbreeding, potentially guiding efforts in breeding, breed improvement, and conservation.

Cyclic loading failed to promote growth in a pig model of midfacial hypoplasia

Yucatan miniature pigs, often used as large animal models in clinical research, are distinguished by a breed-specific midfacial hypoplasia with anterior crossbite. Although this deformity can be corrected by distraction osteogenesis, a less invasive method is desirable. We chose a mechanical cyclic stimulation protocol that has been successful in enhancing sutural growth in small animals and in a pilot study on standard pigs. Yucatan minipigs (n = 14) were obtained in pairs, with one of each pair randomly assigned to sham or loaded groups. All animals had loading implants installed on the right nasal and frontal bones and received labels for cell proliferation and mineral apposition. After a week of healing and under anesthesia, experimental animals received cyclic tensile loads (2.5 Hz, 30 min) delivered to the right nasofrontal suture daily for 5 days. Sutural strains were recorded at the final session for experimental animals. Sham animals received the same treatment except without loading or strain gauge placement. In contrast to pilot results on standard pigs, the treatment did not produce the expected sutural widening and increased growth. Although sutures were not fused and strains were in the normal range, the targeted right nasofrontal suture was narrowed rather than widened, with no statistically significant changes in sutural cell proliferation, mineral apposition, or vascularity. In general, Yucatan minipig sutures were more vascular than those of standard pigs and also tended to have more proliferating cells. In conclusion, either because the sutures themselves are abnormal or because of growth restrictions elsewhere in the skull, this cyclic loading protocol was unable to produce the desired response of sutural widening and growth. This treatment, effective in normal animals, did not improve naturally occurring midfacial hypoplasia in Yucatan minipigs.

Phylogenetic Relations and High-Altitude Adaptation in Wild Boar (Sus scrofa), Identified Using Genome-Wide Data

The Qinghai-Tibet Plateau (QTP) wild boar is an excellent model for investigating high-altitude adaptation. In this study, we analyzed genome-wide data from 93 wild boars compiled from various studies worldwide, including the QTP, southern and northern regions of China, Europe, Northeast Asia, and Southeast Asia, to explore their phylogenetic patterns and high-altitude adaptation based on genome-wide selection signal analysis and run of homozygosity (ROH) estimation. The findings demonstrate the alignment between the phylogenetic associations among wild boars and their geographical location. An ADMIXTURE analysis indicated a relatively close genetic relationship between QTP and southern Chinese wild boars. Analyses of the fixation index and cross-population extended haplotype homozygosity between populations revealed 295 candidate genes (CDGs) associated with high-altitude adaptation, such as TSC2, TELO2, SLC5A1, and SLC5A4. These CDGs were significantly overrepresented in pathways such as the mammalian target of rapamycin signaling and Fanconi anemia pathways. In addition, 39 ROH islands and numerous selective CDGs (e.g., SLC5A1, SLC5A4, and VCP), which are implicated in glucose metabolism and mitochondrial function, were discovered in QTP wild boars. This study not only assessed the phylogenetic history of QTP wild boars but also advanced our comprehension of the genetic mechanisms underlying the adaptation of wild boars to high altitudes.

Convergent dwarfism consequences of minipigs under independent artificial selections

BACKGROUND

Currently, diverse minipigs have acquired a common dwarfism phenotype through independent artificial selections. Characterizing the population and genetic diversity in minipigs is important to unveil genetic mechanisms regulating their body sizes and effects of independent artificial selections on those genetic mechanisms. However, full understanding for the genetic mechanisms and phenotypic consequences in minipigs still lag behind.

RESULTS

Here, using whole genome sequencing data of 41 pig breeds, including eight minipigs, we identified a large genomic diversity in a minipig population compared to other pig populations in terms of population structure, demographic signatures, and selective signatures. Selective signatures reveal diverse biological mechanisms related to body size in minipigs. We also found evidence for neural development mechanism as a minipig-specific body size regulator. Interestingly, selection signatures within those mechanisms containing neural development are also highly different among minipig breeds. Despite those large genetic variances, PLAG1, CHM, and ESR1 are candidate key genes regulating body size which experience different differentiation directions in different pig populations.

CONCLUSIONS

These findings present large variances of genetic structures, demographic signatures, and selective signatures in the minipig population. They also highlight how different artificial selections with large genomic diversity have shaped the convergent dwarfism.

Genome-wide detection of selection signatures in Jianli pigs reveals novel cis-regulatory haplotype in EDNRB associated with two-end black coat color

BACKGROUND

Jianli pig, a renowned indigenous breed in China, has the characteristics of a two-end black (TEB) coat color, excellent meat quality, strong adaptability and increased prolificacy. However, there is limited information available regarding the genetic diversity, population structure and genomic regions under selection of Jianli pig. On the other hand, the genetic mechanism of TEB coat color has remained largely unknown.

RESULTS

In this study, the whole genome resequencing of 30 Jianli pigs within a context of 153 individuals representing 13 diverse breeds was performed. The population structure analysis revealed that Jianli pigs have close genetic relationships with the Tongcheng pig breed, their geographical neighbors. Three methods (observed heterozygosity, expected heterozygosity, and runs of homozygosity) implied a relatively high level of genetic diversity and, a low inbreeding coefficient in Jianli compared with other pigs. We used Fst and XP-EHH to detect the selection signatures in Jianli pigs compared with Asian wild boar. A total of 451 candidate genes influencing meat quality (CREBBP, ADCY9, EEPD1 and HDAC9), reproduction (ESR1 and FANCA), and coat color (EDNRB, MITF and MC1R), were detected by gene annotation analysis. Finally, to fine-map the genomic region for the two-end black (TEB) coat color phenotype in Jianli pigs, we performed three signature selection methods between the TEB coat color and no-TEB coat color pig breeds. The current study, further confirmed that the EDNRB gene is a candidate gene for TEB color phenotype found in Chinese pigs, including Jinhua pigs, and the haplotype harboring 25 SNPs in the EDNRB gene may promote the formation of TEB coat color. Further ATAC-seq and luciferase reporter assays of these regions suggest that the 25-SNPs region was a strong candidate causative mutation that regulates the TEB coat color phenotype by altering enhancer function.

CONCLUSION

Our results advanced the understanding of the genetic mechanism behind artificial selection, and provided further resources for the protection and breeding improvement of Jianli pigs.

Tissue-engineered vocal fold replacement in swine: Methods for functional and structural analysis

We have developed a cell-based outer vocal fold replacement (COVR) as a potential therapy to improve voice quality after vocal fold (VF) injury, radiation, or tumor resection. The COVR consists of multipotent human adipose-derived stem cells (hASC) embedded within a three-dimensional fibrin scaffold that resembles vocal fold epithelium and lamina propria layers. Previous work has shown improved wound healing in rabbit studies. In this pilot study in pigs, we sought to develop methods for large animal implantation and phonatory assessment. Feasibility, safety, and structural and functional outcomes of the COVR implant are described. Of eight pigs studied, six animals underwent COVR implantation with harvest between 2 weeks and 6 months. Recovery of laryngeal tissue structure was assessed by vibratory and histologic analyses. Recovery of voice function was assessed by investigating acoustic parameters that were derived specifically for pigs. Results showed improved lamina propria qualities relative to an injured control animal at 6 months. Acoustic parameters reflected voice worsening immediately after surgery as expected; acoustics displayed clear voice recovery in the animal followed for 6 months after COVR. These methods form the basis for a larger-scale long-term pre-clinical safety and efficacy study.

The nasal septum and midfacial growth

The nasal septum is the only element of the chondrocranium which never completely ossifies. The persistence of this nonarticular cartilage has given rise to a variety of theories concerning cranial mechanics and growth of the midface. Previously, using pigs, we demonstrated that the septum is not a strut supporting the snout and that septal growth seems capable of stretching the overlying nasofrontal suture, a major contributor to snout elongation. Here we investigate whether abnormalities of the septum are implicated in cases of midfacial hypoplasia, in which growth of the midface is inadequate. Mild midfacial hypoplasia is common in domestic pig breeds and often severe in the Yucatan minipig, a popular laboratory breed. Normal-snouted and midfacial hypoplastic heads of standard (farm mixed breed) and minipigs ranging in age from perinatal to 12 months were dissected, imaged by CT, and/or prepared for histology. Even at birth, Yucatan minipigs with midfacial hypoplasia exhibited greater caudal ossification than normal; the ventral cartilaginous sphenoidal "tail" was diminished or missing. In addition, cells that morphologically appeared to have divided recently were less numerous than in newborn standard pigs. Juvenile Yucatan minipigs lacked caudal cartilaginous growth zones almost completely. In standard newborns, the ventral caudal septum was more replicative than the dorsal, but this trend was not seen in Yucatan newborns. In conclusion, accelerated maturation of the caudal septum was associated with midfacial hypoplasia, a further indication that the septum, particularly its ventral portion, is important for midfacial elongation.

Genome-Wide Association Study of Growth Traits in a Four-Way Crossbred Pig Population

Growth traits are crucial economic traits in the commercial pig industry and have a substantial impact on pig production. However, the genetic mechanism of growth traits is not very clear. In this study, we performed a genome-wide association study (GWAS) based on the specific-locus amplified fragment sequencing (SLAF-seq) to analyze ten growth traits on 223 four-way intercross pigs. A total of 227,921 highly consistent single nucleotide polymorphisms (SNPs) uniformly dispersed throughout the entire genome were used to conduct GWAS. A total of 53 SNPs were identified for ten growth traits using the mixed linear model (MLM), of which 18 SNPs were located in previously reported quantitative trait loci (QTL) regions. Two novel QTLs on SSC4 and SSC7 were related to average daily gain from 30 to 60 kg (ADG30-60) and body length (BL), respectively. Furthermore, 13 candidate genes (ATP5O, GHRHR, TRIM55, EIF2AK1, PLEKHA1, BRAP, COL11A2, HMGA1, NHLRC1, SGSM1, NFATC2, MAML1, and PSD3) were found to be associated with growth traits in pigs. The GWAS findings will enhance our comprehension of the genetic architecture of growth traits. We suggested that these detected SNPs and corresponding candidate genes might provide a biological foundation for improving the growth and production performance of pigs in swine breeding.

Integration of multi-omics data reveals cis-regulatory variants that are associated with phenotypic differentiation of eastern from western pigs

Background

The genetic mechanisms that underlie phenotypic differentiation in breeding animals have important implications in evolutionary biology and agriculture. However, the contribution of cis-regulatory variants to pig phenotypes is poorly understood. Therefore, our aim was to elucidate the molecular mechanisms by which non-coding variants cause phenotypic differences in pigs by combining evolutionary biology analyses and functional genomics.

Results

We obtained a high-resolution phased chromosome-scale reference genome with a contig N50 of 18.03 Mb for the Luchuan pig breed (a representative eastern breed) and profiled potential selective sweeps in eastern and western pigs by resequencing the genomes of 234 pigs. Multi-tissue transcriptome and chromatin accessibility analyses of these regions suggest that tissue-specific selection pressure is mediated by promoters and distal cis-regulatory elements. Promoter variants that are associated with increased expression of the lysozyme (LYZ) gene in the small intestine might enhance the immunity of the gastrointestinal tract and roughage tolerance in pigs. In skeletal muscle, an enhancer-modulating single-nucleotide polymorphism that is associated with up-regulation of the expression of the troponin C1, slow skeletal and cardiac type (TNNC1) gene might increase the proportion of slow muscle fibers and affect meat quality.

Conclusions

Our work sheds light on the molecular mechanisms by which non-coding variants shape phenotypic differences in pigs and provides valuable resources and novel perspectives to dissect the role of gene regulatory evolution in animal domestication and breeding.

Pathological Evaluation of Porcine Circovirus 2d (PCV2d) Strain and Comparative Evaluation of PCV2d and PCV2b Inactivated Vaccines against PCV2d Infection in a Specific Pathogen-Free (SPF) Yucatan Miniature Pig Model

Porcine circovirus type 2 (PCV2) is an economically important swine pathogen that causes porcine circovirus-associated diseases (PCVADs). The objective of this study was to evaluate the use of specific pathogen-free Yucatan miniature pigs (YMPs) as an experimental model for PCV2d challenge and vaccine assessment because PCV2-negative pigs are extremely rare in conventional swine herds in Korea. In the first experiment, every three pigs were subjected to PCV2d field isolate or mock challenge. During three weeks of experiments, the PCV2d infection group exhibited clinical outcomes of PCVAD with high viral loads, lymphoid depletion, and detection of PCV2d antigens in lymphoid organs by immunohistochemistry. In the second experiment, three groups of pigs were challenged with PCV2d after immunization for three weeks: a nonvaccinated group (three pigs), a PCV2b-Vac group vaccinated with a commercial PCV2b-based inactivated vaccine SuiShot^® Circo-ONE (five pigs), and a PCV2d-Vac group vaccinated with an experimental PCV2d-based inactivated vaccine (five pigs). During the three weeks of the challenge period, nonvaccinated pigs showed similar clinical outcomes to those observed in the PCV2d infection group from the first experiment. In contrast, both the PCV2b and PCV2d vaccinations produced good levels of protection against PCV2d challenge, as evidenced by reduced viral loads, improved growth performance, high virus-neutralizing antibody titers, and less development of PCV2-associated pathological lesions. Taken together, these data suggest that YMPs could be an alternative model for PCV2 challenge experiments, and these animals displayed typical clinical and pathological features and characteristics of protective immunity induced by the vaccines that were consistent with those resulting from PCV2 infections in conventional pigs.

Identifying selection signatures for litter size in Guangxi Bama Xiang pigs

Litter size is an important economic trait in pig production. However, the genetic mechanisms underlying varying litter size in Guangxi Bama Xiang pigs remain unknown. To identify selection signatures for litter size in Guangxi Bama Xiang pigs, we obtained 297 Illumina PorcineSNP50 BeadChip array data and the average born number (ABN) from parity one to nine in Guangxi Bama Xiang pigs. Fixation index (Fst) methods were used to identify the selection signature of the litter size, and three phenotypic gradient differential population pairs (according to the ABN) in individuals were used to reduce the false positives of signature selections. Single nucleotide polymorphisms (SNPs) were identified in the VEGFA promoter and exons. The general linear model was used to analyse the differences in distinct genotypes after they were typed using three-round multiplex PCR technology. Finally, the transcriptome factor and CpG island in the VEGFA promoter were predicted. A total of 328, 328 and 317 significant loci were identified in the 1st, 2nd and 3rd population pairs, respectively. After removing the false positives, 25 SNPs were defined as the selection signatures in relation to litter size. Ten (VEGFA, USP49, USP25, SRPK1, SLC26A8, RPL10A, PPARD, MAPK14, HMGA1 and CHRDL2) out of 52 genes in the selection regions were annotated as the candidate genes of litter size, respectively, VEGFA. There were no SNPs in the VEGFA exon region, but we obtained three SNPs (rs786889605, rs343769603 and rs323942424) in the VEGFA promoter regions. The ABN in CC was significantly higher than that in TT in rs786889605, and the ABN in TT was significantly lower than that in GG in rs323942424. Meanwhile, the mutation of the VEGFA promoter result in the loss of Sp1 and NF-1 and the formation of Oct-1. In summary, we obtained ten candidate genes, and two mutations in the VEGFA promoter that could be important potential molecular biomarkers for litter size in Bama Xiang pigs.

Identification of Body Size Determination Related Candidate Genes in Domestic Pig Using Genome-Wide Selection Signal Analysis

This study aimed to identify the genes related to the body size of pigs by conducting genome-wide selection analysis (GWSA). We performed a GWSA scan on 50 pigs belonging to four small-bodied pig populations (Diannan small-eared pig, Bama Xiang pig, Wuzhishan pig, and Jeju black pig from South Korea) and 124 large-bodied pigs. We used the genetic parameters of the pairwise fixation index (FST) and π ratio (case/control) to screen candidate genome regions and genes related to body size. The results revealed 47,339,509 high-quality SNPs obtained from 174 individuals, while 280 interacting candidate regions were obtained from the top 1% signal windows of both parameters, along with 187 genes (e.g., ADCK4, AMDHD2, ASPN, ASS1, and ATP6V0C). The results of the candidate gene (CG) annotation showed that a series of CGs (e.g., MSTN, LTBP4, PDPK1, PKMYT1, ASS1, and STAT6) was enriched into the gene ontology terms. Moreover, molecular pathways, such as the PI3K-Akt, HIF-1, and AMPK signaling pathways, were verified to be related to body development. Overall, we identified a series of key genes that may be closely related to the body size of pigs, further elucidating the heredity basis of body shape determination in pigs and providing a theoretical reference for molecular breeding.

Multi-omics approaches for comprehensive analysis and understanding of the immune response in the miniature pig breed

The porcine immune system has an important role in pre-clinical studies together with understanding the biological response mechanisms before entering into clinical trials. The size distribution of the Korean minipig is an important feature that make this breed ideal for biomedical research and safe practice in post clinical studies. The extremely tiny (ET) minipig serves as an excellent model for various biomedical research studies, but the comparatively frail and vulnerable immune response to the environment over its Large (L) size minipig breed leads to additional after born care. To overcome this pitfall, comparative analysis of the genomic regions under selection in the L type breed could provide a better understanding at the molecular level and lead to the development of an enhanced variety of ET type minipig. In this study, we utilized whole genome sequencing (WGS) to identify traces of artificial selection and integrated them with transcriptome data generated from blood samples to find strongly selected and differentially expressed genes of interest. We identified a total of 35 common genes among which 7 were differentially expressed and showed selective sweep in the L type over the ET type minipig breed. The stabilization of these genes were further confirmed using nucleotide diversity analysis, and these genes could serve as potential biomarkers for the development of a better variety of ET type pig breed.

Yucatan Minipig Knee Meniscus Regional Biomechanics and Biochemical Structure Support its Suitability as a Large Animal Model for Translational Research

Knee meniscus injuries are the most frequent causes of orthopedic surgical procedures in the U.S., motivating tissue engineering attempts and the need for suitable animal models. Despite extensive use in cardiovascular research and the existence of characterization data for the menisci of farm pigs, the farm pig may not be a desirable preclinical model for the meniscus due to rapid weight gain. Minipigs are conducive to in vivo experiments due to their slower growth rate than farm pigs and similarity in weight to humans. However, characterization of minipig knee menisci is lacking. The objective of this study was to extensively characterize structural and functional properties within different regions of both medial and lateral Yucatan minipig knee menisci to inform this model’s suitability as a preclinical model for meniscal therapies. Menisci measured 23.2–24.8 mm in anteroposterior length (33–40 mm for human), 7.7–11.4 mm in width (8.3–14.8 mm for human), and 6.4–8.4 mm in peripheral height (5–7 mm for human). Per wet weight, biochemical evaluation revealed 23.9–31.3% collagen (COL; 22% for human) and 1.20–2.57% glycosaminoglycans (GAG; 0.8% for human). Also, per dry weight, pyridinoline crosslinks (PYR) were 0.12–0.16% (0.12% for human) and, when normalized to collagen content, reached as high as 1.45–1.96 ng/µg. Biomechanical testing revealed circumferential Young’s modulus of 78.4–116.2 MPa (100–300 MPa for human), circumferential ultimate tensile strength (UTS) of 18.2–25.9 MPa (12–18 MPa for human), radial Young’s modulus of 2.5–10.9 MPa (10–30 MPa for human), radial UTS of 2.5–4.2 MPa (1–4 MPa for human), aggregate modulus of 157–287 kPa (100–150 kPa for human), and shear modulus of 91–147 kPa (120 kPa for human). Anisotropy indices ranged from 11.2–49.4 and 6.3–11.2 for tensile stiffness and strength (approximately 10 for human), respectively. Regional differences in mechanical and biochemical properties within the minipig medial meniscus were observed; specifically, GAG, PYR, PYR/COL, radial stiffness, and Young’s modulus anisotropy varied by region. The posterior region of the medial meniscus exhibited the lowest radial stiffness, which is also seen in humans and corresponds to the most prevalent location for meniscal lesions. Overall, similarities between minipig and human menisci support the use of minipigs for meniscus translational research.

Assessment of Heterozygosity and Genome-Wide Analysis of Heterozygosity Regions in Two Duroc Pig Populations

Heterozygosity can effectively reflect the diverse models of population structure and demographic history. However, the genomic distribution of heterozygotes and the correlation between regions of heterozygosity (runs of heterozygosity, ROHet) and phenotypes are largely understudied in livestock. The objective of this study was to identify ROHet in the Duroc pig genome, and investigate the relationships between ROHet and eight important economic traits. Here, we genotyped 3,770 American Duroc (S21) and 2,096 Canadian Duroc (S22) pigs using 50 K single nucleotide polymorphism array to analyze heterozygosity. A total of 145,010 and 84,396 ROHets were characterized for S21 and S22 populations, respectively. ROHet segments were mostly enriched in 1–2 Mb length classification (75.48% in S21 and 72.25% in S22). The average genome length covered by ROHet was 66.53 ± 12.20 Mb in S21 and 73.32 ± 13.77 Mb in S22 pigs. Additionally, we detected 20 and 13 ROHet islands in S21 and S22 pigs. Genes in these genomic regions were mainly involved in the biological processes of immunity and reproduction. Finally, the genome-wide ROHet-phenotypes association analysis revealed that 130 ROHets of S21 and 84 ROHets of S22 were significantly associated with eight economic traits. Among the candidate genes in the significant ROHet regions, 16 genes related to growth, metabolism, and meat quality were considered as candidate genes for important economic traits of pigs. This work preliminarily explores the effect of heterozygosity-rich regions in the pig genome on production performance and provides new insights for subsequent research on pig genetic improvement.

Detecting Selection in Multiple Populations by Modeling Ancestral Admixture Components

One of the most powerful and commonly used approaches for detecting local adaptation in the genome is the identification of extreme allele frequency differences between populations. In this article, we present a new maximum likelihood method for finding regions under positive selection. It is based on a Gaussian approximation to allele frequency changes and it incorporates admixture between populations. The method can analyze multiple populations simultaneously and retains power to detect selection signatures specific to ancestry components that are not representative of any extant populations. Using simulated data, we compare our method to related approaches, and show that it is orders of magnitude faster than the state-of-the-art, while retaining similar or higher power for most simulation scenarios. We also apply it to human genomic data and identify loci with extreme genetic differentiation between major geographic groups. Many of the genes identified are previously known selected loci relating to hair pigmentation and morphology, skin, and eye pigmentation. We also identify new candidate regions, including various selected loci in the Native American component of admixed Mexican-Americans. These involve diverse biological functions, such as immunity, fat distribution, food intake, vision, and hair development.

Comparative Proteomics and Phosphoproteomics Analysis Reveal the Possible Breed Difference in Yorkshire and Duroc Boar Spermatozoa

Sperm cells are of unique elongated structure and function, the development of which is tightly regulated by the existing proteins and the posttranslational modifications (PTM) of these proteins. Based on the phylogenetic relationships of various swine breeds, Yorkshire boar is believed to be distinctly different from Duroc boar. The comprehensive differential proteomics and phosphoproteomics profilings were performed on spermatozoa from both Yorkshire and Duroc boars. By both peptide and PTM peptide quantification followed by statistical analyses, 167 differentially expressed proteins were identified from 1,745 proteins, and 283 differentially expressed phosphopeptides corresponding to 102 unique differentially phosphorylated proteins were measured from 1,140 identified phosphopeptides derived from 363 phosphorylated proteins. The representative results were validated by Western blots. Pathway enrichment analyses revealed that majority of differential expression proteins and differential phosphorylation proteins were primarily concerned with spermatogenesis, male gamete generation, sperm motility, energy metabolism, cilium morphogenesis, axonemal dynein complex assembly, sperm–egg recognition, and capacitation. Remarkably, axonemal dynein complex assembly related proteins, such as SMCP, SUN5, ODF1, AKAP3, and AKAP4 that play a key regulatory role in the sperm physiological functions, were significantly higher in Duroc spermatozoa than that of Yorkshire. Furthermore, phosphorylation of sperm-specific proteins, such as CABYR, ROPN1, CALM1, PRKAR2A, and PRKAR1A, participates in regulation of the boar sperm motility mainly through the cAMP/PKA signal pathway in different breeds, demonstrating that protein phosphorylation may be an important mechanism underlying the sperm diversity. Protein–protein interaction analysis revealed that the 14 overlapped proteins between differential expression proteins and differential phosphorylation proteins potentially played a key role in sperm development and motility of the flagellum, including the proteins ODF1, SMCP, AKAP4, FSIP2, and SUN5. Taken together, these physiologically and functionally differentially expressed proteins (DEPs) and differentially expressed phosphorylated proteins (DPPs) may constitute the proteomic backgrounds between the two different boar breeds. The validation will be performed to delineate the roles of these PTM proteins as modulators of Yorkshire and Duroc boar spermatozoa.

Variability in porcine microRNA genes and its association with mRNA expression and lipid phenotypes

BACKGROUND

Mature microRNAs (miRNAs) play an important role in repressing the expression of a wide range of mRNAs. The presence of polymorphic sites in miRNA genes and their corresponding 3'UTR binding sites can disrupt canonical conserved miRNA-mRNA pairings, and thus modify gene expression patterns. However, to date such polymorphic sites in miRNA genes and their association with gene expression phenotypes and complex traits are poorly characterized in pigs.

RESULTS

By analyzing whole-genome sequences from 120 pigs and wild boars from Europe and Asia, we identified 285 single nucleotide polymorphisms (SNPs) that map to miRNA loci, and 109,724 SNPs that are located in predicted 7mer-m8 miRNA binding sites within porcine 3'UTR. In porcine miRNA genes, SNP density is reduced compared with their flanking non-miRNA regions. By sequencing the genomes of five Duroc boars, we identified 12 miRNA SNPs that were subsequently genotyped in their offspring (N = 345, Lipgen population). Association analyses of miRNA SNPs with 38 lipid-related traits and hepatic and muscle microarray expression phenotypes recorded in the Lipgen population were performed. The most relevant detected association was between the genotype of the rs319154814 (G/A) SNP located in the apical loop of the ssc-miR-326 hairpin precursor and PPP1CC mRNA levels in the liver (q-value = 0.058). This result was subsequently confirmed by qPCR (P-value = 0.027). The rs319154814 (G/A) genotype was also associated with several fatty acid composition traits.

CONCLUSIONS

Our findings show a reduced variability of porcine miRNA genes, which is consistent with strong purifying selection, particularly in the seed region that plays a critical role in miRNA binding. Although it is generally assumed that SNPs mapping to the seed region are those with the most pronounced consequences on mRNA expression, we show that a SNP mapping to the apical region of ssc-miR-326 is significantly associated with hepatic mRNA levels of the PPP1CC gene, one of its predicted targets. Although experimental confirmation of such an interaction is reported in humans but not in pigs, this result highlights the need to further investigate the functional effects of miRNA polymorphisms that are located outside the seed region on gene expression in pigs.

Porcine models of acute kidney injury

Pigs represent a potentially attractive model for medical research. Similar body size and physiological patterns of kidney injury that more closely mimic those described in humans make larger animals attractive for experimentation. Using larger animals, including pigs, to investigate the pathogenesis of acute kidney injury (AKI) also serves as an experimental bridge, narrowing the gap between clinical disease and preclinical discoveries. This article compares the advantages and disadvantages of large versus small AKI animal models and provides a comprehensive overview of the development and application of porcine models of AKI induced by clinically relevant insults, including ischemia-reperfusion, sepsis, and nephrotoxin exposure. The primary focus of this review is to evaluate the use of pigs for AKI studies by current investigators, including areas where more information is needed.

Integration of multi-omics approaches for functional characterization of muscle related selective sweep genes in Nanchukmacdon

Pig as a food source serves daily dietary demand to a wide population around the world. Preference of meat depends on various factors with muscle play the central role. In this regards, selective breeding abled us to develop "Nanchukmacdon" a pig breeds with an enhanced variety of meat and high fertility rate. To identify genomic regions under selection we performed whole-genome resequencing, transcriptome, and whole-genome bisulfite sequencing from Nanchukmacdon muscles samples and used published data for three other breeds such as Landrace, Duroc, Jeju native pig and analyzed the functional characterization of candidate genes. In this study, we present a comprehensive approach to identify candidate genes by using multi-omics approaches. We performed two different methods XP-EHH, XP-CLR to identify traces of artificial selection for traits of economic importance. Moreover, RNAseq analysis was done to identify differentially expressed genes in the crossed breed population. Several genes (UGT8, ZGRF1, NDUFA10, EBF3, ELN, UBE2L6, NCALD, MELK, SERP2, GDPD5, and FHL2) were identified as selective sweep and differentially expressed in muscles related pathways. Furthermore, nucleotide diversity analysis revealed low genetic diversity in Nanchukmacdon for identified genes in comparison to related breeds and whole-genome bisulfite sequencing data shows the critical role of DNA methylation pattern in identified genes that leads to enhanced variety of meat. This work demonstrates a way to identify the molecular signature and lays a foundation for future genomic enabled pig breeding.

Selective sweep analysis reveals extensive parallel selection traits between large white and Duroc pigs

In the process of pig genetic improvement, different commercial breeds have been bred for the same purpose, improving meat production. Most of the economic traits, such as growth and fertility, have been selected similarly despite the discrepant selection pressure, which is known as parallel selection. Here, 28 whole-genome sequencing data of Danish large white pigs with an approximately 25-fold depth each were generated, resulting in about 12 million high-quality SNPs for each individual. Combined with the sequencing data of 27 Duroc and 23 European wild boars, we investigated the parallel selection of Danish large white and Duroc pigs using two complementary methods, Fst and iHS. In total, 67 candidate regions were identified as the signatures of parallel selection. The genes in candidate regions of parallel selection were mainly associated with sensory perception, growth rate, and body size. Further functional annotation suggested that the striking consistency of the terms may be caused by the polygenetic basis of quantitative traits, and revealing the complex genetic basis of parallel selection. Besides, some unique terms were enriched in population-specific selection regions, such as the limb development-related terms enriched in Duroc-specific selection regions, suggesting unique selections of breed-specific selected traits. These results will help us better understand the parallel selection process of different breeds. Moreover, we identified several potential causal SNPs that may contribute to the pig genetic breeding process.

A porcine model of phenylketonuria generated by CRISPR/Cas9 genome editing

Phenylalanine hydroxylase-deficient (PAH-deficient) phenylketonuria (PKU) results in systemic hyperphenylalaninemia, leading to neurotoxicity with severe developmental disabilities. Dietary phenylalanine (Phe) restriction prevents the most deleterious effects of hyperphenylalaninemia, but adherence to diet is poor in adult and adolescent patients, resulting in characteristic neurobehavioral phenotypes. Thus, an urgent need exists for new treatments. Additionally, rodent models of PKU do not adequately reflect neurocognitive phenotypes, and thus there is a need for improved animal models. To this end, we have developed PAH-null pigs. After selection of optimal CRISPR/Cas9 genome-editing reagents by using an in vitro cell model, zygote injection of 2 sgRNAs and Cas9 mRNA demonstrated deletions in preimplantation embryos, with embryo transfer to a surrogate leading to 2 founder animals. One pig was heterozygous for a PAH exon 6 deletion allele, while the other was compound heterozygous for deletions of exon 6 and of exons 6-7. The affected pig exhibited hyperphenylalaninemia (2000-5000 μM) that was treatable by dietary Phe restriction, consistent with classical PKU, along with juvenile growth retardation, hypopigmentation, ventriculomegaly, and decreased brain gray matter volume. In conclusion, we have established a large-animal preclinical model of PKU to investigate pathophysiology and to assess new therapeutic interventions.

Evolutionary and Functional Analysis of Korean Native Pig Using Single Nucleotide Polymorphisms

Time and cost-effective production of next-generation sequencing data has enabled the performance of population-scale comparative and evolutionary studies for various species, which are essential for obtaining the comprehensive insight into molecular mechanisms underlying species- or breed-specific traits. In this study, the evolutionary and functional analysis of Korean native pig (KNP) was performed using single nucleotide polymorphism (SNP) data by comparative and population genomic approaches with six different mammalian species and five pig breeds. We examined the evolutionary history of KNP SNPs, and the specific genes of KNP based on the uniqueness of non-synonymous SNPs among the used species and pig breeds. We discovered the evolutionary trajectory of KNP SNPs within the used mammalian species as well as pig breeds. We also found olfaction-associated functions that have been characterized and diversified during evolution, and quantitative trait loci associated with the unique traits of KNP. Our study provides new insight into the evolution of KNP and serves as a good example for a better understanding of domestic animals in terms of evolution and domestication using the combined approaches of comparative and population genomics.

No excessive mutations in transcription activator-like effector nuclease-mediated α-1,3-galactosyltransferase knockout Yucatan miniature pigs

OBJECTIVE

Specific genomic sites can be recognized and permanently modified by genome editing. The discovery of endonucleases has advanced genome editing in pigs, attenuating xenograft rejection and cross-species disease transmission. However, off-target mutagenesis caused by these nucleases is a major barrier to putative clinical applications. Furthermore, off-target mutagenesis by genome editing has not yet been addressed in pigs.

METHODS

Here, we generated genetically inheritable α-1,3-galactosyltransferase (GGTA1) knockout Yucatan miniature pigs by combining transcription activator-like effector nuclease (TALEN) and nuclear transfer. For precise estimation of genomic mutations induced by TALEN in GGTA1 knockout pigs, we obtained the whole-genome sequence of the donor cells for use as an internal control genome.

RESULTS

In-depth whole-genome sequencing analysis demonstrated that TALEN-mediated GGTA1 knockout pigs had a comparable mutation rate to homologous recombination-treated pigs and wild-type strain controls. RNA sequencing analysis associated with genomic mutations revealed that TALEN-induced off-target mutations had no discernable effect on RNA transcript abundance.

CONCLUSION

Therefore, TALEN appears to be a precise and safe tool for generating genome-edited pigs, and the TALEN-mediated GGTA1 knockout Yucatan miniature pigs produced in this study can serve as a safe and effective organ and tissue resource for clinical applications.

In silico approaches to discover the functional impact of non-synonymous single nucleotide polymorphisms in selective sweep regions of the Landrace genome

OBJECTIVE

The aim of this study was to discover the functional impact of non-synonymous single nucleotide polymorphisms (nsSNPs) that were found in selective sweep regions of the Landrace genome.

METHODS

Whole-genome re-sequencing data were obtained from 40 pigs, including 14 Landrace, 16 Yorkshire, and 10 wild boars, which were generated with the Illumina HiSeq 2000 platform. The nsSNPs in the selective sweep regions of the Landrace genome were identified, and the impacts of these variations on protein function were predicted to reveal their potential association with traits of the Landrace breed, such as reproductive capacity.

RESULTS

Total of 53,998 nsSNPs in the mapped regions of pigs were identified, and among them, 345 nsSNPs were found in the selective sweep regions of the Landrace genome which were reported previously. The genes featuring these nsSNPs fell into various functional categories, such as reproductive capacity or growth and development during the perinatal period. The impacts of amino acid sequence changes by nsSNPs on protein function were predicted using two in silico SNP prediction algorithms, i.e., sorting intolerant from tolerant and polymorphism phenotyping v2, to reveal their potential roles in biological processes that might be associated with the reproductive capacity of the Landrace breed.

CONCLUSION

The findings elucidated the domestication history of the Landrace breed and illustrated how Landrace domestication led to patterns of genetic variation related to superior reproductive capacity. Our novel findings will help understand the process of Landrace domestication at the genome level and provide SNPs that are informative for breeding.

Application of partial least squares in exploring the genome selection signatures between populations

Natural and artificial selection have led to substantial variation in the phenotypic traits of different populations. Therefore, there is a need to develop methods that are based on cross-population comparisons to discover loci related to specific traits. Here, we suggested a strategy to detect the genome selection signatures between populations based on the partial least squares (PLS) theory. Using the binary population indicator as the response variable in the PLS analysis, alleles under selection between populations were identified from the first PLS component. We explored the theory behind the PLS analysis to reveal its usefulness in detecting the loci under selection. Through the simulation study, the results showed that the PLS method had a better performance than the F_ST and EigenGWAS methods. In addition, by using the real data hapmap3, we found that rs11150606 in PRSS53 gene and rs1800414 in OCA2 gene were under selection between East Asian populations and three other populations, including African, American, and European populations. We concluded that this strategy was easily carried out and might supplement for the deficiency of the EigenGWAS method in some cases. To facilitate the application of this method, we developed an R script that is freely accessible at http://klab.sjtu.edu.cn/PLS/ .

Genome analysis of Yucatan miniature pigs to assess their potential as biomedical model animals

Objective

Pigs share many physiological, anatomical and genomic similarities with humans, which make them suitable models for biomedical researches. Understanding the genetic status of Yucatan miniature pigs (YMPs) and their association with human diseases will help to assess their potential as biomedical model animals. This study was performed to identify non-synonymous single nucleotide polymorphisms (nsSNPs) in selective sweep regions of the genome of YMPs and present the genetic nsSNP distributions that are potentially associated with disease occurrence in humans.

Methods

nsSNPs in whole genome resequencing data from 12 YMPs were identified and annotated to predict their possible effects on protein function. Sorting intolerant from tolerant (SIFT) and polymorphism phenotyping v2 analyses were used, and gene ontology (GO) network and Kyoto encyclopedia of genes and genomes (KEGG) pathway analyses were performed.

Results

The results showed that 8,462 genes, encompassing 72,067 nsSNPs were identified, and 118 nsSNPs in 46 genes were predicted as deleterious. GO network analysis classified 13 genes into 5 GO terms (p<0.05) that were associated with kidney development and metabolic processes. Seven genes encompassing nsSNPs were classified into the term associated with Alzheimer’s disease by referencing the genetic association database. The KEGG pathway analysis identified only one significantly enriched pathway (p<0.05), hsa04080: Neuroactive ligand-receptor interaction, among the transcripts.

Conclusion

The number of deleterious nsSNPs in YMPs was identified and then these variants-containing genes in YMPs data were adopted as the putative human diseases-related genes. The results revealed that many genes encompassing nsSNPs in YMPs were related to the various human genes which are potentially associated with kidney development and metabolic processes as well as human disease occurrence.

Genomic signatures reveal selection of characteristics within and between Meishan pig populations

The Chinese Meishan pig breed is well known for its high prolificacy. Moreover, this breed can be divided into three types based on their body size: big Meishan, middle Meishan (MMS) and small Meishan (SMS) pigs. Few studies have reported on the genetic signatures of Meishan pigs, particularly on a genome-wide scale. Exploring for genetic signatures could be quite valuable for revealing the genetic architecture of phenotypic variation. Thus, we performed research in two parts based on the genome reducing and sequencing data of 143 Meishan pigs (74 MMS pigs, 69 SMS pigs). First, we detected the selection signatures among all Meishan pigs studied using the relative extended haplotype homozygosity test. Second, we detected the selection signatures between MMS and SMS pigs using the cross-population extended haplotype homozygosity and F_ST methods. A total of 111 398 SNPs were identified from the sequenced genomes. In the population analysis, the most significant genes were associated with the mental development (RGMA), reproduction (HDAC4, FOXL2) and lipid metabolism (ACACB). From the cross-population analysis, we detected genes related to body weight (SPDEF, PACSIN1) in both methods. We suggest that rs341373351, located within the PACSIN1 gene, might be the causal variant. This study may have achieved consistency between selection signatures and characteristics within and between Meishan pig populations. These findings can provide insight into investigating the molecular background of high prolificacy and body size in pig.

Targeting regenerative exosomes to myocardial infarction using cardiac homing peptide

Rationale: Cardiac stem cell-derived exosomes have been demonstrated to promote cardiac regeneration following myocardial infarction in preclinical studies. Recent studies have used intramyocardial injection in order to concentrate exosomes in the infarct. Though effective in a research setting, this method is not clinically appealing due to its invasive nature. We propose the use of a targeting peptide, cardiac homing peptide (CHP), to target intravenously-infused exosomes to the infarcted heart. Methods: Exosomes were conjugated with CHP through a DOPE-NHS linker. Ex vivo targeting was analyzed by incubating organ sections with the CHP exosomes and analyzing with fluorescence microscopy. In vitro assays were performed on neonatal rat cardiomyocytes and H9C2 cells. For the animal study, we utilized an ischemia/reperfusion rat model. Animals were treated with either saline, scramble peptide exosomes, or CHP exosomes 24 h after surgery. Echocardiography was performed 4 h after surgery and 21 d after surgery. At 21 d, animals were sacrificed, and organs were collected for analysis. Results: By conjugating the exosomes with CHP, we demonstrate increased retention of the exosomes within heart sections ex vivo and in vitro with neonatal rat cardiomyocytes. In vitro studies showed improved viability, reduced apoptosis and increased exosome uptake when using CHP-XOs. Using an animal model of ischemia/reperfusion injury, we measured the heart function, infarct size, cellular proliferation, and angiogenesis, with improved outcomes with the CHP exosomes. Conclusions: Our results demonstrate a novel method for increasing delivery of for treatment of myocardial infarction. By targeting exosomes to the infarcted heart, there was a significant improvement in outcomes with reduced fibrosis and scar size, and increased cellular proliferation and angiogenesis.

PigVar: a database of pig variations and positive selection signatures

Pigs are excellent large-animal models for medical research and a promising organ donor source for transplant patients. Next-generation sequencing technology has yielded a dramatic increase in the volume of genomic data for pigs. However, the limited amount of variation data provided by dbSNP, and non-congruent criteria used for calling variation, present considerable hindrances to the utility of this data. We used a uniform pipeline, based on GATK, to identify non-redundant, high-quality, whole-genome SNPs from 280 pigs and 6 outgroup species. A total of 64.6 million SNPs were identified in 280 pigs and 36.8 million in the outgroups. We then used LUMPY to identify a total of 7 236 813 structural variations (SVs) in 211 pigs. Positively selected loci were identified through five statistical tests of different evolutionary attributes of the SNPs. Combining the non-redundant variations and the evolutionary selective scores, we built the first pig-specific variation database, PigVar (http://www.ibiomedical.net/pigvar/), which is a web-based open-access resource. PigVar collects parameters of the variations including summary lists of the locations of the variations within protein-coding and long intergenic non-coding RNA (lincRNA) genes, whether the SNPs are synonymous or non-synonymous, their ancestral and derived states, geographic sampling locations, as well as breed information. The PigVar database will be kept operational and updated to facilitate medical research using the pig as model and agricultural research including pig breeding.

Database URL: http://www.ibiomedical.net/pigvar/

Convergent and divergent genetic changes in the genome of Chinese and European pigs

Since 10,000 BC, continuous human selection has led to intense genetic and phenotypic changes in pig (Sus scrofa) domestication. Through whole genome analysis of 257 individuals, we demonstrated artificial unidirectional and bidirectional selection as the primary force to shape the convergent and divergent changes between Chinese domestic pigs (CHD) and European domestic pigs (EUD). We identified 31 genes in unidirectional selection regions that might be related to fundamental domestication requirements in pigs. And these genes belong predominantly to categories related to the nervous system, muscle development, and especially to metabolic diseases. In addition, 35 genes, representing different breeding preference, were found under bidirectional selection for the distinct leanness and reproduction traits between CHD and EUD. The convergent genetic changes, contributing physical and morphological adaption, represent the common concerns on pig domestication. And the divergent genetic changes reflect distinct breeding goals between Chinese and European pigs. Using ITPR3, AHR and NMU as examples, we explored and validated how the genetic variations contribute to the phenotype changes.

The Yucatan Minipig Temporomandibular Joint Disc Structure-Function Relationships Support Its Suitability for Human Comparative Studies

Frequent involvement of the disc in temporomandibular joint (TMJ) disorders warrants attempts to tissue engineer TMJ disc replacements. Physiologically, a great degree of similarity is seen between humans and farm pigs (FPs), but the pig's rapid growth confers a significant challenge for in vivo experiments. Minipigs have a slower growth rate and are smaller than FPs, but minipig TMJ discs have yet to be fully characterized. The objective of this study was to determine the suitability of the minipig for TMJ studies by extensive structural and functional characterization. The properties of minipig TMJ discs closely reproduced previously reported morphological, biochemical, and biomechanical values of human and FP discs. The width/length dimension ratio of the minipig TMJ disc was 1.95 (1.69 for human and 1.94 for FP). The biochemical evaluation revealed, on average per wet weight, 24.3% collagen (22.8% for human and 24.9% for FP); 0.8% glycosaminoglycan (GAG; 0.5% for human and 0.4% for FP); and 0.03% DNA (0.008% for human and 0.02% for FP). Biomechanical testing revealed, on average, compressive relaxation modulus of 50 kPa (37 kPa for human and 32 kPa for FP), compressive instantaneous modulus of 1121 kPa (1315 kPa for human and 1134 kPa for FP), and coefficient of viscosity of 13 MPa·s (9 MPa·s for human and 3 MPa·s for FP) at 20% strain. These properties also varied topographically in accordance to those of human and FP TMJ discs. Anisotropy, quantified by bidirectional tensile testing and histology, again was analogous among minipig, human, and FP TMJ discs. The minipig TMJ's ginglymoarthrodial nature was verified through cone beam computer tomography. Collectively, the similarities between minipig and human TMJ discs support the use of minipig as a relevant model for TMJ research; considering the practical advantages conferred by its growth rate and size, the minipig may be a preferred model over FP.

MicroRNA in sperm from Duroc, Landrace and Yorkshire boars

Sperm contain microRNAs (miRNAs), which may have roles in epigenetic control. Regarding phylogenetic relationships among various swine breeds, Yorkshire and Landrace, are considered phenotypically and genetically very similar, but distinctly different from Duroc. The objective of the present study was to compare abundance of boar sperm miRNAs in these three breeds. Overall, 252 prioritized miRNAs were investigated using real-time PCR; relative expression of miRNAs in sperm was similar in Yorkshire and Landrace boars, but significantly different compared to Duroc. Seventeen miRNAs (hsa-miR-196a-5p, hsa-miR-514a-3p, hsa-miR-938, hsa-miR-372-3p, hsa-miR-558, hsa-miR-579-3p, hsa-miR-595, hsa-miR-648, hsa-miR-524-3p, hsa-miR-512-3p, hsa-miR-429, hsa-miR-639, hsa-miR-551a, hsa-miR-624-5p, hsa-miR-585-3p, hsa-miR-508-3p and hsa-miR-626) were down-regulated (P < 0.05; fold regulation ≤-2) in Yorkshire and Landrace sperm, compared to Duroc sperm. Furthermore, three miRNAs (hsa-miR-9-5p, hsa-miR-150-5p, and hsa-miR-99a-5p) were significantly up-regulated in Yorkshire and Landrace sperm compared to Duroc sperm, However, 240 miRNAs were not significantly different (within + 2 fold) between Yorkshire and Landrace sperm. We concluded that miRNAs in sperm were not significantly different between Yorkshire and Landrace boars, but there were significant differences between those two breeds and Duroc boars. Furthermore, integrated target genes for selected down-regulated miRNAs (identified via an in-silico method) appeared to participate in spermatogenesis and sperm functions.

A decade of pig genome sequencing: a window on pig domestication and evolution

Insight into how genomes change and adapt due to selection addresses key questions in evolutionary biology and in domestication of animals and plants by humans. In that regard, the pig and its close relatives found in Africa and Eurasia represent an excellent group of species that enables studies of the effect of both natural and human-mediated selection on the genome. The recent completion of the draft genome sequence of a domestic pig and the development of next-generation sequencing technology during the past decade have created unprecedented possibilities to address these questions in great detail. In this paper, I review recent whole-genome sequencing studies in the pig and closely-related species that provide insight into the demography, admixture and selection of these species and, in particular, how domestication and subsequent selection of Sus scrofa have shaped the genomes of these animals.

Miniature Swine Breeds in Toxicology and Drug Safety Assessments

The use of miniature swine as a nonrodent species in safety assessment has continued to expand for over a decade, and they are becoming routinely used in toxicology and in pharmacology as well as a model for human diseases. Miniature swine models are regularly used for regulatory toxicity studies designed to assess safety of new therapeutic compounds given through different routes of exposure and are used as an alternative model to the canine or the nonhuman primate. Translational preclinical swine study data presented support the current finding that miniature swine are the animal model of choice for assessment of drug absorption, tolerance, and systemic toxicity following systemic exposures. Because research investigators need to be familiar with important anatomic and histopathologic features of the miniature swine in order to place toxicopathologic findings in their proper perspective, clinical and anatomic pathology data from a large number of Sinclair, Hanford, Yucatan, and Göttingen breeds from control groups from a wide variety of studies performed between 2004 and 2014 will be presented, compared, and partially illustrated.