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

Single and longitudinal genome-wide association studies for dairy traits available in goats with three recorded lactations

Milk yield and composition phenotypes are systematically recorded across several lactations in goats, but the majority of genome-wide association studies (GWAS) performed so far have rather ignored the longitudinal nature of such data. Here, we have used two different GWAS approaches to analyse data from three lactations recorded in Murciano-Granadina goats. In Analysis 1, independent GWAS have been carried out for each trait and lactation, while a single longitudinal GWAS, jointly considering all data, has been performed in Analysis 2. In both analyses, genome-wide significant QTL for lactose percentage on chromosome 2 (129.77-131.01 Mb) and for milk protein percentage on the chromosome 6 (74.8-94.6 Mb) casein gene cluster region were detected. In Analysis 1, several QTL were not replicated in all three lactations, possibly due to the existence of lactation-specific genetic determinants. In Analysis 2, we identified several genome-wide significant QTL related to milk yield and protein content that were not uncovered in Analysis 1. The increased number of QTL identified in Analysis 2 suggests that the longitudinal GWAS is particularly well suited for the genetic analysis of dairy traits. Moreover, our data confirm that variability within or close to the casein complex is the main genetic determinant of milk protein percentage in Murciano-Granadina goats.

Host genetics and gut microbiota jointly regulate blood biochemical indicators in chickens

Blood biochemical indicators play a crucial role in assessing an individual's overall health status and metabolic function. In this study, we measured five blood biochemical indicators, including total cholesterol (CHOL), low-density lipoprotein cholesterol (LDL-CH), triglycerides (TG), high-density lipoprotein cholesterol (HDL-CH), and blood glucose (BG), as well as 19 growth traits of 206 male chickens. By integrating host whole-genome information and 16S rRNA sequencing of the duodenum, jejunum, ileum, cecum, and feces microbiota, we assessed the contributions of host genetics and gut microbiota to blood biochemical indicators and their interrelationships. Our results demonstrated significant negative phenotypic and genetic correlations (r =  - 0.20 ~  - 0.67) between CHOL and LDL-CH with growth traits such as body weight, abdominal fat content, muscle content, and shin circumference. The results of heritability and microbiability indicated that blood biochemical indicators were jointly regulated by host genetics and gut microbiota. Notably, the heritability of HDL-CH was estimated to be 0.24, while the jejunal microbiability for BG and TG reached 0.45 and 0.23. Furthermore, by conducting genome-wide association study (GWAS) with the single-nucleotide polymorphism (SNPs), insertion/deletion (indels), and structural variation (SV), we identified RAP2C, member of the RAS oncogene family (RAP2C), dedicator of cytokinesis 11 (DOCK11), neurotensin (NTS) and BOP1 ribosomal biogenesis factor (BOP1) as regulators of HDL-CH, and glycerophosphodiester phosphodiesterase domain containing 5 (GDPD5), dihydrodiol dehydrogenase (DHDH), and potassium voltage-gated channel interacting protein 1 (KCNIP1) as candidate genes of BG. Moreover, our findings suggest that cecal RF39 and Clostridia_UCG_014 may be linked to the regulation of CHOL, and jejunal Streptococcaceae may be involved in the regulation of TG. Additionally, microbial GWAS results indicated that the presence of gut microbiota was under host genetic regulation. Our findings provide valuable insights into the complex interaction between host genetics and microbiota in shaping the blood biochemical profile of chickens. KEY POINTS: • Multiple candidate genes were identified for the regulation of CHOL, HDL-CH, and BG. • RF39, Clostridia_UCG_014, and Streptococcaceae were implicated in CHOL and TG modulation. • The composition of gut microbiota is influenced by host genetics.

Genomic prediction of growth in a commercially, recreationally, and culturally important marine resource, the Australian snapper (Chrysophrys auratus)

