Development and Genome Sequencing of a Laboratory-Inbred Miniature Pig Facilitates Study of Human Diabetic Disease

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.

A chromosome-level genome assembly of the Korean minipig (Sus scrofa)

Recent advancements in sequencing and genome assembly technologies have led to rapid generation of high-quality genome assemblies for various species and breeds. Despite the importance as minipigs an animal model in biomedical research, the construction of high-quality genome assemblies of minipigs still lags behind other pig breeds. To address this problem, we constructed a high-quality chromosome-level genome assembly of the Korean minipig (KMP) utilizing multiple different types of sequencing reads and reference genomes. The KMP assembly included 19 chromosome-level sequences with a total length of 2.52 Gb and N50 of 137 Mb. Comparative analyses with the pig reference genome (Sscrofa11.1) demonstrated comparable contiguity and completeness of the KMP assembly. Additionally, genome annotation analyses identified 22,666 protein-coding genes and repetitive elements occupying 40.10% of the genome. The KMP assembly and genome annotation provide valuable resources that can contribute to various future research on minipig and other pig breeds.

Microbiome-gut-brain axis contributes to patients and Bama miniature pigs with acute large ischemic stroke

Acute large hemispheric infarction (ALHI) is an overwhelming emergency with a great challenge of gastrointestinal dysfunction clinically. Here, we initially proposed delayed bowel movements as the clinical phenotype of strike to gut-brain axis (GBA) in ALHI patients by epidemiological analysis of 499 acute ischemic stroke (AIS) patients. 1H NMR-based metabolomics revealed that AIS markedly altered plasma global metabolic profiling of patients compared with healthy controls. Risk factors of strike on GBA were the National Institutes of Health Stroke Scale (NIHSS) score ≥ 5 and stroke onset time ≤ 24 h. As a result, first defecating time after admission to the hospital ≥2 days could be considered as a potential risk factor for strike on GBA. Subsequently, the ALHI Bama miniature (BM) pig model with acute symptomatic seizure was successfully established by ligation of the left ascending pharyngeal artery combined with local air injection. Clinical phenotypes of brain necrosis such as hemiplegia were examined with brain diffusion-weighted imaging (DWI) and pathological diagnosis. In addition to global brain injury and inflammation, we also found that ALHI induced marked alterations of intestinal barrier integrity, the gut microbial community, and microbiota-derived metabolites including serotonin and neurotransmitters in both plasma and multiple brain tissues of BM pigs. These findings revealed that microbiota-gut-brain axis highly contributed to the occurrence and development of ALHI.

Swine models in translational research and medicine

Swine are increasingly studied as animal models of human disease. The anatomy, size, longevity, physiology, immune system, and metabolism of swine are more like humans than traditional rodent models. In addition, the size of swine is preferred for surgical placement and testing of medical devices destined for humans. These features make swine useful for biomedical, pharmacological, and toxicological research. With recent advances in gene-editing technologies, genetic modifications can readily and efficiently be made in swine to study genetic disorders. In addition, gene-edited swine tissues are necessary for studies testing and validating xenotransplantation into humans to meet the critical shortfall of viable organs versus need. Underlying all of these biomedical applications, the knowledge of husbandry, background diseases and lesions, and biosecurity needs are important for productive, efficient, and reproducible research when using swine as a human disease model for basic research, preclinical testing, and translational studies.

Post-meiotic mechanism of facultative parthenogenesis in gonochoristic whiptail lizard species

Facultative parthenogenesis (FP) has historically been regarded as rare in vertebrates, but in recent years incidences have been reported in a growing list of fish, reptile, and bird species. Despite the increasing interest in the phenomenon, the underlying mechanism and evolutionary implications have remained unclear. A common finding across many incidences of FP is either a high degree of homozygosity at microsatellite loci or low levels of heterozygosity detected in next-generation sequencing data. This has led to the proposal that second polar body fusion following the meiotic divisions restores diploidy and thereby mimics fertilization. Here, we show that FP occurring in the gonochoristic Aspidoscelis species A. marmoratus and A. arizonae results in genome-wide homozygosity, an observation inconsistent with polar body fusion as the underlying mechanism of restoration. Instead, a high-quality reference genome for A. marmoratus and analysis of whole-genome sequencing from multiple FP and control animals reveals that a post-meiotic mechanism gives rise to homozygous animals from haploid, unfertilized oocytes. Contrary to the widely held belief that females need to be isolated from males to undergo FP, females housed with conspecific and heterospecific males produced unfertilized eggs that underwent spontaneous development. In addition, offspring arising from both fertilized eggs and parthenogenetic development were observed to arise from a single clutch. Strikingly, our data support a mechanism for facultative parthenogenesis that removes all heterozygosity in a single generation. Complete homozygosity exposes the genetic load and explains the high rate of congenital malformations and embryonic mortality associated with FP in many species. Conversely, for animals that develop normally, FP could potentially exert strong purifying selection as all lethal recessive alleles are purged in a single generation.

