A whole genome sequence association study for puberty in a large Duroc × Erhualian F2 population

Genomic prediction and genome-wide association studies for additive and dominance effects for body composition traits using 50 K and imputed high-density SNP genotypes in Yunong-black pigs

Body composition traits are complex traits controlled by minor genes and, in hybrid populations, are impacted by additive and nonadditive effects. We aimed to identify candidate genes and increase the accuracy of genomic prediction of body composition traits in crossbred pigs by including dominance genetic effects. Genomic selection (GS) and genome-wide association studies were performed on seven body composition traits in 807 Yunong-black pigs using additive genomic models (AM) and additive-dominance genomic models (ADM) with an imputed high-density single nucleotide polymorphism (SNP) array and the Illumina Porcine SNP50 BeadChip. The results revealed that the additive heritabilities estimated for AM and ADM using the 50 K SNP data ranged from 0.20 to 0.34 and 0.11 to 0.30, respectively. However, the ranges of additive heritability for AM and ADM in the imputed data ranged from 0.20 to 0.36 and 0.12 to 0.30, respectively. The dominance variance accounted for 23% and 27% of the total variance for the 50 K and imputed data, respectively. The accuracy of genomic prediction improved by 5% on average for 50 K and imputed data when dominance effect were considered. Without the dominance effect, the accuracies for 50 K and imputed data were 0.35 and 0.38, respectively, and 0.41 and 0.43, respectively, upon considering it. A total of 12 significant SNP and 16 genomic regions were identified in the AM, and 14 significant SNP and 21 genomic regions were identified in the ADM for both the 50 K and imputed data. There were five overlapping SNP in the 50 K and imputed data. In the AM, a significant SNP (CNC10041568) was found in both body length and backfat thickness traits, which was in the PLAG1 gene strongly and significantly associated with body length and backfat thickness in pigs. Moreover, a significant SNP (CNC10031356) with a heterozygous dominant genotype was present in the ADM. Furthermore, several functionally related genes were associated with body composition traits, including MOS, RPS20, LYN, TGS1, TMEM68, XKR4, SEMA4D and ARNT2. These findings provide insights into molecular markers and GS breeding for the Yunong-black pigs.

Genomic Scanning of Inbreeding Depression for Litter Size in Two Varieties of Iberian Pigs

Inbreeding depression is expected to be more pronounced in fitness-related traits, such as pig litter size. Recent studies have suggested that the genetic determinism of inbreeding depression may be heterogeneous across the genome. Therefore, the objective of this study was to conduct a genomic scan of the whole pig autosomal genome to detect the genomic regions that control inbreeding depression for litter size in two varieties of Iberian pigs (Entrepelado and Retinto). The datasets consisted of 2069 (338 sows) and 2028 (327 sows) records of litter size (Total Number Born and Number Born Alive) for the Entrepelado and Retinto varieties. All sows were genotyped using the Geneseek GGP PorcineHD 70 K chip. We employed the Unfavorable Haplotype Finder software to extract runs of homozygosity (ROHs) and conducted a mixed-model analysis to identify highly significant differences between homozygous and heterozygous sows for each specific ROH. A total of eight genomic regions located on SSC2, SSC5, SSC7, SSC8, and SSC13 were significantly associated with inbreeding depression, housing some relevant genes such as FSHR, LHCGR, CORIN, AQP6, and CEP120.

Genome-wide association studies for additive and dominance effects for body composition traits in commercial crossbred Piétrain pigs