Growth is one of the most important traits of an organism. For exploited species, this trait has ecological and evolutionary consequences as well as economical and conservation significance. Rapid changes in growth rate associated with anthropogenic stressors have been reported for several marine fishes, but little is known about the genetic basis of growth traits in teleosts. We used reduced genome representation data and genome-wide association approaches to identify growth-related genetic variation in the commercially, recreationally, and culturally important Australian snapper (Chrysophrys auratus, Sparidae). Based on 17,490 high-quality single-nucleotide polymorphisms and 363 individuals representing extreme growth phenotypes from 15,000 fish of the same age and reared under identical conditions in a sea pen, we identified 100 unique candidates that were annotated to 51 proteins. We documented a complex polygenic nature of growth in the species that included several loci with small effects and a few loci with larger effects. Overall heritability was high (75.7%), reflected in the high accuracy of the genomic prediction for the phenotype (small vs large). Although the single-nucleotide polymorphisms were distributed across the genome, most candidates (60%) clustered on chromosome 16, which also explains the largest proportion of heritability (16.4%). This study demonstrates that reduced genome representation single-nucleotide polymorphisms and the right bioinformatic tools provide a cost-efficient approach to identify growth-related loci and to describe genomic architectures of complex quantitative traits. Our results help to inform captive aquaculture breeding programs and are of relevance to monitor growth-related evolutionary shifts in wild populations in response to anthropogenic pressures.

High-Altitude Genetic Selection and Genome-Wide Association Analysis of Yield-Related Traits in Elymus sibiricus L. Using SLAF Sequencing

The genetic adaptations to harsh climatic conditions in high altitudes and genetic basis of important agronomic traits are poorly understood in Elymus sibiricus L. In this study, an association population of 210 genotypes was used for population structure, selective sweep analysis, and genome-wide association study (GWAS) based on 88,506 single nucleotide polymorphisms (SNPs). We found 965 alleles under the natural selection of high altitude, which included 7 hub genes involved in the response to UV, and flavonoid and anthocyanin biosynthetic process based on the protein-protein interaction (PPI) analysis. Using a mixed linear model (MLM), the GWAS test identified a total of 1,825 significant loci associated with 12 agronomic traits. Based on the gene expression data of two wheat cultivars and the PPI analysis, we finally identified 12 hub genes. Especially, in plant height traits, the top hub gene (TOPLESS protein) encoding auxins and jasmonic acid signaling pathway, shoot apical meristem specification, and xylem and phloem pattern formation was highly overexpressed. These genes might play essential roles in controlling the growth and development of E. sibiricus. Therefore, this study provides fundamental insights relevant to hub genes and will benefit molecular breeding and improvement in E. sibiricus and other Elymus species.

PARPs in lipid metabolism and related diseases

PARPs and tankyrases (TNKS) represent a family of 17 proteins. PARPs and tankyrases were originally identified as DNA repair factors, nevertheless, recent advances have shed light on their role in lipid metabolism. To date, PARP1, PARP2, PARP3, tankyrases, PARP9, PARP10, PARP14 were reported to have multi-pronged connections to lipid metabolism. The activity of PARP enzymes is fine-tuned by a set of cholesterol-based compounds as oxidized cholesterol derivatives, steroid hormones or bile acids. In turn, PARPs modulate several key processes of lipid homeostasis (lipotoxicity, fatty acid and steroid biosynthesis, lipoprotein homeostasis, fatty acid oxidation, etc.). PARPs are also cofactors of lipid-responsive nuclear receptors and transcription factors through which PARPs regulate lipid metabolism and lipid homeostasis. PARP activation often represents a disruptive signal to (lipid) metabolism, and PARP-dependent changes to lipid metabolism have pathophysiological role in the development of hyperlipidemia, obesity, alcoholic and non-alcoholic fatty liver disease, type II diabetes and its complications, atherosclerosis, cardiovascular aging and skin pathologies, just to name a few. In this synopsis we will review the evidence supporting the beneficial effects of pharmacological PARP inhibitors in these diseases/pathologies and propose repurposing PARP inhibitors already available for the treatment of various malignancies.

