Genetic determinants for intramuscular fat content and water-holding capacity in mice selected for high muscle mass

Identification of potential biomarkers associated with meat tenderness in Hanwoo (Korean cattle): An expression quantitative trait loci analysis

Meat tenderness is considered the most important trait contributing to beef quality, level of consumer satisfaction, willingness to pay premium prices and industry profit. Genomic selection method would be helpful for genetic improvement of traits with low heritability and that are difficult to measure. The identification of core genes can aid genomic selection for complex traits with low heritability that are difficult to measure. We performed statistical analysis of associations between longissimus dorsi muscle tenderness and gene expression in 20 Hanwoo cattle, using Warner-Bratzler shear force and RNAseq data, respectively. We found a total of 166 core genes, from which six (ASAP1, CAPN5, ELN, SUMF2, TTC8 and MGAT4A) were regulated by 16 cis-expression quantitative trait loci (eQTL) SNPs. Notably, we found that a cis-eQTL SNP of the ELN gene contained an MFZ-1 binding site in its putative promoter region. These findings provide useful information for genomic prediction of beef tenderness in Hanwoo cattle.

A human adipose tissue cell-type transcriptome atlas

The importance of defining cell-type-specific genes is well acknowledged. Technological advances facilitate high-resolution sequencing of single cells, but practical challenges remain. Adipose tissue is composed primarily of adipocytes, large buoyant cells requiring extensive, artefact-generating processing for separation and analysis. Thus, adipocyte data are frequently absent from single-cell RNA sequencing (scRNA-seq) datasets, despite being the primary functional cell type. Here, we decipher cell-type-enriched transcriptomes from unfractionated human adipose tissue RNA-seq data. We profile all major constituent cell types, using 527 visceral adipose tissue (VAT) or 646 subcutaneous adipose tissue (SAT) samples, identifying over 2,300 cell-type-enriched transcripts. Sex-subset analysis uncovers a panel of male-only cell-type-enriched genes. By resolving expression profiles of genes differentially expressed between SAT and VAT, we identify mesothelial cells as the primary driver of this variation. This study provides an accessible method to profile cell-type-enriched transcriptomes using bulk RNA-seq, generating a roadmap for adipose tissue biology.

Genotyping and Whole-Genome Resequencing of Welsh Sheep Breeds Reveal Candidate Genes and Variants for Adaptation to Local Environment and Socioeconomic Traits

Background

Advances in genetic tools applied to livestock breeding has prompted research into the previously neglected breeds adapted to harsh local environments. One such group is the Welsh mountain sheep breeds, which can be farmed at altitudes of 300 m above sea level but are considered to have a low productive value because of their poor wool quality and small carcass size. This is contrary to the lowland breeds which are more suited to wool and meat production qualities, but do not fare well on upland pasture. Herein, medium-density genotyping data from 317 individuals representing 15 Welsh sheep breeds were used alongside the whole-genome resequencing data of 14 breeds from the same set to scan for the signatures of selection and candidate genetic variants using haplotype- and SNP-based approaches.

Results

Haplotype-based selection scan performed on the genotyping data pointed to a strong selection in the regions of GBA3, PPARGC1A, APOB, and PPP1R16B genes in the upland breeds, and RNF24, PANK2, and MUC15 in the lowland breeds. SNP-based selection scan performed on the resequencing data pointed to the missense mutations under putative selection relating to a local adaptation in the upland breeds with functions such as angiogenesis (VASH1), anti-oxidation (RWDD1), cell stress (HSPA5), membrane transport (ABCA13 and SLC22A7), and insulin signaling (PTPN1 and GIGFY1). By contrast, genes containing candidate missense mutations in the lowland breeds are related to cell cycle (CDK5RAP2), cell adhesion (CDHR3), and coat color (MC1R).

Conclusion

We found new variants in genes with potentially functional consequences to the adaptation of local sheep to their environments in Wales. Knowledge of these variations is important for improving the adaptative qualities of UK and world sheep breeds through a marker-assisted selection.

Evolutionary and Functional Analysis of Korean Native Pig Using Single Nucleotide Polymorphisms

