A genome-wide scan for quantitative trait loci affecting serum glucose and lipids in a White Duroc × Erhualian intercross F2 population

Genome-wide association study of 17 serum biochemical indicators in a chicken F2 resource population

BACKGROUND

Serum biochemical indicators are often regarded as direct reflections of animal metabolism and health. The molecular mechanisms underlying serum biochemical indicators metabolism of chicken (Gallus Gallus) have not been elucidated. Herein, we performed a genome-wide association study (GWAS) to identify the variation associated with serum biochemical indicators. The aim of this research was to broaden the understanding of the serum biochemical indicators in chickens.

RESULTS

A GWAS of serum biochemical indicators was carried out on 734 samples from an F2 Gushi× Anka chicken population. All chickens were genotyped by sequencing, 734 chickens and 321,314 variants were obtained after quality control. Based on these variants, a total of 236 single-nucleotide polymorphisms (SNPs) on 9 chicken chromosomes (GGAs) were identified to be significantly (-log10(P) > 5.72) associated with eight of seventeen serum biochemical indicators. Ten novel quantitative trait locis (QTLs) were identified for the 8 serum biochemical indicator traits of the F2 population. Literature mining revealed that the ALPL, BCHE, GGT2/GGT5 genes at loci GGA24, GGA9 and GGA15 might affect the alkaline phosphatase (AKP), cholinesterase (CHE) and γ-glutamyl transpeptidase (GGT) traits, respectively.

CONCLUSION

The findings of the present study may contribute to a better understanding of the molecular mechanisms of chicken serum biochemical indicator regulation and provide a theoretical basis for chicken breeding programs.

The Effect of BSCL2 Gene on Fat Deposition Traits in Pigs

BSCL2 encodes seipin, a transmembrane endoplasmic reticulum protein associated with lipodystrophy and severe metabolic complications, including diabetes and hepatic steatosis. In pigs, BSCL2 expression increases during adipocyte differentiation. In the present study, we identified significant gene variants associated with fat deposition (FD)-related processes based on subcutaneous fat tissue RNA-seq data. In the association study, to prove our hypothesis, three Polish pig breeds were included: Złotnicka White (ZW, n = 72), Polish Landrace (PL, n = 201), and Polish Large White (PLW, n = 169). Based on variant calling analysis and χ² tests, BSCL2 mutations showing significantly different genotype/allele distribution between high- and low-fat pigs were selected for a comprehensive association study. Four interesting BSCL2 variants (rs346079334, rs341493267, rs330154033, and rs81333153) belonging to downstream and missense mutations were investigated. Our study showed a significant decrease in minor allele frequency for two BSCL2 variants (rs346079334 and rs341493267) in PL pigs in 2020-2021. In ZW, BSCL2 mutations significantly affected loin and ham fats, meat redness, and growth performance traits, such as feed conversion and daily feed intake. Similar observations were noted for PLW and PL, where BSCL2 mutations influenced fat depositions and meat traits, such as loin eye area, loin mass and fat, carcass yield, and growth performance traits. Based on the observation in pigs, our study supports the theory that BSCL2 expressed in subcutaneous fat is involved in the FD process.

Genome-Wide Association Analysis for Blood Lipid Traits Measured in Three Pig Populations Reveals a Substantial Level of Genetic Heterogeneity

Serum lipids are associated with myocardial infarction and cardiovascular disease in humans. Here we dissected the genetic architecture of blood lipid traits by applying genome-wide association studies (GWAS) in 1,256 pigs from Laiwu, Erhualian and Duroc × (Landrace × Yorkshire) populations, and a meta-analysis of GWAS in more than 2,400 pigs from five diverse populations. A total of 22 genomic loci surpassing the suggestive significance level were detected on 11 pig chromosomes (SSC) for six blood lipid traits. Meta-analysis of GWAS identified 5 novel loci associated with blood lipid traits. Comparison of GWAS loci across the tested populations revealed a substantial level of genetic heterogeneity for porcine blood lipid levels. We further evaluated the causality of nine polymorphisms nearby or within the APOB gene on SSC3 for serum LDL-C and TC levels. Of the 9 polymorphisms, an indel showed the most significant association with LDL-C and TC in Laiwu pigs. But the significant association was not identified in the White Duroc × Erhualian F2 resource population, in which the QTL for LDL-C and TC was also detected on SSC3. This indicates that population-specific signals may exist for the SSC3 QTL. Further investigations are warranted to validate this assumption.