The role of ADP-ribose metabolism in metabolic regulation, adipose tissue differentiation, and metabolism

In this review, Szanto et al. summarize the metabolic regulatory roles of PARP enzymes and their associated pathologies.

Poly(ADP-ribose) polymerases (PARPs or ARTDs), originally described as DNA repair factors, have metabolic regulatory roles. PARP1, PARP2, PARP7, PARP10, and PARP14 regulate central and peripheral carbohydrate and lipid metabolism and often channel pathological disruptive metabolic signals. PARP1 and PARP2 are crucial for adipocyte differentiation, including the commitment toward white, brown, or beige adipose tissue lineages, as well as the regulation of lipid accumulation. Through regulating adipocyte function and organismal energy balance, PARPs play a role in obesity and the consequences of obesity. These findings can be translated into humans, as evidenced by studies on identical twins and SNPs affecting PARP activity.

Detection of homozygous genotypes for a putatively lethal recessive mutation in the porcine argininosuccinate synthase 1 (ASS1) gene

The sequencing of the pig genome revealed the existence of homozygous individuals for a nonsense mutation in the argininosuccinate synthase 1 (ASS1) gene (rs81212146, c.944T>A, L315X). Paradoxically, an AA homozygous genotype for this polymorphism is expected to abolish the function of the ASS1 enzyme that participates in the urea cycle, leading to citrullinemia, hyperammonemia, coma and death. Sequencing of five Duroc boars that sired a population of 350 Duroc barrows revealed the segregation of the c.944T>A polymorphism, so we aimed to investigate its phenotypic consequences. Genotyping of this mutation in the 350 Duroc barrows revealed the existence of seven individuals homozygous (AA) for the nonsense mutation. These AA pigs had a normal weight despite the fact that mild citrullinemia often involves impaired growth. Sequencing of the region surrounding the mutation in TT, TA and AA individuals revealed that the A substitution in the second position of the codon (c.944T>A) is in complete linkage disequilibrium with a C replacement (c.943T>C) in the first position of the codon. This second mutation would compensate for the potentially damaging effect of the c.944T>A replacement. In fact, this is the most probable reason why pigs with homozygous AA genotypes at the 944 site of the ASS1 coding region are alive. Our results illustrate the complexities of predicting the consequences of nonsense mutations on gene function and phenotypes, not only because of annotation issues but also owing to the existence of genetic mechanisms that sometimes limit the penetrance of highly harmful mutations.

Polymorphisms of the cryptochrome 2 and mitoguardin 2 genes are associated with the variation of lipid-related traits in Duroc pigs

The genetic factors determining the phenotypic variation of porcine fatness phenotypes are still largely unknown. We investigated whether the polymorphism of eight genes (MIGA2, CRY2, NPAS2, CIART, ARNTL2, PER1, PER2 and PCK1), which display differential expression in the skeletal muscle of fasted and fed sows, is associated with the variation of lipid and mRNA expression phenotypes in Duroc pigs. The performance of an association analysis with the GEMMA software demonstrated that the rs330779504 SNP in the MIGA2 gene is associated with LDL concentration at 190 days (LDL₂, corrected P-value = 0.057). Moreover, the rs320439526 SNP of the CRY2 gene displayed a significant association with stearic acid content in the longissimus dorsi muscle (LD C18:0, corrected P-value = 0.015). Both SNPs were also associated with the mRNA levels of the corresponding genes in the gluteus medius skeletal muscle. From a biological perspective these results are meaningful because MIGA2 protein plays an essential role in mitochondrial fusion, a process tightly connected with the energy status of the cell, while CRY2 is a fundamental component of the circadian clock. However, inclusion of these two SNPs in chromosome-wide association analyses demonstrated that they are not located at the peaks of significance for the two traits under study (LDL₂ for rs330779504 and LD C18:0 for rs320439526), thus implying that these two SNPs do not have causal effects.

