- Research article
- Open Access
Adrenal cortex expression quantitative trait loci in a German Holstein × Charolais cross
BMC Geneticsvolume 17, Article number: 135 (2016)
The Erratum to this article has been published in BMC Genetics 2016 17:148
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.
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.
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.
eQTL analyses and the integration of proteomic and metabolomics data in genetic analyses contributed to the understanding of genes and molecular networks of complex traits in humans [1, 2] and farm animals [3–5]. In cattle most studies investigating the genetic background of complex traits such as milk performance or reproduction used quantitative trait loci (QTL) mapping and genome-wide association studies. Recent advances of high throughput technologies promoted the utilization of transcriptomic, proteomic or metabolomics data in genetic analyses .
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. Stress-related effects on productivity [7–9], health [10, 11] and reproduction [12, 13] in cattle are not fully understood. However, insight into molecular networks affecting these traits in cattle could contribute to the understanding of the mutual relationship between stress and obesity or stress related effects on health in humans.
Adrenal cortex is a promising target to investigate these effects, due to its impact on stress response and regulation of essential metabolic pathways . Human diseases such as Addison’s disease, the Cushing’s syndrome or Conn’s syndrome that are related to aberrations in glucocorticoid or mineralocorticoid levels further accentuate its relevance [15–17]. In cattle first observations indicating the importance of the adrenal gland for economically important traits were related to reproduction and lactation [18, 19]. Recent studies in cattle and pig showed that genetic variants in the glucocorticoid receptor and mineralocorticoid receptors are associated with cortisol concentrations and meat quality traits, for example [20, 21].
Therefore we investigated the genetic background of adrenocortical gene expression in a F2 cattle population (SEGFAM). The SEGFAM population is a cross of Holstein cows that were selected for high milk yield and Charolais bulls that were selected for meat production for a long time. By crossing the two breeds a high genetic and phenotypic variance in regard to milk and meat production traits was expected due to segregation of different QTL alleles in the F2 cross . The population was established to investigate the genetic and physiological background of nutrient transformation types, and first studies showed a huge phenotypic variance in regard to growth at the onset of puberty [4, 23], milk yield  and temperament [25, 26]. By identifying eQTL in this population we assumed to gain insight into the genetic regulation of gene expression of different nutrition types that could provide new targets for further research of genes that have an effect on economically important traits such as milk production and meat quality traits in cattle.
The aims of this study were first to identify eQTL in the adrenal cortex in cattle and secondly, to evaluate whether associated genes and molecular pathways are known to affect economically important traits in cattle. The analysed data comprised genotypes of 37,204 SNPs and expression values for 10,986 transcripts. In order to account for systemic effects and to de-correlate the expression data a sire-dam model was used during pre-processing. To account for relatedness the complete pedigree data comprising the founder and F1 generation was used. Then we performed a genome-wide eQTL analysis using an adaptive linear rank test  on the residuals of the sire-dam model and the genotype data. Finally, we compared the results of this study with previous studies within the same population and performed Ingenuity® Pathway Analysis (IPA, QIAGEN Redwood City, www.qiagen.com/ingenuity) and literature reviews to identify the molecular pathways involved and evaluate potential targets for further research.
Expression quantitative trait loci
In total, 505 SNPs were identified to be associated with the expression of 129 transcripts. A complete list of eQTL including the chromosomal position and annotation of the SNP and the gene encoding the transcript is provided in Additional file 1. In addition to Benjamini-Hochberg adjusted p-values of the adaptive linear rank test we also included the unadjusted p-values of adaptive linear rank test, Kruskal Wallis-test and ANOVA for completeness. The number of eQTL per chromosome is shown in Fig. 1. Identified associations and the corresponding genomic position for both SNP and the gene encoding the transcript are shown in Fig. 2. eQTL were identified on all chromosomes except on chromosomes X, 16, 24 and 28. The highest number of associations was observed on chromosome 23 comprising 55 SNPs associated with 11 transcripts. Apparently, in the case of most of these associations (482 of 505) both the SNP and the gene encoding the transcript were located on the same chromosome (cis associations). Only for 17 associations SNPs and genes were located on different chromosomes (trans associations). For the remaining six associations the position of the SNP was unknown.
Genomic distances between SNPs and genes for all 482 cis associations are shown in Fig. 3. Most of these SNPs (343 of 482) were located in close vicinity (less than 5 MB) to the gene encoding the transcript. On the other hand, the expression of LYPD6B (Chr. 2), GRIP1 (Chr. 5), SNCB (Chr. 7), COQ3 (Chr. 9) and ALDH5A1 (Chr.23), for example, was associated with 21, 15, 35, 28 and 20 SNPs, respectively. These SNPs were located up to 30 MB up- or downstream of the gene encoding the transcript. Considering that we investigated a F2 population that showed considerable linkage disequilibrium ( and Additional file 2) a clear discrimination between local and distant or cis and trans eQTL was hindered for these eQTL [29, 30]. Therefore, we assumed a cis or chromosome specific effect for all eQTL for which the SNP and the associated gene were located on the same chromosome. For the 17 eQTL for which the SNP and the gene encoding the transcript were located on different chromosomes a trans-regulatory effect was assumed.
Ingenuity® pathway analysis
In order to gain insight into molecular and cellular functions and regulatory processes, in which the genes encoding the 129 identified transcripts are involved in, Ingenuity® pathway analysis was performed. The analysis indicated that the genes are mainly involved in molecular and cellular function comprising free radical scavenging, cellular compromise, cell morphology, lipid metabolism and small molecule biochemistry. A complete list of the top five diseases and bio functions including the related genes is provided in Table 1. The results of the Ingenuity® Pathway upstream analyses are provided in Tables 2 and 3. Prolactin was the most important upstream regulator. The most important causal networks were related to ESRRG, GNRH1 and the two retinoic acid receptors RARG and RARB as master regulators (Table 3). Most of the affected diseases and bio functions as well as the most important causal networks and the most important upstream regulators are related to a relatively small number of genes that are involved in many of the identified pathways and causal networks, such as CYCS, LHCGR, TF, VIM, XDH, CYP27A1 and NQO1.
In order to identify molecular pathways that potentially influence the phenotypic variation of economically important traits in cattle we investigated the genetic background of adrenocortical gene expression in an experimental F2 population of a Charolais × German Holstein cross. We identified 505 SNPs associated with the expression of 129 genes. We also showed that some of these genes are involved in molecular pathways affecting economically important traits in cattle. In comparison with eQTL mapping studies in other species the proportion of cis eQTL (482) compared with trans eQTL (17) observed in our study is high. This is in line with the observation that the number of identified local or cis eQTL is related to sample size, extent of linkage disequilibrium, size of the window applied for cis eQTL and effect size of eQTL which is assumed to be stronger for cis eQTL in comparison with trans eQTL [29–31]. A reason for the high number of cis eQTL observed in our study might be that we assumed a cis or chromosome specific effect for all associations for which SNP and gene were located on the same chromosome. We decided to use this definition for cis associations for two reasons. Firstly, the extent of linkage disequilibrium observed in the studied population is vast ( and Additional file 2) and, secondly, the relative genomic position of SNP and associated gene does not necessarily provide information about the underlying regulatory variation . However, even a window size of 1 MB which is often used to differentiate cis and trans effects would result in a higher number of cis eQTL (314) in comparison with trans eQTL (191). In addition, the extent of linkage disequilibrium and the high number of transcripts that are associated with more than one SNP in close vicinity to the gene encoding the transcript also indicated that the identified SNPs are not necessarily the causal variants but are in linkage disequilibrium to the causal polymorphisms and that the identified associations rather indicate genomic regions harbouring the causal variants affecting the transcript abundance.
Ingenuity® pathway analyses of genes encoding affected transcripts revealed that some of the genes are involved in “Diseases and Bio Functions” like lipid metabolism, inflammatory response and reproductive system development and function. This indicates that the identified eQTL potentially contribute to the phenotypic variation of economically important traits in cattle. For some of the genes there is evidence that they have an effect on reproduction, milk performance, growth or health traits, based on genetic and molecular studies in cattle and other farm animals. The following sections comprise details for selected genes. In addition, comparisons with studies within the SEGFAM population indicate that some of the genes affected by eQTL contribute to the phenotypic variation of behaviour, temperament and growth at the onset of puberty. Since the effects of adrenocortical gene expression on most of the economically important traits in cattle still have to be investigated, for other genes only the contribution to disease phenotypes in human or mice or the involvement in specific metabolic pathways like reproductive system development and function suggests that they might contribute to the phenotypic variation of economically important traits in cattle. In addition, the sheer existence of genomic regions affecting the adrenocortical transcript abundance of genes that are known to influence economically important traits in cattle is no evidence that the adrenocortical transcript abundance is causal, but similar effects might be present in other organs and tissues  providing insight into molecular mechanisms that contribute to the phenotypic variation in cattle. Subsequently we will discuss some of the genes that were highlighted by the functional analyses and that have been discussed to be involved in the phenotypic manifestation in cattle in other studies.
Genes involved in reproduction and other traits
Vimentin (VIM) is a type III intermediate filament that is expressed in cells of mesenchymal origin and has key functions in the formation of the cytoskeletal network, organelle positioning, cell migration and adhesions as well as in cell signaling . VIM expression has been shown to be required for the development and establishment of the embryo in cattle. It was reported to be differentially expressed between nuclear transfer and by in vitro fertilization produced embryos [33, 34]. Nowadays in vitro production of embryos is common and approximately 15% of all embryos are produced by in vitro technologies in cattle [35, 36].
The luteinizing hormone/chorigonadotropin receptor (LHCGR) is a G-protein-coupled receptor that is mainly known for binding luteinizing hormone and chorionic gonadotropin and mediating their cellular actions. In humans inactivating mutations in the LHCGR gene have been shown to cause infertility in both male and female . In addition, LHCGR is discussed to be involved in the manifestation of the pregnancy induced Cushing’s Syndrome due to adrenocortical LHCGR overexpression that may lead to adrenocortical hyperfunction [38–40]. In cattle, polymorphisms within the LHCGR gene have been shown to be associated with calving interval, days to first service and milk yield  as well as with superovulation traits [42–44] that are important for generating a high number of eggs for the production of embryos [35, 36].
The cytochrome P450, family 27, subfamily A, polypeptide 1 (Cyp27A1) is a key-enzyme in Vitamin D3 metabolism  and is required for bile acid synthesis [46, 47]. In cattle, Vitamin D3 supplementation has been shown to affect dry-matter intake and beef tenderness . In a previous study in the SEGFAM population, metabolites that belong to the sub class “Vitamine D3 and derivates” of Sterol Lipids were important for the classification of temperament types based on prefrontal cortex and serum metabolites . In addition, likewise to LHCGR, Cyp27A1 has also been discussed to be important for reproduction  but has also been shown to affect milk yield and somatic cell score , an important production and an important health trait in dairy cattle, respectively.
Genes involved in milk production and other traits
The xanthine dehydrogenase (XDH) is a xanthine oxidoreductase that belongs to the group of molybdenum iron-sulfur flavin hydroxylases. XDH contributes to the detoxification of endogenous or xenobiotic compounds and is involved in the oxidative metabolism of purines [51, 52]. In cattle XDH has been shown to be upregulated during the lactation cycle due to its contribution to lipid droplet formation . In addition, XDH protein abundance has been shown to increase in the first nine days of lactation in milk serum , and  suggested XDH to be a candidate gene for milk production and mastitis susceptibility based on studies in mouse and cattle.
Transferrin (TF) together with ferritin is important for acquisition, transport and storage of iron which is essential for many metabolic processes like the synthesis of hemoglobin for oxygen transport . In cattle, TF has been shown to be highly polymorphic  and a study in Chinese native cattle indicated that TF polymorphisms are associated with protein yield, 305-day milk yield and mastitis susceptibility. In addition, they showed that TF mRNA expression was higher in mastitis affected in comparison with unaffected mammary tissue .
Genes potentially affecting behavior and temperament
The aldehyde dehydrogenase 5 family, member A1 (ALDH5A1) is a semialdehyde dehydrogenase that is important for the metabolism of γ-aminobutyric acid (GABA) the major inhibitory neurotransmitter in the brain . A differential expression of ALDH5A1 has been reported in Angus cattle selected for high or low residual feed intake . Increased levels of 2,4-dihidroxy-butanoic acid have been observed in the urine of patients suffering from semialdehyde dehydrogenase deficiency  and an accumulation of gamma-hydroxybutyric acid (GHB) and GABA in the central nervous system was reported in ALDH5A1−/− mice . 2,4-dihidroxy-butanoic was priviously identified to be a potential biomarker in Alzheimer’s disease  and a temperament type related abundance of 2,4-dihidroxy-butanoic acid as well as of γ-aminobutyric acid acid has also been observed in temperament type specific metabolite profiles of the prefrontal cortex of animals with extrem temperament types deriving from the SEGFAM population . Interestingly ALDH5A1 is located on chromosome 23 in the genomic region were bovine MHC genes are located. In this genomic region we identified a high number of SNPs that were not only associated with the expression of ALDH5A1 but also with the expression of some MHC genes like the MHC class I heavy chain BOLA and the major histocompatibility complex BOLA-DQA1. This genomic region has been reported to be associated with the antibody mediated immune response in Canadian Holstein cows .
Another study analyzed the impact of genetic variants on behaviour characteristics assessed in an open field and novel object test within the SEGFAM population . We identified three SNPs (ARS-BFGL-NGS-2942/rs110027993, ARS-BFGL-NGS-98658/rs109674592, ARS-BFGL-NGS-27299/rs109313646) that affect the expression of NME6 and RAB32 in this study. The same SNPs were previously identified to be associated with the inactivity in an open field test and the explorative behaviour in a novel object test, respectively .
The NME/NM23 nucleoside diphosphate kinase 6 (NME6) is a nucleoside diphosphate kinase. Little is known about its function besides that it is assumed to play a role in cell growth and cell cycle progression [65, 66]. In addition, a recent study in pig indicated that NME6 is located in a QTL region associated with back fat .
RAB32, member RAS oncogene family (RAB32) belongs to the RAB family of small GTP-binding proteins. It is involved in the biogenesis of lysosome-related organelles like melanosomes . Similar to NME6 we could not find any additional information other than a suggested role in lipid metabolism .
Ingenuity® upstream analysis
The upstream analysis revealed PRL, ESRRG and GNRH1 as upstream or master regulators, providing insight into regulatory mechanisms that are potentially affected by eQTL.
Prolactin (PRL) is an anterior pituitary protein hormone that is not only important for reproduction but is also involved in the control of behaviours . In cattle especially the effects of PRL in regard to lactation and milk production have been investigated. In the last years several studies in cattle investigated the inhibition of PRL-release on milk yield and showed that the inhibition was associated with a decrease in milk yield additionally indicating an important role of PRL in the maintenance of lactation [71–73]. The SEGFAM population showed remarkable differences in their milk production .
Gonadotropin-releasing hormone 1 (GNRH1) encodes a protein that is secreted and cleaved to form the luteinizing hormone-releasing hormone and the prolactin release-inhibiting factor that regulate LH and FSH gene expression . LH and FSH both have been early recognized to be involved in the onset of puberty in cattle and today Gonadotropin-releasing hormone is used to synchronize cows for timed artificial insemination especially in beef cattle [75, 76]. In addition, GnRH signaling was previously identified as a key modifier of differential growth at the onset of puberty in the SEGFAM population .
The Estrogen-Related Receptor Gamma (ESRRG) that belongs to the estrogen receptor-related receptor family has recently been shown to be important in the fibrate-mediated regulation of lipid metabolism genes in a human ApoA-I transgenic mice model . Interestingly ESRRG was also identified to be associated with the onset of puberty in cattle .
In the present study we investigated the genetic background of adrenocortical gene expression and identified eQTL affecting the expression levels of genes that have been reported to influence economically important traits in cattle. In addition, some of the genes and molecular pathways affected by eQTL were previously shown to contribute to the phenotypic variation of behaviour, temperament and growth at the onset of puberty in the SEGFAM population. Therefore, eQTL analysis appears to be a useful approach to gain insight into the molecular and genetic background of complex traits in cattle. Additional studies are needed to identify the causal polymorphisms and underlying molecular mechanisms (cis or trans) and the effects of genetic variation on gene expression in other tissues.
Animals and tissue samples
This study comprised 145 F2 cows from the experimental resource cross population (SEGFAM) that was initiated at the Leibniz Institute for Farm Animal Biology (FBN) in the P0 generation by crossing five purebred Charolais sires to purebred German Holstein cows . The population was established by multiple ovulation and embryo transfer, which resulted in five half-sib and full-sib families from which 49, 35, 52, 8 and 1 cows were included in this study, respectively. The cows from the fourth and fifth family were maternally related to animals in the first and third family, respectively. All animals were housed in a loose housing barn under identical environmental and feeding conditions at the Leibniz Institute for Farm Animal Biology (FBN) in Dummerstorf, Germany . The animals were slaughtered at day 30 after parturition in their second lactation following a standardized protocol. Besides other tissues mammary gland tissue and the adrenal gland were immediately taken after slaughter for DNA and RNA extraction, respectively. The adrenal gland was further dissected to separate the adrenal cortex from the medulla. All tissue samples were cut in small pieces, snap-frozen in liquid nitrogen and stored at −80 °C or in liquid nitrogen.
DNA extraction and SNP genotyping
SNP genotyping was accomplished using Illumina® BovineSNP50 Beadchip v1.0 and v2.0. Genotype calling was performed with the Genotyping Module of the GenomeStudio V2011.1 Software (Illumina®). For this purpose genomic DNA was extracted from mammary gland tissue with the QIAmp DNA Mini Kit (QIAGEN, Hilden, Germany) and prepared for hybridization on the Illumina® BovineSNP50 Beadchip following manufacturer’s instructions. Genotype calling and quality control was performed according to  for v1.0 and v2.0 data separately. After quality control and filtering the v1.0 and v2.0 datasets were merged and all samples and SNPs with more than 10% missing genotypes and a minor allele frequency of less than 5% were excluded from further analyses. To identify inconsistencies between recorded genotypes and pedigree information, the software PedCheck  and a larger dataset including genotypes of P0 sires and all (male and female) F1 and F2 animals was used . The final dataset comprised 37,204 SNPs of 145 animals. Chromosomal positions of SNPs are taken from  and are based on genome assembly UMD3.1.
RNA extraction and microarray hybridization
Total RNA was isolated from the adrenal cortex of 145 animals using TRI Reagent (Sigma, Taufkirchen, Germany). After DNase I treatment the RNA was further purified with the RNeasy Kit (QIAGEN, Hilden, Germany) following the manufacturers recommendations. RNA was quantified using the NanoDrop ND-1000 spectrophotometer (Peqlab, Erlangen, Germany). Integrity of RNA was checked by running 1 μg of RNA on a 1% agarose gel. Absence of DNA contamination was verified by PCR amplifying fragments of GAPDH with RNA as a template. For hybridization, 500 ng of total RNA were amplified using Ambion WT Expression Kit (Affymetrix). Subsequently, 5.5 μg of the resulting cDNA was fragmented and labeled using the Affymetrix Terminal Labeling Kit. The fragmented cDNA was hybridized to the microarray using the Affymetrix Hyb-WashStain Kit and Affymetrix standard protocols. Fluidic station protocol was FS450_0001.
For expression profiling the custom GeneChip Bovine Gene v1 Array was used. Affymetrix designed the GeneChip Bovine Gene v1 Array based on Ensemble and RefSeq predictions for the Genome Bos Taurus Build 4.0. The design was targeted to develop a whole Genome Expression Array with approximately 25 probes per transcript distributed over the whole length of each transcript. In total 194,712 probe sets targeting at approximately 24,000 transcripts are implemented on the array. In addition standard Affymetrix controls for hybridization and labeling efficiency as well as for non-specific binding are included on the array. Gene-level analysis was performed using Affymetrix® Expression Console™ Software. Microarray raw data were preprocessed using the RMA algorithm . Quality control was performed in accordance to . In addition, the detection above background algorithm was used to identify expressed genes. Probe sets that were present in 75% of all samples and transcripts with at least 50% of all probe sets present were included in the analyses. The final dataset comprised expression values for 10,986 transcripts of 145 animals. For annotation the NetAffx  release 34 chromosomal positions of the transcripts were used that are based on genome assembly UMD3.1. In addition the Ingenuity® gene annotation is provided including information about transcript cluster ids that target the same gene. The microarray data are deposited at Gene Expression Omnibus database  (GEO: GSE75371).
After quality control and filtering the expression data was further pre-processed to account for systemic effects. Furthermore, the statistical test used for eQTL analysis required the data to be decorrelated. Since we expected gametic phase disequilibrium for some of the loci, we refrained from estimating kinship from SNP data but preferred the use of kinship estimates from pedigree information. Therefore the following sire-dam model was applied in ASReml  to account for fixed effects of season and year and age at slaughter as well as for relatedness:
where y i is the expression value for transcript y of animal i (i = 1,…, 145), s j is the fixed effect of year and season (j = 1,…, 25) at slaughter, x i is the age in days of animal i at slaughter, a if is the additive-genetic effect of the sire of animal i, a im is the additive-genetic effect of the dam of animal i and ε i is the random residual effect. To estimate the additive-genetic effects the complete pedigree data comprising 38 P0 and 52 F1 animals was used. The residuals comprise the Mendelian sampling effect as an uncorrelated genetic component as well as a random deviation. For further analysis the residuals were considered as de-correlated expression data.
For eQTL analysis an adaptive linear rank test  was applied. The test was optimized to meet the challenge of non-normally distributed expression data in eQTL studies. This two stage procedure first calculates selector statistics for skewness and tail length of the de-correlated expression data. Then an appropriate linear rank test depending on the selector statistics is applied to identify differences in the expression levels per genotype group. Thus, for each pair of SNP and transcript a linear rank test is selected that is supposed to have large statistical power for the respective actual distribution. The available tests are Kruskal-Wallis test , median test , long tails test , short tails test , right skewness and left skewness test . To account for multiple testing we applied a global significance level of 0.1 for Benjamini-Hochberg adjusted p-values  (considering all 374 million performed tests), which is approximately equivalent to a significance level of 1.3x10−7 for unadjusted p-values. This global false discovery rate (FDR) threshold of 0.1 is less conservative for trans associations than usual, but it is very conservative for cis associations, since it corresponds to a FDR threshold of approximately 0.0035 for an eQTL analysis in which only cis associations (approximately 13 million) are considered for the calculation of the Benjamini-Hochberg adjusted p-values. This was done because the majority of associations was expected between SNPs and transcripts for which the encoding gene of the transcript was located on the same chromosome in close vicinity to the SNP.
Ingenuity® pathway analysis
The functional and upstream analyses were performed through the use of QIAGEN’s Ingenuity® Pathway Analysis (IPA, QIAGEN Redwood City, www.qiagen.com/ingenuity). The purpose of the functional analysis was to determine if the 129 identified genes are statistically significant over-represented (more genes than expected by chance) for a specific molecular or cellular function in comparison to a reference-set of genes. Similar to the functional analysis, the Ingenuity® upstream analysis was performed to measure whether there is a statistically significant overlap between the 129 identified genes and the genes that are regulated by a specific transcription regulator. For both tests IPA uses a Fisher’s Exact Test to determine the significance of the over-representation or overlap in comparison to a reference set. As reference set for the functional analysis the Ingenuity® Knowledge Base was chosen, because the reference network that is used to identify the upstream and master (upstream and downstream) regulators in the upstream analysis is based on the Ingenuity® Knowledge Base. Reported p-values are Benjamini-Hochberg adjusted.
Analysis of variance
Expression quantitative trait loci
False discovery rate
Polymerase chain reaction
Quantitative trait loci
Single nucleotide polymorphism
Cazier JB, Kaisaki PJ, Argoud K, Blaise BJ, Veselkov K, Ebbels TMD, et al. Untargeted metabolome quantitative trait locus mapping associates variation in urine glycerate to mutant glycerate kinase. J Proteome Res. 2012;11:631–42.
Portelli MA, Siedlinski M, Stewart CE, Postma DS, Nieuwenhuis MA, Vonk JM, et al. Genome-wide protein QTL mapping identifies human plasma kallikrein as a post-translational regulator of serum uPAR levels. FASEB J. 2014;28:923–34.
Ponsuksili S, Murani E, Brand B, Schwerin M, Wimmers K. Integrating expression profiling and whole-genome association for dissection of fat traits in a porcine model. J Lipid Res. 2011;52:668–78.
Widmann P, Reverter A, Fortes MR, Weikard R, Suhre K, Hammon H, et al. A systems biology approach using metabolomic data reveals genes and pathways interacting to modulate divergent growth in cattle. BMC Genomics. 2013;14:798.
Yang Y, Bu D, Zhao X, Sun P, Wang J, Zhou L. Proteomic Analysis of Cow, Yak, Buffalo, Goat and Camel milk whey proteins: quantitative differential expression patterns. J Proteome Res. 2013;12:1660–7.
Kadarmideen HN. Genomics to systems biology in animal and veterinary sciences: progress, lessons and opportunities. Livest Sci. 2014;166:232–48.
Hemsworth PH, Coleman GJ, Barnett JL, Borg S. Relationships between human-animal interactions and productivity of commercial dairy cows. J Anim Sci. 2000;78:2821–31.
Müller R, von Keyserlingk MAG. Consistency of flight speed and its correlation to productivity and to personality in Bos taurus beef cattle. Appl Anim Behav Sci. 2006;99:193–204.
Sutherland MA, Rogers AR, Verkerk GA. The effect of temperament and responsiveness towards humans on the behavior, physiology and milk production of multi-parous dairy cows in a familiar and novel milking environment. Physiol Behav. 2012;107:329–37.
Burdick N, Randel R, Carroll J, Welsh T. Interactions between temperament, stress, and immune function in cattle. Int J Zoonoses. 2011;2011:1–9.
Burdick NC, Banta JP, Neuendorff DA, White JC, Vann RC, Laurenz JC, et al. Interrelationships among growth, endocrine, immune, and temperament variables in neonatal Brahman calves. J Anim Sci. 2009;87:3202–10.
Kasimanickam R, Asay M, Schroeder S, Kasimanickam V, Gay JM, Kastelic JP, et al. Calm temperament improves reproductive performance of beef cows. Reprod Domest Anim. 2014;49:1063–7.
Cooke RF, Arthington JD, Araujo DB, Lamb GC. Effects of acclimation to human interaction on performance, temperament, physiological responses, and pregnancy rates of Brahman-crossbred cows. J Anim Sci. 2009;87:4125–32.
Nussey S, Whitehead S. The adrenal gland. In: Endocrinology: an integrated approach. Oxford: BIOS Scientific Publishers; 2001.
Boscaro M, Barzon L, Fallo F, Sonino N. Cushing's syndrome. Lancet. 2001;357:783–91.
Piaditis G, Markou A, Papanastasiou L, Androulakis II, Kaltsas G. Progress in aldosteronism: a review of the prevalence of primary aldosteronism in pre-hypertension and hypertension. Eur J Endocrinol. 2015;172:R191–203.
Ten S, New M, Maclaren N. Addison`s Disease 2001. J Clin Endocrinol Metab. 2001;86:2909–22.
Wagner WC, Oxenreider SL. Adrenal function in the cow. Diurnal changes and the effects of lactation and neurohypophyseal hormones. J Anim Sci. 1972;34:630–5.
Wagner WC, Hansel W. Reproductive physiology of the post partum cow. J Reprod Fertil. 1969;18:493–500.
Poleti MD, DeRijk RH, Rosa AF, Moncau CT, Oliveira PS, Coutinho LL, et al. Genetic variants in glucocorticoid and mineralocorticoid receptors are associated with concentrations of plasma cortisol, muscle glycogen content, and meat quality traits in male Nellore cattle. Domest Anim Endocrinol. 2015;51:105–13.
Reyer H, Ponsuksili S, Wimmers K, Murani E. Association of N-terminal domain polymorphisms of the porcine glucocorticoid receptor with carcass composition and meat quality traits. Anim Genet. 2014;45:125–9.
Kühn C, Bellmann O, Voigt J, Wegner J, Guiard V, Ender K. An experimental approach for studying the genetic and physiological background of nutrient transformation in cattle with respect to nutrient secretion and accretion type. Arch Tierz. 2002;45:317–30.
Weikard R, Altmaier E, Suhre K, Weinberger K, Hammon H, Albrecht E, et al. Metabolomic profiles indicate distinct physiological pathways affected by two loci with major divergent effect on Bos taurus growth and lipid deposition. Physiol Genomics. 2010;42A:79–88.
Hammon HM, Metges CC, Schulz A, Junghans P, Steinhoff J, Schneider F, et al. Differences in milk production, glucose metabolism, and carcass composition of 2 Charolais X Holstein F2 families derived from reciprocal paternal and maternal grandsire crosses. J Dairy Sci. 2010;93:3007–18.
Brand B, Hadlich F, Brandt B, Schauer N, Graunke KL, Langbein J, et al. Temperament type specific metabolite profiles of the prefrontal cortex and serum in cattle. PLoS One. 2015;10:e0125044.
Graunke LK, Nürnberg G, Repsilber D, Puppe B, Langbein J. Describing temperament in an ungulate: a multidimensional approach. PLoS One. 2013;8:e74579.
Szymczak S, Scheinhardt MO, Zeller T, Wild PS, Blankenberg S, Ziegler A. Adaptive linear rank tests for eQTL studies. Stat Med. 2013;32:524–37.
Friedrich J, Brand B, Ponsuksili S, Graunke KL, Langbein J, Knaust J, et al. Detection of genetic variants affecting cattle behaviour and their impact on milk production: a genome-wide association study. Anim Genet. 2015;47:12–8.
Rockman MV, Kruglyak L. Genetics of global gene expression. Nat Rev Genet. 2006;7:862–72.
Williams RBH, Chan EKF, Cowley MJ, Little PFR. The influence of genetic variation on gene expression. Genome Res. 2007;17:1707–16.
Holloway B, Luck S, Beatty M, Rafalski JA, Li B. Genome-wide expression quantitative trait loci (eQTL) analysis in maize. BMC Genomics. 2011;12:336.
Ivaska J, Pallari HM, Nevo J, Eriksson JE. Novel functions of vimentin in cell adhesion, migration, and signaling. Exp Cell Res. 2007;313:2050–62.
Pfister-Genskow M, Myers C, Childs LA, Lacson JC, Patterson T, Betthauser JM, et al. Identification of differentially expressed genes in individual bovine preimplantation embryos produced by nuclear transfer: improper reprogramming of genes required for development. Biol Reprod. 2005;72:546–55.
Maddox-Hyttel P, Alexopoulos NI, Vajta G, Lewis I, Rogers P, Cann L, et al. Immunohistochemical and ultrastructural characterization of the initial post-hatching development of bovine embryos. Reproduction. 2003;125:607–23.
Mapletoft RJ, Hasler JF. Assisted reproductive technologies in cattle: a review. Rev Sci Tech. 2005;24:393.
Rodriguez-Martinez H. Assisted reproductive techniques for cattle breeding in developing countries: a critical appraisal of their value and limitations. Reprod Domest Anim. 2012;47:21–6.
Layman LC. Human gene mutations causing infertility. J Med Genet. 2002;39:153–61.
Chui MH, Ozbey NC, Ezzat S, Kapran Y, Erbil Y, Asa S. Case report: Adrenal LH/hCG receptor overexpression and gene amplification causing pregnancy-induced Cushing's syndrome. Endocr Pathol. 2009;20:256–61.
Bernichtein S, Alevizaki M, Huhtaniemi I. Is the adrenal cortex a target for gonadotropins? Trends Endocrinol Metab. 2008;19:231–8.
Carlson HE. Human adrenal cortex hyperfunction due to LH/hCG. Mol Cell Endocrinol. 2007;269:46–50.
Hastings N, Donn S, Derecka K, Flint AP, Woolliams JA. Polymorphisms within the coding region of the bovine luteinizing hormone receptor gene and their association with fertility traits. Anim Genet. 2006;37:583–5.
Yang WC, Tang KQ, Li SJ, Chao LM, Yang LG. Polymorphisms of the bovine luteinizing hormone/choriogonadotropin receptor (LHCGR) gene and its association with superovulation traits. Mol Biol Rep. 2012;39:2481–7.
Luo W, Gumen A, Haughian JM, Wiltbank MC. The role of luteinizing hormone in regulating gene expression during selection of a dominant follicle in cattle. Biol Reprod. 2011;84:369–78.
Yu Y, Pang Y, Zhao H, Xu X, Wu Z, An L, et al. Association of a missense mutation in the luteinizing hormone/choriogonadotropin receptor gene (LHCGR) with superovulation traits in Chinese Holstein heifers. J Anim Sci Biotechnol. 2012;3:1–5.
Sawada N, Sakaki T, Ohta M, Inouye K. Metabolism of vitamin D3 by human CYP27A1. Biochem Biophys Res Commun. 2000;273:977–84.
Goodwin B, Gauthier KC, Umetani M, Watson MA, Lochansky MI, Jon LC, et al. Identification of bile acid precursors as endogenous ligands for the nuclear xenobiotic pregnane X receptor. Proc Natl Acad Sci U S A. 2003;100:223–8.
Chiang JYL. Bile acids: regulation of synthesis. J Lipid Res. 2009;50:1955–66.
Karges K, Brooks JC, Gill DR, Breazile JE, Owens FN, Morgan JB. Effects of supplemental vitamin D3 on feed intake, carcass characteristics, tenderness, and muscle properties of beef steers. J Anim Sci. 2001;79:2844–50.
Bauersachs S, Ulbrich SE, Zakhartchenko V, Minten M, Reichenbach M, Reichenbach HD, et al. The endometrium responds differently to cloned versus fertilized embryos. Proc Natl Acad Sci U S A. 2009;106:5681–6.
Pimentel ECG, Bauersachs S, Tietze M, Simianer H, Tetens J, Thaller G, et al. Exploration of relationships between production and fertility traits in dairy cattle via association studies of SNPs within candidate genes derived by expression profiling. Anim Genet. 2011;42:251–62.
Xu P, Huecksteadt TP, Hoidal JR. Molecular cloning and characterization of the human xanthine dehydrogenase gene (XDH). Genomics. 1996;34:173–80.
Weseler AR, Bast A. Oxidative stress and vascular function: implications for pharmacologic treatments. Curr Hypertens Rep. 2010;12:154–61.
Bionaz M, Loor J. Gene networks driving bovine milk fat synthesis during the lactation cycle. BMC Genomics. 2008;9:366.
Zhang L, Boeren S, Hageman JA, van Hooijdonk T, Vervoort J, Hettinga K. Bovine milk proteome in the first 9 days: protein interactions in maturation of the immune and digestive system of the newborn. PLoS One. 2015;10:e0116710.
Ogorevc J, Kunej T, Razpet A, Dovc P. Database of cattle candidate genes and genetic markers for milk production and mastitis. Anim Genet. 2009;40:832–51.
Ponka P, Beaumont C, Richardson DR. Function and regulation of transferrin and ferritin. Semin Hematol. 1998;35:35–54.
Lukac D, Vidovic V, Nemeš Z, Stupar M, Popovic-Vranješ A. Genotypic frequencies of the ß-lactoglobulin, k-casein and transferrin in Serbian Holstein-Friesian dairy cattle. Mljekarstvo/Dairy. 2013;63:203–10.
Ju Z, Li Q, Huang J, Hou M, Li R, Li J, et al. Three novel SNPs of the bovine Tf gene in Chinese native cattle and their associations with milk production traits. Genet Mol Res. 2011;10:340–52.
McCormick DA. GABA as an inhibitory neurotransmitter in human cerebral cortex. J Neurophysiol. 1989;62:1018–27.
Chen Y, Gondro C, Quinn K, Herd RM, Parnell PF, Vanselow B. Global gene expression profiling reveals genes expressed differentially in cattle with high and low residual feed intake. Anim Genet. 2011;42:475–90.
Brown GK, Cromby CH, Manning NJ, Pollitt RJ. Urinary organic acids in succinic semialdehyde dehydrogenase deficiency: evidence of a-oxidation of 4-hydroxybutyric acid, interaction of succinic semialdehyde with pyruvate dehydrogenase and possible secondary inhibition of mitochondrial ß-oxidation. J Inherit Metab Dis. 1987;10:367–75.
Knerr I, Pearl PL, Bottiglieri T, Carter Snead O, Jakobs C, Gibson KM. Therapeutic concepts in succinate semialdehyde dehydrogenase (SSADH; ALDH5a1) deficiency (gamma-hydroxybutyric aciduria). Hypotheses evolved from 25 years of patient evaluation, studies in Aldh5a1−/− mice and characterization of gamma-hydroxybutyric acid pharmacology. J Inherit Metab Dis. 2007;30:279–94.
Oresic M, Hyotylainen T, Herukka SK, Sysi-Aho M, Mattila I, Seppanan-Laakso T, et al. Metabolome in progression to Alzheimer's disease. Transl Psychiatry. 2011;1:e57.
Thompson-Crispi K, Sargolzaei M, Ventura R, Abo-Ismail M, Miglior F, Schenkel F, et al. A genome-wide association study of immune response traits in Canadian Holstein cattle. BMC Genomics. 2014;15:559.
Tsuiki H, Nitta M, Furuya A, Hanai N, Fujiwara T, Inagaki M, et al. A novel human nucleoside diphosphate (NDP) kinase, Nm23-H6, localizes in mitochondria and affects cytokinesis. J Cell Biochem. 2000;76:254–69.
Desvignes T, Pontarotti P, Fauvel C, Bobe J. Nme protein family evolutionary history, a vertebrate perspective. BMC Evol Biol. 2009;9:256.
Onteru SK, Gorbach DM, Young JM, Garrick DJ, Dekkers JCM, Rothschild MF. Whole genome association studies of residual feed intake and related traits in the pig. PLoS One. 2013;8:e61756.
Bultema JJ, Boyle JA, Malenke PB, Martin FE, Dell'Angelica EC, Cheney RE, et al. Myosin vc interacts with Rab32 and Rab38 proteins and works in the biogenesis and secretion of melanosomes. J Biol Chem. 2014;289:33513–28.
Wang C, Liu Z, Huang X. Rab32 is important for autophagy and lipid storage in Drosophila. PLoS One. 2012;7:e32086.
Freeman ME, Kanyicska B, Lerant A, Nagy G. Prolactin: structure, function, and regulation of secretion. Physiol Rev. 2000;80:1523–631.
Lacasse P, Lollivier V, Bruckmaier RM, Boisclair YR, Wagner GF, Boutinaud M. Effect of the prolactin-release inhibitor quinagolide on lactating dairy cows. J Dairy Sci. 2011;94:1302–9.
Boutinaud M, Lollivier V, Finot L, Bruckmaier RM, Lacasse P. Mammary cell activity and turnover in dairy cows treated with the prolactin-release inhibitor quinagolide and milked once daily. J Dairy Sci. 2012;95:177–87.
Ollier S, Zhao X, Lacasse P. Effects of feed restriction and prolactin-release inhibition at drying-off on susceptibility to new intramammary infection in cows. J Dairy Sci. 2015;98:221–8.
Thompson IR, Kaiser UB. GnRH pulse frequency-dependent differential regulation of LH and FSH gene expression. Mol Cell Endocrinol. 2014;385:28–35.
Schams D, Schallenberger E, Gombe S, Karg H. Endocrine patterns associated with puberty in male and female cattle. J Reprod Fertil Suppl. 1980;30:103–10.
Day ML. State of the art of GnRH-based timed AI in beef cattle. Anim Reprod. 2015;12:473–8.
Sanoudou D, Duka A, Drosatos K, Hayes KC, Zannis VI. Role of Esrrg in the fibrate-mediated regulation of lipid metabolism genes in human ApoA-I transgenic mice. Pharmacogenomics J. 2010;10:165–79.
Fortes MRS, Reverter A, Zhang Y, Collis E, Nagaraj SH, Jonsson NN, et al. Association weight matrix for the genetic dissection of puberty in beef cattle. Proc Natl Acad Sci U S A. 2010;107:13642–7.
Illumina, Infinium Genotyping Data Analysis. http://www.illumina.com/Documents/products/technotes/technote_infinium_genotyping_data_analysis.pdf. Accessed 1 Oct 2016.
O'Connell JR, Weeks DE. PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet. 1998;63:259–66.
Knaust J, Hadlich F, Weikard R, Kuehn C. Epistatic interactions between at least three loci determine the “rat-tail” phenotype in cattle. Genet Sel Evol. 2016;48:26.
Fadista J, Bendixen C. Genomic position mapping discrepancies of commercial SNP chips. PLoS One. 2012;7:e31025.
Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP. Summaries of Affymetrix GeneChip probe level data. Nucl Acids Res. 2003;31:e15.
Affymetrix, Quality Assessment of Exon and Gene Arrays. http://media.affymetrix.com/support/technical/whitepapers/exon_gene_arrays_qa_whitepaper.pdf. Accessed 1 Oct 2016.
Liu G, Loraine AE, Shigeta R, Cline M, Cheng J, Valmeekam V, et al. NetAffx: Affymetrix probesets and annotations. Nucl Acids Res. 2003;31:82–6.
Gene Expression Omnibus (GEO). http://www.ncbi.nlm.nih.gov/projects/geo/. Accessed 1 Oct 2016.
Gilmour AR, Gogel BJ, Cullis BR, Thompson R. ASReml User Guide Release 2.0. VSN International Ltd, Hemel Hempstead, HP1 1ES, UK. 2006.
Kruskal WH. A nonparametric test for the several sample problem. Ann Math Stat. 1952;23:525–40.
Mood AM, Graybill FA, Boes DC. Introduction to the Theory of Statistics. 3rd ed. Singapore: McGraw-Hill; 1974.
Büning H. Robust analysis of variance. J Appl Stat. 1997;24:319–32.
Gastwirth JL. Percentile modifications of two sample rank tests. J Am Stat Assoc. 1965;60:1127–41.
Hogg RV, Fisher DM, Randles RH. A two sample adaptive distribution-free test. J Am Stat Assoc. 1975;70:656–61.
Benjamini Y, Hochberg Y. Controlling the False Discovery Rate - A practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995;57:289–300.
We thank our colleagues at the FBN Dummerstorf, who helped in the generation and care of the SEGFAM F2 resource population for their continuous support of our work. Skillful technical assistance of A. Jugert, J. Bittner, N. Gentz and H. Tychsen is thankfully acknowledged.
This work is part of the project “PHENOMICS - a system biological approach of genotype-phenotype-distinction in the context of farm animal performance, health and welfare in cattle and pig” (0315536A (BB, JF, DZ), 0315536F (AZ) and 0315536G) (http://www.phaenomics.auf.uni-rostock.de/en-home.html) funded by the German Federal Ministry of Education and Research (BMBF). The funders had no role in study design, data collection, analysis and interpretation, preparation of the manuscript or decision to publish.
Availability of data and materials
All relevant data are within the paper and its supporting Information files. The gene expression data is available at Gene expression Omnibus (GEO, ) public repository (GEO accession number: GSE75371).
JF, DZ, BB and NR were responsible for generating and preprocessing the data. MOS and AZ performed the eQTL analysis. MS, SP, NR and AZ conceived and designed the study and analyses. BB and MOS drafted the manuscript. MS, SP, NR and AZ critically reviewed the manuscript. BB and MOS contributed equally to the work presented in this paper. All authors read and approved the manuscript.
The positions of BB, JF, DZ (0315536A) and MOS (0315536F) were financed by the German Federal Ministry of Education and Research (BMBF) for the conduct of the study with awards to MS (0315536A), AZ (0315536F) and NR (0315536G). The other authors declare that they have no competing interests.
Consent for publication
The experimental procedures were carried out according to the animal care guidelines of the State Mecklenburg-Vorpommern, Germany and were approved by the Landesamt für Landwirtschaft, Lebensmittelsicherheit und Fischerei Mecklenburg-Vorpommern (Reference number: LVL M-V/310-4/7221.3-2.1-017/03).
The original version of this article was revised to fix an error in Equation 8.
An erratum to this article is available at http://dx.doi.org/10.1186/s12863-016-0450-x.
List of eQTL. A complete list of eQTL comprising the chromosomal position and annotation of the SNPs and transcripts as well as the results of the adaptive linear rank test, the Kruskal Wallis-test and ANOVA are provided including the p-values for all tests and the corresponding Benjamini-Hochberg adjusted p-values for the adaptive linear rank test. (XLSX 75 kb)
Linkage disequilibrium decay plot. Genome-wide linkage disequilibrium (LD) decay plot for 145 cows of a F2 resource population deriving from a cross between Charolais and German Holstein founder breeds based on 37,204 polymorphic single nucleotide polymorphism (SNP) markers. The pairwise linkage disequilibrium (r2) between all pairs of SNPs on each chromosome was calculated and plotted against the genetic distance of the two markers. (PDF 521 kb)