A genome-wide association analysis for porcine serum lipid traits reveals the existence of age-specific genetic determinants

Background

The genetic determinism of blood lipid concentrations, the main risk factor for atherosclerosis, is practically unknown in species other than human and mouse. Even in model organisms, little is known about how the genetic determinants of lipid traits are modulated by age-specific factors. To gain new insights into this issue, we have carried out a genome-wide association study (GWAS) for cholesterol (CHOL), triglyceride (TRIG) and low (LDL) and high (HDL) density lipoprotein concentrations measured in Duroc pigs at two time points (45 and 190 days).

Results

Analysis of data with mixed-model methods (EMMAX, GEMMA, GenABEL) and PLINK showed a low positional concordance between trait-associated regions (TARs) for serum lipids at 45 and 190 days. Besides, the proportion of phenotypic variance explained by SNPs at these two time points was also substantially different. The four analyses consistently detected two regions on SSC3 (124 Mb, CHOL and LDL at 190 days) and SSC6 (135 Mb, CHOL and TRIG at 190 days) with highly significant effects on the porcine blood lipid profile. Moreover, we have found that SNP variation within SSC3, SSC6, SSC10, SSC13 and SSC16 TARs is associated with the expression of several genes mapping to other chromosomes and related to lipid metabolism.

Conclusions

Our data demonstrate that the effects of genomic determinants influencing lipid concentrations in pigs, as well as the amount of phenotypic variance they explain, are influenced by age-related factors.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-758) contains supplementary material, which is available to authorized users.

Evaluation of the causality of the low-density lipoprotein receptor gene (LDLR) for serum lipids in pigs

A significant quantitative trait locus (QTL) for low-density lipoprotein cholesterol (LDL-C) and total cholesterol (TC) was identified around the LDLR gene on chromosome 2 (SSC2) in a White Duroc × Erhualian F2 resource population and Sutai pigs in our previous study. However, in previous reports, the causality of LDLR with serum lipids is controversial in pigs. To systematically assess the causality of LDLR with serum lipids, association analyses were successively performed in three populations: Sutai pigs, a White Duroc × Erhualian F2 resource population and a Duroc × (Landrace × Large White) population. We first performed a haplotype-based association study with 60K SNP genotyping data and evidenced the significant association with LDL-C and TC around the LDLR gene region. We also found that there is more than one QTL for LDL-C and TC on SSC2. Then, we evaluated the causalities of two missense mutations, c.1812C>T and c.1520A>G, with LDL-C and TC. We revealed that the c.1812C>T SNP showed the strongest association with LDL-C (P = 5.40 × 10(-11) ) and TC (P = 3.64 × 10(-8) ) and explained all the QTL effect in Sutai pigs. Haplotype analysis found that two missense SNPs locate within a 1.93-Mb haplotype block. One major haplotype showed the strongest significant association with LDL-C (P = 4.62 × 10(-18) ) and TC (P = 1.06 × 10(-9) ). However, the c.1812C>T SNP was not identified in the White Duroc × Erhualian intercross, and the association of c.1520A>G with both LDL-C and TC did not achieve significance in this F2 population, suggesting population heterogeneity. Both missense mutations were identified in the Duroc × (Landrace × Large White) population and showed significant associations with LDL-C and TC. Our data give evidence that the LDLR gene should be a candidate causative gene for LDL-C and TC in pigs, but heterogeneity exists in different populations.