Genome-wide association studies for 30 haematological and blood clinical-biochemical traits in Large White pigs reveal genomic regions affecting intermediate phenotypes

Haematological and clinical-biochemical parameters are considered indicators of the physiological/health status of animals and might serve as intermediate phenotypes to link physiological aspects to production and disease resistance traits. The dissection of the genetic variability affecting these phenotypes might be useful to describe the resilience of the animals and to support the usefulness of the pig as animal model. Here, we analysed 15 haematological and 15 clinical-biochemical traits in 843 Italian Large White pigs, via three genome-wide association scan approaches (single-trait, multi-trait and Bayesian). We identified 52 quantitative trait loci (QTLs) associated with 29 out of 30 analysed blood parameters, with the most significant QTL identified on porcine chromosome 14 for basophil count. Some QTL regions harbour genes that may be the obvious candidates: QTLs for cholesterol parameters identified genes (ADCY8, APOB, ATG5, CDKAL1, PCSK5, PRL and SOX6) that are directly involved in cholesterol metabolism; other QTLs highlighted genes encoding the enzymes being measured [ALT (known also as GPT) and AST (known also as GOT)]. Moreover, the multivariate approach strengthened the association results for several candidate genes. The obtained results can contribute to define new measurable phenotypes that could be applied in breeding programs as proxies for more complex traits.

Poly(ADP-ribose) polymerase-2 is a lipid-modulated modulator of muscular lipid homeostasis

There is a growing body of evidence that poly(ADP-ribose) polymerase-2 (PARP2), although originally described as a DNA repair protein, has a widespread role as a metabolic regulator. We show that the ablation of PARP2 induced characteristic changes in the lipidome. The silencing of PARP2 induced the expression of sterol regulatory element-binding protein-1 and -2 and initiated de novo cholesterol biosynthesis in skeletal muscle. Increased muscular cholesterol was shunted to muscular biosynthesis of dihydrotestosterone, an anabolic steroid. Thus, skeletal muscle fibers in PARP2^-/- mice were stronger compared to those of their wild-type littermates. In addition, we detected changes in the dynamics of the cell membrane, suggesting that lipidome changes also affect the biophysical characteristics of the cell membrane. In in silico and wet chemistry studies, we identified lipid species that can decrease the expression of PARP2 and potentially phenocopy the genetic abruption of PARP2, including artificial steroids. In view of these observations, we propose a new role for PARP2 as a lipid-modulated regulator of lipid metabolism.

Haplotypes on pig chromosome 3 distinguish metabolically healthy from unhealthy obese individuals

We have established a pig resource population specifically designed to elucidate the genetics involved in development of obesity and obesity related co-morbidities by crossing the obesity prone Göttingen Minipig breed with two lean production pig breeds. In this study we have performed genome wide association (GWA) to identify loci with effect on blood lipid levels. The most significantly associated single nucleotide polymorphisms (SNPs) were used for linkage disequilibrium (LD) and haplotype analyses. Three separate haploblocks which influence the ratio between high density lipoprotein cholesterol and total cholesterol (HDL-C/CT), triglycerides (TG) and low density lipoprotein cholesterol (LDL-C) levels respectively were identified on Sus Scrofa chromosome 3 (SSC3). Large additive genetic effects were found for the HDL-C/CT and LDL-C haplotypes. Haplotypes segregating from Göttingen Minipigs were shown to impose a positive effect on blood lipid levels. Thus, the genetic profile of the Göttingen Minipig breed seems to support a phenotype comparable to the metabolic healthy obese (MHO) phenotype in humans.

Pleiotropy, constraint, and modularity in the evolution of life histories: insights from genomic analyses

