Genome-wide association and biological pathway analysis for milk-fat composition in Danish Holstein and Danish Jersey cattle.

BACKGROUND: The milk fat profile of the Danish Holstein (DH) and Danish Jersey (DJ) show clear differences. Identification of the genomic regions, genes and biological pathways underlying the milk fat biosynthesis will improve the understanding of the biology underlying bovine milk fat production and may provide new possibilities to change the milk fat composition by selective breeding. In this study a genome wide association scan (GWAS) in the DH and DJ was performed for a detailed milk fatty acid (FA) profile using the HD bovine SNP array and subsequently a biological pathway analysis based on the SNP data was performed. RESULTS: The GWAS identified in total 1,233 SNPs (FDR < 0.10) spread over 18 chromosomes for nine different FA traits for the DH breed and 1,122 SNPs (FDR < 0.10) spread over 26 chromosomes for 13 different FA traits were detected for the DJ breed. Of these significant SNPs, 108 SNP markers were significant in both DH and DJ (C14-index, BTA26; C16, BTA14; fat percentage (FP), BTA14). This was supported by an enrichment test. The QTL on BTA14 and BTA26 represented the known candidate genes DGAT and SCD. In addition we suggest ACSS3 to be a good candidate gene for the QTL on BTA5 for C10:0 and C15:0. In addition, genetic correlations between the FA traits within breed showed large similarity across breeds. Furthermore, the biological pathway analysis revealed that fat digestion and absorption (KEGG04975) plays a role for the traits FP, C14:1, C16 index and C16:1. CONCLUSION: There was a clear similarity between the underlying genetics of FA in the milk between DH and DJ. This was supported by the fact that there was substantial overlap between SNPs for FP, C14 index, C14:1, C16 index and C16:1. In addition genetic correlations between FA showed a similar pattern across DH and DJ. Furthermore the biological pathway analysis suggested that fat digestion and absorption KEGG04975 is important for the traits FP, C14:1, C16 index and C16:1.

Exploring causal networks of bovine milk fatty acids in a multivariate mixed model context.

BACKGROUND: Knowledge regarding causal relationships among traits is important to understand complex biological systems. Structural equation models (SEM) can be used to quantify the causal relations between traits, which allow prediction of outcomes to interventions applied to such a network. Such models are fitted conditionally on a causal structure among traits, represented by a directed acyclic graph and an Inductive Causation (IC) algorithm can be used to search for causal structures. The aim of this study was to explore the space of causal structures involving bovine milk fatty acids and to select a network supported by data as the structure of a SEM. RESULTS: The IC algorithm adapted to mixed models settings was applied to study 14 correlated bovine milk fatty acids, resulting in an undirected network. The undirected pathway from C4:0 to C12:0 resembled the de novo synthesis pathway of short and medium chain saturated fatty acids. By using prior knowledge, directions were assigned to that part of the network and the resulting structure was used to fit a SEM that led to structural coefficients ranging from 0.85 to 1.05. The deviance information criterion indicated that the SEM was more plausible than the multi-trait model. CONCLUSIONS: The IC algorithm output pointed towards causal relations between the studied traits. This changed the focus from marginal associations between traits to direct relationships, thus towards relationships that may result in changes when external interventions are applied. The causal structure can give more insight into underlying mechanisms and the SEM can predict conditional changes due to such interventions.

Efficiency of genomic selection using Bayesian multi-marker models for traits selected to reflect a wide range of heritabilities and frequencies of detected quantitative traits loci in mice.

