Taking advantage from phenotype variability in a local animal genetic resource: identification of genomic regions associated with the hairless phenotype in Casertana pigs

Casertana is an endangered autochthonous pig breed (raised in south‐central Italy) that is considered to be the descendant of the influential Neapolitan pig population that was used to improve British breeds in the 19th century. Casertana pigs are characterized by a typical, almost complete, hairless phenotype, even though a few Casertana pigs are normal haired. In this work, using Illumina PorcineSNP60 BeadChip data, we carried out a genome‐wide association study and an F_ST analysis with this breed by comparing animals showing the classical hairless phenotype (n = 81) versus pigs classified as haired (n = 15). Combining the results obtained with the two approaches, we identified two significant regions: one on porcine chromosome (SSC) 7 and one on SSC15. The SSC7 region contains the forkhead box N3 (FOXN3) gene, the most plausible candidate gene of this region, considering that mutations in another gene of the same family (forkhead box N1; Foxn1 or FOXN1) are responsible for the nude locus in rodents and alopecia in humans. Another potential candidate gene, rho guanine nucleotide exchange factor 10 (ARHGEF10), is located in the SSC15 region. FOXN3 and ARHGEF10 have been detected as differentially expressed in androgenetic and senescent alopecia respectively. This study on an autochthonous pig breed contributes to shed some light on novel genes potentially involved in hair development and growth and demonstrates that local animal breeds can be valuable genetic resources for disclosing genetic factors affecting unique traits, taking advantage of phenotype variability segregating in small populations.     

Identification of genes related to intramuscular fat independent of backfat thickness in Duroc pigs using single‐step genome‐wide association

Intramuscular fat (IMF) is an important meat‐quality trait of pigs, which influences pork’s shearing force, hydraulics, tenderness and juicy flavor. However, to achieve a higher percentage of lean meat, pigs with lower backfat thickness (BF) are intensively selected for, which may lead to a reduction in pork quality. Therefore, the objective of this study was to locate loci that affect IMF without changing BF. A single‐step GWAS was performed on 950 Duroc pigs genotyped by a 50K SNP chip in order to detect genomic variants relevant to IMF and BF. The significant SNPs detected were afterwards divided into a BF subset (seven SNPs), an IMF subset (11 SNPs) and a subset of both traits (12 SNPs), according to their P‐value and LD. After SNP and QTL annotation, our results indicated that SSC1: 167938652, 166363826, 164829874 and 167171587 might be associated with IMF without changing BF. In the subset of both traits, we found that the combined effect of ALGA0006602 (SSC1: 159538854) and 12784636 (SSC1: 160773437) might improve the IMF without changing BF. Our gene annotation result showed that TLE3, ITGA11, SMAD6, PAQR5 and [RNF152^A/G × MC4R^A/A] genes might affect IMF independently of BF. We believe that the SNPs and genes identified in this study will be valuable for the future molecular breeding of IMF in Duroc pigs.     

Quantitative gene expression analysis on chromosome 6 between Korean native pigs and Yorkshire breeds for fat deposition

Pig chromosome 6 (SSC6) has been reported to have QTL affecting backfat thickness (BFT) and intramuscular fat (IMF). A human-pig comparative map covering 18 autosomes with the highest resolution has been constructed and based on this map SSC6 has conserved syntenicgroups with human chromosome (HSA) 16, 19, 1, and 18. In this study, the pig Affy elements mapped to the SSC6 were analyzed, and the differentially expressed genes in three tissues (liver, backfat and loin muscle) between Yorkshire and Korean Native Pigs (KNP) were collected, in particular those genes located in the internal between markers SW1355 and SW1823 where a quantitative trait loci (QTL) affecting the intramuscular fat content (IMF) have been detected in multiple pig populations. The genes listed here may offer information for further study the candidate genes affecting these QTL on the expression level.     

Association of two porcine reproductive and respiratory syndrome virus (PRRSV) receptor genes,CD163 and SN with immune traits