Genetic dissection of blood lipid traits by integrating genome-wide association study and gene expression profiling in a porcine model

BACKGROUND

Serum concentrations of total cholesterol (TC), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C) and triglycerides (TG) are highly heritable traits that are used clinically to evaluate risk for cardiovascular disease in humans. In this study, we applied a genome-wide association study (GWAS) in 1,075 pigs from two populations and gene expression studies on 497 liver samples to dissect the genetic basis of serum lipids in a pig model.

RESULTS

We totally identified 8, 5, 2 and 3 genomic loci harboring 109 SNPs that were significantly associated with LDL-C, TC, TG and the ratio of HDL-C/LDL-C in two experimental populations, respectively. In the F2 population, the most prominent SNP was identified at the SSC3: 124,769,847 bp where APOB is the well-known candidate gene. However, in the Sutai population, the most number of significant SNPs was identified at SSC2: 64.97-82.22 Mb where LDLR was identified as the candidate gene. Furthermore, we firstly reported 4 novel genomic loci in pigs harboring the LDL-C-associated SNPs. We also observed obvious population heterogeneity in the two tested populations. Through whole-genome gene expression analysis, we detected 718 trait-correlated expressions. Many of these transcripts correspond to candidate genes for blood lipids in humans. The GWAS mapped 120 cis-eQTLs and 523 trans-eQTLs for these transcripts. One gene encoding the transcript gnl|UG|Ssc#S35330332 stands out to be an important candidate gene for LDL-C by an integrative analysis of GWAS, eQTL and trait-associated expression.

CONCLUSIONS

We identified the genomic regions or candidate genes associated with blood lipids by an integrative analysis of GWAS, QTT and eQTL mapping in pigs. The findings would benefit the further identification of the causative genes for blood lipid traits in both pigs and humans.

Association of the rs7395662 SNP in the MADD-FOLH1 and several environmental factors with serum lipid levels in the Mulao and Han populations

BACKGROUND

The rs7395662 single nucleotide polymorphism (SNP) in the MADD-FOLH1 has been associated with serum lipid traits, but the results are inconsistent in different populations. The present study was undertaken to investigate the association of rs7395662 SNP and several environmental factors with serum lipid levels in the Guangxi Mulao and Han populations.

METHOD

A total of 721 subjects of Mulao and 727 subjects of Han Chinese were randomly selected from our previous stratified randomized samples. Genotyping of the SNP was performed by polymerase chain reaction and restriction fragment length polymorphism combined with gel electrophoresis, and confirmed by direct sequencing.

RESULTS

Serum apolipoprotein (Apo) B levels were higher in Mulao than in Han (P < 0.01). The allelic and genotypic frequencies in Han were different between males and females (P < 0.05 for each), but there was no difference between Mulao and Han or between Mulao males and females. The levels of low-density lipoprotein cholesterol (LDL-C) and ApoB in Mulao females were different among the genotypes (P < 0.05), the G allele carriers had higher LDL-C and ApoB levels than the G allele non-carriers. The levels of total cholesterol (TC), triglyceride (TG), LDL-C and ApoB in Han males and TC, TG and high-density lipoprotein cholesterol (HDL-C) in Han females were different among the genotypes (P < 0.05-0.01), the subjects with GG genotype in Han males had higher TC, TG, and ApoB and lower LDL-C levels than the subjects with AA or AG genotype, and the G allele carriers in Han females had lower TC and HDL-C levels than the G allele non-carriers. The levels of LDL-C and ApoB in Mulao females were correlated with the genotypes (P < 0.05 for each). The levels of HDL-C and ApoAI in Han males and HDL-C in Han females were correlated with genotypes (P < 0.05-0.001). Serum lipid parameters were also correlated with several environmental factors in both ethnic groups (P < 0.05-0.01).

CONCLUSION

The association of rs7395662 SNP and serum lipid levels is different between the Mulao and Han populations, and between males and females in both ethnic groups.

A comprehensive survey of copy number variation in 18 diverse pig populations and identification of candidate copy number variable genes associated with complex traits