Time and cost-effective production of next-generation sequencing data has enabled the performance of population-scale comparative and evolutionary studies for various species, which are essential for obtaining the comprehensive insight into molecular mechanisms underlying species- or breed-specific traits. In this study, the evolutionary and functional analysis of Korean native pig (KNP) was performed using single nucleotide polymorphism (SNP) data by comparative and population genomic approaches with six different mammalian species and five pig breeds. We examined the evolutionary history of KNP SNPs, and the specific genes of KNP based on the uniqueness of non-synonymous SNPs among the used species and pig breeds. We discovered the evolutionary trajectory of KNP SNPs within the used mammalian species as well as pig breeds. We also found olfaction-associated functions that have been characterized and diversified during evolution, and quantitative trait loci associated with the unique traits of KNP. Our study provides new insight into the evolution of KNP and serves as a good example for a better understanding of domestic animals in terms of evolution and domestication using the combined approaches of comparative and population genomics.

Integrating Genome-Wide CNVs Into QTLs and High Confidence GWAScore Regions Identified Positional Candidates for Sheep Economic Traits

Copy number variations (CNVs) are important source of genetic variation, which can affect diverse economic traits through a variety of mechanisms. In addition, genome scan can identify many quantitative trait loci (QTLs) for the economic traits, while genome-wide association studies (GWAS) can localize genetic variants associated with the phenotypic variations. Here, we developed a method called GWAScore which collected GWAS summary data to identify potential candidates, and integrated CNVs into QTLs and high confidence GWAScore regions to detect crucial CNV markers for sheep growth traits. We got 197 candidate genes which were overlapping with the candidate CNVs. Some crucial genes (MYLK3, TTC29, HERC6, ABCG2, RUNX1, etc.) showed significantly elevated GWAScore peaks than other candidate genes. In this study, we developed the GWAScore method to excavate the potential value of candidate genes as markers for the sheep molecular breeding.

Genomic regions associated with muscularity in beef cattle differ in five contrasting cattle breeds

BACKGROUND

Linear type traits, which reflect the muscular characteristics of an animal, could provide insight into how, in some cases, morphologically very different animals can yield the same carcass weight. Such variability may contribute to differences in the overall value of the carcass since primal cuts vary greatly in price; such variability may also hinder successful genome-based association studies. Therefore, the objective of our study was to identify genomic regions that are associated with five muscularity linear type traits and to determine if these significant regions are common across five different breeds. Analyses were carried out using linear mixed models on imputed whole-genome sequence data in each of the five breeds, separately. Then, the results of the within-breed analyses were used to conduct an across-breed meta-analysis per trait.

RESULTS

We identified many quantitative trait loci (QTL) that are located across the whole genome and associated with each trait in each breed. The only commonality among the breeds and traits was a large-effect pleiotropic QTL on BTA2 that contained the MSTN gene, which was associated with all traits in the Charolais and Limousin breeds. Other plausible candidate genes were identified for muscularity traits including PDE1A, PPP1R1C and multiple collagen and HOXD genes. In addition, associated (gene ontology) GO terms and KEGG pathways tended to differ between breeds and between traits especially in the numerically smaller populations of Angus, Hereford, and Simmental breeds. Most of the SNPs that were associated with any of the traits were intergenic or intronic SNPs located within regulatory regions of the genome.

CONCLUSIONS

The commonality between the Charolais and Limousin breeds indicates that the genetic architecture of the muscularity traits may be similar in these breeds due to their similar origins. Conversely, there were vast differences in the QTL associated with muscularity in Angus, Hereford, and Simmental. Knowledge of these differences in genetic architecture between breeds is useful to develop accurate genomic prediction equations that can operate effectively across breeds. Overall, the associated QTL differed according to trait, which suggests that breeding for a morphologically different (e.g. longer and wider versus shorter and smaller) more efficient animal may become possible in the future.

Development and validation of a small SNP panel for feed efficiency in beef cattle