CD163 and sialoadhesin (SN) were reported as two essential receptors for the porcine reproductive and respiratory syndrome virus. To investigate the relationship between these two genes and porcine immunity, we assigned porcine CD163 and SN respectively to SSC5q21-q24 and SSC17q23 by IMpRH. Expression profiles revealed that CD163 and SN were ubiquitously expressed in ten tissues, and were expressed highly in lymph gland, spleen and liver, which implied the potential functions of CD163 and SN in immunity. Moreover, a single nucleotide polymorphism (SNP) c.3534C>T was found in 3′-UTR of the CD163 gene and association analysis showed that this gene was significantly associated with the IgG content in blood (P < 0.05). A novel missense mutation c.878A>G located in exon4 of the SN gene which caused the amino acid transition from histidine to arginine was detected, and it was significantly associated with the WBC count in the peripheral blood (P < 0.05). These results provided fundamental evidence for CD163 and SN as two functional candidate genes affecting immunity in pigs.     

过表达PDI对毕赤酵母分泌淀粉样蛋白的影响——以胱抑素为例

[目的] 本研究旨在结合酵母菌蛋白质二硫键异构酶（protein disulfide isomerase,PDI）与其底物蛋白鸡胱抑素C （chicken cystatin C,cC）在酵母中的共表达,理解PDI影响外源蛋白合成与表达的调控规律。运用转录组深度测序技术（RNA-Seq）筛选差异基因,调取并鉴定影响cC表达的关键基因,为解析外源蛋白高效表达机制,改造工程菌株提供理论支撑。[方法] 以巴斯德毕赤酵母GS115、GS115-cC为出发菌株,采用电转的方法将携带PDI编码基因的载体pPIC3.5K转入到GS115/GS115-cC菌株,使其在菌株中过表达,研究过表达PDI对cC表达的影响。采用RNA-Seq深度测序方法,研究重组毕赤酵母基因表达差异情况。并结合KEGG注释结果对数据进行分析,挑选差异显著表达基因进行验证,初步明确其在蛋白表达调控方面的功能。[结果] 本研究通过构建过表达PDI重组毕赤酵母菌株,使得外源蛋白cC的表达量显著增加。利用RNA-seq技术分析过表达PDI菌株与正常菌株的差异,最终筛选了373个差异表达基因,其中有122个差异基因注释到KEGG生物通路,包括12个基因注释到蛋白质转运和分解代谢途径,21个基因注释到蛋白质折叠分选和降解途径,以及24个基因参与蛋白质的翻译途径等。[结论] 在毕赤酵母中过表达PDI能显著增加外源蛋白cC的表达量。通过对过表达与正常表达PDI的毕赤酵母基因的表达谱分析,初步确定了其中一些转录情况变化显著的基因,明确了它们参与的细胞途径和信号通路,为改造具有高效率表达淀粉样蛋白的酵母菌株奠定基础。     

Molecular cloning of WIF1 and HMGA2 reveals ear‐preferential expression while uncovering a missense mutation associated with porcine ear size in WIF1

Considerable diversity exists in porcine ear size, which is an important morphological feature of pig breeds. Previously, we localized four crucial candidate genes—high mobility group AT‐hook 2 (HMGA2), LEM domain‐containing 3 (LEMD3), methionine sulfoxide reductase B3 (MSRB3) and Wnt inhibitory factor 1 (WIF1)—on Sus Scrofa chromosome 5 affecting porcine ear size, then cloned LEMD3 and MSBR3. In this study, we performed rapid amplification of cDNA ends to obtain full‐length cDNA sequences of 2338‐bp WIF1 and 2998‐bp HMGA2. Using quantitative real‐time PCR, we revealed that WIF1 expression was highest in ear cartilage of 60‐day‐old pigs and that this is therefore a better candidate gene for ear size than HMGA2. We further screened coding sequence variants in both genes and identified only one missense mutation (WIF1:c.1167C>G) in a conserved epidermal growth factor‐like domain from the mammalian WIF1 protein. The protein‐altering mutation was significantly associated with ear size across the Large White × Minzhu hybrid and Beijing Black pig populations. When WIF1:c.1167C>G was included as fixed effect in the model to re‐run a genome‐wide association study in the Large White × Minzhu intercross population the P‐value of the peak SNP on SSC5 from re‐running the genome‐wide association study dropped from 2.45E‐12 to 7.33E‐05. Taken together, the WIF1:c.1167C>G could be an important mutation associated with ear size. Our findings provide helpful information for further studies of the molecular mechanisms controlling porcine ear size.     

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₂ 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₂ 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.     

