A first-generation metric linkage disequilibrium map of bovine chromosome 6

Genetic structure analysis in several populations of cattle using SNP genotypes

Parameters such as effective population size (Ne), runs of homozygosity (ROH), and inbreeding based on ROH (FROH) can give new insight into the level of genetic diversity for the population under selection. This research aimed to measure the extent of linkage disequilibrium (LD), effective population size (Ne), Haplotype Block Structure, and runs of homozygosity (ROHs) in several populations of cattle using SNP genotypes. In this study, that the average r² decreased with the increasing distance of SNP pairs. A general decrease in Ne can be seen for all four populations, indicating a loss of genetic diversity. The Iranian Holstein had the lowest level of genomic inbreeding at an ROH of 1, 5, 10 Mb, while the French Holstein had the highest. The maximum number of ROH is seen at a distance of less than 1 Mb, and the lowest number of ROH is seen at a distance of 10 Mb. The number of ROH decreases with increasing distance due to the increased recombination rate. This is a concern as an increase in inbreeding leads to a reduction in the effective population size, which was also evident in the study populations.

Comparative accuracies of genetic values predicted for economically important milk traits, genome-wide association, and linkage disequilibrium patterns of Canadian Holstein cows

Genomic selection methodologies and genome-wide association studies use powerful statistical procedures that correlate large amounts of high-density SNP genotypes and phenotypic data. Actual 305-d milk (MY), fat (FY), and protein (PY) yield data on 695 cows and 76,355 genotyping-by-sequencing-generated SNP marker genotypes from Canadian Holstein dairy cows were used to characterize linkage disequilibrium (LD) structure of Canadian Holstein cows. Also, the comparison of pedigree-based BLUP, genomic BLUP (GBLUP), and Bayesian (BayesB) statistical methods in the genomic selection methodologies and the comparison of Bayesian ridge regression and BayesB statistical methods in the genome-wide association studies were carried out for MY, FY, and PY. Results from LD analysis revealed that as marker distance decreases, LD increases through chromosomes. However, unexpected high peaks in LD were observed between marker pairs with larger marker distances on all chromosomes. The GBLUP and BayesB models resulted in similar heritability estimates through 10-fold cross-validation for MY and PY; however, the GBLUP model resulted in higher heritability estimates than BayesB model for FY. The predictive ability of GBLUP model was significantly lower than that of BayesB for MY, FY, and PY. Association analyses indicated that 28 high-effect markers and markers on Bos taurus autosome 14 located within 6 genes (DOP1B, TONSL, CPSF1, ADCK5, PARP10, and GRINA) associated significantly with FY.

The breeding structure for the small ruminant resources in India

Intense selection for a few desired traits has resulted in reduction of the effective population size (Ne) in most of the plant and livestock populations across the world. The objective of the research was to assess the impact of Ne on the genetic architecture of the population in a simulated data with variable Ne for general population under selection. Along with this, the estimate of Ne and its ratio to adult breeding population (N_B) in the census data of small ruminants of India were also investigated. Results indicated that the average inbreeding ([Formula: see text]) decreases with increase in Ne; similarly, increase in [Formula: see text] per generation was highest in population with lowest Ne. Correlation of estimated breeding value (EBV) with true breeding value (TBV) was not much affected with effective population size. An effective number of chromosome segments (Me) in the populations under selection were significantly affected by magnitude of Ne, with linear positive relation between Ne and Me. Results on livestock census data revealed that all the sheep and goat breeds have sufficiently large Ne based on derived and actual census data. The median for ratio of effective population size to adult census size in sheep breeds was 0.120 and for goat breeds was 0.131. Karnah and Poonchi sheep shares the status of endangered breeds due to a smaller number of breeding female population and hence need attention for conservation. The Ne was large in sheep and goat due to less selection pressure as a result of low coverage of breed improvement programs, availability of large number of breeding males, and absence of artificial insemination (AI) in the field flocks. The estimates of Ne and its ratio to the adult census size (N_B) excluded several factors such as fluctuating population size and overlapping generations. Study revealed introspection from most of the industrial breeding programs on the issue of Ne for populations under selection. Similarly, in small ruminants, large Ne indicates huge genetic diversity and scope of improvement in the productivity in near future.

The utility of low-density genotyping for imputation in the Thoroughbred horse

Background

Despite the dramatic reduction in the cost of high-density genotyping that has occurred over the last decade, it remains one of the limiting factors for obtaining the large datasets required for genomic studies of disease in the horse. In this study, we investigated the potential for low-density genotyping and subsequent imputation to address this problem.

Results

