# Required properties for markers used to calculate unbiased estimates of the genetic correlation between populations

Yvonne C. J. Wientjes, Mario P. L. Calus, Pascal Duenk, Piter Bijma

► To cite this version:

Yvonne C. J. Wientjes, Mario P. L. Calus, Pascal Duenk, Piter Bijma. Required properties for markers used to calculate unbiased estimates of the genetic correlation between populations. Genetics Selection Evolution, 2018, 50 (1), pp.65. 10.1186/s12711-018-0434-6. http-02405290

HAL Id: hal-02405290

https://hal.science/hal-02405290

Submitted on 11 Dec 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Required properties for markers used to calculate unbiased estimates of the genetic correlation between populations

Yvonne C. J. Wientjes*, Mario P. L. Calus, Pascal Duenk and Piter Bijma

**Abstract**

**Background:** Generally, populations differ in terms of environmental and genetic factors, which can create differences in allele substitution effects between populations. Therefore, a single genotype may have different additive genetic values in different populations. The correlation between the two additive genetic values of a single genotype in two populations is known as the additive genetic correlation between populations and thus, can differ from 1. Our objective was to investigate whether differences in linkage disequilibrium (LD) and allele frequencies of markers and causal loci between populations affect the bias of the estimated genetic correlation. We simulated two populations that were separated by 50 generations and differed in LD pattern between markers and causal loci, as measured by the LD-statistic *r*. We used a high marker density to represent a high consistency of LD between populations, and lower marker densities to represent situations with a lower consistency of LD between populations. Markers and causal loci were selected to have either similar or different allele frequencies in the two populations.

**Results:** Our results show that genetic correlations were underestimated only slightly when the difference in allele frequencies between the two populations was similar for the markers and the causal loci. A lower marker density, representing a lower consistency of LD between populations, had only a minor effect on the underestimation of the genetic correlation. When the difference in allele frequencies between the two populations was not similar for markers and causal loci, genetic correlations were severely underestimated. This bias occurred because the markers did not predict accurately the relationships at causal loci.

**Conclusions:** For an unbiased estimation of the genetic correlation between populations, the markers should accurately predict the relationships at the causal loci. To achieve this, it is essential that the difference in allele frequencies between populations is similar for markers and causal loci. Our results show that differences in LD phase between causal loci and markers across populations have little effect on the estimated genetic correlation.

**Background**

Alleles in different populations are often expressed in different environments and genetic backgrounds. Because of genotype by environment interactions and non-additive genetic effects, these differences result in different allele substitution effects between populations [1–3]. In addition, the set of loci that underlie a trait can differ between populations. Therefore, a single genotype may have different additive genetic values in different populations [2,

4]. For each population, the additive genetic value is the product of the genotype, which is measured as the allele count at each locus, multiplied by the allele substitution effects for that population. The additive genetic correlation between two populations is the correlation between the two additive genetic values of a single genotype in the two populations and may differ considerably from 1.

Knowledge of the genetic correlation between populations helps to understand the differences and similarities in genetic architecture of complex traits between populations [5, 6]. For both genomic prediction and genome-wide association studies, combining information from populations is an attractive approach to increase the

<footnote>*Correspondence: yvonne.wientjes@wur.nl
Animal Breeding and Genomics, Wageningen University and Research,
6700 AH Wageningen, The Netherlands</footnote>
accuracy of estimated genetic values or the power to identify quantitative trait loci. This is especially the case when the number of individuals with marker genotypes and phenotypes in a population is limited. For both genomic prediction and genome-wide association studies, the genetic correlation between populations determines the added benefit of combining information from multiple populations [5, 7, 8]. Therefore, the genetic correlation between populations is an important parameter in human studies e.g. [5, 9], as well as in animal and plant breeding e.g. [10, 11].

To estimate a genetic correlation between two populations, it is essential to know the relationships between individuals from the two populations. Traditionally, relationships between individuals are based on pedigree information, which generally is only available within a population. The current availability of genome-wide marker panels has opened up new opportunities to estimate genetic correlations between populations of distantly related individuals, such as between breeds e.g. [10, 12], lines [13], sub-populations e.g. [11], or ethnicities e.g. [5, 9]. Genetic correlations between populations can be estimated using methods based on genomic relationships [10], random regression on marker genotypes [14, 15], or summary statistics of genome-wide association studies [6, 16]. Wientjes et al. [17] showed that it is possible to obtain an unbiased estimate of the genetic correlation from genomic relationships based on causal loci.

Because causal loci are generally unknown, genomic relationships have to be based on marker information. It is expected that this creates bias in the estimates of genetic correlation, because the strength and phase of linkage disequilibrium (LD) between causal loci and markers differ between populations in humans [18], livestock [19, 20], and plants [21, 22]. Due to imperfect LD between causal loci and markers, not all the genetic variance is explained by the markers, which can further distort the estimation of genetic correlations [16, 23]. In contrast to the expectations, a simulation study with populations that had different LD patterns showed that the estimated genetic correlation between populations based on marker information was only slightly biased [7]. This contrast indicates that the impact of differences in LD patterns between populations on the estimated genetic correlation remains unclear.

The objective of this study was to investigate whether differences in LD patterns between populations and differences in allele frequencies of markers and/or causal loci between populations affect bias of the estimated genetic correlation. We simulated two populations that were separated by 50 generations using scenarios that differed in consistency of LD between populations, defined

as the correlation in LD phase between the two populations, and in allele frequencies of markers and causal loci between the populations. We used marker-based relationship matrices to estimate the genetic correlation.

## Methods

### Population structure

Two populations were simulated using the QMSim software [24]. The simulations were set-up with two characteristics: (1) the two populations have different LD patterns, as measured by the LD statistic *r*; and (2) a large number of loci segregate in the last generation with part of these (> 200,000) having similar allele frequencies in both populations, and another part (> 200,000) having different allele frequencies in each population. We simulated a historical population for 211 generations (Fig. 1). The first generation (generation − 211) contained 300 individuals. During the following 100 generations (generations − 211 to − 112), population size decreased gradually to 50 individuals to create LD. From generation -111 to generation − 12, population size increased gradually to 300 individuals and was kept constant for the next 10 generations (generations − 11 to − 2). In the last generation of the historical population (generation − 1), population size increased to 1800 individuals. In the entire historical population, generations were discrete, there was no selection, mating was at random, and the male to female ratio was 1:5. The effective population size (*N*<sub>e</sub>) of the entire historical population was ~79, which was approximated as the harmonic mean of *N*<sub>e</sub> in each generation calculated as: $\frac{1}{N_e} = \frac{1}{4N_m} + \frac{1}{4N_f}$, where *N*<sub>m</sub> is the number of males and *N*<sub>f</sub> the number of females used to create the next generation [25].