Multiomic study of the protective mechanism of Persicaria capitata (Buch.-Ham. ex D.Don) H.Gross against streptozotocin-induced diabetic nephropathy in Guizhou miniature pigs

BACKGROUND

Persicaria capitata (Buch.-Ham. ex D.Don) H.Gross (P. capitata, PCB), a traditional drug of the Miao people in China, is potential traditional drug used for the treatment of diabetic nephropathy (DN).

PURPOSE

The purpose of this study is to investigate the function of P. capitata and clarify its protective mechanism against DN.

METHODS

We induced DN in the Guizhou miniature pig with injections of streptozotocin, and P. capitata was added to the pigs' diet to treat DN. In week 16, all the animals were slaughtered, samples were collected, and the relative DN indices were measured. 16S rRNA sequencing, metagenomics, metabolomics, RNA sequencing, and proteomics were used to explore the protective mechanism of P. capitata against DN.

RESULTS

Dietary supplementation with P. capitata significantly reduced the extent of the disease, not only in term of the relative disease indices but also in hematoxylin-eosin-stained tissues. A multiomic analysis showed that two microbes (Clostridium baratii and Escherichia coli), five metabolites (oleic acid, linoleic acid, 4-phenylbutyric acid, 18-β-glycyrrhetinic acid, and ergosterol peroxide), four proteins (ENTPD5, EPHX1, ARVCF and TREH), four important mRNAs (encoding ENTPD5, EPHX1, ARVCF, and TREH), six lncRNAs (TCONS_00024194, TCONS_00085825, TCONS_00006937, TCONS_00070981, TCONS_00074099, and TCONS_00097913), and two circRNAs (novel_circ_0001514 and novel_circ_0017507) are all involved in the protective mechanism of P. capitata against DN.

CONCLUSIONS

Our results provide multidimensional theoretical support for the study and application of P. capitata.

Assembly of novel sequences for Chinese domestic pigs reveals new genes and regulatory variants providing new insights into their diversity

There is an increasing understanding that a reference genome representing an individual cannot capture all the gene repertoire of a species. Here, we conduct a population-scale missing sequences detection of Chinese domestic pigs using whole-genome sequencing data from 534 individuals. We identify 132.41 Mb of sequences absent in the reference assembly, including eight novel genes. In particular, the breeds spread in Chinese high-altitude regions perform significantly different frequencies of new sequences in promoters than other breeds. Furthermore, we dissect the role of non-coding variants and identify a novel sequence inserted in the 3'UTR of the FMO3 gene, which may be associated with the intramuscular fat phenotype. This novel sequence could be a candidate marker for meat quality. Our study provides a comprehensive overview of the missing sequences in Chinese domestic pigs and indicates that this dataset is a valuable resource for understanding the diversity and biology of pigs.

Pig pangenome graph reveals functional features of non-reference sequences

BACKGROUND

The reliance on a solitary linear reference genome has imposed a significant constraint on our comprehensive understanding of genetic variation in animals. This constraint is particularly pronounced for non-reference sequences (NRSs), which have not been extensively studied.

RESULTS

In this study, we constructed a pig pangenome graph using 21 pig assemblies and identified 23,831 NRSs with a total length of 105 Mb. Our findings revealed that NRSs were more prevalent in breeds exhibiting greater genetic divergence from the reference genome. Furthermore, we observed that NRSs were rarely found within coding sequences, while NRS insertions were enriched in immune-related Gene Ontology terms. Notably, our investigation also unveiled a close association between novel genes and the immune capacity of pigs. We observed substantial differences in terms of frequencies of NRSs between Eastern and Western pigs, and the heat-resistant pigs exhibited a substantial number of NRS insertions in an 11.6 Mb interval on chromosome X. Additionally, we discovered a 665 bp insertion in the fourth intron of the TNFRSF19 gene that may be associated with the ability of heat tolerance in Southern Chinese pigs.

CONCLUSIONS

Our findings demonstrate the potential of a graph genome approach to reveal important functional features of NRSs in pig populations.

