Metabolomic profiles indicate distinct physiological pathways affected by two loci with major divergent effect on Bos taurus growth and lipid deposition

Investigation of cosmopolitan and local Italian beef cattle breeds uncover common patterns of heterozygosity

The analysis of livestock heterozygosity is less common compared to the study of homozygous patterns. Heterozygous-Rich Regions (HRRs) may harbor significant loci for functional traits such as immune response, survival rate, and fertility. For this reason, this study was conducted to investigate and characterize the heterozygosity patterns of four beef cattle breeds, which included two cosmopolitan breeds (Limousine and Charolaise) and two local breeds (Sarda and Sardo Bruna). Our analysis identified regions with a high degree of heterozygosity using a consecutive runs approach, the Tajima D test, nucleotide diversity estimation, and Hardy Weinberg equilibrium test. These regions exhibited recurrent heterozygosity peaks and were consistently found on specific chromosomes across all breeds, specifically autosomes 15, 16, 20, and 23. The cosmopolitan and Sardo Bruna breeds also displayed peaks on autosomes 2 and 21, respectively. Thirty-five top runs shared by more than 25% of the populations were identified. These genomic fragments encompassed 18 genes, two of which are directly linked to male fertility, while four are associated with lactation. Two other genes play roles in survival and immune response. Our study also detected a region related to growth and carcass traits in Limousine breed. Our analysis of heterozygosity-rich regions revealed particular segments of the cattle genome linked to various functional traits. It appears that balancing selection is occurring in specific regions within the four examined breeds, and unexpectedly, they are common across cosmopolitan and local breeds. The genes identified hold potential for applications in breeding programs and conservation studies to investigate the phenotypes associated with these heterozygous genotypes. In addition, Tajima D test, Nucleotide diversity, and Hardy Weinberg equilibrium test confirmed the presence of heterozygous fragments found with Heterozygous-Rich Regions analysis.

Patterns of enrichment and acceleration in evolutionary rates of promoters suggest a role of regulatory regions in cetacean gigantism

BACKGROUND

Cetaceans (whales, porpoises, and dolphins) are a lineage of aquatic mammals from which some species became giants. Only recently, gigantism has been investigated from the molecular point of view. Studies focused mainly on coding regions, and no data on the influence of regulatory regions on gigantism in this group was available. Accordingly, we investigated the molecular evolution of non-coding regulatory regions of genes already described in the literature for association with size in mammals, focusing mainly on the promoter regions. For this, we used Ciiider and phyloP tools. Ciiider identifies significantly enriched transcription factor binding sites, and phyloP estimates the molecular evolution rate of the promoter.

RESULTS

We found evidence of enrichment of transcription binding factors related to large body size, with distinct patterns between giant and non-giant cetaceans in the IGFBP7 and NCAPG promoters, in which repressive agents are present in small cetaceans and those that stimulate transcription, in giant cetaceans. In addition, we found evidence of acceleration in the IGF2, IGFBP2, IGFBP7, and ZFAT promoters.

CONCLUSION

Our results indicate that regulatory regions may also influence cetaceans' body size, providing candidate genes for future research to understand the molecular basis of the largest living animals.

The molecular evolution of genes previously associated with large sizes reveals possible pathways to cetacean gigantism

Cetaceans are a group of aquatic mammals with the largest body sizes among living animals, including giant representatives such as blue and fin whales. To understand the genetic bases of gigantism in cetaceans, we performed molecular evolutionary analyses on five genes (GHSR, IGF2, IGFBP2, IGFBP7, and EGF) from the growth hormone/insulin-like growth factor axis, and four genes (ZFAT, EGF, LCORL, and PLAG1) previously described as related to the size of species evolutionarily close to cetaceans, such as pigs, cows, and sheep. Our dataset comprised 19 species of cetaceans, seven of which are classified as giants because they exceed 10 m in length. Our results revealed signs of positive selection in genes from the growth hormone/insulin-like growth factor axis and also in those related to body increase in cetacean-related species. In addition, pseudogenization of the EGF gene was detected in the lineage of toothless cetaceans, Mysticeti. Our results suggest the action of positive selection on gigantism in genes that act both in body augmentation and in mitigating its consequences, such as cancer suppression when involved in processes such as division, migration, and cell development control.

Exploring and Identifying Candidate Genes and Genomic Regions Related to Economically Important Traits in Hanwoo Cattle

The purpose of the current review was to explore and summarize different studies concerning the detection and characterization of candidate genes and genomic regions associated with economically important traits in Hanwoo beef cattle. Hanwoo cattle, the indigenous premium beef cattle of Korea, were introduced for their marbled fat, tenderness, characteristic flavor, and juiciness. To date, there has been a strong emphasis on the genetic improvement of meat quality and yields, such as backfat thickness (BFT), marbling score (MS), carcass weight (CW), eye muscle area (EMA), and yearling weight (YW), as major selection criteria in Hanwoo breeding programs. Hence, an understanding of the genetics controlling these traits along with precise knowledge of the biological mechanisms underlying the traits would increase the ability of the industry to improve cattle to better meet consumer demands. With the development of high-throughput genotyping, genomewide association studies (GWAS) have allowed the detection of chromosomal regions and candidate genes linked to phenotypes of interest. This is an effective and useful tool for accelerating the efficiency of animal breeding and selection. The GWAS results obtained from the literature review showed that most positional genes associated with carcass and growth traits in Hanwoo are located on chromosomes 6 and 14, among which LCORL, NCAPG, PPARGC1A, ABCG2, FAM110B, FABP4, DGAT1, PLAG1, and TOX are well known. In conclusion, this review study attempted to provide comprehensive information on the identified candidate genes associated with the studied traits and genes enriched in the functional terms and pathways that could serve as a valuable resource for future research in Hanwoo breeding programs.

Genome-wide association analysis of milk production, somatic cell score, and body conformation traits in Holstein cows

Milk production and body conformation traits are critical economic traits for dairy cows. To understand the basic genetic structure for those traits, a genome wide association study was performed on milk yield, milk fat yield, milk fat percentage, milk protein yield, milk protein percentage, somatic cell score, body form composite index, daily capacity composite index, feed, and leg conformation traits, based on the Illumina Bovine HD100k BeadChip. A total of 57, 12 and 26 SNPs were found to be related to the milk production, somatic cell score and body conformation traits in the Holstein cattle. Genes with pleiotropic effect were also found in this study. Seven significant SNPs were associated with multi-traits and were located on the PLEC, PLEKHA5, TONSL, PTGER4, and LCORL genes. In addition, some important candidate genes, like GPAT3, CEBPB, AGO2, SLC37A1, and FNDC3B, were found to participate in fat metabolism or mammary gland development. These results can be used as candidate genes for milk production, somatic cell score, and body conformation traits of Holstein cows, and are helpful for further gene function analysis to improve milk production and quality.

Genome-wide Association Study for Carcass Primal Cut Yields Using Single-step Bayesian Approach in Hanwoo Cattle

The importance of meat and carcass quality is growing in beef cattle production to meet both producer and consumer demands. Primal cut yields, which reflect the body compositions of carcass, could determine the carcass grade and, consequently, command premium prices. Despite its importance, there have been few genome-wide association studies on these traits. This study aimed to identify genomic regions and putative candidate genes related to 10 primal cut traits, including tenderloin, sirloin, striploin, chuck, brisket, top round, bottom round, shank, flank, and rib in Hanwoo cattle using a single-step Bayesian regression (ssBR) approach. After genomic data quality control, 43,987 SNPs from 3,745 genotyped animals were available, of which 3,467 had phenotypic records for the analyzed traits. A total of 16 significant genomic regions (1-Mb window) were identified, of which five large-effect quantitative trait loci (QTLs) located on chromosomes 6 at 38–39 Mb, 11 at 21–22 Mb, 14 at 6–7 Mb and 26–27 Mb, and 19 at 26–27 Mb were associated with more than one trait, while the remaining 11 QTLs were trait-specific. These significant regions were harbored by 154 genes, among which TOX, FAM184B, SPP1, IBSP, PKD2, SDCBP, PIGY, LCORL, NCAPG, and ABCG2 were noteworthy. Enrichment analysis revealed biological processes and functional terms involved in growth and lipid metabolism, such as growth (GO:0040007), muscle structure development (GO:0061061), skeletal system development (GO:0001501), animal organ development (GO:0048513), lipid metabolic process (GO:0006629), response to lipid (GO:0033993), metabolic pathways (bta01100), focal adhesion (bta04510), ECM–receptor interaction (bta04512), fat digestion and absorption (bta04975), and Rap1 signaling pathway (bta04015) being the most significant for the carcass primal cut traits. Thus, identification of quantitative trait loci regions and plausible candidate genes will aid in a better understanding of the genetic and biological mechanisms regulating carcass primal cut yields.