The objective of this study was to develop and validate a customized cost-effective single nucleotide polymorphism (SNP) panel for genetic improvement of feed efficiency in beef cattle. The SNPs identified in previous association studies and through extensive analysis of candidate genomic regions and genes, were screened for their functional impact and allele frequency in Angus and Hereford breeds used as validation candidates for the panel. Association analyses were performed on genotypes of 159 SNPs from new samples of Angus (n = 160), Hereford (n = 329), and Angus-Hereford crossbred (n = 382) cattle using allele substitution and genotypic models in ASReml. Genomic heritabilities were estimated for feed efficiency traits using the full set of SNPs, SNPs associated with at least one of the traits (at P ≤ 0.05 and P < 0.10), as well as the Illumina bovine 50K representing a widely used commercial genotyping panel. A total of 63 SNPs within 43 genes showed association (P ≤ 0.05) with at least one trait. The minor alleles of SNPs located in the GHR and CAST genes were associated with decreasing effects on residual feed intake (RFI) and/or RFI adjusted for backfat (RFIf), whereas minor alleles of SNPs within MKI67 gene were associated with increasing effects on RFI and RFIf. Additionally, the minor allele of rs137400016 SNP within CNTFR was associated with increasing average daily gain (ADG). The SNPs genotypes within UMPS, SMARCAL, CCSER1, and LMCD1 genes showed significant over-dominance effects whereas other SNPs located in SMARCAL1, ANXA2, CACNA1G, and PHYHIPL genes showed additive effects on RFI and RFIf. Gene enrichment analysis indicated that gland development, as well as ion and cation transport are important physiological mechanisms contributing to variation in feed efficiency traits. The study revealed the effect of the Jak-STAT signaling pathway on feed efficiency through the CNTFR, OSMR, and GHR genes. Genomic heritability using the 63 significant (P ≤ 0.05) SNPs was 0.09, 0.09, 0.13, 0.05, 0.05, and 0.07 for ADG, dry matter intake, midpoint metabolic weight, RFI, RFIf, and backfat, respectively. These SNPs contributed to genetic variation in the studied traits and thus can potentially be used or tested to generate cost-effective molecular breeding values for feed efficiency in beef cattle.

The effect of QTL-rich region polymorphisms identified by targeted DNA-seq on pig production traits

The aim of the present study was to analyse the effect of PLCD4, PECR, FN1 and PNKD mutations on pig productive traits and tested the usefulness of targeted enrichment DNA sequencing method as tool for preselection of genetic markers. The potential genetic markers for pig productive traits were identified by using targeted enrichment DNA sequencing of chromosome 15 region that is QTL-rich. The selected mutations were genotyped by using HRM, Sanger sequencing and PCR-ACRS methods. The association study was performed by using GLM model in SAS program and included over 500 pigs of 5 populations maintained in Poland. The variation (C/T) of PLCD4 gene affected feed conversion, intramuscular fat and water exudation. The PNKD mutations were associated with texture parameters measured after cooking. In turn, the variation rs792423408 (C/T) in the FN1 gene influenced toughness measured in semimembranosus muscle and growth traits that was observed particularly in Duroc breed. Summarizing, the investigated gene variants delivered valuable information that could be used during developing SNP microarray for genomic estimated breeding value procedure in pigs. Moreover, the study showed that the TEDNA-seq method could be used to preselect the molecular markers associated with pig traits. However, the further association study that included large number animal populations is necessary.

Fine mapping on chromosome 13q32-34 and brain expression analysis implicates MYO16 in schizophrenia

We previously reported linkage of schizophrenia and schizoaffective disorder to 13q32-34 in the European descent Afrikaner population from South Africa. The nature of genetic variation underlying linkage peaks in psychiatric disorders remains largely unknown and both rare and common variants may be contributing. Here, we examine the contribution of common variants located under the 13q32-34 linkage region. We used densely spaced SNPs to fine map the linkage peak region using both a discovery sample of 415 families and a meta-analysis incorporating two additional replication family samples. In a second phase of the study, we use one family-based data set with 237 families and independent case-control data sets for fine mapping of the common variant association signal using HapMap SNPs. We report a significant association with a genetic variant (rs9583277) within the gene encoding for the myosin heavy-chain Myr 8 (MYO16), which has been implicated in neuronal phosphoinositide 3-kinase signaling. Follow-up analysis of HapMap variation within MYO16 in a second set of Afrikaner families and additional case-control data sets of European descent highlighted a region across introns 2-6 as the most likely region to harbor common MYO16 risk variants. Expression analysis revealed a significant increase in the level of MYO16 expression in the brains of schizophrenia patients. Our results suggest that common variation within MYO16 may contribute to the genetic liability to schizophrenia.

Single nucleotide polymorphisms for feed efficiency and performance in crossbred beef cattle

Background

This study was conducted to: (1) identify new SNPs for residual feed intake (RFI) and performance traits within candidate genes identified in a genome wide association study (GWAS); (2) estimate the proportion of variation in RFI explained by the detected SNPs; (3) estimate the effects of detected SNPs on carcass traits to avoid undesirable correlated effects on these economically important traits when selecting for feed efficiency; and (4) map the genes to biological mechanisms and pathways. A total number of 339 SNPs corresponding to 180 genes were tested for association with phenotypes using a single locus regression (SLRM) and genotypic model on 726 and 990 crossbred animals for feed efficiency and carcass traits, respectively.

Results

