Genomic selection pressure discovery using site-frequency spectrum and reduced local variability statistics in Pakistani Dera-Din-Panah goat

BACKGROUND

Population geneticists have long sought to comprehend various selection traces accumulated in the goat genome due to natural or human driven artificial selection through breeding practices, which led the wild animals to domestication, so understanding evolutionary process may helpful to utilize the full genetic potential of goat genome.

METHODS AND RESULTS

As a step forward to pinpoint the selection signals in Pakistani Dera-Din-Panah (DDP) goat, whole-genome pooled sequencing (n = 12) was performed, and 618,236,192 clean paired-end reads were mapped against ARS1 reference goat assembly. Five different selection signature statistics were applied using four site-frequency spectrum (SFS) methods (Tajima's D ([Formula: see text]), Fay and Wu's H ([Formula: see text]), Zeng's E ([Formula: see text]), [Formula: see text]) and one reduced local variability approach named pooled heterozygosity ([Formula: see text]). The under-selection regions were annotated with significant threshold values of [Formula: see text]≥4.7, [Formula: see text]≥6, [Formula: see text]≥2.5, Pool-HMM ≥ 12, and [Formula: see text]≥5 that resulted in accumulative 364 candidate gene hits. The highest genomic selection signals were observed on Chr. 4, 6, 10, 12, 15, 16, 18, 20, and 27 and harbor ADAMTS6, CWC27, RELN, MYCBP2, FGF14, STIM1, CFAP74, GNB1, CALML6, TMEM52, FAM149A, NADK, MMP23B, OPN3, FH, MFHAS1, KLKB1, RRM1, KMO, SPEF2, F11, KIT, KMO, ERI1, ATP8B4, and RHOG genes. Next, the validation of our captured genomic hits was also performed by more than one applied statistics which harbor meat production, immunity, and reproduction associated genes to strengthen our hypothesis of under-selection traits in this Pakistani goat breed. Furthermore, common candidate genes captured by more than one statistical method were subjected to gene ontology and KEGG pathway analysis to get insights of particular biological processes associated with this goat breed.

CONCLUSION

Current perception of genomic architecture of DDP goat provides a better understanding to improve its genetic potential and other economically important traits of medium to large body size, milk, and fiber production by updating the genomic insight driven breeding strategies to boost the livestock and agriculture-based economy of the country.

Genomic prediction in pigs using data from a commercial crossbred population: insights from the Duroc x (Landrace x Yorkshire) three-way crossbreeding system

BACKGROUND

Genomic selection is widely applied for genetic improvement in livestock crossbreeding systems to select excellent nucleus purebred (PB) animals and to improve the performance of commercial crossbred (CB) animals. Most current predictions are based solely on PB performance. Our objective was to explore the potential application of genomic selection of PB animals using genotypes of CB animals with extreme phenotypes in a three-way crossbreeding system as the reference population. Using real genotyped PB as ancestors, we simulated the production of 100,000 pigs for a Duroc x (Landrace x Yorkshire) DLY crossbreeding system. The predictive performance of breeding values of PB animals for CB performance using genotypes and phenotypes of (1) PB animals, (2) DLY animals with extreme phenotypes, and (3) random DLY animals for traits of different heritabilities ([Formula: see text] = 0.1, 0.3, and 0.5) was compared across different reference population sizes (500 to 6500) and prediction models (genomic best linear unbiased prediction (GBLUP) and Bayesian sparse linear mixed model (BSLMM)).

RESULTS

Using a reference population consisting of CB animals with extreme phenotypes showed a definite predictive advantage for medium- and low-heritability traits and, in combination with the BSLMM model, significantly improved selection response for CB performance. For high-heritability traits, the predictive performance of a reference population of extreme CB phenotypes was comparable to that of PB phenotypes when the effect of the genetic correlation between PB and CB performance ([Formula: see text]) on the accuracy obtained with a PB reference population was considered, and the former could exceed the latter if the reference size was large enough. For the selection of the first and terminal sires in a three-way crossbreeding system, prediction using extreme CB phenotypes outperformed the use of PB phenotypes, while the optimal design of the reference group for the first dam depended on the percentage of individuals from the corresponding breed that the PB reference data comprised and on the heritability of the target trait.

CONCLUSIONS

A commercial crossbred population is promising for the design of the reference population for genomic prediction, and selective genotyping of CB animals with extreme phenotypes has the potential for maximizing genetic improvement for CB performance in the pig industry.

Association of ghrelin gene polymorphisms with slaughter traits in pig

It has been hypothesized that mutations in the ghrelin gene in pigs may play a similar role as in humans and may be associated with obesity. The aim of study was to analyse the polymorphisms of ghrelin gene and to evaluate its effect on pigs’ carcass traits. The effect of c.-93A>G, 4428T>C and g.4486C>T polymorphisms at the ghrelin gene on slaughter performance were analysed in 346 gilts represented by three breeds (Polish Landrace, Duroc, Pietrain). Animals were fattened from 30 to 100 (±2.5) kg body weight. After slaughter, the carcasses were chilled for 24 hours (4ºC), weighted and the right half-carcasses were dissected and evaluated. A number of data were obtained including: meat weight in primary cuts, weight of ham, backfat thickness and carcass yield. From breeding and production point of view, the favourable results were obtained for pigs with the GG genotype at the c.-93A>G locus, characterized by better carcass results than those with the AA genotype, e.g. higher ham weight and lower average backfat thickness. In pigs with the TT genotype at the g.4428T>C locus, we found lower mean backfat thickness than pigs with the CC genotype.

Mutations in melanocortin-4 receptor: From fish to men

Melanocortin-4 receptor (MC4R), expressed abundantly in the hypothalamus, is a critical regulator of energy homeostasis, including both food intake and energy expenditure. Shortly after the publication in 1997 of the Mc4r knockout phenotypes in mice, including increased food intake and severe obesity, the first mutations in MC4R were reported in humans in 1998. Studies in the subsequent two decades have established MC4R mutation as the most common monogenic form of obesity, especially in early-onset severe obesity. Studies in animals, from fish to mammals, have established the conserved physiological roles of MC4R in all vertebrates in regulating energy balance. Drug targeting MC4R has been recently approved for treating morbid genetic obesity. How the MC4R can be exploited for animal production is highly worthy of active investigation.

α/β-Hydrolase Domain (ABHD) Inhibitors as New Potential Therapeutic Options against Lipid-Related Diseases

Much of the experimental evidence in the literature has linked altered lipid metabolism to severe diseases such as cancer, obesity, cardiovascular pathologies, diabetes, and neurodegenerative diseases. Therefore, targeting key effectors of the dysregulated lipid metabolism may represent an effective strategy to counteract these pathological conditions. In this context, α/β-hydrolase domain (ABHD) enzymes represent an important and diversified family of proteins, which are involved in the complex environment of lipid signaling, metabolism, and regulation. Moreover, some members of the ABHD family play an important role in the endocannabinoid system, being designated to terminate the signaling of the key endocannabinoid regulator 2-arachidonoylglycerol. This Perspective summarizes the research progress in the development of ABHD inhibitors and modulators: design strategies, structure-activity relationships, action mechanisms, and biological studies of the main ABHD ligands will be highlighted.

Genome-Wide Association Study of Meat Quality Traits in a Three-Way Crossbred Commercial Pig Population

Meat quality is an important trait for pig-breeding programs aiming to meet consumers' demands. Geneticists must improve meat quality based on their understanding of the underlying genetic mechanisms. Previous studies showed that most meat-quality indicators were low-to-moderate heritability traits; therefore, improving meat quality using conventional techniques remains a challenge. Here, we performed a genome-wide association study of meat-quality traits using the GeneSeek Porcine SNP50K BeadChip in 582 crossbred Duroc × (Landrace × Yorkshire) commercial pigs (249 males and 333 females). Meat conductivity, marbling score, moisture, meat color, pH, and intramuscular fat (IMF) content were investigated. The genome-wide association study was performed using both fixed and random model Circulating Probability Unification (FarmCPU) and a mixed linear model (MLM) with the rMVP software. The genomic heritability of the studied traits ranged from 0.13 ± 0.07 to 0.55 ± 0.08 for conductivity and meat color, respectively. Thirty-two single-nucleotide polymorphisms (SNPs) were identified for meat quality in the crossbred pigs using both FarmCPU and MLM. Among the detected SNPs, five, nine, seven, four, six, and five were significantly associated with conductivity, IMF, marbling score, meat color, moisture, and pH, respectively. Several candidate genes for meat quality were identified in the detected genomic regions. These findings will contribute to the ongoing improvement of meat quality, meeting consumer demands and improving the economic outlook for the swine industry.

Runs of homozygosity provide a genome landscape picture of inbreeding and genetic history of European autochthonous and commercial pig breeds

ROHs are long stretches of DNA homozygous at each polymorphic position. The proportion of genome covered by ROHs and their length are indicators of the level and origin of inbreeding. Frequent common ROHs within the same population define ROH islands and indicate hotspots of selection. In this work, we investigated ROHs in a total of 1131 pigs from 20 European local pig breeds and in three cosmopolitan breeds, genotyped with the GGP Porcine HD Genomic Profiler. plink software was used to identify ROHs. Size classes and genomic inbreeding parameters were evaluated. ROH islands were defined by evaluating different thresholds of homozygous SNP frequency. A functional overview of breed-specific ROH islands was obtained via over-representation analyses of GO biological processes. Mora Romagnola and Turopolje breeds had the largest proportions of genome covered with ROH (~1003 and ~955 Mb respectively), whereas Nero Siciliano and Sarda breeds had the lowest proportions (~207 and 247 Mb respectively). The highest proportion of long ROH (>16 Mb) was in Apulo-Calabrese, Mora Romagnola and Casertana. The largest number of ROH islands was identified in the Italian Landrace (n = 32), Cinta Senese (n = 26) and Lithuanian White Old Type (n = 22) breeds. Several ROH islands were in regions encompassing genes known to affect morphological traits. Comparative ROH structure analysis among breeds indicated the similar genetic structure of local breeds across Europe. This study contributed to understanding of the genetic history of the investigated pig breeds and provided information to manage these pig genetic resources.

TBC1D1 interacting proteins, VPS13A and VPS13C, regulate GLUT4 homeostasis in C2C12 myotubes

Proteins involved in the spaciotemporal regulation of GLUT4 trafficking represent potential therapeutic targets for the treatment of insulin resistance and type 2 diabetes. A key regulator of insulin- and exercise-stimulated glucose uptake and GLUT4 trafficking is TBC1D1. This study aimed to identify proteins that regulate GLUT4 trafficking and homeostasis via TBC1D1. Using an unbiased quantitative proteomics approach, we identified proteins that interact with TBC1D1 in C2C12 myotubes including VPS13A and VPS13C, the Rab binding proteins EHBP1L1 and MICAL1, and the calcium pump SERCA1. These proteins associate with TBC1D1 via its phosphotyrosine binding (PTB) domains and their interactions with TBC1D1 were unaffected by AMPK activation, distinguishing them from the AMPK regulated interaction between TBC1D1 and AMPKα1 complexes. Depletion of VPS13A or VPS13C caused a post-transcriptional increase in cellular GLUT4 protein and enhanced cell surface GLUT4 levels in response to AMPK activation. The phenomenon was specific to GLUT4 because other recycling proteins were unaffected. Our results provide further support for a role of the TBC1D1 PTB domains as a scaffold for a range of Rab regulators, and also the VPS13 family of proteins which have been previously linked to fasting glycaemic traits and insulin resistance in genome wide association studies.

Whole genome re-sequencing reveals recent signatures of selection in three strains of farmed Nile tilapia (Oreochromis niloticus)

Nile tilapia belongs to the second most cultivated group of fish in the world, mainly because of its favorable characteristics for production. Genetic improvement programs and domestication process of Nile tilapia may have modified the genome through selective pressure, leaving signals that can be detected at the molecular level. In this work, signatures of selection were identified using genome-wide SNP data, by two haplotype-based (iHS and Rsb) and one F_ST based method. Whole-genome re-sequencing of 326 individuals from three strains (A, B and C) of farmed tilapia maintained in Brazil and Costa Rica was carried out using Illumina HiSeq 2500 technology. After applying conventional SNP-calling and quality-control filters, ~ 1.3 M high-quality SNPs were inferred and used as input for the iHS, Rsb and F_ST based methods. We detected several candidate genes putatively subjected to selection in each strain. A considerable number of these genes are associated with growth (e.g. NCAPG, KLF3, TBC1D1, TTN), early development (e.g. FGFR3, PFKFB3), and immunity traits (e.g. NLRC3, PIGR, MAP1S). These candidate genes represent putative genomic landmarks that could be associated to traits of biological and commercial interest in farmed Nile tilapia.

Genome-wide association studies for the number of teats and teat asymmetry patterns in Large White pigs

The number of teats is a morphological trait that influences the mothering ability of the sows and thus their reproduction performances. In this study, we carried out GWASs for the total number of teats and other 12 related parameters in 821 Italian Large White heavy pigs. All pigs were genotyped with the Illumina PorcineSNP60 BeadChip array. For four investigated parameters (total number of teats, the number of teats of the left line, the number of teats of the right line and the maximum number of teats comparing the two sides), significant markers were identified on SSC7, in the region of the vertnin (VRTN) gene. Significant markers for the numbers of posterior teats and the absolute difference between anterior and posterior teat numbers were consistently identified on SSC6. The most significant SNP for these parameters was an intron variant in the TOX high mobility group box family member 3 (TOX3) gene. For the other four parameters (absolute difference between the two sides; anterior teats; the ratio between the posterior and the anterior number of teats; and the absence or the presence of extra teats) only suggestively significant markers were identified on several other chromosomes. This study further supported the role of the VRTN gene region in affecting the recorded variability of the number of teats in the Italian Large White pig population and identified a genomic region potentially affecting the biological mechanisms controlling the developmental programme of morphological features in pigs.

Identification of genetic loci associated with growth traits at weaning in yak through a genome-wide association study

A multilocus GWAS was performed to explore the genetic architecture of four growth traits in yak. In total, 354 female yaks for which measurements of body weight (BW), withers height (WH), body length (BL) and chest girth (CG) at weaning were available underwent genotyping with the Illumina BovineHD BeadChip (770K). After quality control, we retained 98 688 SNPs and 354 animals for GWAS analysis. We identified seven, 18, seven and nine SNPs (corresponding to seven, 17, seven and eight candidate genes) associated with BW, WH, BL and CG at weaning respectively. Interestingly, most of these candidate genes were reported to be involved in growth-related processes such as muscle formation, lipid deposition, feed efficiency, carcass composition and development of the central and peripheral nervous system. Our results offer novel insight into the molecular architecture underpinning yak growth traits. Further functional analyses are thus warranted to explore the molecular mechanisms whereby these genes affect these traits of interest.

An association analysis for 14 candidate genes mapping to meat quality quantitative trait loci in a Duroc pig population reveals that the ATP1A2 genotype is highly associated with muscle electric conductivity

In previous GWAS carried out in a Duroc commercial line (Lipgen population), we detected on pig chromosomes 3, 4 and 14 several QTL for gluteus medius muscle redness (GM a*), electric conductivity in the longissimus dorsi muscle (LD CE) and vaccenic acid content in the LD muscle (LD C18:1 n - 7), respectively. We have genotyped, in the Lipgen population, 19 SNPs mapping to 14 genes located within these QTL. Subsequently, association analyses have been performed. After correction for multiple testing, two SNPs in the TGFBRAP1 (rs321173745) and SELENOI (rs330820437) genes were associated with GM a*, whereas ACADSB (rs81449951) and GPR26 (rs343087568) genotypes displayed significant associations with LD vaccenic content. Moreover, the polymorphisms located at the ATP1A2 (rs344748241), ATP8B2 (rs81382410) and CREB3L4 (rs321278469 and rs330133789) genes showed significant associations with LD CE. We made a second round of association analyses including the SNPs mentioned above as well as other SNPs located in the chromosomes to which they map. After performing a correction for multiple testing, the only association that remained significant at the chromosome-wide level was that between the ATP1A2 genotype and LD CE. From a functional point of view, this association is meaningful because this locus encodes a subunit of the Na⁺ /K⁺ -ATPase responsible for maintaining an electrochemical gradient across the plasma membrane.

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.

Sequence analysis and structure prediction of ABHD16A and the roles of the ABHD family members in human disease

Abhydrolase domain containing 16A (ABHD16A) is a member of the α/β hydrolase domain-containing (ABHD) protein family and is expressed in a variety of animal cells. Studies have shown that ABHD16A has acylglycerol lipase and phosphatidylserine lipase activities. Its gene location in the main histocompatibility complex (MHC) III gene cluster suggests that this protein may participate in the immunomodulation of the body. The results of studies investigating nearly 20 species of ABHDs reveal that the ABHD proteins are key factors in metabolic regulation and disease occurrence and development. In this paper, we summarize the related progress regarding the function of ABHD16A and other ABHD proteins. A prediction of the active sites and structural domains of ABHD16A and an analysis of the amino acid sites are included. Moreover, we analysed the amino acid sequences of the ABHD16A molecules in different species and provide an overview of the related functions and diseases associated with these proteins. The functions and diseases related to ABHD are systematically summarized and highlighted. Future research directions for studies investigating the functions and mechanisms of these proteins are also suggested. Further studies investigating the function of ABHD proteins may further confirm their positions as important determinants of lipid metabolism and related diseases.

Isoform-specific AMPK association with TBC1D1 is reduced by a mutation associated with severe obesity

AMP-activated protein kinase (AMPK) is a key regulator of cellular and systemic energy homeostasis which achieves this through the phosphorylation of a myriad of downstream targets. One target is TBC1D1 a Rab-GTPase-activating protein that regulates glucose uptake in muscle cells by integrating insulin signalling with that promoted by muscle contraction. Ser²³⁷ in TBC1D1 is a target for phosphorylation by AMPK, an event which may be important in regulating glucose uptake. Here, we show AMPK heterotrimers containing the α1, but not the α2, isoform of the catalytic subunit form an unusual and stable association with TBC1D1, but not its paralogue AS160. The interaction between the two proteins is direct, involves a dual interaction mechanism employing both phosphotyrosine-binding (PTB) domains of TBC1D1 and is increased by two different pharmacological activators of AMPK (AICAR and A769962). The interaction enhances the efficiency by which AMPK phosphorylates TBC1D1 on its key regulatory site, Ser²³⁷. Furthermore, the interaction is reduced by a naturally occurring R125W mutation in the PTB1 domain of TBC1D1, previously found to be associated with severe familial obesity in females, with a concomitant reduction in Ser²³⁷ phosphorylation. Our observations provide evidence for a functional difference between AMPK α-subunits and extend the repertoire of protein kinases that interact with substrates via stabilisation mechanisms that modify the efficacy of substrate phosphorylation.

A comparative analysis of label-free liquid chromatography-mass spectrometry liver proteomic profiles highlights metabolic differences between pig breeds

The liver is a complex organ governing several physiological processes that define biological mechanisms affecting growth, feed efficiency and performance traits in all livestock species, including pig. Proteomics may contribute to a better understanding of the relationship between liver functions and complex production traits in pigs and to characterize this species as biomedical model. This study applied, for the first time, a label‐free liquid chromatography-mass spectrometry (LC‐MS) proteomic approach to compare the liver proteome profiles of two important heavy pig breeds, Italian Duroc and Italian Large White. Liver specimens were collected (after slaughtering) from performance tested pigs of these two breeds, raised in standard conditions. The label‐free LC‐MS method captured a total of 501 proteins of which 200 were subsequently considered in the between breeds comparison. A statistical pipeline based on the sparse Partial Least Squares Discriminant Analysis (sPLS-DA), coupled with stability and significance tests, was applied for the identification of up or down regulated proteins between breeds. This analysis revealed a total of 25 proteins clearly separating Italian Duroc and Italian Large White pigs. Among the top proteins differentiating the two breeds, 3-ketoacyl-CoA thiolase, mitochondrial (ACAA2) and histone H2B type 2-F (HIST2H2BF) were up-regulated in Italian Duroc pigs and carboxylesterase 3 (CES3) and ketohexokinase (KHK) were up-regulated in Italian Large White pigs. Fatty acid synthase (FASN), involved in fatty acid metabolism and encoded by a gene located in a QTL region for fatty acid composition, was up-regulated in Italian Large White pigs. The in silico protein interaction analysis showed that 16 of these proteins were connected in one big module. Bioinformatic functional analysis indicated that differentially expressed proteins were involved in several biological processes related to the metabolism of lipids, amino-acids, carbohydrates, cofactors and antibiotics/drugs, suggesting that these functions might distinguish Italian Duroc and Italian Large White pigs. This pilot comparative proteomic analysis of the porcine liver highlighted several biological factors that could determine the peculiar production potentials of these two heavy pig breeds, derived by their different genetic backgrounds.

Feed-efficient pigs exhibit molecular patterns allowing a timely circulation of hormones and nutrients

Feed efficiency (FE) is a measure of the rate between feed intake and body weight gain and is subject to constant progress in pigs, based on extensive performance tests and analyses of physiological parameters. However, endocrine regulatory circuits that comprise the sensation and perception of intrinsic requirements and appropriate systemic responses have not yet been fully elucidated. It is hypothesized that the gut-brain axis, which is a network of hierarchical anterior regulatory tissues, contributes largely to variations in FE. Therefore, full-sib pigs with extreme residual feed intake values were assigned to experimental groups of high and low FE. Relevant hormones, minerals, and metabolites including fatty acid profiles were analyzed in serum to assess postprandial conditions. Transcriptome profiles were deduced from intestinal (duodenum, jejunum, ileum) and neuroendocrine tissues (hypothalamus). Serum analyses of feed-efficient animals showed an increased content of the incretin GIP, calcium, magnesium, β-hydroxybutyric acid, and fat compared with low-FE pigs. Complementary expression profiles in intestinal tissues indicate a modulated permeability and host-microbe interaction in FE-divergent pigs. Transcriptomic analyses of the hypothalamus showed that differences between the FE groups in appetite and satiety regulation are less pronounced. However, hypothalamic abundance of transcripts like ADCY7, LHCGR, and SLC2A7 and molecular signatures in local and systemic tissue sites indicate that increased allocation and circulation of energy equivalents, minerals, and hormones are promoted in feed-efficient animals. Overall, patterns of gastrointestinal hormones and gene expression profiles identified host-microbiota interaction, intestinal permeability, feed intake regulation, and energy expenditure as potential mechanisms affecting FE in pigs.

Regulation of RabGAPs involved in insulin action

Rab (Ras-related proteins in brain) GTPases are key proteins responsible for a multiplicity of cellular trafficking processes. Belonging to the family of monomeric GTPases, they are regulated by cycling between their active GTP-bound and inactive GDP-bound conformations. Despite possessing a slow intrinsic GTP hydrolysis activity, Rab proteins rely on RabGAPs (Rab GTPase-activating proteins) that catalyze GTP hydrolysis and consequently inactivate the respective Rab GTPases. Two related RabGAPs, TBC1D1 and TBC1D4 (=AS160) have been described to be associated with obesity-related traits and type 2 diabetes in both mice and humans. Inactivating mutations of TBC1D1 and TBC1D4 lead to substantial changes in trafficking and subcellular distribution of the insulin-responsive glucose transporter GLUT4, and to subsequent alterations in energy substrate metabolism. The activity of the RabGAPs is controlled through complex phosphorylation events mediated by protein kinases including AKT and AMPK, and by putative regulatory interaction partners. However, the dynamics and downstream events following phosphorylation are not well understood. This review focuses on the specific role and regulation of TBC1D1 and TBC1D4 in insulin action.

Genome-wide scan reveals population stratification and footprints of recent selection in Nelore cattle

Background

This study aimed at (1) assessing the genomic stratification of experimental lines of Nelore cattle that have experienced different selection regimes for growth traits, and (2) identifying genomic regions that have undergone recent selection. We used a sample of 763 animals genotyped with the Illumina BovineHD BeadChip, among which 674 animals originated from two lines that are maintained under directional selection for increased yearling body weight and 89 animals from a control line that is maintained under stabilizing selection.

Results

Multidimensional analysis of the genomic dissimilarity matrix and admixture analysis revealed a substantial level of population stratification between the directional selection lines and the stabilizing selection control line. Two of the three tests used to detect selection signatures (F_ST, XP-EHH and iHS) revealed six candidate regions with indications of selection, which strongly indicates truly positive signals. The set of identified candidate genes included several genes with roles that are functionally related to growth metabolism, such as COL14A1, CPT1C, CRH, TBC1D1, and XKR4.

Conclusions

The current study identified genetic stratification that resulted from almost four decades of divergent selection in an experimental Nelore population, and highlighted autosomal genomic regions that present patterns of recent selection. Our findings provide a basis for a better understanding of the metabolic mechanism that underlies the growth traits, which are modified by selection for yearling body weight.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0381-2) contains supplementary material, which is available to authorized users.