Global Analysis of Differentially Expressed Genes and Proteins in the Wheat Callus Infected by Agrobacterium tumefaciens

Agrobacterium-mediated plant transformation is an extremely complex and evolved process involving genetic determinants of both the bacteria and the host plant cells. However, the mechanism of the determinants remains obscure, especially in some cereal crops such as wheat, which is recalcitrant for Agrobacterium-mediated transformation. In this study, differentially expressed genes (DEGs) and differentially expressed proteins (DEPs) were analyzed in wheat callus cells co-cultured with Agrobacterium by using RNA sequencing (RNA-seq) and two-dimensional electrophoresis (2-DE) in conjunction with mass spectrometry (MS). A set of 4,889 DEGs and 90 DEPs were identified, respectively. Most of them are related to metabolism, chromatin assembly or disassembly and immune defense. After comparative analysis, 24 of the 90 DEPs were detected in RNA-seq and proteomics datasets simultaneously. In addition, real-time RT-PCR experiments were performed to check the differential expression of the 24 genes, and the results were consistent with the RNA-seq data. According to gene ontology (GO) analysis, we found that a big part of these differentially expressed genes were related to the process of stress or immunity response. Several putative determinants and candidate effectors responsive to Agrobacterium mediated transformation of wheat cells were discussed. We speculate that some of these genes are possibly related to Agrobacterium infection. Our results will help to understand the interaction between Agrobacterium and host cells, and may facilitate developing efficient transformation strategies in cereal crops.     

Genome-Wide Association Study Singles Out SCD and LEPR as the Two Main Loci Influencing Intramuscular Fat Content and Fatty Acid Composition in Duroc Pigs

Intramuscular fat (IMF) content and fatty acid composition affect the organoleptic quality and nutritional value of pork. A genome-wide association study was performed on 138 Duroc pigs genotyped with a 60k SNP chip to detect biologically relevant genomic variants influencing fat content and composition. Despite the limited sample size, the genome-wide association study was powerful enough to detect the association between fatty acid composition and a known haplotypic variant in SCD (SSC14) and to reveal an association of IMF and fatty acid composition in the LEPR region (SSC6). The association of LEPR was later validated with an independent set of 853 pigs using a candidate quantitative trait nucleotide. The SCD gene is responsible for the biosynthesis of oleic acid (C18:1) from stearic acid. This locus affected the stearic to oleic desaturation index (C18:1/C18:0), C18:1, and saturated (SFA) and monounsaturated (MUFA) fatty acids content. These effects were consistently detected in gluteus medius, longissimus dorsi, and subcutaneous fat. The association of LEPR with fatty acid composition was detected only in muscle and was, at least in part, a consequence of its effect on IMF content, with increased IMF resulting in more SFA, less polyunsaturated fatty acids (PUFA), and greater SFA/PUFA ratio. Marker substitution effects estimated with a subset of 65 animals were used to predict the genomic estimated breeding values of 70 animals born 7 years later. Although predictions with the whole SNP chip information were in relatively high correlation with observed SFA, MUFA, and C18:1/C18:0 (0.48–0.60), IMF content and composition were in general better predicted by using only SNPs at the SCD and LEPR loci, in which case the correlation between predicted and observed values was in the range of 0.36 to 0.54 for all traits. Results indicate that markers in the SCD and LEPR genes can be useful to select for optimum fatty acid profiles of pork.     

SNPs associated with body weight and backfat thickness in two pig breeds identified by a genome-wide association study

