High-resolution physical mapping and construction of a porcine contig spanning the intramuscular fat content QTL

Genome-Wide Association Study Reveals Additive and Non-Additive Effects on Growth Traits in Duroc Pigs

Growth rate plays a critical role in the pig industry and is related to quantitative traits controlled by many genes. Here, we aimed to identify causative mutations and candidate genes responsible for pig growth traits. In this study, 2360 Duroc pigs were used to detect significant additive, dominance, and epistatic effects associated with growth traits. As a result, a total number of 32 significant SNPs for additive or dominance effects were found to be associated with various factors, including adjusted age at a specified weight (AGE), average daily gain (ADG), backfat thickness (BF), and loin muscle depth (LMD). In addition, the detected additive significant SNPs explained 2.49%, 3.02%, 3.18%, and 1.96% of the deregressed estimated breeding value (DEBV) variance for AGE, ADG, BF, and LMD, respectively, while significant dominance SNPs could explain 2.24%, 13.26%, and 4.08% of AGE, BF, and LMD, respectively. Meanwhile, a total of 805 significant epistatic effects SNPs were associated with one of ADG, AGE, and LMD, from which 11 sub-networks were constructed. In total, 46 potential genes involved in muscle development, fat deposition, and regulation of cell growth were considered as candidates for growth traits, including CD55 and NRIP1 for AGE and ADG, TRIP11 and MIS2 for BF, and VRTN and ZEB2 for LMD, respectively. Generally, in this study, we detected both new and reported variants and potential candidate genes for growth traits of Duroc pigs, which might to be taken into account in future molecular breeding programs to improve the growth performance of pigs.

Identification and Characterization of Circular RNAs in Association With the Deposition of Intramuscular Fat in Aohan Fine-Wool Sheep

Aohan fine-wool sheep (AFWS) is a high-quality fine-wool sheep breed that supplies wool and meat. Research is needed on the molecular mechanism behind intramuscular fat (IMF) deposition that greatly improves mutton quality. The widely expressed non-coding RNA is physiologically used in roles such as competitive endogenous RNA (ceRNA) that includes circular RNAs (circRNAs). Although circRNAs were studied in many fields, little research was devoted to IMF in sheep. We used the longissimus dorsi muscle of 2 and 12-month-old AWFS as research material to identify circRNAs related to IMF deposition in these sheep by RNA-seq screening for differentially expressed circRNAs in the two age groups. A total of 11,565 candidate circRNAs were identified, of which the 104 differentially expressed circRNAs in the two age groups were analyzed. Enrichment analysis was performed using Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes. The enriched pathways included lipid transport (GO:0006869), negative regulation of canonical Wnt signaling pathway (GO:0090090), fat digestion and absorption (ko04975), and sphingolipid metabolism (ko00600). The differentially expressed circRNAs included ciRNA455, circRNA9086, circRNA7445, circRNA4557, and others. The source genes involved in these pathways might regulate IMF deposition. We used the TargetScan and miRanda software for interaction analysis, and a network diagram of circRNA-miRNA interactions was created. CircRNA455-miR-127, circRNA455-miR-29a, circRNA455-miR-103, circRNA4557-mir149-5p, and circRNA2440-mir-23a might be involved in the IMF deposition process. The targeting relationship of circRNA4557-miR-149-5p was verified by a dual-luciferase reporter assay. The RT-qPCR results of seven randomly selected circRNAs were consistent with the sequencing results. This study provides additional information on circRNA regulation of IMF deposition in AFWS and is a useful resource for future research on this sheep breed.

Genomic structure, expression and association study of the porcine FSD2