Genome‑wide association study and genomic prediction for growth traits in yellow-plumage chicken using genotyping-by-sequencing

Background

Growth traits are of great importance for poultry breeding and production and have been the topic of extensive investigation, with many quantitative trait loci (QTL) detected. However, due to their complex genetic background, few causative genes have been confirmed and the underlying molecular mechanisms remain unclear, thus limiting our understanding of QTL and their potential use for the genetic improvement of poultry. Therefore, deciphering the genetic architecture is a promising avenue for optimising genomic prediction strategies and exploiting genomic information for commercial breeding. The objectives of this study were to: (1) conduct a genome-wide association study to identify key genetic factors and explore the polygenicity of chicken growth traits; (2) investigate the efficiency of genomic prediction in broilers; and (3) evaluate genomic predictions that harness genomic features.

Results

We identified five significant QTL, including one on chromosome 4 with major effects and four on chromosomes 1, 2, 17, and 27 with minor effects, accounting for 14.5 to 34.1% and 0.2 to 2.6% of the genomic additive genetic variance, respectively, and 23.3 to 46.7% and 0.6 to 4.5% of the observed predictive accuracy of breeding values, respectively. Further analysis showed that the QTL with minor effects collectively had a considerable influence, reflecting the polygenicity of the genetic background. The accuracy of genomic best linear unbiased predictions (BLUP) was improved by 22.0 to 70.3% compared to that of the conventional pedigree-based BLUP model. The genomic feature BLUP model further improved the observed prediction accuracy by 13.8 to 15.2% compared to the genomic BLUP model.

Conclusions

A major QTL and four minor QTL were identified for growth traits; the remaining variance was due to QTL effects that were too small to be detected. The genomic BLUP and genomic feature BLUP models yielded considerably higher prediction accuracy compared to the pedigree-based BLUP model. This study revealed the polygenicity of growth traits in yellow-plumage chickens and demonstrated that the predictive ability can be greatly improved by using genomic information and related features.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00672-9.

Dietary Energy Level Impacts the Performance of Donkeys by Manipulating the Gut Microbiome and Metabolome

Considerable evidence suggests that dietary energy levels and gut microbiota are pivotal for animal health and productivity. However, little information exists about the correlations among dietary energy level, performance, and the gut microbiota and metabolome of donkeys. The objective of this study was to investigate the mechanisms by which dietary energy content dictates the growth performance by modulating the intestinal microbiome and metabolome of donkeys. Thirty-six nine-month-old male Dezhou donkeys with similar body weights were randomly assigned to two groups fed low- or high-energy diets (LE or HE). The results showed that donkeys fed HE had increased (p < 0.05) the average daily gain (ADG) and feed efficiency (G/F) compared with those that received LE diet. The gut microbiota in both groups was dominated by the phyla Firmicutes and Bacteroidetes regardless of the dietary energy level. However, feeding HE to donkeys significantly decreased (p < 0.05) the ratio of Firmicutes to Bacteroidetes (F/B). Compared to the LE group, feeding HE specifically increased the abundances of unidentified_Prevotellaceae (p = 0.02) while decreasing the richness of unidentified_Ruminococcaceae (p = 0.05). Compared to the LE group, feeding the HE diet significantly (p < 0.05) upregulated certain metabolic pathways involving the aspartate metabolism and the urea cycle. In addition, the increased bacteria and metabolites in the HE-fed group exhibited a positive correlation with improved growth performance of donkeys. Taken together, feeding the HE diet increased the richness of Prevotellaceae and upregulated growth-related metabolic pathways, which may have contributed to the ameliorated growth performance of donkeys. Thus, it is a recommendable dietary strategy to feed HE diets to fattening donkeys for superior product performance and feed efficiency.

Genomic analyses reveal distinct genetic architectures and selective pressures in Chinese donkeys

Donkey (Equus asinus) is an important livestock animal in China because of its draft and medicinal value. After a long period of natural and artificial selection, the variety and phenotype of donkeys have become abundant. We clarified the genetic and demographic characteristics of Chinese domestic donkeys and the selection pressures by analyzing 78 whole genomes from 12 breeds. According to population structure, most Chinese domestic donkeys showed a dominant ancestral type. However, the Chinese donkeys still represented a significant geographical distribution trend. In the selective sweep, gene annotation, functional enrichment, and differential expression analyses between large and small donkey groups, we identified selective signals, including NCAPG and LCORL, which are related to rapid growth and large body size. Our findings elucidate the evolutionary history and formation of different donkey breeds and provide theoretical insights into the genetic mechanism underlying breed characteristics and molecular breeding programs of donkey clades.

Identification of candidate genes associated with slaughter traits in F2 chicken population using genome-wide association study

Slaughter traits are crucial economic traits of chickens. We performed a GWAS to discover critical loci and candidate genes for 21 slaughter traits in an F2 chicken population resulting from crossing Luxi gamecocks and recessive white feather broilers. We found some SNPs and genes which were significantly associated with keel length, head length, body slope length, bilateral leg weight without shin, bilateral foot weight, subcutaneous fat thickness, heart weight, muscular stomach weight and glandular stomach weight. This study provides references for further investigation of slaughter traits and molecular breeding in chicken.

Metabolomics Comparison of Hanwoo (Bos taurus coreanae) Biofluids Using Proton Nuclear Magnetic Resonance Spectroscopy

The aim of this study was to identify the metabolomic profiles of rumen fluid, serum, and urine from Hanwoo (Bos taurus coreanae), using proton nuclear magnetic resonance (1H-NMR) spectroscopy. In all, 189, 110, and 188 metabolites were identified in rumen fluid, serum, and urine, and 107, 49, and 99 were quantified, respectively. Organic acids, carbohydrates, and aliphatic acyclic compound metabolites were present at the highest concentrations in rumen fluid, serum, and urine, respectively. In addition, acetate, glucose, and urea were the most highly concentrated individual metabolites in rumen fluid, serum, and urine, respectively. In all, 77 metabolites were commonly identified, and 19 were quantified across three biofluids. Metabolic pathway analysis showed that the common quantified metabolites could provide relevant information about three main metabolic pathways, phenylalanine, tyrosine, and tryptophan biosynthesis; caffeine metabolism; and histidine metabolism. These results can be useful as reference values for future metabolomic research on Hanwoo biofluids in Korea.

Genome-wide association studies reveal novel loci associated with carcass and body measures in beef cattle

Genomic selection has an essential role in the livestock economy by increasing selection productivity. Genomics provides a mechanism to increase the rate of genetic gain using marker-assisted selection. Various quantitative trait loci (QTL) associated with body, carcass and meat quality traits in beef cattle have been found. It is widely accepted that QTL traits in livestock species are regulated by several genes and factors from the environment. Genome-wide association studies (GWAS) are a powerful approach in identifying QTL and to establish genomic regions harboring the genes and polymorphisms associated with specific characteristics in beef cattle. Due to their impact on economic returns, growth, carcass and meat quality traits of cattle are frequently used as essential criteria in selection in breeding programs., GWAS has been used in beef cattle breeding and genetic program and some progress has been made. Furthermore, numerous genes and markers related to productivity traits in beef cattle have been found. This review summarizes the advances in the use of GWAS in beef cattle production and outlines the associations with growth, carcass, and meat quality.

The LCORL Locus is under Selection in Large-Sized Pakistani Goat Breeds