The last generation of the historical population was divided randomly into two equally-sized populations (A and B) of 900 individuals. In the next generation, the size of both populations was increased to 1800 individuals and was kept constant for the following 40 generations (generations 1–40). These reasonably large population sizes limited the drift of allele frequencies. The number of offspring was set to 10 and selection was at random, such that the number of selected offspring per individual followed approximately a Poisson distribution, as assumed in the Wright-Fisher model of genetic drift. During the last 10 generations (generations 41–50), population size decreased to 120 individuals in each population to increase the extent of LD in each population, and the number of offspring was set to 20. Similar to the historical population, these generations were discrete, there was no selection, mating was at random, and the male to female ratio was 1:5. All individuals from the last generation (2000) were used for the analyses.
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Schematic representation of the simulated population structure</figcaption></figure>

## Genome size

A genome of 10 chromosomes of one Morgan each was simulated. This genome size was a balance between the computational effort necessary for the analyses and the

variation in relationships between family members [26]. Each chromosome contained 300,000 randomly spaced loci, with a recurrent mutation rate of 0.00005 in the historical population. In the last generation of the historical
population, segregating loci were selected and mutation was stopped. The chosen population size and mutation rate resulted in a U-shaped allele frequency distribution of segregating loci in the two populations, as common in real populations.

In the last generation (generation 50), markers and 2000 causal loci were selected from all segregating loci. Three marker panels were constructed: a high-density panel (HDP) with 200,000 markers, a medium-density panel (MDP) with 20,000 markers, and a low-density panel (LDP) with 2000 markers. Each of the smaller marker panels was a subset from the larger marker panels.

Markers and causal loci were selected to have either similar or different allele frequencies in populations A and B. For both approaches, three selection criteria were used; namely (1) the segregation in one or both populations, (2) the absolute difference in allele frequency between population A ($p_A$) and population B ($p_B$), and (3) the difference in $2p(1-p)$ between populations A and B, which is a measure of the difference in variance explained by a locus when allele substitution effects are the same in the two populations. The last criterion was mainly effective for loci with a low allele frequency, since an apparently small difference in allele frequency can result in a relatively large difference in variance explained for those loci. This criterion was used to ensure that the proportion of genetic variance explained by a locus was more or less similar in the two populations when populations had similar allele frequencies.

As a first step, markers were selected from the segregating loci. To select markers with similar allele frequencies in the two populations, (1) loci had to segregate in both populations, (2) $|p_A - p_B|$ should be less than 0.14, and (3) $|[2p_A(1-p_A) - 2p_B(1-p_B)]| / [2\tilde{p}_{AB}(1 - \tilde{p}_{AB})]$ should be less than 2, where $\tilde{p}_{AB}$ is the average of $p_A$ and $p_B$. To select markers with different allele frequencies in the two populations, (1) loci had to segregate in at least one population, (2) $|p_A - p_B|$ should be more than 0.14, and (3) $|[2p_A(1-p_A) - 2p_B(1-p_B)]| / [2\tilde{p}_{AB}(1 - \tilde{p}_{AB})]$ should be more than 1. The cut-off values were chosen to either minimize or maximize the difference in allele frequencies between the populations, while ensuring that enough loci in each replicate met the criteria. Our aim was to select marker panels with a uniform distribution of allele frequencies to reflect commercially available marker chips [27–30]. For this step, the loci that met the criteria were divided into 50 bins based on average allele frequency over the two populations (i.e., allele frequencies of bin 1 ranged from 0 to 0.02, of bin 2 from 0.02 to 0.04, etc.), and from each bin an equal number of loci was randomly selected. When the number of loci was too

small in the two extreme bins (0.00–0.02, and 0.98–1.00), the bins were combined with the neighbouring bin.

As a second step, causal loci were selected from the segregating loci that were not selected as markers. To select causal loci, the same criteria and cut-off values were used as for markers, with one exception. In all scenarios, causal loci did not have to segregate in both populations, since some causal loci are known to be at least partly population-specific [31]. Causal loci were selected randomly from all simulated loci that met the criteria, and therefore the pattern of their allele frequencies followed an approximate U-shaped distribution as expected for causal loci [32, 33].

## LD patterns and consistency of LD between populations

The LD pattern and consistency of LD between the populations was investigated. Within each population and between all causal loci and markers less than 10 cM apart, the parameter $r$ was calculated [34]:

$$r = \frac{(f_{11}f_{22} - f_{12}f_{21})}{\sqrt{f_{11}f_{22}f_{12}f_{21}}}$$

where $f_{11}$ is the haplotype frequency with allele 1 at the first locus and allele 1 at the second locus, $f_{22}$, $f_{12}$ and $f_{21}$ are frequencies of the other possible haplotypes, $f_1$ and $f_2$ are the frequencies of allele 1 and allele 2 at the first locus, and $f_1$ and $f_2$ are the frequencies of allele 1 and allele 2 at the second locus. The LD pattern within each population was represented by the average $r^2$ for intervals of 0.1 cM between the markers. The consistency of LD between the two simulated populations A and B was calculated as the correlation between $r$ values of the two populations for intervals of 0.1 cM, following De Roos et al. [35]. The consistency of LD between populations is known to decrease when the genomic distance between markers and causal loci increases [35]. This indicates that the populations were expected to have a higher consistency of LD between markers and causal loci when using HDP markers than when using MDP or LDP markers. In this way, the analysis with the HDP markers represented a situation with the highest consistency of LD between populations, and the analysis with LDP markers represented a situation with the lowest consistency of LD between populations.

## Phenotypes

For each causal locus, allele substitution effects were sampled from a bivariate normal distribution, with a mean of 0, a standard deviation of 1, and a correlation between the populations of 1, 0.8, 0.6, 0.4, 0.2 or 0. For each individual, its allele counts for the causal loci (coded as 0, 1,
and 2) were multiplied by the corresponding allele substitution effects and the results were summed over loci to calculate the additive genetic value (AGV) of the individual. The AGV were scaled to a mean of 0 and a variance of 1 across all individuals. Since allele substitution effects were sampled independently from allele frequency, the correlation between AGV of populations A and B (i.e., genetic correlation) was similar to the correlation between allele substitution effects (i.e., 1, 0.8, 0.6, 0.4, 0.2 or 0). A normally distributed environmental effect was sampled for each individual to obtain a heritability of 0.3 in each population. Phenotypes of all 2000 individuals in generation 50 were computed by summing the AGV and the environmental effects.

The simulation of phenotypes was replicated 50 times for each scenario. For each replicate, markers and causal loci were selected for three scenarios; (1) similar allele frequencies between populations for both markers and causal loci, (2) similar allele frequencies between populations for markers and different allele frequencies between populations for causal loci, and (3) different allele frequencies for both markers and causal loci. Within each scenario, phenotypes were simulated for each of the six genetic correlations. Scripts and seeds to simulate the data are in Additional file 2.

## Estimation of the genetic correlation