Background

Copy number variation (CNV) is a major source of structural variants and has been commonly identified in mammalian genome. It is associated with gene expression and may present a major genetic component of phenotypic diversity. Unlike many other mammalian genomes where CNVs have been well annotated, studies of porcine CNV in diverse breeds are still limited.

Result

Here we used Porcine SNP60 BeadChip and PennCNV algorithm to identify 1,315 putative CNVs belonging to 565 CNV regions (CNVRs) in 1,693 pigs from 18 diverse populations. Total 538 out of 683 CNVs identified in a White Duroc × Erhualian F₂ population fit Mendelian transmission and 6 out of 7 randomly selected CNVRs were confirmed by quantitative real time PCR. CNVRs were non-randomly distributed in the pig genome. Several CNV hotspots were found on pig chromosomes 6, 11, 13, 14 and 17. CNV numbers differ greatly among different pig populations. The Duroc pigs were identified to have the most number of CNVs per individual. Among 1,765 transcripts located within the CNVRs, 634 genes have been reported to be copy number variable genes in the human genome. By integrating analysis of QTL mapping, CNVRs and the description of phenotypes in knockout mice, we identified 7 copy number variable genes as candidate genes for phenotypes related to carcass length, backfat thickness, abdominal fat weight, length of scapular, intermuscle fat content of logissimus muscle, body weight at 240 day, glycolytic potential of logissimus muscle, mean corpuscular hemoglobin, mean corpuscular volume and humerus diameter.

Conclusion

We revealed the distribution of the unprecedented number of 565 CNVRs in pig genome and investigated copy number variable genes as the possible candidate genes for phenotypic traits. These findings give novel insights into porcine CNVs and provide resources to facilitate the identification of trait-related CNVs.

QTL analysis of clinical-chemical traits in an F₂ intercross between Landrace and Korean native pigs

Clinical-chemical traits are essential when examining the health status of individuals. The aim of this study was to identify quantitative trait loci (QTL) and the associated positional candidate genes affecting clinical-chemical traits in a reciprocal F(2) intercross between Landrace and Korean native pigs. Following an overnight fast, 25 serum phenotypes related to clinical-chemical traits (e.g., hepatic function parameters, renal function parameters, electrolyte, lipids) were measured in >970 F(2) progeny. All experimental samples were subjected to genotyping analysis using 165 microsatellite markers located across the genome. We identified eleven genome-wide significant QTL in six chromosomal regions (SSC 2, 7, 8, 13, 14, and 15) and 59 suggestive QTL in 17 chromosomal regions (SSC 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13, 14, 15, 16, 17, and 18). We also observed significant effects of reciprocal crosses on some of the traits, which would seem to result from maternal effect, QTL on sex chromosomes, imprinted genes, or genetic difference in mitochondrial DNA. The role of genomic imprinting in clinical-chemical traits also was investigated. Genome-wide analysis revealed a significant evidence for an imprinted QTL in SSC4 affecting serum amylase levels. Additionally, a series of bivariate linkage analysis provided strong evidence that QTL in SSC 2, 13, 15, and 18 have a pleiotropic effect on clinical-chemical traits. In conclusion, our study detected both novel and previously reported QTL influencing clinical-chemical traits in pigs. The identified QTL together with the positional candidate genes identified here could play an important role in elucidating the genetic structure of clinical-chemical phenotype variation in humans and swine.

Sex-specific association of rs16996148 SNP in the NCAN/CILP2/PBX4 and serum lipid levels in the Mulao and Han populations

BACKGROUND

The association of rs16996148 single nucleotide polymorphism (SNP) in NCAN/CILP2/PBX4 and serum lipid levels is inconsistent. Furthermore, little is known about the association of rs16996148 SNP and serum lipid levels in the Chinese population. We therefore aimed to detect the association of rs16996148 SNP and several environmental factors with serum lipid levels in the Guangxi Mulao and Han populations.