Integrated Analysis of the Transcriptome and Microbial Diversity in the Intestine of Miniature Pig Obesity Model

Obesity, a key contributor to metabolic disorders, necessitates an in-depth understanding of its pathogenesis and prerequisites for prevention. Guangxi Bama miniature pig (GBM) offers an apt model for obesity-related studies. In this research, we used transcriptomics and 16S rRNA gene sequencing to discern the differentially expressed genes (DEGs) within intestinal (jejunum, ileum, and colon) tissues and variations in microbial communities in intestinal contents of GBM subjected to normal diets (ND) and high-fat, high-carbohydrate diets (HFHCD). After a feeding duration of 26 weeks, the HFHCD-fed experimental group demonstrated notable increases in backfat thickness, BMI, abnormal blood glucose metabolism, and blood lipid levels alongside the escalated serum expression of pro-inflammatory factors and a marked decline in intestinal health status when compared to the ND group. Transcriptomic analysis revealed a total of 1669 DEGs, of which 27 had similar differences in three intestinal segments across different groups, including five immune related genes: COL6A6, CYP1A1, EIF2AK2, NMI, and LGALS3B. Further, we found significant changes in the microbiota composition, with a significant decrease in beneficial bacterial populations within the HFHCD group. Finally, the results of integrated analysis of microbial diversity with transcriptomics show a positive link between certain microbial abundance (Solibacillus, norank_f__Saccharimonadaceae, Candidatus_Saccharimonas, and unclassified_f__Butyricicoccaceae) and changes in gene expression (COL6A6 and NMI). Overall, HFHCD appears to co-contribute to the initiation and progression of obesity in GBM by aggravating inflammatory responses, disrupting immune homeostasis, and creating imbalances in intestinal flora.

Young SINEs in pig genomes impact gene regulation, genetic diversity, and complex traits

Transposable elements (TEs) are a major source of genetic polymorphisms and play a role in chromatin architecture, gene regulatory networks, and genomic evolution. However, their functional role in pigs and contributions to complex traits are largely unknown. We created a catalog of TEs (n = 3,087,929) in pigs and found that young SINEs were predominantly silenced by histone modifications, DNA methylation, and decreased accessibility. However, some transcripts from active young SINEs showed high tissue-specificity, as confirmed by analyzing 3570 RNA-seq samples. We also detected 211,067 dimorphic SINEs in 374 individuals, including 340 population-specific ones associated with local adaptation. Mapping these dimorphic SINEs to genome-wide associations of 97 complex traits in pigs, we found 54 candidate genes (e.g., ANK2 and VRTN) that might be mediated by TEs. Our findings highlight the important roles of young SINEs and provide a supplement for genotype-to-phenotype associations and modern breeding in pigs.

Pangenome obtained by long-read sequencing of 11 genomes reveal hidden functional structural variants in pigs

Long-read sequencing (LRS) facilitates both the genome assembly and the discovery of structural variants (SVs). Here, we built a graph-based pig pangenome by incorporating 11 LRS genomes with an average of 94.01% BUSCO completeness score, revealing 206-Mb novel sequences. We discovered 183,352 nonredundant SVs (63% novel), representing 12.12% of the reference genome. By genotyping SVs in an additional 196 short-read sequencing samples, we identified thousands of population stratified SVs. Particularly, we detected 7,568 Tibetan specific SVs, some of which demonstrate significant population differentiation between Tibetan and low-altitude pigs, which might be associated with the high-altitude hypoxia adaptation in Tibetan pigs. Further integrating functional genomic data, the most promising candidate genes within the SVs that might contribute to the high-altitude hypoxia adaptation were discovered. Overall, our study generates a benchmark pangenome resource for illustrating the important roles of SVs in adaptive evolution, domestication, and genetic improvement of agronomic traits in pigs.

Generation of a genetically modified pig model with CREBRFR457Q variant