Using the haplotype phasing and imputation program, BEAGLE, it is possible to impute genotypes from low- to high-density (50K) in the Thoroughbred horse with reasonable to high accuracy. Analysis of the sources of variation in imputation accuracy revealed dependence both on the minor allele frequency of the single nucleotide polymorphisms (SNPs) being imputed and on the underlying linkage disequilibrium structure. Whereas equidistant spacing of the SNPs on the low-density panel worked well, optimising SNP selection to increase their minor allele frequency was advantageous, even when the panel was subsequently used in a population of different geographical origin. Replacing base pair position with linkage disequilibrium map distance reduced the variation in imputation accuracy across SNPs. Whereas a 1K SNP panel was generally sufficient to ensure that more than 80% of genotypes were correctly imputed, other studies suggest that a 2K to 3K panel is more efficient to minimize the subsequent loss of accuracy in genomic prediction analyses. The relationship between accuracy and genotyping costs for the different low-density panels, suggests that a 2K SNP panel would represent good value for money.

Conclusions

Low-density genotyping with a 2K SNP panel followed by imputation provides a compromise between cost and accuracy that could promote more widespread genotyping, and hence the use of genomic information in horses. In addition to offering a low cost alternative to high-density genotyping, imputation provides a means to combine datasets from different genotyping platforms, which is becoming necessary since researchers are starting to use the recently developed equine 70K SNP chip. However, more work is needed to evaluate the impact of between-breed differences on imputation accuracy.

Maximum likelihood estimation of linkage disequilibrium in half-sib families

Maximum likelihood methods for the estimation of linkage disequilibrium between biallelic DNA-markers in half-sib families (half-sib method) are developed for single and multifamily situations. Monte Carlo computer simulations were carried out for a variety of scenarios regarding sire genotypes, linkage disequilibrium, recombination fraction, family size, and number of families. A double heterozygote sire was simulated with recombination fraction of 0.00, linkage disequilibrium among dams of δ=0.10, and alleles at both markers segregating at intermediate frequencies for a family size of 500. The average estimates of δ were 0.17, 0.25, and 0.10 for Excoffier and Slatkin (1995), maternal informative haplotypes, and the half-sib method, respectively. A multifamily EM algorithm was tested at intermediate frequencies by computer simulation. The range of the absolute difference between estimated and simulated δ was between 0.000 and 0.008. A cattle half-sib family was genotyped with the Illumina 50K BeadChip. There were 314,730 SNP pairs for which the sire was a homo-heterozygote with average estimates of r2 of 0.115, 0.067, and 0.111 for half-sib, Excoffier and Slatkin (1995), and maternal informative haplotypes methods, respectively. There were 208,872 SNP pairs for which the sire was double heterozygote with average estimates of r2 across the genome of 0.100, 0.267, and 0.925 for half-sib, Excoffier and Slatkin (1995), and maternal informative haplotypes methods, respectively. Genome analyses for all possible sire genotypes with 829,042 tests showed that ignoring half-sib family structure leads to upward biased estimates of linkage disequilibrium. Published inferences on population structure and evolution of cattle should be revisited after accommodating existing half-sib family structure in the estimation of linkage disequilibrium.

Functional characterisation of the bovine neuropeptide Y gene promoter and evaluation of the transcriptional activities of promoter haplotypes

Neuropeptide Y (NPY) is a potent orexigenic agent. The molecular mechanisms underlying the regulation of bovine NPY gene expression by its promoter region is currently unknown. The objectives of this research were to: (i) identify the SNPs in the promoter region of the bovine NPY gene, (ii) investigate the effects of these SNPs by measuring promoter transcriptional activities of different bovine NPY promoter haplotypes and; (iii) identify the minimal promoter region (MPR) required for basal activity of the NPY gene in vitro. Seventeen SNPs were identified in the promoter region. Of these, 14 affected putative transcription factors binding motifs including a TATA binding protein factor at -20, GC-Box factors SP1 at -170 and GATA binding motifs at -120 and -347. The SNPs were assigned to five major haplotypes (BtNPY_H1-5), of which BtNPY_H5 had maximum transcriptional activity. The region extending to -134 nt was identified as the MPR. This MPR was confirmed by the identification of a putative TATA box (-29 nt) and two SP1/GC binding sites (-94 and -118 nt), within this region. However, promoter expression was significantly enhanced when the construct contained the -614 to -1019 nt region. In conclusion, a number of SNPs characterised in the bovine NPY promoter especially those affecting the transcription factor binding sites, enhancer and repressor regions have the potential to affect NPY gene expression. Natural variation exists in the promoter region of the bovine NPY gene, which should be further explored for selection of energetic efficiency in cattle.

