Identification and characterization of structural variants related to meat quality in pigs using chromosome-level genome assemblies

BACKGROUND

Many studies have been performed to identify various genomic loci and genes associated with the meat quality in pigs. However, the full genetic architecture of the trait still remains unclear in part because of the lack of accurate identification of related structural variations (SVs) which resulted from the shortage of target breeds, the limitations of sequencing data, and the incompleteness of genome assemblies. The recent generation of a new pig breed with superior meat quality, called Nanchukmacdon, and its chromosome-level genome assembly (the NCMD assembly) has provided new opportunities.

RESULTS

By applying assembly-based SV calling approaches to various genome assemblies of pigs including Nanchukmacdon, the impact of SVs on meat quality was investigated. Especially, by checking the commonality of SVs with other pig breeds, a total of 13,819 Nanchukmacdon-specific SVs (NSVs) were identified, which have a potential effect on the unique meat quality of Nanchukmacdon. The regulatory potentials of NSVs for the expression of nearby genes were further examined using transcriptome- and epigenome-based analyses in different tissues.

CONCLUSIONS

Whole-genome comparisons based on chromosome-level genome assemblies have led to the discovery of SVs affecting meat quality in pigs, and their regulatory potentials were analyzed. The identified NSVs will provide new insights regarding genetic architectures underlying the meat quality in pigs. Finally, this study confirms the utility of chromosome-level genome assemblies and multi-omics analysis to enhance the understanding of unique phenotypes.

A chromosome-level genome assembly of the Korean crossbred pig Nanchukmacdon (Sus scrofa)

As plentiful high-quality genome assemblies have been accumulated, reference-guided genome assembly can be a good approach to reconstruct a high-quality assembly. Here, we present a chromosome-level genome assembly of the Korean crossbred pig called Nanchukmacdon (the NCMD assembly) using the reference-guided assembly approach with short and long reads. The NCMD assembly contains 20 chromosome-level scaffolds with a total size of 2.38 Gbp (N50: 138.77 Mbp). Its BUSCO score is 93.1%, which is comparable to the pig reference assembly, and a total of 20,588 protein-coding genes, 8,651 non-coding genes, and 996.14 Mbp of repetitive elements are annotated. The NCMD assembly was also used to close many gaps in the pig reference assembly. This NCMD assembly and annotation provide foundational resources for the genomic analyses of pig and related species.

The AIRN lncRNA is imprinted and paternally expressed in pigs

Genomic imprinting plays critical roles during the development of mammalian species and underlying epigenetic mechanisms frequently involve long non-coding RNAs (lncRNAs). The paternal transcription of the antisense Igf2r RNA noncoding (Airn) is responsible for paternal silencing of the mouse insulin-like growth factor 2 receptor (Igf2r) gene and maternal Igf2r expression. Although the corresponding maternal DNA methylation imprint is conserved in humans and pigs, the orthologous AIRN lncRNA has been identified in humans but not in pigs. Here, we aimed to examine imprinted allelic expression of the porcine AIRN lncRNA along with a corresponding differentially methylated region (DMR) and to analyze allelic expression of AIRN and IGF2R in pigs. By comparing parthenogenetic and control porcine embryos, we identified a maternally methylated DMR and a significantly higher expression of AIRN lncRNA in control embryos (P < 0.05) indicating its paternal expression. Further analyses revealed that the expression of AIRN lncRNA was enriched in the pig brain and its subregions, and it was monoallelically expressed; whereas, IGF2R was expressed biallelically suggesting an absence of allele-specific transcriptional regulation. Our findings will lead to further investigations into the role of the imprinted porcine AIRN lncRNA during pig development.

PSXIII-13 Effect of Dietary Supplementation of Crushed Sasa Quelpaertensis Nakai on Growth Performance, Carcass Characteristics and Meat Quality in Landrace Pigs

This study was conducted to investigate the characteristics and physiological changes of growth performance by adding crushed Sasa quelpaertensis Nakai(SQN) to pig feed. The reason for using SQN is that SQN grows naturally under trees or at the foot of mountains and inhibits the growth of other plants, thereby preventing biodiversity maintenance. A total of 48 Landrace pigs with an average body weight of 45.11±6.16kg were assigned to four treatments with different levels of crushed SQN addition. The treatments were basal diet (control), basal diet supplemented with 2% SQN(T1), 4% SQN(T2) and 4% SQN(T2). The compound feeds were given ad libitum for 100 days. During the finishing period, the weight gain values in the four treatment groups were 111.37±3.66 kg, 109.40±3.27 kg, 106.10±3.28 kg and 103.50±3.27 kg, respectively. There was no significant difference in end-stage body weight, but the daily weight gain was 0.66±0.02kg, 0.64±0.02kg, 0.61±0.02kg, and 0.58±0.02kg, indicating that T1 was significantly (p&lt; 0.05) greater than that of other treatment groups. Carcass characteristics were tested at 180 days of age, and the thickness of back fat in each treatment group was 21.75±1.48mm, 20.40±1.32mm, 18.60±1.32mm, and 19.20±1.32mm. The SQN-added group showed a tendency to thin back fat. In the physicochemical analysis of pig sirloin, the fat content was also significantly decreased in the SQN-added group to 3.76±0.38%, 3.52±0.39%, 2.44±0.39% and 2.79±0.39% by treatment group. In the fatty acid composition, there was a difference in linolenic acid (C18:3n3), but there was no significant difference overall. According to the results of this study, it was difficult to observe anti-obesity effect of crushd SQN because it was not digested in pigs, a unit animal. However, if the SQN extract is used, the possibility of using it as a pig feed has been confirmed, and it is judged that it can be used as a reference for the development of feed additives.