Goat domestication and human selection for valued traits have formed diverse breeds with characteristic phenotypes. This process led to the fixation of causative genetic variants controlling breed-specific traits within regions of reduced genetic diversity-so-called "selection signatures". We previously reported an analysis of selection signatures based on pooled whole-genome sequencing data of 20 goat breeds and bezoar goats. In the present study, we reanalyzed the data and focused on a subset of eight Pakistani goat breeds (Angora, Barbari, Beetal, Dera Din Panah, Kamori, Nachi, Pahari, Teddy). We identified 749 selection signatures based on reduced heterozygosity in these breeds. A search for signatures that are shared across large-sized goat breeds revealed that five medium-to-large-sized Pakistani goat breeds had a common selection signature on chromosome 6 in a region harboring the LCORL gene, which has been shown to modulate height or body size in several mammalian species. Fine-mapping of the region confirmed that all five goat breeds with the selection signature were nearly fixed for the same haplotype in a ~191 kb region spanning positions 37,747,447-37,938,449. From the pool sequencing data, we identified a frame-shifting single base insertion into an isoform-specific exon of LCORL as a potential candidate causal variant mediating the size-increasing effect. If this preliminary result can be confirmed in independent replication studies, genotyping of this variant might be used to improve breeding programs and the selection for stature in goats.

Biological Network Approach for the Identification of Regulatory Long Non-Coding RNAs Associated With Metabolic Efficiency in Cattle

Background: Genomic regions associated with divergent livestock feed efficiency have been found predominantly outside protein coding sequences. Long non-coding RNAs (lncRNA) can modulate chromatin accessibility, gene expression and act as important metabolic regulators in mammals. By integrating phenotypic, transcriptomic, and metabolomic data with quantitative trait locus data in prioritizing co-expression network analyses, we aimed to identify and functionally characterize lncRNAs with a potential key regulatory role in metabolic efficiency in cattle. Materials and Methods: Crossbred animals (n = 48) of a Charolais x Holstein F₂-population were allocated to groups of high or low metabolic efficiency based on residual feed intake in bulls, energy corrected milk in cows and intramuscular fat content in both genders. Tissue samples from jejunum, liver, skeletal muscle and rumen were subjected to global transcriptomic analysis via stranded total RNA sequencing (RNAseq) and blood plasma samples were used for profiling of 640 metabolites. To identify lncRNAs within the indicated tissues, a project-specific transcriptome annotation was established. Subsequently, novel transcripts were categorized for potential lncRNA status, yielding a total of 7,646 predicted lncRNA transcripts belonging to 3,287 loci. A regulatory impact factor approach highlighted 92, 55, 35, and 73 lncRNAs in jejunum, liver, muscle, and rumen, respectively. Their ensuing high regulatory impact factor scores indicated a potential regulatory key function in a gene set comprising loci displaying differential expression, tissue specificity and loci overlapping with quantitative trait locus regions for residual feed intake or milk production. These were subjected to a partial correlation and information theory analysis with the prioritized gene set. Results and Conclusions: Independent, significant and group-specific correlations (|r| > 0.8) were used to build a network for the high and the low metabolic efficiency group resulting in 1,522 and 1,732 nodes, respectively. Eight lncRNAs displayed a particularly high connectivity (>100 nodes). Metabolites and genes from the partial correlation and information theory networks, which each correlated significantly with the respective lncRNA, were included in an enrichment analysis indicating distinct affected pathways for the eight lncRNAs. LncRNAs associated with metabolic efficiency were classified to be functionally involved in hepatic amino acid metabolism and protein synthesis and in calcium signaling and neuronal nitric oxide synthase signaling in skeletal muscle cells.

Genetic mapping of distal femoral, stifle, and tibial radiographic morphology in dogs with cranial cruciate ligament disease

Cranial cruciate ligament disease (CCLD) is a complex trait. Ten measurements were made on orthogonal distal pelvic limb radiographs of 161 pure and mixed breed dogs with, and 55 without, cranial cruciate partial or complete ligament rupture. Dogs with CCLD had significantly smaller infrapatellar fat pad width, higher average tibial plateau angle, and were heavier than control dogs. The first PC weightings captured the overall size of the dog’s stifle and PC2 weightings reflected an increasing tibial plateau angle coupled with a smaller fat pad width. Of these dogs, 175 were genotyped, and 144,509 polymorphisms were used in a genome-wide association study with both a mixed linear and a multi-locus model. For both models, significant (p_genome <3.46×10⁻⁷ for the mixed and< 6.9x10^-8 for the multilocus model) associations were found for PC1, tibial diaphyseal length and width, fat pad base length, and femoral and tibial condyle width at LCORL, a known body size-regulating locus. Other body size loci with significant associations were growth hormone 1 (GH1), which was associated with the length of the fat pad base and the width of the tibial diaphysis, and a region on CFAX near IRS4 and ACSL4 in the multilocus model. The tibial plateau angle was associated significantly with a locus on CFA10 in the linear mixed model with nearest candidate genes BET1 and MYH9 and on CFA08 near candidate genes WDHD1 and GCH1. MYH9 has a major role in osteoclastogenesis. Our study indicated that tibial plateau slope is associated with CCLD and a compressed infrapatellar fat pad, a surrogate for stifle osteoarthritis. Because of the association between tibial plateau slope and CCLD, and pending independent validation, these candidate genes for tibial plateau slope may be tested in breeds susceptible to CCLD before they develop disease or are bred.

Ten years of the horse reference genome: insights into equine biology, domestication and population dynamics in the post-genome era

The horse reference genome from the Thoroughbred mare Twilight has been available for a decade and, together with advances in genomics technologies, has led to unparalleled developments in equine genomics. At the core of this progress is the continuing improvement of the quality, contiguity and completeness of the reference genome, and its functional annotation. Recent achievements include the release of the next version of the reference genome (EquCab3.0) and generation of a reference sequence for the Y chromosome. Horse satellite‐free centromeres provide unique models for mammalian centromere research. Despite extremely low genetic diversity of the Y chromosome, it has been possible to trace patrilines of breeds and pedigrees and show that Y variation was lost in the past approximately 2300 years owing to selective breeding. The high‐quality reference genome has led to the development of three different SNP arrays and WGSs of almost 2000 modern individual horses. The collection of WGS of hundreds of ancient horses is unique and not available for any other domestic species. These tools and resources have led to global population studies dissecting the natural history of the species and genetic makeup and ancestry of modern breeds. Most importantly, the available tools and resources, together with the discovery of functional elements, are dissecting molecular causes of a growing number of Mendelian and complex traits. The improved understanding of molecular underpinnings of various traits continues to benefit the health and performance of the horse whereas also serving as a model for complex disease across species.

Retinol binding protein 4 abundance in plasma and tissues is related to body fat deposition in cattle

Retinol binding protein 4 (RBP4) facilitates the transport of retinol in the body but is also an adipokine and fatty acid transporter. Our study was aimed at investigating the associations between RBP4 abundance and fat deposition in cattle. Blood samples of 246 crossbred bulls were taken at 8 months of age and at slaughter at 18 months of age for the determination of RBP4, hormone levels, and fatty acid composition. Significant correlations between plasma RBP4 abundance at 8 months of age and carcass traits at 18 months of age were detected (e.g., r = 0.3; P < 0.001 to carcass fat). Furthermore, RBP4 abundances in the plasma and subcutaneous fat were higher (P < 0.05) in bulls with increased fat deposition, whereas the liver RBP4 expression was not (P > 0.05). Retinol binding protein 4 was immunohistochemically localized in or close to adipocytes within muscle and adipose tissue and in liver stellate cells but not in hepatocytes. Overall, our results indicate that increased RBP4 levels were associated with increased fat deposition and altered fatty acid composition, but not with altered glucose tolerance, in crossbred bulls. Moreover, our results suggest that adipose-tissue-derived RBP4 may contribute to the circulating RBP4 level.

A genome-wide scan for diversifying selection signatures in selected horse breeds