Multicellular organisms display an enormous range of life history (LH) strategies and present an evolutionary conundrum; despite strong natural selection, LH traits are characterized by high levels of genetic variation. To understand the evolution of life histories and maintenance of this variation, the specific phenotypic effects of segregating alleles and the genetic networks in which they act need to be elucidated. In particular, the extent to which LH evolution is constrained by the pleiotropy of alleles contributing to LH variation is generally unknown. Here, we review recent empirical results that shed light on this question, with an emphasis on studies employing genomic analyses. While genome-scale analyses are increasingly practical and affordable, they face limitations of genetic resolution and statistical power. We describe new research approaches that we believe can produce new insights and evaluate their promise and applicability to different kinds of organisms. Two approaches seem particularly promising: experiments that manipulate selection in multiple dimensions and measure phenotypic and genomic response and analytical approaches that take into account genome-wide associations between markers and phenotypes, rather than applying a traditional marker-by-marker approach.

Joint QTL mapping and gene expression analysis identify positional candidate genes influencing pork quality traits

Meat quality traits have an increasing importance in the pig industry because of their strong impact on consumer acceptance. Herewith, we have combined phenotypic and microarray expression data to map loci with potential effects on five meat quality traits recorded in the longissimus dorsi (LD) and gluteus medius (GM) muscles of 350 Duroc pigs, i.e. pH at 24 hours post-mortem (pH₂₄), electric conductivity (CE) and muscle redness (a*), lightness (L*) and yellowness (b*). We have found significant genome-wide associations for CE of LD on SSC4 (~104 Mb), SSC5 (~15 Mb) and SSC13 (~137 Mb), while several additional regions were significantly associated with meat quality traits at the chromosome-wide level. There was a low positional concordance between the associations found for LD and GM traits, a feature that reflects the existence of differences in the genetic determinism of meat quality phenotypes in these two muscles. The performance of an eQTL search for SNPs mapping to the regions associated with meat quality traits demonstrated that the GM a* SSC3 and pH₂₄ SSC17 QTL display positional concordance with cis-eQTL regulating the expression of several genes with a potential role on muscle metabolism.

Metabolic roles of poly(ADP-ribose) polymerases

Poly(ADP-ribosyl)ation (PARylation) is an evolutionarily conserved reaction that had been associated with numerous cellular processes such as DNA repair, protein turnover, inflammatory regulation, aging or metabolic regulation. The metabolic regulatory tasks of poly(ADP-ribose) polymerases (PARPs) are complex, it is based on the regulation of metabolic transcription factors (e.g. SIRT1, nuclear receptors, SREBPs) and certain cellular energy sensors. PARP over-activation can cause damage to mitochondrial terminal oxidation, while the inhibition of PARP-1 or PARP-2 can induce mitochondrial oxidation by enhancing the mitotropic tone of gene transcription and signal transduction. These PARP-mediated processes impact on higher order metabolic regulation that modulates lipid metabolism, circadian oscillations and insulin secretion and signaling. PARP-1, PARP-2 and PARP-7 are related to metabolic diseases such as diabetes, alcoholic and non-alcoholic fatty liver disease (AFLD, NAFLD), or on a broader perspective to Warburg metabolism in cancer or the metabolic diseases accompanying aging.

Identification and characterization of novel single nucleotide polymorphism markers for fat deposition in muscle tissue of pigs using amplified fragment length polymorphism

Objective

This study was conducted to identify and evaluate the effective single nucleotide polymorphism (SNP) markers for fat deposition in the longissimus dorsi muscles of pigs using the amplified fragment length polymorphism (AFLP) approach.

Methods

Sixty-four selective primer combinations were used to identify the AFLP markers in the 20 highest- and 20 lowest-intramuscular fat (IMF) content phenotypes. Five AFLP fragments were converted into simple codominant SNP markers. These SNP markers were tested in terms of their association with IMF content and fatty acid (FA) composition traits in 620 commercially crossbred pigs.

Results

The SSC7 g.4937240C>G marker showed an association with IMF content (p<0.05). The SSC9 g.5496647_5496662insdel marker showed a significant association with IMF content and arachidonic levels (p<0.05). The SSC10 g.71225134G>A marker revealed an association with palmitoleic and ω9 FA levels (p<0.05), while the SSC17 g.61976696G>T marker showed a significant association with IMF content and FA levels of palmitoleic, eicosenoic, arachidonic, monounsaturated fatty acids, and ω9 FA levels. However, no significant association of SSC8 g.47338181G>A was observed with any IMF and FA levels in this study.

