Genome-wide mapping of 10 calving and fertility traits in Holstein dairy cattle with special regard to chromosome 18.

Over the last decades, a dramatic decrease in reproductive performance has been observed in Holstein cattle and fertility problems have become the most common reason for a cow to leave the herd. The premature removal of animals with high breeding values results in both economic and breeding losses. For efficient future Holstein breeding, the identification of loci associated with low fertility is of major interest and thus constitutes the aim of this study. To reach this aim, a genome-wide combined linkage disequilibrium and linkage analysis (cLDLA) was conducted using data on the following 10 calving and fertility traits in the form of estimated breeding values: days from first service to conception of heifers and cows, nonreturn rate on d 56 of heifers and cows, days from calving to first insemination, days open, paternal and maternal calving ease, paternal and maternal stillbirth. The animal data set contained 2,527 daughter-proven Holstein bulls from Germany that were genotyped with Illumina's BovineSNP50 BeadChip (Illumina Inc., San Diego, CA). For the cLDLA, 41,635 sliding windows of 40 adjacent single nucleotide polymorphisms (SNP) were used. At each window midpoint, a variance component analysis was executed using ASReml. The underlying mixed linear model included random quantitative trait locus (QTL) and polygenic effects. We identified 50 genome-wide significant QTL. The most significant peak was detected for direct calving ease at 59,179,424 bp on chromosome 18 (BTA18). Next, a mixed-linear model association (MLMA) analysis was conducted. A comparison of the cLDLA and MLMA results with special regard to BTA18 showed that the genome-wide most significant SNP from the MLMA was associated with the same trait and located on the same chromosome at 57,589,121 bp (i.e., about 1.5 Mb apart from the cLDLA peak). The results of 5 different cLDLA and 2 MLMA models, which included the fixed effects of either SNP or haplotypes, suggested that the cLDLA method outperformed the MLMA in accuracy and precision. The haplotype-based cLDLA method allowed for a more precise mapping and the definition of ancestral and derived QTL alleles, both of which are essential for the detection of underlying quantitative trait nucleotides.

Genomic selection in the German Landrace population of the Bavarian herdbook.

The aim of our study was to compare different validation methods with respect to their impact on validation results and to evaluate the feasibility of genomic selection in the German Landrace population of the Bavarian herdbook. For this purpose, a sample of 337 boars and 1,676 sows was genotyped with the Illumina PorcineSNP60 BeadChip. Conventional BLUP breeding values for fertility, growth, carcass, and quality traits were deregressed and used as phenotypes in genomic BLUP. The resulting genomic breeding values were also blended with information from the full conventional breeding value estimation to include information from nongenotyped parents. Subsequent validation used forward prediction, realized reliabilities, and theoretical reliabilities. The results indicate that the validation methods showed a relatively large effect on in the displayed reliability levels in our study: forward prediction reliabilities were found to be much lower than the conventional parent-average reliabilities whereas corresponding realized and theoretical reliabilities were found substantially greater. Theoretical reliabilities appear to be the most consistent validation approach tested in our study, because they avoid the use of proxy variables. Generally, our results suggest a substantial potential for a genomic selection implementation for the Bavarian herdbook by using both sows and boars. Theoretical genomic reliabilities of direct genomic values of selection candidates were, on average, 31 to 36% greater than the conventional parent average reliabilities. However, the inclusion of residual information from conventional breeding values had only a marginal effect on reliabilities.

Genome-wide association mapping of milk production traits in Braunvieh cattle.