Obesity is among the strongest risk factors for type 2 diabetes (T2D). The CREBRF missense allele rs373863828 (p. Arg457Gln, p. R457Q) is associated with increased body mass index but reduced risk of T2D in people of Pacific ancestry. To investigate the functional consequences of the CREBRF variant, we introduced the corresponding human mutation R457Q into the porcine genome. The CREBRF^R457Q pigs displayed dramatically increased fat deposition, which was mainly distributed in subcutaneous adipose tissue other than visceral adipose tissue. The CREBRF^R457Q variant promoted preadipocyte differentiation. The increased differentiation capacity of precursor adipocytes conferred pigs the unique histological phenotype that adipocytes had a smaller size but a greater number in subcutaneous adipose tissue (SAT) of CREBRF^R457Q variant pigs. In addition, in SAT of CREBRF^R457Q pigs, the contents of the peroxidative metabolites 4-hydroxy-nonenal and malondialdehyde were significantly decreased, while the activity of antioxidant enzymes, such as glutathione peroxidase, superoxide dismutase, and catalase, was increased, which was in accordance with the declined level of the reactive oxygen species (ROS) in CREBRF^R457Q pigs. Together, these data supported a causal role of the CREBRF^R457Q variant in the pathogenesis of obesity, partly via adipocyte hyperplasia, and further suggested that reduced oxidative stress in adipose tissue may mediate the relative metabolic protection afforded by this variant despite the related obesity.

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.

Resveratrol Inhibits Proliferation and Differentiation of Porcine Preadipocytes by a Novel LincRNA-ROFM/miR-133b/AdipoQ Pathway

Resveratrol (RES) has a wide range of biological and pharmacological activities with various health benefits for humans as a food additive. In animal production, RES has been considered a potential functional feed additive for producing high-quality pork. Long noncoding RNAs (lncRNAs) have emerged as essential regulators of fat metabolism, and phytochemicals can regulate fat metabolism through lncRNA. However, it is unclear whether RES can improve back-fat thickness by regulating lncRNA. In this study, we identified a novel lncRNA, which was named a long intergenic non-protein coding RNA, a regulator of fat metabolism (LincRNA-ROFM), from our previous lncRNA sequencing data. LincRNA-ROFM can inhibit adipocyte proliferation and differentiation. In-depth analyses showed that LincRNA-ROFM acts as a molecular sponge for miR-133b, and adiponectin (AdipoQ) is a direct target of miR-133b in porcine preadipocytes. In addition, the expression of LincRNA-ROFM was positively correlated with AdipoQ. RES can promote the expression of LincRNA-ROFM by PPARα and C/EBPα. Altogether, our research showed that LincRNA-ROFM acts as a ceRNA to sequester miR-133b and is upregulated by RES, leading to heightened AdipoQ expression, and thus decreased adipocyte proliferation and differentiation, which reduces back-fat thickness of pigs. Taken together, the RES/LincRNA-ROFM/miR-133b/AdipoQ regulatory network preliminarily explains the mechanism of action of RES in inhibiting fat deposition, which provides new insight into the downstream mechanism of RES inhibition of fat deposits by regulating the lncRNA.

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.

Blood metabolic and physiological profiles of Bama miniature pigs at different growth stages

Background

Bama miniature pigs aged between six (6 M) and twelve months (12 M) are usually used in human medical research as laboratory pigs. However, the difference in serum metabolic profiles from 6 to 12 M-old pigs remains unclear. This study aimed to identify the metabolic and physiological profiles present in the blood to further explain changes in Bama miniature pig growth. We collected blood samples from 6 M-, eight-month- (8 M-), ten-month- (10 M-), and 12 M-old healthy Guangxi Bama miniature pigs. A total of 20 blood physiological indices (BPIs) were measured: seven for white blood cells, eight for red blood cells, and five for platelet indices. Liquid chromatography and mass spectrometry-based non-targeted metabolomic approaches were used to analyze the difference in metabolites. The associations between the differences were calculated using Spearman correlations with Benjamini–Hochberg adjustment. The 100 most abundant differential metabolites were selected for analysis of their metabolic profiles.

Results

There were no significant differences in BPIs at different ages, but the mid cell ratio and red blood cell number increased with age. Seven BPIs in Bama miniature pigs were closer to human BPIs than to mouse BPIs. A total of 14 and 25 significant differential metabolites were identified in 6 M vs. 12 M and 8 M vs. 12 M, respectively. In total, 9 and 18 amino acids and their derivatives showed significantly lower concentrations in 6 M- and 8 M-old pigs than in 12 M-old pigs. They were identified as the core significantly different metabolites between the age groups 6 M vs. 12 M and 8 M vs. 12 M. Half of the enriched pathways were the amino acids metabolism pathways. The concentration of six amino acids (dl-tryptophan, phenylacetylglycine, muramic acid, N-acetylornithine, l(−)-pipecolinic acid, and creatine) and their derivatives increased with age. A total of 61 of the top 100 most abundant metabolites were annotated. The metabolic profiles contained 14 amino acids and derivatives, six bile acids and derivatives, 19 fatty acids and derivatives, and 22 others. The concentrations of fatty acids and derivatives were found to be inversely proportional to those of amino acids and derivatives.