BACKGROUND: Genomic selection uses dense single nucleotide polymorphisms (SNP) markers to predict breeding values, as compared to conventional evaluations which estimate polygenic effects based on phenotypic records and pedigree information. The objective of this study was to compare polygenic, genomic and combined polygenic-genomic models, including mixture models (labelled according to the percentage of genotyped SNP markers considered to have a substantial effect, ranging from 2.5% to 100%). The data consisted of phenotypes and SNP genotypes (10,946 SNPs) of 2,188 mice. Various growth, behavioural and physiological traits were selected for the analysis to reflect a wide range of heritabilities (0.10 to 0.74) and numbers of detected quantitative traits loci (QTL) (1 to 20) affecting those traits. The analysis included estimation of variance components and cross-validation within and between families. RESULTS: Genomic selection showed a high predictive ability (PA) in comparison to traditional polygenic selection, especially for traits of moderate heritability and when cross-validation was between families. This occurred although the proportion of genomic variance of traits using genomic models was 22 to 33% smaller than using polygenic models. Using a 2.5% mixture genomic model, the proportion of genomic variance was 79% smaller relative to the polygenic model. Although the proportion of variance explained by the markers was reduced further when a smaller number of SNPs was assumed to have a substantial effect on the trait, PA of genomic selection for most traits was little affected. These low mixture percentages resulted in improved estimates of single SNP effects. Genomic models implemented for traits with fewer QTLs showed even lower PA than the polygenic models. CONCLUSIONS: Genomic selection generally performed better than traditional polygenic selection, especially in the context of between family cross-validation. Reducing the number of markers considered to affect the trait did not significantly change PA for most traits, particularly in the case of within family cross-validation, but increased the number of markers found to be associated with QTLs. The underlying number of QTLs affecting the trait has an effect on PA, with a smaller number of QTLs resulting in lower PA using the genomic model compared to the polygenic model.

Population and systems genetics analyses of cortisol in pigs divergently selected for stress.

This study presents a systems genetic analysis on the physiology of cortisol in mice and pigs with an aim to show the potential of a comprehensive computational approach to quickly identify candidate genes and avoid a costly whole-genome quantitative trait locus (QTL) mapping. Population genetics analyses were performed on measurements of cortisol from a pig selection experiment. Expression QTL were mapped and gene networks were built using gene expressions for Crhr1 (corticotrophin-releasing hormone receptor) gene and single nucleotide polymorphisms from public mouse data. Results from mouse data were used to infer potential candidate regulatory genes involved in pig cortisol regulation, using a comparative or translational systems genetics approach. The pig data used were from a 10-yr divergent genetic selection experiment, providing data on 417 individuals. Population genetics analysis showed that cortisol is highly genetically determined with heritabilities of 0.40-0.70. Furthermore, a major gene with an additive effect of 86 ng/ml is segregating. Genetical-genomics investigations revealed two trans-acting eQTL for Crhr1 gene expression on chromosomes 2 and 13. Candidate gene search under trans-eQTL peaks yielded 63 genes for Crhr1 expression phenotypes. Functional links for Crhr1 genes with other genes/proteins in the gene network using mouse data were shown for the first 10 statistically significant genes involved. Results show translational or comparative systems genetics approaches reduce costs and time in large-scale genetics and "-omics" investigations. This is the first study to report a strong genetic basis for cortisol physiology using a systems approach.

Biochemical pathways analysis of microarray results: regulation of myogenesis in pigs.

BACKGROUND: Combining microarray results and biological pathway information will add insight into biological processes. Pathway information is widely available in databases through the internet. Mammalian muscle formation has been previously studied using microarray technology in pigs because these animals are an interesting animal model for muscle formation due to selection for increased muscle mass. Results indicated regulation of the expression of genes involved in proliferation and differentiation of myoblasts, and energy metabolism. The aim of the present study was to analyse microarrays studying myogenesis in pigs. It was necessary to develop methods to search biochemical pathways databases. RESULTS: PERL scripts were developed that used the names of the genes on the microarray to search databases. Synonyms of gene names were added to the list by searching the Gene Ontology database. The KEGG database was searched for pathway information using this updated gene list. The KEGG database returned 88 pathways. Most genes were found in a single pathway, but others were found in up to seven pathways. Combining the pathways and the microarray information 21 pathways showed sufficient information content for further analysis. These pathways were related to regulation of several steps in myogenesis and energy metabolism. Pathways regulating myoblast proliferation and muscle fibre formation were described. Furthermore, two networks of pathways describing the formation of the myoblast cytoskeleton and regulation of the energy metabolism during myogenesis were presented. CONCLUSION: Combining microarray results and pathways information available through the internet provide biological insight in how the process of porcine myogenesis is regulated.