Loss of Monoallelic Expression of IGF2 in the Adult Liver Via Alternative Promoter Usage and Chromatin Reorganization

In mammals, genomic imprinting operates via gene silencing mechanisms. Although conservation of the imprinting mechanism at the H19/IGF2 locus has been generally described in pigs, tissue-specific imprinting at the transcript level, monoallelic-to-biallelic conversion, and spatio-temporal chromatin reorganization remain largely uninvestigated. Here, we delineate spatially regulated imprinting of IGF2 transcripts, age-dependent hepatic mono- to biallelic conversion, and reorganization of topologically associating domains at the porcine H19/IGF2 locus for better translation to human and animal research. Whole-genome bisulfite sequencing (WGBS) and RNA sequencing (RNA-seq) of normal and parthenogenetic porcine embryos revealed the paternally hypermethylated H19 differentially methylated region and paternal expression of IGF2. Using a polymorphism-based approach and omics datasets from chromatin immunoprecipitation sequencing (ChIP–seq), whole-genome sequencing (WGS), RNA-seq, and Hi-C, regulation of IGF2 during development was analyzed. Regulatory elements in the liver were distinguished from those in the muscle where the porcine IGF2 transcript was monoallelically expressed. The IGF2 transcript from the liver was biallelically expressed at later developmental stages in both pigs and humans. Chromatin interaction was less frequent in the adult liver compared to the fetal liver and skeletal muscle. The duration of genomic imprinting effects within the H19/IGF2 locus might be reduced in the liver with biallelic conversion through alternative promoter usage and chromatin remodeling. Our integrative omics analyses of genome, epigenome, and transcriptome provided a comprehensive view of imprinting status at the H19/IGF2 cluster.

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.

Comparative methylation and RNA-seq expression analysis in CpG context to identify genes involved in Backfat vs. Liver diversification in Nanchukmacdon Pig

BACKGROUND

DNA methylation and demethylation at CpG islands is one of the main regulatory factors that allow cells to respond to different stimuli. These regulatory mechanisms help in developing tissue without affecting the genomic composition or undergoing selection. Liver and backfat play important roles in regulating lipid metabolism and control various pathways involved in reproductive performance, meat quality, and immunity. Genes inside these tissue store a plethora of information and an understanding of these genes is required to enhance tissue characteristics in the future generation.

RESULTS

A total of 16 CpG islands were identified, and they were involved in differentially methylation regions (DMRs) as well as differentially expressed genes (DEGs) of liver and backfat tissue samples. The genes C7orf50, ACTB and MLC1 in backfat and TNNT3, SIX2, SDK1, CLSTN3, LTBP4, CFAP74, SLC22A23, FOXC1, GMDS, GSC, GATA4, SEMA5A and HOXA5 in the liver, were categorized as differentially-methylated. Subsequently, Motif analysis for DMRs was performed to understand the role of the methylated motif for tissue-specific differentiation. Gene ontology studies revealed association with collagen fibril organization, the Bone Morphogenetic Proteins (BMP) signaling pathway in backfat and cholesterol biosynthesis, bile acid and bile salt transport, and immunity-related pathways in methylated genes expressed in the liver.

CONCLUSIONS

In this study, to understand the role of genes in the differentiation process, we have performed whole-genome bisulfite sequencing (WGBS) and RNA-seq analysis of Nanchukmacdon pigs. Methylation and motif analysis reveals the critical role of CpG islands and transcriptional factors binding site (TFBS) in guiding the differential patterns. Our findings could help in understanding how methylation of certain genes plays an important role and can be used as biomarkers to study tissue specific characteristics.

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.

PSIII-11 Effect of Breast cancer metastasis-suppressor 1-like (BRMS1L) genotype on teat number traits in Jeju native pigs (JNP) and Landrace

The purpose of this study was to compare the changes of teat number traits of Jeju native pig (JNP) and Landrace varieties according to the breast cancer metastasis-suppressor 1-like (BRMS1L) genotype. The number of teat varies according to the single-base polymorphism with a G or A base 1,087 from the start codon ATG in the exon 1 region of the BRMS1L gene. The total teat number was examined at birth for 28 JNPs and 72 Landraces, and the BRMS1L genotype was analyzed by Pyrosequencing. As a result, the genotyping frequency of JNP was 1 A / A type, 9 A / G type, and 18 G / G type, whereas, on the contrary, Landrace had 70 A / A type, 2 A / C type, and C / C was identified as 0 heads. The total teat number was between the two varieties was 13.32 ± 0.95 and 14.51 ± 1.03, respectively (P &lt; 0.001). According to genotypes in Jeju native pig breeds, the total number of nipples was 15 A / A, 14.0 ± 0.82 A / G, and 12.89 ± 0.74 G / G (P &lt; 0.018). Conversely, Landrace was 14.54 ± 1.02 A / A and 13.5 ± 0.5 A / G. As a result of the above study, the frequency of varieties and genotypes was found to have a significant effect on the difference in nipple count. Based on the results of this study, if the improvement of Jeju native pigs is to be carried out, it is expected that the number of nipples will increase and the mammalian capacity will increase.

A Comparison between Hi-C and 10X Genomics Linked Read Sequencing for Whole Genome Phasing in Hanwoo Cattle