The pattern of linkage disequilibrium in German Holstein cattle

This study presents a second generation of linkage disequilibrium (LD) map statistics for the whole genome of the Holstein-Friesian population, which has a four times higher resolution compared with that of the maps available so far. We used DNA samples of 810 German Holstein-Friesian cattle genotyped by the Illumina Bovine SNP50K BeadChip to analyse LD structure. A panel of 40 854 (75.6%) markers was included in the final analysis. The pairwise r(2) statistic of SNPs up to 5 Mb apart across the genome was estimated. A mean value of r(2) = 0.30 +/- 0.32 was observed in pairwise distances of <25 kb and it dropped to 0.20 +/- 0.24 at 50-75 kb, which is nearly the average inter-marker space in this study. The proportion of SNPs in useful LD (r(2) > or = 0.25) was 26% for the distance of 50 and 75 kb between SNPs. We found a lower level of LD for SNP pairs at the distance < or =100 kb than previously thought. Analysis revealed 712 haplo-blocks spanning 4.7% of the genome and containing 8.0% of all SNPs. Mean and median block length were estimated as 164 +/- 117 kb and 144 kb respectively. Allele frequencies of the SNPs have a considerable and systematic impact on the estimate of r(2). It is shown that minimizing the allele frequency difference between SNPs reduces the influence of frequency on r(2) estimates. Analysis of past effective population size based on the direct estimates of recombination rates from SNP data showed a decline in effective population size to N(e) = 103 up to approximately 4 generations ago. Systematic effects of marker density and effective population size on observed LD and haplotype structure are discussed.

Extensive linkage disequilibrium in a wild bird population

Knowledge about the extent and patterns of linkage disequilibrium (LD) can provide important insights into demographic processes and strategies to identify the genetic basis of complex phenotypes in wild populations. However, data on the extent and patterns of LD from non-model vertebrate species from the wild are still scarce. We conducted so far the most extensive and detailed examination of LD in a pedigreed wild bird population using genotypes from 97 autosomal and 6 gonosomal microsatellites and a recently established linkage map of Siberian jays (Perisoreus infaustus). Analysis of syntenic marker pairs showed high levels of LD that extended over tens of centimorgans or several megabases and generally decayed as an increasing function of intermarker distance. In addition, significant LD was also very common between nonsyntenic markers. Patterns of LD varied across different linkage groups possibly because of the differences in chromosomal structure (macro-, micro-, and Z-chromosome). In particular, the level of LD was significantly lower on the Z-chromosome than on the autosomes at comparable genetic distances. In general, the high levels and extent of LD in this population are likely owing to its relatively small size, significant intrapopulation genetic structure, and occurrence of inbreeding. Whatever the cause, the long-range LD between syntenic loci suggests that LD mapping of phenotypic traits in this population using low-density markers maps is feasible. However, the frequent occurrence of LD between nonsyntenic markers suggests that the combined use of linkage and LD methods is needed to reduce the likelihood of false-positive associations between marker loci and traits of ecological and evolutionary interest.

A genome wide survey of SNP variation reveals the genetic structure of sheep breeds

The genetic structure of sheep reflects their domestication and subsequent formation into discrete breeds. Understanding genetic structure is essential for achieving genetic improvement through genome-wide association studies, genomic selection and the dissection of quantitative traits. After identifying the first genome-wide set of SNP for sheep, we report on levels of genetic variability both within and between a diverse sample of ovine populations. Then, using cluster analysis and the partitioning of genetic variation, we demonstrate sheep are characterised by weak phylogeographic structure, overlapping genetic similarity and generally low differentiation which is consistent with their short evolutionary history. The degree of population substructure was, however, sufficient to cluster individuals based on geographic origin and known breed history. Specifically, African and Asian populations clustered separately from breeds of European origin sampled from Australia, New Zealand, Europe and North America. Furthermore, we demonstrate the presence of stratification within some, but not all, ovine breeds. The results emphasize that careful documentation of genetic structure will be an essential prerequisite when mapping the genetic basis of complex traits. Furthermore, the identification of a subset of SNP able to assign individuals into broad groupings demonstrates even a small panel of markers may be suitable for applications such as traceability.

Linkage disequilibrium in the North American Holstein population

