- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
Gestation length (GL) plays an important role in piglet maturation of major organs and development of body, while the genetic molecular markers of GL have not been extensively identified. In this study, according to the 5,662 effective records of 3,072 sows, the heritability and repeatability of GL were estimated through the dmuai of DMU Version 6.5.1 with a repeatability model, namely,
Introduction
Gestation length (GL) is the interval (days) of sows from insemination to farrowing, which belongs to the reproduction traits. The GL of pigs is well documented as 3 months, 3 weeks, and 3 days (114 days). GL was influenced by the genotypes of piglets and sow (Rydhmer et al., 2008). Within the last few days of gestation, the piglet maturation of major organs and development of body mass occurs (Ullrey et al., 1965), herein, increasing GL might be utilized to improve piglet development and survivability. George Charbonneau, a veterinarian, talked about “Keeping Large Litters Alive” and shared the attached two charts at the Shakespeare Swine Seminar (https://www.nationalhogfarmer.com/business/gestation-length-and-inducing-what-do-we-know), showing that sows allowed to farrow at day 118 had pigs born at an average weight of 3.7 pounds versus day 113 at 2.9 pounds (Supplementary Figure S1A) and indicating that earlier a sow farrows the lower the level of Immunoglobulin G colostrum in the milk for the pigs (Supplementary Figure S1B). The gain of GL was negatively correlated with litter size, while was positively correlated with the average pig weight (Omtvedt et al., 1965; Leigh, 1981). Selection for prolonged gestation might probably improve piglet survival after birth and piglet growth, while it might result in more stillbirths (Rydhmer et al., 2008). Hence, it is very important to understand the genetic mechanism of pregnancy days in pigs.
GL is not measurable until sexual maturity, and traditional breeding technologies are limited. Genomic selection improving the efficiency and accuracy of breeding is a form of marker-assisted selection in which genetic markers covering the entire genome are used (Yang et al., 2020). Some of the early reports proposed quantitative trait loci (QTL) and novel genetic variants associated with pig GL, while the limited number of DNA markers and F2 populations with large linkage disequilibrium blocks resulted in large genomic regions associated with GL (Wilkie et al., 1999; Chen et al., 2010). There were some potential QTLs, variants, and genes associated with GL of swine (Onteru et al., 2012; Hidalgo et al., 2016; Wang et al., 2018; See et al., 2019), while the relevant research was still lacking.
In this study, we obtained 5,662 valid GL records from 3,097 pigs, and adopted a strict repeatability model to estimate the heritability of GL. To identify the genetic SNPs and genes associated with GL, we conducted the GWAS analyses according to the GGP Porcine 50K Chip, imputed SNPs and performed the Fst analysis. Further, we carried out the functional enrichments to annotate these candidate genes. This research will provide valuable genetic resources for the future understanding of GL regulation mechanisms and breeding work.
Materials and methods
Animals and phenotype collection
A total of 3,072 pigs were used in the present study, and they were from two farms in Liaoning and Shanxi Provinces, China. The GL, breed, farm, parity, and total litter birth weight were collected for these sows, and 5,662 records were shown in Supplementary Table S1. The pedigrees of these sows were traced to four generations with 5,786 individuals to be used to estimate the genetic parameters of GL.
Genetic parameter estimation
By the dmuai of DMU Version 6.5.1, we estimated the variance and covariance components of GL based on a repeatability model, which was described as follows:
in which,
Further, we calculated the heritability and repeatability of GL by the following formulas:
where
Quality control of genotype data
Of all 3,072 sows with GL phenotypes, 906 individuals were used for GWAS. The ear samples of these 906 sows were collected and the DNA was isolated from them by a commercially available kit (Q1Aamp DNA Mini Kit, QIAGEN, Germany). Then these pigs were genotyped with the GeneSeek Genomic Profiler (GGP) Porcine 50K Chip (Illumina, San Diego, CA, United States), which contained 50,697 SNPs. With the PLINK, we performed the quality control by removing the SNPs with Minor Allele Frequencies (MAF) < 0.05, call rate <95%, and a deviation from Hardy-Weinberg equilibrium (HWE, p < 0.000001). Finally, we obtained 34,609 SNPs and 906 animals for GWAS, in which, all the individuals had a high SNP detection rate (>90%). All SNP positions were annotated according to the pig genome assembly Sscrofa 11.1.
Genotype imputation
To improve the identification of QTL, the whole-genome sequence data is expected to use (van den Berg et al., 2019), while sequencing thousands of individuals is expensive. In this study, we conducted the imputation for the GGP Porcine 50K Chip after quality control according to the online software PHARP v1 (http://alphaindex.zju.edu.cn/PHARP/index.php), which included 2,012 haplotypes, 34 million SNPs of autosomes, and 71 pig breeds (Wang et al., 2022). After imputation, we carried out the quality control (R2 ≥ 0.8 and MAF = 0.05), and obtained 9,212,179 SNPs, in which, the average R2 was 0.9248. Further, we conducted the quality control and GWAS according to the filter and model used in the front GGP Porcine 50K Chip, where 9,043,409 SNPs and 906 individuals were utilized.
Genome-wide association study
We conducted the GWAS by GCTA software with the model as follows:
where
The results of GWAS were shown by Manhattan and quantile-quantile (Q-Q) plots. A significant SNP was determined if the raw p-value <0.05/n, in which, n was the number of SNPs. Further, the threshold of suggestive significant association was calculated according to 1/n.
F-statistics analysis
In population genetics, Fst is an index to measure whether the actual frequency of genotypes in a population deviates from the theoretical proportion of genetic balance. In recent years, the number of litter sizes has increased, and the GL has been obvious for a long. To evaluate the degree of genetic differentiation of individuals with the different GL, Fst was carried out through the vcftools software. We calculated the Fst value of each SNP between GL 114 days and 115/116/117/118 days based on the following formulas:
and
In which,
The SNPs involved in the top 1% (Fst values) were selected as the core SNPs associated with the GL trait.
Candidate genes and functional annotation
By BioMart of Ensembl database, the candidate genes within 500 kb of significant SNPs were screened based on the pig reference genome (Sscrofa11.1). Further, we investigated the functions of these genes through the KOBAS (http://kobas.cbi.pku.edu.cn/kobas3/genelist/), and the significance threshold was corrected p-value < 0.05.
Results
Phenotypes and genetic parameters of gestation length
GL phenotype distribution and descriptive statistics were shown in Figure 1. The average GL was 115.97 days in this population, and the maximum and minimum values were 120 and 112, respectively. By estimation, the additive genetic variance V(a), permanent environmental effect variance V(p), and residual effect variance V(e) were detected, namely, 0.3014, 0.1594, and 1.4299, respectively. Correspondingly, the heritability and repeatability were calculated:
Identification of significant SNPs associated with GL by genome-wide association study
Through the GWAS based on the GGP Porcine 50K Chip, four genome-wide significant SNPs were identified (p < 1.44E-06 (0.05/34,609), Table 1 and Figures 2A,B), including WU_10.2_8_32052217 and WU_10.2_3_124292040 located on (Sus scrofa Chromosome (SSC) 3, and ALGA0073634 and WU_10.2_13_204954425 located on SSC13). According to the suggestive significant threshold p < 2.89E-05 (1/34,609), additional five SNPs were detected (Table 1 and Figures 2A,B), namely, WU_10.2_3_124371187 and WU_10.2_3_124478097 located on SSC3, H3GA0016351 and WU_10.2_5_111419379 located on SSC5, and CASI0009756 located on SSC13.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1046423/fgene-13-1046423-HTML-r1/image_m/fgene-13-1046423-g002.jpg)
FIGURE 2. Manhanttan and qq plots of GWAS for GL. The solid and dotted lines indicated significant and suggestive significant thresholds, respectively. (A,B). GWAS results based on the GGP Porcine 50K Chip. (C,D). GWAS results were based on the SNPs imputed.
Genome-wide association study based on the imputation genotype data
Further, the genotype was imputed by PHARP v1 to be used for a GWAS, and a total of 53 genome-wide significant SNPs [p < 5.52E-09 (0.05/9,043,409)] on pig chromosomes SSC13 were identified (Figures 2C,D; Table 2), including SNP ALGA0073634. Based on the suggestive significant threshold p < 1.11E-07 (1/9,043,409), 112 SNPs were detected to be associated with GL, containing 111 located on SSC13 and one located on SSC7 (Figures 2C,D; Table 2).
Detection of score SNPs associated with gestation length by F-statistics analysis
GLs of 114, 115, 116, 117, and 118 days were involved in 309, 263, 116, 117, and 49 individuals, respectively. Based on the GGP Porcine 50K Chip data, we calculated the Fst values of GL between 114 days and 115/116/117/118 days for each SNP, and the FST density distribution maps were constructed (Figure 3). The average Fst values of the four groups were 0.0041, 0.0279, 0.0341, and 0.0376, indicating that there was a differentiation between individuals with GL of 114 days and 115/116/117/118 days. In addition, the group of 114 and 118 days had the highest degree of differentiation than other groups.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1046423/fgene-13-1046423-HTML-r1/image_m/fgene-13-1046423-g003.jpg)
FIGURE 3. Distribution of Fst density on chromosome. (A) Analysis between 114 and 115 days. (B) Analysis between 114 and 116 days. (C) Analysis between 114 and 117 days. (D) Analysis between 114 and 118 days. The black lines indicated 320 SNPs in the top 1% (FST > 0.0317, 0.1489, 0.1748 and 0.2108 for 114 and 115/116/117/118 days, respectively).
For each group, 320 SNPs in the top 1% (FST > 0.0317, 0.1489, 0.1748, and 0.2108 for 114 and 115/116/117/118 days, respectively) were the core SNPs in the region where the signal occurred (Figure 3), a total of 885 core SNPs were identified (Supplementary Table S2). Out of them, SSC13 had the largest number of core SNPs with 119, followed by SSC7 (116 core SNPs) and SSC9 (89 core SNPs). The highest FST value (FST = 0.4035) was located on SSC1. According to the chromosomal location and FST value of the selection signal, the distribution of 885 core SNPs for each group was marked in Figure 4, which showed the differences in the distribution of selection signals related to GL across the genome. Combined with GWAS results, four SNPs, WU_10.2_8_32052217, WU_10.2_3_124292040, WU_10.2_3_124371187 and WU_10.2_3_124478097 located on SSC3, were co-identified by GWAS and Fst analysis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1046423/fgene-13-1046423-HTML-r1/image_m/fgene-13-1046423-g004.jpg)
FIGURE 4. The maps of 320 selection signatures on pig chromosome. (A) Analysis between 114 and 115 days. (B) Analysis between 114 and 116 days. (C) Analysis between 114 and 117 days. (D) Analysis between 114 and 118 days.
Screening candidate genes associated with gestation length and functional annotation
The genes located within a 500 kb distance on either side of the significant SNPs were considered candidates for GL. In total, 4,588 candidate genes were detected, of which, 27 were identified by GWAS of GGP Porcine 50K Chip data (Table 1), 12 were identified by GWAS of imputation SNPs (Table 2), and 4,574 were screened by core SNPs of Fst analysis (Supplementary Table S2). There were 20 genes simultaneously confirmed by GWAS and Fst analysis, namely HUNK, ARHGDIB, ERP27, RERG, NEDD9, TMEM170B, SCAF4, SOD1, TIAM1, ENSSSCG00000048838, ENSSSCG00000047227, EDN1, HIVEP1, ENSSSCG00000043944, LRATD1, ENSSSCG00000048577, ENSSSCG00000042932, ENSSSCG00000041405, ENSSSCG00000045589, and ADTRP.
To research the functions of candidate genes, we performed functional annotation by KOBAS. 12 genes from GWAS of GGP Porcine 50K Chip data were significantly involved in four KEGG pathways and 84 GO terms, including protein and lipid metabolism, regulation of multicellular organism growth, embryo implantation, and ovarian follicle development (Supplementary Table S3). Seven candidates from GWAS of imputation SNPs were strongly involved in 11 KEGG and 109 GO enrichments, which were mainly associated with body and organ signaling pathways, cellular metabolism, regulation of neuron apoptotic process, and heart development (Supplementary Table S3). 2,518 genes from Fst results heavily participated in 24 KEGG and 128 GO enrichments, containing cell differentiation, early endosome, and carbohydrate metabolic process etc. (Supplementary Table S3).
Discussion
GL is closely related to sow reproductive traits, such as litter size, herein deciphering the genetic mechanism of GL has important implications for genomic selective breeding. In the present study, we obtained the GL trait of pigs from a separate farm and estimated the heritability based on the 5,662 records of 3,072 sows. A total of 906 individuals were genotyped by 50K SNP Chip, and used for the GWAS by GCTA. Also, the genotypes were imputed to 9,212,179 SNPs for a further GWAS. In addition, we conducted the Fst analysis to detect more genetic markers for the GL trait. For these SNPs and candidate genes, we searched their enrichment annotation to provide preliminary functional analysis. While, the biological mechanism of these genetic associations between the SNPs or genes with GL still needs to be verified in future research.
The heritability of GL was estimated by a repeatability model with the pooling data of three breeds and was shown with
By the GWAS, nine and 112 SNPs significantly associated with GL were identified by GGP Porcine 50K Chip and imputed SNPs, respectively, and 34 candidate genes were screened. Among them, five genes located on SSC13 were identified at the same time by two association analyses. ENSSSCG00000039409 was a novel gene, and its function study was insufficient until now. KRTAP11-1 was a hair keratin-associated protein, as a possible crucial element for the physical properties of hair shafts, and was related to the embryonic hair development of Foxn1nu mice (Fujimoto et al., 2014). SCAF4 encoded the member of the serine/arginine-rich splicing factor family, which was an antiterminator during transcription (Gregersen et al., 2019). SCAF4 was important for correct usage of polyA sites for mRNA termination, and might be associated with the variable neurodevelopmental disorder (Fliedner et al., 2020). SOD1 was involved in fetal loss of women (Zhang et al., 2019) and muscle development (Abati et al., 2020). TIAM1 was part of the AMPK-Tiam1-Rac1 signaling pathway, which mediates contraction-stimulated glucose uptake in skeletal muscle cells and tissues (Yue et al., 2021). The lack of the Tiam1/Rac1 signal would contribute to the inhibition of embryo implantation in mice (Ma et al., 2015).
Fst is a measure of population divergence, and it is widespread to quantify the genetic distance between populations and to assess differentiation at individual SNPs by Fst (Bhatia et al., 2013). In this study, we observed the Fst values between 114 days and 115/116/117/118 days for each SNP, and the degree of differentiation was gradually increased from the groups of 114 and 115 days to 114 and 118 days, which meant the highest divergence between GL of 114 and 118 days. Combined with GWAS results, 20 genes were simultaneously identified through GWAS and Fst analysis, out of them, SCAF4, SOD1, and TIAM1 were detected in two GWASs, and ENSSSCG00000048838, ENSSSCG00000047227, ENSSSCG00000043944, ENSSSCG00000048577, ENSSSCG00000042932, ENSSSCG00000041405, and ENSSSCG00000045589 were the novel genes. HUNK was involved in transferase activity, transfer phosphorus-containing groups, and protein tyrosine kinase activity, and was thought to play a role in breast cancer metastasis (Dilday et al., 2020; Williams et al., 2020). ARHGDIB, RERG, and TMEM170B were also related to breast cancer (Li et al., 2018; Huang et al., 2020; Wang et al., 2020), implying they might participate in the development of mammary glands. ERP27 was a redox-inactive member of the PDI family, and enhanced receptor stability and downstream signaling maintain the stability of the organism (Kober et al., 2013). NEDD9 was found to regulate the developing embryonic nervous system (Knutson and Clagett-Dame, 2015), and HIVEP1 was involved in the obesity of newborns (Rizzo et al., 2020). ADTRP played a critical role in vascular development and vessel integrity and function (Lupu et al., 2021), and also specifies primitive myelopoiesis and definitive hematopoiesis (Luo et al., 2020). LRATD1 was the substrate protein of NMT1/2 (Su et al., 2021) and involved in cell morphogenesis and cell motility. EDN1 was linked to the local hemodynamic response of uteroplacental, and was a key player in the pathophysiology of hypoxia-induced fetal growth restriction (Thaete et al., 2007), and its signaling was crucial for the formation of the embryonic pharyngeal arches and their skeletal derivatives (Sasaki et al., 2013).
Conclusion
In the present study, the heritability 0.1594 of GL was shown, and a total of 1,002 SNPs were detected to impact GL phenotype by GWAS analyses of GGP Porcine 50K Chip and imputed SNPs, and Fst analysis. Further, 4,588 genes were considered candidates for GL, and their functions were confirmed, in which, 20 genes simultaneously confirmed by GWAS and Fst analysis were the most suggestive candidates.
Data availability statement
All raw data generated in this study were provided in this study and public repository (https://figshare.com/s/edda38a1c99aa7ab7ae0).
Ethics statement
The animal study was reviewed and approved by the Animal Ethics Committee of the Institute of Animal Science, Chinese Academy of Agricultural Sciences (Beijing, China).
Author contributions
LW conceived the study. LS and HL collected and analyzed the data. LS wrote the original draft. LW and LS performed the review and editing. All authors have read and agreed to the published version of the manuscript.
Funding
This research was supported by the funding of Agricultural Science and Technology Innovation Project (ASTIP-IAS02).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer XH declared a shared affiliation with the author(s) to the handling editor at the time of review.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.1046423/full#supplementary-material
References
Abati, E., Bresolin, N., Comi, G., and Corti, S. (2020). Silence superoxide dismutase 1 (SOD1): A promising therapeutic target for amyotrophic lateral sclerosis (ALS). Expert Opin. Ther. Targets 24 (4), 295–310. doi:10.1080/14728222.2020.1738390
Bhatia, G., Patterson, N., Sankararaman, S., and Price, A. L. (2013). Estimating and interpreting FST: The impact of rare variants. Genome Res. 23 (9), 1514–1521. doi:10.1101/gr.154831.113
Chen, C. Y., Guo, Y. M., Zhang, Z. Y., Ren, J., and Huang, L. S. (2010). A whole genome scan to detect quantitative trait loci for gestation length and sow maternal ability related traits in a White Duroc × Erhualian F2 resource population. Animal 4 (6), 861–866. doi:10.1017/S1751731110000169
Crump, R. E., Haley, C. S., Thompson, R., and Mercer, J. (1997). Individual animal model estimates of genetic parameters for reproduction traits of Landrace pigs performance tested in a commercial nucleus herd. Anim. Sci. 65, 285–290. doi:10.1017/s1357729800016593
Dilday, T., Ramos, N., and Yeh, E. (2020). HUNK signaling in metastatic breast cancer. Oncoscience 7 (5-6), 30–33. doi:10.18632/oncoscience.504
Fliedner, A., Kirchner, P., Wiesener, A., van de Beek, I., Waisfisz, Q., van Haelst, M., et al. (2020). Variants in SCAF4 cause a neurodevelopmental disorder and are associated with impaired mRNA processing. Am. J. Hum. Genet. 107 (3), 544–554. doi:10.1016/j.ajhg.2020.06.019
Fujimoto, S., Takase, T., Kadono, N., Maekubo, K., and Hirai, Y. (2014). Krtap11-1, a hair keratin-associated protein, as a possible crucial element for the physical properties of hair shafts. J. Dermatol. Sci. 74 (1), 39–47. doi:10.1016/j.jdermsci.2013.12.006
Gregersen, L. H., Mitter, R., Ugalde, A. P., Nojima, T., Proudfoot, N. J., Agami, R., et al. (2019). SCAF4 and SCAF8, mRNA anti-terminator proteins. Cell 177 (7), 1797–1813. doi:10.1016/j.cell.2019.04.038
Hanenberg, E. H. A. T., Knol, E. F., and Merks, J. W. M. (2001). Estimates of genetic parameters for reproduction traits at different parities in Dutch Landrace pigs. Livest. Prod. Sci. 69 (2), 179–186. doi:10.1016/s0301-6226(00)00258-x
Hidalgo, A. M., Lopes, M. S., Harlizius, B., and Bastiaansen, J. W. M. (2016). Genome-wide association study reveals regions associated with gestation length in two pig populations. Anim. Genet. 47 (2), 223–226. doi:10.1111/age.12395
Hollema, B. L., Zwiers, S., and Hermesch, S. (2020). Genetic parameters for haemoglobin levels in sows and piglets as well as sow reproductive performance and piglet survival. Animal 14 (4), 688–696. doi:10.1017/S1751731119002532
Huang, L., Tang, X., Shi, X., and Su, L. (2020). miR-532-5p promotes breast cancer proliferation and migration by targeting RERG. Exp. Ther. Med. 19 (1), 400–408. doi:10.3892/etm.2019.8186
Imboonte, N., and Kuhaaudomlarp, P. (2012). Genetic associations between stillbirth, total number of piglets born and gestation length in a commercial pig farm. Thai J. Vet. Med. 42 (2), 165–172.
Knutson, D. C., and Clagett-Dame, M. (2015). A complex RARE is required for the majority of Nedd9 embryonic expression. Transgenic Res. 24 (1), 123–134. doi:10.1007/s11248-014-9825-9
Kober, F. X., Koelmel, W., Kuper, J., Drechsler, J., Mais, C., Hermanns, H. M., et al. (2013). The crystal structure of the protein-disulfide isomerase family member ERp27 provides insights into its substrate binding capabilities. J. Biol. Chem. 288 (3), 2029–2039. doi:10.1074/jbc.M112.410522
Leigh, A. O. (1981). Factors affecting the gestation period of pigs in Nigeria. Trop. Anim. Health Prod. 13 (2), 87–93. doi:10.1007/BF02237898
Li, M., Han, Y., Zhou, H., Li, X., Lin, C., Zhang, E., et al. (2018). Transmembrane protein 170B is a novel breast tumorigenesis suppressor gene that inhibits the Wnt/β-catenin pathway. Cell Death Dis. 9 (2), 91. doi:10.1038/s41419-017-0128-y
Luo, C., Pook, E., Wang, F., Archacki, S. R., Tang, B., Zhang, W., et al. (2020). ADTRP regulates TFPI expression via transcription factor POU1F1 involved in coronary artery disease. Gene 753, 144805. doi:10.1016/j.gene.2020.144805
Lupu, C., Patel, M. M., and Lupu, F. (2021). Insights into the functional role of ADTRP (Androgen-Dependent TFPI-regulating protein) in health and disease. Int. J. Mol. Sci. 22 (9), 4451. doi:10.3390/ijms22094451
Ma, H-L., Gong, F., Tang, Y., Li, X., Li, X., Yang, X., et al. (2015). Inhibition of endometrial tiam1/rac1 signals induced by miR-22 up-regulation leads to the failure of embryo implantation during the implantation window in pregnant mice. Biol. Reprod. 92 (6), 152. doi:10.1095/biolreprod.115.128603
Nagy, I., Curik, I., Farkas, J., Csato, L., and Csornyei, Z. (2009). Bayesian inference of genetic parameters on litter size and gestation length in Hungarian Landrace and Hungarian Large White pigs. Ital. J. Anim. Sci. 8, 68–70. doi:10.4081/ijas.2009.s3.68
Ogawa, S., Konta, A., Kimata, M., Ishii, K., Uemoto, Y., and Satoh, M. J. A. S. J. (2019). Estimation of genetic parameters for farrowing traits in purebred Landrace and Large White pigs. Animal Sci. J. 90 (1), 23–28. doi:10.1111/asj.13120
Omtvedt, I. T., Stanislaw, C. M., and Whatley, J. A. (1965). Relationship of gestation length, age and weight at breeding, and gestation gain to sow productivity at farrowing. J. Anim. Sci. 24 (2), 531–535. doi:10.2527/jas1965.242531x
Onteru, S. K., Fan, B., Du, Z. Q., Garrick, D. J., Stalder, K. J., and Rothschild, M. F. (2012). A whole-genome association study for pig reproductive traits. Anim. Genet. 43 (1), 18–26. doi:10.1111/j.1365-2052.2011.02213.x
Rizzo, H. E., Escaname, E. N., Alana, N. B., Lavender, E., Gelfond, J., Fernandez, R., et al. (2020). Maternal diabetes and obesity influence the fetal epigenome in a largely Hispanic population. Clin. Epigenetics 12 (1), 34. doi:10.1186/s13148-020-0824-9
Rydhmer, L., Lundeheim, N., and Canario, L. (2008). Genetic correlations between gestation length, piglet survival and early growth. Livest. Sci. 115, 287–293. doi:10.1016/j.livsci.2007.08.014
Sasaki, M. M., Nichols, J. T., and Kimmel, C. B. (2013). edn1 and hand2 Interact in early regulation of pharyngeal arch outgrowth during zebrafish development. PLoS One 8 (6), e67522. doi:10.1371/journal.pone.0067522
See, G. M., Trenhaile-Grannemann, M. D., Spangler, M. L., Ciobanu, D. C., and Mote, B. E. (2019). A genome-wide association study for gestation length in swine. Anim. Genet. 50 (5), 539–542. doi:10.1111/age.12822
Su, D., Kosciuk, T., Yang, M., Price, I. R., and Lin, H. (2021). Binding affinity determines substrate specificity and enables discovery of substrates for N-myristoyltransferases. ACS Catal. 11 (24), 14877–14883. doi:10.1021/acscatal.1c03330
Thaete, L. G., Jilling, T., Synowiec, S., Khan, S., and Neerhof, M. G. (2007). Expression of endothelin 1 and its receptors in the hypoxic pregnant rat. Biol. Reprod. 77 (3), 526–532. doi:10.1095/biolreprod.107.061820
Ullrey, D. E., Sprague, J. I., Becker, D. E., and Miller, E. R. (1965). Growth of the swine fetus. J. Anim. Sci. 24, 711–717. doi:10.2527/jas1965.243711x
van den Berg, S., Vandenplas, J., van Eeuwijk, F. A., Bouwman, A. C., Lopes, M. S., and Veerkamp, R. F. (2019). Imputation to whole-genome sequence using multiple pig populations and its use in genome-wide association studies. Genet. Sel. Evol. 51, 2. doi:10.1186/s12711-019-0445-y
Wang, X., Bi, X., Huang, X., Wang, B., Guo, Q., and Wu, Z. (2020). Systematic investigation of biomarker-like role of ARHGDIB in breast cancer. Cancer Biomark. 28 (1), 101–110. doi:10.3233/CBM-190562
Wang, Y., Ding, X., Tan, Z., Xing, K., Yang, T., Pan, Y., et al. (2018). Genome-wide association study for reproductive traits in a Large White pig population. Anim. Genet. 49 (2), 127–131. doi:10.1111/age.12638
Wang, Z., Zhang, Z., Chen, Z., Sun, J., Cao, C., Wu, F., et al. (2022). Author correction: Pharp: A pig haplotype reference panel for genotype imputation. Sci. Rep. 12 (1), 13964. doi:10.1038/s41598-022-18078-y
Wilkie, P. J., Paszek, A. A., Beattie, C. W., Alexander, L. J., Wheeler, M. B., and Schook, L. B. (1999). A genomic scan of porcine reproductive traits reveals possible quantitative trait loci (QTLs) for number of corpora lutea. Mamm. Genome 10 (6), 573–578. doi:10.1007/s003359901047
Williams, C. B., Phelps-Polirer, K., Dingle, I. P., Williams, C. J., Rhett, M. J., Eblen, S. T., et al. (2020). HUNK phosphorylates EGFR to regulate breast cancer metastasis. Oncogene 39 (5), 1112–1124. doi:10.1038/s41388-019-1046-5
Yang, A. Q., Chen, B., Ran, M. L., Yang, G. M., and Zeng, C. (2020). The application of genomic selection in pig cross breeding. Yi Chuan 42 (2), 145–152. doi:10.16288/j.yczz.19-253
Yue, Y., Zhang, C., Zhao, X., Liu, S., Lv, X., Zhang, S., et al. (2021). Tiam1 mediates Rac1 activation and contraction-induced glucose uptake in skeletal muscle cells. Faseb J. 35 (2), e21210. doi:10.1096/fj.202001312R
Zhang, Y., Zhao, W., Xu, H., Hu, M., Guo, X., Jia, W., et al. (2019). Hyperandrogenism and insulin resistance-induced fetal loss: Evidence for placental mitochondrial abnormalities and elevated reactive oxygen species production in pregnant rats that mimic the clinical features of polycystic ovary syndrome. J. Physiol. 597 (15), 3927–3950. doi:10.1113/JP277879
Keywords: pig gestation length, genetic parameter estimation, GWAS, Fst, genetic molecular marker
Citation: Shi L, Li H and Wang L (2023) Genetic parameters estimation and genome molecular marker identification for gestation length in pigs. Front. Genet. 13:1046423. doi: 10.3389/fgene.2022.1046423
Received: 16 September 2022; Accepted: 28 November 2022;
Published: 05 January 2023.
Edited by:
Shahin Eghbalsaied, Islamic Azad University, IranReviewed by:
Xiaoyun He, Institute of Animal Sciences (CAAS), ChinaCong Li, Northwest A&F University, China
Guiyan Chu, Northwest A&F University, China
Copyright © 2023 Shi, Li and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lixian Wang, iaswlx@263.net
†These authors have contributed equally to this work