Until recently, genome-scale phasing was limited due to the short read sizes of sequence data. Though the use of long-read sequencing can overcome this limitation, they require extensive error correction. The emergence of technologies such as 10X genomics linked read sequencing and Hi-C which uses short-read sequencers along with library preparation protocols that facilitates long-read assemblies have greatly reduced the complexities of genome scale phasing. Moreover, it is possible to accurately assemble phased genome of individual samples using these methods. Therefore, in this study, we compared three phasing strategies which included two sample preparation methods along with the Long Ranger pipeline of 10X genomics and HapCut2 software, namely 10X-LG, 10X-HapCut2, and HiC-HapCut2 and assessed their performance and accuracy. We found that the 10X-LG had the best phasing performance amongst the method analyzed. They had the highest phasing rate (89.6%), longest adjusted N50 (1.24 Mb), and lowest switch error rate (0.07%). Moreover, the phasing accuracy and yield of the 10X-LG stayed over 90% for distances up to 4 Mb and 550 Kb respectively, which were considerably higher than 10X-HapCut2 and Hi-C Hapcut2. The results of this study will serve as a good reference for future benchmarking studies and also for reference-based imputation in Hanwoo.

A functional regulatory variant of MYH3 influences muscle fiber-type composition and intramuscular fat content in pigs

Muscle development and lipid accumulation in muscle critically affect meat quality of livestock. However, the genetic factors underlying myofiber-type specification and intramuscular fat (IMF) accumulation remain to be elucidated. Using two independent intercrosses between Western commercial breeds and Korean native pigs (KNPs) and a joint linkage-linkage disequilibrium analysis, we identified a 488.1-kb region on porcine chromosome 12 that affects both reddish meat color (a*) and IMF. In this critical region, only the MYH3 gene, encoding myosin heavy chain 3, was found to be preferentially overexpressed in the skeletal muscle of KNPs. Subsequently, MYH3-transgenic mice demonstrated that this gene controls both myofiber-type specification and adipogenesis in skeletal muscle. We discovered a structural variant in the promotor/regulatory region of MYH3 for which Q allele carriers exhibited significantly higher values of a* and IMF than q allele carriers. Furthermore, chromatin immunoprecipitation and cotransfection assays showed that the structural variant in the 5'-flanking region of MYH3 abrogated the binding of the myogenic regulatory factors (MYF5, MYOD, MYOG, and MRF4). The allele distribution of MYH3 among pig populations worldwide indicated that the MYH3 Q allele is of Asian origin and likely predates domestication. In conclusion, we identified a functional regulatory sequence variant in porcine MYH3 that provides novel insights into the genetic basis of the regulation of myofiber type ratios and associated changes in IMF in pigs. The MYH3 variant can play an important role in improving pork quality in current breeding programs.

Analysis of allele-specific expression using RNA-seq of the Korean native pig and Landrace reciprocal cross

Objective

We tried to analyze allele-specific expression in the pig neocortex using bioinformatic analysis of high-throughput sequencing results from the parental genomes and offspring transcriptomes from reciprocal crosses between Korean Native and Landrace pigs.

Methods

We carried out sequencing of parental genomes and offspring transcriptomes using next generation sequencing. We subsequently carried out genome scale identification of single nucleotide polymorphisms (SNPs) in two different ways using either individual genome mapping or joint genome mapping of the same breed parents that were used for the reciprocal crosses. Using parent-specific SNPs, allele-specifically expressed genes were analyzed.

Results

Because of the low genome coverage (~4×) of the sequencing results, most SNPs were non-informative for parental lineage determination of the expressed alleles in the offspring and were thus excluded from our analysis. Consequently, 436 SNPs covering 336 genes were applicable to measure the imbalanced expression of paternal alleles in the offspring. By calculating the read ratios of parental alleles in the offspring, we identified seven genes showing allele-biased expression (p<0.05) including three previously reported and four newly identified genes in this study.

Conclusion

The newly identified allele-specifically expressing genes in the neocortex of pigs should contribute to improving our knowledge on genomic imprinting in pigs. To our knowledge, this is the first study of allelic imbalance using high throughput analysis of both parental genomes and offspring transcriptomes of the reciprocal cross in outbred animals. Our study also showed the effect of the number of informative animals on the genome level investigation of allele-specific expression using RNA-seq analysis in livestock species.

Ectopic Overexpression of Porcine Myh1 Increased in Slow Muscle Fibers and Enhanced Endurance Exercise in Transgenic Mice

Myosin heavy chain (MyHC) isoforms consist of Myh7, Myh2, Myh1, and Myh4, which are expressed in skeletal muscle tissues during postnatal development. These genes influence the contraction–relaxation activity in skeletal muscles and are involved in determining muscle composition such as the proportion of fast-to-slow and/or slow-to-fast fiber types. Among them, Myh1 is associated with skeletal muscle contraction and is involved in both slow-to-fast and fast-to-slow transition. However, the muscle transition mechanism is not well understood. For this study, we first produced porcine Myh1 transgenic (TG) mice to study whether the ectopic expressed porcine Myh1 gene had any effects on muscle composition, especially on slow-type muscle components. Our results showed that the factors associated with slow muscles, such as Myh7, Myoglobin, Troponin (slow-type units), and cytochrome C, were highly expressed in the quadriceps muscles of Myh1 transgenic mice. Furthermore, the ectopic porcine MYH1 protein was located only in the slow-type muscle fibers of the quadriceps muscles in Myh1 transgenic mice. In physical endurance tests, Myh1 transgenic mice ran longer and further on a treadmill than wild-type (WT) mice. These data fully supported our hypothesis that Myh1 is associated with slow muscle composition, with overexpression of Myh1 in muscle tissues possibly being a new key in modulating muscle fiber types. Our study provides a better understanding of muscle composition metabolism, physical mobility, and genetic factors in muscle fatigue.