Evidence of a major gene from Bayesian segregation analyses of liability to osteochondral diseases in pigs.

Bayesian segregation analyses were used to investigate the mode of inheritance of osteochondral lesions (osteochondrosis, OC) in pigs. Data consisted of 1163 animals with OC and their pedigrees included 2891 animals. Mixed-inheritance threshold models (MITM) and several variants of MITM, in conjunction with Markov chain Monte Carlo methods, were developed for the analysis of these (categorical) data. Results showed major genes with significant and substantially higher variances (range 1.384-37.81), compared to the polygenic variance (sigmau2). Consequently, heritabilities for a mixed inheritance (range 0.65-0.90) were much higher than the heritabilities from the polygenes. Disease allele frequencies range was 0.38-0.88. Additional analyses estimating the transmission probabilities of the major gene showed clear evidence for Mendelian segregation of a major gene affecting osteochondrosis. The variants, MITM with informative prior on sigmau2, showed significant improvement in marginal distributions and accuracy of parameters. MITM with a "reduced polygenic model" for parameterization of polygenic effects avoided convergence problems and poor mixing encountered in an "individual polygenic model." In all cases, "shrinkage estimators" for fixed effects avoided unidentifiability for these parameters. The mixed-inheritance linear model (MILM) was also applied to all OC lesions and compared with the MITM. This is the first study to report evidence of major genes for osteochondral lesions in pigs; these results may also form a basis for underpinning the genetic inheritance of this disease in other animals as well as in humans.

Investigation of major gene for milk yield, milking speed, dry matter intake, and body weight in dairy cattle.

The main aim of this study was to determine if there exist any major gene for milk yield (MY), milking speed (MS), dry matter intake (DMI), and body weight (BW) recorded at various stages of lactation in first-lactation dairy cows (2543 observations from 320 cows) kept at the research farm of the Swiss Federal Institute of Technology between April 1994 and April 2004. Data were modelled based a simple repeatability covariance structure and analysed by using Bayesian segregation analyses. Gibbs sampling was used to make statistical inferences on posterior distributions; inferences were based on a single run of the Markov chain for each trait with 500,000 samples, with each 10th sample collected because of the high correlation among the samples. The posterior mean (+/-SD) of major gene variance was 2.61 (+/-2.46) for MY, 0.83 (+/-1.26) for MS, 4.37 (+/-2.34) for DMI, and 2056.43 (+/-665.67) for BW. Highest posterior density regions for 3 of the 4 traits did not include 0 (except MS), which supported the evidence for major gene. With additional tests for agreement with Mendelian transmission probabilities, we could only confirm the existence of a major gene for MY, but not for MS, DMI, and BW. Expected Mendelian transmission probabilities and their model fits were also compared.

The use of a genetic-counselling program by Dutch breeders for four hereditary health problems in boxer dogs.