METHOD

A total of 712 subjects of Mulao nationality and 736 participants of Han nationality were randomly selected from our stratified randomized cluster samples. Genotyping of the rs16996148 SNP was performed by polymerase chain reaction and restriction fragment length polymorphism combined with gel electrophoresis, and then confirmed by direct sequencing.

RESULTS

The levels of apolipoprotein (Apo) B were higher in Mulao than in Han (P < 0.001). The frequencies of G and T alleles were 87.2% and 12.8% in Mulao, and 89.9% and 10.1% in Han (P <0.05); respectively. The frequencies of GG, GT and TT genotypes were 76.0%, 22.5% and 1.5% in Mulao, and 81.2%, 17.4% and 1.4% in Han (P <0.05); respectively. There were no significant differences in the genotypic and allelic frequencies between males and females in both ethnic groups. The levels of HDL-C, ApoAI, and the ratio of ApoAI to ApoB in Mulao were different between the GG and GT/TT genotypes in males but not in females (P < 0.01 for all), the subjects with GT/TT genotypes had higher serum levels of HDL-C, ApoAI, and the ratio of ApoAI to ApoB than the subjects with GG genotype. The levels of TC, TG, LDL-C, ApoAI, and ApoB in Han were different between the GG and GT/TT genotypes in males but not in females (P < 0.05-0.001), the T allele carriers had higher serum levels of TC, TG, LDL-C, ApoAI, and ApoB than the T allele noncarriers. The levels of HDL-C, ApoAI, and the ratio of ApoAI to ApoB in Mulao were correlated with the genotypes in males (P < 0.05-0.01) but not in females. The levels of TC, TG, HDL-C, LDL-C, ApoAI and ApoB in Han were associated with the genotypes in males (P < 0.05-0.001) but not in females. Serum lipid parameters were also correlated with several enviromental factors in both ethnic groups (P < 0.05-0.001).

CONCLUSIONS

The genotypic and allelic frequencies of rs16996148 SNP and the associations of the SNP and serum lipid levels are different in the Mulao and Han populations. Sex (male)-specific association of rs16996148 SNP in the NCAN/CILP2/PBX4 and serum lipid levels is also observed in the both ethnic groups.

A global view of porcine transcriptome in three tissues from a full-sib pair with extreme phenotypes in growth and fat deposition by paired-end RNA sequencing

Background

Elucidation of the pig transcriptome is essential for interpreting functional elements of the genome and understanding the genetic architecture of complex traits such as fat deposition, metabolism and growth.

Results

Here we used massive parallel high-throughput RNA sequencing to generate a high-resolution map of the porcine mRNA and miRNA transcriptome in liver, longissimus dorsi and abdominal fat from two full-sib F₂ hybrid pigs with segregated phenotypes on growth, blood physiological and biochemical parameters, and fat deposition. We obtained 8,508,418-10,219,332 uniquely mapped reads that covered 78.0% of the current annotated transcripts and identified 48,045-122,931 novel transcript fragments, which constituted 17,085-29,499 novel transcriptional active regions in six tested samples. We found that about 18.8% of the annotated genes showed alternative splicing patterns, and alternative 3' splicing is the most common type of alternative splicing events in pigs. Cross-tissue comparison revealed that many transcriptional events are tissue-differential and related to important biological functions in their corresponding tissues. We also detected a total of 164 potential novel miRNAs, most of which were tissue-specifically identified. Integrated analysis of genome-wide association study and differential gene expression revealed interesting candidate genes for complex traits, such as IGF2, CYP1A1, CKM and CES1 for heart weight, hemoglobin, pork pH value and serum cholesterol, respectively.

Conclusions

This study provides a global view of the complexity of the pig transcriptome, and gives an extensive new knowledge about alternative splicing, gene boundaries and miRNAs in pigs. Integrated analysis of genome wide association study and differential gene expression allows us to find important candidate genes for porcine complex traits.

Association between plasma metabolites and gene expression profiles in five porcine endocrine tissues