Linkage disequilibrium was estimated using 7119 single nucleotide polymorphism markers across the genome and 200 animals from the North American Holstein cattle population. The analysis of maternally inherited haplotypes revealed strong linkage disequilibrium (r(2) > 0.8) in genomic regions of approximately 50 kb or less. While linkage disequilibrium decays as a function of genomic distance, genomic regions within genes showed greater linkage disequilibrium and greater variation in linkage disequilibrium compared with intergenic regions. Identification of haplotype blocks could characterize the most common haplotypes. Although maximum haplotype block size was over 1 Mb, mean block size was 26-113 kb by various definitions, which was larger than that observed in humans ( approximately 10 kb). Effective population size of the dairy cattle population was estimated from linkage disequilibrium between single nucleotide polymorphism marker pairs in various haplotype ranges. Rapid reduction of effective population size of dairy cattle was inferred from linkage disequilibrium in recent generations. This result implies a loss of genetic diversity because of the high rate of inbreeding and high selection intensity in dairy cattle. The pattern observed in this study indicated linkage disequilibrium in the current dairy cattle population could be exploited to refine mapping resolution. Changes in effective population size during past generations imply a necessity of plans to maintain polymorphism in the Holstein population.

A second generation genetic map for rainbow trout (Oncorhynchus mykiss)

Background

Genetic maps characterizing the inheritance patterns of traits and markers have been developed for a wide range of species and used to study questions in biomedicine, agriculture, ecology and evolutionary biology. The status of rainbow trout genetic maps has progressed significantly over the last decade due to interest in this species in aquaculture and sport fisheries, and as a model research organism for studies related to carcinogenesis, toxicology, comparative immunology, disease ecology, physiology and nutrition. We constructed a second generation genetic map for rainbow trout using microsatellite markers to facilitate the identification of quantitative trait loci for traits affecting aquaculture production efficiency and the extraction of comparative information from the genome sequences of model fish species.

Results

A genetic map ordering 1124 microsatellite loci spanning a sex-averaged distance of 2927.10 cM (Kosambi) and having 2.6 cM resolution was constructed by genotyping 10 parents and 150 offspring from the National Center for Cool and Cold Water Aquaculture (NCCCWA) reference family mapping panel. Microsatellite markers, representing pairs of loci resulting from an evolutionarily recent whole genome duplication event, identified 180 duplicated regions within the rainbow trout genome. Microsatellites associated with genes through expressed sequence tags or bacterial artificial chromosomes produced comparative assignments with tetraodon, zebrafish, fugu, and medaka resulting in assignments of homology for 199 loci.

Conclusion

The second generation NCCCWA genetic map provides an increased microsatellite marker density and quantifies differences in recombination rate between the sexes in outbred populations. It has the potential to integrate with cytogenetic and other physical maps, identifying paralogous regions of the rainbow trout genome arising from the evolutionarily recent genome duplication event, and anchoring a comparative map with the zebrafish, medaka, tetraodon, and fugu genomes. This resource will facilitate the identification of genes affecting traits of interest through fine mapping and positional cloning of candidate genes.

The concordance test emerges as a powerful tool for identifying quantitative trait nucleotides: lessons from BTA6 milk yield QTL

The lack of conventions for confirming the discovery of quantitative trait nucleotides in livestock was evidenced by the proposals of mutations in two different genes (SPP1 and ABCG2) as the underlying functional mutation for a major quantitative trait locus (QTL) for milk concentration on bovine chromosome 6 (BTA6). Of these conflicting candidates, SPP1 was excluded by follow-up studies and by the data described here. A simple test for concordance of the zygosity state between QTL segregation status and the candidate polymorphism was shown, in this case, to be a critical step towards establishing the proof. If a given sample effectively represents the genetic variation across the QTL region, haplotype-based concordance may further enhance the functionality and resolution power of this test, allowing identification of the causative gene.

Linkage disequilibrium compared between five populations of domestic sheep

BACKGROUND

The success of genome-wide scans depends on the strength and magnitude of linkage disequilibrium (LD) present within the populations under investigation. High density SNP arrays are currently in development for the sheep genome, however little is known about the behaviour of LD in this livestock species. This study examined the behaviour of LD within five sheep populations using two LD metrics, D' and x2'. Four economically important Australian sheep flocks, three pure breeds (White Faced Suffolk, Poll Dorset, Merino) and a crossbred population (Merino x Border Leicester), along with an inbred Australian Merino museum flock were analysed.

RESULTS

