An association study using imputed whole-genome sequence data identifies novel significant loci for growth-related traits in a Duroc × Erhualian F2 population

Glycyrrhizic Acid and Compound Probiotics Supplementation Alters the Intestinal Transcriptome and Microbiome of Weaned Piglets Exposed to Deoxynivalenol

Deoxynivalenol (DON) is a widespread mycotoxin that affects the intestinal health of animals and humans. In the present study, we performed RNA-sequencing and 16S rRNA sequencing in piglets after DON and glycyrrhizic acid and compound probiotics (GAP) supplementation to determine the changes in intestinal transcriptome and microbiota. Transcriptome results indicated that DON exposure altered intestinal gene expression involved in nutrient transport and metabolism. Genes related to lipid metabolism, such as PLIN1, PLIN4, ADIPOQ, and FABP4 in the intestine, were significantly decreased by DON exposure, while their expressions were significantly increased after GAP supplementation. KEGG enrichment analysis showed that GAP supplementation promoted intestinal digestion and absorption of proteins, fats, vitamins, and other nutrients. Results of gut microbiota composition showed that GAP supplementation significantly improved the diversity of gut microbiota. DON exposure significantly increased Proteobacteria, Actinobacteria, and Bacillus abundances and decreased Firmicutes, Lactobacillus, and Streptococcus abundances; however, dietary supplementation with GAP observably recovered their abundances to normal. In addition, predictive functions by PICRUSt analysis showed that DON exposure decreased lipid metabolism, whereas GAP supplementation increased immune system. This result demonstrated that dietary exposure to DON altered the intestinal gene expressions related to nutrient metabolism and induced disturbances of intestinal microbiota, while supplementing GAP to DON-contaminated diets could improve intestinal health for piglets.

GWAS of Reproductive Traits in Large White Pigs on Chip and Imputed Whole-Genome Sequencing Data

Total number born (TNB), number of stillborn (NSB), and gestation length (GL) are economically important traits in pig production, and disentangling the molecular mechanisms associated with traits can provide valuable insights into their genetic structure. Genotype imputation can be used as a practical tool to improve the marker density of single-nucleotide polymorphism (SNP) chips based on sequence data, thereby dramatically improving the power of genome-wide association studies (GWAS). In this study, we applied Beagle software to impute the 50 K chip data to the whole-genome sequencing (WGS) data with average imputation accuracy (R2) of 0.876. The target pigs, 2655 Large White pigs introduced from Canadian and French lines, were genotyped by a GeneSeek Porcine 50K chip. The 30 Large White reference pigs were the key ancestral individuals sequenced by whole-genome resequencing. To avoid population stratification, we identified genetic variants associated with reproductive traits by performing within-population GWAS and cross-population meta-analyses with data before and after imputation. Finally, several genes were detected and regarded as potential candidate genes for each of the traits: for the TNB trait: NOTCH2, KLF3, PLXDC2, NDUFV1, TLR10, CDC14A, EPC2, ORC4, ACVR2A, and GSC; for the NSB trait: NUB1, TGFBR3, ZDHHC14, FGF14, BAIAP2L1, EVI5, TAF1B, and BCAR3; for the GL trait: PPP2R2B, AMBP, MALRD1, HOXA11, and BICC1. In conclusion, expanding the size of the reference population and finding an optimal imputation strategy to ensure that more loci are obtained for GWAS under high imputation accuracy will contribute to the identification of causal mutations in pig breeding.

Genetic parameters and genome-wide association for milk production traits and somatic cell score in different lactation stages of Shanghai Holstein population

The aim of this study was to investigate the genetic parameters and genetic architectures of six milk production traits in the Shanghai Holstein population. The data used to estimate the genetic parameters consisted of 1,968,589 test-day records for 305,031 primiparous cows. Among the cows with phenotypes, 3,016 cows were genotyped with Illumina Bovine SNP50K BeadChip, GeneSeek Bovine 50K BeadChip, GeneSeek Bovine LD BeadChip v4, GeneSeek Bovine 150K BeadChip, or low-depth whole-genome sequencing. A genome-wide association study was performed to identify quantitative trait loci and genes associated with milk production traits in the Shanghai Holstein population using genotypes imputed to whole-genome sequences and both fixed and random model circulating probability unification and a mixed linear model with rMVP software. Estimated heritabilities (h2) varied from 0.04 to 0.14 for somatic cell score (SCS), 0.07 to 0.22 for fat percentage (FP), 0.09 to 0.27 for milk yield (MY), 0.06 to 0.23 for fat yield (FY), 0.09 to 0.26 for protein yield (PY), and 0.07 to 0.35 for protein percentage (PP), respectively. Within lactation, genetic correlations for SCS, FP, MY, FY, PY, and PP at different stages of lactation estimated in random regression model were ranged from -0.02 to 0.99, 0.18 to 0.99, 0.04 to 0.99, 0.04 to 0.99, 0.01 to 0.99, and 0.33 to 0.99, respectively. The genetic correlations were highest between adjacent DIM but decreased as DIM got further apart. Candidate genes included those related to production traits (DGAT1, MGST1, PTK2, and SCRIB), disease-related (LY6K, COL22A1, TECPR2, and PLCB1), heat stress–related (ITGA9, NDST4, TECPR2, and HSF1), and reproduction-related (7SK and DOCK2) genes. This study has shown that there are differences in the genetic mechanisms of milk production traits at different stages of lactation. Therefore, it is necessary to conduct research on milk production traits at different stages of lactation as different traits. Our results can also provide a theoretical basis for subsequent molecular breeding, especially for the novel genetic loci.

