Identification of nonsense variants in the genomes of 15 Murciano-Granadina bucks and analysis of their segregation in parent-offspring trios.

Nonsense variants can inactivate gene function by causing the synthesis of truncated proteins or by inducing nonsense mediated decay of messenger RNAs. The occurrence of such variants in the genomes of livestock species is modulated by multiple demographic and selective factors. Even though nonsense variants can have causal effects on embryo lethality, abortions, and disease, their genomic distribution and segregation in domestic goats have not been characterized in depth yet. In this work, we have sequenced the genomes of 15 Murciano-Granadina bucks with an average coverage of 32.92 × ± 1.45 × . Bioinformatic analysis revealed 947 nonsense variants consistently detected with SnpEff and Ensembl-VEP. These variants were especially abundant in the 3'end of the protein-coding regions. Genes related to olfactory perception, ATPase activity coupled to transmembrane movement of substances, defense to virus, hormonal response, and sensory perception of taste were particularly enriched in nonsense variants. Seventeen nonsense variants expected to have harmful effects on fitness were genotyped in parent-offspring trios. We observed that several nonsense variants predicted to be lethal based on mouse knockout data did not have such effect, a finding that could be explained by the existence of multiple mechanisms counteracting lethality. These findings demonstrate that predicting the effects of putative nonsense variants on fitness is extremely challenging. As a matter of fact, such a goal could only be achieved by generating a high quality telomere-to-telomere goat reference genome combined with carefully curated annotation and functional testing of promising candidate variants.

Comparing the mRNA expression profile and the genetic determinism of intramuscular fat traits in the porcine gluteus medius and longissimus dorsi muscles.

BACKGROUND: Intramuscular fat (IMF) content and composition have a strong impact on the nutritional and organoleptic properties of porcine meat. The goal of the current work was to compare the patterns of gene expression and the genetic determinism of IMF traits in the porcine gluteus medius (GM) and longissimus dorsi (LD) muscles. RESULTS: A comparative analysis of the mRNA expression profiles of the pig GM and LD muscles in 16 Duroc pigs with available microarray mRNA expression measurements revealed the existence of 106 differentially expressed probes (fold-change > 1.5 and q-value < 0.05). Amongst the genes displaying the most significant differential expression, several loci belonging to the Hox transcription factor family were either upregulated (HOXA9, HOXA10, HOXB6, HOXB7 and TBX1) or downregulated (ARX) in the GM muscle. Differences in the expression of genes with key roles in carbohydrate and lipid metabolism (e.g. FABP3, ORMDL1 and SLC37A1) were also detected. By performing a GWAS for IMF content and composition traits recorded in the LD and GM muscles of 350 Duroc pigs, we identified the existence of one region on SSC14 (110-114 Mb) displaying significant associations with C18:0, C18:1(n-7), saturated and unsaturated fatty acid contents in both GM and LD muscles. Moreover, we detected several genome-wide significant associations that were not consistently found in both muscles. Further studies should be performed to confirm whether these associations are muscle-specific. Finally, the performance of an eQTL scan for 74 genes, located within GM QTL regions and with available microarray measurements of gene expression, made possible to identify 14 cis-eQTL regulating the expression of 14 loci, and six of them were confirmed by RNA-Seq. CONCLUSIONS: We have detected significant differences in the mRNA expression patterns of the porcine LD and GM muscles, evidencing that the transcriptomic profile of the skeletal muscle tissue is affected by anatomical, metabolic and functional factors. A highly significant association with IMF composition on SSC14 was replicated in both muscles, highlighting the existence of a common genetic determinism, but we also observed the existence of a few associations whose magnitude and significance varied between LD and GM muscles.

Identification of protein-damaging mutations in 10 swine taste receptors and 191 appetite-reward genes.

BACKGROUND: Taste receptors (TASRs) are essential for the body's recognition of chemical compounds. In the tongue, TASRs sense the sweet and umami and the toxin-related bitter taste thus promoting a particular eating behaviour. Moreover, their relevance in other organs is now becoming evident. In the intestine, they regulate nutrient absorption and gut motility. Upon ligand binding, TASRs activate the appetite-reward circuitry to signal the nervous system and keep body homeostasis. With the aim to identify genetic variation in the swine TASRs and in the genes from the appetite and the reward pathways, we have sequenced the exons of 201 TASRs and appetite-reward genes from 304 pigs belonging to ten breeds, wild boars and to two phenotypically extreme groups from a F2 resource with data on growth and fat deposition. RESULTS: We identified 2,766 coding variants 395 of which were predicted to have a strong impact on protein sequence and function. 334 variants were present in only one breed and at predicted alternative allele frequency (pAAF) ≥ 0.1. The Asian pigs and the wild boars showed the largest proportion of breed specific variants. We also compared the pAAF of the two F2 groups and found that variants in TAS2R39 and CD36 display significant differences suggesting that these genes could influence growth and fat deposition. We developed a 128-variant genotyping assay and confirmed 57 of these variants. CONCLUSIONS: We have identified thousands of variants affecting TASRs as well as genes involved in the appetite and the reward mechanisms. Some of these genes have been already associated to taste preferences, appetite or behaviour in humans and mouse. We have also detected indications of a potential relationship of some of these genes with growth and fat deposition, which could have been caused by changes in taste preferences, appetite or reward and ultimately impact on food intake. A genotyping array with 57 variants in 31 of these genes is now available for genotyping and start elucidating the impact of genetic variation in these genes on pig biology and breeding.

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.