Association of single nucleotide polymorphisms in fat metabolism candidate genes with fatty acid profiles of muscle and subcutaneous fat in heavy pigs

Dietary and organoleptic qualities of pork products are largely influenced by the profiles of fatty acids (FAs) of meat and fat. The objective of this work was to investigate the potential associations between FA profile in subcutaneous adipose (back fat and leg fat) and muscular (longissimus dorsi and biceps femoris) tissues in heavy pigs (n = 129, 9-month-old) with single nucleotide polymorphisms (SNPs) in six candidate genes involved in fat metabolism: Stearoyl-CoA desaturase (SCD), Diacylglycerol acyltransferase 1 and 2 (DGAT1 &DGAT2), Microsomal Triglyceride Transfer Protein (MTTP), Fatty Acid Synthase (FASN) and Heart- fatty acid binding protein (H-FABP). Preliminary results suggest a putative association between MTTP, DGAT2 and FASN and the FA content in both fat and meat, whereas between DGAT1, SCD and H-FABP the association was found in adipose tissue only. However, the effect of the analyzed genes, needs to be verified in a larger and better characterized pig population to support the hypothesized associations with FA content.

Genome-wide association study for feed efficiency and growth traits in U.S. beef cattle

Background

Single nucleotide polymorphism (SNP) arrays for domestic cattle have catalyzed the identification of genetic markers associated with complex traits for inclusion in modern breeding and selection programs. Using actual and imputed Illumina 778K genotypes for 3887 U.S. beef cattle from 3 populations (Angus, Hereford, SimAngus), we performed genome-wide association analyses for feed efficiency and growth traits including average daily gain (ADG), dry matter intake (DMI), mid-test metabolic weight (MMWT), and residual feed intake (RFI), with marker-based heritability estimates produced for all traits and populations.