Growth and fat deposition are important economic traits due to the influence on production in pigs. In this study, a dataset of 1200 pigs with 345,570 SNPs genotyped by sequencing (GBS) was used to conduct a GWAS with single-marker regression method to identify SNPs associated with body weight and backfat thickness (BFT) and to search for candidate genes in Landrace and Yorkshire pigs. A total of 27 and 13 significant SNPs were associated with body weight and BFT, respectively. In the region of 149.85–149.89 Mb on SSC6, the SNP (SSC6: 149876737) for body weight and the SNP (SSC6: 149876507) for BFT were in the same locus region (a gap of 230 bp). Two SNPs were located in the DOCK7 gene, which is a protein-coding gene that plays an important role in pigmentation. Two SNPs located on SSC8: 54567459 and SSC11: 33043081 were found to overlap weight and BFT; however, no candidate gene was found in these regions. In addition, based on other significant SNPs, two positional candidate genes, NSRP1 and CADPS, were proposed to influence weight. In conclusion, this is the first study report using GBS data to identify the significant SNPs for weight and BFT. A total of four particularly interesting SNPs and one potential candidate genes (DOCK7) were found for these traits in domestic pigs. This study improves our knowledge to better understand the complex genetic architecture of weight and BFT, but further validation studies of these candidate loci and genes are recommended in pigs     

Molecular characterization of the porcine <Emphasis Type="Italic">JHDM1A</Emphasis> gene associated with average daily gain: evaluation its role in skeletal muscle development and growth

JHDM1A, a member of the JHDM (JmjC-domain-containing histone demethylase) family, plays an central role in gene silencing, cell cycle, cell growth and cancer development through histone H3K36 demethylation modification. Here reported the cloning, expression, chromosomal location and association analysis with growth traits of porcine JHDM1A gene. Sequence analysis showed that the porcine JHDM1A gene encodes 1,162 amino acids and contains JmjC, F-box, and CXXC zinc-finger domains, which coding sequence and deduced protein shares 91 and 99% similarity with human JHDM1A, respectively. Spatio-Temporal expression analysis indicated that the mRNA expression of porcine JHDM1A had significantly higher levels in the middle (65 days) and later (90 days) period’s embryo skeletal muscle than that of 33 days, and showed a ubiquitously expression but with the highest abundance in kidney, lung and liver of an adult pig. Radiation hybrid mapping and the following linkage mapping data indicate that JHDM1A maps to 2p17 region of pig chromosome 2 (SSC2). Allele frequency differences were detected in different pig breeds and an association study was performed with a SNP within 3′UTR. The results showed that there is a tendency for allele frequencies to differ between the fast growth breeds (Yorkshire) and slow growth pig breeds (Qingping pigs, Yushan Black pigs, Erhualian pigs and Dahuabai pigs). The association analysis using a Berkshire × Yorkshire F₂ population indicated that the C224G polymorphism had a highly significant association with average daily gain on test (P < 0.01), a trend association with average back fat thickness (P < 0.07), and significant associations (P < 0.01) on percent of average drip loss, Fiber Type II Ratio, muscle shear force and average lactate content in μmol/g. This study provides the first evidence that JHDM1A is differentially expressed in porcine embryonic skeletal muscle and associated with meat growth and quality traits.     

The chromosomal localization,expression pattern and polymorphism analysis of porcine <Emphasis Type="Italic">FSCN1</Emphasis> gene differently expressed from LongSAGE library

Fascin homologue 1 (FSCN1) has established roles in cell adhesion, motility, and cell–cell interactions. Our LongSAGE analysis suggested that FSCN1 was potentially differentially expressed in prenatal pig skeletal muscle. We have cloned the genomic DNA and mRNA sequence of FSCN1 gene and mapped it to SSC3p16-p17. The FSCN1 gene was differently expressed during prenatal skeletal muscle development and exhibited different expression pattern between Tongcheng and Landrace pigs. In Tongcheng pigs, FSCN1 expression was similar at 33 and 65 days post conception (dpc), and then sharply increased to a peak at 90 dpc. In Landrace pigs, however, expression increased between 33 and 65 dpc, peaked at 65 dpc, and was down-regulated thereafter. Significantly different expression levels between Tongcheng and Landrace were observed at 65 and 90 dpc. In postnatal pigs, it was strongly expressed only in the brain, but weakly in skeletal muscle and other tissues. We initially identified 32 SNPs through genomic DNA of FSCN1 gene. Association analysis suggested that the 6840^C/T mutation was significantly associated with the age at market weight (AGE) (p = 0.0004), average day gain from birth to market (ADG1) (p = 0.0002), and average day gain at testing period (ADG2) (p < 0.0001). Our study suggested that FSCN1 gene plays an in prenatal skeletal muscle development and was a candidate gene for meat production trait.     