Genome scan linkage analysis identifies a major quantitative trait loci for fatty acid composition in longissimus dorsi muscle in an F2 intercross between Landrace and Korean native pigs

OBJECTIVE

This study was conducted to locate quantitative trait loci (QTL) influencing fatty acid (FA) composition in a large F₂ intercross between Landrace and Korean native pigs.

METHODS

Eighteen FA composition traits were measured in more than 960 F₂ progeny. All experimental animals were genotyped with 165 microsatellite markers located throughout the pig autosomes.

RESULTS

We detected 112 QTLs for the FA composition; Forty seven QTLs reached the genome-wide significant threshold. In particular, we identified a cluster of highly significant QTLs for FA composition on SSC12. QTL for polyunsaturated fatty acid on pig chromosome 12 (F-value = 97.2 under additive and dominance model, nominal p-value 3.6×10^-39) accounted for 16.9% of phenotypic variance. In addition, four more QTLs for C18:1, C18:2, C20:4, and monounsaturated fatty acids on the similar position explained more than 10% of phenotypic variance.

CONCLUSION

Our findings of a major QTL for FA composition presented here could provide helpful information to locate causative variants to improve meat quality traits in pigs.

Association of a missense mutation in the positional candidate gene glutamate receptor-interacting protein 1 with backfat thickness traits in pigs

OBJECTIVE

Previously, we reported quantitative trait loci (QTLs) affecting backfat thickness (BFT) traits on pig chromosome 5 (SW1482-SW963) in an F2 intercross population between Landrace and Korean native pigs. The aim of this study was to evaluate glutamate receptor-interacting protein 1 (GRIP1) as a positional candidate gene underlying the QTL affecting BFT traits.

METHODS

Genotype and phenotype analyses were performed using the 1,105 F₂ progeny. A mixed-effect linear model was used to access association between these single nucleotide polymorphism (SNP) markers and the BFT traits in the F₂ intercross population.

RESULTS

Highly significant associations of two informative SNPs (c.2442 T>C, c.3316 C>G [R1106G]) in GRIP1 with BFT traits were detected. In addition, the two SNPs were used to construct haplotypes that were also highly associated with the BFT traits.

CONCLUSION

The SNPs and haplotypes of the GRIP1 gene determined in this study can contribute to understand the genetic structure of BFT traits in pigs.

Association of the Single Nucleotide Polymorphisms in RUNX1, DYRK1A, and KCNJ15 with Blood Related Traits in Pigs

The aim of this study was to detect positional candidate genes located within the support interval (SI) regions based on the results of red blood cell, mean corpuscular volume (MCV), and mean corpuscular hemoglobin quantitative trait locus (QTL) in Sus scrofa chromosome 13, and to verify the correlation between specific single-nucleotide polymorphisms (SNPs) located in the exonic region of the positional candidate gene and the three genetic traits. The flanking markers of the three QTL SI regions are SW38 and S0215. Within the QTL SI regions, 44 genes were located, and runt-related transcription factor 1, dual-specificity tyrosine-(Y)-phosphorylation regulated kinase 1A (DYRK1A), and potassium inwardly-rectifying channel, subfamily J, member 15 KCNJ15-which are reported to be related to the hematological traits and clinical features of Down syndrome-were selected as positional candidate genes. The ten SNPs located in the exonic region of the three genes were detected by next generation sequencing. A total of 1,232 pigs of an F₂ resource population between Landrace and Korean native pigs were genotyped. To investigate the effects of the three genes on each genotype, a mixed-effect model which is the considering family structure model was used to evaluate the associations between the SNPs and three genetic traits in the F₂ intercross population. Among them, the MCV level was highly significant (nominal p = 9.8×10^-9) in association with the DYRK1A-SNP1 (c.2989 G<A), and explained 4.9% of the phenotypic variance. However, since linkage disequilibrium is quite extensive in an F₂ intercross, our approach has limited power to distinguish one particular positional candidate gene from a QTL region.

Whole-transcriptome analysis of mouse adipose tissue in response to short-term caloric restriction

Caloric restriction (CR) has been shown to extend the lifespan of many species by improving cellular function and organismal health. Additionally, fat reduction by CR may play an important role in lengthening lifespan and preventing severe age-related diseases. Interestingly, CR induced the greatest transcriptome change in the epididymal fat of mice in our study. In this transcriptome analysis, we identified and categorized 446 genes that correlated with CR level. We observed down-regulation of several signaling pathways, including insulin/insulin-like growth factor 1 (insulin/IGF-1), epidermal growth factor (EGF), transforming growth factor beta (TGF-β), and canonical wingless-type mouse mammary tumor virus integration site (Wnt). Many genes related to structural features, including extracellular matrix structure, cell adhesion, and the cytoskeleton, were down-regulated, with a strong correlation to the degree of CR. Furthermore, genes related to the cell cycle and adipogenesis were down-regulated. These biological processes are well-identified targets of insulin/IGF-1, EGF, TGF-β, and Wnt signaling. In contrast, genes involved in specific metabolic processes, including the tricarboxylic acid cycle and the electron transport chain were up-regulated. We performed in silico analysis of the promoter sequences of CR-responsive genes and identified two associated transcription factors, Paired-like homeodomain 2 (Pitx2) and Paired box gene 6 (Pax6). Our results suggest that strict regulation of signaling pathways is critical for creating the optimal energy homeostasis to extend lifespan.