Results

Moderate and/or large-effect QTL were detected for all traits in all populations, as jointly defined by the estimated proportion of variance explained (PVE) by marker effects (PVE ≥ 1.0%) and a nominal P-value threshold (P ≤ 5e-05). Lead SNPs with PVE ≥ 2.0% were considered putative evidence of large-effect QTL (n = 52), whereas those with PVE ≥ 1.0% but < 2.0% were considered putative evidence for moderate-effect QTL (n = 35). Identical or proximal lead SNPs associated with ADG, DMI, MMWT, and RFI collectively supported the potential for either pleiotropic QTL, or independent but proximal causal mutations for multiple traits within and between the analyzed populations. Marker-based heritability estimates for all investigated traits ranged from 0.18 to 0.60 using 778K genotypes, or from 0.17 to 0.57 using 50K genotypes (reduced from Illumina 778K HD to Illumina Bovine SNP50). An investigation to determine if QTL detected by 778K analysis could also be detected using 50K genotypes produced variable results, suggesting that 50K analyses were generally insufficient for QTL detection in these populations, and that relevant breeding or selection programs should be based on higher density analyses (imputed or directly ascertained).

Conclusions

Fourteen moderate to large-effect QTL regions which ranged from being physically proximal (lead SNPs ≤ 3Mb) to fully overlapping for RFI, DMI, ADG, and MMWT were detected within and between populations, and included evidence for pleiotropy, proximal but independent causal mutations, and multi-breed QTL. Bovine positional candidate genes for these traits were functionally conserved across vertebrate species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3754-y) contains supplementary material, which is available to authorized users.