The additive genetic correlation between populations was estimated using the following bivariate model:

$$ \begin{bmatrix} \mathbf{y}_A \\ \mathbf{y}_B \end{bmatrix} = \begin{bmatrix} \mathbf{x}_A & \mathbf{0} \\ \mathbf{0} & \mathbf{x}_B \end{bmatrix} \begin{bmatrix} \mu_A \\ \mu_B \end{bmatrix} + \begin{bmatrix} \mathbf{Z}_A & \mathbf{0} \\ \mathbf{0} & \mathbf{Z}_B \end{bmatrix} \begin{bmatrix} \mathbf{a}_A \\ \mathbf{a}_B \end{bmatrix} + \begin{bmatrix} \mathbf{e}_A \\ \mathbf{e}_B \end{bmatrix}, $$

where $\mathbf{y}_k$ is a vector of phenotypes for population k ($k = A, B$), $\mathbf{x}_k$ is an incidence vector relating phenotypes to the mean in population k ($\mu_k$), $\mathbf{Z}_k$ is an incidence matrix relating phenotypes to estimated additive genetic values ($\begin{bmatrix} \mathbf{a}_A \\ \mathbf{a}_B \end{bmatrix} \sim N(\begin{bmatrix} \mathbf{0} \\ \mathbf{0} \end{bmatrix}, \begin{bmatrix} \sigma_A^2 & \sigma_{AB} \\ \sigma_{AB} & \sigma_B^2 \end{bmatrix} \otimes \begin{bmatrix} \mathbf{G}_{AA} & \mathbf{G}_{AB} \\ \mathbf{G}_{BA} & \mathbf{G}_{BB} \end{bmatrix})$) with $\otimes$ representing the Kronecker product function, and $\mathbf{e}_k$ are vectors with independent residual effects. Genetic and residual variances were estimated using restricted maximum likelihood (REML). The first analyses were performed using the ASReml software [36]. For the scenarios analysed later, we switched to MTG2 [37] to reduce computation time. We verified that the estimated variance components were identical using both programs.

The genomic relationship matrix (G) between all individuals was calculated as [17]:

$$ \begin{aligned} \mathbf{G} &= \begin{bmatrix} \mathbf{G}_{AA} & \mathbf{G}_{AB} \\ \mathbf{G}_{BA} & \mathbf{G}_{BB} \end{bmatrix} \\ &= \begin{bmatrix} \frac{\mathbf{W}_A \mathbf{W}_A'}{\sum \sqrt{2p_{Ai}(1-p_{Ai})}} & \frac{\mathbf{W}_A \mathbf{W}_B'}{\sqrt{\sum \sqrt{2p_{Ai}(1-p_{Ai})} \sqrt{2p_{Bi}(1-p_{Bi})}}} \\ \frac{\mathbf{W}_B \mathbf{W}_A'}{\sqrt{\sum \sqrt{2p_{Ai}(1-p_{Ai})} \sqrt{2p_{Bi}(1-p_{Bi})}}} & \frac{\mathbf{W}_B \mathbf{W}_B'}{\sum \sqrt{2p_{Bi}(1-p_{Bi})}} \end{bmatrix}, \end{aligned} $$

where $\mathbf{W}_k$ is a matrix with centered allele counts of all individuals from population k, and $p_{ki}$ is the allele frequency for locus i in population k. Centered allele counts were calculated as $g_{ijk} - 2p_{ki}$, where $g_{ijk}$ is the allele count of locus i for individual j from population k, coded as 0, 1 or 2. This G defines the relationships as standardized covariances between the genetic values of individuals [17]. In all scenarios and for all replicates, we calculated G using allele counts of (1) causal loci, (2) HDP markers, (3) MDP markers, or (4) LDP markers in both populations.

The relationships at causal loci are the true relationships for that trait, which are approximated when using markers. Marker-based relationships are subject to sampling error, since markers are a subset of the genome and in imperfect LD with the causal loci. A way to account for this sampling error is by regressing G towards the pedigree relationship matrix (A) [32, 38, 39], which is expected to reduce bias of the estimated variance components [32]. To investigate the effect of this regression, G matrices based on the three marker panels were regressed towards A and used for the scenarios with a correlation of 0.8 or 0.4.

Before regressing G towards A, the inbreeding level of each within-population block in G was rescaled to the inbreeding level in A, following [39]:

$$ \mathbf{G}^* = (1 - \bar{\mathbf{F}}_k)\mathbf{G} + 2\bar{\mathbf{F}}_k\mathbf{J}, $$

where $\bar{\mathbf{F}}_k$ is the average inbreeding coefficient of all individuals of population k based on the pedigree, and J is a matrix of 1s. The rescaled $\mathbf{G}^*$ was regressed towards A following [32, 38]:

$$ \hat{\mathbf{G}} = \mathbf{A} + b(\mathbf{G}^* - \mathbf{A}) \quad \text{with} \quad b = \frac{\mathrm{Var}(\mathbf{G}^* - \mathbf{A})}{\mathrm{Var}(\mathbf{G}^* - \mathbf{A}) + \frac{1}{n}}, $$

where n is the number of markers. To set-up A, the pedigree of the last 10 generations was used, such that between-population A relationships were zero. The regression was done separately within each population per bin of pedigree relationships (< 0.10, 0.10–0.25, 0.25–0.50, > 0.5) and between populations, since regression coefficients are higher for stronger pedigree relationships
[40, 41]. For the diagonal elements, only the inbreeding coefficients in $G^*$ were regressed towards $A$ [32]. Regression coefficients were all close to 1 for higher marker density panels (>0.99 for HDP and >0.97 for MDP). For LDP markers, regression coefficients were lower, i.e. ~0.84 for between-population relationships, ~0.89, ~0.91, ~0.94 and ~0.96 for the four bins of within-population relationships, and ~0.93 for inbreeding coefficients.

# Results

## Characteristics of the simulations

The criteria for selecting markers and causal loci resulted in clear differences between the scenarios with similar and different allele frequencies in the two populations (Fig. 2) and [see Additional file 1]. As designed in the simulations, the distribution of allele frequencies was uniform for markers and U-shaped for causal loci as expected in livestock populations. Therefore, the percentage of causal loci with a minor allele frequency lower than 0.05 was higher (on average 33% in each population) than the percentage of markers with a minor allele frequency lower than 0.05 (on average only 15% in each population). The decay of LD was similar in both populations (Fig. 3), with a stronger decay as distances between the loci in the 0–2 cM interval increased. The correlation of LD phase between the populations, representing the consistency of LD phase between the populations, decreased rapidly at short distances (0–5 cM), and fluctuated around zero at distances longer than 5 cM. As designed in the simulations, the extent and consistency of LD between populations were comparable to those in chicken and pig populations [20, 42–44].

## Proportion of explained variance