Conclusion

These findings suggest high levels of MID cell ratio, red blood count, and amino acids in 12 M-old pigs as indicators for improved body function over time in Bama miniature pigs, similar to those in human development. This makes the pig a more suitable medical model organism than the mouse. The results of this study are limited to the characteristics of blood metabolism in the inbred Bama miniature pigs, and the effects of impacting factors such as breed, age, sex, health status and nutritional level should be considered when studying other pig populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40813-022-00278-7.

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.

Genome-Wide Association Studies, Runs of Homozygosity Analysis, and Copy Number Variation Detection to Identify Reproduction-Related Genes in Bama Xiang Pigs

Litter size and teat number are economically important traits in the porcine industry. However, the genetic mechanisms influencing these traits remain unknown. In this study, we analyzed the genetic basis of litter size and teat number in Bama Xiang pigs and evaluated the genomic inbreeding coefficients of this breed. We conducted a genome-wide association study to identify runs of homozygosity (ROH), and copy number variation (CNV) using the novel Illumina PorcineSNP50 BeadChip array in Bama Xiang pigs and annotated the related genes in significant single nucleotide polymorphisms and common copy number variation region (CCNVR). We calculated the ROH-based genomic inbreeding coefficients (FROH) and the Spearman coefficient between FROH and reproduction traits. We completed a mixed linear model association analysis to identify the effect of high-frequency copy number variation (HCNVR; over 5%) on Bama Xiang pig reproductive traits using TASSEL software. Across eight chromosomes, we identified 29 significant single nucleotide polymorphisms, and 12 genes were considered important candidates for litter-size traits based on their vital roles in sperm structure, spermatogenesis, sperm function, ovarian or follicular function, and male/female infertility. We identified 9,322 ROHs; the litter-size traits had a significant negative correlation to FROH. A total of 3,317 CNVs, 24 CCNVR, and 50 HCNVR were identified using cnvPartition and PennCNV. Eleven genes related to reproduction were identified in CCNVRs, including seven genes related to the testis and sperm function in CCNVR1 (chr1 from 311585283 to 315307620). Two candidate genes (NEURL1 and SH3PXD2A) related to reproduction traits were identified in HCNVR34. The result suggests that these genes may improve the litter size of Bama Xiang by marker-assisted selection. However, attention should be paid to deter inbreeding in Bama Xiang pigs to conserve their genetic diversity.

The de novo assembly of a European wild boar genome revealed unique patterns of chromosomal structural variations and segmental duplications

The rapid progress of sequencing technology has greatly facilitated the de novo genome assembly of pig breeds. However, the assembly of the wild boar genome is still lacking, hampering our understanding of chromosomal and genomic evolution during domestication from wild boars into domestic pigs. Here, we sequenced and de novo assembled a European wild boar genome (ASM2165605v1) using the long‐range information provided by 10× Linked‐Reads sequencing. We achieved a high‐quality assembly with contig N50 of 26.09 Mb. Additionally, 1.64% of the contigs (222) with lengths from 107.65 kb to 75.36 Mb covered 90.3% of the total genome size of ASM2165605v1 (~2.5 Gb). Mapping analysis revealed that the contigs can fill 24.73% (93/376) of the gaps present in the orthologous regions of the updated pig reference genome (Sscrofa11.1). We further improved the contigs into chromosome level with a reference‐assistant scaffolding method. Using the ‘assembly‐to‐assembly’ approach, we identified intra‐chromosomal large structural variations (SVs, length >1 kb) between ASM2165605v1 and Sscrofa11.1 assemblies. Interestingly, we found that the number of SV events on the X chromosome deviated significantly from the linear models fitting autosomes (R² > 0.64, p < 0.001). Specifically, deletions and insertions were deficient on the X chromosome by 66.14 and 58.41% respectively, whereas duplications and inversions were excessive on the X chromosome by 71.96 and 107.61% respectively. We further used the large segmental duplications (SDs, >1 kb) events as a proxy to understand the large‐scale inter‐chromosomal evolution, by resolving parental‐derived relationships for SD pairs. We revealed a significant excess of SD movements from the X chromosome to autosomes (p < 0.001), consistent with the expectation of meiotic sex chromosome inactivation. Enrichment analyses indicated that the genes within derived SD copies on autosomes were significantly related to biological processes involving nervous system, lipid biosynthesis and sperm motility (p < 0.01). Together, our analyses of the de novo assembly of ASM2165605v1 provides insight into the SVs between European wild boar and domestic pig, in addition to the ongoing process of meiotic sex chromosome inactivation in driving inter‐chromosomal interaction between the sex chromosome and autosomes.

