Expression quantitative trait loci in sheep liver and muscle contribute to variations in meat traits

Background

Variants that regulate transcription, such as expression quantitative trait loci (eQTL), have shown enrichment in genome-wide association studies (GWAS) for mammalian complex traits. However, no study has reported eQTL in sheep, although it is an important agricultural species for which many GWAS of complex meat traits have been conducted. Using RNA sequence data produced from liver and muscle from 149 sheep and imputed whole-genome single nucleotide polymorphisms (SNPs), our aim was to dissect the genetic architecture of the transcriptome by associating sheep genotypes with three major molecular phenotypes including gene expression (geQTL), exon expression (eeQTL) and RNA splicing (sQTL). We also examined these three types of eQTL for their enrichment in GWAS of multi-meat traits and fatty acid profiles.

Results

Whereas a relatively small number of molecular phenotypes were significantly heritable (h² > 0, P < 0.05), their mean heritability ranged from 0.67 to 0.73 for liver and from 0.71 to 0.77 for muscle. Association analysis between molecular phenotypes and SNPs within ± 1 Mb identified many significant cis-eQTL (false discovery rate, FDR < 0.01). The median distance between the eQTL and transcription start sites (TSS) ranged from 68 to 153 kb across the three eQTL types. The number of common variants between geQTL, eeQTL and sQTL within each tissue, and the number of common variants between liver and muscle within each eQTL type were all significantly (P < 0.05) larger than expected by chance. The identified eQTL were significantly (P < 0.05) enriched in GWAS hits associated with 56 carcass traits and fatty acid profiles. For example, several geQTL in muscle mapped to the FAM184B gene, hundreds of sQTL in liver and muscle mapped to the CAST gene, and hundreds of sQTL in liver mapped to the C6 gene. These three genes are associated with body composition or fatty acid profiles.

Conclusions

We detected a large number of significant eQTL and found that the overlap of variants between eQTL types and tissues was prevalent. Many eQTL were also QTL for meat traits. Our study fills a gap in the knowledge on the regulatory variants and their role in complex traits for the sheep model.

A method for implementing methane breeding values in Australian dairy cattle

Context There has been a lot of interest in recent years in developing estimated breeding values (EBVs) to reduce methane emissions from the livestock sector. However, while a major limitation is the availability of high-quality methane phenotypes measured on individual animals required to develop these EBVs, it has been recognised that selecting for improved efficiency of milk production, longevity, feed efficiency and fertility may be an effective strategy to genetically reduce methane emissions in dairy cows. Aim Applying carbon dioxide equivalents (CO2-eq) weights to these EBVs, we hypothesise that it is possible to develop a genetic tool to reduce greenhouse-gas emissions (GHG). Methods We calculated the effect of an EBV unit change in each trait in the Balanced Performance Index on CO2-eq emissions per cow per year. The estimated environmental weights were used to calculate a prototype index of CO2-eq emissions. The final set of EBVs selected for inclusion in the GHG subindex were milk volume, fat yield and protein yield, survival and feed saved, as these traits had an independent effect on emissions. Feed saved is the Australian feed efficiency trait. A further modification was to include a direct methane trait in the GHG subindex, which is a more direct genomic evaluation of methane estimated from measured methane data, calculated as the difference between actual and predicted emissions, for example, a residual methane EBV. Key results The accuracy of the GHG subindex (excluding residual methane EBV) is ~0.50, calculated as the correlation between the index and gross methane (using 3-day mean gross methane phenotypes corrected for fixed effects, such as batch and parity and adjusting for the heritability). The addition of the residual methane EBV had a minimal effect with a correlation of 0.99 between the indexes. This was likely to be due to limited availability of methane phenotypes, resulting in residual methane EBVs with low reliabilities. Conclusions We expect that as more methane data becomes available and the accuracy of the residual methane trait increases, the two GHG subindexes will become differentiated. When the GHG subindex estimates are applied to bull EBVs, it can be seen that selecting for bulls that are low emitters of GHG can be achieved with a small compromise in the BPI of ~20 BPI units (standard deviation of BPI = 100). Implications Therefore, selection for more sustainable dairy cattle, both economic and environmental, may be promptly implemented until sufficient data are collected on methane.

Population structure and history of the Welsh sheep breeds determined by whole genome genotyping

Background

One of the most economically important areas within the Welsh agricultural sector is sheep farming, contributing around £230 million to the UK economy annually. Phenotypic selection over several centuries has generated a number of native sheep breeds, which are presumably adapted to the diverse and challenging landscape of Wales. Little is known about the history, genetic diversity and relationships of these breeds with other European breeds. We genotyped 353 individuals from 18 native Welsh sheep breeds using the Illumina OvineSNP50 array and characterised the genetic structure of these breeds. Our genotyping data were then combined with, and compared to, those from a set of 74 worldwide breeds, previously collected during the International Sheep Genome Consortium HapMap project.

Results

Model based clustering of the Welsh and European breeds indicated shared ancestry. This finding was supported by multidimensional scaling analysis (MDS), which revealed separation of the European, African and Asian breeds. As expected, the commercial Texel and Merino breeds appeared to have extensive co-ancestry with most European breeds. Consistently high levels of haplotype sharing were observed between native Welsh and other European breeds. The Welsh breeds did not, however, form a genetically homogeneous group, with pairwise F_ST between breeds averaging 0.107 and ranging between 0.020 and 0.201. Four subpopulations were identified within the 18 native breeds, with high homogeneity observed amongst the majority of mountain breeds. Recent effective population sizes estimated from linkage disequilibrium ranged from 88 to 825.

Conclusions

Welsh breeds are highly diverse with low to moderate effective population sizes and form at least four distinct genetic groups. Our data suggest common ancestry between the native Welsh and European breeds. These findings provide the basis for future genome-wide association studies and a first step towards developing genomics assisted breeding strategies in the UK.

Electronic supplementary material

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