An imputation-based genome-wide association study for growth and fatness traits in Sujiang pigs

Sujiang pigs are a synthetic breed derived from Jiangquhai, Fengjing, and Duroc pigs. In this study, we sequenced the genome of 62 pigs with a coverage depth of 10× to 20×, including 27 Sujiang and 35 founder breed pigs, and we collected 360 global pigs' genome sequence data from public databases including 39 Duroc pigs. We obtained a high-quality variant dataset of 365 Sujiang pigs by imputing the porcine 80 K single nucleotide polymorphism (SNP) Beadchip to the whole-genome scale with a total of 422 pigs as a reference panel. A dataset of 365 imputated Sujiang pigs was used to perform single-trait genome-wide association study (GWAS) and meta-analyses for growth and fatness traits. Single-trait GWAS identified 1 907, 18, and 14 SNPs surpassing the suggestively significant threshold for backfat thickness, chest circumference, and chest width, respectively. Meta-analyses identified 2 400 genome-wide significant SNPs and 520 suggestively significant SNPs for backfat thickness and chest circumference, and 719 genome-wide significant SNPs and 1 225 suggestively significant SNPs for all seven traits. According to the meta-analysis of backfat thickness and chest circumference, a remarkable region of 2.69 Mb on Sus scrofa chromosome 4 containing FAM110B, IMPAD1, LYN, MOS, PENK, PLAG1, SDR16C5 and XKR4 was identified as a candidate region. The haplotype heat map of the 2.69 Mb region verified that Sujiang pigs were derived from Duroc and Chinese indigenous pigs, especially Jiangquhai pigs. The Kruskal-Wallis test showed that haplotypes of the 2.69 Mb region significantly affected backfat thickness and chest circumference traits. We then focused on PLAG1, an important growth-related gene, and identified two synonymous SNPs with obvious differences among different breeds in the PLAG1 gene. We then performed genotyping of 365 Sujiang, 150 Duroc, 95 Jiangquhai, and 100 Fengjing pigs to confirm the above result and verified that the two variants significantly affected phenotypes of growth and fatness traits. Our findings not only provide insights into the genetic architecture of porcine growth and fatness traits but also provide potential markers for selective breeding of these traits in Sujiang pigs.

Application of Genetic, Genomic and Biological Pathways in Improvement of Swine Feed Efficiency

Despite the significant improvement of feed efficiency (FE) in pigs over the past decades, feed costs remain a major challenge for producers profitability. Improving FE is a top priority for the global swine industry. A deeper understanding of the biology underlying FE is crucial for making progress in genetic improvement of FE traits. This review comprehensively discusses the topics related to the FE in pigs including: measurements, genetics, genomics, biological pathways and the advanced technologies and methods involved in FE improvement. We first provide an update of heritability for different FE indicators and then characterize the correlations of FE traits with other economically important traits. Moreover, we present the quantitative trait loci (QTL) and possible candidate genes associated with FE in pigs and outline the most important biological pathways related to the FE traits in pigs. Finally, we present possible ways to improve FE in swine including the implementation of genomic selection, new technologies for measuring the FE traits, and the potential use of genome editing and omics technologies.

An association study on imputed whole-genome resequencing from high-throughput sequencing data for body traits in crossbred pigs

Body traits are important economic factors in the pig industry. Genome-wide association studies (GWASs) have been widely applied using high-density genotype data to detect QTL in pigs. The aim of the present study was to detect the genetic variants significantly associated with body traits in crossbred pigs using the Illumina Porcine SNP50 BeadChip and imputed whole-genome sequence data. A set of seven body traits - body length, body height, chest circumference, cannon bone circumference, leg buttock circumference, back fat thickness and loin muscle depth - were measured. Moderate to high heritabilities were obtained for most traits (from 0.14 to 0.46), and significant genetic and phenotypic correlations among them were observed. GWAS identified 714 significantly associated SNPs located at 39 regions on all autosomes for body traits, and a total of seven functionally related candidate genes: PIK3CD, HOXA, PCGF2, CHST11, COL2A1, BMI1 and OSR2. Functional enrichment analysis revealed that candidate genes were enriched in the estrogen signaling pathway, embryonic skeletal system morphogenesis and embryonic skeletal system development. These results aim to uncover the genetic mechanisms underlying body development and marker-assisted selection programs focusing on body traits in pigs.