Strong evidence of associations for RFI were located on chromosomes 8, 15, 16, 18, 19, 21, and 28. The strongest association with RFI (P = 0.0017) was found with a newly discovered SNP located on BTA 8 within the ELP3 gene. SNPs rs41820824 and rs41821600 on BTA 16 within the gene HMCN1 were strongly associated with RFI (P = 0.0064 and P = 0.0033, respectively). A SNP located on BTA 18 within the ZNF423 gene provided strong evidence for association with RFI (P = 0.0028). Genomic estimated breeding values (GEBV) from 98 significant SNPs were moderately correlated (0.47) to the estimated breeding values (EBVs) from a mixed animal model. The significant (P < 0.05) SNPs (98) explained 26% of the genetic variance for RFI. In silico functional analysis for the genes suggested 35 and 39 biological processes and pathways, respectively for feed efficiency traits.

Conclusions

This study identified several positional and functional candidate genes involved in important biological mechanisms associated with feed efficiency and performance. Significant SNPs should be validated in other populations to establish their potential utilization in genetic improvement programs.

Molecular heterogeneities of adipose depots - potential effects on adipose-muscle cross-talk in humans, mice and farm animals

Adipose tissue is considered as a major endocrine organ that secretes numerous proteins called adipokines. The heterogeneous nature of adipose tissue in different parts of the body suggests respective heterogeneity of proteomes and secretomes. This review consolidates knowledge from recent studies targeting the diversity of different adipose depots affecting the pattern of secreted adipokines and discusses potential consequences for the cross-talk between adipose and skeletal muscle in humans, rodent models and farm animals. Special attention is paid to muscle-associated fat depots like inter- and intramuscular fat that become focus of attention in the context of the rather new notion of skeletal muscle as a major endocrine organ. Understanding the complexity of communication between adipocytes and skeletal muscle cells will allow developing strategies for improvement of human health and for sustainable production of high quality meat.

Genomic imprinting and genetic effects on muscle traits in mice

BACKGROUND

Genomic imprinting refers to parent-of-origin dependent gene expression caused by differential DNA methylation of the paternally and maternally derived alleles. Imprinting is increasingly recognized as an important source of variation in complex traits, however, its role in explaining variation in muscle and physiological traits, especially those of commercial value, is largely unknown compared with genetic effects.

RESULTS

We investigated both genetic and genomic imprinting effects on key muscle traits in mice from the Berlin Muscle Mouse population, a key model system to study muscle traits. Using a genome scan, we first identified loci with either imprinting or genetic effects on phenotypic variation. Next, we established the proportion of phenotypic variation explained by additive, dominance and imprinted QTL and characterized the patterns of effects. In total, we identified nine QTL, two of which show large imprinting effects on glycogen content and potential, and body weight. Surprisingly, all imprinting patterns were of the bipolar type, in which the two heterozygotes are different from each other but the homozygotes are not. Most QTL had pleiotropic effects and explained up to 40% of phenotypic variance, with individual imprinted loci accounting for 4-5% of variation alone.

CONCLUSION

Surprisingly, variation in glycogen content and potential was only modulated by imprinting effects. Further, in contrast to general assumptions, our results show that genomic imprinting can impact physiological traits measured at adult stages and that the expression does not have to follow the patterns of paternal or maternal expression commonly ascribed to imprinting effects.

Quantitative trait loci segregating in crosses between New Hampshire and White Leghorn chicken lines: III. Fat deposition and intramuscular fat content

In this study, a genome scan was performed to detect genomic loci that affect fat deposition in white adipose tissues and muscles in 278 F (2) males of reciprocal crosses between the genetically and phenotypically extreme inbred chicken lines New Hampshire (NHI) and White Leghorn (WL77). Genome-wide highly significant quantitative trait loci (QTL) influencing fat deposition in white adipose tissues were found on GGA2 and 4. The peak QTL positions for different visceral and subcutaneous white adipose tissues were located between 41.4 and 112.4 Mb on GGA2 and between 76.2 and 78.7 Mb on GGA4, which explained 4.2-10.4% and 4.3-11.6% respectively of the phenotypic F (2) variances. Contrary to our expectations, the QTL allele descending from the lean line WL77 on GGA4 led to increased fat deposition. We suggest a transgressive action of the obesity allele only if it is not in the genetic background of the line WL77. Additional highly significant loci for subcutaneous adipose tissue mass were identified on GGA12 and 15. For intramuscular fat content, a suggestive QTL was located on GGA14. The analysed crosses provide a valuable resource for further fine mapping of fatness genes and subsequent gene discovery.