The genetic differentiation of the current horse population was evolutionarily created by natural or artificial selection which shaped the genomes of individual breeds in several unique ways. The availability of high throughput genotyping methods created the opportunity to study this genetic variation on a genome-wide level allowing detection of genome regions divergently selected between separate breeds as well as among different horse types sharing similar phenotypic features. In this study, we used the population differentiation index (FST) that is generally used for measuring locus-specific allele frequencies variation between populations, to detect selection signatures among six horse breeds maintained in Poland. These breeds can be classified into three major categories, including light, draft and primitive horses, selected mainly in terms of type (utility), exterior, performance, size, coat color and appearance. The analysis of the most pronounced selection signals found in this study allowed us to detect several genomic regions and genes connected with processes potentially important for breed phenotypic differentiation and associated with energy homeostasis during physical effort, heart functioning, fertility, disease resistance and motor coordination. Our results also confirmed previously described association of loci on ECA3 (spanning LCORL and NCAPG genes) and ECA11 (spanning LASP1 gene) with the regulation of body size in our draft and primitive (small size) horses. The efficiency of the applied FST-based approach was also confirmed by the identification of a robust selection signal in the blue dun colored Polish Konik horses at the locus of TBX3 gene, which was previously shown to be responsible for dun coat color dilution in other horse breeds. FST-based method showed to be efficient in detection of diversifying selection signatures in the analyzed horse breeds. Especially pronounced signals were observed at the loci responsible for fixed breed-specific features. Several candidate genes under selection were proposed in this study for traits selected in separate breeds and horse types, however, further functional and comparative studies are needed to confirm and explain their effect on the observed genetic diversity of the horse breeds.

Reducing the interval of a growth QTL on chromosome 4 in laying hens

In our previous research, we identified a QTL with an interval of 3.4 Mb for growth on chicken chromosome (GGA) 4 in an advanced intercross population of an initial cross between the New Hampshire inbred line (NHI) and the White Leghorn inbred line (WL77). In the current study, an association analysis was performed in a population of purebred white layers (WLA) with White Leghorn origin. Genotypic data of 130 SNPs within the previously identified 3.4-Mb region were obtained using a 60K SNP chip. In total, 24 significant SNPs (LOD ≥ 4.44) on GGA4 were detected for daily weigh gain from 8 to 14 weeks and two SNPs (LOD ≥ 4.80) for body weight at 14 weeks. The QTL interval was reduced by 1.9 Mb to an interval of 1.5 Mb (74.6-76.1 Mb) that harbors 15 genes. Furthermore, to identify additional loci for chicken growth, a genome-wide association study (GWAS) was carried out in a WLA population. The GWAS identified an additional QTL on GGA6 for body weight at six weeks (19.8-21.2 Mb). Our findings showed that by using a WLA population we were able to further reduce the QTL confidence interval previously detected using a NHI × WL77 advanced intercross population.

Livestock metabolomics and the livestock metabolome: A systematic review