Fat depth (FD) and muscle depth (MD) are economically important traits and used to estimate carcass lean content (LMP), which is one of the main breeding objectives in pig breeding programmes. We assessed the genetic architectures of body composition traits for additive and dominance effects in commercial crossbred Piétrain pigs using both 50 K array and sequence genotypes. We first performed a genome-wide association study (GWAS) using single-marker association analysis with a false discovery rate of 0.1. Then, we estimated the additive and dominance effects of the most significant variant in the quantitative trait loci (QTL) regions. It was investigated whether the use of whole-genome sequence (WGS) will improve the QTL detection (both additive and dominance) with a higher power compared with lower density SNP arrays. Our results showed that more QTL regions were detected by WGS compared with 50 K array (n = 54 vs. n = 17). Of the novel associated regions associated with FD and LMP and detected by WGS, the most pronounced peak was on SSC13, situated at ~116-118, 121-127 and 129-134 Mbp. Additionally, we found that only additive effects contributed to the genetic architecture of the analysed traits and no significant dominance effects were found for the tested SNPs at QTL regions, regardless of panel density. The associated SNPs are located in or near several relevant candidate genes. Of these genes, GABRR2, GALR1, RNGTT, CDH20 and MC4R have been previously reported as being associated with fat deposition traits. However, the genes on SSC1 (ZNF292, ORC3, CNR1, SRSF12, MDN1, TSHZ1, RELCH and RNF152) and SSC18 (TTC26 and KIAA1549) have not been reported previously to our best knowledge. Our current findings provide insights into the genomic regions influencing composition traits in Piétrain pigs.

Genome-Wide Association Studies, Runs of Homozygosity Analysis, and Copy Number Variation Detection to Identify Reproduction-Related Genes in Bama Xiang Pigs

Litter size and teat number are economically important traits in the porcine industry. However, the genetic mechanisms influencing these traits remain unknown. In this study, we analyzed the genetic basis of litter size and teat number in Bama Xiang pigs and evaluated the genomic inbreeding coefficients of this breed. We conducted a genome-wide association study to identify runs of homozygosity (ROH), and copy number variation (CNV) using the novel Illumina PorcineSNP50 BeadChip array in Bama Xiang pigs and annotated the related genes in significant single nucleotide polymorphisms and common copy number variation region (CCNVR). We calculated the ROH-based genomic inbreeding coefficients (FROH) and the Spearman coefficient between FROH and reproduction traits. We completed a mixed linear model association analysis to identify the effect of high-frequency copy number variation (HCNVR; over 5%) on Bama Xiang pig reproductive traits using TASSEL software. Across eight chromosomes, we identified 29 significant single nucleotide polymorphisms, and 12 genes were considered important candidates for litter-size traits based on their vital roles in sperm structure, spermatogenesis, sperm function, ovarian or follicular function, and male/female infertility. We identified 9,322 ROHs; the litter-size traits had a significant negative correlation to FROH. A total of 3,317 CNVs, 24 CCNVR, and 50 HCNVR were identified using cnvPartition and PennCNV. Eleven genes related to reproduction were identified in CCNVRs, including seven genes related to the testis and sperm function in CCNVR1 (chr1 from 311585283 to 315307620). Two candidate genes (NEURL1 and SH3PXD2A) related to reproduction traits were identified in HCNVR34. The result suggests that these genes may improve the litter size of Bama Xiang by marker-assisted selection. However, attention should be paid to deter inbreeding in Bama Xiang pigs to conserve their genetic diversity.

An imputed whole-genome sequence-based GWAS approach pinpoints causal mutations for complex traits in a specific swine population

Sequencing-based genome-wide association studies (GWAS) have facilitated the identification of causal associations between genetic variants and traits in diverse species. However, it is cost-prohibitive for the majority of research groups to sequence a large number of samples. Here, we carried out genotype imputation to increase the density of single nucleotide polymorphisms in a large-scale Swine F₂ population using a reference panel including 117 individuals, followed by a series of GWAS analyses. The imputation accuracies reached 0.89 and 0.86 for allelic concordance and correlation, respectively. A quantitative trait nucleotide (QTN) affecting the chest vertebrate was detected directly, while the investigation of another QTN affecting the residual glucose failed due to the presence of similar haplotypes carrying wild-type and mutant allelesin the reference panel used in this study. A high imputation accuracy was confirmed by Sanger sequencing technology for the most significant loci. Two candidate genes, CPNE5 and MYH3, affecting meat-related traits were proposed. Collectively, we illustrated four scenarios in imputation-based GWAS that may be encountered by researchers, and our results will provide an extensive reference for future genotype imputation-based GWAS analyses in the future.

The distribution of runs of homozygosity in the genome of river and swamp buffaloes reveals a history of adaptation, migration and crossbred events

Background