Background

Endocrine tissues play a fundamental role in maintaining homeostasis of plasma metabolites such as non-esterified fatty acids and glucose, the levels of which reflect the energy balance or the health status of animals. However, the relationship between the transcriptome of endocrine tissues and plasma metabolites has been poorly studied.

Methods

We determined the blood levels of 12 plasma metabolites in 27 pigs belonging to five breeds, each breed consisting of both females and males. The transcriptome of five endocrine tissues i.e. hypothalamus, adenohypophysis, thyroid gland, gonads and backfat tissues from 16 out of the 27 pigs was also determined. Sex and breed effects on the 12 plasma metabolites were investigated and associations between genes expressed in the five endocrine tissues and the 12 plasma metabolites measured were analyzed. A probeset was defined as a quantitative trait transcript (QTT) when its association with a particular metabolic trait achieved a nominal P value < 0.01.

Results

A larger than expected number of QTT was found for non-esterified fatty acids and alanine aminotransferase in at least two tissues. The associations were highly tissue-specific. The QTT within the tissues were divided into co-expression network modules enriched for genes in Kyoto Encyclopedia of Genes and Genomes or gene ontology categories that are related to the physiological functions of the corresponding tissues. We also explored a multi-tissue co-expression network using QTT for non-esterified fatty acids from the five tissues and found that a module, enriched in hypothalamus QTT, was positioned at the centre of the entire multi-tissue network.

Conclusions

These results emphasize the relationships between endocrine tissues and plasma metabolites in terms of gene expression. Highly tissue-specific association patterns suggest that candidate genes or gene pathways should be investigated in the context of specific tissues.

Detection of quantitative trait loci affecting serum cholesterol, LDL, HDL, and triglyceride in pigs

Background

Serum lipids are associated with many serious cardiovascular diseases and obesity problems. Many quantitative trait loci (QTL) have been reported in the pig mostly for performance traits but very few for the serum lipid traits. In contrast, remarkable numbers of QTL are mapped for serum lipids in humans and mice. Therefore, the objective of this research was to investigate the chromosomal regions influencing the serum level of the total cholesterol (CT), triglyceride (TG), high density protein cholesterol (HDL) and low density protein cholesterol (LDL) in pigs. For this purpose, a total of 330 animals from a Duroc × Pietrain F2 resource population were phenotyped for serum lipids using ELISA and were genotyped by using 122 microsatellite markers covering all porcine autosomes for QTL study in QTL Express. Blood sampling was performed at approximately 175 days before slaughter of the pig.

Results

Most of the traits were correlated with each other and were influenced by average daily gain, slaughter date and age. A total of 18 QTL including three QTL with imprinting effect were identified on 11 different porcine autosomes. Most of the QTL reached to 5% chromosome-wide (CW) level significance including a QTL at 5% experiment-wide (GW) and a QTL at 1% GW level significance. Of these QTL four were identified for both the CT and LDL and two QTL were identified for both the TG and LDL. Moreover, three chromosomal regions were detected for the HDL/LDL ratio in this study. One QTL for HDL on SSC2 and two QTL for TG on SSC11 and 17 were detected with imprinting effect. The highly significant QTL (1% GW) was detected for LDL at 82 cM on SSC1, whereas significant QTL (5% GW) was identified for HDL/LDL on SSC1 at 87 cM. Chromosomal regions with pleiotropic effects were detected for correlated traits on SSC1, 7 and 12. Most of the QTL identified for serum lipid traits correspond with the previously reported QTL for similar traits in other mammals. Two novel QTL on SSC16 for HDL and HDL/LDL ratio and an imprinted QTL on SSS17 for TG were detected in the pig for the first time.

Conclusion

The newly identified QTL are potentially involved in lipid metabolism. The results of this work shed new light on the genetic background of serum lipid concentrations and these findings will be helpful to identify candidate genes in these QTL regions related to lipid metabolism and serum lipid concentrations in pigs.