Genome-wide annotation of protein-coding genes in pig

BACKGROUND

There is a need for functional genome-wide annotation of the protein-coding genes to get a deeper understanding of mammalian biology. Here, a new annotation strategy is introduced based on dimensionality reduction and density-based clustering of whole-body co-expression patterns. This strategy has been used to explore the gene expression landscape in pig, and we present a whole-body map of all protein-coding genes in all major pig tissues and organs.

RESULTS

An open-access pig expression map ( www.rnaatlas.org ) is presented based on the expression of 350 samples across 98 well-defined pig tissues divided into 44 tissue groups. A new UMAP-based classification scheme is introduced, in which all protein-coding genes are stratified into tissue expression clusters based on body-wide expression profiles. The distribution and tissue specificity of all 22,342 protein-coding pig genes are presented.

CONCLUSIONS

Here, we present a new genome-wide annotation strategy based on dimensionality reduction and density-based clustering. A genome-wide resource of the transcriptome map across all major tissues and organs in pig is presented, and the data is available as an open-access resource ( www.rnaatlas.org ), including a comparison to the expression of human orthologs.

Measurement and Destruction of Porcine Endogenous Retrovirus in the Chinese Bama Minipig

Porcine hepatocytes are widely used in bioartificial liver (BAL) systems for the treatment of liver failure, and Chinese Bama minipigs (BMPs) are extensively used for animal experiments in the field of medicine in China. The genome of porcine endogenous retroviruses (PERVs) has not yet been accurately quantified, posing a threat to their clinical application because they act as a source of cells. In this study, we used genome sequence data from BMPs to predict PERV copies and their distribution. We validated and quantified the identified PERV copies and subtypes across different BMP individuals and tissues using quantitative real-time polymerase chain reaction and droplet digital polymerase chain reaction, respectively, and found that the BMP genome contains only 11 to 21 PERV copies. Notably, they lack the C subtype, which is a relatively safe cell source. Moreover, we applied CRISPR/Cas9 technology to knock out the pol fragment of PERVs in primary renal fibroblasts (PRFs) from BMPs and obtain PERV-destructed cells. Overall, our results lay a foundation for obtaining PERV-destructed BMPs as a safe source of hepatocytes for BALs for future applications.

Long-read assembly of the Chinese indigenous Ningxiang pig genome and identification of genetic variations in fat metabolism among different breeds

Advances in long-read sequencing technology and genome assembly provide an opportunity to improve the pig genome and reveal the full range of structural variations (SVs) between local Chinese and European pigs. To date, little is known about the genomes of some unique Chinese indigenous breeds, such as the Ningxiang pig. Here, we report the sequencing and assembly of a highly contiguous Ningxiang pig genome (NX) via an integration of PacBio single-molecule real-time sequencing, Illumina next-generation sequencing, BioNano optical mapping and Hi-C (chromosome conformation capture) approaches. The assembled genome comprises 2.44 Gb with a contig N50 of 26.1 Mb and 418 contigs in total. These contigs are organized into 121 scaffolds with a scaffold N50 of 139.0 Mb. More than 99.1% of the assembled sequence could be localized to 19 pseudochromosomes and is annotated with 20,914 protein-coding genes and 34.04% repetitive sequences. Comparisons between the NX and European Duroc assemblies revealed many SVs in genes involved in the immune system, nervous system, lipid metabolism and environmental adaptation. The genetic variants include 47 Chinese domestic pig-specific SVs and the associated 74 genes may contribute to the differences in domestic traits compared to European pigs. Moreover, single nucleotide polymorphisms (SNPs) identified from whole genome resequencing data of 73 Chinese pigs, representing 17 geographically isolated breeds, showed their specific genetic variations, population structure and evolutionary patterns. Finally, we explore transcriptional regulation in the first intron of the MYL4 gene, as the genomic SV (281-bp deletion) in Ningxiang pig promotes its subcutaneous fat compared to European pig breeds. This work identifies a set of Asian-specific SVs and SNPs, which will be important resources for modern pig breeding and genetic conservation.

Genome-wide single nucleotide polymorphism array and whole-genome sequencing reveal the inbreeding progression of Banna minipig inbred line