The construction of a haplotype reference panel using extremely low coverage whole genome sequences and its application in genome-wide association studies and genomic prediction in Duroc pigs

Extremely low coverage whole genome sequencing (lcWGS) is an economical technique to obtain high-density single nucleotide polymorphisms (SNPs). Here, we explored the feasibility of constructing a haplotype reference panel (lcHRP) using lcWGS and evaluated the effects of lcHRP through a genome-wide association study (GWAS) and genomic prediction in pigs. A total of 297 and 974 Duroc pigs were genotyped using lcWGS and a 50 K SNP array, respectively. We obtained 19,306,498 SNPs using lcWGS with an accuracy of 0.984. With the help of lcHRP, the accuracy of imputation from the SNP array to lcWGS was 0.922. Compared to the SNP array findings, those from the imputation-based GWAS identified more signals across four traits. With the integration of the top 1% imputation-based GWAS findings as genomic features, the accuracies of genomic prediction was improved by 6.0% to 13.2%. This study showed the great potential of lcWGS in pigs' molecular breeding.

Genome-Wide Association Study Based on Random Regression Model Reveals Candidate Genes Associated with Longitudinal Data in Chinese Simmental Beef Cattle

Simple Summary

Genome-wide association study (GWAS) has become the main approach for detecting functional genes that affects complex traits. For growth traits, the conventional GWAS method can only deal with the single-record traits observed at specific time points, rather than the longitudinal traits measured at multiple time points. Previous studies have reported the random regression model (RRM) for longitudinal data could overcome the limitation of the traditional GWAS model. Here, we present an association analysis based on RRM (GWAS-RRM) for 808 Chinese Simmental beef cattle at four stages of age. Ultimately, 37 significant single-nucleotide polymorphisms (SNPs) and several important candidate genes were screened to be associated with the body weight. Enrichment analysis showed these genes were significantly enriched in the signaling transduction pathway and lipid metabolism. This study not only offers a further understanding of the genetic basis for growth traits in beef cattle, but also provides a robust analytics tool for longitudinal traits in various species.

Abstract

Body weight (BW) is an important longitudinal trait that directly described the growth gain of bovine in production. However, previous genome-wide association study (GWAS) mainly focused on the single-record traits, with less attention paid to longitudinal traits. Compared with traditional GWAS models, the association studies based on the random regression model (GWAS-RRM) have better performance in the control of the false positive rate through considering time-stage effects. In this study, the BW trait data were collected from 808 Chinese Simmental beef cattle aged 0, 6, 12, and 18 months, then we performed a GWAS-RRM to fit the time-varied SNP effect. The results showed a total of 37 significant SNPs were associated with BW. Gene functional annotation and enrichment analysis indicated FGF4, ANGPT4, PLA2G4A, and ITGA5 were promising candidate genes for BW. Moreover, these genes were significantly enriched in the signaling transduction pathway and lipid metabolism. These findings will provide prior molecular information for bovine gene-based selection, as well as facilitate the extensive application of GWAS-RRM in domestic animals.

Accelerated deciphering of the genetic architecture of agricultural economic traits in pigs using a low-coverage whole-genome sequencing strategy

BACKGROUND

Uncovering the genetic architecture of economic traits in pigs is important for agricultural breeding. However, high-density haplotype reference panels are unavailable in most agricultural species, limiting accurate genotype imputation in large populations. Moreover, the infinitesimal model of quantitative traits implies that weak association signals tend to be spread across most of the genome, further complicating the genetic analysis. Hence, there is a need to develop new methods for sequencing large cohorts without large reference panels.

RESULTS

We describe a Tn5-based highly accurate, cost- and time-efficient, low-coverage sequencing method to obtain 11.3 million whole-genome single-nucleotide polymorphisms in 2,869 Duroc boars at a mean depth of 0.73×. On the basis of these single-nucleotide polymorphisms, a genome-wide association study was performed, resulting in 14 quantitative trait loci (QTLs) for 7 of 21 important agricultural traits in pigs. These QTLs harbour genes, such as ABCD4 for total teat number and HMGA1 for back fat thickness, and provided a starting point for further investigation. The inheritance models of the different traits varied greatly. Most follow the minor-polygene model, but this can be attributed to different reasons, such as the shaping of genetic architecture by artificial selection for this population and sufficiently interconnected minor gene regulatory networks.