Analysis of porcine MUC4 gene as a candidate gene for prolificacy QTL on SSC13 in an Iberian × Meishan F2 population.

BACKGROUND: Reproductive traits, such as prolificacy, are of great interest to the pig industry. Better understanding of their genetic architecture should help to increase the efficiency of pig productivity through the implementation of marker assisted selection (MAS) programmes. RESULTS: The Mucin 4 (MUC4) gene has been evaluated as a candidate gene for a prolificacy QTL described in an Iberian × Meishan (Ib × Me) F2 intercross. For association analyses, two previously described SNPs (DQ124298:g.243A>G and DQ124298:g.344A>G) were genotyped in 347 pigs from the Ib × Me population. QTL for the number of piglets born alive (NBA) and for the total number of piglets born (TNB) were confirmed on SSC13 at positions 44 cM and 51 cM, respectively. The MUC4 gene was successfully located within the confidence intervals of both QTL. Only DQ124298:g.344A>G MUC4 polymorphism was significantly associated with both NBA and TNB (P-value < 0.05) with favourable effects coming from the Meishan origin. MUC4 expression level was determined in F2 sows displaying extreme phenotypes for the number of embryos (NE) at 30-32 days of gestation. Differences in the uterine expression of MUC4 were found between high (NE ≥ 13) and low (NE ≤ 11) prolificacy sows. Overall, MUC4 expression in high prolificacy sows was almost two-fold increased compared with low prolificacy sows. CONCLUSIONS: Our data suggest that MUC4 could play an important role in the establishment of an optimal uterine environment that would increase embryonic survival during pig gestation.

Copy number variation in the porcine genome inferred from a 60 k SNP BeadChip.

BACKGROUND: Recent studies in pigs have detected copy number variants (CNVs) using the Comparative Genomic Hybridization technique in arrays designed to cover specific porcine chromosomes. The goal of this study was to identify CNV regions (CNVRs) in swine species based on whole genome SNP genotyping chips. RESULTS: We used predictions from three different programs (cnvPartition, PennCNV and GADA) to analyze data from the Porcine SNP60 BeadChip. A total of 49 CNVRs were identified in 55 animals from an Iberian x Landrace cross (IBMAP) according to three criteria: detected in at least two animals, contained three or more consecutive SNPs and recalled by at least two programs. Mendelian inheritance of CNVRs was confirmed in animals belonging to several generations of the IBMAP cross. Subsequently, a segregation analysis of these CNVRs was performed in 372 additional animals from the IBMAP cross and its distribution was studied in 133 unrelated pig samples from different geographical origins. Five out of seven analyzed CNVRs were validated by real time quantitative PCR, some of which coincide with well known examples of CNVs conserved across mammalian species. CONCLUSIONS: Our results illustrate the usefulness of Porcine SNP60 BeadChip to detect CNVRs and show that structural variants can not be neglected when studying the genetic variability in this species.

A high throughput genotyping approach reveals distinctive autosomal genetic signatures for European and Near Eastern wild boar.

The lack of a Near Eastern genetic signature in modern European porcine breeds indicates that, although domestic pigs from the Fertile Crescent entered Europe during the Neolithic, they were completely replaced by their European counterparts in a short window of time. Whilst the absence of such genetic signature has been convincingly demonstrated at the mitochondrial level, variation at the autosomal genomes of European and Near Eastern Sus scrofa has not been compared yet. Herewith, we have explored the genetic relationships among 43 wild boar from Europe (N = 21), Near East (N = 19) and Korea (N = 3), and 40 Iberian (N = 16), Canarian (N = 4) and Mangalitza (N = 20) pigs by using a high throughput SNP genotyping platform. After data filtering, 37,167 autosomal SNPs were used to perform population genetics analyses. A multidimensional scaling plot based on genome-wide identity-by-state pairwise distances inferred with PLINK showed that Near Eastern and European wild boar populations are genetically differentiated. Maximum likelihood trees built with TreeMix supported this conclusion i.e. an early population split between Near Eastern and European Sus scrofa was observed. Moreover, analysis of the data with Structure evidenced that the sampled Iberian, Canarian and Mangalitza pigs did not carry any autosomal signature compatible with a Near Eastern ancestry, a finding that agrees well with previous mitochondrial studies.

On growth, fatness, and form: a further look at porcine chromosome 4 in an Iberian x Landrace cross.

A crossed population between Iberian x Landrace pigs consisting of 321 F2, 87 F3, and 85 backcross individuals has been analyzed to refine the number and positions of quantitative trait loci (QTL) affecting shape, growth, fatness, and meat quality traits in SSC4. A multitrait multi-QTL approach has been used. Our results suggest that carcass length and shoulder weight are affected by two loci. The first one, close to the AFABP gene, has a very strong pleiotropic effect on fatness, whereas the second one, in the interval between S0073 and S0214, also affects live weight, although to a lesser extent. This latter QTL would correspond to the FAT1 locus described initially in pigs. It seems that SSC4's loci play an important role in redistributing total weight, and the Landrace allele increases shoulder weight and carcass length much more than ham or total weight. Furthermore, there is also strong evidence of additional loci influencing pH and color in more distant, telomeric positions.