Water buffalo is one of the most important livestock species in the world. Two types of water buffalo exist: river buffalo (Bubalus bubalis bubalis) and swamp buffalo (Bubalus bubalis carabanensis). The buffalo genome has been recently sequenced, and thus a new 90 K single nucleotide polymorphism (SNP) bead chip has been developed. In this study, we investigated the genomic population structure and the level of inbreeding of 185 river and 153 swamp buffaloes using runs of homozygosity (ROH). Analyses were carried out jointly and separately for the two buffalo types.

Results

The SNP bead chip detected in swamp about one-third of the SNPs identified in the river type. In total, 18,116 ROH were detected in the combined data set (17,784 SNPs), and 16,251 of these were unique. ROH were present in both buffalo types mostly detected (~ 59%) in swamp buffalo. The number of ROH per animal was larger and genomic inbreeding was higher in swamp than river buffalo. In the separated datasets (46,891 and 17,690 SNPs for river and swamp type, respectively), 19,760 and 10,581 ROH were found in river and swamp, respectively. The genes that map to the ROH islands are associated with the adaptation to the environment, fitness traits and reproduction.

Conclusions

Analysis of ROH features in the genome of the two water buffalo types allowed their genomic characterization and highlighted differences between buffalo types and between breeds. A large ROH island on chromosome 2 was shared between river and swamp buffaloes and contained genes that are involved in environmental adaptation and reproduction.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00616-3.

Dynamic DNA methylation of ovaries during pubertal transition in gilts

Background

In female mammals, the initiation of puberty, coupling with the dramatically morphological changes in ovaries, indicates the sexual and follicular maturation. Previous studies have suggested that the disrupted DNA methylation results in the delayed puberty. However, to date, the changes in ovarian methylomes during pubertal transition have not been investigated. In this study, using gilts as a pubertal model, the genome-wide DNA methylation were profiled to explore their dynamics during pubertal transition across Pre-, In- and Post-puberty.

Results

During pubertal transition, the follicles underwent maturation and luteinization, coupled with the significant changes in the mRNA expression of DNMT1 and DNMT3a. DNA methylation levels of In-puberty were higher than that of Pre- and Post-puberty at the locations of genes and CpG islands (CGIs). Analysis of the DNA methylation changes identified 12,313, 20,960 and 17,694 differentially methylated CpGs (DMCs) for the comparisons of Pre- vs. In-, In vs. Post-, and Pre- vs. Post-puberty, respectively. Moreover, the CGIs, upstream and exonic regions showed a significant underrepresentation of DMCs, but the CGI shores, CGI shelves, intronic, downstream and intergenic regions showed a significant overrepresentation of DMCs. Furthermore, biological functions of these methylation changes enriched in PI3K-Akt signaling pathway, GnRH signaling pathway, and Insulin secretion, and the mRNA expressions of several genes of these signaling pathway, including MMP2, ESR1, GSK3B, FGF21, IGF1R, and TAC3, were significantly changed across Pre-, In- and Post-puberty in ovaries.

Conclusions

During pubertal transition in gilts, the DNA methylation changes of ovaries were likely to affect the transcription of genes related to PI3K-Akt signaling pathway, GnRH signaling pathway, and Insulin secretion. These observations can provide new insight into the epigenetic mechanism of follicular and sexual maturation during pubertal transition in mammals.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5884-x) contains supplementary material, which is available to authorized users.

Genome-Wide DNA Methylation Analysis of Hypothalamus During the Onset of Puberty in Gilts