Conclusion

Four SNP markers (SSC7 g.4937240C>G, SSC9 g.5496647_5496662insdel, SSC10 g.71225134G>A, and SSC17 g.61976696G>T) were found to be associated with IMF and/or FA content traits in commercially crossbred pigs. These findings provide evidence of the novel SNP markers as being potentially useful for selecting pigs with the desirable IMF content and FA composition.

Association between the DOCK7, PCSK9 and GALNT2 Gene Polymorphisms and Serum Lipid levels

This study was to determine the association between several single nucleotide polymorphisms (SNPs) in the dedicator of cytokinesis 7 (DOCK7), proprotein convertase subtilisin/kexin type 9 (PCSK9) and polypeptide N-acetylgalactosaminyltransferase 2 (GALNT2) and serum lipid levels. Genotyping of 9 SNPs was performed in 881 Jing subjects and 988 Han participants. Allele and genotype frequencies of the detected SNPs were different between the two populations. Several SNPs were associated with triglyceride (TG, rs10889332, rs615563, rs7552841, rs1997947, rs2760537, rs4846913 and rs11122316), high-density lipoprotein (HDL) cholesterol (rs1997947), low-density lipoprotein (LDL) cholesterol (rs1168013 and rs7552841), apolipoprotein (Apo) A1 (rs1997947), ApoB (rs10889332 and rs7552841), and ApoA1/ApoB ratio (rs7552841) in Jing minority; and with TG (rs10889332, rs615563, rs7552841, rs11206517, rs1997947, rs4846913 and rs11122316), HDL cholesterol (rs11206517 and rs4846913), LDL cholesterol (rs1168013), ApoA1 (rs11206517 and rs4846913), ApoB (rs7552841), and ApoA1/ApoB ratio (rs4846913) in Han nationality. Strong linkage disequilibria were noted among the SNPs. The commonest haplotype was G-C-G-C-T-G-C-C-G (>10%). The frequencies of C-C-G-C-T-G-T-C-G, G-C-A-C-T-G-C-C-G, G-C-G-C-T-A-C-C-A, G-C-G-C-T-G-C-C-A, G-C-G-C-T-G-T-C-A haplotypes were different between the two populations. Haplotypes could explain much more serum lipid variation than any single SNP alone especially for TG. Differences in lipid profiles between the two populations might partially attribute to these SNPs and their haplotypes.

Functional and association studies of the cholesteryl ester transfer protein (CETP) gene in a Wannan Black pig model

Some polymorphisms of the human CETP gene are causally and significantly associated with serum lipids levels; however, the information regarding this gene in pigs is sparse. To evaluate the effects of CETP on blood lipid traits and fat deposition in pig, porcine CETP tissue expression patterns were observed by quantitative real-time polymerase chain reaction (qPCR) first. High expression was detected in liver, spleen, gluteus medius (GM) muscle and backfat. A de novo polymorphism (AF333037:g.795C>T) in the intron 1 region of porcine CETP was identified. This polymorphism was further genotyped by direct sequencing of the PCR products of 390 Wannan Black pigs, a Chinese native breed population. Association analyses at 45 and 300 days of age revealed highly significant associations between CETP genotypes and serum lipid traits. Furthermore, this polymorphism was proved to be associated with differences in liver CETP mRNA levels: pigs at 300 days of age with the TT genotype had higher levels than did those with other genotypes (P = 0.021). Additionally, analysis at 300 days of age showed that GM CETP mRNA expression correlated positively with serum lipids levels as well as with carcass backfat thickness and intramuscular fat content in GM. These results indicate that CETP is involved in serum, adipose and muscle lipid metabolism in pigs. The mechanisms underlying such relationships and their functional implications are worthy of further research.

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.