Production of Cloned Korean Native Pig by Somatic Cell Nuclear Transfer

The Korean native pig (KNP) have been considered as animal models for animal biotechnology research because of their relatively small body size and their presumably highly inbred status due to the closed breeding program. However, little is reported about the use of KNP for animal biotechnology researches. This study was performed to establish the somatic cell nuclear transfer (SCNT) protocol for the production of swine leukocyte antigens (SLA) homotype-defined SCNT KNP. The ear fibroblast cells originated from KNP were cultured and used as donor cell. After thawing, the donor cells were cultured for 1 hour with 15 μM roscovitine prior to the nuclear transfer. The numbers of reconstructed and parthenogenetic embryos transferred were 98 ± 35.2 and 145 ± 11.2, respectively. The pregnancy and delivery rate were 3/5 (60%) and 2/5 (40%). One healthy SLA homotype-defined SCNT KNP was successfully generated. The recipient-based individual cloning efficiency ranged from 0.65 to 1.08%. Taken together, it can be postulated that the methodological establishment of the production of SLA homotype-defined cloned KNP can be applied to the generation of transgenic cloned KNP as model animals for human disease and xenotransplantation researches.

Two cases of meningocele and meningoencephalocele in Jeju native pigs

BACKGROUND

Meningocele and meningoencephalocele of the skull are congenital deformities. Various species, such as pigs, dogs, and cats, are susceptible to congenital meningocele and meningoencephalocele and the incidence is higher in large white and landrace pigs.

CASE PRESENTATION

In this study, swelling was observed in the fontanel areas of the median planes of the skull cap in two female piglets of the same litter. Gross clinical examination, neurological examination, computed tomography (CT), and magnetic resonance imaging (MRI) were conducted on the symptomatic piglets. The gross clinical and neurological examinations revealed no specific findings, except for the swellings. According to the CT results, the length of the defect on the sagittal section of the skull was 4.7 mm in case 1 and 20.62 mm in case 2. Connected flow between the skull swellings and the cerebrospinal fluid (CSF) of the lateral ventricles was observed, and partial herniation was identified in case 2. On MRI, CSF with high T2 signals was identified in the arachnoid spaces between the cerebrum and the cerebellum in the two cases, which is consistent with intracranial hypertension. The size of the swelling formed in the parietal bones was 1.6 × 1.1 × 1.8 cm(3) (case 1) and 1.2 × 1.38 × 1.7 cm(3) (case 2). The increase in intracranial pressure was more obvious in case 2 than in case 1, and was accompanied by posterior displacements of the mesencephalon and cerebellum.

CONCLUSIONS

Case 1 was diagnosed as meningocele resulting from meningeal herniation and case 2 was diagnosed as meningoencephalocele caused by brain tissue herniation.

Genome-wide QTL analysis of meat quality-related traits in a large F2 intercross between Landrace and Korean native pigs

Background

We conducted a genome-wide linkage analysis to identify quantitative trait loci (QTL) that influence meat quality-related traits in a large F₂ intercross between Landrace and Korean native pigs. Thirteen meat quality-related traits of the m. longissimus lumborum et thoracis were measured in more than 830 F₂ progeny. All these animals were genotyped with 173 microsatellite markers located throughout the pig genome, and the GridQTL program based on the least squares regression model was used to perform the QTL analysis.

Results

We identified 23 genome-wide significant QTL in eight chromosome regions (SSC1, 2, 6, 7, 9, 12, 13, and 16) (SSC for Sus Scrofa) and detected 51 suggestive QTL in the 17 chromosome regions. QTL that affect 10 meat quality traits were detected on SSC12 and were highly significant at the genome-wide level. In particular, the QTL with the largest effect affected crude fat percentage and explained 22.5% of the phenotypic variance (F-ratio = 278.0 under the additive model, nominal P = 5.5 × 10⁻⁵⁵). Interestingly, the QTL on SSC12 that influenced meat quality traits showed an obvious trend for co-localization.

Conclusions

Our results confirm several previously reported QTL. In addition, we identified novel QTL for meat quality traits, which together with the associated positional candidate genes improve the knowledge on the genetic structure that underlies genetic variation for meat quality traits in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0080-6) contains supplementary material, which is available to authorized users.

Genome-wide association analysis to identify SNP markers affecting teat numbers in an F2 intercross population between Landrace and Korean native pigs

Most reproductive traits have low heritability and are greatly affected by environmental factors. Teat number and litter size are traits related to the reproduction ability of pigs. To identify quantitative trait loci (QTLs) for teat number traits, a genome-wide association study (GWAS) was conducted using an F2 intercross between Landrace and Korean native pigs. Genotype analysis was performed using the porcine SNP 60 K beadchip. The GWAS was performed using a mixed-effects model and linear regression approach. When a genome-wide threshold was determined using the Bonferroni method (P = 1.61 × 10(-6)), 38 single nucleotide polymorphism (SNP) markers in pig chromosome 7 (SSC7) were significantly associated with three teat number traits (total teat number, left teat number, and right teat number). Among these, SNPs in 5 genes (HDDC3, LOC100156276, LOC100155863, ANPEP, SCAMP2) were selected for further study based primarily on their statistical significance. A significant association was detected in SCAMP2 g.25280 G>A for total teat number (P = 2.0 × 10(-12)), HDDC3 g.1319 G>A SNP for left teat number (P = 2.3 × 10(-7)), and SCAMP2 g.14198 G>A for right teat number (P = 4.7 × 10(-12)). These results provide valuable information about the selective breeding for desirable teat numbers in pigs.