Although selection of the early age at puberty in gilts will make for a favorable effect on the reproductivity of sow, a large proportion of phenotypic variation in age at puberty of gilts cannot be explained by genetics. Previous studies have implicated hypothalamic DNA methylation in the onset of puberty in mammals. However, the underlying molecular mechanism regarding the regulation of the onset of puberty has remained largely unexplored in gilts. Herein, the genome-scale DNA methylation of hypothalamus was acquired, using the reduced representation bisulfite sequencing, to compare and describe the changes of DNA methylation across Pre-, In- and Post-pubertal gilts. In this study, the average methylation levels of CpGs and CpHs (where H = C, T, or A) in CpG islands- and gene-related regions were gradually decreased in hypothalamic methylomes during the pubertal transition. Comparisons of Pre- vs. In-, In- vs. Post-, and Pre- vs. Post-pubertal stage revealed that there were 85726, 92914, and 100421 differentially methylated CpGs and 5940, 14804, and 16893 differentially methylated CpHs (where H = C, T, or A) in the hypothalamic methylomes. The methylation changes of CpHs were more dynamic than that of CpGs, and methylation changes of CpGs and CpHs were likely to be, respectively, involved in the developmental processes of reproduction and the molecular processes of cellular communications in the hypothalamus. Moreover, methylation changes of CpHs were observed to overrepresent in the quantitative trait loci of age at puberty, and the biological function of these CpH methylation changes was enriched in the pancreas development in gilts. Furthermore, the mRNA levels of several differentially CpG or CpH methylated genes related to the transcription of RNA II polymerase, GnRH signaling pathway, Estrogen signaling pathway, PI3K-AKt signaling pathway, and Insulin signaling pathway, including MAX, MMP2, FGF11, IGF1R, FGF21, and GSK3B, were significantly changed across these pubertal stages in the hypothalamus. These results will help our understanding of how DNA methylation contributes to phenotypic variation of age at puberty.

Cloning and expression of lin-28 homolog B gene in the onset of puberty in Duolang sheep

Objective

Recent studies have demonstrated that lin-28 homolog B (LIN28B)/miRNA let-7 (let-7) plays a role in the regulation of pubertal onset in mammals. However, the role of LIN28B/let-7 in the onset of ovine puberty remains unknown. We cloned the Duolang sheep Lin28B cDNA sequence, detected the expression change of LIN28B, let-7a and let-7g in hypothalamus, pituitary and ovary tissues at three different pubertal stages.

Methods

The reverse transcriptase polymerase chain reaction (RT-PCR) was used to clone the cDNA sequence of LIN28B gene from Duolang sheep and the bioinformatics methods were applied to analyze the amino acid sequence of LIN28B protein. The mRNA expression levels of the LIN28B gene at different pubertal stages were examined by real time RT-PCR.

Results

LIN28B cDNA of Duolang sheep was cloned, and two transcripts were obtained. The amino acid sequence of transcript 1 shares 99.60%, 98.78%, and 94.80% identity with those of goat, wild yak and pig, respectively. Strong LIN28B mRNA expression was detected in the hypothalamus, pituitary, ovary, oviduct and uterus, while moderate expression was found in the liver, kidney, spleen and heart, weak expression was observed in the heart. No expression was found in the lungs. Quantitative real-time PCR (QPCR) and western-blot analysis revealed that the LIN28B was highly expressed in the hypothalamus and ovary at prepuberty stages, and this expression significantly decreased from the prepuberty to puberty stages (p<0.05). Markedly increased levels of mRNA expression were detected in the pituitary from prepuberty to puberty (p<0.05) and then significantly decreased from puberty to postpuberty (p<0.05). The expression levels of let-7a and let-7g showed no significant changes among different pubertal stages (p>0.05).

Conclusion

These results provided a foundation for determining the functions of LIN28B/let-7 and their role in the onset of sheep puberty.

A whole-genome sequence based association study on pork eating quality traits and cooking loss in a specially designed heterogeneous F6 pig population

To determine the genetic basis of pork eating quality traits and cooking loss, we herein performed a genome-wide association study (GWAS) for tenderness, juiciness, oiliness, umami, overall liking and cooking loss by using whole genome sequences of heterogeneous stock F6 pigs which were generated by crossing 4 typical western pig breeds (Duroc, Landrace, Large White and Pietrain) and 4 typical Asian pig breeds (Erhualian, Laiwu, Bamaxiang and Tibetan). We identified 50 associated loci (QTLs) and most of them are novel. Seven loci also showed pleiotropic associations with different traits. In addition, we identified multiple promising candidate genes for these traits, including PAK1 and AQP11 for cooking loss, EP300 for tenderness, SDK1 for juiciness, FITM2 and 5-linked MYH genes for oiliness, and TNNI2 and TNNT3 for overall liking. Our results provide not only a better understanding of the genetic basis for meat quality, but also a potential application in future breeding for these complex traits.