Inbred pigs are promising animal models for biomedical research and xenotransplantation. Established in 1980, the Banna minipig inbred (BMI) line originated from a sow and its own male offspring. It was selected from a small backcountry minority Lahu village, where records show that no other pig breed has ever been introduced. During the inbreeding process, we perfomed extreme inbreeding over 23 consecutive generations using full-sibling or parent-offspring mating. In order to investigate the inbreeding effects in BMI pigs across generations over the past 40 years, in this study we conducted a genome-wide SNP genotyping of the last 10 generations, representing generations 14-23. In total, we genotyped 57,746 SNPs, corresponding to an average decrease in heterozygosity rate of 0.0078 per generation. Furthermore, we were only able to identify 18,216 polymorphic loci with a MAF larger than 0.05, which is substantially lower than the values in previous reports on other pig breeds. In addition, we sequenced the genome of the first pig in the twenty-third generation (inbreeding coefficient 99.28%) to an average coverage of 12.4× to evaluate at the genome level the impact of advanced inbreeding. ROH analysis indicates that BMI pigs have longer ROHs than Wuzhishan and Duroc pigs. Those long ROH regions in BMI pigs are enriched for distinct functions compared with the highly polymorphic regions. Our study reveals a genome-wide allele diversity loss during the progress of inbreeding in BMI pigs and characterizes ROH and polymorphic regions as a result of inbreeding. Overall, our results indicate the successful establishment of the BMI line, which paves the way for further in-depth studies.

Gene cloning, prokaryotic expression and antiviral activities of interferon-αω from Chinese Bama miniature pigs

Porcine interferon (PoIFN) complex represents an ideal model for studying IFN evolution that resulted from viral pressure during domestication. IFN-αω is an emergent subtype of type I IFNs which has been primarily characterized in domestic pigs. In this study, the PoIFN-αω cDNA was cloned from Chinese Bama miniature pigs by RT-PCR, and its tissue expression profile was analyzed by real-time RT-PCR. The cDNA was expressed in Escherichia coli as a His-tagged protein and purified by nickel affinity chromatography. The antiviral activities of recombinant PoIFN-αω (rPoIFN-αω) against four different pig viruses were measured using cytopathic effect (CPE) inhibition assay. Although the PoIFN-αω sequence of Bama miniature pigs was identical to that of domestic pigs, the tissue expression profiles differed significantly between the two pig species. The rPoIFN-αω showed dose-dependent pre-infection antiviral activities against porcine pseudorabies virus, vesicular stomatitis virus and porcine reproductive and respiratory syndrome virus, but not against porcine circovirus type 2. When used as treatment post infection with the three viruses, rPoIFN-αω showed the efficacy in decreasing CPE in the infected cells in a time-dependent manner. Therefore, the expressed rPoIFN-αω could be used as an antiviral agent against pig virus infections.

The genome of the naturally evolved obesity-prone Ossabaw miniature pig

The feral pigs of Ossabaw Island (USA) have an outstanding propensity to obesity and develop complete metabolic syndrome (MetS) upon prolonged high energy dieting. We now report the first high quality genome of the Ossabaw pig with Contig N50 of ∼6.03 Mb, significantly higher than most other published pig genomes. Genomic comparison to Duroc reveals that variations including SNPs, INDELs and one ∼2 Mb inversion identified in Ossabaw pig may be related to its "thrifty" phenotype. Finally, an important positively selected gene (PSG) was found to be LEPR (leptin receptor) containing two positively selected sites which may lead to pseudogenization of this gene with possible significant effects on obesity and inflammation development. This work provides the first complete mapping of a genome representing a naturally 'feast and famine' evolved phenotype of MetS, serving as a blueprint to guide the search for new targets and new biomarkers for obesity comorbidities.

Genome-Wide Characterization and Comparative Analyses of Simple Sequence Repeats among Four Miniature Pig Breeds

Simple Summary

Simple sequence repeats (SSRs) are present at high densities in regulatory elements, suggesting that they may affect gene function and phenotypic traits. Therefore, SSRs can be exploited in marker-assisted selection. In addition, they can be widely used as molecular markers to study genetic diversity, population structure, and evolution. While SSRs have been widely studied in many mammalian species, very little research has focused on genome-wide SSRs of miniature pigs, a small but special group of pigs that express the dwarf phenotype. Based on the SSR-enriched library building and sequencing, about 30,000 novel polymorphic SSRs for four miniature pig breeds were mapped to the Duroc pig reference genome. The four miniature pig breeds had different numbers and types of SSRs and distributions of repeat units. There were 2518 polymorphic SSRs in the intron or exon regions that were common to all four breeds and functional analyses revealed 17 genes that were associated with body size and other genes that were associated with growth and development. In conclusion, the SSRs detected in the miniature pigs in this study may provide useful genetic markers for the selection of farm animals and the polymorphic SSRs provide valuable insights into the determination of mature body size, as well as the immunity, growth and development of animals.