Genetic Marker Discovery in Complex Traits: A Field Example on Fat Content and Composition in Pigs

Among the large number of attributes that define pork quality, fat content and composition have attracted the attention of breeders in the recent years due to their interaction with human health and technological and sensorial properties of meat. In livestock species, fat accumulates in different depots following a temporal pattern that is also recognized in humans. Intramuscular fat deposition rate and fatty acid composition change with life. Despite indication that it might be possible to select for intramuscular fat without affecting other fat depots, to date only one depot-specific genetic marker (PCK1 c.2456C>A) has been reported. In contrast, identification of polymorphisms related to fat composition has been more successful. For instance, our group has described a variant in the stearoyl-coA desaturase (SCD) gene that improves the desaturation index of fat without affecting overall fatness or growth. Identification of mutations in candidate genes can be a tedious and costly process. Genome-wide association studies can help in narrowing down the number of candidate genes by highlighting those which contribute most to the genetic variation of the trait. Results from our group and others indicate that fat content and composition are highly polygenic and that very few genes explain more than 5% of the variance of the trait. Moreover, as the complexity of the genome emerges, the role of non-coding genes and regulatory elements cannot be disregarded. Prediction of breeding values from genomic data is discussed in comparison with conventional best linear predictors of breeding values. An example based on real data is given, and the implications in phenotype prediction are discussed in detail. The benefits and limitations of using large SNP sets versus a few very informative markers as predictors of genetic merit of breeding candidates are evaluated using field data as an example.