Genome-wide association studies for iris pigmentation and heterochromia patterns in Large White pigs

Eye colour genetics have been extensively studied in humans since the rediscovery of Mendel’s laws. This trait was first interpreted using simplistic genetic models but soon it was realised that it is more complex. In this study, we analysed eye colour variability in a Large White pig population (n = 897) and report the results of GWASs based on several comparisons including pigs having four main eye colour categories (three with both pigmented eyes of different brown grades: pale, 17.9%; medium, 14.8%; and dark, 54.3%; another one with both eyes completely depigmented, 3.8%) and heterochromia patterns (heterochromia iridis – depigmented iris sectors in pigmented irises, 3.2%; heterochromia iridum – one whole eye iris of depigmented phenotype and the other eye with the iris completely pigmented, 5.9%). Pigs were genotyped with the Illumina PorcineSNP60 BeadChip and GEMMA was used for the association analyses. The results indicated that SLC45A2 (on chromosome 16, SSC16), EDNRB (SSC11) and KITLG (SSC5) affect the different grades of brown pigmentation of the eyes, the bilateral eye depigmentation defect and the heterochromia iridis defect recorded in this white pig population respectively. These genes are involved in several mechanisms affecting pigmentation. Significant associations for the eye depigmented patterns were also identified for SNPs on two SSC4 regions (including two candidate genes: NOTCH2 and PREX2) and on SSC6, SSC8 and SSC14 (including COL17A1 as candidate gene). This study provided useful information to understand eye pigmentation mechanisms, further valuing the pig as animal model to study complex phenotypes in humans.     

Putative Regulatory Factors Associated with Intramuscular Fat Content

Intramuscular fat (IMF) content is related to insulin resistance, which is an important prediction factor for disorders, such as cardiovascular disease, obesity and type 2 diabetes in human. At the same time, it is an economically important trait, which influences the sensorial and nutritional value of meat. The deposition of IMF is influenced by many factors such as sex, age, nutrition, and genetics. In this study Nellore steers (Bos taurus indicus subspecies) were used to better understand the molecular mechanisms involved in IMF content. This was accomplished by identifying differentially expressed genes (DEG), biological pathways and putative regulatory factors. Animals included in this study had extreme genomic estimated breeding value (GEBV) for IMF. RNA-seq analysis, gene set enrichment analysis (GSEA) and co-expression network methods, such as partial correlation coefficient with information theory (PCIT), regulatory impact factor (RIF) and phenotypic impact factor (PIF) were utilized to better understand intramuscular adipogenesis. A total of 16,101 genes were analyzed in both groups (high (H) and low (L) GEBV) and 77 DEG (FDR 10%) were identified between the two groups. Pathway Studio software identified 13 significantly over-represented pathways, functional classes and small molecule signaling pathways within the DEG list. PCIT analyses identified genes with a difference in the number of gene-gene correlations between H and L group and detected putative regulatory factors involved in IMF content. Candidate genes identified by PCIT include: ANKRD26, HOXC5 and PPAPDC2. RIF and PIF analyses identified several candidate genes: GLI2 and IGF2 (RIF1), MPC1 and UBL5 (RIF2) and a host of small RNAs, including miR-1281 (PIF). These findings contribute to a better understanding of the molecular mechanisms that underlie fat content and energy balance in muscle and provide important information for the production of healthier beef for human consumption.