The number of vertebrae (NV), especially the number of thoracic vertebrae (NTV), varies among pig breeds. The NTV is controlled by vertebral segmentation and the number of somites during embryonic development. Although there is a high correlation between the NTV and NV, studies on a fixed NV have mainly considered the absolute numbers of thoracic vertebrae instead of vertebral segmentation. Therefore, this study aimed to discover variants associated with the NTV by considering the effect of the NV in pigs. The NTV and NV of 542 F2 individuals from a Large White × Minzhu pig crossbreed were recorded. All animals were genotyped for VRTN g.19034 A > C, LTBP2 c.4481A > C, and 37 missense or splice variants previously reported in a 951-kb interval on SSC7 and 147 single nucleotide polymorphisms (SNPs) on SSC14. To identify NTV-associated SNPs, we firstly performed a genome-wide association study (GWAS) using the Q + K (population structure + kinship matrix) model in TASSEL. With the NV as a covariate, the obtained data were used to identify the SNPs with the most significant genome-wide association with the NTV by performing a GWAS on a PorcineSNP60K Genotyping BeadChip. Finally, a conditional GWAS was performed by fixing this SNP. The GWAS showed that 31 SNPs on SSC7 have significant genome-wide associations with the NTV. No missense or splice variants were found to be associated with the NTV significantly. A linkage disequilibrium analysis suggested the existence of quantitative trait loci (QTL) in a 479-Kb region on SSC7, which contained a critical candidate gene FOS for the NTV in pigs. Subsequently, a conditional GWAS was performed by fixing M1GA0010658, the most significant of these SNPs. Two SNPs in BMPR1A were found to have significant genome-wide associations and a significant dominant effect. The leading SNP, S14_87859370, accounted for 3.86% of the phenotypic variance. Our study uncovered that regulation variants in FOS on SSC7 and in BMPR1A on SSC14 might play important roles in controlling the NTV, and thus these genetic factors may be harnessed for increasing the NTV in pigs.

Being in a confined environment causes chronic stress in gestating sows, which is detrimental for sow health, welfare and, consequently, offspring physiology. This study assessed the health and welfare of gestating sows housed in a group housing system compared to individual gestation stalls. After pregnancy was confirmed, experimental sows were divided randomly into two groups: the group housing system (GS), with the electronic sow feeding (ESF) system; or individual stall (IS). The behavior of sows housed in the GS or IS was then compared; throughout pregnancy, GS sows displayed more exploratory behavior, less vacuum chewing, and less sitting behavior (p < 0.05). IS sows showed higher stress hormone levels than GS sows. In particular, at 41 days of gestation, the concentration of the adrenocorticotropic hormone (ACTH) and adrenaline (A) in IS sows was significantly higher than that of GS sows, and the A level of IS sows remained significantly higher at 71 days of gestation (p < 0.01). The lipopolysaccharide (LPS) test was carried out in the weaned piglets of the studied sows. Compared with the offspring of gestating sows housed in GS (PG) or IS (PS), PG experienced a shorter period of high temperature and showed a quicker return to the normal state (p < 0.05). Additionally, their lower levels of stress hormone (p < 0.01) suggest that PG did not suffer from as much stress as PS. These findings suggested that gestating sows housed in GS were more able to carry out their natural behaviors and, therefore, had lower levels of stress and improved welfare. In addition, PG also showed better disease resistance and resilience. These results will provide a research basis for the welfare and breeding of gestating sows.

Background: Our previous genome-wide association study (GWAS) on milk fatty acid traits in Chinese Holstein cows revealed, the SNP, BTB-01556197, was significantly associated with C10:0 at genome-wide level (P = 0.0239). It was located in the down-stream of 5-hydroxytryptamine receptor 1B (HTR1B) gene that has been shown to play an important role in the regulation of fatty acid oxidation. Hence, we considered it as a promising candidate gene for milk fatty acids in dairy cattle. In this study, we aimed to investigate whether the HTR1B gene had significant genetic effects on milk fatty acid traits. Results: We re-sequenced the entire coding region and 3000 bp of 5′ and 3′ flanking regions of HTR1B gene. A total of 13 SNPs was identified, containing one in 5′ flanking region, two in 5′ untranslated region (UTR), two in exon 1, five in 3′ UTR, and three in 3′ flanking region. By performing genotype-phenotype association analysis with SAS9.2 software, we observed that 13 SNPs were significantly associated with medium-chain saturated fatty acids such as C6:0, C8:0 and C10:0 (P < 0.0001 ~ 0.042). With Haploview 4.1 software, linkage disequilibrium (LD) analysis was performed. Two haplotype blocks formed by two and ten SNPs were observed. Haplotype-based association analysis indicated that both haplotype blocks were strongly associated with C6:0, C8:0 and C10:0 as well (P < 0.0001 ~ 0.0071). With regards to the missense mutation in exon 1 (g.17303383G > T) that reduced amino acid change from alanine to serine, we predicted that it altered the secondary structure of HTR1B protein with SOPMA. In addition, we predicted that three SNPs in promoter region, g.17307103A > T, g.17305206 T > G and g.17303761C > T, altered the binding sites of transcription factors (TFs) HMX2, PAX2, FOXP1ES, MIZ1, CUX2, DREAM, and PPAR-RXR by Genomatix. Of them, luciferase assay experiment further confirmed that the allele T of g.17307103A > T significantly increased the transcriptional activity of HTR1B gene than allele A (P = 0.0007). Conclusions: In conclusion, our findings provided first evidence that the HTR1B gene had significant genetic effects on milk fatty acids in dairy cattle.