The fibronectin type III and SPRY domain containing 2 (FSD2) on porcine chromosome 7 is considered a candidate gene for pork quality, since its two domains, which were present in fibronectin and ryanodine receptor. The fibronectin type III and SPRY domains were first identified in fibronectin and ryanodine receptor, respectively, which are candidate genes for meat quality. The aim of this study was to elucidate the genomic structure of FSD2 and functions of single nucleotide polymorphisms (SNPs) within FSD2 that are related to meat quality in pigs. Using a bacterial artificial chromosome clone sequence, we revealed that porcine FSD2 consisted of 13 exons encoding 750 amino acids. In addition, FSD2 was expressed in heart, longissimus dorsi muscle, psoas muscle, and tendon among 23 kinds of porcine tissues tested. A total of ten SNPs, including four missense mutations, were identified in the exonic region of FSD2, and two major haplotypes were obtained based on the SNP genotypes of 633 Berkshire pigs. Both haplotypes were associated significantly with intramuscular fat content (IMF, P < 0.020) and moisture percentage (MP, P < 0.002). Moreover, haplotype 2 was associated with meat color, affecting yellowness (P = 0.002). These haplotype effects were further supported by the alteration of putative protein structures with amino acid substitutions. Taken together, our results suggest that FSD2 haplotypes are involved in regulating meat quality including IMF, MP, and meat color in pigs, and may be used as meaningful molecular makers to identify pigs with preferable pork quality.

KCNK10, a tandem pore domain potassium channel, is a regulator of mitotic clonal expansion during the early stage of adipocyte differentiation

KCNK10, a member of tandem pore domain potassium channel family, gives rise to leak K⁺ currents. It plays important roles in stabilizing the negative resting membrane potential and in counterbalancing depolarization. We previously demonstrated that kcnk10 expression is quickly elevated during the early stage of adipogenesis of 3T3-L1 cells and that reduction of kcnk10 expression inhibits adipocyte differentiation. However, the molecular mechanism of KCNK10 in adipocyte differentiation remains unclear. Here we revealed that kcnk10 is induced by 3-isobutyl-1-methylxanthine, a cyclic nucleotide phosphodiesterase inhibitor and a potent inducer of adipogenesis, during the early stage of adipocyte differentiation. We also demonstrated that KCNK10 functions as a positive regulator of mitotic clonal expansion (MCE), a necessary process for terminal differentiation. The reduction of kcnk10 expression repressed the expression levels of CCAAT/enhancer-binding protein β (C/EBPβ) and C/EBPδ as well as the phosphorylation level of Akt during the early phase of adipogenesis. In addition, knockdown of kcnk10 expression suppressed insulin-induced Akt phosphorylation. These results indicate that KCNK10 contributes to the regulation of MCE through the control of C/EBPβ and C/EBPδ expression and insulin signaling.

Genome-wide QTL mapping of nine body composition and bone mineral density traits in pigs

Background

Since the pig is one of the most important livestock animals worldwide, mapping loci that are associated with economically important traits and/or traits that influence animal welfare is extremely relevant for efficient future pig breeding. Therefore, the purpose of this study was a genome-wide mapping of quantitative trait loci (QTL) associated with nine body composition and bone mineral traits: absolute (Fat, Lean) and percentage (FatPC, LeanPC) fat and lean mass, live weight (Weight), soft tissue X-ray attenuation coefficient (R), absolute (BMC) and percentage (BMCPC) bone mineral content and bone mineral density (BMD).

Methods

Data on the nine traits investigated were obtained by Dual-energy X-ray absorptiometry for 551 pigs that were between 160 and 200 days old. In addition, all pigs were genotyped using Illumina’s PorcineSNP60 Genotyping BeadChip. Based on these data, a genome-wide combined linkage and linkage disequilibrium analysis was conducted. Thus, we used 44 611 sliding windows that each consisted of 20 adjacent single nucleotide polymorphisms (SNPs). For the middle of each sliding window a variance component analysis was carried out using ASReml. The underlying mixed linear model included random QTL and polygenic effects, with fixed effects of sex, housing, season and age.

Results