A whole-genome association study of milk production traits: milk yield, protein yield, fat yield, protein percentage, and fat percentage, was performed on the population of Braunvieh cattle. Five hundred and fifty-four progeny-tested bulls and 36,219 autosomal single nucleotide polymorphism (SNP) markers on 29 Bos taurus autosomes (BTA) were included in the analysis. A principal component analysis was conducted to adjust for the effect of population stratification in the analyzed data set. For the principal component analysis, genome-wide relationships between individuals were calculated. Three different criteria (Horn's test, Kaiser's criterion, and Jolliffe's criterion) were tested to determine the number of significant principal components. Estimation of putative associations between SNP and milk production traits was carried out using a linear regression model in R software (R Foundation for Statistical Computing, Vienna, Austria). Significant principal components, adjusting for population stratification separately for each criterion and family relationships and genotypes at individual SNP were included as fixed effects in the model. The inflation factor λ and quantile-quantile plots were calculated to compare how the different criteria deal with stratification in our mapping population. Based on the analyses on all of the aforementioned criteria, we can conclude that Jolliffe's criterion deals the best with population stratification. Furthermore, significance thresholds for a given genome-wide false discovery rate of 5% were estimated and used for specific traits. Three of the analyzed traits showed genome-wide significant association with SNP. Two SNP had an effect on milk yield on BTA4, 2 SNP affected fat yield on BTA14 and BTA23, and 1 SNP was associated with fat percent on BTA1. Single nucleotide polymorphisms identified in this study as associated with milk production traits will further contribute to the mapping of corresponding quantitative trait loci and investigation of the genes responsible for polymorphisms in milk production traits in dairy cattle. Described comparison of different criteria for determination of significant principal components can provide important information for similar studies in stratified populations.

SNP-based association mapping of the polled gene in divergent cattle breeds.

Naturally, hornless cattle are called polled. Although the POLL locus could be assigned to a c. 1.36-Mb interval in the centromeric region of BTA1, the underlying genetic basis for the polled trait is still unknown. Here, an association mapping design was set up to refine the candidate region of the polled trait for subsequent high-throughput sequencing. The case group comprised 101 homozygous polled animals from nine divergent cattle breeds, the majority represented by Galloway, Angus, Fleckvieh and Holstein Friesian. Additionally, this group included some polled individuals of Blonde d'Aquitaine, Charolais, Hereford, Jersey and Limousin breeds. The control group comprised horned Belgian Blue, Fleckvieh, Holstein Friesian and Illyrian Busa cattle. A genome-wide scan using 49,163 SNPs was performed, which revealed one shared homozygous haplotype block consisting of nine neighbouring SNPs in all polled animals. This segment defines a 381-kb interval on BTA1 that we consider to be the most likely location of the POLL mutation. Our results further demonstrate that the polled-associated haplotype is also frequent in horned animals included in this study, and thus the haplotype as such cannot be used for population-wide genetic testing. The actual trait-associated haplotype may be revealed by using higher-density SNP arrays. For the final identification of the causal mutation, we suggest high-throughput sequencing of the entire candidate region, because the identification of functional candidate genes is difficult owing to the lack of a comparable model.

SNP-based association mapping of Arachnomelia in Fleckvieh cattle.

Bovine arachnomelia is an inherited congenital disorder with malformation mainly of the limbs, the vertebral column and the skull, following a monogenic autosomal recessive heredity. Despite almost identical pathological findings, arachnomelia has previously been mapped to bovine chromosome 23 and 5 in Fleckvieh and Braunvieh respectively. Therefore, this disorder may be an example of locus heterogeneity in cattle. This study aimed to refine the candidate region to allow positional cloning and sequence analyses of candidate genes in Fleckvieh cattle. For that purpose, a case-control association mapping design was set up with a case group of 16 pre-selected affected individuals and a control group of 50 unrelated animals. The subset of affected animals was selected from a total of 129 pathologically confirmed cases due to the occurrence of recombination(s) within a 14.5 cM candidate interval previously mapped to chromosome 23. Six linked microsatellites currently used for indirect gene testing in Fleckvieh were analysed for this purpose. In all selected cases, a genome-wide scan using 44 473 informative SNPs revealed shared segments of homozygosity at 15 adjacent SNPs on chromosome 23. Additional haplotype analysis of 37 carrier bulls confirmed the localization of the arachnomelia locus to a region of 927 kb (13.622-14.549 Mb) containing molybdenum cofactor biosynthesis protein 1 gene, the most likely candidate gene for arachnomelia in Fleckvieh. The number of recombinant haplotypes observed in cases was more than doubled compared with the number of expected recombinations. This remarkably increased mapping resolution and thus illustrates the benefit of pre-selection in association studies.

Confirmation and refinement of a QTL on BTA5 affecting milk production traits in the Fleckvieh dual purpose cattle breed.

We analysed a QTL affecting milk yield (MY), milk protein yield (PY) and milk fat yield (FY) in the dual purpose cattle breed Fleckvieh on BTA5. Twenty-six microsatellite markers covering 135 cM were selected to analyse nine half-sib families containing 605 sons in a granddaughter design. We thereby assigned two new markers to the public linkage map using the CRI-MAP program. Phenotypic records were daughter yield deviations (DYD) originating from the routinely performed genetic evaluations of breeding animals. To determine the position of the QTL, three different approaches were applied: interval mapping (IM), linkage analysis by variance component analysis (LAVC), and combined linkage disequilibrium (LD) and linkage (LDL) analysis. All three methods mapped the QTL in the same marker interval (BM2830-ETH152) with the greatest test-statistic value at 118, 119.33 and 119.33 cM respectively. The positive QTL allele simultaneously increases DYD in the first lactation by 272 kg milk, 7.1 kg milk protein and 7.0 kg milk fat. Although the mapping accuracy and the significance of a QTL effect increased from IM over LAVC to LDL, the confidence interval was large (13, 20 and 24 cM for FY, MY and PY respectively) for the positional cloning of the causal gene. The estimated averages of pair wise marker LD with a distance <5 cM were low (0.107) and reflect the large effective population size of the Fleckvieh subpopulation analysed. This low level of LD suggests a need for increase in marker density in following fine mapping steps.

Mapping of QTL for Body Conformation and Behavior in Cattle.

Genome scans for quantitative trait loci (QTL) in farm animals have concentrated on primary production and health traits, and information on QTL for other important traits is rare. We performed a whole genome scan in a granddaughter design to detect QTL affecting body conformation and behavior in dairy cattle. The analysis included 16 paternal half-sib families of the Holstein breed with 872 sons and 264 genetic markers. The markers were distributed across all 29 autosomes and the pseudoautosomal region of the sex chromosomes with average intervals of 13.9 cM and covering an estimated 3155.5 cM. All families were analyzed jointly for 22 traits using multimarker regression and significance thresholds determined empirically by permutation. QTL that exceeded the experiment-wise significance threshold (5% level) were detected on chromosome 6 for foot angle, teat placement, and udder depth, and on chromosome 29 for temperament. QTL approaching experiment-wise significance (10% level) were located on chromosome 6 for general quality of feet and legs and general quality of udder, on chromosome 13 for teat length, on chromosome 23 for general quality of feet and legs, and on chromosome 29 for milking speed. An additional 51 QTL significant at the 5% chromosome-wise level were distributed over 21 chromosomes. This study provides the first evidence for QTL involved in behavior of dairy cattle and identifies QTL for udder conformation on chromosome 6 that could form the basis of recently reported QTL for clinical mastitis.

Quantitative trait loci mapping of functional traits in the German Holstein cattle population.

A whole-genome scan to detect quantitative trait loci (QTL) for functional traits was performed in the German Holstein cattle population. For this purpose, 263 genetic markers across all autosomes and the pseudoautosomal region of the sex chromosomes were genotyped in 16 granddaughter-design families with 872 sons. The traits investigated were deregressed breedingvalues for maternal and direct effects on dystocia (DYSm, DYSd) and stillbirth (STIm, STId) as well as maternal and paternal effects on nonreturn rates of 90 d (NR90m, NR90p). Furthermore, deregressed breeding values for functional herd life (FHL) and daughter yield deviation for somatic cell count (SCC) were investigated. Weighted multimarker regression analyses across families and permutation tests were applied for the detection of QTL and the calculation of statistical significance. A ten percent genomewise significant QTL was localized for DYSm on chromosome 8 and for SCC on chromosome 18. A further 24 putative QTL exceeding the 5% chromosomewise threshold were detected. On chromosomes 7, 8, 10, 18, and X/Yps, coincidence of QTL for several traits was observed. Our results suggest that loci with influence on udder health may also contribute to genetic variance of longevity. Prior to implementation of these QTL in marker assisted selection programs for functional traits, information about direct and correlated effects of these QTL as well as fine mapping of their chromosomal positions is required.

Mapping of the bovine blood group systems J, N', R', and Z show evidence for oligo-genetic inheritance.

Genes determining the bovine erythrocyte antigens were mapped by linkage analysis. In total 9591 genotypes of 20 grandsire families with 1074 sires from a grand-daughter design were elucidated for the genes determining the erythrocyte antigens EAA, EAB, EAC, EAF, EAJ, EAL, EAM, EAN', EAR', EAS, EAT', and EAZ according to standard paternity testing procedures in the blood typing laboratories. Linkage analyses were performed with 248 microsatellite markers, eight SSCP markers and four polymorphic proteins and enzymes covering the 29 autosomes and the pseudoautosomal region of the sex chromosomes. The number of informative meioses for the blood group systems ranged from 76 to 947. Blood group systems EAM and EAT' were non-informative. Most of the erythrocyte antigen loci showed significant linkage to a single chromosome and were mapped unequivocally. The genes determining erythrocyte antigen EAA, EAB, EAC, EAL, and EAS were mapped to chromosomes 15, 12, 18, 3, and 21, respectively. Lod-score values ranged from 11.43 to 107.83. Moreover, the EAF system could be mapped to chromosome 17. However, the EAN' system previously known as part of the EAF system could be mapped to chromosome 5. In addition, the blood group systems EAJ, the new EAN', EAR', and EAZ, showed significant linkage to microsatellite markers on various chromosomes and also to other blood groups. The appearance of a single blood group system might be therefore either dependent on the existence of other blood group systems or because of an interaction between different loci on various chromosomes as is known in humans and in pigs.

A whole genome scan for differences in recombination rates among three Bos taurus breeds.

Twenty paternal half-sib families of a granddaughter design were genotyped for 265 genetic markers, most of them microsatellites. These were 16 Holstein families, 3 Simmental families, and 1 Brown Swiss family. The number of sires per breed was 872, 170, and 32, respectively. Two-point recombination rates were estimated both jointly for all breeds and each single breed separately. Of 1168 marker intervals, 865 provided estimates for at least two breeds. Differences between breeds were tested by likelihood ratio tests. Four marker intervals, representing three genomic regions on BTA19, BTA24, and BTA27, show a significant impact of the breed at a false discovery rate of 0.23 and indicate a genetic component of observed heterogeneity of recombination. The variability of recombination rates between cattle breeds might not be a common feature of the whole genome, but rather might be restricted to certain chromosomal segments. Thus, attention should be paid to heterogeneities when pooling data of such regions from different breeds.

A QTL for the degree of spotting in cattle shows synteny with the KIT locus on chromosome 6.

The proportion of unpigmented coat on the trunk was determined from photographs of 38 German Simmental and 627 German Holstein bulls distributed over three generations. All 665 animals were members of 18 Holstein and 3 Simmental half-sib families. A Bayesian estimation of heritability yielded a posterior mean of 0.88 and a standard error of 0.08. A quantitative trait loci (QTL) scan over all chromosomes covered by 229 microsatellite marker loci (2926 cM) was performed by fitting a multiple marker regression model to 625 observations from the youngest generation in 18 families. On chromosome 6 a QTL for the proportion of white coat with large effects (experiment-wise error probability < .0001) was found and a less important one on chromosome 3 (chromosome-wise error probability < .009). Chromosome 6 is known to harbor the KIT locus (receptor tyrosinase kinase), which is associated with various depigmentation phenotypes in mice, humans, and pigs. Similarity of phenotypic KIT effects in other species and synteny with the reported QTL suggest that KIT is a serious candidate gene for the degree of spotting in cattle. The results are also discussed with respect to resistance to solar radiation, heat stress, and photosensitization.

A male linkage map of the cattle (Bos taurus) genome.

A male linkage map of the cattle (Bos taurus) genome was constructed using nine large half-sib families. The map consists of 269 loci, of which 249 are microsatellites and 20 are structural genes. Among the 249 microsatellites, 140 are markers selected from other maps and 98 are new assignments. Chromosome assignment were established for 35 new markers by somatic cell hybrid analysis, of which 26 were confirmed by linkage analysis. Genome coverage is 1975 cM contained within terminal markers on all 29 autosomes. The average distance between adjacent loci is 9.7 cM, with 72.1% of the map intervals < or = 15 cM and 4.9% of the intervals > or = 25 cM. The inclusion of mapped markers permitted integration and comparisons with other maps, facilitating the identification of discrepancies in chromosome assignment, gene order, and map distance. The inclusion of Type I and blood group markers in the map was useful for comparative mapping, revealing possible blood group orthologies between humans and cattle. The map generated will serve as a useful tool for comparative mapping, mapping of quantitative trait loci and marker assisted selection.

Isolation and characterization of 45 polymorphic microsatellites from the bovine genome.

A small-insert bovine genomic library was constructed in pBluescript II SK(+) and enriched for microsatellites by selective rescue of single-stranded pBluescript DNA carrying (CA)n/(TG)n tandem repeats. Approximately 50% of the clones in the enriched library contained (CA)n repeats or CA-rich sequences. Sequencing of clones selected for (CA)n repeats resulted in the identification and characterization of 45 (CA)n polymorphic microsatellites. Genotyping in 9 large paternal half-sib families indicated that 40 of these microsatellite markers exhibit autosomal Mendelian inheritance. The numbers of alleles range from 2 to 18, with an average of 6.3 per locus. The polymorphic microsatellite markers we have identified and characterized will contribute to the construction of a high-resolution linkage map of bovine genome.

Genetic mapping of F13A to BTA23 by sperm typing: difference in recombination rate between bulls in the DYA-PRL interval.

Our objective was to extend the linkage and comparative maps of BTA23 by determining whether the structural gene for the A subunit of blood coagulation factor XIII (F13A) is linked to BoLA-DYA, a centromeric marker, or the distally located gene encoding prolactin (PRL). Bovine F13A was mapped relative to DYA and PRL in an experiment that examined segregation of alleles in 176 sperm. Genotyping was performed by PCR-RFLP for all loci, following amplification of the haploid genome by primer extension preamplification. F13A was found to be linked to PRL (theta = 0.314 +/- 0.038). The most likely order is DYA-PRL-F13A (odds > 10(4):1). This result demonstrates conservation of synteny between BTA23 and most of HSA6p. Surprisingly, theta DYA-PRL was 0.310 +/- 0.039, 83.4% greater (P < 0.02) than we found for another bull (Van Eijk et al., Mamm. Genome 4: 113, 1993). The difference in recombination rate in the DYA-PRL interval provides further evidence for an unusual recombination hot spot between the bovine Mhc class IIa and class IIb subregions and suggests that bull-specific maps may be necessary for marker-assisted selection.

Marker-derived phylogeny of European cattle supports demic expansion of agriculture.

Using frequencies of 86 genes from 23 loci of blood group systems, blood and milk proteins, the genetic relationships among 14 cattle breeds including four native Balkan and four synthetic Balkan-Alp breeds were studied. The dendrogram and nonlinear map construction shows a consensus 'Balkan breed cluster', an 'Alp breed cluster', an unstable position of synthetic breeds and well-separated American breeds. Positive partial correlations between genetic distance and time elapsed since introduction of farming while keeping geographical distances constant, and regular patterns over thousands of kilometers indicate that large-scale cattle population movements together with human migration (in the Neolithic age) from the Near East into Europe across the Balkans are the most likely explanation for the genetic distances observed in our data. More recent breed differentiation and selection do not yet blur this initial pattern of European cattle populations.

Variable expression of epitopes on the surface of Mycoplasma gallisepticum demonstrated with monoclonal antibodies.

Twelve monoclonal antibodies (Mabs) against Mycoplasma gallisepticum (Mg) strains F, R, S6(208) and PET2 were used for analysis of epitopes of 22 Mg strains. Six Mabs recognized surface epitopes in the majority of strains, but did not react with variant strains like K 503 and K 703. Two Mabs reacted with epitopes on about 56 kilodalton (kDa) proteins and showing consistent expression on Mg colonies. Three Mabs recognized three different variable surface epitopes associated with about 67 kDa proteins and one Mab variable epitope on about 33 and 80 kDa proteins. Two-dimensional immunoblotting showed considerable differences in the charge of proteins bearing variable surface epitopes in different Mg strains. Subcloning of four low passage Mg strains using Mabs for screening populations that derived from a single colony with defined surface epitopes showed that some colonies may switch surface epitopes associated with 67 and 80 kDa proteins. This switching was reversible and generated subpopulations of Mg expressing different combinations of surface epitopes. Phenotypic switching of epitopes probably occurs also in vivo and may be the mechanism enabling Mg to evade the host immune response.

Order of bovine DRB3, DYA, and PRL determined by sperm typing.

The order and recombination fractions (theta) between the bovine major histocompatibility complex DRB3, DYA, and prolactin (PRL) genes were determined by typing of 254 sperm from a triply heterozygous bull. A recently developed method, primer extension preamplification (PEP), was used to amplify the bovine sperm genome prior to amplification of specific loci by the polymerase chain reaction (PCR). At least 28 copies of the DRB3, PRL, or DYA gene were obtained from 50 cycles of PEP. For sperm typing, alleles of each locus were discriminated by restriction endonuclease cleavage of PCR products and polyacrylamide gel electrophoresis of the restriction fragments. The most likely gene order is PRL-DRB3-DYA, with theta = 0.025 (+/- 0.012) and theta = 0.150 (+/- 0.024), respectively. The odds are 128:1 in favor of this order in comparison with the second most likely order DRB3-PRL-DYA. Our results demonstrate the power of sperm typing in concert with PEP for multilocus gene mapping.