Our group developed a genetic-counselling program for boxer-dog breeders in The Netherlands, using data for cryptorchidism (uni- and/or bilateral), epilepsy, knee-problems (including ligament rupture, fractured or ruptured meniscus, severe osteo-arthrosis of the knee, or a combination of these disorders), and schisis (including cheiloschisis, palatoschisis, or cheilopalatoschisis). We transformed the estimated breeding values (EBVs) into odds ratios (ORs), to enable the breeder to compare the risk for each of the traits for a certain dam-sire combination with the average population risk (set at 1). The goal of the study was to evaluate the use of our genetic-counselling program by Dutch breeders of boxer dogs. We asked breeders of the Dutch Boxer Club to send in an application form for genetic-counselling from June 1 to December 1, 2000. Breeders indicated on this application form three desirable sires for their dam (sire 1, sire 2, sire 3) in random order. On the basis of this information, a counselling report was produced which included ORs for the four diseases in litters of the dam-sire combinations indicated on the application form. Together with the counselling report, the breeders received an evaluation form. We received 129 application forms from 70 breeders, and collected 125 evaluations. Of these evaluations, 96 were informative about the influence of the counselling report on sire choice. The most-important criteria used by breeders to select sires were: the exterior characteristics (60%) and known progeny (52%). Although it was the first time breeders could make use of genetic-counselling, 32% of the breeders indicated that the genetic-counselling played a major role in their sire selection. Breeders expressed little difference in importance for the four genetic traits, but there was a tendency to consider epilepsy more than the others. Breeders hesitated to put long-term population interest above short-term personal interest. Nevertheless, the general conclusion of this study is that breeders acknowledge the value of counselling, especially to reduce the frequency of a threatening disease such as epilepsy.

Bayesian multi-QTL mapping for growth curve parameters.

BACKGROUND: Identification of QTL affecting a phenotype which is measured multiple times on the same experimental unit is not a trivial task because the repeated measures are not independent and in most cases show a trend in time. A complicating factor is that in most cases the mean increases non-linear with time as well as the variance. A two- step approach was used to analyze a simulated data set containing 1000 individuals with 5 measurements each. First the measurements were summarized in latent variables and subsequently a genome wide analysis was performed of these latent variables to identify segregating QTL using a Bayesian algorithm. RESULTS: For each individual a logistic growth curve was fitted and three latent variables: asymptote (ASYM), inflection point (XMID) and scaling factor (SCAL) were estimated per individual. Applying an 'animal' model showed heritabilities of approximately 48% for ASYM and SCAL while the heritability for XMID was approximately 24%. The genome wide scan revealed four QTLs affecting ASYM, one QTL affecting XMID and four QTLs affecting SCAL. The size of the QTL differed. QTL with a larger effect could be more precisely located compared to QTL with small effect. The locations of the QTLs for separate parameters were very close in some cases and probably caused the genetic correlation observed between ASYM and XMID and SCAL respectively. None of the QTL appeared on chromosome five. CONCLUSIONS: Repeated observations on individuals were affected by at least nine QTLs. For most QTL a precise location could be determined. The QTL for the inflection point (XMID) was difficult to pinpoint and might actually exist of two closely linked QTL on chromosome one.

From genetical genomics to systems genetics: potential applications in quantitative genomics and animal breeding.

This article reviews methods of integration of transcriptomics (and equally proteomics and metabolomics), genetics, and genomics in the form of systems genetics into existing genome analyses and their potential use in animal breeding and quantitative genomic modeling of complex traits. Genetical genomics or the expression quantitative trait loci (eQTL) mapping method and key findings in this research are reviewed. Various procedures and potential uses of eQTL mapping, global linkage clustering, and systems genetics are illustrated using actual analysis on recombinant inbred lines of mice with data on gene expression (for diabetes- and obesity-related genes), pathway, and single nucleotide polymorphism (SNP) linkage maps. Experimental and bioinformatics difficulties and possible solutions are discussed. The main uses of this systems genetics approach in quantitative genomics were shown to be in refinement of the identified QTL, candidate gene and SNP discovery, understanding gene-environment and gene-gene interactions, detection of candidate regulator genes/eQTL, discriminating multiple QTL/eQTL, and detection of pleiotropic QTL/eQTL, in addition to its use in reconstructing regulatory networks. The potential uses in animal breeding are direct selection on heritable gene expression measures, termed "expression assisted selection," and genetical genomic selection of both QTL and eQTL based on breeding values of the respective genes, termed "expression-assisted evaluation."