Using a Bonferroni-corrected genome-wide significance threshold of P < 0.001, significant peaks were identified for all traits except BMCPC. Overall, we identified 72 QTL on 16 chromosomes, of which 24 were significantly associated with one trait only and the remaining with more than one trait. For example, a QTL on chromosome 2 included the highest peak across the genome for four traits (Fat, FatPC, LeanPC and R). The nearby gene, ZNF608, is known to be associated with body mass index in humans and involved in starvation in Drosophila, which makes it an extremely good candidate gene for this QTL.

Conclusions

Our QTL mapping approach identified 72 QTL, some of which confirmed results of previous studies in pigs. However, we also detected significant associations that have not been published before and were able to identify a number of new and promising candidate genes, such as ZNF608.

Electronic supplementary material

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

Galaxy tools to study genome diversity

Background

Intra-species genetic variation can be used to investigate population structure, selection, and gene flow in non-model vertebrates; and due to the plummeting costs for genome sequencing, it is now possible for small labs to obtain full-genome variation data from their species of interest. However, those labs may not have easy access to, and familiarity with, computational tools to analyze those data.

Results

We have created a suite of tools for the Galaxy web server aimed at handling nucleotide and amino-acid polymorphisms discovered by full-genome sequencing of several individuals of the same species, or using a SNP genotyping microarray. In addition to providing user-friendly tools, a main goal is to make published analyses reproducible. While most of the examples discussed in this paper deal with nuclear-genome diversity in non-human vertebrates, we also illustrate the application of the tools to fungal genomes, human biomedical data, and mitochondrial sequences.

Conclusions

This project illustrates that a small group can design, implement, test, document, and distribute a Galaxy tool collection to meet the needs of a particular community of biologists.

Genome-wide association study of meat quality traits in a White Duroc×Erhualian F2 intercross and Chinese Sutai pigs

Thousands of QTLs for meat quality traits have been identified by linkage mapping studies, but most of them lack precise position or replication between populations, which hinder their application in pig breeding programs. To localize QTLs for meat quality traits to precise genomic regions, we performed a genome-wide association (GWA) study using the Illumina PorcineSNP60K Beadchip in two swine populations: 434 Sutai pigs and 933 F2 pigs from a White Duroc×Erhualian intercross. Meat quality traits, including pH, color, drip loss, moisture content, protein content and intramuscular fat content (IMF), marbling and firmness scores in the M. longissimus (LM) and M. semimembranosus (SM) muscles, were recorded on the two populations. In total, 127 chromosome-wide significant SNPs for these traits were identified. Among them, 11 SNPs reached genome-wise significance level, including 1 on SSC3 for pH, 1 on SSC3 and 3 on SSC15 for drip loss, 3 (unmapped) for color a*, and 2 for IMF each on SSC9 and SSCX. Except for 11 unmapped SNPs, 116 significant SNPs fell into 28 genomic regions of approximately 10 Mb or less. Most of these regions corresponded to previously reported QTL regions and spanned smaller intervals than before. The loci on SSC3 and SSC7 appeared to have pleiotropic effects on several related traits. Besides them, a few QTL signals were replicated between the two populations. Further, we identified thirteen new candidate genes for IMF, marbling and firmness, on the basis of their positions, functional annotations and reported expression patterns. The findings will contribute to further identification of the causal mutation underlying these QTLs and future marker-assisted selection in pigs.

Segregation of regulatory polymorphisms with effects on the gluteus medius transcriptome in a purebred pig population

Background

The main goal of the present study was to analyse the genetic architecture of mRNA expression in muscle, a tissue with an outmost economic importance for pig breeders. Previous studies have used F₂ crosses to detect porcine expression QTL (eQTL), so they contributed with data that mostly represents the between-breed component of eQTL variation. Herewith, we have analysed eQTL segregation in an outbred Duroc population using two groups of animals with divergent fatness profiles. This approach is particularly suitable to analyse the within-breed component of eQTL variation, with a special emphasis on loci involved in lipid metabolism.

Methodology/Principal Findings