RNA-Seq technology is widely used to analyze global changes in the transcriptome and investigate the influence on relevant phenotypic traits. Beijing Black pigs show differences in growth rate and meat quality compared to western pig breeds. However, the molecular mechanisms responsible for such phenotypic differences remain unknown. In this study, longissimus dorsi muscles from Beijing Black and Yorkshire pigs were used to construct RNA libraries and perform RNA-seq. Significantly different expressions were observed in 1051 mRNAs, 322 lncRNAs, and 82 circRNAs. GO and KEGG pathway annotation showed that differentially expressed mRNAs participated in skeletal muscle development and fatty acid metabolism, which determined the muscle-related traits. To explore the regulatory role of lncRNAs, the cis and trans-target genes were predicted and these lncRNAswere involved in the biological processes related to skeletal muscle development and fatty acid metabolismvia their target genes. CircRNAs play a ceRNA role by binding to miRNAs. Therefore, the potential miRNAs of differentially expressed circRNAs were predicted and interaction networks among circRNAs, miRNAs, and key regulatory mRNAs were constructed to illustrate the function of circRNAs underlying skeletal muscle development and fatty acid metabolism. This study provides new clues for elucidating muscle phenotypic variation in pigs.

Porcine epidemic diarrhea virus (PEDV) causes porcine epidemic diarrhea (PED), a highly infectious disease, which has resulted in huge economic losses for the pig industry. To date, the pathogenic and immune response mechanism was not particularly clear. The purpose of this study was to investigate the pathogenic and immune responses of pigs infected with PEDV.In this study, 12 Min pigs were randomly selected without taking colostrum. At 3 days old, eight piglets were infected with 1 mL of PEDV solution (10 TCID50/ml), and the remaining four piglets were handled by 1 mL of 0.9% normal saline. Within the age of 7 days old, four piglets died and were considered as the death group. Correspondingly, four alive individuals were classified into the resistance group. Tissues of the duodenum, jejunum, ileum, colon, cecum, and rectum of piglets in the three groups were collected to measure the PEDV content. Additionally, the jejunum was used for the measurements and analyses of Hematoxylin-eosinstaining (HE), immunohistochemical sections, and transcriptomics. The phenotypes of Min piglets infected with PEDV showed that the viral copy numbers and jejunal damage had significant differences between the death and resistance groups. We also observed the transcriptome of the jejunum, and the differentially expressed (DE) analysis observed 6,585 DE protein-coding genes (PCGs), 3,188 DE long non-coding RNAs (lncRNAs), and 350 DE microRNAs (miRNAs), which were mainly involved in immune response and metabolic pathways. Furthermore, the specific expressed molecules for each group were identified, and 97 PCGs,108 lncRNAs, and 51 miRNAs were included in the ceRNA-regulated networks. By weighted gene co-expression network analysis (WGCNA) and transcription factor (TF) prediction, 27 significant modules and 32 significant motifs (E-value < 0.05) annotated with 519 TFs were detected. Of these TFs, 53 were DE PCGs. In summary, the promising key PCGs, lncRNAs, and miRNAs related to the pathogenic and immunological response of pigs infected with PEDV were detected and provided new insights into the pathogenesis of PEDV.

Sow milk production ability is an important limiting factor impacting suboptimal growth and the survival of piglets. Through pig genetic improvement, litter sizes have been increased. Larger litters need more suckling mammary glands, which results in increased milk from the lactating sow. Hence, there is much significance to exploring sow lactation performance. For milk production ability, it is not practical to directly measure the milk yield, we used litter weight gain (LWG) throughout sow lactation as an indicator. In this study, we estimated the heritability of LWG, namely, 0.18 ± 0.07. We then performed a GWAS, and detected seven significant SNPs, namely, Sus scrofa Chromosome (SSC) 2: ASGA0010040 (p = 7.73E-11); SSC2:MARC0029355 (p = 1.30E-08), SSC6: WU_10.2_6_65751151 (p = 1.32E-10), SSC7: MARC0058875 (p = 4.99E-09), SSC10: WU_10.2_10_49571394 (p = 6.79E-08), SSC11: M1GA0014659 (p = 1.19E-07), and SSC15: MARC0042106 (p = 1.16E-07). We performed the distribution of phenotypes corresponding to the genotypes of seven significant SNPs and showed that ASGA0010040, MARC0029355, MARC0058875, WU_10.2_10_49571394, M1GA0014659, and MARC0042106 had extreme phenotypic values that corresponded to the homozygous genotypes, while the intermediate values corresponded to the heterozygous genotypes. We screened for flanking regions ± 200 kb nearby the seven significant SNPs, and identified 38 genes in total. Among them, 28 of the candidates were involved in lactose metabolism, colostrum immunity, milk protein, and milk fat by functional enrichment analysis. Through the combined analysis between 28 candidate genes and transcriptome data of the sow mammary gland, we found nine commons (ANO3, MUC15, DISP3, FBXO6, CLCN6, HLA-DRA, SLA-DRB1, SLA-DQB1, and SLA-DQA1). Furthermore, by comparing the chromosome positions of the candidate genes with the quantitative trait locus (QTLs) as previously reported, a total of 17 genes were found to be within 0.86–94.02 Mb of the reported QTLs for sow milk production ability, in which, NAV2 was found to be located with 0.86 Mb of the QTL region ssc2: 40936355. In conclusion, we identified seven significant SNPs located on SSC2, 6, 7, 10, 11, and 15, and propose 28 candidate genes for the ability to produce milk in Yorkshire pigs, 10 of which were key candidates.