Association study between single nucleotide polymorphisms in porcine genes and pork quality traits for fresh consumption and processing into Italian dry-cured ham

Single nucleotide polymorphisms (SNPs) of six genes (TTN, PRKAG3, CAST, CTSB, CTSF, and MYPN), known for associations with carcass and meat quality traits, post mortem proteolysis, were screened in a commercial crossed population of 368 heavy pigs (Large White x Landrace)×Duroc, reared according to the rules of Italian Protected Designation of Origin (PDO) dry-cured ham. Carcass, longissimus thoracis et lumborum muscle (LTL), and green ham traits were obtained after slaughtering, main weight losses of dry-cured hams were collected during processing. The results showed the impact of CAST variants on carcass weight, of CTSF on LTL tenderness, ham weight and fatness, of PRKAG3 and TTN on ultimate pH, hamweight. This study, while confirming significant associations between SNPs of genes and qualitative traits of carcass, longissimus and ham, supports CTSF as candidate gene suitable for fresh consumption purpose (tenderness of longissimus at 24h post mortem), and for dry-cured ham processing (higher thickness of ham subcutaneous fat).

Deconstructing the pig sex metabolome: Targeted metabolomics in heavy pigs revealed sexual dimorphisms in plasma biomarkers and metabolic pathways

Metabolomics has opened new possibilities to investigate metabolic differences among animals. In this study, we applied a targeted metabolomic approach to deconstruct the pig sex metabolome as defined by castrated males and entire gilts. Plasma from 545 performance-tested Italian Large White pigs (172 castrated males and 373 females) sampled at about 160 kg live weight were analyzed for 186 metabolites using the Biocrates AbsoluteIDQ p180 Kit. After filtering, 132 metabolites (20 AA, 11 biogenic amines, 1 hexose, 13 acylcarnitines, 11 sphingomyelins, 67 phosphatidylcholines, and 9 lysophosphatidylcholines) were retained for further analyses. The multivariate approach of the sparse partial least squares discriminant analysis was applied, together with a specifically designed statistical pipeline, that included a permutation test and a 10 cross-fold validation procedure that produced stability and effect size statistics for each metabolite. Using this approach, we identified 85 biomarkers (with metabolites from all analyzed chemical families) that contributed to the differences between the 2 groups of pigs ( < 0.05 at the stability statistic test). All acylcarnitines and almost all biogenic amines were higher in castrated males than in gilts. Metabolites involved in tryptophan catabolism had the largest differences (i.e., delta = 20% for serotonin) between castrated males (higher) and gilts (lower). The level of several AA (Ala, Arg, Gly, His, Lys, Ser, Thr, and Trp) was higher in gilts (delta was from approximately 1.0 to approximately 4.8%) whereas products of AA catabolism (taurine, 2-aminoadipic acid, and methionine sulfoxide) were higher in castrated males (delta was approximately 5.0-6.0%), suggesting a metabolic shift in castrated males toward energy storage and lipid production. Similar general patterns were observed for most sphingomyelins, phosphatidylcholines, and lysophosphatidylcholines. Metabolomic pathway analysis and pathway enrichment identified several differences between the 2 sexes. This metabolomic overview opened new clues on the biochemical mechanisms underlying sexual dimorphism that, on one hand, might explain differences in terms of economic traits between castrated male pigs and entire gilts and, on the other hand, could strengthen the pig as a model to define metabolic mechanisms related to fat deposition.

Next generation semiconductor based sequencing of bitter taste receptor genes in different pig populations and association analysis using a selective DNA pool-seq approach

Taste perception in animals affects feed intake and may influence production traits. In particular, bitter is sensed by receptors encoded by the family of TAS2R genes. In this research, using a DNA pool-seq approach coupled with next generation semiconductor based target resequencing, we analysed nine porcine TAS2R genes (TAS2R1, TAS2R3, TAS2R4, TAS2R7, TAS2R9, TAS2R10, TAS2R16, TAS2R38 and TAS2R39) to identify variability and, at the same time, estimate single nucleotide polymorphism (SNP) allele frequencies in several populations and testing differences in an association analysis. Equimolar DNA pools were prepared for five pig breeds (Italian Duroc, Italian Landrace, Pietrain, Meishan and Casertana) and wild boars (5-10 individuals each) and for two groups of Italian Large White pigs with extreme and divergent back fat thickness (50 + 50 pigs). About 1.8 million reads were obtained by sequencing amplicons generated from these pools. A total of 125 SNPs were identified, of which 37 were missense mutations. Three of them (p.Ile53Phe and p.Trp85Leu in TAS2R4; p.Leu37Ser in TAS2R39) could have important effects on the function of these bitter taste receptors, based on in silico predictions. Variability in wild boars seems lower than that in domestic breeds potentially as a result of selective pressure in the wild towards defensive bitter taste perception. Three SNPs in TAS2R38 and TAS2R39 were significantly associated with back fat thickness. These results may be important to understand the complexity of taste perception and their associated effects that could be useful to develop nutrigenetic approaches in pig breeding and nutrition.

Analysis of the porcine APOA2 gene expression in liver, polymorphism identification and association with fatty acid composition traits

APOA2 is a protein implicated in triglyceride, fatty acid and glucose metabolism. In pigs, the APOA2 gene is located on pig chromosome 4 (SSC4) in a QTL region affecting fatty acid composition, fatness and growth traits. In this study, we evaluated APOA2 as a candidate gene for meat quality traits in an Iberian × Landrace backcross population. The APOA2:c.131T>A polymorphism, located in exon 3 of APOA2 and determining a missense mutation, was associated with the percentage of hexadecenoic acid [C16:1(n-9)], linoleic acid [C18:2(n-6)], α-linolenic acid [C18:3(n-3)], dihomo-gamma-linolenic acid [C20:3(n-6)] and polyunsaturated fatty acids (PUFAs) in backfat. Furthermore, this SNP was associated with the global mRNA expression levels of APOA2 in liver and was used as a marker to determine allelic expression imbalance by pyrosequencing. We determined an overexpression of the T allele in heterozygous samples with a mean ratio of 2.8 (T/A), observing a high variability in the allelic expression among individuals. This result suggests that complex regulatory mechanisms, beyond a single polymorphism (e.g. epigenetic effects or multiple cis-acting polymorphisms), may be regulating APOA2 gene expression.

Dynamic gene expression profiles during postnatal development of porcine subcutaneous adipose

A better understanding of the control of lipogenesis is of critical importance for both human and animal physiology. This requires a better knowledge of the changes of gene expression during the process of adipose tissue development. Thus, the objective of the current study was to determine the effects of development on subcutaneous adipose tissue gene expression in growing and adult pigs. Here, we present a comprehensive investigation of mRNA transcriptomes in porcine subcutaneous adipose tissue across four developmental stages using digital gene expression profiling. We identified 3,274 differential expressed genes associated with oxidative stress, immune processes, apoptosis, energy metabolism, insulin stimulus, cell cycle, angiogenesis and translation. A set of universally abundant genes (ATP8, COX2, COX3, ND1, ND2, SCD and TUBA1B) was found across all four developmental stages. This set of genes may play important roles in lipogenesis and development. We also identified development-related gene expression patterns that are linked to the different adipose phenotypes. We showed that genes enriched in significantly up-regulated profiles were associated with phosphorylation and angiogenesis. In contrast, genes enriched in significantly down-regulated profiles were related to cell cycle and cytoskeleton organization, suggesting an important role for these biological processes in adipose growth and development. These results provide a resource for studying adipose development and promote the pig as a model organism for researching the development of human obesity, as well as being used in the pig industry.