GeneChip Porcine Genome arrays (Affymetrix) were used to determine the mRNA expression levels of gluteus medius samples from 105 Duroc barrows. A whole-genome eQTL scan was carried out with a panel of 116 microsatellites. Results allowed us to detect 613 genome-wide significant eQTL unevenly distributed across the pig genome. A clear predominance of trans- over cis-eQTL, was observed. Moreover, 11 trans-regulatory hotspots affecting the expression levels of four to 16 genes were identified. A Gene Ontology study showed that regulatory polymorphisms affected the expression of muscle development and lipid metabolism genes. A number of positional concordances between eQTL and lipid trait QTL were also found, whereas limited evidence of a linear relationship between muscle fat deposition and mRNA levels of eQTL regulated genes was obtained.

Conclusions/Significance

Our data provide substantial evidence that there is a remarkable amount of within-breed genetic variation affecting muscle mRNA expression. Most of this variation acts in trans and influences biological processes related with muscle development, lipid deposition and energy balance. The identification of the underlying causal mutations and the ascertainment of their effects on phenotypes would allow gaining a fundamental perspective about how complex traits are built at the molecular level.

Joint analysis of quantitative trait loci and major-effect causative mutations affecting meat quality and carcass composition traits in pigs

Background

Detection of quantitative trait loci (QTLs) affecting meat quality traits in pigs is crucial for the design of efficient marker-assisted selection programs and to initiate efforts toward the identification of underlying polymorphisms. The RYR1 and PRKAG3 causative mutations, originally identified from major effects on meat characteristics, can be used both as controls for an overall QTL detection strategy for diversely affected traits and as a scale for detected QTL effects. We report on a microsatellite-based QTL detection scan including all autosomes for pig meat quality and carcass composition traits in an F2 population of 1,000 females and barrows resulting from an intercross between a Pietrain and a Large White-Hampshire-Duroc synthetic sire line. Our QTL detection design allowed side-by-side comparison of the RYR1 and PRKAG3 mutation effects seen as QTLs when segregating at low frequencies (0.03-0.08), with independent QTL effects detected from most of the same population, excluding any carrier of these mutations.

Results

Large QTL effects were detected in the absence of the RYR1 and PRKGA3 mutations, accounting for 12.7% of phenotypic variation in loin colour redness CIE-a* on SSC6 and 15% of phenotypic variation in glycolytic potential on SSC1. We detected 8 significant QTLs with effects on meat quality traits and 20 significant QTLs for carcass composition and growth traits under these conditions. In control analyses including mutation carriers, RYR1 and PRKAG3 mutations were detected as QTLs, from highly significant to suggestive, and explained 53% to 5% of the phenotypic variance according to the trait.

Conclusions

Our results suggest that part of muscle development and backfat thickness effects commonly attributed to the RYR1 mutation may be a consequence of linkage with independent QTLs affecting those traits. The proportion of variation explained by the most significant QTLs detected in this work is close to the influence of major-effect mutations on the least affected traits, but is one order of magnitude lower than effect on variance of traits primarily affected by these causative mutations. This suggests that uncovering physiological traits directly affected by genetic polymorphisms would be an appropriate approach for further characterization of QTLs.

RNA interference-mediated knockdown of the mouse gene encoding potassium channel subfamily K member 10 inhibits hormone-induced differentiation of 3T3-L1 preadipocytes

Intramuscular fat (IMF) is an economically important trait of domestic meat animals; thus, it is important to identify the factors that influence the IMF content. In this study, we identified the gene associated with adipogenesis from all the positional candidate genes located in the quantitative trait loci (QTL) for IMF content on porcine chromosome 7. We analyzed the expression of the abovementioned genes during differentiation of mouse 3T3-L1 preadipocytes by using real-time polymerase chain reaction (PCR). Total cellular RNA was extracted before and 6, 12, 36, and 48 h and 4, 6, and 8d after treatment with standard hormonal inducers of differentiation-insulin, dexamethasone, and 3-isobutyl-1-methylxanthine (IBMX). Six hours after induction, potassium channel subfamily K member 10 (KCNK10) gene expression in the preadipocytes was found to be 100-fold greater than that at the baseline; this expression declined until day 4 after the induction. Moreover, knockdown of the KCNK10 gene by transfection with short-hairpin RNA (shRNA) significantly decreased triacylglycerol accumulation on day 8 after the induction. An RNA interference study revealed that KCNK10 knockdown inhibited the differentiation of 3T3-L1 cells. These results indicate that KCNK10 plays an important role in the early stages of preadipocyte differentiation.