The Possibility of TBC1D21 as a Candidate Gene for Teat Numbers in Pigs

Based on a quantitative traits locus (QTL) study using a F₂ intercross between Landrace and Korean native pigs, a significant QTL affecting teat numbers in SSC7 was identified. The strong positional candidate gene, TBC1D21, was selected due to its biological function for epithelial mesenchymal cell development. Sequence analysis revealed six single nucleotide polymorphisms (SNPs) in the TBC1D21 gene. Among these, two SNP markers, one silent mutation (SNP01) for g.13,050A>G and one missense mutation (SNP04) for c.829A>T (S277C), were genotyped and they showed significant associations with teat number traits (p value = 6.38E-05 for SNP01 and p value = 1.06E-07 for SNP04 with total teat numbers). Further functional validation of these SNPs could give valuable information for understanding the teat number variation in pigs.

Characterization of insertional variation of porcine endogenous retroviruses in six different pig breeds

Pigs may need to be exploited as xenotransplantation donors due to the shortage of human organs, tissues and cells. Porcine endogenous retroviruses (PERVs) are a significant obstacle to xenotransplantation because they can infect human cells in vitro and have the potential for transmission of unexpected pathogens to humans. In this research, 101 pigs, including four commercial breeds (23 Berkshire, 13 Duroc, 22 Landrace and 14 Yorkshire pigs), one native breed (19 Korean native pigs) and one miniature breed (10 NIH miniature pigs) were used to investigate insertional variations for 11 PERV loci (three PERV-A, six PERV-B and two PERV-C). Over 60% of the pigs harbored one PERV-A (907F8) integration and five PERV-B (B3-3G, B3-7G, 742H1, 1155D9 and 465D1) integrations. However, two PERV-A loci (A1-6C and 1347C1) and one PERV-B locus (B3-7F) were absent in Duroc pigs. Moreover, two PERV-C loci (C2-6C and C4-2G) only existed in Korean native pigs and NIH miniature pigs. The results suggest that PERV insertional variations differ among pig breeds as well as among individuals within a breed. Also, the results presented here can be used for the selection of animals that do not have specific PERV integration for xenotransplantation research.

Comparison of live performance and meat quality parameter of cross bred (korean native black pig and landrace) pigs with different coat colors

Five hundred and forty crossbred (Korean native black pig×Landrace) F₂ were selected at a commercial pig farm and then divided into six different coat color groups: (A: Black, B: White, C: Red, D: White spot in black, E: Black spot in white, F: Black spot in red). Birth weight, 21st d weight, 140th d weight and carcass weight varied among the different coat color groups. D group (white spot in black coat) showed a significantly higher body weight at each weigh (birth weight, 140th d weight and carcass weight) than did the other groups, whereas the C group (red coat color) showed a significantly lower body weight at finishing stage (140th d weight and carcass weight) compared to other groups. Meat quality characteristics, shear force, cooking loss and meat color were not significantly different among the different coat color groups, whereas drip loss was significantly higher in F than in other groups. Most blood characteristics were not significantly different among the different groups, except for the red blood cells.

Genome-wide association study identifies quantitative trait loci affecting hematological traits in an F2 intercross between Landrace and Korean native pigs

Changes affecting the status of health and robustness can bring about physiological alterations including hematological parameters in swine. To identify quantitative trait loci (QTL) associated with eight hematological traits (one leukocyte trait, six erythrocyte traits and one platelet trait), we conducted a genome-wide association study using the PorcineSNP60K BeadChip in a resource population derived from an intercross between Landrace and Korean native pigs. A total of 36 740 SNPs from 816 F2 progeny were analyzed for each blood-related trait after filtering for quality control. Data were analyzed by the genome-wide rapid association using mixed model and regression (GRAMMAR) approach. A total of 257 significant SNPs (P < 1.36 × 10(-6) ) on SSC3, 6, 8, 13 and 17 were identified for blood-related traits in this study. Interestingly, the genomic region between 17.9 and 130 Mb on SSC8 was found to be significantly associated with red blood cell, mean corpuscular volume and mean corpuscular hemoglobin. Our results include the identification of five significant SNPs within five candidate genes (KIT, IL15, TXK, ARAP2 and ERG) for hematopoiesis. Further validation of these identified SNPs could give valuable information for understanding the variation of hematological traits in pigs.

Important role of the C-terminal region of pig aldo-keto reductase family 1 member C1 in the NADPH-dependent reduction of steroid hormones

The NADPH-dependent reduction activities of two paralogous pig AKR1C1s with and without 19 additional amino acid residues in C-terminus were evaluated against steroid hormones including 5alpha-dihydrotestosterone, testosterone, progesterone, androstenedione and 5alpha-androstane-3.17-dione, which act as substrates of the AKR1C1S. Among the hormones, the AKR1C1s exhibited the highest activity against 5alpha-dihydrotestosterone and the lowest activity against testosterone and progesterone. Furthermore, the AKR1C1s showed the largest differential activities against; 5alpha-dihydrotestosterone, but no such change of activities was found against progestrone and testosterone. These results suggest that the C-terminal region of AKR1C1 plays an important effect in the reduction activities of pig AKR1C1. Thus, the differential activities of two AKR1C1 paralogs observed in the present study provide important insights in understanding the molecular evolution.

Swine liver methylomes of Berkshire, Duroc and Landrace breeds by MeDIPS