The proportion of the phenotypic variance explained by the markers, known as the genomic heritability [45], was close to the simulated heritability for the scenarios with HDP and MDP markers and slightly lower than the simulated heritability for the scenarios with LDP markers (estimated: ~0.29; simulated: 0.3). This implies that genetic variances were estimated accurately regardless of the marker panel used.

## Estimated genetic correlation

With relationships based on causal loci, all estimated genetic correlations were unbiased, irrespective of whether causal loci had similar or different allele frequencies in the two populations (Fig. 4). This was expected based on our previous results [17].

With relationships based on markers, all estimated genetic correlations were slightly to severely biased. The bias was very small when the difference in allele frequencies between the two populations was similar for the markers and the causal loci. For example, when marker-based relationships were not regressed towards the pedigree relationships, genetic correlations were underestimated by only ~2.5% for HDP, ~3% for MDP, and ~11% for LDP markers (Fig. 4a, c). The bias was much larger when markers had similar allele frequencies in the two populations and causal loci had different allele frequencies, i.e. when the difference in allele frequencies between the two populations differed between markers and causal loci (Fig. 4b; ~28% for HDP, ~30% for MDP, and ~41% for LDP markers). It should be noted that the

<figure><img src="image_4.png" /><figcaption>**Fig. 2** Allele frequencies of markers in two populations using two approaches to select markers. For one random replicate, allele frequencies of markers from both populations are plotted against each other for markers that are selected to have **a** similar allele frequencies or **b** different allele frequencies in the two populations</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Extent of LD in two populations and correlation of LD phase between the populations as a function of distance. The average LD (r²) between causal loci and markers for both populations, and the correlation of LD-phase (correlation of r) between the populations, as a function of distance between causal loci and markers for one random replicate</figcaption></figure>

distribution of allele frequencies was always uniform for markers and always U-shaped for causal loci.

Across all scenarios, regressing **G** towards **A** had only a small effect on the estimated genetic correlation (Fig. 5). At a high marker density, regressing **G** towards **A** lowered the estimated genetic correlation. Therefore, the underestimation for HDP and MDP markers increased from ~4 to ~9% when the difference in allele frequencies between the populations was similar for markers and causal loci, and from ~28 to ~32% when this was not the case. In contrast, regressing **G** towards **A** resulted in higher estimated genetic correlations at low marker density. For LDP markers, the underestimation decreased from ~12 to ~8% when the difference in allele frequencies between the populations was similar for markers and causal loci and from ~41 to ~38% when this was not the case. Thus, regressing **G** towards **A** was only beneficial for the estimation of the genetic correlation between populations when the marker density was low.

In general, standard errors of the mean across replicates for the estimated genetic correlation were small for all scenarios (~0.02), and tended to be slightly larger for lower true genetic correlations. Moreover, standard errors were slightly larger when the difference in allele frequencies between populations was not similar for markers and causal loci (Fig. 4b vs a, c). Regression of **G** towards **A** had no effect on the standard error.

## Genomic relationships

Estimates of the genetic variance are biased when the regression of true relationships on marker-based relationships is not equal to 1 [38]. We investigated whether this could explain the underestimation of the genetic correlation by considering the genomic relationships at the causal loci as the true relationships for that trait. In Figs. 6 and 7, we plotted the relationships at the causal loci versus the unregressed relationships at the markers for one of the replicates. The regression coefficients for within-population genomic relationships were close to 1, and were only slightly lower when causal loci had different allele frequencies (Fig. 7) compared to similar allele frequencies (Fig. 6) in the two populations. This means that the within-population relationships at the markers can predict quite accurately the relationships at the causal loci.

Regression coefficients of between-population relationships deviated more from 1, especially at low marker density. When the difference in allele frequencies between the populations was similar for markers and causal loci, the regression coefficients were equal to ~0.8 for HDP and MDP and 0.67 for LDP markers (Fig. 6). This means that the relationships at the markers led to an over-prediction of the relationships at the causal loci. When the difference in allele frequencies between the populations was not similar for markers and causal loci, regression coefficients of between-population relationships
a Populations with similar allele frequencies of both markers and causal loci

<figure><img src="image_6.png" /></figure>

b Populations with similar allele frequencies of markers and different allele frequencies of causal loci

<figure><img src="image_7.png" /></figure>

c Populations with different allele frequencies of both markers and causal loci

<figure><img src="image_8.png" /></figure>

**Fig. 4** Estimated genetic correlations between populations without regressing the genomic relationship matrix. The average estimated genetic correlation (± standard error of the mean) at different simulated genetic correlations for the scenario in which a markers and causal loci have similar allele frequencies in the two populations, b markers have similar and causal loci different allele frequencies in the two populations, or c markers and causal loci have different allele frequencies in the two populations, when the genomic relationship matrix is either based on the genotypes of causal loci (2000), HDP (200,000), MDP (20,000), or LDP (2000) markers without regression towards the pedigree relationship matrix. Standard errors were calculated as the standard deviation over replicates divided by the square root of the number of replicates
<figure><img src="image_9.png" /><figcaption>**Fig. 5** Estimated genetic correlations between populations with regression of the genomic relationship matrix. The average estimated genetic correlation (± standard error of the mean) at a simulated genetic correlation of a 0.8 or b 0.4 for the three scenarios with HDP (200,000), MDP (20,000), or LDP (2000) markers and regression of G towards A. Standard errors were calculated as the standard deviation over replicates divided by the square root of the number of replicates</figcaption></figure>

were equal to ~0.30 (Fig. 7). Thus, the over-prediction of between-population relationships using markers was much larger when the difference in allele frequencies between the populations was not similar for markers and causal loci.

The correlation between the relationships at the causal loci and at the markers, i.e., the accuracy of the marker-based relationships, decreased when the density of the markers decreased (Figs. 6, 7). When the difference in allele frequencies between the populations was similar for markers and causal loci, the correlation for within-population relationships was ~0.91 for HDP and MDP, and ~0.88 for LDP markers. The correlation for between-population relationships was ~0.70 for HDP and MDP, and 0.60 for LDP markers. The correlation between relationships at causal loci

and at markers was much lower when the difference in allele frequencies between the populations was not similar for markers and causal loci (within-population relationships: ~0.66 for HDP and MDP, ~0.63 for LDP; between-population relationships: ~0.09 for HDP and MDP, ~0.08 for LDP).

## Discussion

Our objective was to investigate whether differences in LD patterns between populations and differences in allele frequencies of markers and/or causal loci between populations affect bias of the estimated genetic correlation. We simulated two populations that differed in LD pattern between markers and causal loci, as measured by the LD-statistic *r*. Our results show that when the difference in allele frequencies between the two populations is similar
<figure><img src="image_10.png" /></figure>

for markers and causal loci, estimated genetic correlations are only slightly underestimated using markers. When the difference in allele frequencies between the two populations was not similar for markers and causal loci, genetic correlations were severely underestimated. Differences in LD and allele frequencies of causal loci between populations had only a very slight effect on the precision of the estimated genetic correlation.