Abstract

Simple sequence repeats (SSRs) are commonly used as molecular markers in research on genetic diversity and discrimination among taxa or breeds because polymorphisms in these regions contribute to gene function and phenotypically important traits. In this study, we investigated genome-wide characteristics, repeat units, and polymorphisms of SSRs using sequencing data from SSR-enriched libraries created from Wuzhishan (WZS), Bama (BM), inbred Luchuan (LC) and Zangxiang (ZX) miniature pig breeds. The numbers and types of SSRs, distributions of repeat units and polymorphic SSRs varied among the four breeds. Compared to the Duroc pig reference genome, 2518 polymorphic SSRs were unique and common to all four breeds and functional annotation revealed that they may affect the coding and regulatory regions of genes. Several examples, such as FGF23, MYF6, IGF1R, and LEPROT, are associated with growth and development in pigs. Three of the polymorphic SSRs were selected to confirm the polymorphism and the corresponding alleles through fluorescence polymerase chain reaction (PCR) and capillary electrophoresis. Together, this study provides useful insights into the discovery, characteristics and distribution of SSRs in four pig breeds. The polymorphic SSRs, especially those common and unique to all four pig breeds, might affect associated genes and play important roles in growth and development.

Identification of Multigene Biomarker for Shrimp White Feces Syndrome by Full-Length Transcriptome Sequencing

The pacific white shrimp, Litopenaeus vannamei, with the largest shrimp industry production in the world, is currently threatened by a severe disease, white feces syndrome (WFS), which cause devastating losses globally, while its causal agents remain largely unknown. Herein, compared to the Control shrimp by metagenomic analysis, we firstly investigated that the altered functions of intestinal microbial community in WFS shrimp were the enrichment of bacterial chemotaxis and flagellar assembly pathways, hinting at a potential role of pathogenic bacteria for growth and development, which might be related to WFS occurrence. Single-molecule real-time (SMRT) sequencing was to further identify the gene structure and gene regulation for more clues in WFS aetiology. Totally 50,049 high quality transcripts were obtained, capturing 39,995 previously mapped and 10,054 newly detected transcripts, which were annotated to 30,554 genes. A total of 158 differentially expressed genes (DEGs) were characterized in WFS shrimp. These DEGs were strongly associated with various immune related genes that regulated the expression of multiple antimicrobial peptides (e.g., antilipopolysaccharide factors, penaeidins, and crustin), which were further experimentally validated using quantitative PCR on transcript level. Collectively, multigene biomarkers were identified to be closely associated with WFS, especially those functional alterations in microbial community and the upregulated immune related gene with antibacterial activities. Our finding not only inspired our cogitation on WFS aetiology from both microbial and host immune response perspectives with combined metagenomic and full-length transcriptome sequencing, but also provided valuable information for enhancing shrimp aquaculture.

Porcine models for studying complications and organ crosstalk in diabetes mellitus

The worldwide prevalence of diabetes mellitus and obesity is rapidly increasing not only in adults but also in children and adolescents. Diabetes is associated with macrovascular complications increasing the risk for cardiovascular disease and stroke, as well as microvascular complications leading to diabetic nephropathy, retinopathy and neuropathy. Animal models are essential for studying disease mechanisms and for developing and testing diagnostic procedures and therapeutic strategies. Rodent models are most widely used but have limitations in translational research. Porcine models have the potential to bridge the gap between basic studies and clinical trials in human patients. This article provides an overview of concepts for the development of porcine models for diabetes and obesity research, with a focus on genetically engineered models. Diabetes-associated ocular, cardiovascular and renal alterations observed in diabetic pig models are summarized and their similarities with complications in diabetic patients are discussed. Systematic multi-organ biobanking of porcine models of diabetes and obesity and molecular profiling of representative tissue samples on different levels, e.g., on the transcriptome, proteome, or metabolome level, is proposed as a strategy for discovering tissue-specific pathomechanisms and their molecular key drivers using systems biology tools. This is exemplified by a recent study providing multi-omics insights into functional changes of the liver in a transgenic pig model for insulin-deficient diabetes mellitus. Collectively, these approaches will provide a better understanding of organ crosstalk in diabetes mellitus and eventually reveal new molecular targets for the prevention, early diagnosis and treatment of diabetes mellitus and its associated complications.