Gene-environment interactions in complex diseases: genetic models and methods for QTL mapping in multiple half-sib populations.

An interval quantitative trait locus (QTL) mapping method for complex polygenic diseases (as binary traits) showing QTL by environment interactions (QEI) was developed for outbred populations on a within-family basis. The main objectives, within the above context, were to investigate selection of genetic models and to compare liability or generalized interval mapping (GIM) and linear regression interval mapping (RIM) methods. Two different genetic models were used: one with main QTL and QEI effects (QEI model) and the other with only a main QTL effect (QTL model). Over 30 types of binary disease data as well as six types of continuous data were simulated and analysed by RIM and GIM. Using table values for significance testing, results show that RIM had an increased false detection rate (FDR) for testing interactions which was attributable to scale effects on the binary scale. GIM did not suffer from a high FDR for testing interactions. The use of empirical thresholds, which effectively means higher thresholds for RIM for testing interactions, could repair this increased FDR for RIM, but such empirical thresholds would have to be derived for each case because the amount of FDR depends on the incidence on the binary scale. RIM still suffered from higher biases (15-100% over- or under-estimation of true values) and high standard errors in QTL variance and location estimates than GIM for QEI models. Hence GIM is recommended for disease QTL mapping with QEI. In the presence of QEI, the model including QEI has more power (20-80% increase) to detect the QTL when the average QTL effect is small (in a situation where the model with a main QTL only is not too powerful). Top-down model selection is proposed in which a full test for QEI is conducted first and then the model is subsequently simplified. Methods and results will be applicable to human, plant and animal QTL mapping experiments.

Efficiency of population structures for mapping of Mendelian and imprinted quantitative trait loci in outbred pigs using variance component methods.

In a simulation study different designs for a pure line pig population were compared for efficiency of mapping QTL using the variance component method. Phenotypes affected by a Mendelian QTL, a paternally expressed QTL, a maternally expressed QTL or by a QTL without an effect were simulated. In all alternative designs 960 progeny were phenotyped. Given the limited number of animals there is an optimum between the number of families and the family size. Estimation of Mendelian and parentally expressed QTL is more efficient in a design with large family sizes. Too small a number of sires should be avoided to minimize chances of sires to be non-segregating. When a large number of families is used, the number of haplotypes increases which reduces the accuracy of estimating the QTL effect and thereby reduces the power to show a significant QTL and to correctly position the QTL. Dense maps allow for smaller family size due to exploitation of LD-information. Given the different possible modes of inheritance of the QTL using 8 to 16 boars, two litters per dam was optimal with respect to determining significance and correct location of the QTL for a data set consisting of 960 progeny. The variance component method combining linkage disequilibrium and linkage analysis seems to be an appropriate choice to analyze data sets which vary in marker density and which contain complex family structures.

Selection against genetic defects in conservation schemes while controlling inbreeding.

We studied different genetic models and evaluation systems to select against a genetic disease with additive, recessive or polygenic inheritance in genetic conservation schemes. When using optimum contribution selection with a restriction on the rate of inbreeding (DeltaF) to select against a disease allele, selection directly on DNA-genotypes is, as expected, the most efficient strategy. Selection for BLUP or segregation analysis breeding value estimates both need 1-2 generations more to halve the frequency of the disease allele, while these methods do not require knowledge of the disease mutation at the DNA level. BLUP and segregation analysis methods were equally efficient when selecting against a disease with single gene or complex polygene inheritance, i.e. knowledge about the mode of inheritance of the disease was not needed for efficient selection against the disease. Smaller schemes or schemes with a more stringent restriction on DeltaF needed more generations to halve the frequency of the disease alleles or the fraction of diseased animals. Optimum contribution selection maintained DeltaF at its predefined level, even when selection of females was at random. It is argued that in the investigated small conservation schemes with selection against a genetic defect, control of DeltaF is very important.