Using a methyl-DNA immunoprecipitation technique in combination with next-generation deep sequencing, we conducted comprehensive DNA methylation profiling of liver genomes from three pig breeds: Berkshire, Duroc and Landrace. The profiles revealed that the distribution patterns of methylation signals along the genome are conserved among the three pig breeds. Specifically, many signals in coding genes were found in introns, and most signals in the repetitive elements were identified in non-long terminal repeat (LTR) retrotransposons such as long and short interspersed repetitive elements, implying a significant association with alternative splicing and expression of retrotransposable elements respectively. Differentially methylated regions among the three pig breeds were identified in the non-LTR retrotransposons, suggesting that they may lead to differential retrotransposable element activity. Altogether, this study provides advanced swine methylome data and valuable resources for understanding the function of DNA methylation in the evolutionary divergence of different pig breeds.

Worldwide genetic relationships of pigs as inferred from X chromosome SNPs

The phylogeography of the porcine X chromosome has not been studied despite the unique characteristics of this chromosome. Here, we genotyped 59 single nucleotide polymorphisms (SNPs) in 312 pigs from around the world, representing 39 domestic breeds and wild boars in 30 countries. Overall, widespread commercial breeds showed the highest heterozygosity values, followed by African and American populations. Structuring, as inferred from FST and analysis of molecular variance, was consistently larger in the non-pseudoautosomal (NPAR) than in the pseudoautosomal regions (PAR). Our results show that genetic relationships between populations can vary widely between the NPAR and the PAR, underscoring the fact that their genetic trajectories can be quite different. NPAR showed an increased commercial-like genetic component relative to the PAR, probably because human selection processes to obtain individuals with high productive parameters were mediated by introgressing boars rather than sows.

QTL analysis of clinical-chemical traits in an F₂ intercross between Landrace and Korean native pigs

Clinical-chemical traits are essential when examining the health status of individuals. The aim of this study was to identify quantitative trait loci (QTL) and the associated positional candidate genes affecting clinical-chemical traits in a reciprocal F(2) intercross between Landrace and Korean native pigs. Following an overnight fast, 25 serum phenotypes related to clinical-chemical traits (e.g., hepatic function parameters, renal function parameters, electrolyte, lipids) were measured in >970 F(2) progeny. All experimental samples were subjected to genotyping analysis using 165 microsatellite markers located across the genome. We identified eleven genome-wide significant QTL in six chromosomal regions (SSC 2, 7, 8, 13, 14, and 15) and 59 suggestive QTL in 17 chromosomal regions (SSC 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13, 14, 15, 16, 17, and 18). We also observed significant effects of reciprocal crosses on some of the traits, which would seem to result from maternal effect, QTL on sex chromosomes, imprinted genes, or genetic difference in mitochondrial DNA. The role of genomic imprinting in clinical-chemical traits also was investigated. Genome-wide analysis revealed a significant evidence for an imprinted QTL in SSC4 affecting serum amylase levels. Additionally, a series of bivariate linkage analysis provided strong evidence that QTL in SSC 2, 13, 15, and 18 have a pleiotropic effect on clinical-chemical traits. In conclusion, our study detected both novel and previously reported QTL influencing clinical-chemical traits in pigs. The identified QTL together with the positional candidate genes identified here could play an important role in elucidating the genetic structure of clinical-chemical phenotype variation in humans and swine.

Novel alternative splicing by exon skipping in KIT associated with whole-body roan in an intercrossed population of Landrace and Korean Native pigs

The KIT locus has been suggested to be a strong candidate region linked with whole-body roan in the F(2) population produced by intercrosses between Landrace and Korean Native pigs. In this manuscript, we report the finding of a novel alternative splicing event in the porcine KIT gene that results in the skipping of exon 5 in the I(Rn) allele. KIT mRNAs that lack exon 5 were identified in the large intestine and skin, suggesting that the mechanism responsible for the skipping of exon 5 may be tissue specific. A U(26) repeat in intron 5 showed complete linkage (LOD = 11.8) with the roan phenotype and absolute association with the black phenotype of the Korean Native pig (KNP) population samples, inferring that the repeat pattern may alter the complementary base-pairing-mediated looping-out of introns 4 and 5, which may mediate the exon 5-skipping event. Although the sample size in our study was relatively small, we speculate that the R3 allele containing the U(26) repeat is a causative element for the roan phenotype via alternative control of the exon skipping in our roan pedigree.

QTL analysis of white blood cell, platelet and red blood cell-related traits in an F2 intercross between Landrace and Korean native pigs

Haematological traits play important roles in disease resistance and defence functions. The objective of this study was to locate quantitative trait loci (QTL) and the associated positional candidate genes influencing haematological traits in an F(2) intercross between Landrace and Korean native pigs. Eight blood-related traits (six erythrocyte traits, one leucocyte trait and one platelet trait) were measured in 816 F(2) progeny. All experimental animals were genotyped with 173 informative microsatellite markers located throughout the pig genome. We report that nine chromosomes harboured QTL for the baseline blood parameters: genomic regions on SSC 1, 4, 5, 6, 8, 9, 11, 13 and 17. Eight of twenty identified QTL reached genome-wide significance. In addition, we evaluated the KIT locus, an obvious candidate gene locus affecting variation in blood-related traits. Using dense single nucleotide polymorphism marker data on SSC 8 and the marker-assisted association test, the strong association of the KIT locus with blood phenotypes was confirmed. In conclusion, our study identified both previously reported and novel QTL affecting baseline haematological parameters in pigs. Additionally, the positional candidate genes identified here could play an important role in elucidating the genetic architecture of haematological phenotype variation in swine and in humans.