Metabolomics uses advanced analytical chemistry techniques to comprehensively measure large numbers of small molecule metabolites in cells, tissues and biofluids. The ability to rapidly detect and quantify hundreds or even thousands of metabolites within a single sample is helping scientists paint a far more complete picture of system-wide metabolism and biology. Metabolomics is also allowing researchers to focus on measuring the end-products of complex, hard-to-decipher genetic, epigenetic and environmental interactions. As a result, metabolomics has become an increasingly popular “omics” approach to assist with the robust phenotypic characterization of humans, crop plants and model organisms. Indeed, metabolomics is now routinely used in biomedical, nutritional and crop research. It is also being increasingly used in livestock research and livestock monitoring. The purpose of this systematic review is to quantitatively and objectively summarize the current status of livestock metabolomics and to identify emerging trends, preferred technologies and important gaps in the field. In conducting this review we also critically assessed the applications of livestock metabolomics in key areas such as animal health assessment, disease diagnosis, bioproduct characterization and biomarker discovery for highly desirable economic traits (i.e., feed efficiency, growth potential and milk production). A secondary goal of this critical review was to compile data on the known composition of the livestock metabolome (for 5 of the most common livestock species namely cattle, sheep, goats, horses and pigs). These data have been made available through an open access, comprehensive livestock metabolome database (LMDB, available at http://www.lmdb.ca). The LMDB should enable livestock researchers and producers to conduct more targeted metabolomic studies and to identify where further metabolome coverage is needed.

Searching for new loci and candidate genes for economically important traits through gene-based association analysis of Simmental cattle

Single-marker genome-wide association study (GWAS) is a convenient strategy of genetic analysis that has been successful in detecting the association of a number of single-nucleotide polymorphisms (SNPs) with quantitative traits. However, analysis of individual SNPs can only account for a small proportion of genetic variation and offers only limited knowledge of complex traits. This inadequacy may be overcome by employing a gene-based GWAS analytic approach, which can be considered complementary to the single-SNP association analysis. Here we performed an initial single-SNP GWAS for bone weight (BW) and meat pH value with a total of 770,000 SNPs in 1141 Simmental cattle. Additionally, 21836 cattle genes collected from the Ensembl Genes 83 database were analyzed to find supplementary evidence to support the importance of gene-based association study. Results of the single SNP-based association study showed that there were 11 SNPs significantly associated with bone weight (BW) and two SNPs associated with meat pH value. Interestingly, all of these SNPs were located in genes detected by the gene-based association study.

Adrenal cortex expression quantitative trait loci in a German Holstein × Charolais cross

Background

The importance of the adrenal gland in regard to lactation and reproduction in cattle has been recognized early. Caused by interest in animal welfare and the impact of stress on economically important traits in farm animals the adrenal gland and its function within the stress response is of increasing interest. However, the molecular mechanisms and pathways involved in stress-related effects on economically important traits in farm animals are not fully understood. Gene expression is an important mechanism underlying complex traits, and genetic variants affecting the transcript abundance are thought to influence the manifestation of an expressed phenotype. We therefore investigated the genetic background of adrenocortical gene expression by applying an adaptive linear rank test to identify genome-wide expression quantitative trait loci (eQTL) for adrenal cortex transcripts in cattle.

Results

A total of 10,986 adrenal cortex transcripts and 37,204 single nucleotide polymorphisms (SNPs) were analysed in 145 F2 cows of a Charolais × German Holstein cross. We identified 505 SNPs that were associated with the abundance of 129 transcripts, comprising 482 cis effects and 17 trans effects. These SNPs were located on all chromosomes but X, 16, 24 and 28. Associated genes are mainly involved in molecular and cellular functions comprising free radical scavenging, cellular compromise, cell morphology and lipid metabolism, including genes such as CYP27A1 and LHCGR that have been shown to affect economically important traits in cattle.

Conclusions

In this study we showed that adrenocortical eQTL affect the expression of genes known to contribute to the phenotypic manifestation in cattle. Furthermore, some of the identified genes and related molecular pathways were previously shown to contribute to the phenotypic variation of behaviour, temperament and growth at the onset of puberty in the same population investigated here. We conclude that eQTL analysis appears to be a useful approach providing insight into the molecular and genetic background of complex traits in cattle and will help to understand molecular networks involved.

Electronic supplementary material

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

Population structure of eleven Spanish ovine breeds and detection of selective sweeps with BayeScan and hapFLK

The goals of the current work were to analyse the population structure of 11 Spanish ovine breeds and to detect genomic regions that may have been targeted by selection. A total of 141 individuals were genotyped with the Infinium 50 K Ovine SNP BeadChip (Illumina). We combined this dataset with Spanish ovine data previously reported by the International Sheep Genomics Consortium (N = 229). Multidimensional scaling and Admixture analyses revealed that Canaria de Pelo and, to a lesser extent, Roja Mallorquina, Latxa and Churra are clearly differentiated populations, while the remaining seven breeds (Ojalada, Castellana, Gallega, Xisqueta, Ripollesa, Rasa Aragonesa and Segureña) share a similar genetic background. Performance of a genome scan with BayeScan and hapFLK allowed us identifying three genomic regions that are consistently detected with both methods i.e. Oar3 (150–154 Mb), Oar6 (4–49 Mb) and Oar13 (68–74 Mb). Neighbor-joining trees based on polymorphisms mapping to these three selective sweeps did not show a clustering of breeds according to their predominant productive specialization (except the local tree based on Oar13 SNPs). Such cryptic signatures of selection have been also found in the bovine genome, posing a considerable challenge to understand the biological consequences of artificial selection.

Genome-wide association study for carcass traits, fatty acid composition, chemical composition, sugar, and the effects of related candidate genes in Japanese Black cattle

We performed a genome-wide association study (GWAS) and candidate gene analysis to: (i) evaluate the effectiveness of the GWAS in our small population by performing GWAS for carcass weight (CW) and fatty acid composition; (ii) detect novel candidate regions affecting non-CW carcass traits, chemical composition and sugar; and (iii) evaluate the association of the candidate genes previously detected in CW and fatty acid composition with other economically important traits. A total of 574 Japanese Black cattle and 40 657 Single nucleotide polymorphisms were used. In addition, candidate gene analyses were performed to evaluate the association of three CW-related genes and two fatty acid-related genes with carcass traits, fatty acid composition, chemical composition and sugar. The significant regions with the candidate genes were detected for CW and fatty acid composition, and these results showed that a significant region would be detectable despite the small sample size. The novel candidate regions were detected on BTA23 for crude protein and on BTA19 for fructose. CW-related genes associated with the rib-eye area and fatty acid composition were identified, and fatty acid-related genes had no relationship with other traits. Moreover, the favorable allele of CW-related genes had an unfavorable effect on fatty acid composition.

Uncovering Adaptation from Sequence Data: Lessons from Genome Resequencing of Four Cattle Breeds

Detecting the molecular basis of adaptation is one of the major questions in population genetics. With the advance in sequencing technologies, nearly complete interrogation of genome-wide polymorphisms in multiple populations is becoming feasible in some species, with the expectation that it will extend quickly to new ones. Here, we investigate the advantages of sequencing for the detection of adaptive loci in multiple populations, exploiting a recently published data set in cattle (Bos taurus). We used two different approaches to detect statistically significant signals of positive selection: a within-population approach aimed at identifying hard selective sweeps and a population-differentiation approach that can capture other selection events such as soft or incomplete sweeps. We show that the two methods are complementary in that they indeed capture different kinds of selection signatures. Our study confirmed some of the well-known adaptive loci in cattle (e.g., MC1R, KIT, GHR, PLAG1, NCAPG/LCORL) and detected some new ones (e.g., ARL15, PRLR, CYP19A1, PPM1L). Compared to genome scans based on medium- or high-density SNP data, we found that sequencing offered an increased detection power and a higher resolution in the localization of selection signatures. In several cases, we could even pinpoint the underlying causal adaptive mutation or at least a very small number of possible candidates (e.g., MC1R, PLAG1). Our results on these candidates suggest that a vast majority of adaptive mutations are likely to be regulatory rather than protein-coding variants.

NCAPG is differentially expressed during longissimus muscle development and is associated with growth traits in Chinese Qinchuan beef cattle

Based on RNA-seq analysis, we recently found that the bovine NCAPG (non-SMC condensin I complex, subunit G) gene is differentially expressed during development of the longissimus muscle. In the present study, we validated this result and, using quantitative real-time PCR analysis, identified two adjacent genes, LCORL and DCAF16, that are more abundant in fetal muscle tissue; further analysis of tissue-specific expression patterns indicated high abundance of NCAPG in muscle. Since no polymorphisms were detected in a previous study of Qinchuan cattle, we extended our investigation to examine the occurrence of single-nucleotide polymorphisms (SNPs) in the NCAPG gene. Three SNPs, i.e., one located in the intron region (g47747: T > G), a synonymous mutation (g52535: A > G) and a missense mutation (g53208: T > G) that leads to a change in the amino acid of interest (pIle442Met), were detected in a population of Qinchuan beef cattle (n = 300). Association analysis showed that these SNPs were significantly associated with the growth traits of Qinchuan beef cattle. Our results indicate that the bovine NCAPG gene may be involved in the development of the longissimus muscle. These polymorphisms in the NCAPG gene may be useful for marker-assisted selection of optimal body size in Qinchuan beef cattle.

A Non-Synonymous HMGA2 Variant Decreases Height in Shetland Ponies and Other Small Horses

The identification of quantitative trait loci (QTL) such as height and their underlying causative variants is still challenging and often requires large sample sizes. In humans hundreds of loci with small effects control the heritable portion of height variability. In domestic animals, typically only a few loci with comparatively large effects explain a major fraction of the heritability. We investigated height at withers in Shetland ponies and mapped a QTL to ECA 6 by genome-wide association (GWAS) using a small cohort of only 48 animals and the Illumina equine SNP70 BeadChip. Fine-mapping revealed a shared haplotype block of 793 kb in small Shetland ponies. The HMGA2 gene, known to be associated with height in horses and many other species, was located in the associated haplotype. After closing a gap in the equine reference genome we identified a non-synonymous variant in the first exon of HMGA2 in small Shetland ponies. The variant was predicted to affect the functionally important first AT-hook DNA binding domain of the HMGA2 protein (c.83G>A; p.G28E). We assessed the functional impact and found impaired DNA binding of a peptide with the mutant sequence in an electrophoretic mobility shift assay. This suggests that the HMGA2 variant also affects DNA binding in vivo and thus leads to reduced growth and a smaller stature in Shetland ponies. The identified HMGA2 variant also segregates in several other pony breeds but was not found in regular-sized horse breeds. We therefore conclude that we identified a quantitative trait nucleotide for height in horses.

Genome-wide association study of body weight in Australian Merino sheep reveals an orthologous region on OAR6 to human and bovine genomic regions affecting height and weight

Background

Body weight (BW) is an important trait for meat production in sheep. Although over the past few years, numerous quantitative trait loci (QTL) have been detected for production traits in cattle, few QTL studies have been reported for sheep, with even fewer on meat production traits. Our objective was to perform a genome-wide association study (GWAS) with the medium-density Illumina Ovine SNP50 BeadChip to identify genomic regions and corresponding haplotypes associated with BW in Australian Merino sheep.

Methods

A total of 1781 Australian Merino sheep were genotyped using the medium-density Illumina Ovine SNP50 BeadChip. Among the 53 862 single nucleotide polymorphisms (SNPs) on this array, 48 640 were used to perform a GWAS using a linear mixed model approach. Genotypes were phased with hsphase; to estimate SNP haplotype effects, linkage disequilibrium blocks were identified in the detected QTL region.

Results

Thirty-nine SNPs were associated with BW at a Bonferroni-corrected genome-wide significance threshold of 1 %. One region on sheep (Ovis aries) chromosome 6 (OAR6) between 36.15 and 38.56 Mb, included 13 significant SNPs that were associated with BW; the most significant SNP was OAR6_41936490.1 (P = 2.37 × 10⁻¹⁶) at 37.69 Mb with an allele substitution effect of 2.12 kg, which corresponds to 0.248 phenotypic standard deviations for BW. The region that surrounds this association signal on OAR6 contains three genes: leucine aminopeptidase 3 (LAP3), which is involved in the processing of the oxytocin precursor; NCAPG non-SMC condensin I complex, subunit G (NCAPG), which is associated with foetal growth and carcass size in cattle; and ligand dependent nuclear receptor corepressor-like (LCORL), which is associated with height in humans and cattle.

Conclusions

The GWAS analysis detected 39 SNPs associated with BW in sheep and a major QTL region was identified on OAR6. In several other mammalian species, regions that are syntenic with this region have been found to be associated with body size traits, which may reflect that the underlying biological mechanisms share a common ancestry. These findings should facilitate the discovery of causative variants for BW and contribute to marker-assisted selection.

Electronic supplementary material

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

PLAG1 and NCAPG-LCORL in livestock

A recent progress on stature genetics has revealed simple genetic architecture in livestock animals in contrast to that in humans. PLAG1 and/or NCAPG‐LCORL, both of which are known as a locus for adult human height, have been detected for association with body weight/height in cattle and horses, and for selective sweep in dogs and pigs. The findings indicate a significant impact of these loci on mammalian growth or body size and usefulness of the natural variants for selective breeding. However, association with an unfavorable trait, such as late puberty or risk for a neuropathic disease, was also reported for the respective loci, indicating an importance to discriminate between causality and association. Here I review the recent findings on quantitative trait loci (QTL) for stature in livestock animals, mainly focusing on the PLAG1 and NCAPG‐LCORL loci. I also describe our recent efforts to identify the causative variation for the third major locus for carcass weight in Japanese Black cattle.

Metabolic biomarkers for chronic kidney disease

Chronic kidney disease (CKD) is an increasingly recognized burden for patients and health care systems with high (and growing) global incidence and prevalence, significant mortality, and disproportionately high treatment costs. Yet, the available diagnostic tools are either impractical in clinical routine or have serious shortcomings impeding a well-informed disease management although optimized treatment strategies with proven benefits for the patients have become available. Advances in bioanalytical technologies have facilitated studies that identified genomic, proteomic, and metabolic biomarker candidates, and confirmed some of them in independent cohorts. This review summarizes the CKD-related markers discovered so far, and focuses on compounds and pathways, for which there is quantitative data, substantiating evidence from translational research, and a mechanistic understanding of the processes involved. Also, multiparametric marker panels have been suggested that showed promising diagnostic and prognostic performance in initial analyses although the data basis from prospective trials is very limited. Large-scale studies, however, are underway and will provide the information for validating a set of parameters and discarding others. Finally, the path from clinical research to a routine application is discussed, focusing on potential obstacles such as the use of mass spectrometry, and the feasibility of obtaining regulatory approval for targeted metabolomics assays.

Loci associated with adult stature also affect calf birth survival in cattle

Background

Understanding the underlying pleiotropic relationships among quantitative traits is necessary in order to predict correlated responses to artificial selection. The availability of large-scale next-generation sequence data in cattle has provided an opportunity to examine whether pleiotropy is responsible for overlapping QTL in multiple economic traits. In the present study, we examined QTL affecting cattle stillbirth, calf size, and adult stature located in the same genomic region.

Results

A genome scan using imputed whole genome sequence variants revealed one QTL with large effects on the service sire calving index (SCI), and body conformation index (BCI) at the same location (~39 Mb) on chromosome 6 in Nordic Red cattle. The targeted region was analyzed for SCI and BCI component traits. The QTL peak included LCORL and NCAPG genes, which had been reported to influence fetal growth and adult stature in several species. The QTL exhibited large effects on calf size and stature in Nordic Red cattle. Two deviant haplotypes (HAP1 and HAP2) were resolved which increased calf size at birth, and affected adult body conformation. However, the haplotypes also resulted in increased calving difficulties and calf mortality due to increased calf size at birth. Haplotype locations overlapped, however linkage disequilibrium (LD) between the sites was low, suggesting that two independent mutations were responsible for similar effects. The difference in prevalence between the two haplotypes in Nordic Red subpopulations suggested independent origins in different populations.

Conclusions

Results of our study identified QTL with large effects on body conformation and service sire calving traits on chromosome 6 in cattle. We present robust evidence that variation at the LCORL and NCAPG locus affects calf size at birth and adult stature. We suggest the two deviant haplotypes within the QTL were due to two independent mutations.

Electronic supplementary material

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

Systems biology analysis merging phenotype, metabolomic and genomic data identifies Non-SMC Condensin I Complex, Subunit G (NCAPG) and cellular maintenance processes as major contributors to genetic variability in bovine feed efficiency

Feed efficiency is a paramount factor for livestock economy. Previous studies had indicated a substantial heritability of several feed efficiency traits. In our study, we investigated the genetic background of residual feed intake, a commonly used parameter of feed efficiency, in a cattle resource population generated from crossing dairy and beef cattle. Starting from a whole genome association analysis, we subsequently performed combined phenotype-metabolome-genome analysis taking a systems biology approach by inferring gene networks based on partial correlation and information theory approaches. Our data about biological processes enriched with genes from the feed efficiency network suggest that genetic variation in feed efficiency is driven by genetic modulation of basic processes relevant to general cellular functions. When looking at the predicted upstream regulators from the feed efficiency network, the Tumor Protein P53 (TP53) and Transforming Growth Factor beta 1 (TGFB1) genes stood out regarding significance of overlap and number of target molecules in the data set. These results further support the hypothesis that TP53 is a major upstream regulator for genetic variation of feed efficiency. Furthermore, our data revealed a significant effect of both, the Non-SMC Condensin I Complex, Subunit G (NCAPG) I442M (rs109570900) and the Growth /differentiation factor 8 (GDF8) Q204X (rs110344317) loci, on residual feed intake and feed conversion. For both loci, the growth promoting allele at the onset of puberty was associated with a negative, but favorable effect on residual feed intake. The elevated energy demand for increased growth triggered by the NCAPG 442M allele is obviously not fully compensated for by an increased efficiency in converting feed into body tissue. As a consequence, the individuals carrying the NCAPG 442M allele had an additional demand for energy uptake that is reflected by the association of the allele with increased daily energy intake as observed in our study.

QTLs associated with dry matter intake, metabolic mid-test weight, growth and feed efficiency have little overlap across 4 beef cattle studies

Background

The identification of genetic markers associated with complex traits that are expensive to record such as feed intake or feed efficiency would allow these traits to be included in selection programs. To identify large-effect QTL, we performed a series of genome-wide association studies and functional analyses using 50 K and 770 K SNP genotypes scored in 5,133 animals from 4 independent beef cattle populations (Cycle VII, Angus, Hereford and Simmental × Angus) with phenotypes for average daily gain, dry matter intake, metabolic mid-test body weight and residual feed intake.

Results

A total of 5, 6, 11 and 10 significant QTL (defined as 1-Mb genome windows with Bonferroni-corrected P-value <0.05) were identified for average daily gain, dry matter intake, metabolic mid-test body weight and residual feed intake, respectively. The identified QTL were population-specific and had little overlap across the 4 populations. The pleiotropic or closely linked QTL on BTA 7 at 23 Mb identified in the Angus population harbours a promising candidate gene ACSL6 (acyl-CoA synthetase long-chain family member 6), and was the largest effect QTL associated with dry matter intake and mid-test body weight explaining 10.39% and 14.25% of the additive genetic variance, respectively. Pleiotropic or closely linked QTL associated with average daily gain and mid-test body weight were detected on BTA 6 at 38 Mb and BTA 7 at 93 Mb confirming previous reports. No QTL for residual feed intake explained more than 2.5% of the additive genetic variance in any population. Marker-based estimates of heritability ranged from 0.21 to 0.49 for residual feed intake across the 4 populations.

Conclusions

This GWAS study, which is the largest performed for feed efficiency and its component traits in beef cattle to date, identified several large-effect QTL that cumulatively explained a significant percentage of additive genetic variance within each population. Differences in the QTL identified among the different populations may be due to differences in power to detect QTL, environmental variation, or differences in the genetic architecture of trait variation among breeds. These results enhance our understanding of the biology of growth, feed intake and utilisation in beef cattle.

Large-effect pleiotropic or closely linked QTL segregate within and across ten US cattle breeds

Background

The availability of high-density SNP assays including the BovineSNP50 (50 K) enables the identification of novel quantitative trait loci (QTL) and improvement of the resolution of the locations of previously mapped QTL. We performed a series of genome-wide association studies (GWAS) using 50 K genotypes scored in 18,274 animals from 10 US beef cattle breeds with observations for twelve body weights, calving ease and carcass traits.

Results

A total of 159 large-effects QTL (defined as 1-Mb genome windows explaining more than 1% of additive genetic variance) were identified. In general, more QTL were identified in analyses with bigger sample sizes. Four large-effect pleiotropic or closely linked QTLs located on BTA6 at 37–42 Mb (primarily at 38 Mb), on BTA7 at 93 Mb, on BTA14 at 23–26 Mb (primarily at 25 Mb) and on BTA20 at 4 Mb were identified in more than one breed. Several breed-specific large-effect pleiotropic or closely linked QTL were also identified. Some identified QTL regions harbor genes known to have large effects on a variety of traits in cattle such as PLAG1 and MSTN and others harbor promising candidate genes including NCAPG, ARRDC3, ERGIC1, SH3PXD2B, HMGA2, MSRB3, LEMD3, TIGAR, SEPT7, and KIRREL3. Gene ontology analysis revealed that genes involved in ossification and in adipose tissue development were over-represented in the identified pleiotropic QTL. Also, the MAPK signaling pathway was identified as a common pathway affected by the genes located near the pleiotropic QTL.

Conclusions

This largest GWAS ever performed in beef cattle, led us to discover several novel across-breed and breed-specific large-effect pleiotropic QTL that cumulatively account for a significant percentage of additive genetic variance (e.g. more than a third of additive genetic variance of birth and mature weights; and calving ease direct in Hereford). These results will improve our understanding of the biology of growth and body composition in cattle.

Electronic supplementary material

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

Locus characterization and gene expression of bovine FNDC5: is the myokine irisin relevant in cattle?

The transmembrane protein FNDC5 was recently characterized as precursor of an exercise induced myokine named irisin. Previous studies found a relationship between circulating irisin levels and muscle mass in humans. Consequently, we tested the hypothesis whether FNDC5/irisin is involved in the modulation of body composition in cattle. Since information on the bovine FNDC5 locus was scarce, we characterized the gene experimentally as prerequisite for these investigations. We provide here a revised and extended gene model for bovine FNDC5. Although similarly organized like the human and murine loci, a higher variability was observed at transcript level in the bovine locus. FNDC5 mRNA was abundant in bovine skeletal muscle and was detected at lower levels in adipose tissue and liver. There were no expression differences between two groups of bulls highly different in muscularity and adiposity. Full-length FNDC5 protein (25 kDa) was present in bovine skeletal muscle independent of muscularity. Neither FNDC5 nor its putatively secreted peptide irisin were found in circulation of bulls. In contrast, we demonstrated that FNDC5 (25 kDa) and irisin (12 kDa) were present in murine skeletal muscle and that irisin was circulating in murine serum. This indicates fundamental differences in the regulation of FNDC5 and irisin between rodents and cattle.

Adipose and muscle tissue gene expression of two genes (NCAPG and LCORL) located in a chromosomal region associated with cattle feed intake and gain

A region on bovine chromosome 6 has been implicated in cattle birth weight, growth, and length. Non-SMC conodensin I complex subunit G (NCAPG) and ligand dependent nuclear receptor corepressor-like protein (LCORL) are positional candidate genes within this region. Previously identified genetic markers in both genes were associated with average daily gain (ADG) and average daily feed intake (ADFI) in a crossbred population of beef steers. These markers were also associated with hot carcass weight, ribeye area and adjusted fat thickness suggesting that they may have a role in lean muscle growth and/or fat deposition. The purpose of this study was to determine whether the transcript abundance of either of these genes in cattle adipose and muscle tissue was associated with variation in feed intake and average daily gain phenotypes. Transcript abundance for NCAPG and LCORL in adipose and muscle tissue was measured in heifers (adipose only), cows and steers using real-time polymerase chain reaction. In the adipose tissue from cows and heifers, a negative correlation between LCORL transcript abundance and ADFI were detected (P = 0.05). In the muscle tissue from cows, transcript abundance of NCAPG was associated with ADG (r = 0.26; P = 0.009). A positive correlation between LCORL transcript abundance from muscle tissue of steers and ADFI was detected (P = 0.04). LCORL protein levels in the muscle of steers were investigated and were associated with ADFI (P = 0.01). These data support our earlier genetic associations with ADFI and ADG within this region and represent the potential for biological activity of these genes in the muscle and adipose tissues of beef cattle; however, they also suggest that sex, age and/or nutrition-specific interactions may affect the expression of NCAPG and LCORL in these tissues.

A systems biology approach using metabolomic data reveals genes and pathways interacting to modulate divergent growth in cattle

BACKGROUND

Systems biology enables the identification of gene networks that modulate complex traits. Comprehensive metabolomic analyses provide innovative phenotypes that are intermediate between the initiator of genetic variability, the genome, and raw phenotypes that are influenced by a large number of environmental effects. The present study combines two concepts, systems biology and metabolic analyses, in an approach without prior functional hypothesis in order to dissect genes and molecular pathways that modulate differential growth at the onset of puberty in male cattle. Furthermore, this integrative strategy was applied to specifically explore distinctive gene interactions of non-SMC condensin I complex, subunit G (NCAPG) and myostatin (GDF8), known modulators of pre- and postnatal growth that are only partially understood for their molecular pathways affecting differential body weight.

RESULTS

Our study successfully established gene networks and interacting partners affecting growth at the onset of puberty in cattle. We demonstrated the biological relevance of the created networks by comparison to randomly created networks. Our data showed that GnRH (Gonadotropin-releasing hormone) signaling is associated with divergent growth at the onset of puberty and revealed two highly connected hubs, BTC and DGKH, within the network. Both genes are known to directly interact with the GnRH signaling pathway. Furthermore, a gene interaction network for NCAPG containing 14 densely connected genes revealed novel information concerning the functional role of NCAPG in divergent growth.

CONCLUSIONS

Merging both concepts, systems biology and metabolomic analyses, successfully yielded new insights into gene networks and interacting partners affecting growth at the onset of puberty in cattle. Genetic modulation in GnRH signaling was identified as key modifier of differential cattle growth at the onset of puberty. In addition, the benefit of our innovative concept without prior functional hypothesis was demonstrated by data suggesting that NCAPG might contribute to vascular smooth muscle contraction by indirect effects on the NO pathway via modulation of arginine metabolism. Our study shows for the first time in cattle that integration of genetic, physiological and metabolomics data in a systems biology approach will enable (or contribute to) an improved understanding of metabolic and gene networks and genotype-phenotype relationships.

Association of height and pubertal timing with lipoprotein subclass profile: exploring the role of genetic and environmental effects

OBJECTIVES

Little is known about the relationship between growth and lipoprotein profile. We aimed to analyze common genetic and environmental factors in the association of height from late childhood to adulthood and pubertal timing with serum lipid and lipoprotein subclass profile.

METHODS

A longitudinal cohort of Finnish twin pairs (FinnTwin12) was analyzed using self-reported height at 11-12, 14, 17 years and measured stature at adult age (21-24 years). Data were available for 719 individual twins including 298 complete pairs. Serum lipids and lipoprotein subclasses were measured by proton nuclear magnetic resonance spectroscopy. Multivariate variance component models for twin data were fitted. Cholesky decomposition was used to partition the phenotypic covariation among traits into additive genetic and unique environmental correlations.

RESULTS

In men, the strongest associations for both adult height and puberty were observed with total cholesterol, low-density lipoprotein cholesterol, intermediate-density lipoprotein cholesterol, and low-density lipoprotein particle subclasses (max. r = -0.19). In women, the magnitude of the correlations was weaker (max. r = -0.13). Few associations were detected between height during adolescence and adult lipid profile. Early onset of puberty was related to an adverse lipid profile, but delayed pubertal development in girls was associated with an unfavorable profile, as well. All associations were mediated mainly by additive genetic factors, but unique environmental effects cannot be disregarded.

CONCLUSIONS

Early puberty and shorter adult height relate to higher concentrations of atherogenic lipids and lipoprotein particles in early adulthood. Common genetic effects behind these phenotypes substantially contribute to the observed associations.

A genome-wide association study indicates LCORL/NCAPG as a candidate locus for withers height in German Warmblood horses

A genome-wide association scan for loci affecting withers height was conducted in 782 German Warmblood stallions, which were genotyped using the Illumina EquineSNP50 Bead Chip. A principal components approach was applied to correct for population structure. The analysis revealed a single major QTL on ECA3 explaining ~18 per cent of the phenotypic variance, which is in concordance with recent reports from other horse populations. The LCORL/NCAPG locus represents a strong candidate gene for this QTL. This locus is among a small number that have consistently been identified to influence human height in several large meta-analyses. Furthermore, a mutation within the NCAPG gene was found to affect growth and body frame size in cattle. Together with the results of this study in German Warmbloods, these findings strongly indicate LCORL/NCAPG as a candidate locus for withers height in horses. Further studies are, however, needed to confirm this.

Comparison of the effects explained by variations in the bovine PLAG1 and NCAPG genes on daily body weight gain, linear skeletal measurements and carcass traits in Japanese Black steers from a progeny testing program

The c.1326T>G single nucleotide polymorphism (SNP) in the NCAPG gene, which leads to an amino acid change of Ile442 to Met442, was previously identified as a candidate causative variation for a bovine carcass weight quantitative trait loci (QTL) on chromosome 6, which was associated with linear skeletal measurement gains and daily body weight gain at puberty. Recently, we identified the stature quantitative trait nucleotides (QTNs) in the PLAG1-CHCHD7 intergenic region as the causative variations for another carcass weight QTL on chromosome 14. This study aimed to compare the effects of the two QTL on growth and carcass traits using 768 Japanese Black steers from a progeny testing program and to determine whether a genetic interaction was present between them. The FJX_250879 SNP representing the stature QTL was associated with linear skeletal measurements and average daily body weight gain at early and late periods during adolescence. A genetic interaction between FJX_250879 and NCAPG c.1326T>G was detected only for body and rump lengths. Both were associated with increased carcass weight and Longissimus muscle area, and NCAPG c.1326T>G was also associated with reduced subcutaneous fat thickness and increased carcass yield estimate. These results will provide useful information to improve carcass weight in Japanese Black cattle.

Investigation of the genetic architecture of a bone carcass weight QTL on BTA6

A previous analysis of an F(2) /Backcross Charolais × Holstein cross population identified the presence of a highly significant QTL on chromosome 6 (BTA6) affecting the proportion of bone in the carcass. Two closely linked QTL affected birth weight (BW) and body length at birth (BBL). In this report, the marker density around the QTL on BTA6 was increased, adding four additional microsatellite markers across the chromosome and 46 SNPs within the target QTL confidence interval. Of the SNPs, 26 were in positional candidate genes and the remaining 20 provided an even distribution of markers in the target QTL region. As a bone-related trait, the sum of the bone weight for all the left fore- and hindquarter joints of the carcass was analysed. We also studied the BW and BBL. Analyses of the data substantially reduced the QTL confidence interval. No strong evidence was found that the QTL for the three traits studied are different, and we conclude that the results are consistent with a single pleiotropic QTL influencing the three traits, with the largest effects on the proportion of bone in the carcass. The analyses also suggest that none of the SNPs tested is the sole causative variant of the QTL effects. Specifically, the SNP in the NCAPG gene previously reported as a causal mutation for foetal growth and carcass traits in other cattle populations was excluded as the causal mutation for the QTL reported here. Polymorphisms located in other previously identified candidate genes including SPP1, ABCG2, IBSP, MEPE and PPARGC1A were also excluded. The results suggest that SNP51_BTA-119876 is the polymorphism in strongest linkage disequilibrium with the causal mutation(s). Further research is required to identify the causal variant(s) associated with this bone-related QTL.

Association, effects and validation of polymorphisms within the NCAPG - LCORL locus located on BTA6 with feed intake, gain, meat and carcass traits in beef cattle

BACKGROUND

In a previously reported genome-wide association study based on a high-density bovine SNP genotyping array, 8 SNP were nominally associated (P ≤ 0.003) with average daily gain (ADG) and 3 of these were also associated (P ≤ 0.002) with average daily feed intake (ADFI) in a population of crossbred beef cattle. The SNP were clustered in a 570 kb region around 38 Mb on the draft sequence of bovine chromosome 6 (BTA6), an interval containing several positional and functional candidate genes including the bovine LAP3, NCAPG, and LCORL genes. The goal of the present study was to develop and examine additional markers in this region to optimize the ability to distinguish favorable alleles, with potential to identify functional variation.

RESULTS

Animals from the original study were genotyped for 47 SNP within or near the gene boundaries of the three candidate genes. Sixteen markers in the NCAPG-LCORL locus displayed significant association with both ADFI and ADG even after stringent correction for multiple testing (P ≤ 005). These markers were evaluated for their effects on meat and carcass traits. The alleles associated with higher ADFI and ADG were also associated with higher hot carcass weight (HCW) and ribeye area (REA), and lower adjusted fat thickness (AFT). A reduced set of markers was genotyped on a separate, crossbred population including genetic contributions from 14 beef cattle breeds. Two of the markers located within the LCORL gene locus remained significant for ADG (P ≤ 0.04).

CONCLUSIONS

Several markers within the NCAPG-LCORL locus were significantly associated with feed intake and body weight gain phenotypes. These markers were also associated with HCW, REA and AFT suggesting that they are involved with lean growth and reduced fat deposition. Additionally, the two markers significant for ADG in the validation population of animals may be more robust for the prediction of ADG and possibly the correlated trait ADFI, across multiple breeds and populations of cattle.

Association of an ACSL1 gene variant with polyunsaturated fatty acids in bovine skeletal muscle

Background

The intramuscular fat deposition and the fatty acid profiles of beef affect meat quality. High proportions of unsaturated fatty acids are related to beef flavor and are beneficial for the nutritional value of meat. Moreover, a variety of clinical and epidemiologic studies showed that particularly long-chain omega-3 fatty acids from animal sources have a positive impact on human health and disease.

Results

To screen for genetic factors affecting fatty acid profiles in beef, we initially performed a microsatellite-based genome scan in a F₂ Charolais × German Holstein resource population and identified a quantitative trait locus (QTL) for fatty acid composition in a region on bovine chromosome 27 where previously QTL affecting marbling score had been detected in beef cattle populations. The long-chain acyl-CoA synthetase 1 (ACSL1) gene was identified as the most plausible functional and positional candidate gene in the QTL interval due to its direct impact on fatty acid metabolism and its position in the QTL interval. ACSL1 is necessary for synthesis of long-chain acyl-CoA esters, fatty acid degradation and phospholipid remodeling. We validated the genomic annotation of the bovine ACSL1 gene by in silico comparative sequence analysis and experimental verification. Re-sequencing of the complete coding, exon-flanking intronic sequences, 3' untranslated region (3'UTR) and partial promoter region of the ACSL1 gene revealed three synonymous mutations in exons 6, 7, and 20, six noncoding intronic gene variants, six polymorphisms in the promoter region, and four variants in the 3' UTR region. The association analysis identified the gene variant in intron 5 of the ACSL1 gene (c.481-233A>G) to be significantly associated with the relative content of distinct fractions and ratios of fatty acids (e.g., n-3 fatty acids, polyunsaturated, n-3 long-chain polyunsaturated fatty acids, trans vaccenic acid) in skeletal muscle. A tentative association of the ACSL1 gene variant with intramuscular fat content indicated that an indirect effect on fatty acid composition via modulation of total fat content of skeletal muscle cannot be excluded.