Genetics of Adiposity in Large Animal Models for Human Obesity-Studies on Pigs and Dogs

The role of domestic mammals in the development of human biomedical sciences has been widely documented. Among these model species the pig and dog are of special importance. Both are useful for studies on the etiology of human obesity. Genome sequences of both species are known and advanced genetic tools [eg, microarray SNP for genome wide association studies (GWAS), next generation sequencing (NGS), etc.] are commonly used in such studies. In the domestic pig the accumulation of adipose tissue is an important trait, which influences meat quality and fattening efficiency. Numerous quantitative trait loci (QTLs) for pig fatness traits were identified, while gene polymorphisms associated with these traits were also described. The situation is different in dog population. Generally, excessive accumulation of adipose tissue is considered, similar to humans, as a complex disease. However, research on the genetic background of canine obesity is still in its infancy. Between-breed differences in terms of adipose tissue accumulation are well known in both animal species. In this review we show recent advances of studies on adipose tissue accumulation in pigs and dogs, and their potential importance for studies on human obesity.

Rab GAPs AS160 and Tbc1d1 play nonredundant roles in the regulation of glucose and energy homeostasis in mice

The related Rab GTPase-activating proteins (Rab GAPs) AS160 and Tbc1d1 regulate the trafficking of the glucose transporter GLUT4 that controls glucose uptake in muscle and fat cells and glucose homeostasis. AS160- and Tbc1d1-deficient mice exhibit different adipocyte- and skeletal muscle-specific defects in glucose uptake, GLUT4 expression and trafficking, and glucose homeostasis. A recent study analyzed male mice with simultaneous deletion of AS160 and Tbc1d1 (AS160(-/-)/Tbc1d1(-/-) mice). Herein, we describe abnormalities in male and female AS160(-/-)/Tbc1d1(-/-) mice on another strain background. We confirm the earlier observation that GLUT4 expression and glucose uptake defects of single-knockout mice join in AS160(-/-)/Tbc1d1(-/-) mice to affect all skeletal muscle and adipose tissues. In large mixed fiber-type skeletal muscles, changes in relative basal GLUT4 plasma membrane association in AS160(-/-) and Tbc1d1(-/-) mice also combine in AS160(-/-)/Tbc1d1(-/-) mice. However, we found different glucose uptake abnormalities in isolated skeletal muscles and adipocytes than reported previously, resulting in different interpretations of how AS160 and Tbc1d1 regulate GLUT4 translocation to the cell surface. In support of a larger role for AS160 in glucose homeostasis, in contrast with the previous study, we find similarly impaired glucose and insulin tolerance in AS160(-/-)/Tbc1d1(-/-) and AS160(-/-) mice. However, in vivo glucose uptake abnormalities in AS160(-/-)/Tbc1d1(-/-) skeletal muscles differ from those observed previously in AS160(-/-) mice, indicating additional defects due to Tbc1d1 deletion. Similar to AS160- and Tbc1d1-deficient mice, AS160(-/-)/Tbc1d1(-/-) mice show sex-specific abnormalities in glucose and energy homeostasis. In conclusion, our study supports nonredundant functions for AS160 and Tbc1d1.

Twenty years of artificial directional selection have shaped the genome of the Italian Large White pig breed

In this study, we investigated at the genome-wide level if 20 years of artificial directional selection based on boar genetic evaluation obtained with a classical BLUP animal model shaped the genome of the Italian Large White pig breed. The most influential boars of this breed (n = 192), born from 1992 (the beginning of the selection program of this breed) to 2012, with an estimated breeding value reliability of >0.85, were genotyped with the Illumina Porcine SNP60 BeadChip. After grouping the boars in eight classes according to their year of birth, filtered single nucleotide polymorphisms (SNPs) were used to evaluate the effects of time on genotype frequency changes using multinomial logistic regression models. Of these markers, 493 had a PBonferroni  < 0.10. However, there was an increasing number of SNPs with a decreasing level of allele frequency changes over time, representing a continuous profile across the genome. The largest proportion of the 493 SNPs was on porcine chromosome (SSC) 7, SSC2, SSC8 and SSC18 for a total of 204 haploblocks. Functional annotations of genomic regions, including the 493 shifted SNPs, reported a few Gene Ontology terms that might underly the biological processes that contributed to increase performances of the pigs over the 20 years of the selection program. The obtained results indicated that the genome of the Italian Large White pigs was shaped by a directional selection program derived by the application of methodologies assuming the infinitesimal model that captured a continuous trend of allele frequency changes in the boar population.

Deletion of the Rab GAP Tbc1d1 modifies glucose, lipid, and energy homeostasis in mice

Tbc1d1 is a Rab GTPase-activating protein (GAP) implicated in regulating intracellular retention and cell surface localization of the glucose transporter GLUT4 and thus glucose uptake in a phosphorylation-dependent manner. Tbc1d1 is most abundant in skeletal muscle but is expressed at varying levels among different skeletal muscles. Previous studies with male Tbc1d1-deficient (Tbc1d1(-/-)) mice on standard and high-fat diets established a role for Tbc1d1 in glucose, lipid, and energy homeostasis. Here we describe similar, but also additional abnormalities in male and female Tbc1d1(-/-) mice. We corroborate that Tbc1d1 loss leads to skeletal muscle-specific and skeletal muscle type-dependent abnormalities in GLUT4 expression and glucose uptake in female and male mice. Using subcellular fractionation, we show that Tbc1d1 controls basal intracellular GLUT4 retention in large skeletal muscles. However, cell surface labeling of extensor digitorum longus muscle indicates that Tbc1d1 does not regulate basal GLUT4 cell surface exposure as previously suggested. Consistent with earlier observations, female and male Tbc1d1(-/-) mice demonstrate increased energy expenditure and skeletal muscle fatty acid oxidation. Interestingly, we observe sex-dependent differences in in vivo phenotypes. Female, but not male, Tbc1d1(-/-) mice have decreased body weight and impaired glucose and insulin tolerance, but only male Tbc1d1(-/-) mice show increased lipid clearance after oil gavage. We surmise that similar changes at the tissue level cause differences in whole-body metabolism between male and female Tbc1d1(-/-) mice and between male Tbc1d1(-/-) mice in different studies due to variations in body composition and nutrient handling.

Polymorphism of fat metabolism genes as candidate markers for meat quality and production traits in heavy pigs

High meat quality is required for dry-cured ham production, which quality depends on meat fat quantity and composition. The aim was to study the polymorphisms of six genes involved in fat metabolism, namely, Stearoyl-CoA desaturase (SCD), Diacylglycerol acyltransferase 1 and 2 (DGAT1 &DGAT2), Microsomal triglyceride transfer protein (MTTP), Fatty acid synthase (FASN) and Heart fatty acid binding protein (H-FABP) in two traditional (Large White or Duroc × (Landrace × Large White)) and two industrial hybrids (Goland and Danbred), which are used for dry-cured ham production. Significant associations of SCD and MTTP were found with carcass weight. DGAT2 was associated with back fat thickness and L* fat colour (objective white colour score). Several genes (DGAT2, MTTP &FASN) were associated with weight loss during salting, first step in dry-cured ham production, affecting final yield. Finally, MTTP was associated with shear force. Our findings suggest that the SCD, DGAT2, MTTP and FASN polymorphisms are associated with quality of heavy pig meat products.

Reduced Representation Libraries from DNA Pools Analysed with Next Generation Semiconductor Based-Sequencing to Identify SNPs in Extreme and Divergent Pigs for Back Fat Thickness

