# A meta-analysis including pre-selected sequence variants associated with seven traits in three French dairy cattle populations

Andrew G. Marete, Bernt Guldbrandtsen, Mogens S. Lund, Sebastien Fritz, Goutam Sahana, Didier Boichard

► To cite this version:

Andrew G. Marete, Bernt Guldbrandtsen, Mogens S. Lund, Sebastien Fritz, Goutam Sahana, et al. A meta-analysis including pre-selected sequence variants associated with seven traits in three French dairy cattle populations. Frontiers in Genetics, 2018, 9, pp.522. 10.3389/fgene.2018.00522. hal-02625793

HAL Id: hal-02625793

https://hal.inrae.fr/hal-02625793

Submitted on 26 May 2020

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>

# A Meta-Analysis Including Pre-selected Sequence Variants Associated With Seven Traits in Three French Dairy Cattle Populations

Andrew G. Marete¹,²*, Bernt Guldbrandtsen², Mogens S. Lund², Sébastien Fritz¹,³, Goutam Sahana² and Didier Boichard¹

¹ UMR GABI, INRA, AgroParisTech, Université Paris Saclay, 78350 Jouy en Josas, France, ² Center for Quantitative Genetics and Genomics, Aarhus University, Aarhus, Denmark, ³ ALLICE, Paris, France

A within-breed genome-wide association study (GWAS) is useful when identifying the QTL that segregates in a breed. However, an across-breed meta-analysis can be used to increase the power of identification and precise localization of QTL that segregate in multiple breeds. Precise localization will allow including QTL information from other breeds in genomic prediction due to the persistence of the linkage phase between the causal variant and the marker. This study aimed to identify and confirm QTL detected in within-breed GWAS through a meta-analysis in three French dairy cattle breeds. A set of sequence variants selected based on their functional annotations were imputed into 50 k genotypes for 46,732 Holstein, 20,096 Montbeliarde, and 11,944 Normande cows to identify QTL for milk production, the success rate at insemination of cows (fertility) and stature. We conducted within-breed GWAS followed by across-breed meta-analysis using a weighted Z-scores model on the GWAS summary data (i.e., *P*-values, effect direction, and sample size). After Bonferroni correction, the GWAS result identified 21,956 significantly associated SNP (*P*<sub>FWER</sub> < 0.05), while meta-analysis result identified 9,604 significant SNP (*P*<sub>FWER</sub> < 0.05) associated with the phenotypes. The meta-analysis identified 36 QTL for milk yield, 48 QTL for fat yield and percentage, 29 QTL for protein yield and percentage, 13 QTL for fertility, and 16 QTL for stature. Some of these QTL were not significant in the within-breed GWAS. Some previously identified causal variants were confirmed, e.g., BTA14:1802265 (fat percentage, *P* = 1.5 × 10<sup>-760</sup>; protein percentage, *P* = 7.61 × 10<sup>-348</sup>) both mapping the DGAT1-K232A mutation and BTA14:25006125 (*P* = 8.58 × 10<sup>-140</sup>) mapping *PLAG1* gene was confirmed for stature in Montbeliarde. New QTL lead SNP shared between breeds included the intronic variant rs109205829 (*NFIB* gene), and the intergenic variant rs41592357 (1.38 Mb upstream of the *CNTN6* gene and 0.65 Mb downstream of the *CNTN4* gene). Rs110425867 (*ZFAT* gene) was the top variant associated with fertility, and new QTL lead SNP included rs109483390 (0.1 Mb upstream of the *TNFAIP3* gene and 0.07 Mb downstream of *PERP* gene), and rs42412333 (0.45 Mb downstream of the *RPL10L* gene). An across-breed meta-analysis had greater power to detect QTL as opposed to a within breed GWAS. The QTL detected here can be incorporated in routine genomic predictions.

**Keywords:** candidate genes, candidate SNP, GWAS, meta-analysis, sequence imputation

## OPEN ACCESS

**Edited by:**
Tad Stewart Sonstegard,
Recombinetics, United States

**Reviewed by:**
Dan Nonneman,
Agricultural Research Service (USDA),
United States
Luca Fontanesi,
Università degli Studi di Bologna, Italy

**\*Correspondence:**
Andrew G. Marete
didier.boichard@inra.fr

**Specialty section:**
This article was submitted to Livestock Genomics, a section of the journal Frontiers in Genetics

**Received:** 03 June 2018
**Accepted:** 16 October 2018
**Published:** 06 November 2018

**Citation:**
Marete AG, Guldbrandtsen B, Lund MS, Fritz S, Sahana G and Boichard D (2018) A Meta-Analysis Including Pre-selected Sequence Variants Associated With Seven Traits in Three French Dairy Cattle Populations. Front. Genet. 9:522. doi: 10.3389/fgene.2018.00522
# INTRODUCTION

Ideally, the magnitude of estimated effects on the quantitative trait of interest could be used to rank single nucleotide polymorphisms (SNP) for a functional genomic study to identify causal variants. The physical locations of the associated SNP on the genome are then flagged and mined for causative variants underlying quantitative trait loci (QTL). Causal variants reported for QTL in previous dairy cattle studies include polymorphisms causing variation in milk production, fertility, and embryo mortality (Grisart et al., 2002; Hoff et al., 2017; Michot et al., 2017; Bouwman et al., 2018). Most analyses testing genomic regions, like those above, perform within-breed GWAS (Chang et al., 2018). Even so, our knowledge regarding causal variants in these genomic regions remains limited because the GWAS analyses yield similar P-values for many adjacent SNP variants, a consequence of linkage disequilibrium (LD). Even with a multitude of GWAS results, strong LD prevents distinguishing the actual causal variant from linked markers (Goddard and Hayes, 2009). A meta-analysis can be used to improve the resolution of QTL detection and identify causal variants provided that LD is conserved at short distances across breeds (van den Berg et al., 2016). The main advantage of a meta-analysis is that it allows simultaneous analysis of many breeds by combining GWAS summary statistics across populations, thereby increasing power to detect QTL (van den Berg et al., 2016; Bouwman et al., 2018).

In this paper, we report an across-breed meta-analysis of GWAS summary statistics based on imputed pre-selected sequence variant genotypes from 78,772 cattle from three dairy breeds and for seven traits. The meta-analysis across breeds allowed us to identify 142 QTL for milk production, stature, and fertility.

# MATERIALS AND METHODS

## Studied Population and Phenotypes

We studied phenotypes from three French dairy cattle breeds. No ethical approval was necessary since the data used was from existing databases. Phenotypes were from cows born between 2007 and 2013 in the French dairy production regions and were defined as follows: (1) Production traits including milk, protein, and fat yields were obtained from test-day records expressed as 305 d yield (kg). Fat content (g/kg) was calculated as 1,000*(fat yield/milk yield). Protein content (g/kg) was calculated as 1,000* (protein yield/milk yield). (2) Fertility was defined as success/failure at each insemination of lactating cows and recorded as a success (=1) or a failure (=0). (3) Stature was measured as the vertical distance from the plank to the sacrum.

These traits were analyzed using the French national evaluation model (Boichard et al., 2012b) to obtain Yield Deviations (YD) (VanRaden and Wiggans, 1991). A YD can be defined as the mean performance adjusted for all environmental effects, including the permanent environment for cows with multiple records. The model varied according to the trait, especially concerning environmental factors. We retained traits

for genotyped cows and this included data from 46,732 Holstein cows, 20,096 Montbeliarde cows, and 11,944 Normande cows.

## Genotyping and Imputation

Cows were genotyped using the Illumina Infinium® BovineSNP50 BeadChip (50 K, Illumina, San Diego, CA), BovineLD BeadChip (Boichard et al., 2012a) or one of the first four versions of the EuroG10k SNP chip (Boichard et al., 2018). The EuroG10k SNP chip is composed of two parts: (1) Between 7,000 and 8,500 generic (and supposedly neutral) SNP from BovineLD Genotyping BeadChip v.2 (Boichard et al., 2012a). (2) Whole genome sequence variants-a custom part of up to 7,232 SNP selected from sequence data as part of 1,000 Bull Genomes Project Run 4 (Daetwyler et al., 2014) based on their functional annotation. The full description of the EuroG10k chip and quality control procedure was previously described by Marete et al. (2018a) and Boichard et al. (2018). Cows with phenotypes were rather old animals and were genotyped with the 50 k chip. Younger animals including cows without phenotypes were genotyped with the EuroG10k chip. To get complete marker information across animals, and to include the sequence variants in the analysis, we run two successive steps of imputation for the entire populations (males and females, with and without phenotypes, see Figure 1). First, we imputed all 50 k markers for all animals using animals with 50 k genotypes as a reference. Then we imputed sequence variants using the animals genotyped with the EuroG10k chip as the reference population. The imputation was done within breed using the FImpute software (Sargolzaei et al., 2014). After imputation, 48,576 SNP distributed across 29 Bos taurus autosomes (BTA) remained. Of these, 42,967 were from the 50 k chip, and 5,609 sequence variants selected from French cattle populations. Imputation accuracy was estimated using both concordance rate and allelic squared correlation ($r^2$), both for animals genotyped with 50 k only and those genotyped with EuroG10k chip only. In both instances, we randomly masked 20% of the genotyped markers on each chromosome and 15% of the cows per breed, then performed an imputation and compared the imputed genotypes with true genotypes by estimating a Pearson correlation coefficient. For the 50 k SNP, the average concordance was >0.98, and average allelic $r^2$ was >0.97, and for the sequence variants, the average concordance was >0.95, and average $r^2$ was >0.96 for all breeds.

## Association Studies, Meta-Analysis, and QTL Heterogeneity

First, 48,576 SNP (pre-selected sequence variants and 50 k SNP) were tested individually for association with each trait within the breed. The GCTA software (Yang et al., 2011) was used to fit a mixed linear regression model to test associations between a SNP and the trait. For any given trait, the fitted model was

$$y = 1\mu + x\beta + g + \varepsilon \quad (1)$$

where y was a vector of phenotypes (yield deviation) for all cows, 1 was a row vector of 1s, μ was the mean, x was the vector of true or imputed allele dosages for the SNP, β was the
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Genotype imputation workflow. Table showing number of animals genotyped with various chips is included in **Supplementary Table 1**.</figcaption></figure>

fixed allele substitution effect of the SNP, g was the vector of random additive genetic effects with g ~ N(0, Gσ²g), G being the genomic relationship matrix, σ²g was the variance explained by all the SNPs, ε was a vector of random residual effects with ε ~ N(0, Iσ²ε), σ²ε was the error variance, and I was an identity matrix. The variance of y was var(y) = Gσ²g + Iσ²ε. The gjk term of G matrix was estimated by

$$g_{jk} = \frac{1}{w} \sum_{i=1}^{w} \frac{(z_{ij} - 2p_i)(z_{ik} - 2p_i)}{2p_i(1 - p_i)} \quad (2)$$

where w was the total number of SNP, zij and zik were numbers of copies of the reference allele for the ith SNP in the jth and kth cows, respectively, and pi was the frequency of the reference allele estimated from the marker data (VanRaden, 2008).

To control the type I error rate and assuming all tests were independent (a conservative assumption due to LD among SNP), a Bonferroni genome-wide correction was applied for the total number of tests for each trait i.e., the nominal type I error rate (α = 5%) was divided by the number of SNP (w = 48,576) to obtain a genome-wide error rate threshold of t = 1.03 × 10⁻⁶. Any SNP whose probability of observing the test statistics was below this value was considered genome-wide significant.

The within-breed analyses were followed by a meta-analysis using weighted Z-score model as implemented in METAL software (Willer et al., 2010). The weighted Z-score model used P-values and directions of effect estimates and weights individual GWAS based on the sample size to compute a Z-score, i.e.,

$$Z_b = \Phi^{-1} \left( 1 - \frac{P_b}{2} \right) * \Delta b \quad (3)$$

where Zb was the Z-score for breed b, Φ was the standard normal cumulative distribution, Pb was the P-value, Δb was the direction of the SNP effect estimated within breed b. Overall Z-score was calculated as:

$$Z = \frac{\sum_b Z_b w_b}{\sqrt{\sum_b w_b^2}} \quad (4)$$

where w_b = √N_b and N_b was the sample size for breed b. The overall P-value was calculated as:

$$P = 2\Phi(-|Z|) \quad (5)$$

We evaluated heterogeneity of the effect sizes across breeds using Cochran's Q-test (Cochran, 1954) as implemented in the METAL software. A QTL was deemed to segregate across populations if the heterogeneity P-value was ≤ 0.05. The functional consequence of significantly associated variants was predicted using the Variant Effect Predictor tool from Ensembl Genome Browser 90 (McLaren et al., 2016). We classified variants that (i) mapped to a gene, (ii) <5-kb to known genes, and, (iii) >5-kb from any coding region. In case of multiple variants representing the same gene, we kept the variant associated with most traits and with the lowest P-value as representative for that gene.

## Identification of QTL That Segregates Within and Across Breeds

Since it is challenging to identify the causative variant, we aimed at determining the QTL region that most likely harbor the causal variant. First, for each chromosome, we subjectively inspected significantly associated variants within the genomic regions. Secondly, we defined the within-breed QTL region as 1 Mb windows and the SNP with the lowest P-value designated as the lead SNP. A QTL was selected as significant in more than one breed if any of the SNP in the QTL region of the other breed had a P-value lower than the genome-wide threshold (P < 1.03 × 10⁻⁶). The UMD3.1 (Zimin et al., 2009) bovine genome annotation was used to annotate genes centering within 1 Mb intervals on the lead SNP. We compared these annotations to known QTL for bovine milk production, fertility and stature traits using QTLdb (Hu et al., 2016) and literature review.
# RESULTS

## Association Studies

Overall, there were 21,956 SNP significantly associated with traits in three French cattle breeds. Forty-three percent of these associated SNP were from Holstein, 24% from Montbeliarde, 14% from Normande. Seven percent of these SNP showed association in both Holstein and Montbeliarde, 4% in Holstein and Normande, 5% in Montbeliarde and Holstein, and 3% in all three breeds. Production traits and stature had more significant hits across breeds, compared to fertility. As presented in **Figures 2A–E**, the overlap of significant tests in any two-breed combination and for all traits was most evident between Holstein and Montbeliarde (1,621 SNP) and least evident between Montbeliarde and Normande (1,082 SNP). Seven hundred and one SNP were significantly associated with the same trait in all three breeds.

As presented in **Figures 2F–I**, some of the imputed sequence variants had lower *P*-values compared with the 50 k variants *P*-values. Although more than 13% of the pre-selected sequence variants had MAF between 0.5 and 5%, only 2 QTL lead SNP had MAF lower than 5%. We did not detect QTL with MAF <5% in Normande. Among the pre-selected sequence variants, 1,478 were significant (*P* < 1.03 × 10⁻⁶).

## Within Breed GWAS And Multi-Breed Meta-Analysis

### Production Traits

#### Milk yield