Conclusions

The initial QTL analysis suggested the ACSL1 gene as a positional and functional candidate gene for fatty acid composition in bovine skeletal muscle. The findings of subsequent association analyses indicate that ACSL1 or a separate gene in close proximity might play a functional role in mediating the lipid composition of beef.

From brain to food: analysis of phosphatidylcholins, lyso-phosphatidylcholins and phosphatidylcholin-plasmalogens derivates in Alzheimer's disease human post mortem brains and mice model via mass spectrometry

Alzheimer's disease (AD) is a devastating neurodegenerative disorder characterized by extracellular senile plaques mainly consisting of Aβ, a 40-42 amino acid long peptide, and intracellular neurofibrillary tangles, accompanied by an excessive loss of synapses. Recently evidence accumulated that nutrition, especially polyunsaturated fatty acids, influences AD pathogenesis. Especially mid-life food habits with the consumption of specific fatty acids (FA) appear to influence the disease risk. The timely separation between food intake and disease makes a direct correlation with detailed analysis of eating habits combined with accurate food analysis nearly unattainable. A possible solution to circumvent these difficulties is to investigate the FA composition in human post mortem brain. In this study we focused on the main phospholipids phosphatidylcholin (PC), phosphatidylcholin-plasmalogen (PC-PL) and lyso-phosphatidylcholin (lyso-PC) in AD brains compared to control brains. Frontal cortices, temporal cortices and cerebellum of 30 AD (mean 78 years) and 14 control aged matched brains (mean 77.4 years) as well as APP transgenic mice compared to control mice were analyzed using an AB Sciex 4000 Qtrap mass spectrometer utilizing a FIA MS/MS method. PC, PC-PL and lyso-PC metabolites were analyzed in respect to saturation level and FA composition. As expected, the majority of the lipid species showed no significant differences, but interestingly a few species revealed a highly significant reduction in AD brains. These FAs are potential candidates for further food analysis in respect to AD pathology. Additionally, we show that the method applied with multiple reaction monitoring (MRM) used for this study is suitable for semi quantitative analysis of small amounts (10 μl) of brain tissue.