The aim of this study was to identify single nucleotide polymorphisms (SNPs) that could be associated with back fat thickness (BFT) in pigs. To achieve this goal, we evaluated the potential and limits of an experimental design that combined several methodologies. DNA samples from two groups of Italian Large White pigs with divergent estimating breeding value (EBV) for BFT were separately pooled and sequenced, after preparation of reduced representation libraries (RRLs), on the Ion Torrent technology. Taking advantage from SNAPE for SNPs calling in sequenced DNA pools, 39,165 SNPs were identified; 1/4 of them were novel variants not reported in dbSNP. Combining sequencing data with Illumina PorcineSNP60 BeadChip genotyping results on the same animals, 661 genomic positions overlapped with a good approximation of minor allele frequency estimation. A total of 54 SNPs showing enriched alleles in one or in the other RRLs might be potential markers associated with BFT. Some of these SNPs were close to genes involved in obesity related phenotypes.

A retrospective analysis of allele frequency changes of major genes during 20 years of selection in the Italian Large White pig breed

In this study, we investigated whether a selection programme based on boar genetic evaluation obtained with a classical BLUP animal model can change allele frequencies in a pig population. All Italian Large White boars born from 1992 to 2012 with estimated breeding value reliability >0.85 (n = 200) were selected among all boars of this breed. Boars were genotyped with markers in major genes (IGF2 intron3-g.3072G>A, MC4R p.D298N, VRTN PRE1 insertion, PRKAG3 p.I199V and FTO g.276T>G). Genotyping data were analysed grouping boars in eight classes according to their year of birth. To evaluate the influence of time on allele frequencies of the genotyped markers, multinomial logistic regression models were computed. Four of five polymorphic sites (IGF2, MC4R, VRTN and FTO) showed significant (p < 0.01) changes in allele frequencies over time due to a progressive and continuous increase of one allele (associated with higher lean meat content, higher average daily gain and favourable feed: gain ratio) and, consequently, decrease of the other one, following the directional selection of the selection programme of this pig breed. The retrospective analysis that was carried out in Italian Large White boars suggests that selection based on methodologies assuming the infinitesimal model is able to modify in a quite short period of time allele frequencies in major genes, increasing the frequency of alleles explaining a relevant (non-infinitesimal) fraction of the overall genetic variability for production traits.

Genome-wide QTL mapping of nine body composition and bone mineral density traits in pigs

Background

Since the pig is one of the most important livestock animals worldwide, mapping loci that are associated with economically important traits and/or traits that influence animal welfare is extremely relevant for efficient future pig breeding. Therefore, the purpose of this study was a genome-wide mapping of quantitative trait loci (QTL) associated with nine body composition and bone mineral traits: absolute (Fat, Lean) and percentage (FatPC, LeanPC) fat and lean mass, live weight (Weight), soft tissue X-ray attenuation coefficient (R), absolute (BMC) and percentage (BMCPC) bone mineral content and bone mineral density (BMD).

Methods

Data on the nine traits investigated were obtained by Dual-energy X-ray absorptiometry for 551 pigs that were between 160 and 200 days old. In addition, all pigs were genotyped using Illumina’s PorcineSNP60 Genotyping BeadChip. Based on these data, a genome-wide combined linkage and linkage disequilibrium analysis was conducted. Thus, we used 44 611 sliding windows that each consisted of 20 adjacent single nucleotide polymorphisms (SNPs). For the middle of each sliding window a variance component analysis was carried out using ASReml. The underlying mixed linear model included random QTL and polygenic effects, with fixed effects of sex, housing, season and age.

Results

Using a Bonferroni-corrected genome-wide significance threshold of P < 0.001, significant peaks were identified for all traits except BMCPC. Overall, we identified 72 QTL on 16 chromosomes, of which 24 were significantly associated with one trait only and the remaining with more than one trait. For example, a QTL on chromosome 2 included the highest peak across the genome for four traits (Fat, FatPC, LeanPC and R). The nearby gene, ZNF608, is known to be associated with body mass index in humans and involved in starvation in Drosophila, which makes it an extremely good candidate gene for this QTL.

Conclusions

Our QTL mapping approach identified 72 QTL, some of which confirmed results of previous studies in pigs. However, we also detected significant associations that have not been published before and were able to identify a number of new and promising candidate genes, such as ZNF608.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0068-2) contains supplementary material, which is available to authorized users.

Biochemical and pharmacological characterization of the human lymphocyte antigen B-associated transcript 5 (BAT5/ABHD16A)

Background

Human lymphocyte antigen B-associated transcript 5 (BAT5, also known as ABHD16A) is a poorly characterized 63 kDa protein belonging to the α/β-hydrolase domain (ABHD) containing family of metabolic serine hydrolases. Its natural substrates and biochemical properties are unknown.

Methodology/Principal Findings

Amino acid sequence comparison between seven mammalian BAT5 orthologs revealed that the overall primary structure was highly (≥95%) conserved. Activity-based protein profiling (ABPP) confirmed successful generation of catalytically active human (h) and mouse (m) BAT5 in HEK293 cells, enabling further biochemical characterization. A sensitive fluorescent glycerol assay reported hBAT5-mediated hydrolysis of medium-chain saturated (C14∶0), long-chain unsaturated (C18∶1, C18∶2, C20∶4) monoacylglycerols (MAGs) and 15-deoxy-Δ^12,14-prostaglandin J₂-2-glycerol ester (15d-PGJ₂-G). In contrast, hBAT5 possessed only marginal diacylglycerol (DAG), triacylglycerol (TAG), or lysophospholipase activity. The best MAG substrates were 1-linoleylglycerol (1-LG) and 15d-PGJ₂-G, both exhibiting low-micromolar K_m values. BAT5 had a neutral pH optimum and showed preference for the 1(3)- vs. 2-isomers of MAGs C18∶1, C18∶2 and C20∶4. Inhibitor profiling revealed that β-lactone-based lipase inhibitors were nanomolar inhibitors of hBAT5 activity (palmostatin B > tetrahydrolipstatin > ebelactone A). Moreover, the hormone-sensitive lipase inhibitor C7600 (5-methoxy-3-(4-phenoxyphenyl)-3H-[1], [3], [4]oxadiazol-2-one) was identified as a highly potent inhibitor (IC₅₀ 8.3 nM). Phenyl and benzyl substituted analogs of C7600 with increased BAT5 selectivity were synthesized and a preliminary SAR analysis was conducted to obtain initial insights into the active site dimensions.

Conclusions/Significance

This study provides an initial characterization of BAT5 activity, unveiling the biochemical and pharmacological properties with in vitro substrate preferences and inhibitor profiles. Utilization of glycerolipid substrates and sensitivity to lipase inhibitors suggest that BAT5 is a genuine lipase with preference for long-chain unsaturated MAGs and could in this capacity regulate glycerolipid metabolism in vivo as well. This preliminary SAR data should pave the way towards increasingly potent and BAT5-selective inhibitors.

Identification of genomic regions associated with feed efficiency in Nelore cattle

Background

Feed efficiency is jointly determined by productivity and feed requirements, both of which are economically relevant traits in beef cattle production systems. The objective of this study was to identify genes/QTLs associated with components of feed efficiency in Nelore cattle using Illumina BovineHD BeadChip (770 k SNP) genotypes from 593 Nelore steers. The traits analyzed included: average daily gain (ADG), dry matter intake (DMI), feed-conversion ratio (FCR), feed efficiency (FE), residual feed intake (RFI), maintenance efficiency (ME), efficiency of gain (EG), partial efficiency of growth (PEG) and relative growth rate (RGR). The Bayes B analysis was completed with Gensel software parameterized to fit fewer markers than animals. Genomic windows containing all the SNP loci in each 1 Mb that accounted for more than 1.0% of genetic variance were considered as QTL region. Candidate genes within windows that explained more than 1% of genetic variance were selected by putative function based on DAVID and Gene Ontology.