## Estimation of genetic correlations using marker-based relationships

Estimates of the genetic variance and heritability are known to be biased when the regression coefficient of the true relationships on the marker-based relationships is not equal to 1, i.e., when $E(G_{causal loci} | G_{markers}) \neq G_{markers}$ [32, 38, 46]. When this regression coefficient is less than 1, relationships at the markers show too much variation,
<figure><img src="image_11.png" /></figure>

which results in an underestimation of the genetic variance. Yang et al. [32] argued that a regression coefficient less than 1 can be due to two effects: (1) sampling error on the relationships because the number of markers is finite; and (2) a difference in the distribution of allele frequencies between causal loci and markers. In all our scenarios, the number of markers was finite and the

distribution of allele frequencies differed between causal loci and markers. However, within populations, the estimated genomic heritability [45] was close to the simulated trait heritability for all scenarios. This suggests that the number of markers used was sufficient to constrain the sampling error on within-population relationships to an acceptable level, and that our estimated genetic
variances were affected only slightly by the difference in the distribution of allele frequencies between causal loci and markers. Thus, the underestimation of the genetic correlation between populations is not a result of biased estimates of the genetic variance.

The relative sampling error due to the use of a finite number of markers is much larger for between-population relationships than for within-population relationships, because more markers are needed to accurately estimate the small between-population relationships [38]. Moreover, we showed that the accuracy of predicting the between-population relationships at the causal loci using markers depended on whether the difference in allele frequency of causal loci between populations was reflected by the markers. These two factors can result in an underestimation of the genetic covariance between populations, which can explain the slight underestimation of the genetic correlation in the scenarios in which the difference in allele frequencies between the two populations was similar for markers and causal loci and the larger underestimation in the scenarios in which this was not the case. The higher sampling error on between-population relationships can also explain the larger underestimation of the genetic correlation for the LDP markers than for the HDP and MDP markers. Thus, to estimate the genetic correlation between populations, it is important that the difference in allele frequencies between the populations is similar for markers and causal loci and that the number of markers is sufficiently large.

The additive genetic correlation between populations is defined as the correlation between the two additive genetic values of a single genotype, measured as the allele count at each causal locus, in the two populations. This correlation is equal to the correlation between allele substitution effects when allele substitution effects are independent of allele frequency. This independency was used in our simulated phenotypes and also to set-up the G matrix, as it is an implicit assumption in Method 1 of VanRaden [47]. When this assumption is not met, the genetic correlation is no longer equal to the correlation between allele substitution effects, and violation of this assumption in the set-up of G may result in biased estimates of the genetic correlation.

## Regression of the maker-based relationships

Regressing **G** towards **A** is a way of correcting the marker-based relationships for the sampling error due to a finite number of markers [39]. The regression was strongest for LDP markers and reduced the underestimation of the genetic correlation. These results agree with the findings that regressing **G** towards **A** is important when the number of markers is small [32] and supports our statement that relationships at LDP markers were affected by

sampling error. However, regressing **G** towards **A** slightly increased the underestimation of the genetic correlation with HDP and MDP markers. The reason for this is not clear. It might be that the regression of **G** towards **A** not only reduces the sampling error, but also amplifies the effect of the difference in the distribution of allele frequencies between causal loci and markers.

In our study, regressing **G** towards **A** was slightly detrimental for the estimation of the genetic correlation when using HDP (200,000) or MDP (20,000) markers, with regression coefficients being only slightly less than 1, and it was beneficial when using LDP (2000) markers with regression coefficients being considerably less than 1. The simulated genome represented about one-third of the genome of livestock species such as cattle and chicken [48, 49]. This suggests that regressing **G** could be detrimental when using a genome-wide total of 60,000 or more markers in livestock. Note that this number of markers will depend on the consistency of LD between populations. Between-population relationships are all closer to zero when the LD pattern is less consistent between populations [50]. Such weaker relationships generally require more markers to reduce their relative sampling error to an acceptable level [32]. Hence, we think that the regression coefficients may be a better indicator for deciding whether or not to regress **G**; when all regression coefficients are close to 1, e.g., higher than 0.95, it is probably better not to regress **G** towards **A** when estimating the genetic correlation between populations.