Establishment of a resource population of SLA haplotype-defined Korean native pigs

The highly polymorphic porcine major histocompatibility complex (MHC), or the swine leukocyte antigens (SLA), has been repeatedly associated with variations in swine immune response to pathogens and vaccines as well as with production traits. The SLA antigens are also important targets for immunological recognition of foreign tissue grafts. We recently established a resource population of Korean native pigs as models for human transplantation and xenotransplantation research. In this study, 115 animals derived from three generations of the Korean native pigs were genotyped for three SLA class I (SLA-2, SLA-3 and SLA-1) and three SLA class II loci (DRB1, DQB1, DQA) using PCR with sequence-specific primers (PCR-SSP) at the allele group resolution. A total of seven SLA haplotypes (Lr-5.34, Lr-7.23, Lr-31.13, Lr-56.23, Lr-56.30, Lr-59.1, Lr-65.34), comprising six unique class I and five unique class II haplotypes, were characterized in the founding animals. Class I haplotype Lr-65.0 and class II haplotype Lr-0.34 were novel; and together with Lr-56.0 these haplotypes appeared to be breed-specific. In the progeny population, Lr-7.23 and Lr-56.30 appeared to be the most prevalent haplotypes with frequencies of 34.7% and 31.6%, respectively; the overall homozygosity was 27.4%. This resource population of SLA-defined Korean native pigs will be useful as large animal models for various transplantation and xenotransplantation experiments, as well as for dissecting the roles of SLA proteins in swine disease resistance and production traits.

The robust phylogeny of Korean wild boar (Sus scrofa coreanus) using partial D-loop sequence of mtDNA

In order to elucidate the precise phylogenetic relationships of Korean wild boar (Sus scrofa coreanus), a partial mtDNA D-loop region (1,274 bp, NC_000845 nucleotide positions 16576-1236) was sequenced among 56 Korean wild boars. In total, 25 haplotypes were identified and classified into four distinct subgroups (K1 to K4) based on Bayesian phylogenetic analysis using Markov chain Monte Carlo methods. An extended analysis, adding 139 wild boars sampled worldwide, confirmed that Korean wild boars clearly belong to the Asian wild boar cluster. Unexpectedly, the Myanmarese/Thai wild boar population was detected on the same branch as Korean wild boar subgroups K3 and K4. A parsimonious median-joining network analysis including all Asian wild boar haplotypes again revealed four maternal lineages of Korean wild boars, which corresponded to the four Korean wild boar subgroups identified previously. In an additional analysis, we supplemented the Asian wild boar network with 34 Korean and Chinese domestic pig haplotypes. We found only one haplotype, C31, that was shared by Chinese wild, Chinese domestic and Korean domestic pigs. In contrast to our expectation that Korean wild boars contributed to the gene pool of Korean native pigs, these data clearly suggest that Korean native pigs would be introduced from China after domestication from Chinese wild boars.

Integrating Y-chromosome, mitochondrial, and autosomal data to analyze the origin of pig breeds

We have investigated the origin of swine breeds through the joint analysis of mitochondrial, microsatellite, and Y-chromosome polymorphisms in a sample of pigs and wild boars with a worldwide distribution. Genetic differentiation between pigs and wild boars was remarkably weak, likely as a consequence of a sustained gene flow between both populations. The analysis of nuclear markers evidenced the existence of a close genetic relationship between Near Eastern and European wild boars making it difficult to infer their relative contributions to the gene pool of modern European breeds. Moreover, we have shown that European and Far Eastern pig populations have contributed maternal and paternal lineages to the foundation of African and South American breeds. Although West African pigs from Nigeria and Benin exclusively harbored European alleles, Far Eastern and European genetic signatures of similar intensity were detected in swine breeds from Eastern Africa. This region seems to have been a major point of entry of livestock species in the African continent as a result of the Indian Ocean trade. Finally, South American creole breeds had essentially a European ancestry although Asian Y-chromosome and mitochondrial haplotypes were found in a few Nicaraguan pigs. The existence of Spanish and Portuguese commercial routes linking Asia with America might have favored the introduction of Far Eastern breeds into this continent.

Genome-wide identification of quantitative trait loci in a cross between Hampshire and Landrace II: meat quality traits

Background

Meat quality traits are important in pig breeding programs, but they are difficult to include in a traditional selection program. Marker assisted selection (MAS) of meat quality traits is therefore of interest in breeding programs and a Quantitative Trait Locus (QTL) analysis is the key to identifying markers that can be used in MAS. In this study, Landrace and Hampshire intercross and backcross families were used to investigate meat quality traits. Hampshire pigs are commonly used as the sire line in commercial pig breeding. This is the first time a pedigree including Hampshire pigs has been used for a QTL analysis of meat quality traits.

Results

In total, we analyzed 39 meat quality traits and identified eight genome-wide significant QTL peaks in four regions: one on chromosome 3, two on chromosome 6 and one on chromosome 16. At least two of the QTLs do not appear to have been detected in previous studies. On chromosome 6 we identified QTLs for water content in M. longissimus dorsi (LD), drip loss in LD and post mortem pH decline in LD. On chromosomes 3 and 16 we identified previously undetected QTLs for protein content in LD and for freezing and cooking loss respectively.

Conclusion

We identified at least two new meat quality trait QTLs at the genome-wide significance level. We detected two QTLs on chromosome 6 that possibly coincide with QTLs detected in other studies. We were also able to exclude the C1843T mutation in the ryanodine receptor (RYR1) as a causative mutation for one of the chromosome 6 QTLs in this cross.