Results

Thirty-six QTL (1-Mb SNP window) were identified on chromosomes 1, 2, 3, 5, 6, 7, 8, 9, 10, 12, 14, 15, 16, 18, 19, 20, 21, 22, 24, 25 and 26 (UMD 3.1). The amount of genetic variance explained by individual QTL windows for feed efficiency traits ranged from 0.5% to 9.07%. Some of these QTL minimally overlapped with previously reported feed efficiency QTL for Bos taurus. The QTL regions described in this study harbor genes with biological functions related to metabolic processes, lipid and protein metabolism, generation of energy and growth. Among the positional candidate genes selected for feed efficiency are: HRH4, ALDH7A1, APOA2, LIN7C, CXADR, ADAM12 and MAP7.

Conclusions

Some genomic regions and some positional candidate genes reported in this study have not been previously reported for feed efficiency traits in Bos indicus. Comparison with published results indicates that different QTLs and genes may be involved in the control of feed efficiency traits in this Nelore cattle population, as compared to Bos taurus cattle.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-014-0100-0) contains supplementary material, which is available to authorized users.

Copy number variants in Italian Large White pigs detected using high-density single nucleotide polymorphisms and their association with back fat thickness

The aim of this study was to identify copy number variants (CNVs) in Italian Large White pigs and test them for association with back fat thickness (BFT). Within a population of 12 000 performance-tested pigs, two groups of animals with extreme and divergent BFT estimated breeding values (EBVs; 147 with negative and 150 with positive EBVs) were genotyped with the Illumina Porcine SNP60 BeadChip. CNVs were detected with PENNCNV software. We identified a total of 4146 CNV events in 170 copy number variation regions (CNVRs) located on 15 porcine autosomes. Validation of detected CNVRs was carried out (i) by comparing CNVRs already detected by other studies and (ii) by semiquantitative fluorescent multiplex (SQFM) PCR of a few CNVRs. Most of CNVRs detected in Italian Large White pigs (71.2%) were already reported in other pig breeds/populations, and 82.1% of the CNV events detected by PENNCNV were confirmed by SQFM PCR. For each CNVR, we compared the occurrence of CNV events between the pigs of the high and low BFT EBV tails. Sixteen regions showed significance at P < 0.10, and seven were significant at P < 0.05 but were not significant after Bonferroni correction (Fisher's exact test). These results indicated that CNVs could explain a limited fraction of the genetic variability of fat deposition in Italian Large White pigs. However, it was interesting to note that one of these CNVRs encompassed the ZPLD1 gene. In humans, a rare CNV event including this gene is associated with obesity. Studies identifying CNVs in pigs could assist in elucidating the genetic mechanisms underlying human obesity.

A genomewide association study for average daily gain in Italian Large White pigs

Average daily gain is an important target trait in pig breeding programs. In this study we performed a genomewide association study for ADG in Italian Large White pigs using a selective genotyping approach. Two extreme and divergent groups of Italian Large White pigs (number 190 + 190) were selected among a population of about 10,000 performance tested gilts (EBV for ADG in the 2 groups were -30 ± 14 g and 81 ± 12 g, respectively) and genotyped with the Illumina PorcineSNP60 BeadChip. Association analysis was performed treating the pigs of the 2 extreme groups as cases and controls after correction for family-based stratification. A total of 127 SNP resulted significantly associated with ADG (P nominal value [P(raw)] < 2.0 × 10(-7), P < 0.01 Bonferroni corrected [P(Bonferroni)] < 0.01, false discovery rate < 7.76 × 10(-5)). Another 102 SNP were suggestively associated with the target trait (P(raw) between 2.0 × 10(-7) and 2.02 × 10(-6), P(Bonferroni) < 0.10, false discovery rate < 4.19 × 10(-4)). These SNP were located on all autosomes and on porcine chromosome (SSC) X. The largest number of SNP within this list was on SSC5 (n = 42), SSC7 (34), SSC6 (30), SSC4 (23), and SSC16 (16). These chromosomes were richer in significant or suggestively significant markers than expected (P < 0.001). A quite high number of these SNP (n = 23) were associated with backfat thickness in a previous genomewide association study performed in the same pig population, confirming the negative correlation between the 2 traits. Two or more SNP targeted the same gene: IGSF3 and HS2ST1 (SSC4), OTOGL (SSC5), FTO region (SSC6), and MYLK4 and MCUR1 (SSC7). Other regions that were associated with ADG in previous candidate gene studies (e.g., MC4R on SSC1, IGF2 and LDHA on SSC2, MUC4 on SSC13) 1) included markers with P(raw) < 0.01 that, however, did not pass the stringent threshold of significance adopted in this study or 2) could not be tested because not assigned to the Sscrofa10.2 genome version. Functional annotation of the significant regions using Gene Ontology suggested that many and complex processes at different levels are involved in affecting ADG, indicating the complexity of the genetic factors controlling this ultimate phenotype. The obtained results may contribute to understand the genetic mechanisms determining ADG that could open new perspectives to improve selection efficiency in this breed.

Galaxy tools to study genome diversity

Background

Intra-species genetic variation can be used to investigate population structure, selection, and gene flow in non-model vertebrates; and due to the plummeting costs for genome sequencing, it is now possible for small labs to obtain full-genome variation data from their species of interest. However, those labs may not have easy access to, and familiarity with, computational tools to analyze those data.

Results

We have created a suite of tools for the Galaxy web server aimed at handling nucleotide and amino-acid polymorphisms discovered by full-genome sequencing of several individuals of the same species, or using a SNP genotyping microarray. In addition to providing user-friendly tools, a main goal is to make published analyses reproducible. While most of the examples discussed in this paper deal with nuclear-genome diversity in non-human vertebrates, we also illustrate the application of the tools to fungal genomes, human biomedical data, and mitochondrial sequences.

Conclusions

This project illustrates that a small group can design, implement, test, document, and distribute a Galaxy tool collection to meet the needs of a particular community of biologists.

Identification of immune-related SNPs in the transcriptome of Mytilus chilensis through high-throughput sequencing

Single nucleotide polymorphisms (SNPs) identified in coding regions represent a useful tool for understanding the immune response against pathogens and stressful environmental conditions. In this study, a SNPs database was generated from transcripts involved in the innate immune response of the mussel Mytilus chilensis. The SNPs were identified through hemocytes transcriptome sequencing from 18 individuals, and SNPs mining was performed in 225,336 contigs, yielding 20,306 polymorphisms associated to immune-related genes. Classification of identified SNPs was based on different pathways of the immune response for Mytilus sp. A total of 28 SNPs were identified in the Toll-like receptor pathway and included 5 non-synonymous polymorphisms; 19 SNPs were identified in the apoptosis pathway and included 3 non-synonymous polymorphisms; 35 SNPs were identified in the Ubiquitin-mediated proteolysis pathway and included 4 non-synonymous variants; and 54 SNPs involved in other molecular functions related to the immune response, such as molecular chaperones, antimicrobial peptides, and genes that interacts with marine toxins were also identified. The molecular markers identified in this work could be useful for novel studies, such as those related to associations between high-resolution molecular markers and functional response to pathogen agents.