Short range LD (0 - 5 cM) was observed in all five populations, however the persistence with increasing distance and magnitude of LD varied considerably between populations. Average LD (x2') for markers spaced up to 20 cM exceeded the non-syntenic average within the White Faced Suffolk, Poll Dorset and Macarthur Merino. LD decayed faster within the Merino and Merino x Border Leicester, with LD below or consistent with observed background levels. Using marker-marker LD as a guide to the behaviour of marker-QTL LD, estimates of minimum marker spacing were made. For a 95% probability of detecting QTL, a microsatellite marker would be required every 0.1 - 2.5 centimorgans, depending on the population used.

CONCLUSION

Sheep populations were selected which were inbred (Macarthur Merino), highly heterogeneous (Merino) or intermediate between these two extremes. This facilitated analysis and comparison of LD (x2') between populations. The strength and magnitude of LD was found to differ markedly between breeds and aligned closely with both observed levels of genetic diversity and expectations based on breed history. This confirmed that breed specific information is likely to be important for genome wide selection and during the design of successful genome scans where tens of thousands of markers will be required.

Linkage disequilibrium and signatures of selection on chromosomes 19 and 29 in beef and dairy cattle

The objective of this study was to quantify the extent of linkage disequilibrium (LD) on bovine chromosomes 19 and 29 and to study the pattern of selection signatures in beef and dairy breeds (Angus and Holstein) of Bos taurus. The extent of LD was estimated for 370 and 186 single nucleotide polymorphism markers on BTA19 and 29 respectively using the square of the correlation coefficient (r(2)) among alleles at pairs of loci. A comparison of the extent of LD found that the decline of LD followed a similar pattern in both breeds. We observed long-range LD and found that LD dissipates to background levels at a locus separation of about 20 Mb on both chromosomes. Along each chromosome, patterns of LD were variable in both breeds. We find that a minimum of 30 000 informative and evenly spaced markers would be required for whole-genome association studies in cattle. In addition, we have identified chromosomal regions that show some evidence of selection for economically important traits in Angus and Holstein cattle. The results of this study are of importance for the design and application of association studies.

Linkage disequilibrium and persistence of phase in Holstein-Friesian, Jersey and Angus cattle

When a genetic marker and a quantitative trait locus (QTL) are in linkage disequilibrium (LD) in one population, they may not be in LD in another population or their LD phase may be reversed. The objectives of this study were to compare the extent of LD and the persistence of LD phase across multiple cattle populations. LD measures r and r(2) were calculated for syntenic marker pairs using genomewide single-nucleotide polymorphisms (SNP) that were genotyped in Dutch and Australian Holstein-Friesian (HF) bulls, Australian Angus cattle, and New Zealand Friesian and Jersey cows. Average r(2) was approximately 0.35, 0.25, 0.22, 0.14, and 0.06 at marker distances 10, 20, 40, 100, and 1000 kb, respectively, which indicates that genomic selection within cattle breeds with r(2) >or= 0.20 between adjacent markers would require approximately 50,000 SNPs. The correlation of r values between populations for the same marker pairs was close to 1 for pairs of very close markers (<10 kb) and decreased with increasing marker distance and the extent of divergence between the populations. To find markers that are in LD with QTL across diverged breeds, such as HF, Jersey, and Angus, would require approximately 300,000 markers.

Extent of linkage disequilibrium in Holstein cattle in North America

The success of fine-scale mapping and genomic selection depends mainly on the strength of linkage disequilibrium (LD) between markers and causal mutations. With Lewontin's measure of LD (known as D'), high levels of LD that extend over several million base pairs have been reported in livestock. However, this measure of LD can be strongly biased upward by small samples and by low allele frequencies. The aim of this study was to characterize the level and extent of LD in Holstein cattle in North America (Canada and the United States for purposes of this study) by using the squared correlation of the alleles at 2 loci (r(2)). The Affymetrix MegAllele GeneChip Bovine Mapping 10K single nucleotide polymorphism (SNP) array was used to genotype 821 bulls, from which 497 were used in the analysis of the extent of LD. A total of 5,564 SNP were used after filtering out SNP with more than 5% of Mendelian inconsistencies, with more than 20% missing genotypes, or with a minor allele frequency of less than 10%. Analysis of syntenic pairs revealed that useful LD (measured as r(2) > 0.3) occurred at distances shorter than 100 kb. Linkage disequilibrium decayed very rapidly, within a few hundred kilobase pairs. In addition, no substantial LD between unlinked loci was found. Using a sliding window analysis, we observed an irregular pattern of LD across the genome. These findings suggest that to capture useful LD, which is required for whole-genome fine mapping and genomic selection, a denser SNP map would be needed.

Extent of genome-wide linkage disequilibrium in Australian Holstein-Friesian cattle based on a high-density SNP panel

BACKGROUND

The extent of linkage disequilibrium (LD) within a population determines the number of markers that will be required for successful association mapping and marker-assisted selection. Most studies on LD in cattle reported to date are based on microsatellite markers or small numbers of single nucleotide polymorphisms (SNPs) covering one or only a few chromosomes. This is the first comprehensive study on the extent of LD in cattle by analyzing data on 1,546 Holstein-Friesian bulls genotyped for 15,036 SNP markers covering all regions of all autosomes. Furthermore, most studies in cattle have used relatively small sample sizes and, consequently, may have had biased estimates of measures commonly used to describe LD. We examine minimum sample sizes required to estimate LD without bias and loss in accuracy. Finally, relatively little information is available on comparative LD structures including other mammalian species such as human and mouse, and we compare LD structure in cattle with public-domain data from both human and mouse.

RESULTS

We computed three LD estimates, D', Dvol and r2, for 1,566,890 syntenic SNP pairs and a sample of 365,400 non-syntenic pairs. Mean D' is 0.189 among syntenic SNPs, and 0.105 among non-syntenic SNPs; mean r2 is 0.024 among syntenic SNPs and 0.0032 among non-syntenic SNPs. All three measures of LD for syntenic pairs decline with distance; the decline is much steeper for r2 than for D' and Dvol. The value of D' and Dvol are quite similar. Significant LD in cattle extends to 40 kb (when estimated as r2) and 8.2 Mb (when estimated as D'). The mean values for LD at large physical distances are close to those for non-syntenic SNPs. Minor allelic frequency threshold affects the distribution and extent of LD. For unbiased and accurate estimates of LD across marker intervals spanning < 1 kb to > 50 Mb, minimum sample sizes of 400 (for D') and 75 (for r2) are required. The bias due to small samples sizes increases with inter-marker interval. LD in cattle is much less extensive than in a mouse population created from crossing inbred lines, and more extensive than in humans.

CONCLUSION

For association mapping in Holstein-Friesian cattle, for a given design, at least one SNP is required for each 40 kb, giving a total requirement of at least 75,000 SNPs for a low power whole-genome scan (median r2 > 0.19) and up to 300,000 markers at 10 kb intervals for a high power genome scan (median r2 > 0.62). For estimation of LD by D' and Dvol with sufficient precision, a sample size of at least 400 is required, whereas for r2 a minimum sample of 75 is adequate.

Associations between 2 paternal casein haplotypes and milk yield traits of Swiss Fleckvieh cattle

Associations between casein haplotypes and milk yield traits of offspring from 5 Swiss Fleckvieh AI test bulls were investigated. The analysis was performed by using a daughter design, where each daughter inherited either paternal haplotype B-A1-A-A or B-A2-A-A for alleles of alpha s1-, beta-, alpha s2- and kappa-casein genes. The substitution effects of paternal CSN2 A1 versus A2 on protein yield deviations (YDs) were significant (P < 0.05), whereas their effects on milk and fat YDs were not. The paternal substitution effects of the CSN2 A1 versus the A2 allele on protein YDs within the 5 sires did not reach the significance level. This is due to the contrary allele substitution effect of a sire compared to the other 4 sires. The effects of maternal haplotypes on milk, protein and fat YDs were not significant. However, it is noteworthy that the effects of haplotypes with a low frequency in the population deviate largely from the most frequent haplotype B-A2-A-A. The effects of beta-lactoglobulin (BLG) genotypes were significant for protein YDs but not for milk and fat YDs. The association between the paternal CSN2 A1 and A2 alleles and milk protein YDs within sires but not milk and fat YDs indicate an interaction, which might be a consequence of CSN2 heterogeneity or a closely linked gene that is contributing to the estimated effects.

The origin of selection signatures on bovine chromosome 6

The extent and pattern of linkage disequilibrium (LD) between closely spaced markers contain information about population history, including past population size and selection history. Selection signatures can be identified by comparing the LD surrounding a putative selected allele at a locus to the putative non-selected allele. In livestock populations, locations of selection signatures identified in this way should be correlated with QTL affecting production traits, as the populations have been under strong artificial selection for these traits. We used a dense SNP map of bovine chromosome 6 to characterize the pattern of LD on this chromosome in Norwegian Red cattle, a breed which has been strongly selected for milk production. The pattern of LD was generally consistent with strong selection in regions containing QTL affecting milk production traits, including a strong selection signature in a region containing a mutation known to affect milk production. The results demonstrate that in livestock populations, the origin of selection signatures will often be QTL for livestock production traits, and illustrate the value of selection signatures in uncovering new mutations with potential effects on quantitative traits.

Enhancing genetic mapping of complex genomes through the design of highly-multiplexed SNP arrays: application to the large and unsequenced genomes of white spruce and black spruce

Background

To explore the potential value of high-throughput genotyping assays in the analysis of large and complex genomes, we designed two highly multiplexed Illumina bead arrays using the GoldenGate SNP assay for gene mapping in white spruce (Picea glauca [Moench] Voss) and black spruce (Picea mariana [Mill.] B.S.P.).

Results

Each array included 768 SNPs, identified by resequencing genomic DNA from parents of each mapping population. For white spruce and black spruce, respectively, 69.2% and 77.1% of genotyped SNPs had valid GoldenGate assay scores and segregated in the mapping populations. For each of these successful SNPs, on average, valid genotyping scores were obtained for over 99% of progeny. SNP data were integrated to pre-existing ALFP, ESTP, and SSR markers to construct two individual linkage maps and a composite map for white spruce and black spruce genomes. The white spruce composite map contained 821 markers including 348 gene loci. Also, 835 markers including 328 gene loci were positioned on the black spruce composite map. In total, 215 anchor markers (mostly gene markers) were shared between the two species. Considering lineage divergence at least 10 Myr ago between the two spruces, interspecific comparison of homoeologous linkage groups revealed remarkable synteny and marker colinearity.

Conclusion

The design of customized highly multiplexed Illumina SNP arrays appears as an efficient procedure to enhance the mapping of expressed genes and make linkage maps more informative and powerful in such species with poorly known genomes. This genotyping approach will open new avenues for co-localizing candidate genes and QTLs, partial genome sequencing, and comparative mapping across conifers.

Genomewide analysis of secretory activation in mouse models

Mouse models have been widely used to elucidate the biology of mammary gland development and secretory activation. Recent advances in the availability of genomic resources for mice will generate a renewed effort to define the genetic basis of lactation phenotypes and help identify candidate gene pathways. Specific aspects of these advances are relevant to the dairy industry and may provide a rationale for improving milk production in the dairy cow. Differences are evident in mammary gland morphology and various characteristics of milk production of inbred mouse strains, but few studies have undertaken any systematic phenotypic analysis of the different inbred strains of mice for lactation performance. Whole genome association analysis using recent strain-specific genotype data and detailed phenotype measurements from available inbred strains, along with transcript profiling of divergent inbred strains for lactation performance, provides a valuable approach to identify putative candidate genes and associated pathways underlying dairy QTL intervals. Here we discuss the utility of integrating mouse phenomic and genomic resources for understanding secretory activation in the mammary gland.

Whole genome linkage disequilibrium maps in cattle

Background

Bovine whole genome linkage disequilibrium maps were constructed for eight breeds of cattle. These data provide fundamental information concerning bovine genome organization which will allow the design of studies to associate genetic variation with economically important traits and also provides background information concerning the extent of long range linkage disequilibrium in cattle.

Results

Linkage disequilibrium was assessed using r² among all pairs of syntenic markers within eight breeds of cattle from the Bos taurus and Bos indicus subspecies. Bos taurus breeds included Angus, Charolais, Dutch Black and White Dairy, Holstein, Japanese Black and Limousin while Bos indicus breeds included Brahman and Nelore. Approximately 2670 markers spanning the entire bovine autosomal genome were used to estimate pairwise r² values. We found that the extent of linkage disequilibrium is no more than 0.5 Mb in these eight breeds of cattle.

Conclusion

Linkage disequilibrium in cattle has previously been reported to extend several tens of centimorgans. Our results, based on a much larger sample of marker loci and across eight breeds of cattle indicate that in cattle linkage disequilibrium persists over much more limited distances. Our findings suggest that 30,000–50,000 loci will be needed to conduct whole genome association studies in cattle.

Genetic and haplotypic structure in 14 European and African cattle breeds

To evaluate and compare the extent of LD in cattle, 1536 SNPs, mostly localized on BTA03, were detected in silico from available sequence data using two different methods and genotyped on samples from 14 distinct breeds originating from Europe and Africa. Only 696 SNPs could be validated, confirming the importance of trace-quality information for the in silico detection. Most of the validated SNPs were informative in several breeds and were used for a detailed description of their genetic structure and relationships. Results obtained were in agreement with previous studies performed on microsatellite markers and using larger samples. In addition, the majority of the validated SNPs could be mapped precisely, reaching an average density of one marker every 311 kb. This allowed us to analyze the extent of LD in the different breeds. Decrease of LD with physical distance across breeds revealed footprints of ancestral LD at short distances (<10 kb). As suggested by the haplotype block structure, these ancestral blocks are organized, within a breed, into larger blocks of a few hundred kilobases. In practice, such a structure similar to that already reported in dogs makes it possible to develop a chip of <300,000 SNPs, which should be efficient for mapping purposes in most cattle breeds.

A validated whole-genome association study of efficient food conversion in cattle

The genetic factors that contribute to efficient food conversion are largely unknown. Several physiological systems are likely to be important, including basal metabolic rate, the generation of ATP, the regulation of growth and development, and the homeostatic control of body mass. Using whole-genome association, we found that DNA variants in or near proteins contributing to the background use of energy of the cell were 10 times as common as those affecting appetite and body-mass homeostasis. In addition, there was a genic contribution from the extracellular matrix and tissue structure, suggesting a trade-off between efficiency and tissue construction. Nevertheless, the largest group consisted of those involved in gene regulation or control of the phenotype. We found that the distribution of micro-RNA motifs was significantly different for the genetic variants associated with residual feed intake than for the genetic variants in total, although the distribution of promoter sequence motifs was not different. This suggests that certain subsets of micro-RNA are more important for the regulation of this trait. Successful validation depended on the sign of the allelic association in different populations rather than on the strength of the initial association or its size of effect.

A primary assembly of a bovine haplotype block map based on a 15,036-single-nucleotide polymorphism panel genotyped in holstein-friesian cattle

Analysis of data on 1000 Holstein-Friesian bulls genotyped for 15,036 single-nucleotide polymorphisms (SNPs) has enabled genomewide identification of haplotype blocks and tag SNPs. A final subset of 9195 SNPs in Hardy-Weinberg equilibrium and mapped on autosomes on the bovine sequence assembly (release Btau 3.1) was used in this study. The average intermarker spacing was 251.8 kb. The average minor allele frequency (MAF) was 0.29 (0.05-0.5). Following recent precedents in human HapMap studies, a haplotype block was defined where 95% of combinations of SNPs within a region are in very high linkage disequilibrium. A total of 727 haplotype blocks consisting of > or =3 SNPs were identified. The average block length was 69.7 +/- 7.7 kb, which is approximately 5-10 times larger than in humans. These blocks comprised a total of 2964 SNPs and covered 50,638 kb of the sequence map, which constitutes 2.18% of the length of all autosomes. A set of tag SNPs, which will be useful for further fine-mapping studies, has been identified. Overall, the results suggest that as many as 75,000-100,000 tag SNPs would be needed to track all important haplotype blocks in the bovine genome. This would require approximately 250,000 SNPs in the discovery phase.

Genome-wide genetic diversity of Holstein Friesian cattle reveals new insights into Australian and global population variability, including impact of selection

Past breeding strategies for dairy cattle have been very effective in producing rapid genetic gain to achieve industry targets and raise profitability. Such gains have been largely facilitated by intense selection of sires combined with the use of artificial insemination. However, this practice can potentially limit the level of genetic diversity through inbreeding and selection plateaus. The rate of inbreeding in Australia is increasing, primarily as a result of semen importation from a small number of prominent bulls from the USA. The effect of this genetic influx in the Australian dairy cattle population is poorly understood both in terms of diversity and local adaptation/divergence. This study uses 845 genome-wide SNP genetic markers and 431 bulls to characterize the level of genetic diversity and genetic divergence within the Australian and international Holstein Friesian dairy population. No significant differences in genetic diversity (as measured by heterozygosity [H(o)] and allelic richness [A]) were observed over the 25-year time period (1975-1999) for bulls used in Australia. The importation of foreign semen into Australia has increased the effective population size until it was in effect a sub-sample of the global population. Our data indicate that most individuals are equally closely related to one another, regardless of country of origin and year of birth. In effect, the global population can be considered as one single population unit. These results indicate that inbreeding, genetic drift and selection has had little effect at reducing genetic diversity and differentiating the Australian Holstein Friesian population at a genome-wide level.

Exploiting large scale computing to construct high resolution linkage disequilibrium maps of the human genome

Linkage disequilibrium (LD) maps increase power and precision in association mapping, define optimal marker spacing and identify recombination hot-spots and regions influenced by natural selection. Phase II of HapMap provides approximately 2.8-fold more single nucleotide polymorphisms (SNPs) than phase I for constructing higher resolution maps. LDMAP-cluster, is a parallel program for rapid map construction in a Linux environment used here to construct genome-wide LD maps with >8.2 million SNPs from the phase II data.

AVAILABILITY

The LD maps, LDMAP-cluster and documentation are available from: http://www.som.soton.ac.uk/research/geneticsdiv/epidemiology/LDMAP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.