We observed both known and novel QTL in all three cattle breeds. The most significant QTL lead SNP from meta-analysis are presented in **Table 1**. We observed QTL on 16 chromosomes where Holstein had 23 QTL, Montbeliarde had 10 QTL, Normande had 15 QTL. The most significant lead SNP was located at *Bos taurus autosome* (BTA) 14:1801116 (rs109421300, *P* = 9.7 × 10⁻²⁸³), an intronic variant near the causal variants in the *DGAT1* gene. Other significant QTL were observed at BTA5:50804085 (*P* = 1.3 × 10⁻²², *intergenic*), rs109205829 (BTA8:30088133, *P* = 2.3 × 10⁻²¹, *NFIB* gene), and BTA20:32670639 (*P* = 4.0 × 10⁻⁴¹, intergenic flanking 3' end of the *GHR* gene).

A QTL on BTA8 at 29–30 Mb with the same lead SNP at BTA8:30088133 (rs109205829) was observed in both Holstein and Montbeliarde. This lead SNP is an intronic variant within the *NFIB* gene, and the meta-analysis had a lower *P*-value (*P* = 2.3 × 10⁻²²) than either of the single breed GWAS (*P* = 8.4 × 10⁻¹⁴). Another QTL on BTA22 at 23–26 Mb with a lead SNP at BTA22:24219999 (rs41592357) was observed in all three breeds. Rs41592357 is an intergenic variant 1.38 Mb upstream of the *CNTN6* gene and 0.65 Mb downstream of the *CNTN4* gene.

<figure><img src="image_4.png" /></figure>
<table><caption>TABLE 1 | Lead SNP variants at 36 QTL for milk yield.</caption><thead><tr><th rowspan="2">QTL</th><th rowspan="2">Reference SNP<br>cluster ID</th><th rowspan="2">CHR</th><th rowspan="2">BP</th><th rowspan="2">Allele</th><th rowspan="2">Nearby<br>gene</th><th rowspan="2">P<sub>meta</sub></th><th rowspan="2">Effect direction (and evidence of<br>heterogeneity) across three<br>populations<sup>1</sup></th><th colspan="2">MON</th><th colspan="2">NOR</th><th colspan="2">HOL</th></tr><tr><th>b(se)</th><th>p</th><th>b(se)</th><th>p</th><th>b(se)</th><th>p</th></tr></thead><tbody><tr><td>1</td><td>rs29015637</td><td>2</td><td>9810493</td><td>a/g</td><td>—</td><td>4.2e-09</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>2</td><td>rs41615145</td><td>2</td><td>45026036</td><td>t/c</td><td>—</td><td>4.5e-17</td><td>++</td><td>++(1.5e-02)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>3</td><td>rs42829305</td><td>2</td><td>70680985</td><td>t/g</td><td>—</td><td>1.8e-12</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>4</td><td>rs43709763</td><td>3</td><td>20031136</td><td>a/c</td><td>CTSS</td><td>3.0e-15</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>5</td><td>rs29023272</td><td>3</td><td>50804085</td><td>a/g</td><td>RPL5</td><td>1.3e-22</td><td>++</td><td>++(5.9e-04)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>6</td><td>rs42737927</td><td>3</td><td>113967407</td><td>t/c</td><td>UGT1A1</td><td>1.3e-11</td><td>++</td><td>++(6.5e-03)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>7</td><td>rs109960400</td><td>7</td><td>19676040</td><td>t/c</td><td>RANBP3</td><td>4.8e-08</td><td>++</td><td>++(8.2e-02)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>8</td><td>rs110220494</td><td>7</td><td>46109256</td><td>a/g</td><td>AFF4</td><td>3.9e-10</td><td>++</td><td>++(2.2e-02)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>9</td><td>rs109983693</td><td>7</td><td>59590759</td><td>t/c</td><td>—</td><td>8.6e-12</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>10</td><td>rs109205829</td><td>8</td><td>30083133</td><td>c/g</td><td>NFIB</td><td>2.3e-21</td><td>++</td><td>++(5.4e-01)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>11</td><td>rs42254183</td><td>8</td><td>65107050</td><td>a/g</td><td>—</td><td>2.8e-15</td><td>++</td><td>++(9.5e-05)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>12</td><td>rs134292120</td><td>8</td><td>102635397</td><td>a/t</td><td>ZNF483</td><td>6.3e-13</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>13</td><td>rs41612115</td><td>9</td><td>51985736</td><td>t/c</td><td>—</td><td>9.1e-16</td><td>++</td><td>++(2.3e-01)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>14</td><td>rs110315613</td><td>10</td><td>13737253</td><td>t/c</td><td>—</td><td>5.5e-12</td><td>++</td><td>++(4.5e-05)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>15</td><td>rs110053896</td><td>10</td><td>50918709</td><td>a/g</td><td>FAM81A</td><td>2.1e-15</td><td>++</td><td>++(1.1e-03)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>16</td><td>rs109750969</td><td>10</td><td>91086389</td><td>a/g</td><td>—</td><td>1.3e-15</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>17</td><td>rs109209415</td><td>11</td><td>101301047</td><td>t/c</td><td>LAMC3</td><td>2.9e-15</td><td>++</td><td>++(5.8e-01)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>18</td><td>rs110434838</td><td>13</td><td>40837170</td><td>a/g</td><td>KLZ</td><td>9.3e-18</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>19</td><td>rs109421300</td><td>14</td><td>1801116</td><td>t/c</td><td>DGAT1</td><td>9.7e-283</td><td>++</td><td>++(7.7e-43)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>20</td><td>rs41657812</td><td>14</td><td>5952697</td><td>a/g</td><td>—</td><td>1.2e-29</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>21</td><td>rs41629306</td><td>14</td><td>40812291</td><td>a/c</td><td>HNF4G</td><td>1.2e-16</td><td>++</td><td>++(9.3e-04)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>22</td><td>rs42221575</td><td>15</td><td>60250434</td><td>a/g</td><td>—</td><td>5.4e-19</td><td>++</td><td>++(9.1e-02)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>23</td><td>rs381479851</td><td>16</td><td>1819778</td><td>a/g</td><td>PLEKHA6</td><td>9.7e-14</td><td>++</td><td>++(6.1e-01)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>24</td><td>rs41603770</td><td>16</td><td>28319407</td><td>t/c</td><td>—</td><td>1.9e-13</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>25</td><td>rs10127719</td><td>16</td><td>42843345</td><td>t/c</td><td>DRAXIN</td><td>1.6e-15</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>26</td><td>rs109008971</td><td>16</td><td>63549240</td><td>t/c</td><td>KIAA1614</td><td>9.8e-14</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>27</td><td>rs41666403</td><td>17</td><td>24539327</td><td>a/g</td><td>—</td><td>1.9e-15</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>28</td><td>rs41628722</td><td>17</td><td>60612635</td><td>a/c</td><td>C17H12orf49</td><td>1.6e-13</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>29</td><td>rs419115550</td><td>19</td><td>451642235</td><td>a/g</td><td>GUC1</td><td>2.2e-18</td><td>++</td><td>++(2.2e-07)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>30</td><td>rs109153933</td><td>19</td><td>54820324</td><td>t/c</td><td>—</td><td>8.1e-17</td><td>++</td><td>++(4.0e-06)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>31</td><td>rs436825729</td><td>20</td><td>32670659</td><td>a/g</td><td>FBXO4</td><td>4.0e-41</td><td>??</td><td>(-1.0e+00)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>32</td><td>rs110840223</td><td>23</td><td>30279220</td><td>a/g</td><td>ZKSCAN8</td><td>3.2e-15</td><td>??</td><td>(+4.1e-02)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>33</td><td>rs110836892</td><td>23</td><td>34231434</td><td>c/g</td><td>—</td><td>1.7e-12</td><td>++</td><td>++(3.0e-01)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>34</td><td>rs42094442</td><td>26</td><td>21800857</td><td>a/g</td><td>—</td><td>2.4e-11</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>35</td><td>rs43060720</td><td>26</td><td>30146092</td><td>a/g</td><td>—</td><td>5.2e-15</td><td>++</td><td>++(8.9e-02)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>36</td><td>rs41606803</td><td>26</td><td>44892777</td><td>t/c</td><td>CTBP2</td><td>2.3e-14</td><td>++</td><td>++(8.0e-04)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr></tbody></table>

<footnote n="1">¹Populations in order are MON, NOR, and HOL, “+” and “-” denote positive and negative substitution effects of the alternate allele. “?” indicates that the variant did not segregate in the respective population. The P-value of Cochran's Q-test for heterogeneity of the effect sizes across breeds is given in parentheses, and is significant if P < 0.05.</footnote>
<table><caption>TABLE 2 | Lead SNP variants at 48 QTL for fat yield and percentage in milk.</caption><thead><tr><th>QTL</th><th>Reference SNP<br>cluster ID</th><th>CHR</th><th>BP</th><th>Allele</th><th>Nearby gene</th><th>P<sub>meta</sub></th><th>Effect direction (and evidence of<br>heterogeneity) across three<br>populations<sup>a</sup></th><th>b(se)</th><th>MON</th><th>NOR</th><th>p</th><th>b(se)</th><th>p</th><th>b(se)</th><th>p</th><th>HOL</th></tr></thead><tbody><tr><td>1</td><td>rs42243014</td><td>2</td><td>14950298</td><td>t/c</td><td>NEUROD1</td><td>7.50e-14</td><td>++ + (1.4e-01)</td><td>-0.15(0.03)</td><td>7.40e-08</td><td>-0.14(0.04)</td><td>1.90e-03</td><td>-0.13(0.03)</td><td>4.00e-06</td><td></td><td></td><td></td></tr><tr><td>2</td><td>rs110410153</td><td>2</td><td>26212397</td><td>t/c</td><td>MYO3B</td><td>2.40e-13</td><td>-- -- (4.7e-01)</td><td>-0.07(0.03)</td><td>5.60e-03</td><td>-0.09(0.04)</td><td>1.00e-02</td><td>-0.14(0.02)</td><td>1.60e-10</td><td></td><td></td><td></td></tr><tr><td>3</td><td>rs29025329</td><td>2</td><td>44947955</td><td>a/g</td><td>NM</td><td>2.10e-17</td><td>++ + (3.4e-01)</td><td>-1.50(0.39)</td><td>1.50e-04</td><td>-1.09(0.49)</td><td>2.50e-02</td><td>-2.05(0.28)</td><td>1.30e-13</td><td></td><td></td><td></td></tr><tr><td>4</td><td>rs110253060</td><td>2</td><td>109873484</td><td>a/g</td><td>—</td><td>1.10e-16</td><td>++ + (3.9e-06)</td><td>1.42(0.48)</td><td>2.90e-03</td><td>-0.38(0.52)</td><td>4.70e-01</td><td>4.23(0.46)</td><td>4.40e-20</td><td></td><td></td><td></td></tr><tr><td>5</td><td>rs110900415</td><td>3</td><td>7416886</td><td>a/c</td><td>NOS1AP</td><td>4.20e-14</td><td>++ + (2.6e-03)</td><td>-0.92(0.37)</td><td>1.30e-02</td><td>-0.25(0.47)</td><td>5.90e-01</td><td>-2.17(0.27)</td><td>2.60e-15</td><td></td><td></td><td></td></tr><tr><td>6</td><td>rs110824611</td><td>3</td><td>11040167</td><td>t/c</td><td>OR6N2</td><td>2.30e-29</td><td>++ + (2.6e-01)</td><td>-0.16(0.03)</td><td>4.60e-09</td><td>-0.14(0.05)</td><td>4.10e-03</td><td>-0.25(0.03)</td><td>1.20e-20</td><td></td><td></td><td></td></tr><tr><td>7</td><td>rs42370628</td><td>3</td><td>63136013</td><td>t/c</td><td>—</td><td>1.10e-22</td><td>?? + (1.9e+00)</td><td>-1.18(0.59)</td><td>4.40e-02</td><td>1.63(0.43)</td><td>1.60e-04</td><td>3.00(0.31)</td><td>1.10e-22</td><td></td><td></td><td></td></tr><tr><td>8</td><td>rs109919360</td><td>3</td><td>113845303</td><td>a/g</td><td>USP40</td><td>1.80e-12</td><td>++ + (1.6e-05)</td><td>-1.63(0.41)</td><td>6.50e-05</td><td>1.08(0.70)</td><td>1.20e-01</td><td>-2.17(0.30)</td><td>2.60e-13</td><td></td><td></td><td></td></tr><tr><td>9</td><td>rs209323908</td><td>5</td><td>88807577</td><td>a/c</td><td>ABCC9</td><td>5.60e-16</td><td>-- -- (1.5e-01)</td><td>-2.55(0.46)</td><td>3.30e-08</td><td>-2.39(0.65)</td><td>2.20e-04</td><td>-1.99(0.40)</td><td>5.10e-07</td><td></td><td></td><td></td></tr><tr><td>10</td><td>rs416553559</td><td>5</td><td>94570828</td><td>a/g</td><td>—</td><td>3.50e-31</td><td>-- -- (4.3e-07)</td><td>0.19(0.07)</td><td>4.00e-03</td><td>0.32(0.20)</td><td>1.10e-01</td><td>0.81(0.07)</td><td>3.20e-35</td><td></td><td></td><td></td></tr><tr><td>11</td><td>rs109168860</td><td>7</td><td>20325702</td><td>a/c</td><td>KDM4B</td><td>3.70e-16</td><td>-- -- (4.6e-01)</td><td>0.09(0.03)</td><td>4.70e-04</td><td>0.10(0.04)</td><td>1.50e-02</td><td>0.14(0.02)</td><td>1.70e-12</td><td></td><td></td><td></td></tr><tr><td>12</td><td>rs110304043</td><td>7</td><td>45135080</td><td>t/c</td><td>GRIN3B</td><td>6.10e-15</td><td>-- -- (1.9e-01)</td><td>-0.93(0.40)</td><td>1.80e-02</td><td>-1.72(0.50)</td><td>5.70e-04</td><td>-1.97(0.29)</td><td>8.00e-12</td><td></td><td></td><td></td></tr><tr><td>13</td><td>rs109939693</td><td>7</td><td>59590759</td><td>t/c</td><td>—</td><td>2.80e-17</td><td>-- -- (4.6e-01)</td><td>-2.97(0.93)</td><td>1.40e-03</td><td>-1.69(0.48)</td><td>4.10e-04</td><td>-2.24(0.32)</td><td>1.20e-12</td><td></td><td></td><td></td></tr><tr><td>14</td><td>rs110116017</td><td>7</td><td>102889266</td><td>t/c</td><td>—</td><td>7.20e-10</td><td>-2.(-1.5e-02)</td><td>-0.05(0.04)</td><td>1.80e-01</td><td>0.00(0.04)</td><td>9.20e-01</td><td>-0.15(0.02)</td><td>8.70e-11</td><td></td><td></td><td></td></tr><tr><td>15</td><td>rs425889512</td><td>8</td><td>31586227</td><td>c/g</td><td>LURAP1L</td><td>5.60e-15</td><td>++ + (2.0e-01)</td><td>0.15(0.03)</td><td>2.40e-07</td><td>0.29(0.08)</td><td>2.00e-04</td><td>0.12(0.02)</td><td>1.10e-06</td><td></td><td></td><td></td></tr><tr><td>16</td><td>rs109004526</td><td>8</td><td>76473229</td><td>a/c</td><td>AQP7</td><td>1.50e-13</td><td>-- -- (2.4e-02)</td><td>-0.68(0.40)</td><td>8.90e-02</td><td>-1.74(0.80)</td><td>2.90e-02</td><td>-3.49(0.47)</td><td>1.60e-13</td><td></td><td></td><td></td></tr><tr><td>17</td><td>rs416162829</td><td>8</td><td>93896073</td><td>a/g</td><td>—</td><td>1.10e-31</td><td>++ + (2.4e-02)</td><td>0.10(0.03)</td><td>1.60e-04</td><td>0.17(0.04)</td><td>4.00e-05</td><td>0.25(0.02)</td><td>1.80e-26</td><td></td><td></td><td></td></tr><tr><td>18</td><td>rs41593345</td><td>9</td><td>45503296</td><td>t/g</td><td>PREP</td><td>4.30e-18</td><td>-- -- (1.4e-03)</td><td>-0.04(0.03)</td><td>1.80e-01</td><td>-0.12(0.03)</td><td>5.90e-04</td><td>-0.18(0.02)</td><td>5.30e-18</td><td></td><td></td><td></td></tr><tr><td>19</td><td>rs41612115</td><td>9</td><td>51985736</td><td>t/c</td><td>—</td><td>5.10e-21</td><td>++ + (1.4e-02)</td><td>1.36(0.46)</td><td>2.90e-03</td><td>1.48(0.65)</td><td>2.20e-02</td><td>3.32(0.36)</td><td>8.50e-20</td><td></td><td></td><td></td></tr><tr><td>20</td><td>rs41579916</td><td>9</td><td>79191994</td><td>c/g</td><td>—</td><td>3.80e-15</td><td>++ + (3.8e-02)</td><td>-0.06(0.02)</td><td>1.60e-02</td><td>-0.13(0.07)</td><td>6.10e-02</td><td>-0.16(0.02)</td><td>1.50e-14</td><td></td><td></td><td></td></tr><tr><td>21</td><td>rs110588142</td><td>10</td><td>46388294</td><td>a/t</td><td>HERC1</td><td>2.40e-29</td><td>++ + (8.8e-14)</td><td>-0.23(0.44)</td><td>6.00e-01</td><td>0.69(0.08)</td><td>7.00e-18</td><td>0.31(0.03)</td><td>3.30e-26</td><td></td><td></td><td></td></tr><tr><td>22</td><td>rs41612429</td><td>11</td><td>30220330</td><td>t/c</td><td>FBX011</td><td>1.40e-19</td><td>-- -- (1.8e-01)</td><td>-0.14(0.02)</td><td>7.50e-10</td><td>-0.10(0.03)</td><td>4.70e-03</td><td>-0.14(0.02)</td><td>3.30e-10</td><td></td><td></td><td></td></tr><tr><td>23</td><td>rs41668653</td><td>11</td><td>63467507</td><td>t/c</td><td>RAB1A</td><td>7.50e-48</td><td>++ + (1.6e-03)</td><td>0.11(0.03)</td><td>2.10e-05</td><td>0.24(0.04)</td><td>1.70e-11</td><td>0.30(0.02)</td><td>8.20e-37</td><td></td><td></td><td></td></tr><tr><td>24</td><td>rs110436636</td><td>11</td><td>103289035</td><td>a/g</td><td>—</td><td>4.10e-47</td><td>++ + (1.7e-55)</td><td>0.24(0.02)</td><td>2.00e-24</td><td>-0.31(0.03)</td><td>2.20e-19</td><td>0.34(0.02)</td><td>9.70e-62</td><td></td><td></td><td></td></tr><tr><td>25</td><td>rs109657055</td><td>12</td><td>69106382</td><td>t/c</td><td>—</td><td>1.60e-28</td><td>++ + (3.5e-06)</td><td>0.17(0.03)</td><td>7.80e-10</td><td>-0.01(0.04)</td><td>7.90e-01</td><td>0.21(0.02)</td><td>1.00e-25</td><td></td><td></td><td></td></tr><tr><td>26</td><td>rs108961732</td><td>13</td><td>46433697</td><td>t/c</td><td>ADARB2</td><td>5.50e-14</td><td>++ + (6.7e-06)</td><td>-0.15(0.02)</td><td>1.50e-10</td><td>0.04(0.04)</td><td>2.40e-01</td><td>-0.13(0.02)</td><td>7.30e-10</td><td></td><td></td><td></td></tr><tr><td>27</td><td>rs109421300</td><td>14</td><td>1801116</td><td>t/c</td><td>DGAT1</td><td>2.20e-284</td><td>-- -- (6.1e-42)</td><td>15.24(1.95)</td><td>5.80e-15</td><td>9.12(0.84)</td><td>3.20e-27</td><td>12.77(0.35)</td><td>5.30e-287</td><td></td><td></td><td></td></tr><tr><td>28</td><td>rs109234250</td><td>14</td><td>1802265</td><td>a/g</td><td>DGAT1</td><td>2.6e-322</td><td>++ + (6.3e-53)</td><td>4.13(0.13)</td><td>6.30e-234</td><td>2.57(0.06)</td><td>0.00e+00</td><td>2.82(0.03)</td><td>0.00e+00</td><td></td><td></td><td></td></tr><tr><td>29</td><td>rs417178954</td><td>14</td><td>21129363</td><td>a/g</td><td>PRKDC</td><td>1.30e-25</td><td>++ + (2.6e-08)</td><td>0.04(0.02)</td><td>7.60e-02</td><td>0.04(0.03)</td><td>2.40e-01</td><td>0.26(0.02)</td><td>2.80e-32</td><td></td><td></td><td></td></tr><tr><td>30</td><td>rs109007040</td><td>14</td><td>67443766</td><td>t/c</td><td>STK3</td><td>1.90e-75</td><td>++ + (1.2e-10)</td><td>-0.16(0.02)</td><td>2.20e-11</td><td>-0.09(0.04)</td><td>1.50e-02</td><td>-0.46(0.03)</td><td>2.40e-74</td><td></td><td></td><td></td></tr><tr><td>31</td><td>rs110249976</td><td>15</td><td>53166998</td><td>t/c</td><td>FCHSD2</td><td>1.90e-44</td><td>-- -- (3.2e-09)</td><td>-0.34(0.06)</td><td>1.30e-07</td><td>-0.02(0.03)</td><td>5.20e-01</td><td>-0.47(0.03)</td><td>8.30e-47</td><td></td><td></td><td></td></tr><tr><td>32</td><td>rs42331218</td><td>15</td><td>62236286</td><td>t/c</td><td>—</td><td>4.10e-19</td><td>-- -- (1.7e-01)</td><td>-1.61(0.54)</td><td>2.80e-03</td><td>-1.56(0.47)</td><td>8.70e-04</td><td>-2.46(0.31)</td><td>1.70e-15</td><td></td><td></td><td></td></tr><tr><td>33</td><td>rs378003729</td><td>16</td><td>24591986</td><td>t/c</td><td>—</td><td>3.90e-14</td><td>++ (1.3e-04)</td><td>-2.47(0.37)</td><td>2.10e-11</td><td>0.07(0.55)</td><td>9.00e-01</td><td>-1.53(0.28)</td><td>4.00e-08</td><td></td><td></td><td></td></tr><tr><td>34</td><td>rs41628722</td><td>17</td><td>60612635</td><td>a/c</td><td>—</td><td>3.80e-16</td><td>++ - (6.5e-05)</td><td>1.09(0.38)</td><td>4.40e-03</td><td>-0.09(0.48)</td><td>8.50e-01</td><td>3.00(0.34)</td><td>1.40e-18</td><td></td><td></td><td></td></tr><tr><td>35</td><td>rs41874204</td><td>18</td><td>38471116</td><td>a/g</td><td>ZFHX3</td><td>4.90e-11</td><td>++ - (7.9e-07)</td><td>-0.10(0.45)</td><td>8.30e-01</td><td>0.04(0.47)</td><td>9.20e-01</td><td>-2.35(0.28)</td><td>3.10e-17</td><td></td><td></td><td></td></tr></tbody></table>

(Continued)
<table><caption>TABLE 2 | Continued</caption><thead><tr><th>QTL</th><th>Reference SNP<br>cluster ID</th><th>CHR</th><th>BP</th><th>Allele</th><th>Nearby gene</th><th>P<sub>meta</sub></th><th>Effect direction (and evidence of<br>heterogeneity) across three<br>populations<sup>1</sup></th><th>b(se)</th><th>p</th><th>NOR<br>b(se)</th><th>p</th><th>HOL<br>b(se)</th><th>p</th></tr></thead><tbody><tr><td>36</td><td>rs109882115</td><td>18</td><td>58067310</td><td>a/g</td><td>ENSBTAG00000 9.00e-28</td><td></td><td>++(7.0e-03)</td><td>0.20(0.04)</td><td>1.50e-07</td><td>0.06(0.04)</td><td>1.30e-01</td><td>0.24(0.02)</td><td>2.20e-23</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td>039491</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>37</td><td>rs109751680</td><td>19</td><td>429376859</td><td>a/g</td><td>STAT5B</td><td>1.60e-37</td><td>-+(2.5e-18)</td><td>0.08(0.03)</td><td>4.80e-03</td><td>-0.04(0.04)</td><td>2.50e-01</td><td>0.34(0.02)</td><td>3.50e-53</td></tr><tr><td>38</td><td>rs41921177</td><td>19</td><td>51326750</td><td>a/g</td><td>CCDC67</td><td>6.80e-25</td><td>- - -(1.2e-03)</td><td>5.26(0.87)</td><td>1.40e-09</td><td>0.58(0.91)</td><td>5.30e-01</td><td>2.69(0.30)</td><td>1.00e-19</td></tr><tr><td>39</td><td>rs207867568</td><td>20</td><td>31129820</td><td>t/c</td><td>ENSBTAG00000 6.60e-166</td><td></td><td>??+(1.0e+00)</td><td>-0.64(0.55)</td><td>2.40e-01</td><td>-0.01(0.04)</td><td>8.80e-01</td><td>-0.84(0.03)</td><td>6.60e-166</td></tr><tr><td>40</td><td>rs109262355</td><td>20</td><td>35249040</td><td>t/c</td><td>033187</td><td>166</td><td></td><td>-0.36(0.48)</td><td>4.50e-01</td><td>-1.50(0.58)</td><td>9.50e-03</td><td>-2.56(0.29)</td><td>1.60e-18</td></tr><tr><td>41</td><td>rs137102117</td><td>23</td><td>27450989</td><td>a/g</td><td>FYB</td><td>3.50e-16</td><td>- - -(1.3e-04)</td><td>0.38(0.05)</td><td>7.60e-15</td><td>0.59(0.06)</td><td>4.20e-20</td><td>0.14(0.24)</td><td>5.60e-01</td></tr><tr><td>42</td><td>rs110212032</td><td>24</td><td>37652964</td><td>a/g</td><td>LY6G5B</td><td>4.20e-20</td><td>?+(1.0e+00)</td><td>0.10(0.02)</td><td>3.80e-05</td><td>0.01(0.03)</td><td>8.40e-01</td><td>0.18(0.02)</td><td>2.60e-18</td></tr><tr><td>43</td><td>rs41659834</td><td>26</td><td>14814944</td><td>t/c</td><td>—</td><td>6.40e-19</td><td>++ +(8.0e-04)</td><td>-1.36(0.40)</td><td>6.30e-04</td><td>-1.99(0.47)</td><td>2.10e-05</td><td>-2.05(0.31)</td><td>2.10e-11</td></tr><tr><td>44</td><td>rs42090237</td><td>26</td><td>20290497</td><td>t/c</td><td>MYOF</td><td>1.30e-17</td><td>++ +(4.3e-01)</td><td>0.08(0.03)</td><td>1.90e-03</td><td>0.24(0.05)</td><td>6.60e-06</td><td>0.15(0.02)</td><td>2.40e-13</td></tr><tr><td>45</td><td>rs109330072</td><td>26</td><td>32442734</td><td>t/g</td><td>GOT1</td><td>3.00e-19</td><td>++ +(2.0e-01)</td><td>-4.00(1.01)</td><td>7.50e-05</td><td>-5.09(4.75)</td><td>2.80e-01</td><td>-2.54(0.41)</td><td>4.40e-10</td></tr><tr><td>46</td><td>rs109877780</td><td>27</td><td>26628005</td><td>t/g</td><td>—</td><td>5.10e-13</td><td>- - -(1.7e-01)</td><td>-0.10(0.03)</td><td>9.30e-05</td><td>-0.13(0.04)</td><td>3.00e-03</td><td>-0.25(0.02)</td><td>1.50e-33</td></tr><tr><td>47</td><td>rs109599512</td><td>27</td><td>36117365</td><td>a/g</td><td>GOLGA7</td><td>2.00e-72</td><td>++ +(4.7e-04)</td><td>0.21(0.03)</td><td>9.30e-16</td><td>0.39(0.04)</td><td>8.70e-25</td><td>0.29(0.02)</td><td>4.50e-38</td></tr><tr><td>48</td><td>rs42181608</td><td>29</td><td>29254185</td><td>a/g</td><td>PKNOX2</td><td>3.70e-13</td><td>- - -(1.7e-03)</td><td>-0.62(0.41)</td><td>1.30e-01</td><td>-1.80(0.46)</td><td>1.10e-04</td><td>-1.95(0.30)</td><td>9.10e-11</td></tr></tbody></table>

Some QTL overlapped between Holstein and Montbeliarde, some with same lead SNP, but others with different lead SNP, e.g., BTA15 at 57–61 Mb and BTA25 at 10–11 Mb. The lead SNP for the QTL on BTA15 were rs110049689 (BTA15:57333896, *P* = 1.5 × 10⁻⁷, *MYO7A* gene) in Montbeliarde, and rs41690133 (BTA15:56842162, *P* = 2.4 × 10⁻¹², intergenic) in Holstein. On BTA25, the lead SNP was rs42062121 (BTA25:11015593, SNX29 gene) in both Holstein (*P* = 1.1 × 10⁻⁷) and Montbeliarde (*P* = 2.5 × 10⁻⁶).

### Fat yield and percentage.

We observed QTL on 19 chromosomes where Holstein had 26 QTL, Montbeliarde had 12 QTL, Normande had 5 QTL. Meta-analysis confirmed BTA14:1802265 (*P* = 1.5 × 10⁻⁷⁶⁰, *DGAT1* gene). Other QTL lead SNP whose meta-analysis *P*-values exceeded the GWAS significance threshold include: rs110824611 (BTA3:11040167, *P* = 1.2 × 10⁻⁷, *OR6K6* gene), rs41616289 (BTA8:93896073, *P* = 2.9 × 10⁻⁷, *intergenic*), rs109394729 (BTA10:46486647, *P* = 1.9 × 10⁻⁷, *USP3* gene), rs41668653 (BTA11:63467507, *P* = 1.0 × 10⁻¹⁰, *intergenic*), rs110249976 (BTA15:53166998, *P* = 4.8 × 10⁻⁹, *FCHSD2* gene), rs41921177 (BTA19:51326750, *P* = 8.0 × 10⁻¹⁴, *CCDC57* gene), and rs109599512 (BTA27:36117365, *P* = 9.4 × 10⁻¹⁷, *GOLGA7* gene).

We observed QTL overlapping in Holstein and Montbeliarde and in Holstein and Normande, but with different lead SNP. The two lead SNP for QTL significantly associated in both Holstein and Montbeliarde were BTA12:15173534 (*P* = 3.1 × 10⁻¹³, *NUFIP1* gene, 14.5–15.6 Mb), and rs41921177 (BTA19:51326750, *P* = 1.1 × 10⁻¹⁹, *CCDC57* gene, 51–52 Mb). The three lead SNP for QTL significantly associated in both Holstein and Normande were rs109394729 (BTA10:46486647, *USP3* gene, 43–47 Mb), rs109599512 (BTA27:36117365, *GOLGA7* gene, and 35.5–38.5 Mb), and rs41668653 (BTA11:63467507, *intergenic*, 61.2–64.9 Mb).

### Protein yield and percentage.

We observed QTL on 22 chromosomes where Holstein had 31 QTL, Montbeliarde had 18 QTL, Normande had 22 QTL. We confirmed already known causal variants, BTA6:87199843 (*P* = 9.5 × 10⁻³¹, CSN2 gene), and BTA14:1802265 (*P* = 8.0 × 10⁻⁷⁸, *DGAT1* gene). As presented in **Table 4**, new QTL lead SNP whose meta-analysis *P*-value exceeded the within breed GWAS significance levels included: rs41593345 (BTA9:45303296, *P* = 4.3 × 10⁻¹⁸, *PREP* gene), rs41668653 (BTA11:63467507, *P* = 7.5 × 10⁻⁴⁸, *RAB1A* gene), and rs109882115 (BTA18:58067310, *P* = 9.0 × 10⁻²⁸, *CEACAM18* gene). Two QTL overlaps were observed between Holstein and Montbeliarde, on BTA6 (85–88 Mb) with the most significant lead SNP at BTA6:87199843 (*HSTN*), and on BTA21 with lead SNP. rs29011638 (BTA21:37684315, *P* = 3.8 × 10⁻⁷, *intergenic*). We observed two QTL overlapping between Holstein and Normande on BTA7 (13–14 Mb) with the most significant SNP, rs41568613 (BTA7:13526016, *P* = 7.7 × 10⁻⁶, *intergenic*), and on BTA26 (21.5–24.5 Mb) with lead SNP, rs29014382 (BTA26:22339074, *P* = 1.8 × 10⁻⁹, *intergenic*). We observed two QTL overlapping among all three breeds on BTA16 (1–2 Mb) with lead SNP, rs42450080 (BTA16:1608132, *P* = 1.2

<footnote n="1">¹Populations in order are MON, NOR, and HOL; “+” and “-” denote positive and negative substitution effects of the alternate allele. “?” indicates that the variant did not segregate in the respective population. The P-value of Cochran's Q-test for heterogeneity of the effect sizes across breeds is given in parentheses, and is significant if P < 0.05.</footnote>
<table><caption>TABLE 3 | Lead SNP variants at 29 QTL for protein yield and percentage in milk.</caption><thead><tr><th rowspan="2">QTL</th><th rowspan="2">Reference SNP cluster ID</th><th rowspan="2">CHR</th><th rowspan="2">BP</th><th rowspan="2">Allele</th><th rowspan="2">Nearby gene</th><th rowspan="2">P<sub>meta</sub></th><th rowspan="2">Effect direction (and evidence of heterogeneity) across three populations<sup>1</sup></th><th colspan="2">MON</th><th colspan="2">NOR</th><th colspan="2">HOL</th></tr><tr><th>b(se)</th><th>p</th><th>b(se)</th><th>p</th><th>b(se)</th><th>p</th></tr></thead><tbody><tr><td>1</td><td>rs29025329</td><td>2</td><td>44947955</td><td>a/g</td><td>NMI</td><td>2.94e-18</td><td>++ + (9.6e-02)</td><td>-0.82(0.32)</td><td>1.0e-02</td><td>-1.28(0.36)</td><td>3.6e-04</td><td>-1.61(0.21)</td><td>4.8e-15</td></tr><tr><td>2</td><td>rs42889305</td><td>2</td><td>70680985</td><td>t/g</td><td>—</td><td>1.15e-18</td><td>— + (3.2e-04)</td><td>0.36(0.33)</td><td>2.8e-01</td><td>1.52(0.43)</td><td>4.5e-04</td><td>1.84(0.20)</td><td>3.0e-19</td></tr><tr><td>3</td><td>rs110532777</td><td>3</td><td>21640254</td><td>t/c</td><td>POLR3C</td><td>8.05e-11</td><td>— + (3.8e-02)</td><td>0.49(0.30)</td><td>1.0e-01</td><td>0.52(0.35)</td><td>1.5e-01</td><td>1.38(0.21)</td><td>3.3e-11</td></tr><tr><td>4</td><td>rs208464932</td><td>3</td><td>48855688</td><td>a/g</td><td>SLC44A3</td><td>1.17e-13</td><td>?? + (1.0e+00)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>2.02(0.27)</td><td>1.2e-13</td></tr><tr><td>5</td><td>rs42737927</td><td>3</td><td>113907407</td><td>t/c</td><td>UGT1A1</td><td>3.31e-12</td><td>++ + (4.3e-05)</td><td>-0.49(0.34)</td><td>1.4e-01</td><td>0.05(0.35)</td><td>8.9e-01</td><td>-1.73(0.21)</td><td>3.6e-16</td></tr><tr><td>6</td><td>rs43703011</td><td>6</td><td>87181619</td><td>t/c</td><td>CSN2</td><td>9.84e-31</td><td>— + (6.4e-03)</td><td>-1.92(0.42)</td><td>5.9e-06</td><td>-3.71(0.51)</td><td>2.2e-13</td><td>-2.67(0.32)</td><td>1.2e-16</td></tr><tr><td>7</td><td>rs41609745</td><td>7</td><td>65358446</td><td>t/c</td><td>—</td><td>1.89e-14</td><td>++ + (2.4e-04)</td><td>0.11(0.32)</td><td>7.3e-01</td><td>1.52(0.39)</td><td>8.6e-05</td><td>1.57(0.20)</td><td>1.0e-14</td></tr><tr><td>8</td><td>rs109205829</td><td>8</td><td>30088133</td><td>c/g</td><td>NFIB</td><td>1.79e-17</td><td>++ + (3.2e-01)</td><td>1.22(0.32)</td><td>1.2e-04</td><td>1.69(0.36)</td><td>2.7e-06</td><td>1.29(0.21)</td><td>7.4e-10</td></tr><tr><td>9</td><td>rs42254183</td><td>8</td><td>65107050</td><td>a/g</td><td>—</td><td>5.25e-12</td><td>++ + (1.2e-03)</td><td>-0.26(0.43)</td><td>5.5e-01</td><td>-1.13(0.53)</td><td>3.4e-02</td><td>-2.20(0.29)</td><td>6.7e-14</td></tr><tr><td>10</td><td>rs109630763</td><td>8</td><td>103526477</td><td>t/c</td><td>—</td><td>1.85e-12</td><td>— + (1.5e-03)</td><td>-5.37(1.26)</td><td>2.0e-05</td><td>0.20(0.35)</td><td>5.6e-01</td><td>-1.79(0.27)</td><td>3.0e-11</td></tr><tr><td>11</td><td>rs41612115</td><td>9</td><td>51985736</td><td>t/c</td><td>—</td><td>1.64e-22</td><td>++ + (4.8e-02)</td><td>1.32(0.37)</td><td>3.8e-04</td><td>1.19(0.48)</td><td>1.3e-02</td><td>2.46(0.27)</td><td>1.1e-19</td></tr><tr><td>12</td><td>rs41647455</td><td>10</td><td>15333419</td><td>a/g</td><td>—</td><td>1.22e-16</td><td>— + (1.1e-07)</td><td>-0.32(0.46)</td><td>4.9e-01</td><td>-0.23(0.38)</td><td>5.3e-01</td><td>-2.12(0.21)</td><td>1.8e-23</td></tr><tr><td>13</td><td>rs43625971</td><td>10</td><td>47689396</td><td>t/c</td><td>TLN2</td><td>1.47e-13</td><td>++ + (9.7e-03)</td><td>-0.36(0.32)</td><td>2.7e-01</td><td>-1.25(0.35)</td><td>3.7e-04</td><td>-1.70(0.24)</td><td>1.5e-12</td></tr><tr><td>14</td><td>rs109186903</td><td>10</td><td>94083525</td><td>t/c</td><td>—</td><td>2.51e-18</td><td>?? + (1.0e+00)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>2.23(0.26)</td><td>2.5e-18</td></tr><tr><td>15</td><td>rs41684057</td><td>13</td><td>17104043</td><td>t/c</td><td>PRKCQ</td><td>1.48e-10</td><td>++ + (8.7e-02)</td><td>1.14(0.83)</td><td>1.7e-01</td><td>1.19(0.47)</td><td>1.1e-02</td><td>1.30(0.21)</td><td>8.6e-10</td></tr><tr><td>16</td><td>rs110991150</td><td>13</td><td>40815971</td><td>c/g</td><td>—</td><td>3.29e-20</td><td>++ + (4.0e-02)</td><td>1.59(0.46)</td><td>5.1e-04</td><td>7.27(3.66)</td><td>4.7e-02</td><td>2.64(0.30)</td><td>4.1e-18</td></tr><tr><td>17</td><td>rs109421300</td><td>14</td><td>1801116</td><td>t/c</td><td>HSF1</td><td>8.03e-78</td><td>++ + (8.2e-11)</td><td>-8.08(1.59)</td><td>3.7e-07</td><td>-2.70(0.62)</td><td>1.5e-05</td><td>-4.91(0.26)</td><td>3.1e-78</td></tr><tr><td>18</td><td>rs41617796</td><td>14</td><td>37670101</td><td>t/c</td><td>—</td><td>3.11e-11</td><td>++ + (1.6e-01)</td><td>0.63(0.31)</td><td>4.0e-02</td><td>0.74(0.40)</td><td>6.2e-02</td><td>1.71(0.27)</td><td>2.4e-10</td></tr><tr><td>19</td><td>rs109020671</td><td>14</td><td>60893281</td><td>a/g</td><td>—</td><td>2.25e-12</td><td>++ + (4.4e-01)</td><td>-0.81(0.33)</td><td>1.5e-02</td><td>-1.25(0.43)</td><td>3.8e-03</td><td>-1.29(0.21)</td><td>1.5e-09</td></tr><tr><td>20</td><td>rs42221575</td><td>15</td><td>60250434</td><td>a/g</td><td>—</td><td>2.21e-14</td><td>++ + (1.4e-02)</td><td>-0.63(0.46)</td><td>1.7e-01</td><td>-1.28(0.39)</td><td>1.1e-03</td><td>-1.49(0.20)</td><td>1.7e-13</td></tr><tr><td>21</td><td>rs109053739</td><td>15</td><td>77590401</td><td>g/c,t</td><td>CKAP5</td><td>6.60e-13</td><td>++ + (3.1e-03)</td><td>-2.24(0.64)</td><td>4.6e-04</td><td>0.12(0.63)</td><td>8.4e-01</td><td>-1.93(0.27)</td><td>9.7e-13</td></tr><tr><td>22</td><td>rs41628722</td><td>17</td><td>60612635</td><td>a/c</td><td>—</td><td>1.48e-16</td><td>— + (2.8e-03)</td><td>0.94(0.31)</td><td>2.7e-03</td><td>0.26(0.36)</td><td>4.7e-01</td><td>2.12(0.25)</td><td>4.7e-17</td></tr><tr><td>23</td><td>rs110031071</td><td>19</td><td>51767413</td><td>a/g</td><td>PDE6G</td><td>1.59e-09</td><td>— + (5.7e-02)</td><td>-0.36(0.31)</td><td>2.3e-01</td><td>-1.34(0.36)</td><td>1.6e-04</td><td>-1.06(0.21)</td><td>2.6e-07</td></tr><tr><td>24</td><td>rs436825729</td><td>20</td><td>32670639</td><td>a/g</td><td>—</td><td>3.66e-304</td><td>?? + (1.0e+00)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>-0.59(0.01)</td><td>0.00</td></tr><tr><td>25</td><td>rs41592357</td><td>22</td><td>24219999</td><td>t/c</td><td>—</td><td>5.63e-16</td><td>— + (4.3e-01)</td><td>-1.22(0.35)</td><td>4.6e-04</td><td>-1.50(0.35)</td><td>1.8e-05</td><td>-1.24(0.20)</td><td>1.5e-09</td></tr><tr><td>26</td><td>rs385558887</td><td>23</td><td>29349996</td><td>a/g</td><td>—</td><td>1.20e-17</td><td>?? + (1.2e-07)</td><td>0.08(0.31)</td><td>7.9e-01</td><td>—</td><td>—</td><td>2.81(0.28)</td><td>9.0e-24</td></tr><tr><td>27</td><td>rs41624917</td><td>26</td><td>15338386</td><td>a/g</td><td>PLCE1</td><td>2.09e-11</td><td>++ + (5.0e-03)</td><td>0.87(0.39)</td><td>2.6e-02</td><td>0.14(0.51)</td><td>7.9e-01</td><td>1.46(0.21)</td><td>1.2e-12</td></tr><tr><td>28</td><td>rs109839618</td><td>26</td><td>42403875</td><td>a/g</td><td>TACC2</td><td>8.56e-13</td><td>++ + (1.6e-02)</td><td>-0.64(0.31)</td><td>3.9e-02</td><td>-0.47(0.40)</td><td>2.4e-01</td><td>-1.91(0.26)</td><td>2.2e-13</td></tr><tr><td>29</td><td>rs41686102</td><td>29</td><td>41989397</td><td>t/c</td><td>SLC22A6</td><td>1.38e-51</td><td>++ + (1.1e-16)</td><td>-0.03(0.02)</td><td>3.3e-02</td><td>-0.06(0.03)</td><td>3.0e-02</td><td>-2.3(0.01)</td><td>1.4e-65</td></tr></tbody></table>

<footnote n="1">¹Populations in order are MON, NOR, and HOL, “+” and “-” denote positive and negative substitution effects of the alternate allele. “?” indicates that the variant did not segregate in the respective population. The P-value of Cochran's Q-test for heterogeneity of the effect sizes across breeds is given in parentheses, and is significant if P < 0.05.</footnote>
<table><caption>TABLE 4 | Lead SNP variants at 13 QTL for success rate of insemination of lactating cows.</caption><thead><tr><th rowspan="2">QTL</th><th rowspan="2">Reference SNP<br>cluster ID</th><th rowspan="2">CHR</th><th rowspan="2">BP</th><th rowspan="2">Allele</th><th rowspan="2">Nearby<br>gene</th><th rowspan="2">P<sub>meta</sub></th><th rowspan="2">Effect direction (and evidence of<br>heterogeneity) across three<br>populations<sup>1</sup></th><th colspan="2">MON</th><th colspan="2">NOR</th><th colspan="2">HOL</th></tr><tr><th>b(se)</th><th>p</th><th>b(se)</th><th>p</th><th>b(se)</th><th>p</th></tr></thead><tbody><tr><td>1</td><td>rs109193035</td><td>2</td><td>24147156</td><td>a/g</td><td>ITGA6</td><td>1.07e-06</td><td>+?+(3.3e-01)</td><td>0.02(0.01)</td><td>6.3e-02</td><td>0.01(0.01)</td><td>4.5e-01</td><td>0.07(0.02)</td><td>3.9e-06</td></tr><tr><td>2</td><td>rs29010833</td><td>2</td><td>6688588</td><td>t/g</td><td>—</td><td>1.07e-06</td><td>+ + +(5.5e-01)</td><td>-0.01(0.00)</td><td>1.2e-01</td><td>-0.01(0.00)</td><td>1.9e-02</td><td>-0.01(0.00)</td><td>3.5e-05</td></tr><tr><td>3</td><td>rs110027085</td><td>2</td><td>101307425</td><td>a/g</td><td>—</td><td>3.00e-07</td><td>+ + +(3.6e-01)</td><td>-0.01(0.00)</td><td>1.5e-01</td><td>-0.01(0.01)</td><td>5.9e-02</td><td>-0.02(0.00)</td><td>1.9e-06</td></tr><tr><td>4</td><td>rs41897801</td><td>8</td><td>23964538</td><td>a/g</td><td>MLL73</td><td>6.39e-06</td><td>+ + +(1.7e-01)</td><td>0.01(0.00)</td><td>1.2e-02</td><td>0.00(0.00)</td><td>9.8e-01</td><td>0.01(0.00)</td><td>2.7e-05</td></tr><tr><td>5</td><td>rs425288799</td><td>8</td><td>59245157</td><td>a/c</td><td>—</td><td>6.07e-06</td><td>+ + +(6.5e-01)</td><td>-0.01(0.00)</td><td>9.5e-02</td><td>-0.01(0.01)</td><td>1.4e-02</td><td>-0.01(0.00)</td><td>4.0e-04</td></tr><tr><td>6</td><td>rs109483390</td><td>9</td><td>76868290</td><td>a/g</td><td>—</td><td>5.24e-10</td><td>+ + +(3.1e-01)</td><td>0.01(0.00)</td><td>2.9e-02</td><td>0.01(0.00)</td><td>8.6e-02</td><td>0.02(0.00)</td><td>8.3e-09</td></tr><tr><td>7</td><td>rs42412333</td><td>10</td><td>39278374</td><td>a/g</td><td>—</td><td>1.19e-08</td><td>+ + +(5.1e-01)</td><td>-0.01(0.01)</td><td>4.3e-02</td><td>-0.01(0.00)</td><td>2.7e-03</td><td>-0.01(0.00)</td><td>5.2e-06</td></tr><tr><td>8</td><td>rs109269005</td><td>11</td><td>67347867</td><td>a/g</td><td>ANTXR1</td><td>8.39e-07</td><td>— + − (2.9e-01)</td><td>0.01(0.00)</td><td>1.5e-02</td><td>0.01(0.02)</td><td>5.9e-01</td><td>0.04(0.01)</td><td>6.0e-06</td></tr><tr><td>9</td><td>rs41620105</td><td>11</td><td>82643114</td><td>t/c</td><td>—</td><td>3.93e-07</td><td>+ + +(2.0e-01)</td><td>0.00(0.00)</td><td>3.1e-01</td><td>-0.01(0.00)</td><td>1.4e-02</td><td>-0.01(0.00)</td><td>3.0e-06</td></tr><tr><td>10</td><td>rs110425867</td><td>14</td><td>8264685</td><td>t/c</td><td>ZFAT</td><td>1.63e-11</td><td>+ + +(3.2e-01)</td><td>-0.02(0.01)</td><td>2.8e-04</td><td>-0.01(0.01)</td><td>2.1e-01</td><td>-0.03(0.01)</td><td>9.9e-09</td></tr><tr><td>11</td><td>rs41568388</td><td>15</td><td>39634819</td><td>a/g</td><td>—</td><td>4.94e-07</td><td>+ + +(1.9e-02)</td><td>0.00(0.00)</td><td>5.0e-01</td><td>0.00(0.00)</td><td>4.1e-01</td><td>0.02(0.00)</td><td>1.5e-08</td></tr><tr><td>12</td><td>rs109269484</td><td>20</td><td>40301530</td><td>a/c</td><td>ADAMTS12</td><td>4.76e-05</td><td>+ + +(5.2e-01)</td><td>0.01(0.00)</td><td>2.3e-03</td><td>0.01(0.01)</td><td>1.7e-01</td><td>0.01(0.00)</td><td>9.7e-03</td></tr><tr><td>13</td><td>rs42049541</td><td>24</td><td>29349786</td><td>t/g</td><td>—</td><td>1.44e-06</td><td>+ + +(3.8e-05)</td><td>0.00(0.00)</td><td>3.5e-01</td><td>-0.01(0.00)</td><td>1.5e-01</td><td>0.02(0.00)</td><td>1.8e-10</td></tr></tbody></table>

× 10⁻¹⁰, intergenic), and on BTA23 (50–51 Mb) with lead SNP, rs42507912 (BTA23:51118713, *P* = 1.3 × 10⁻¹¹, *GMDS* gene). QTL overlap between fat and protein was evident, e.g., on BTA11 at 61–64 Mb with its lead SNP at rs41668653 (BTA11:63467507, *P* = 7.2 × 10⁻¹³, intergenic), and on BTA14 at 1.6–1.9 Mb with the most significant lead SNP being BTA14:1802265 (*P* = 2.6 × 10⁻³², *DGAT1* gene).

## Fertility

For the success rate at insemination of lactating cows, we observed QTL on 9 chromosomes where Holstein had 13 QTL, Montbeliarde had 2 QTL, Normande had 1 QTL. As presented in Table 5, the lead SNP for new QTL from meta-analysis that exceeded the genome wide threshold (*P* < 1.03 × 10⁻⁶) included rs109483390 (BTA9:76868290, *P* = 5.2 × 10⁻¹⁰, 0.1 Mb upstream of *TNFAIP3* gene and 0.07 Mb downstream of *PERP* gene), rs42412333 (BTA10:39278374, *P* = -1.2 × 10⁻⁸, 38–40 Mb, 0.45 Mb downstream of *RPL10L* gene), and rs110425867 (BTA14:8264685, *P* = 1.6 × 10⁻¹¹, *ZFAT* gene 8.1–8.3 Mb). Lead SNP and QTL intervals for QTL observed only in Holstein were: rs41589904 (BTA8:76582220, *P* = 7.5 × 10⁻⁹, *UBE2R2* gene, 74–77 Mb), rs41593363 (BTA9:53806121, *P* = 1.3 × 10⁻⁸, 49–55 Mb, and 6.5 kb upstream of *GPR63* gene), rs110543856 (BTA18:48150900, *P* = 1.1 × 10⁻¹⁸, *SIPA1L3* gene, 43–48 Mb), rs110588160 (BTA21:47117318, *P* = 2.6 × 10⁻¹⁴, intergenic, 43–48 Mb), rs29023151 (BTA22:23303686, *P* = 1.6 × 10⁻¹¹, *IL5RA* gene, 22–24 Mb), rs109629413 (BTA24:33624891, *P* = 7.1 × 10⁻¹⁵, *TMEM241* gene, 28–33 Mb), and two, QTL on BTA26: rs43158237 (BTA26:6741237, *P* = 2 × 10⁻¹¹, intergenic, 6.7–6.9 Mb) and rs42096924 (BTA26:33443386, *P* = 2.8 × 10⁻¹¹, intergenic, 24–33 Mb). Lead SNP and QTL intervals for QTL observed only in Montbeliarde were: rs41657531 (BTA9:39816061, 1.1 × 10⁻⁸, *RPF2* gene, 38–44 Mb), and rs42417896 (BTA11:34068419, 4.6 × 10⁻⁶, intergenic, 31–39 Mb). The suggestive fertility QTL observed in Normande was rs41570140 (BTA11:30065164, 1.6 × 10⁻⁴, *FBX011* gene, 27–30 Mb).

## Stature

We identified 16 QTL on 11 chromosomes associated with stature in the three breeds (Table 5). The most significant variant was observed in the Montbeliarde breed and was BTA14:25006125 (*P* = 8.6 × 10⁻¹⁴⁰ *PLAG1* gene). Seven novel QTL for stature displayed heterogeneity across the three breeds, e.g., rs43121344 (BTA7:68281468, intergenic and 0.14 Mb toward the 3'-flanking region of *MRPL22* gene) was lead SNP in all three breeds with *P*<sub>NOR</sub> = 3.0 × 10⁻⁴, *P*<sub>MON</sub> = 4.5 × 10⁻⁸, *P*<sub>HOL</sub> = 1.1 × 10⁻⁹, and *P*<sub>META</sub> = 7.6 × 10⁻¹⁹. Lead SNP for breed specific QTL in Montbeliarde include: rs41645172 (BTA2:109778246, *P* = 3.9 × 10⁻¹⁷, *DOCK10* gene), rs109860141 (BTA10:51558941, *P* = 2.4 × 10⁻²⁰ *FAM63B* gene), BTA14:25006125 (*P* = 8.6 × 10⁻¹⁴⁰, *PLAG1* gene), and BTA20:39917111 (*P* = 5.0 × 10⁻²², *ADAMTS12* gene). Lead SNP for QTL observed in Normande include rs41644660 (BTA26:19460476, *P* = 3.5 × 10⁻⁸, 73 kb downstream of *HPS1* gene). Two QTL lead SNP observed in Holstein include rs110652403 (BTA7:90460692, *P* = 4.8 × 10⁻¹²,

<footnote n="1">¹Populations in order are MON, NOR, and HOL, “+” and “-” denote positive and negative substitution effects of the alternate allele. “?” indicates that the variant did not segregate in the respective population. The P-value of Cochran's Q-test for heterogeneity of the effect sizes across breeds is given in parentheses, and is significant if P < 0.05.</footnote>
<table><caption>TABLE 5 | Lead SNP variants at 16 QTL for stature.</caption><thead><tr><th rowspan="2">QTL</th><th rowspan="2">Reference SNP cluster ID</th><th rowspan="2">CHR</th><th rowspan="2">BP</th><th rowspan="2">Allele</th><th rowspan="2">Nearby gene</th><th rowspan="2">P<sub>meta</sub></th><th rowspan="2">Effect direction (and evidence of heterogeneity) across three populations<sup>1</sup></th><th colspan="2">MON</th><th colspan="2">NOR</th><th colspan="2">HOL</th></tr><tr><th>b(se)</th><th>p</th><th>b(se)</th><th>p</th><th>b(se)</th><th>p</th></tr></thead><tbody><tr><td>1</td><td>rs41645172</td><td>2</td><td>109778246</td><td>t/c</td><td>—</td><td>3.3e-10</td><td>++ + (8.1e-10)</td><td>0.31(0.04)</td><td>3.9e-17</td><td>0.22(0.07)</td><td>2.3e-03</td><td>0.18(0.01)</td><td>2.7e-01</td></tr><tr><td>2</td><td>rs43121344</td><td>7</td><td>68281468</td><td>t/c</td><td>—</td><td>7.6e-19</td><td>— — (4.6e-01)</td><td>0.30(0.05)</td><td>4.5e-08</td><td>0.44(0.12)</td><td>3.0e-05</td><td>0.04(0.01)</td><td>1.1e-09</td></tr><tr><td>3</td><td>rs110652403</td><td>7</td><td>90460692</td><td>t/c</td><td>—</td><td>1.2e-11</td><td>— — (1.3e-02)</td><td>—0.04(0.04)</td><td>3.0e-01</td><td>—0.26(0.01)</td><td>1.5e-02</td><td>—0.05(0.01)</td><td>4.8e-12</td></tr><tr><td>4</td><td>rs109703572</td><td>9</td><td>46739973</td><td>t/c</td><td>—</td><td>1.8e-23</td><td>++ + (6.3e-03)</td><td>—0.09(0.04)</td><td>1.8e-02</td><td>—0.22(0.06)</td><td>7.7e-05</td><td>—0.06(0.01)</td><td>4.9e-21</td></tr><tr><td>5</td><td>rs41664341</td><td>9</td><td>80859295</td><td>a/g</td><td>—</td><td>5.8e-22</td><td>— — (2.6e-02)</td><td>0.15(0.04)</td><td>5.2e-05</td><td>0.10(0.06)</td><td>8.7e-02</td><td>0.06(0.01)</td><td>2.5e-19</td></tr><tr><td>6</td><td>rs109860141</td><td>10</td><td>51558941</td><td>a/g</td><td>FAM63B</td><td>5.9e-16</td><td>++ + (2.6e-09)</td><td>0.34(0.04)</td><td>2.4e-20</td><td>0.22(0.06)</td><td>4.2e-04</td><td>0.02(0.01)</td><td>7.9e-03</td></tr><tr><td>7</td><td>rs110128058</td><td>10</td><td>58318595</td><td>a/g</td><td>GNB5</td><td>3.5e-18</td><td>++ + (3.5e-09)</td><td>0.40(0.04)</td><td>3.6e-22</td><td>0.22(0.07)</td><td>2.5e-03</td><td>0.04(0.01)</td><td>6.4e-04</td></tr><tr><td>8</td><td>rs43430263</td><td>12</td><td>53241975</td><td>t/g</td><td>SLA/N1</td><td>6.5e-16</td><td>++ - (1.4e-04)</td><td>0.19(0.04)</td><td>5.7e-07</td><td>—0.04(0.06)</td><td>4.6e-01</td><td>0.06(0.01)</td><td>3.3e-14</td></tr><tr><td>9</td><td>rs29019889</td><td>13</td><td>33366445</td><td>a/g</td><td>NSU/N6</td><td>6.4e-16</td><td>++ + (9.0e-03)</td><td>0.21(0.03)</td><td>4.5e-10</td><td>0.10(0.10)</td><td>3.3e-01</td><td>0.04(0.01)</td><td>3.5e-09</td></tr><tr><td>10</td><td>rs134215421</td><td>14</td><td>25006125</td><td>t/c</td><td>PLAG1</td><td>8.6e-140</td><td>-?? (1.0e+00)</td><td>1.70(0.07)</td><td>8.6e-140</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>11</td><td>rs41584869</td><td>19</td><td>17569423</td><td>a/g</td><td>SPACA3</td><td>2.4e-11</td><td>— — (1.3e-01)</td><td>0.32(0.06)</td><td>4.7e-07</td><td>0.18(0.07)</td><td>6.0e-03</td><td>0.03(0.01)</td><td>6.9e-05</td></tr><tr><td>12</td><td>rs110291686</td><td>19</td><td>43941099</td><td>a/c</td><td>U2; ARLL4D</td><td>1.7e-15</td><td>— — (3.3e-01)</td><td>0.20(0.04)</td><td>1.4e-07</td><td>0.14(0.06)</td><td>1.7e-02</td><td>0.04(0.01)</td><td>1.3e-08</td></tr><tr><td>13</td><td>rs109082193</td><td>20</td><td>18413770</td><td>a/g</td><td>ERCC8</td><td>1.3e-12</td><td>++ + (7.0e-01)</td><td>-0.19(0.05)</td><td>5.9e-04</td><td>-0.11(0.05)</td><td>3.9e-02</td><td>-0.05(0.01)</td><td>3.2e-09</td></tr><tr><td>14</td><td>rs385826566</td><td>20</td><td>39917111</td><td>a/c</td><td>ADAMTS12</td><td>5.0e-22</td><td>-?? (1.0e+00)</td><td>-0.37(0.04)</td><td>5.0e-22</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>15</td><td>rs133559545</td><td>21</td><td>8262352</td><td>a/g</td><td>IGF1R</td><td>2.4e-17</td><td>++?? (1.0e+00)</td><td>-0.30(0.04)</td><td>2.4e-17</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>16</td><td>rs41644660</td><td>26</td><td>19460476</td><td>t/c</td><td>—</td><td>1.1e-27</td><td>++?? (3.2e-07)</td><td>-0.27(0.05)</td><td>3.5e-08</td><td>-0.76(0.07)</td><td>9.3e-27</td><td>—</td><td>—</td></tr></tbody></table>

<footnote n="1">¹Populations in order are MON, NOR, and HOL, “+” and “-” denote positive and negative substitution effects of the alternate allele. “?” indicates that the variant did not segregate in the respective population. The P-value of Cochran's Q-test for heterogeneity of the effect sizes across breeds is given in parentheses, and is significant if P < 0.05.</footnote>
<figure><img src="image_5.png" /></figure>

**FIGURE 3** | Heterogeneity of SNP association with fat and protein in milk across three cattle breeds. (A) Composite Manhattan plot that shows the association of 43,421 imputed variants including 5,609 pre-selected sequence variants with fat and protein in the meta-analysis. The composite Manhattan plot summarizes the results of the meta-analyses with each dot representing the more significant *P*-value that was observed across both traits. Light green represents sequence variants with *P* < 1 × 10⁻⁶. (B) Quantile-quantile plot of the meta-analyses. Green and brown color represent *P*-value of 43,421 imputed variants for fat and protein, respectively. (C) Overview of 59 QTL that were significant at *P* < 1 × 10⁻⁶ in the meta-analysis and within-breed association studies. Each column represents one of 29 *Bos taurus* autosomes. Row colors are breed specific with the top row being overall meta-analysis QTL. Filled squares indicate that QTL were significant in the respective breeds and chromosome. (D) Allelic substitution effects of 36 QTL on fat and protein standardized with the phenotypic standard deviations. The vertical axis is protein, and the horizontal axis is fat.

0.16 Mb upstream of *MEF2C* gene, and 0.3 Mb downstream of *TMEM161B* gene), and rs43430263 (BTA12:53241975, *P* = 3.3 × 10⁻¹⁴, *SLAIN1* gene).

Meta-Analysis Heterogeneity Among Three Breeds

As an example of the heterogeneity observed from a meta-analysis of the three breeds, we observed peaks from a combined Manhattan plot for fat and protein percentage (Figures 3A,B) and standardized the allelic substitution effects by the phenotypic standard deviation of fat and protein percentage meta-analysis results. We observed heterogeneity (Cochran's Q, *P* < 0.05) for 36 QTL from the meta-analysis in three breeds (Figure 3D). We then overlay observed breed specific QTL for these two traits and Holstein had 26 QTL in 24 chromosomes, Montbeliarde had 17 QTL in 14 chromosomes, and Normande had 16 QTL in 7 chromosomes (Figure 3C).

DISCUSSION

Our meta-analysis of association summary statistics for seven
traits across three French dairy cattle breeds discovered 120 QTL
including 13 QTL that had not been detected at *P* < 1.03 ×
10⁻⁶ in the within-breed analyses. In agreement with previous
studies (Pausch et al., 2017) our results show that combining
GWAS summary data from several breeds increases the power
of association studies. However, for lowly heritable traits, and
for distinct populations, a meta-analysis may not be as robust

as the within-breed GWAS (Raven et al., 2014). For instance,
within-breed GWAS showed associated SNP for fertility had no
overlap between breeds probably due to small effects per locus
(i.e., low power) and possibly the action of many, reasonably
rare recessive lethal alleles. Holstein, with larger sample size,
had associated SNP in strong LD and thus had stronger fertility
predicted effects. In some instances, the meta-analysis increased
the power, and consequently, the confidence interval for some
fertility related QTL became visible, e.g., QTL on BTA2, 8 and 20
(Table 4). However, in other instances, the QTL detection power
reduced when the QTL are private, e.g., the meta-analysis for
fertility lost power compared to the within-breed GWAS when
the Normande breed (with the smallest number of animals) was
included in the analysis, e.g., QTL on BTA15 and 24 (Table 4).

One advantage of a meta-analysis over within-breed GWAS
is that it is expected to result in more precision in QTL position
due to the breakdown of long-range LD due to many generations
of recombination since the breeds were separated (Raven et al.,
2014; Pausch et al., 2016). In this study, we performed a within-
breed GWAS first, followed by meta-analyses of GWAS summary
statistics from three breeds using a weighted Z-score model to
get the *P*-value for each SNP. In our study, the meta-analysis
increased the *P*-value for some variants, probably because the
QTL was not segregating across populations. Meta-analysis
decreased *P*-value for some variants, and consequently, some
QTL peaks became more distinct and narrower. For example, a
region on BTA6:87296809 was strongly associated with protein
yield and mapped near the protein-coding genes *ODAM* and
CSN3 in a 120 kb region. Bovine odontogenic, ameloblast-associated (ODAM) participates in structuring the extracellular matrix to attach epithelial cells to mineralized surfaces thus forming a protective seal that is antagonistic to bacterial invasion (Fouillen et al., 2017). On the other hand, Casein kappa (CSN3) is associated with daily milk yield (Bonfatti et al., 2011; Bartonova et al., 2012) as well as better cheese properties (Di Gregorio et al., 2017).

Strongly associated SNP ($P < 1 \times 10^{-50}$) were identified in all three French dairy breeds. Among these were several SNP that had previously been described such as those close to the DGAT1-K232A mutation on BTA14 (Grisart et al., 2002). These SNP were highly associated with all milk production traits in all three breeds. A very significant sequence variant at BTA5:93948357 (rs209372883) was associated with protein and fat in Holstein and Montbeliarde; it has previously been ascribed to a variant in *MGST1* (microsomal glutathione S-transferase 1; Raven et al., 2016), an upstream intron variant. *MGST1* is also an inflammation response gene which is highly expressed through pregnancy and lactation (Church et al., 2012). Previous studies in Japanese Black cattle (Wang et al., 2005) suggested upregulation of *MGST1* during adipocyte development in the longissimus muscle. Slightly upstream, a highly significant peak was centered within *SLC15A5*, a protein-coding gene. Our results suggest that both *MGST1* and *SLC15A5* may contain variants affecting milk production. Raven et al. (2014) identified a fat yield QTL within 3kbp of *MGST1* using a multi-breed analysis and we report a QTL which was within 2kbp of the genic region of *MGST1* using meta-analysis. Similar studies in Canadian Holstein, and using a similar density chip (50 k), reported a QTL 200 bp from that reported in our study (Li et al., 2010). The same Canadian Holstein study reported several fat yield candidate genes including *SLC2A3* and *GDF3* at 101.7–101.8 Mb and *LRP6*, *EMP1*, and *DUSP16* at 97–98 Mb on BTA5 (Li et al., 2010). Their results agrees with our meta-analysis results. Another very significant QTL for milk yield in Holstein was identified on BTA20 at ~38 Mb. The lead SNP of this QTL was located near the *PRLR* gene which has previously been associated with fertility in Holstein cows (Leyva-Corona et al., 2018; Marete et al., 2018b). The majority of the associated SNP with fertility, a low heritability trait, were 50 k chip variants. For instance, rs110992367, the lead SNP for a 3.3 Mb QTL for the success rate of insemination in Holstein cows was 0.38 Mb from the 5' flanking region of the *G2E3* gene. Meta-analysis heterogeneity for this QTL was observed in all three breeds ($P = 9.5 \times 10^{-3}$), but opposite in the effect direction between Holstein and the other two breeds. The heterogeneity is important because the *G2E3* gene is an essential protein-coding gene that prevents apoptotic death during embryonic development (Brooks et al., 2008). This gene has been linked to apoptosis and cellular stress in US beef cattle (Howard et al., 2014) and somatic cell score in Chinese Holstein (Wang et al., 2015). This gene is a good candidate affecting cow fertility.

All QTL identified by meta-analysis were not significant in all within-breed analysis (Tables 1–6, Figure 3C) probably due to low power to reach significance at $P < 1 \times 10^{-6}$ or lack of segregation in the breed (Sham and Purcell, 2014).

For instance, the variant mapping to the *PLAG1* gene only segregated in the Montbeliarde breed. However, most of the SNP had an effect in the same direction in all breeds for the more significant variants, e.g., BTA18:58067310 (*CEACAM18* gene, fat) and BTA10:51558941 (*FAM63B* gene, stature). This suggests that most detected QTL segregate in all three breeds even though they were not significant in the within-breed analysis. In this study, pre-selected sequence variants were directly genotyped or imputed from a large reference population genotyped with the EuroG10k chip, and thus were known with very high accuracy. This was evident from the high average allelic correlation ($r^2 = 0.97$) between real and imputed genotypes for the three breeds. In conclusion, many QTL segregated across the three French cattle breeds more so for highly heritable traits. The meta-analysis across the breeds using the within-breed association summary data increased the power of QTL detection. This QTL information can be incorporated in routine genomic evaluation in French populations. If an independent population of dairy cattle is genotyped with the EuroG10k chip, then the true causal variant underlying the QTL may be validated in independent population/breed. The next step of this study is to use independent Nordic cattle data to validate the QTL reported here.

## AUTHOR CONTRIBUTIONS

AM analyzed data and wrote the manuscript. BG contributed to the design of the custom part of the EuroG10k chip, interpretation of results and revision of the manuscript. ML worked on the study design, secured project funding and revised the manuscript. SF contributed to the design of the custom part of the EuroG10k chip and the management of the marker information. GS interpreted the results and revised the manuscript. DB worked on the study design, interpretation of results and revision of the manuscript.

## FUNDING

AM is a recipient of an Erasmus Mundus Ph.D. grant funded jointly by European Graduate School in Animal Breeding and Genetics (France), and Center for Genomic Selection in Animals and Plants (GenSAP) supported by Innovation Fund Denmark (grant 0603-00519B).

## ACKNOWLEDGMENTS

We are grateful to French National Animal Breeding Database for providing phenotypic data, Valogene for access to genotypes and INRA for providing study environment. We also acknowledge Chris Hozé for initial genotype preparation.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2018.00522/full#supplementary-material
REFERENCES

Bartonova, P., Vrtkova, I., Kaplanova, K., and Urban, T. (2012). Association between CSN3 and BCO2 gene polymorphisms and milk performance traits in the Czech Fleckvieh cattle breed. *Genet. Mol. Res.* 11, 1058–1063. doi: 10.4238/2012.April.27.4

Boichard, D., Boussaha, M., Capitan, A., Rocha, D., Hozé C., Sanchez, M. P., et al. (2018). "Experience from large scale use of the EuroGenomics custom SNP chip in cattle," in *11th World Congress of Genetics Applied to Livestock Production* (Auckland), 675.

Boichard, D., Chung, H., Dassonneville, R., David, X., Eggen, A. A., Fritz, S. S., et al. (2012a). Design of a bovine low-density snp array optimized for imputation. *PLoS ONE* 7:e34130. doi: 10.1371/journal.pone.0034130

Boichard, D., Guillaume, F., Baur, A., Croiseau, P., Rossignol, M. N., Boscher, M. Y., et al. (2012b). Genomic selection in French dairy cattle. *Anim. Prod. Sci.* 52, 115–120. doi: 10.1071/AN11119

Bonfatti, V., Cecchinato, A., Di Martino, G., De Marchi, M., Gallo, L., and Carnier, P. (2011). Effect of κ-casein B relative content in bulk milk κ-casein on Montasio, Asiago, and Caciotta cheese yield using milk of similar protein composition. *J. Dairy Sci.* 94, 602–613. doi: 10.3168/jds.20110203368

Bouwman, A. C., Daetwyler, H. D., Chamberlain, A. J., Hurtado Ponce, C., Sargolzaei, M., Schenkel, F. S., et al. (2018). Meta-analysis of genome wide association studies for the stature of cattle reveals common genes that regulate size in mammals. *Nat. Gen.* 50, 362–367. doi: 10.1038/s41588-018-0056-5

Brooks, W. S., Helton, E. S., Banerjee, S., Venable, M., Johnson, L., Schoeb, T. R., et al. (2008). G2E3 is a dual function ubiquitin ligase required for early embryonic development. *J. Biol. Chem.* 283, 22304–22315. doi: 10.1074/jbc.m803238200

Chang, M., He, L., and Cai, L. (2018). An overview of genome wide association studies. *Comput. Syst. Biol.* 1754, 97–108. doi: 10.1007/978-1-4939-717-8_6

Church, P. C., Goscinski, A., and Lefèvre, C. (2012). EXP-PAC: Providing comparative analysis and storage of next generation gene expression data. *Genomics* 100, 8–13. doi: 10.1016/j.ygeno.2012.05.007

Cochran, W. G. (1954). The combination of estimates from different experiments. *Biometrics* 10, 101–129. doi: 10.2307/3001666

Daetwyler, H. D., Capitan, A., Pausch, H., Stothard, P., van Binsbergen, R., Brøndum, R. F., et al. (2014). Whole-genome sequencing of 234 bulls facilitates mapping of monogenic and complex traits in cattle. *Nat. Genet.* 46, 858–865. doi: 10.1038/ng.3034

Di Gregorio, P., Di Grigoli, A., Di Trana, A., Alabiso, M., Maniaci, G., Rando, A., et al. (2017). Effects of different genotypes at the CSN3 and LGB loci on milk and cheese-making characteristics of the bovine Cinisara breed. *Int. Dairy J.* 71:43105. doi: 10.1016/j.idairyj.2016.11.001

Fouillen, A., Dos Santos Neves, J., Mary, C., Castonguay, J.-D., Moffatt, P., Baron, C., et al. (2017). Interactions of AMTN, ODAM and SCPPPQ1 proteins of a specialized basal lamina that attaches epithelial cells to tooth mineral. *Sci. Rep.* 24:46683. doi: 10.1038/srep46683

Goddard, M. E., and Hayes, B. J. (2009). Mapping genes for complex traits in domestic animals and their use in breeding programmes. *Nat. Rev. Genet.* 10, 381–391. doi: 10.1038/nrg2575

Grisart, B., Coppieters, W., Farnir, F., Karim, L., Ford, C., Berzi, P., et al. (2002). Positional candidate cloning of a QTL in dairy cattle: identification of a missense mutation in the bovine DGAT1 gene with major effect on milk yield and composition. *Genome Res.* 12, 222–231. doi: 10.1101/gr.224202

Hoff, J. L., Decker, J. E., Schnabel, R. D., and Taylor, J. F. (2017). Candidate lethal haplotypes and causal mutations in Angus cattle. *BMC Genomics.* 18:799. doi: 10.1186/s12864-017-4196-2

Howard, J. T., Kachman, S. D., Snelling, W. M., Pollak, E. J., Ciobanu, D. C., Kuehn, L. A., et al. (2014). Beef cattle body temperature during climatic stress: a genome-wide association study. *Int. J. Biometeorol.* 58, 1665–1672. doi: 10.1007/s00484-013-0773-5

Hu, Z.-L., Park, C. A., Reecy, J. M. (2016). Developmental progress and current status of the Animal QTLdb. *Nucleic Acids Res.* 44, D827–D833. doi: 10.1093/nar/gkv1233

Leyva-Corona, J. C., Reyna-Granados, J. R., Zamorano-Algandar, R., et al. (2018). Polymorphisms within the prolactin and growth hormone/insulin-like growth factor-1 functional pathways associated with fertility traits in Holstein cows raised in a hot-humid climate. *Trop. Anim. Health Prod.* doi: 10.1007/s11250-018-1645-0. [Epub ahead of print].

Li, H., Wang, Z., Moore, S. S., Schenkel, F. S., and Stothard, P. (2010). *Genome-wide scan for positional and functional candidate genes affecting milk production traits in Canadian Holstein Cattle,* in *Proc 9Th WCGALP, Genetics of Trait Complexes: Lactation - Lecture Sessions: 0535* (Leipzig). Available online at: http://www.wcgalp.org

Marete, A., Lund, M. S., Boichard, D., and Ramayo-Caldas, Y. (2018a). A system-based analysis of the genetic determinism of udder conformation and health phenotypes across three French dairy cattle breeds. *PLoS ONE* 13:e0199931. doi: 10.1371/journal.pone.0199931

Marete, A., Sahana, G., Fritz, S., Lefebvre, R., Barbat, A., Lund, M. S., et al. (2018b). Genome-wide association study for milking speed in French Holstein cows. *J. Dairy Sci.* 101, 6205–6219. doi: 10.3168/jds.2017-14067

McLaren, W., Gil, L., Hunt, S. E., Riat, H. S., Ritchie, G. R. S., Thormann, A., et al. (2016). The Ensembl variant effect predictor. *Genome Biol.* 17:122. doi: 10.1186/s13059-016-0974-4

Michot, P., Fritz, S., Barbat, A., Boussaha, M., Deloche, M. C., Grohs, C., et al. (2017). A missense mutation in PFAS (phosphoribosylformylglycinamidine synthase) is likely causal for embryonic lethality associated with the MH1 haplotype in Montbéliarde dairy cattle. *J. Dairy Sci.* 100, 8176–8187. doi: 10.3168/jds.2017-12579

Pausch, H., Emmerling, R., Grendler-Grandl, B., Fries, R., Daetwyler, D., and Goddard, M. E. (2017). Meta-analysis of sequence-based association studies across three cattle breeds reveals 25 QTL for fat and protein percentages in milk at nucleotide resolution. *BMC Genomics* 18:853. doi: 10.1186/s12864-017-4263-8

Pausch, H., Emmerling, R., Schwarzenbacher, H., and Fries, R. (2016). A multi-trait meta-analysis with imputed sequence variants reveals twelve QTL for mammary gland morphology in Fleckvieh cattle. *Genet. Sel. Evol.* 48:14. doi: 10.1186/s12711-016-0190-4

Raven, L. A., Cocks, B. G., and Hayes, B. J. (2014). Multibreed genome wide association can improve precision of mapping causative variants underlying milk production in dairy cattle. *BMC Genomics.* 15:62. doi: 10.1186/1471-2164-15-62

Raven, L. A., Cocks, B. G., Kemper, K. E., Chamberlain, A. J., Jagt, C. J. V., Goddard, M. E., et al. (2016). Targeted imputation of sequence variants and gene expression profiling identifies twelve candidate genes associated with lactation volume, composition and calving interval in dairy cattle. *Mamm. Genome* 27, 81–97. doi: 10.1007/s00335-015-9613-8

Sargolzaei, M., Chesnais, J. P., and Schenkel, F. S. (2014). A new approach for efficient genotype imputation using information from relatives. *BMC Genomics* 15:478. doi: 10.1186/1471-2164-15-478

Sham, P. C., and Purcell, S. M. (2014). Statistical power and significance testing in large-scale genetic studies. *Nat. Rev. Genet.* 15, 335–346. doi: 10.1038/nr g3706

van den Berg, I., Boichard, D., and Lund, M. S. (2016). Comparing power and precision of within-breed and multibreed genome-wide association studies of production traits using whole-genome sequence data for 5 French and Danish dairy cattle breeds. *J. Dairy Sci.* 99, 8932–8945. doi: 10.3168/jds.2016-11073

VanRaden, P. M. (2008). Efficient methods to compute genomic predictions. *J. Dairy Sci.* 91, 4414–44232. doi: 10.3168/jds.2007-0980

VanRaden, P. M., and Wiggans, G. R. (1991). Derivation, calculation, and use of national animal model information. *J. Dairy Sci.* 74, 2737–2746. doi: 10.3168/jds.S0022-0302(91)78453-1

Wang, X., Ma, P., Liu, J., Zhang, Q., Zhang, Y., Ding, X., et al. (2015). Genome-wide association study in Chinese Holstein cows reveal two candidate genes for somatic cell score as an indicator for mastitis susceptibility. *BMC Genet.* 16:111. doi: 10.1186/s12863-015-0263-3

Wang, Y. H., Reverter, A., Mannen, H., Taniguchi, M., Harper, G. S., Oyama, K., et al. (2005). Transcriptional profiling of muscle tissue in growing Japanese
Black cattle to identify genes involved with the development of intramuscular fat. *Aust. J. Exp. Agric.* 45, 809–820. doi: 10.1071/EA05058

Willer, C. J., Li, Y., and Abecasis, G. R. (2010). METAL: fast and efficient meta-analysis of genomewide association scans. *Bioinformatics* 26, 2190–2191. doi: 10.1093/bioinformatics/btq340

Yang, J., Lee, S. H., Goddard, M. E., and Visscher, P. M. (2011). GCTA: a tool for genome-wide complex trait analysis. *Am. J. Hum. Genet.* 88, 76–82. doi: 10.1016/j.ajhg.2010.11.011

Zimin, A. V., Delcher, A. L., Florea, L., Kelley, D. R., Schatz, M. C., Puiu, D., et al. (2009). A whole-genome assembly of the domestic cow. *Genome Biol.* 10:R42. doi: 10.1186/gb-2009-10-4-r42

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2018 Marete, Guldbrandtsen, Lund, Fritz, Sahana and Boichard. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.