The coefficients to regress **G** towards **A** were approximated using the number of markers and the variation in $G_{\text{markers}} - A$, assuming that the sampling error was only a result of using a limited number of markers [38]. To investigate the impact of this approximation, we repeated the analysis by using $b = \frac{\text{Cov}(G_{\text{causal loci}} - A, G_{\text{markers}} - A)}{\text{Var}(G_{\text{markers}} - A)}$ [38] as the regression coefficient to regress **G** towards **A**. This regression requires that the causal loci are known, which generally is not the case. Here, we used this approach to investigate whether the bias in the estimated genetic correlation indeed occurs because the marker-based relationships do not accurately predict the relationships at the causal loci, i.e., $E(G_{\text{causal loci}} | G_{\text{markers}}) \neq G_{\text{markers}}$, and not because of differences in LD between populations. We calculated *b* separately for within- and between-population relationships, using 11 bins based on pedigree relationships within populations ($(-0.05, -0.05-0.10, 0.10-0.15, 0.15-0.20, 0.20-0.25, 0.25-0.30, 0.30-0.35, 0.35-0.40, 0.40-0.50, > 0.50$, self-relationships) and 3 bins based on genomic relationships between populations ($(-0.10, -0.10-0.10, > 0.10$), and used those *b* values to rescale the relationships. As shown in Fig. 8, this rescaling removed almost completely the bias in the estimated genetic correlations based on HDP and MDP markers,
<figure><img src="image_12.png" /><figcaption>**Fig. 8** Estimated genetic correlations between populations after rescaling the marker-based genomic relationship matrix. The average estimated genetic correlation (± standard error of the mean) at different simulated genetic correlations for the scenario in which markers and causal loci have similar allele frequencies in the two populations when the genomic relationship matrix is either based on the genotypes of HDP (200,000), MDP (20,000), or LDP (2000) markers, after rescaling the marker-based relationships using a regression coefficient based on the relationships at causal loci. Standard errors were calculated as the standard deviation over replicates divided by the square root of the number of replicates</figcaption></figure>

but overestimated the genetic correlation when using rescaled relationships based on LDP markers. This might result from the much larger sampling error for LDP markers compared to HDP and MDP markers, which could result in underestimated *b* values. Thus, there appears to be a lower boundary for the number of markers necessary to calculate between-population genomic relationships that can be corrected using regression. Altogether, our results confirm that for an unbiased estimate of the genetic correlation between populations, the regression coefficient of true relationships on marker-based relationships should be equal to 1.

Consistency of LD between populations

When calculating the marker-based relationships, the current generation within each population was used as the base population, since we used current population-specific allele frequencies. This means that between-population relationships were zero on average. When the LD is at least partly consistent between the populations, due to the existence of a recent or distant common ancestor, between-population relationships will show variation around zero [50]. This variation is essential to estimate the genetic correlation between populations, and genetic correlation estimates are more precise when it is larger [51].

We expected that a lower consistency of LD between populations would reduce the estimated genetic correlation between populations, because it reduces the correlation between (apparent) marker effects. Surprisingly, our results showed that estimated genetic correlations were similar with HDP and MDP markers, and only

slightly lower with LDP markers. This can be explained by the potential of marker-based relationships to accurately predict the relationships at the causal loci, which is essential to estimate without bias the genetic (co)variances and the genetic correlation between populations. A lower consistency of LD between populations results in a smaller variation in between-population relationships [38, 50], both at causal loci and at markers. Therefore, the regression coefficient of the relationships at the causal loci on the relationships at the markers may not be greatly affected (Figs. 6, 7; HDP and MDP markers). Hence, the consistency of LD between the populations seems to have little impact on the estimated genetic correlation between populations.

The consistency of LD between populations does affect the correlation between the relationships at the causal loci and at the markers (Figs. 6, 7), i.e., the accuracy of the marker-based relationships. For an unbiased estimate of the genetic correlation between populations, the regression of true relationships on marker-relationships should be equal to 1 and marker-based relationships do not necessarily have to be accurate. This contrasts with the estimation of genetic values, as in genomic prediction, for which relationships must be accurate and must show variation [38]. Thus, an unbiased estimate of the genetic correlation between populations does not guarantee that accurate genomic prediction across populations can be performed.

In our study, the analysis with the HDP markers represented a situation with the highest consistency of LD between populations, and the analysis with LDP markers represented a situation with the lowest consistency of
LD between populations. Therefore, the effects of marker
density and consistency of LD between populations were
confounded. The combined impact of marker density and
consistency of LD appears to be limited, because the bias
in the HDP and MDP scenarios was similar and only a
little stronger in the LDP scenario. The impact of marker
density can be reduced by regressing the genomic rela-
tionships towards the pedigree relationships. If the causal
loci are known and the regression coefficients are calcu-
lated using relationships at causal loci, we showed that
this regression completely removed the bias in estimated
genetic correlations based on HDP and MDP markers.
This suggests that the slight bias in the HDP and MDP
scenario was due to marker density, and that differences
in LD between the populations had almost no effect.

Simulated population versus livestock populations

In order to investigate whether the simulated genome
represented livestock genomes, we compared the dis-
tribution of allele frequencies and LD pattern with real
genomes, as suggested by Daetwyler et al. [52]. The sim-
ulated genome showed a comparable pattern of allele
frequencies between markers and causal loci, and a com-
parable extent and consistency of LD between popula-
tions as shown in chicken and pig populations [20, 28,
30, 42–44]. Thus, our results can be translated directly to
livestock populations if the same marker density is used.
This simulated LD was much higher than that generally
observed in human populations [53, 54]. Since marker
density, and thereby the average LD between causal loci
and nearest marker, had almost no effect on the esti-
mated genetic correlation, it is expected that the simu-
lated LD pattern will not affect the results.

We simulated causal loci that were spread randomly
across the genome, which is not always the case in real
populations. When causal loci are enriched in regions
with either high or low LD, (co)variance estimates can
be over- or underestimated [46, 55]. However, we expect
that the impact of the heterogeneity of LD will be smaller
on the estimated genetic correlation than on the herit-
ability, since differences in LD across the genome affect
both the genetic variance and covariance estimates. This
mechanism may also explain why estimates of the genetic
correlation between traits within a population are less
affected by incomplete LD between causal loci and mark-
ers than genetic variance estimates [56].

As explained above, the genetic correlation is equal
to the correlation between allele substitution effects
when allele substitution effects are independent of
allele frequency. Differences in allele substitution effects
between populations result from non-additive genetic

effects and from differences in allele frequencies, and/
or from genotype by environment interactions [1–3]. To
date, the magnitude of these additive, dominance and
epistatic effects is not well known. Therefore, we chose
to simulate directly the different allele substitution
effects from a bivariate normal distribution, instead of
simulating the underlying non-additive effects.

Contrary to our simulations, selection generally
occurs in livestock populations. Selection creates nega-
tive correlations between causal loci, known as the Bul-
mer effect [57]. However, the impact of the Bulmer
effect on the correlation between loci is very small
because the number of causal loci (n_causal loci) is large
for most breeding goal traits (the average correlation is
at maximum $\frac{-1}{n_{causal loci}-1}$). Moreover, in general, selection
acts on multiple traits, which further reduces the cor-
relation between the causal loci affecting a trait. There-
fore, the Bulmer effect will have only a small effect on
the correlation between loci in one population. Since
selection is within population, the Bulmer effect does
not cause covariances between loci in different popula-
tions. For these reasons, we do not expect that selec-
tion and the Bulmer effect would have a large impact of
on the results of our study.

Furthermore, the Bulmer effect is a transient phenome-
non that depends on the type and intensity of selection.
Hence, the additive genetic variance as affected by the
Bulmer effect does not represent a fundamental biological
property of a population, but it can result from selection
decisions that may fluctuate over time. The biologically
relevant quantity is the genic (co)variance, which is always
twice the Mendelian sampling (co)variance. The relevance
of the genic (co)variance follows from the decomposition
of the additive genetic value of an individual into the
Mendelian sampling deviations of its ancestors [58],
**A** = **c'**m, where **c** is a vector of contributions of ancestors
to the individual (including the individual itself, for which
$c_i = 1$), and **m** is a vector of Mendelian sampling devia-
tions of those ancestors. Hence, the variance of the addi-
tive genetic values equals $V_A = c'\text{var}(\mathbf{m})c$. In the absence
of selection, $\text{var}(\mathbf{m})$ is diagonal and
$V_A = \sum_i c_i^2 \sigma_m^2 = (1^2 + 2 \times \frac{1^2}{2} + 4 \times \frac{1^2}{4} + \dots) \sigma_m^2 = 2\sigma_m^2$, where
$\sigma_m^2$ is the (full) Mendelian sampling variance. Selection
causes $\text{var}(m_i)$ to deviate from $\sigma_m^2$ and creates covariances
between Mendelian sampling deviations of different
ancestors. These deviations are transient, and are eroded
by recombination when selection ceases, so that $V_A$
returns to $2\sigma_m^2$ [57]. This illustrates that the genic (co)vari-
ance is the biologically relevant quantity. We have investi-
gated the estimation of the genic correlation, by focusing
on a population in the absence of selection, so that genic
and genetic correlations are equal.
**Implications**

Generally, marker panels are designed such that the markers have intermediate allele frequencies across multiple populations [27–29]. Hence, markers tend to have a higher average minor allele frequency than causal loci [32, 33]. Moreover, the difference in allele frequencies of causal loci between populations is probably not accurately represented by markers. These factors likely result in underestimated genetic correlations between populations using real data, but the impact of each of the factors requires further research.

**Conclusions**

For an unbiased estimate of the genetic correlation between populations based on marker information, it is important that marker-based relationships accurately predict the relationships at causal loci, i.e., $E(G_{causal loci} | G_{markers}) = G_{markers}$. This is obtained when the difference in allele frequencies between the two populations is similar for the markers and the causal loci, and the number of markers is sufficiently large to constrain the sampling error on between-population relationships to an acceptable level. Our results show that differences in LD phase between causal loci and markers across populations have little effect on the estimated genetic correlation (Additional file 2).

**Additional files**

**Additional file 1.** Allele frequency distribution of markers and causal loci. **Figure S1:** Allele frequency distribution of HDP markers for the two populations. **Figure S2:** Allele frequency distribution of causal loci for the two populations.

**Additional file 2.** Programs and seeds to simulate data. This file contains the input file used for QMSim, the Fortran-programs to select markers and causal loci for the different scenarios, the Fortran-program to simulate phenotypes and the seeds for the different programs in each of the replicates.

**Authors' contributions**

YCWJ, MPLC, PD, and PB (all authors) participated in the design of the study. YCWJ performed the simulations and statistical analyses and wrote the first draft of the paper. YCWJ, MPLC and PB were involved in the interpretation of the results. All authors read and approved the final manuscript.

**Acknowledgements**

Not applicable.

**Competing interests**

The authors declare that they have no competing interests.

**Availability of data and materials**

All data generated or analysed during this study are included in Additional file 2. This file contains the input file used for QMSim, the Fortran-programs to select markers and causal loci for the different scenarios, the Fortran-program to simulate phenotypes and the seeds for the different programs in each of the replicates.

**Consent for publication**

Not applicable.

**Ethics approval and consent to participate**

Not applicable.

**Funding**

This study was financially supported by NWO-TTW and the Breed4Food partners Cobb Europe, CRV, Hendrix Genetics and Topigs Norsvin. The use of the HPC cluster has been made possible by CAT-AgroFood (Shared Research Facilities Wageningen UR).

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Received: 28 June 2018 Accepted: 28 November 2018

Published online: 14 December 2018

**References**

1. Falconer DS. The problem of environment and selection. Am Nat. 1952;86:293–8.

2. Goodnight CJ. Epistasis and the increase in additive genetic variance: implication for phase 1 of Wright's shifting-balance process. Evolution. 1995;49:502–11.

3. Wright S. Evolution and the genetics of populations: the theory of gene frequencies, vol. 2. Chicago: University of Chicago Press; 1969.

4. Wade MJ, Goodnight CJ. The theories of Fisher and Wright in the context of metapopulations: when nature does many small experiments. Evolution. 1998;52:1537–53.

5. de Candia TR, Lee SH, Yang J, Browning BL, Gejman PV, Levinson DF, et al. Additive genetic variation in schizophrenia risk is shared by populations of African and European descent. Am J Hum Genet. 2013;93:463–70.

6. Brown BC. Asian genetic epidemiology network type 2 diabetes consortium, Ye CJ, Price AL, Zaitlen N. Transethnic genetic-correlation estimates from summary statistics. Am J Hum Genet. 2016;99:76–88.

7. Wientjes YCJ, Veerkamp RF, Bijma P, Bovenhuis H, Schrooten C, Calus MPL. Empirical and deterministic accuracies of across-population genomic prediction. Genet Sel Evol. 2015;47:5.

8. Wientjes YCJ, Bijma P, Veerkamp RF, Calus MPL. An equation to predict the accuracy of genomic values by combining data from multiple traits, populations, or environments. Genetics. 2016;202:799–823.

9. Yang L, Neale BM, Liu L, Lee SH, Wray NR, Ji N, et al. Polygenic transmission and complex neuro developmental network for attention deficit hyperactivity disorder: genome-wide association study of both common and rare variants. Am J Med Genet B Neuropsychiatr Genet. 2013;162B:419–30.

10. Karoui S, Carabaño MJ, Díaz C, Legarra A. Joint genomic evaluation of French dairy cattle breeds using multiple-trait models. Genet Sel Evol. 2012;44:39.

11. Lehermeier C, Schön CC, de los Campos G. Assessment of genetic heterogeneity in structured plant populations using multivariate whole-genome regression models. Genetics. 2015;201:323–37.

12. Carillier C, Larroque H, Robert-Granié C. Comparison of joint versus purebred genomic evaluation in the French multi-breed dairy goat population. Genet Sel Evol. 2014;46:67.

13. Huang H, Windig JJ, Vereijken A, Calus MPL. Genomic prediction based on data from three layer lines using non-linear regression models. Genet Sel Evol. 2014;46:75.

14. Krag K, Poulsen NA, Larsen MK, Larsen LB, Janss LL, Buitenhuis B. Genetic parameters for milk fatty acids in Danish Holstein cattle based on SNP markers using a Bayesian approach. BMC Genet. 2013;14:79.

15. Sørensen LP, Janss L, Madsen P, Mark T, Lund MS. Estimation of (co) variances for genomic regions of flexible sizes: application to complex infectious udder diseases in dairy cattle. Genet Sel Evol. 2012;44:18.
16. Bulik-Sullivan B, Finucane HK, Anttila V, Gusev A, Day FR, Loh PR, et al. An atlas of genetic correlations across human diseases and traits. Nat Genet. 2015;47:1236–41.

17. Wientjes YCJ, Bijma P, Vandenplas J, Calus MPL. Multi-population genomic relationships for estimating current genetic variances within and genetic correlations between populations. Genetics. 2017;207:503–15.

18. Sawyer SL, Mukherjee N, Pakstis AJ, Feuk L, Kidd JR, Brookes AJ, et al. Linkage disequilibrium patterns vary substantially among populations. Eur J Hum Genet. 2005;13:677–86.

19. Heifetz EM, Fulton JE, O'Sullivan N, Zhao H, Dekkers JCM, Soller M. Extent and consistency across generations of linkage disequilibrium in commercial layer chicken breeding populations. Genetics. 2005;171:1173–81.

20. Veroneze R, Lopes PS, Guimarães SEF, Silva FF, Lopes MS, Harlizius B, et al. Linkage disequilibrium and haplotype block structure in six commercial pig lines. J Anim Sci. 2013;91:3493–501.

21. Flint-Garcia SA, Thornsberry JM, Buckler ES 4th. Structure of linkage disequilibrium in plants. Annu Rev Plant Biol. 2003;54:357–74.

22. Lehermeier C, Krämer N, Bauer E, Bauland C, Camisan C, Campo L, et al. Usefulness of multiparental populations of maize (Zea mays L.) for genome-based prediction. Genetics. 2014;198:3–16.

23. Gianola D, de los Campos G, Toro MA, Naya H, Schön CC, Sorensen D. Do molecular markers inform about pleiotropy? Genetics. 2015;201:23–9.

24. Sargolzaei M, Schenkel FS. QMSim: a large-scale genome simulator for livestock. Bioinformatics. 2009;25:680–1.

25. Falconer DS, Mackay TFC. Introduction to quantitative genetics. 4th ed. Harlow: Pearson Education Limited; 1996.

26. Hill WG. Variation in genetic identity within kinships. Heredity. 1993;71:652–3.

27. Matsuzaki H, Dong S, Loi H, Di X, Liu G, Hubbell E, et al. Genotyping over 100,000 SNPs on a pair of oligonucleotide arrays. Nat Methods. 2004;1:109–11.

28. Groenen MAM, Megens HJ, Zare Y, Warren WC, Hillier LW, Crooijmans RPMA, et al. The development and characterization of a 60 K SNP chip for chicken. BMC Genomics. 2011;12:274.

29. Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, Heaton MP, et al. Development and characterization of a high density SNP genotyping assay for cattle. PLoS One. 2009;4:e5350.

30. Ramos AM, Crooijmans RPMA, Affara NA, Amaral AJ, Archibald AL, Beever JE, et al. Design of a high density SNP genotyping assay in the pig using SNPs identified and characterized by next generation sequencing technology. PLoS One. 2009;4:e6524.

31. Kemper KE, Hayes BJ, Daetwyler HD, Goddard ME. How old are quantitative trait loci and how widely do they segregate? J Anim Breed Genet. 2015;132:121–34.

32. Yang J, Benyamin B, McEvoy BP, Gordon S, Henders AK, Nyholt DR, et al. Common SNPs explain a large proportion of the heritability for human height. Nat Genet. 2010;42:565–9.

33. Kemper KE, Goddard ME. Understanding and predicting complex traits: knowledge from cattle. Hum Mol Genet. 2012;21:R45–51.

34. Hill WG, Robertson A. Linkage disequilibrium in finite populations. Theor Appl Genet. 1968;38:226–31.

35. de Roos APW, Hayes BJ, Spelman RJ, Goddard ME. Linkage disequilibrium and persistence of phase in Holstein-Friesian, Jersey and Angus cattle. Genetics. 2008;179:1503–12.

36. Gilmour AR, Gogel BJ, Cullis BR, Welham SJ, Thompson R. ASReml user guide release 4.1. Hemel Hempstead: VSN International Ltd; 2015.

37. Lee SH, van der Werf HJH. MTG2: an efficient algorithm for multivariate linear mixed model analysis based on genomic information. Bioinformatics. 2016;32:1420–2.

38. Goddard ME, Hayes BJ, Meuwissen THE. Using the genomic relationship matrix to predict the accuracy of genomic selection. J Anim Breed Genet. 2011;128:409–21.

39. Powell JE, Visscher PM, Goddard ME. Reconciling the analysis of IBD and IBS in complex trait studies. Nat Rev Genet. 2010;11:800–5.

40. Wientjes YCJ, Veerkamp RF, Calus MPL. The effect of linkage disequilibrium and family relationships on the reliability of genomic prediction. Genetics. 2013;193:621–31.

41. Veerkamp RF, Mulder HA, Thompson R, Calus MPL. Genomic and pedigree-based genetic parameters for scarcely recorded traits when some animals are genotyped. J Dairy Sci. 2011;94:4189–97.

42. Badke YM, Bates RO, Ernst CW, Schwab C, Teibel JP. Estimation of linkage disequilibrium in four US pig breeds. BMC Genomics. 2012;13:24.

43. Veroneze R, Bastiaansen JW, Knol EF, Guimarães SE, Silva FF, Harlizius B, et al. Linkage disequilibrium patterns and persistence of phase in purebred and crossbred pig (Sus scrofa) populations. BMC Genet. 2014;15:126.

44. Andreescu C, Avendano S, Brown SR, Hassen A, Lamont SJ, Dekkers JCM. Linkage disequilibrium in related breeding lines of chickens. Genetics. 2007;177:2161–9.

45. de los Campos G, Sorensen D, Gianola D. Genomic heritability: What is it? PLoS Genet. 2015;11:e1005048.

46. Yang J, Bakshi A, Zhu Z, Hemani G, Vinkhuyzen AAE, Lee SH, et al. Genetic variance estimation with imputed variants finds negligible missing heritability for human height and body mass index. Nat Genet. 2015;47:1114–20.

47. VanRaden PM. Efficient methods to compute genomic predictions. J Dairy Sci. 2008;91:4414–23.

48. Ihara N, Takasuga A, Mizoshita K, Takeda H, Sugimoto M, Mizoguchi Y, et al. A comprehensive genetic map of the cattle genome based on 3802 microsatellites. Genome Res. 2004;14:1987–98.

49. Groenen MAM, Wahlberg P, Foglio M, Cheng HH, Megens HJ, Crooijmans RPMA, et al. A high-density SNP-based linkage map of the chicken genome reveals sequence features correlated with recombination rate. Genome Res. 2009;19:510–9.

50. Goddard ME. Genomic selection: prediction of accuracy and maximisation of long term response. Genetica. 2009;136:245–57.

51. Visscher PM, Hemani G, Vinkhuyzen AAE, Chen GB, Lee SH, Wray NR, et al. Statistical power to detect genetic (co)variance of complex traits using SNP data in unrelated samples. PLoS Genet. 2014;10:e1004269.

52. Daetwyler HD, Calus MPL, Pong-Wong R, de los Campos G, Hickey JM. Genomic prediction in animals and plants: simulation of data, validation, reporting, and benchmarking. Genetics. 2013;193:347–65.

53. Pritchard JK, Przeworski M. Linkage disequilibrium in humans: models and data. Am J Hum Genet. 2001;69:1–14.

54. Shifman S, Kuypers J, Kokoris M, Yakir B, Darvasi A. Linkage disequilibrium patterns of the human genome across populations. Hum Mol Genet. 2003;12:771–6.

55. Speed D, Hemani G, Johnson MR, Balding DJ. Improved heritability estimation from genome-wide SNPs. Am J Hum Genet. 2012;91:1011–21.

56. Trzaskowski M, Davis OSP, DeFries JC, Yang J, Visscher PM, Plomin R. DNA evidence for strong genome-wide pleiotropy of cognitive and learning abilities. Behav Genet. 2013;43:267–73.

57. Bulmer MG. The effect of selection on genetic variability. Am Nat. 1971;105:201–11.

58. Thompson R. The estimation of heritability with unbalanced data: ii. Data available on more than two generations. Biometrics. 1977;33:497–504.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_13.png" /></figure>