CONCLUSIONS

Genome-wide association study results for 21 important agricultural traits identified 14 QTLs/genes and showed their genetic architectures, providing guidance for genetic improvement harnessing genomic features. The Tn5-based low-coverage sequencing method can be applied to large-scale genome studies for any species without a good reference panel and can be used for agricultural breeding.

Genome-wide association study for feed efficiency in collective cage-raised rabbits under full and restricted feeding

Feed efficiency (FE) is one of the most economically and environmentally relevant traits in the animal production sector. The objective of this study was to gain knowledge about the genetic control of FE in rabbits. To this end, GWASs were conducted for individual growth under two feeding regimes (full feeding and restricted) and FE traits collected from cage groups, using 114 604 autosome SNPs segregating in 438 rabbits. Two different models were implemented: (1) an animal model with a linear regression on each SNP allele for growth trait; and (2) a two‐trait animal model, jointly fitting the performance trait and each SNP allele content, for FE traits. This last modeling strategy is a new tool applied to GWAS and allows information to be considered from non‐genotyped individuals whose contribution is relevant in the group average traits. A total of 189 SNPs in 17 chromosomal regions were declared to be significantly associated with any of the five analyzed traits at a chromosome‐wide level. In 12 of these regions, 20 candidate genes were proposed to explain the variation of the analyzed traits, including genes such as FTO, NDUFAF6 and CEBPA previously associated with growth and FE traits in monogastric species. Candidate genes associated with behavioral patterns were also identified. Overall, our results can be considered as the foundation for future functional research to unravel the actual causal mutations regulating growth and FE in rabbits.

A Genome-Wide Association Study on Feed Efficiency Related Traits in Landrace Pigs

Feed efficiency (FE) traits in pigs are of utmost economic importance. Genetic improvement of FE related traits in pigs might significantly reduce production cost and energy consumption. Hence, our study aimed at identifying SNPs and candidate genes associated with FE related traits, including feed conversion ratio (FCR), average daily gain (ADG), average daily feed intake (ADFI), and residual feed intake (RFI). A genome-wide association study (GWAS) was performed for the four FE related traits in 296 Landrace pigs genotyped with PorcineSNP50 BeadChip. Two different single-trait methods, single SNP linear model GWAS (LM-GWAS) and single-step GWAS (ssGWAS), were implemented. Our results showed that the two methods showed high consistency with respect to SNP identification. A total of 32 common significant SNPs associated with the four FE related traits were identified. Bioinformatics analysis revealed eight common QTL regions, of which three QTL regions related to ADFI and RFI traits were overlapped. Gene ontology analysis revealed six common candidate genes (PRELID2, GPER1, PDX1, TEX2, PLCL2, ICAM2) relevant for the four FE related traits. These genes are involved in the processes of fat synthesis and decomposition, lipid transport process, insulin metabolism, among others. Our results provide, new insights into the genetic mechanisms and candidate function genes of FE related traits in pigs. However, further investigations to validate these results are warranted.

The SSC15 QTL-Rich Region Mutations Affecting Intramuscular Fat and Production Traits in Pigs

One of the more interesting regions in the pig genome is on chromosome 15 (115,800,000-122,100,000, SSC15, Sus scrofa 11.1) that has high quantitative trait locus (QTL) density associated with fattening, slaughter and meat quality characteristics. The SSC15 region encodes over 80 genes and a few miRNA sequences where potential genetic markers can be found. The goal of the study was to evaluate the effects of SSC15 mutations associated with villin 1 (VIL1), tensin 1 (TNS1), obscurin-like 1 (OBSL1) genes and with one long non-coding RNA (lncRNA) on productive pig traits and to enrich the genetic marker pool in further selection purpose. The potential genetic markers were identified using the targeted enrichment DNA sequencing (TEDNA-seq) of chromosome 15 region. The selected mutations were genotyped by using HRM, PCR and PCRRFLP methods. The association study was performed using the general linear model (GLM) in the sas program that included over 600 pigs of 5 Polish populations. The rs332253419 VIL1 mutation shows a significant effect on intramuscular fat (IMF) content in Duroc population where AA pigs had a 16% higher level than heterozygotes. The IMF content is also affected by the OBSL1 mutation, and the differences between groups are even up to 30%, but it is strongly dependent on breed factor. The OBSL1 mutation also significantly influences the meat yellowness, backfat thickness and pH level. The performed study delivers valuable information that could be highly useful during the development of the high-throughput genotyping method for further selection purposes in pigs. The OBSL1 and VIL1 mutations seem to be the most promising DNA marker showing a high effect on IMF level.