A BAC contig map over the proximal approximately 3.3 Mb region of horse chromosome 21

A total of 207 BAC clones containing 155 loci were isolated and arranged into a map of linearly ordered overlapping clones over the proximal part of horse chromosome 21 (ECA21), which corresponds to the proximal half of the short arm of human chromosome 19 (HSA19p) and part of HSA5. The clones form two contigs - each corresponding to the respective human chromosomes - that are estimated to be separated by a gap of approximately 200 kb. Of the 155 markers present in the two contigs, 141 (33 genes and 108 STS) were generated and mapped in this study. The BACs provide a 4-5x coverage of the region and span an estimated length of approximately 3.3 Mb. The region presently contains one mapped marker per 22 kb on average, which represents a major improvement over the previous resolution of one marker per 380 kb obtained through the generation of a dense RH map for this segment. Dual color fluorescence in situ hybridization on metaphase and interphase chromosomes verified the relative order of some of the BACs and helped to orient them accurately in the contigs. Despite having similar gene order and content, the equine region covered by the contigs appears to be distinctly smaller than the corresponding region in human (3.3 Mb vs. 5.5-6 Mb) because the latter harbors a host of repetitive elements and gene families unique to humans/primates. Considering limited representation of the region in the latest version of the horse whole genome sequence EquCab2, the dense map developed in this study will prove useful for the assembly and annotation of the sequence data on ECA21 and will be instrumental in rapid search and isolation of candidate genes for traits mapped to this region.

Identification of QTL with effects on intramuscular fat content and fatty acid composition in a Duroc x Large White cross

Background

Improving pork quality can be done by increasing intramuscular fat (IMF) content. This trait is influenced by quantitative trait loci (QTL) sought out in different pig populations. Considering the high IMF content observed in the Duroc pig, it was appealing to determine whether favourable alleles at a major gene or QTL could be found. The detection was performed in an experimental F2 Duroc × Large White population first by segregation analysis, then by QTL mapping using additional molecular information.

Results

Segregation analysis provided evidence for a major gene, with a recessive Duroc allele increasing IMF by 1.8% in Duroc homozygous pigs. However, results depended on whether data were normalised or not. After Box-Cox transformation, likelihood ratio was indeed 12 times lower and no longer significant. The QTL detection results were partly consistent with the segregation analysis. Three QTL significant at the chromosome wide level were evidenced. Two QTL, located on chromosomes 13 and 15, showed a high IMF Duroc recessive allele with an overall effect slightly lower than that expected from segregation analysis (+0.4 g/100 g muscle). The third QTL was located on chromosome 1, with a dominant Large White allele inducing high IMF content (+0.5 g/100 g muscle). Additional QTL were detected for muscular fatty acid composition.

Conclusion

The study presented results from two complementary approaches, a segregation analysis and a QTL detection, to seek out genes involved in the higher IMF content observed in the Duroc population. Discrepancies between both methods might be partially explained by the existence of at least two QTL with similar characteristics located on two different chromosomes for which different boars were heterozygous. The favourable and dominant allele detected in the Large White population was unexpected. Obviously, in both populations, the favourable alleles inducing high IMF content were not fixed and improving IMF by fixing favourable alleles using markers can then be applied both in Duroc and LW populations. With QTL affecting fatty acid composition, combining an increase of IMF content enhancing monounsaturated fatty acid percentage would be of great interest.