Linkage mapping bovine EST-based SNP

Biomarker Discovery in Animal Health and Disease: The Application of Post-Genomic Technologies

The causes of many important diseases in animals are complex and multifactorial, which present unique challenges. Biomarkers indicate the presence or extent of a biological process, which is directly linked to the clinical manifestations and outcome of a particular disease. Identifying biomarkers or biomarker profiles will be an important step towards disease characterization and management of disease in animals. The emergence of post-genomic technologies has led to the development of strategies aimed at identifying specific and sensitive biomarkers from the thousands of molecules present in a tissue or biological fluid. This review will summarize the current developments in biomarker discovery and will focus on the role of transcriptomics, proteomics and metabolomics in biomarker discovery for animal health and disease.

An integrated approach to gene discovery and marker development in Atlantic cod (Gadus morhua)

Atlantic cod is a species that has been overexploited by the capture fishery. Programs to domesticate this species are underway in several countries, including Canada, to provide an alternative route for production. Selective breeding programs have been successfully applied in the domestication of other species, with genomics-based approaches used to augment conventional methods of animal production in recent years. Genomics tools, such as gene sequences and sets of variable markers, also have the potential to enhance and accelerate selective breeding programs in aquaculture, and to provide better monitoring tools to ensure that wild cod populations are well managed. We describe the generation of significant genomics resources for Atlantic cod through an integrated genomics/selective breeding approach. These include 158,877 expressed sequence tags (ESTs), a set of annotated putative transcripts and several thousand single nucleotide polymorphism markers that were developed from, and have been shown to be highly variable in, fish enrolled in two selective breeding programs. Our EST collection was generated from various tissues and life cycle stages. In some cases, tissues from which libraries were generated were isolated from fish exposed to stressors, including elevated temperature, or antigen stimulation (bacterial and viral) to enrich for transcripts that are involved in these response pathways. The genomics resources described here support the developing aquaculture industry, enabling the application of molecular markers within selective breeding programs. Marker sets should also find widespread application in fisheries management.

Status of the cattle genome map

During the last 30 years, the cattle genome map has been expanded from 4 genes linked on chromosome X to over 22,000 genes identified in the cattle genome sequence assembly. This progress has been achieved due to numerous projects on linkage and physical mapping of the cattle genome driven by its agricultural and scientific significance. Indeed, the high-resolution mapping and functional analysis of the genome led to the discovery of major quantitative trait loci (QTL) regions and several quantitative trait nucleotides (QTNs), as well as some disease genes in the cow population. In addition, a comparison of the cattle genome to the genomes of other mammals has revealed its unique features gained during the speciation and adaptation. With the development of non-expensive sequencing techniques, the analysis of the cattle genome will shift towards the identification of differences between breeds or individuals within breeds that account for the unique features of each breed. This approach holds promise for the development of effective tools for the marker assistant selection and disease diagnostics in cattle.

Precision of distances and ordering of microsatellite markers in consensus linkage maps of chromosomes 1, 3 and 4 from two reciprocal chicken populations using bootstrap sampling

Some factors complicate comparisons between linkage maps from different studies. This problem can be resolved if measures of precision, such as confidence intervals and frequency distributions, are associated with markers. We examined the precision of distances and ordering of microsatellite markers in the consensus linkage maps of chromosomes 1, 3 and 4 from two F(2) reciprocal Brazilian chicken populations, using bootstrap sampling. Single and consensus maps were constructed. The consensus map was compared with the International Consensus Linkage Map and with the whole genome sequence. Some loci showed segregation distortion and missing data, but this did not affect the analyses negatively. Several inversions and position shifts were detected, based on 95% confidence intervals and frequency distributions of loci. Some discrepancies in distances between loci and in ordering were due to chance, whereas others could be attributed to other effects, including reciprocal crosses, sampling error of the founder animals from the two populations, F(2) population structure, number of and distance between microsatellite markers, number of informative meioses, loci segregation patterns, and sex. In the Brazilian consensus GGA1, locus LEI1038 was in a position closer to the true genome sequence than in the International Consensus Map, whereas for GGA3 and GGA4, no such differences were found. Extending these analyses to the remaining chromosomes should facilitate comparisons and the integration of several available genetic maps, allowing meta-analyses for map construction and quantitative trait loci (QTL) mapping. The precision of the estimates of QTL positions and their effects would be increased with such information.

Assignment of chromosomal locations for unassigned SNPs/scaffolds based on pair-wise linkage disequilibrium estimates

Background

Recent developments of high-density SNP chips across a number of species require accurate genetic maps. Despite rapid advances in genome sequence assembly and availability of a number of tools for creating genetic maps, the exact genome location for a number of SNPs from these SNP chips still remains unknown. We have developed a locus ordering procedure based on linkage disequilibrium (LODE) which provides estimation of the chromosomal positions of unaligned SNPs and scaffolds. It also provides an alternative means for verification of genetic maps. We exemplified LODE in cattle.

Results

The utility of the LODE procedure was demonstrated using data from 1,943 bulls genotyped for 73,569 SNPs across three different SNP chips. First, the utility of the procedure was tested by analysing the masked positions of 1,500 randomly-chosen SNPs with known locations (50 from each chromosome), representing three classes of minor allele frequencies (MAF), namely >0.05, 0.01<MAF ≤ 0.05 and 0.001<MAF ≤ 0.01. The efficiency (percentage of masked SNPs that could be assigned a location) was 96.7%, 30.6% and 2.0%; with an accuracy (the percentage of SNPs assigned correctly) of 99.9%, 98.9% and 33.3% in the three classes of MAF, respectively. The average precision for placement of the SNPs was 914, 3,137 and 6,853 kb, respectively. Secondly, 4,688 of 5,314 SNPs unpositioned in the Btau4.0 assembly were positioned using the LODE procedure. Based on these results, the positions of 485 unordered scaffolds were determined. The procedure was also used to validate the genome positions of 53,068 SNPs placed on Btau4.0 bovine assembly, resulting in identification of problem areas in the assembly. Finally, the accuracy of the LODE procedure was independently validated by comparative mapping on the hg18 human assembly.

Conclusion

The LODE procedure described in this study is an efficient and accurate method for positioning SNPs (MAF>0.05), for validating and checking the quality of a genome assembly, and offers a means for positioning of unordered scaffolds containing SNPs. The LODE procedure will be helpful in refining genome sequence assemblies, especially those being created from next-generation sequencing where high-throughput SNP discovery and genotyping platforms are integrated components of genome analysis.

A gene-based SNP linkage map for pacific white shrimp, Litopenaeus vannamei

Pacific white shrimp (Litopenaeus vannamei) are of particular economic importance to the global shrimp aquaculture industry. However, limited genomics information is available for the penaeid species. We utilized the limited public information available, mainly single nucleotide polymorphisms (SNPs) and expressed sequence tags, to discover markers for the construction of the first SNP genetic map for Pacific white shrimp. In total, 1344 putative SNPs were discovered, and out of 825 SNPs genotyped, 418 SNP markers from 347 contigs were mapped onto 45 sex-averaged linkage groups, with approximate coverage of 2071 and 2130 cm for the female and male maps, respectively. The average-squared correlation coefficient (r(2)), a measure of linkage disequilibrium, for markers located more than 50 cm apart on the same linkage group, was 0.15. Levels of r(2) increased with decreasing inter-marker distance from approximately 80 cm, and increased more rapidly from approximately 30 cm. A QTL for shrimp gender was mapped on linkage group 13. Comparative mapping to model organisms, Daphnia pulex and Drosophila melanogaster, revealed extensive rearrangement of genome architecture for L. vannamei, and that L. vannamei was more related to Daphnia pulex. This SNP genetic map lays the foundation for future shrimp genomics studies, especially the identification of genetic markers or regions for economically important traits.

Refined positioning of a quantitative trait locus affecting somatic cell score on chromosome 18 in the German Holstein using linkage disequilibrium

Combined linkage and linkage disequilibrium analysis (LALD) was conducted to more accurately map a previously reported quantitative trait locus (QTL) affecting somatic cell score on bovine chromosome 18. A grand-daughter design consisting of 6 German Holstein grandsire families with 1,054 progeny-tested genotyped sons was used in this study. Twenty microsatellite markers, 5 single nucleotide polymorphisms, and an erythrocyte antigen marker with an average marker spacing of 1.95 cM were analyzed along a chromosomal segment of 50.80 cM. Variance components were estimated and restricted maximum likelihood test statistics were calculated at the midpoint of each marker interval. The test statistics calculated in single-QTL linkage analysis exceeded the genome-wide significance threshold at several putative QTL positions. Using LALD, we were successful in assigning a genome-wide significant QTL to a confidence interval of 10.8 cM between the markers ILSTS002 and BMS833. The QTL in this marker interval was estimated to be responsible for between 5.89 and 13.86% of the genetic variation in somatic cell score. In contrast to the single-QTL linkage analysis model, LALD analyses with a 2-QTL model confirmed the position of one QTL, but gave no conclusive evidence for the existence or position of a second QTL. Ultimately, the QTL position was narrowed down considerably compared with previous results with a refined confidence interval of less than 11 cM.

A whole genome scan to map QTL for milk production traits and somatic cell score in Canadian Holstein bulls

The detection and mapping of genetic markers linked to quantitative trait loci (QTL) can be utilized to enhance genetic improvement of livestock populations. With the completion of the bovine genome sequence assembly, single nucleotide polymorphisms (SNP) assays spanning the whole bovine genome and research work on large scale identification, validation and analysis of genotypic variation in cattle has become possible. The objective of the present study was to perform a whole genome scan to identify and map QTL affecting milk production traits and somatic cell scores using linkage disequilibrium (LD) regression and 1536 SNP markers. Three and 18 SNP were found to be associated with only milk yield (MY) at a genome and chromosome wise significance (p < 0.05) level respectively. Among the 21 significant SNP, 16 were in a region reported to have QTL for MY in other dairy cattle populations and while the rest five were new QTL finding. Four SNP out of 21 are significant for the milk production traits (MY, fat yield, protein yield (PY), and milk contents) in the present study. Six and nine SNP were associated with PY at a genome and chromosome wise significant (p < 0.05) level respectively. Three and 17 SNP were found to be associated with FY at a genome and chromosome wise significant (p < 0.05) level. Five and seven SNP were mapped with somatic cell score at a genome and chromosome wise significant (p < 0.05) level respectively. The results of this study have revealed QTL for MY, PY, protein percentage, FY, fat percentage, somatic cell score and persistency of milk in the Canadian dairy cattle population. The chromosome regions identified in this study should be further investigated to potentially identify the causative mutations underlying the QTL.

Identification of a two-marker-haplotype on Bos taurus autosome 18 associated with somatic cell score in German Holstein cattle

Background

The somatic cell score (SCS) is implemented in routine sire evaluations in many countries as an indicator trait for udder health. Somatic cell score is highly correlated with clinical mastitis, and in the German Holstein population quantitative trait loci (QTL) for SCS have been repeatedly mapped on Bos taurus autosome 18 (BTA18). In the present study, we report a refined analysis of previously detected QTL regions on BTA18 with the aim of identifying marker and marker haplotypes in linkage disequilibrium with SCS. A combined linkage and linkage disequilibrium approach was implemented, and association analyses of marker genotypes and maternally inherited two-marker-haplotypes were conducted to identify marker and haplotypes in linkage disequilibrium with a locus affecting SCS in the German Holstein population.

Results

We detected a genome-wide significant QTL within marker interval 9 (HAMP_c.366+109G>A - BMS833) in the middle to telomeric region on BTA18 and a second putative QTL in marker interval 12-13 (BB710 - PVRL2_c.392G>A). Association analyses with genotypes of markers flanking the most likely QTL positions revealed the microsatellite marker BMS833 (interval 9) to be associated with a locus affecting SCS within the families investigated. A further analysis of maternally inherited two-marker haplotypes and effects of maternally inherited two-marker-interval gametes indicated haplotype 249-G in marker interval 12-13 (BB710 - PVRL2_c.392G>A) to be associated with SCS in the German Holstein population.

Conclusion

Our results confirmed previous QTL mapping results for SCS and support the hypothesis that more than one locus presumably affects udder health in the middle to telomeric region of BTA18. However, a subsequent investigation of the reported QTL regions is necessary to verify the two-QTL hypothesis and confirm the association of two-marker-haplotype 249-G in marker interval 12-13 (BB710 - PVRL2_c.392G>A) with SCS. For this purpose, higher marker density and multiple-trait and multiple-QTL models are required to narrow down the position of the causal mutation or mutations affecting SCS in German Holstein cattle.

A high density linkage map of the bovine genome

BACKGROUND

Recent technological advances have made it possible to efficiently genotype large numbers of single nucleotide polymorphisms (SNPs) in livestock species, allowing the production of high-density linkage maps. Such maps can be used for quality control of other SNPs and for fine mapping of quantitative trait loci (QTL) via linkage disequilibrium (LD).

RESULTS

A high-density bovine linkage map was constructed using three types of markers. The genotypic information was obtained from 294 microsatellites, three milk protein haplotypes and 6769 SNPs. The map was constructed by combining genetic (linkage) and physical information in an iterative mapping process. Markers were mapped to 3,155 unique positions; the 6,924 autosomal markers were mapped to 3,078 unique positions and the 123 non-pseudoautosomal and 19 pseudoautosomal sex chromosome markers were mapped to 62 and 15 unique positions, respectively. The linkage map had a total length of 3,249 cM. For the autosomes the average genetic distance between adjacent markers was 0.449 cM, the genetic distance between unique map positions was 1.01 cM and the average genetic distance (cM) per Mb was 1.25.

CONCLUSION

There is a high concordance between the order of the SNPs in our linkage map and their physical positions on the most recent bovine genome sequence assembly (Btau 4.0). The linkage maps provide support for fine mapping projects and LD studies in bovine populations. Additionally, the linkage map may help to resolve positions of unassigned portions of the bovine genome.

Construction of genetic linkage maps and comparative genome analysis of catfish using gene-associated markers

A genetic linkage map of the channel catfish genome (N=29) was constructed using EST-based microsatellite and single nucleotide polymorphism (SNP) markers in an interspecific reference family. A total of 413 microsatellites and 125 SNP markers were polymorphic in the reference family. Linkage analysis using JoinMap 4.0 allowed mapping of 331 markers (259 microsatellites and 72 SNPs) to 29 linkage groups. Each linkage group contained 3-18 markers. The largest linkage group contained 18 markers and spanned 131.2 cM, while the smallest linkage group contained 14 markers and spanned only 7.9 cM. The linkage map covered a genetic distance of 1811 cM with an average marker interval of 6.0 cM. Sex-specific maps were also constructed; the recombination rate for females was 1.6 times higher than that for males. Putative conserved syntenies between catfish and zebrafish, medaka, and Tetraodon were established, but the overall levels of genome rearrangements were high among the teleost genomes. This study represents a first-generation linkage map constructed by using EST-derived microsatellites and SNPs, laying a framework for large-scale comparative genome analysis in catfish. The conserved syntenies identified here between the catfish and the three model fish species should facilitate structural genome analysis and evolutionary studies, but more importantly should facilitate functional inference of catfish genes. Given that determination of gene functions is difficult in nonmodel species such as catfish, functional genome analysis will have to rely heavily on the establishment of orthologies from model species.

Confirmation of quantitative trait loci using a low-density single nucleotide polymorphism map for twinning and ovulation rate on bovine chromosome 5

Traditional genetic selection in cattle for traits with low heritability, such as reproduction, has had very little success. With the addition of DNA technologies to the genetic selection toolbox for livestock, the opportunity may exist to improve reproductive efficiency more rapidly in cattle. The US Meat Animal Research Center Production Efficiency Population has 9,186 twinning and 29,571 ovulation rate records for multiple generations of animals, but a significant number of these animals do not have tissue samples available for DNA genotyping. The objectives of this study were to confirm QTL for twinning and ovulation rate previously found on BTA5 and to evaluate the ability of GenoProb to predict genotypic information in a pedigree containing 16,035 animals when using genotypes for 24 SNP from 3 data sets containing 48, 724, or 2,900 animals. Marker data for 21 microsatellites on BTA5 with 297 to 3,395 animals per marker were used in conjunction with each data set of genotyped animals. Genotypic probabilities for females were used to calculate independent variables for regressions of additive, dominance, and imprinting effects. Genotypic regressions were fitted as fixed effects in a 2-trait mixed model analysis by using multiple-trait derivative-free REML. Each SNP was analyzed individually, followed by backward selection fitting all individually significant SNP simultaneously and then removing the least significant SNP until only significant SNP were left. Five significant SNP associations were detected for twinning rate and 3 were detected for ovulation rate. Two of these SNP, 1 for each trait, were significant for imprinting. Additional modeling of paternal and maternal allelic effects confirmed the initial results of imprinting done by contrasting heterozygotes. These results are supported by comparative mapping of mouse and human imprinted genes to this region of bovine chromosome 5.

A whole-genome scan to map quantitative trait loci for conformation and functional traits in Canadian Holstein bulls

Genetic improvement of livestock populations can be achieved through detection and mapping of genetic markers linked to quantitative trait loci (QTL). With the completion of the bovine genome sequence assembly, single nucleotide polymorphism (SNP) assays spanning the whole bovine genome and research work on large-scale identification, validation, and analysis of genotypic variation in cattle has become possible. A total of 462 Canadian Holstein Bulls were used to test the association between SNP and QTL. Single locus linkage disequilibrium regression model was implemented to perform a whole genome scan to identify and map QTL affecting conformation and functional traits. One thousand five hundred thirty-six SNP markers from introns and exons of potential QTL regions for economically important traits across the bovine genome were selected for association analysis. A total of 45 and 151 SNP were found to be associated with 17 conformation and functional traits at a genome- and chromosome-wise significance level, respectively. Among the 196 significant SNP, 169 of them are newly detected in this study, whereas 27 of them have been reported in previous literature and 161 of these were located in genes and are worth further investigating to potentially identify the causative mutations underlying the QTL. The single locus linkage disequilibrium regression method using SNP marker genotypes has proven to be a successful methodology for detecting and mapping QTL in dairy cattle populations.

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.

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.

A physical map of the bovine genome

A new physical map of the bovine genome has been constructed by integrating data from genetic and radiation hybrid maps, and a new bovine BAC map, with the bovine genome draft assembly.

Background

Cattle are important agriculturally and relevant as a model organism. Previously described genetic and radiation hybrid (RH) maps of the bovine genome have been used to identify genomic regions and genes affecting specific traits. Application of these maps to identify influential genetic polymorphisms will be enhanced by integration with each other and with bacterial artificial chromosome (BAC) libraries. The BAC libraries and clone maps are essential for the hybrid clone-by-clone/whole-genome shotgun sequencing approach taken by the bovine genome sequencing project.

Results

A bovine BAC map was constructed with HindIII restriction digest fragments of 290,797 BAC clones from animals of three different breeds. Comparative mapping of 422,522 BAC end sequences assisted with BAC map ordering and assembly. Genotypes and pedigree from two genetic maps and marker scores from three whole-genome RH panels were consolidated on a 17,254-marker composite map. Sequence similarity allowed integrating the BAC and composite maps with the bovine draft assembly (Btau3.1), establishing a comprehensive resource describing the bovine genome. Agreement between the marker and BAC maps and the draft assembly is high, although discrepancies exist. The composite and BAC maps are more similar than either is to the draft assembly.

Conclusion

Further refinement of the maps and greater integration into the genome assembly process may contribute to a high quality assembly. The maps provide resources to associate phenotypic variation with underlying genomic variation, and are crucial resources for understanding the biology underpinning this important ruminant species so closely associated with humans.

Use of linked loci as individuals or haplotypes for marker-assisted breed assignment

The objective of this study was to use simulation to evaluate the benefits of considering haplotypes of loci when linked single nucleotide polymorphisms are used for breed assignment. Three breeds of 10,000 females each were simulated under eight scenarios that differed according to the number of generations separating the breeds, size of breed founder populations and recombination rate between linked loci. Molecular genotypes consisted of 20 groups of three linked loci each. Breed assignment was performed in the final generation and was based on the frequency method. Haplotypes were reconstructed using the expectation-maximization algorithm. Accuracy of breed assignment was based on the frequency of correct breed assignment. Assignment accuracy increased as more genotypes (loci or haplotypes) were considered and more animals were used to estimate genotypic frequencies within breed. For most scenarios, use of haplotypes yielded equal or greater accuracies than when loci were considered independent. The advantage of haplotypes tended to increase as linkage disequilibrium between adjacent loci increased. The greatest advantage for using haplotypes was observed when recombination rate was low (0.001), breeds were separated by few generations (100), and a relatively large number of founder animals (110) was used to form new breeds. In this situation, 90% accuracy of breed assignment was achieved using nine to 14 haplotypes (i.e. 27-42 loci) depending on breed, vs. 39-57 individual loci.

High resolution radiation hybrid maps of bovine chromosomes 19 and 29: comparison with the bovine genome sequence assembly

Background

High resolution radiation hybrid (RH) maps can facilitate genome sequence assembly by correctly ordering genes and genetic markers along chromosomes. The objective of the present study was to generate high resolution RH maps of bovine chromosomes 19 (BTA19) and 29 (BTA29), and compare them with the current 7.1X bovine genome sequence assembly (bovine build 3.1). We have chosen BTA19 and 29 as candidate chromosomes for mapping, since many Quantitative Trait Loci (QTL) for the traits of carcass merit and residual feed intake have been identified on these chromosomes.

Results

We have constructed high resolution maps of BTA19 and BTA29 consisting of 555 and 253 Single Nucleotide Polymorphism (SNP) markers respectively using a 12,000 rad whole genome RH panel. With these markers, the RH map of BTA19 and BTA29 extended to 4591.4 cR and 2884.1 cR in length respectively. When aligned with the current bovine build 3.1, the order of markers on the RH map for BTA19 and 29 showed inconsistencies with respect to the genome assembly. Maps of both the chromosomes show that there is a significant internal rearrangement of the markers involving displacement, inversion and flips within the scaffolds with some scaffolds being misplaced in the genome assembly. We also constructed cattle-human comparative maps of these chromosomes which showed an overall agreement with the comparative maps published previously. However, minor discrepancies in the orientation of few homologous synteny blocks were observed.

Conclusion

The high resolution maps of BTA19 (average 1 locus/139 kb) and BTA29 (average 1 locus/208 kb) presented in this study suggest that by the incorporation of RH mapping information, the current bovine genome sequence assembly can be significantly improved. Furthermore, these maps can serve as a potential resource for fine mapping QTL and identification of causative mutations underlying QTL for economically important traits.

Primary genome scan to identify putative quantitative trait loci for feedlot growth rate, feed intake, and feed efficiency of beef cattle

Feed intake and feed efficiency of beef cattle are economically relevant traits. The study was conducted to identify QTL for feed intake and feed efficiency of beef cattle by using genotype information from 100 microsatellite markers and 355 SNP genotyped across 400 progeny of 20 Angus, Charolais, or Alberta Hybrid bulls. Traits analyzed include feedlot ADG, daily DMI, feed-to-gain ratio [F:G, which is the reciprocal of the efficiency of gain (G:F)], and residual feed intake (RFI). A mixed model with sire as random and QTL effects as fixed was used to generate an F-statistic profile across and within families for each trait along each chromosome, followed by empirical permutation tests to determine significance thresholds for QTL detection. Putative QTL for ADG (chromosome-wise P < 0.05) were detected across families on chromosomes 5 (130 cM), 6 (42 cM), 7 (84 cM), 11 (20 cM), 14 (74 cM), 16 (22 cM), 17 (9 cM), 18 (46 cM), 19 (53 cM), and 28 (23 cM). For DMI, putative QTL that exceeded the chromosome-wise P < 0.05 threshold were detected on chromosomes 1 (93 cM), 3 (123 cM), 15 (31 cM), 17 (81 cM), 18 (49 cM), 20 (56 cM), and 26 (69 cM) in the across-family analyses. Putative across-family QTL influencing F:G that exceeded the chromosome-wise P < 0.05 threshold were detected on chromosomes 3 (62 cM), 5 (129 cM), 7 (27 cM), 11 (16 cM), 16 (30 cM), 17 (81 cM), 22 (72 cM), 24 (55 cM), and 28 (24 cM). Putative QTL influencing RFI that exceeded the chromosome-wise P < 0.05 threshold were detected on chromosomes 1 (90 cM), 5 (129 cM), 7 (22 cM), 8 (80 cM), 12 (89 cM), 16 (41 cM), 17 (19 cM), and 26 (48 cM) in the across-family analyses. In addition, a total of 4, 6, 1, and 8 chromosomes showed suggestive evidence (chromosome-wise, P < 0.10) for putative ADG, DMI, F:G, and RFI QTL, respectively. Most of the QTL detected across families were also detected within families, although the locations across families were not necessarily the locations within families, which is likely because of differences among families in marker informativeness for the different linkage groups. The locations and direction of some of the QTL effects reported in this study suggest potentially favorable pleiotropic effects for the underlying genes. Further studies will be required to confirm these QTL in other populations so that they can be fine-mapped for potential applications in marker-assisted selection and management of beef cattle.

Biomarker discovery in animal health and disease: the application of post-genomic technologies

SUMMARY

The causes of many important diseases in animals are complex and multifactorial, which present unique challenges. Biomarkers indicate the presence or extent of a biological process, which is directly linked to the clinical manifestations and outcome of a particular disease. Identifying biomarkers or biomarker profiles will be an important step towards disease characterization and management of disease in animals. The emergence of post-genomic technologies has led to the development of strategies aimed at identifying specific and sensitive biomarkers from the thousands of molecules present in a tissue or biological fluid. This review will summarize the current developments in biomarker discovery and will focus on the role of transcriptomics, proteomics and metabolomics in biomarker discovery for animal health and disease.

Localization of Period 1 mRNA in the ruminant oocyte and investigations of its role in ovarian function

The clock gene Period 1 (Per1) may be a prolificacy gene, because it localized to the mouse oocyte and Per1-null drosophila shed fewer eggs. Because Per1 mapped to a region of mouse chromosome 11 syntenic to bovine chromosome 19 where a quantitative trait loci (QTL) for ovulation rate existed, we hypothesized that Per1 influenced folliculogenesis and ovulation rate in ruminants. Ovarian cortex was collected at slaughter on days 5, 12, 15, 17, and 20 after estrus for real-time RT-PCR evaluation of Per1 mRNA expression in Dorset (n = 18), Romanov (n = 10), Romanov/Dorset (n = 21), and Composite (n = 22) ewes. Ovarian cortex was also collected from cows selected for increased ovulation rate (n=37) or unselected controls (n = 28) on days 4, 5, and 6 of the estrous cycle for in situ hybridization and real-time RT-PCR. To examine the role of Per1 in early follicular development, ovarian cortex from neonatal calves (n = 5) was cultured for 10 days and Per1 mRNA levels were measured on day 0 and on day 10 of culture. The primers generated a 483bp amplicon with 100% sequence homology to bovine RIGUI-like protein (Per1). In silico mapping of this sequence placed Per1 on bovine chromosome 19; however, it was 20cM from the QTL. Per1 mRNA expression was unaffected by prolificacy, day of the cycle, or pregnancy status in ewes or cows. The riboprobe hybridized to oocytes of bovine preantral and antral follicles. In bovine ovarian cortical cultures on day 0, the tissue contained mostly primordial follicles (5.6+/-0.6 follicles/section); however, after 10 days in culture, the number of primordial follicles per section decreased (0.5 follicles/section) and the number of primary follicles increased as follicles activated (day 0 = 0.5+/- 0.6 versus day 10 = 10.4 +/-0.6 primary follicles/section; P < 0.001). Per1 mRNA did not change over time in culture. We conclude that Per1 mRNA is expressed by ruminant oocytes in preantral and antral follicles; however, its physiological role in mammalian ovarian function remains to be elucidated.

Comparative mapping of sheep chromosome 2q

Sheep chromosome 2q (OAR2q), which is homologous with human chromosome 2q (HSA2q), and cattle chromosome 2 (BTA2), is known to contain several loci contributing to carcass traits. However, the chromosomal rearrangements differentiating these chromosomes among the three species have not yet been determined and thus precise correspondences between the locations of sheep and human genes are not known. Twenty-six genes from HSA2q (2q21.1-->2q36) have been assigned to OAR2q by genetic linkage mapping to refine this area of the sheep genome. Seventy-six genes were initially selected from HSA2q. Sixty-eight percent of the PCR primer sets designed for these genes amplified successfully in sheep, and 34% amplified polymorphic products. Part of the proximal arm of OAR2q was found to be inverted compared with HSA2q. The breakpoint has been localised near the growth differentiation factor 8 gene (GDF8), spanning 380 kb between the positions of the hypothetical protein (FLJ20160) (HSA2:191008944-191075046) and glutaminase (GLS) (HSA2:191453847-191538510) (Build36.1).

Putative quantitative trait loci associated with the probability of contracting infectious bovine keratoconjunctivitis1,2

Infectious bovine keratoconjunctivitis, also known as pinkeye, is an economically important disease in cattle. The objective of this study was to detect QTL associated with infectious bovine keratoconjunctivitis in offspring from a Brahman x Hereford sire. The sire was mated to Hereford, Angus, and F1 cows to produce 288 offspring in 1994 and mated to MARC III ((1/4) Hereford, (1/4) Angus, (1/4) Red Poll, and (1/4) Pinzgauer) cows in 1996 to produce 259 offspring (547 animals total). Infectious bovine keratoconjunctivitis was diagnosed by physical examination in 36 animals of the family. Records included unilateral and bilateral frequency, but not severity. Records were binary: 0 for unaffected and 1 for affected cattle. A putative QTL for infectious bovine keratoconjunctivitis was identified on chromosome 1, with a maximum F-statistic (F = 10.15; P = 0.0015) at centimorgan 79 of the linkage group. The support interval spanned centimorgans 66 to 110. There was also evidence suggesting the presence of a QTL for infectious bovine keratoconjunctivitis on chromosome 20, with a maximum F-statistic (F = 10.35; P = 0.0014) at centimorgan 16 of the linkage group. The support interval ranged from centimorgan 2 to 35. This report provides the initial evidence of QTL for infectious bovine keratoconjunctivitis. Although a candidate gene was identified for one of the regions of interest, further studies are needed to identify the genetic basis of resistance to the disease.

Association of a single nucleotide polymorphism in SPP1 with growth traits and twinning in a cattle population selected for twinning rate

Continued validation of genetic markers for economically important traits is crucial to establishing marker-assisted selection as a tool in the cattle industry. The objective of the current study was to evaluate the association of a SNP (T(9)/T(10)) in the osteopontin gene (SPP1) with growth rate in a large cattle population spanning multiple generations and representing alleles from 12 founding breeds. This population has been maintained at the US Meat Animal Research Center since 1981 and subjected to selection for twinning rate. Phenotypic records for this population included twinning rate and ovulation rate, providing an opportunity to examine the potential effects of SPP1 genotype on reproductive traits. A set of 2,701 animals was geno-typed for the T(9)/T(10) polymorphism at SPP1. The geno-typic data, including previously genotyped markers on chromosome 6 (BTA6), were used in conjunction with pedigree information to estimate genotypic probabilities for all 14,714 animals with phenotypic records. The genotypic probabilities for females were used to calculate independent variables for regressions of additive, dominance, and imprinting effects. Genotypic regressions were fit as fixed effects in a mixed model analysis, in which each trait was analyzed in a 2-trait model where single births were treated as a separate trait from twin births. The association of the SPP1 marker with birth weight (P < 0.006), weaning weight (P < 0.007), and yearling weight (P < 0.003) was consistent with the previously reported effects of SPP1 genotype on yearling weight. Our data supports the conclusion that the SNP successfully tracks functional alleles affecting growth in cattle. The previously undetected effect of the SNP on birth and weaning weight suggests this particular SPP1 marker may explain a portion of the phenotypic variance explained by QTL for birth and HCW on BTA6.

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

We constructed a metric linkage disequilibrium (LD) map of bovine chromosome 6 (BTA6) on the basis of data from 220 SNPs genotyped on 433 Australian dairy bulls. This metric LD map has distances in LD units (LDUs) that are analogous to centimorgans in linkage maps. The LD map of BTA6 has a total length of 8.9 LDUs. Within the LD map, regions of high LD (represented as blocks) and regions of low LD (steps) are observed, when plotted against the integrated map in kilobases. At the most stringent block definition, namely a set of loci with zero LDU increase over the span of these markers, BTA6 comprises 40 blocks, accounting for 41% of the chromosome. At a slightly lower stringency of block definition (a set of loci covering a maximum of 0.2 LDUs on the LD map), up to 81% of BTA6 is spanned by 46 blocks and with 13 steps that are likely to reflect recombination hot spots. The mean swept radius (the distance over which LD is likely to be useful for mapping) is 13.3 Mb, confirming extensive LD in Holstein-Friesian dairy cattle, which makes such populations ideal for whole-genome association studies.

Comparative analysis of the bovine MHC class IIb sequence identifies inversion breakpoints and three unexpected genes

The bovine major histocompatibility complex (MHC) or BoLA is organized differently from typical mammalian MHCs in that a large portion of the class II region, called class IIb, has been transposed to a position near the centromere on bovine chromosome 23. Gene mapping indicated that the rearrangement resulted from a single inversion, but the boundaries and gene content of the inverted segment have not been fully determined. Here, we report the genomic sequence of BoLA IIb. Comparative sequence analysis with the human MHC revealed that the proximal inversion breakpoint occurred approximately 2.5 kb from the 3' end of the glutamate-cysteine ligase, catalytic subunit (GCLC) locus and that the distal breakpoint occurred about 2 kb from the 5' end from a divergent class IIDRbeta-like sequence designated DSB. Gene content, order and orientation of BoLA IIb are consistent with the single inversion hypothesis when compared with the corresponding region of the human class II MHC (HLA class II). Differences with HLA include the presence of a single histone H2B gene located between the proteasome subunit, beta type, 9 (PSMB9) and DMB loci and a duplicated TAP2 with a variant splice site. BoLA IIb spans approximately 450 kb DNA, with 20 apparently intact genes and no obvious pseudogenes. The region contains 227 simple sequence repeats (SSRs) and approximately 167 kb of retroviral-related repetitive DNA. Nineteen of the 20 genes identified in silico are supported by bovine EST data indicating that the functional gene content of BoLA IIb has not been diminished because it has been transposed from the remainder of BoLA genes.

Single nucleotide polymorphisms for integrative mapping in the Turkey (Meleagris gallopavo)

When multiple genetic maps exist for a species, integration of these maps requires a set of common markers be genotyped across the individual mapping populations. In the turkey, three genetic maps based on separate mapping populations are available. In this study, SNP-based markers were developed for integrating the cDNA/RFLP-based map (1) with microsatellite markers of the second-generation turkey genome map (2). Forty-eight primer sets were designed and tested and 33 (69%) correctly amplified turkey genomic DNA by PCR. Putative SNPs were detected in 20 (61%) of the amplified gene fragments, and 10 SNP markers were subsequently genotyped by PCR/RFLP for segregation analysis. Eight SNP markers were incorporated into the turkey genetic map.

A gene-based high-resolution comparative radiation hybrid map as a framework for genome sequence assembly of a bovine chromosome 6 region associated with QTL for growth, body composition, and milk performance traits

Background

A number of different quantitative trait loci (QTL) for various phenotypic traits, including milk production, functional, and conformation traits in dairy cattle as well as growth and body composition traits in meat cattle, have been mapped consistently in the middle region of bovine chromosome 6 (BTA6). Dense genetic and physical maps and, ultimately, a fully annotated genome sequence as well as their mutual connections are required to efficiently identify genes and gene variants responsible for genetic variation of phenotypic traits. A comprehensive high-resolution gene-rich map linking densely spaced bovine markers and genes to the annotated human genome sequence is required as a framework to facilitate this approach for the region on BTA6 carrying the QTL.

Results

Therefore, we constructed a high-resolution radiation hybrid (RH) map for the QTL containing chromosomal region of BTA6. This new RH map with a total of 234 loci including 115 genes and ESTs displays a substantial increase in loci density compared to existing physical BTA6 maps. Screening the available bovine genome sequence resources, a total of 73 loci could be assigned to sequence contigs, which were already identified as specific for BTA6. For 43 loci, corresponding sequence contigs, which were not yet placed on the bovine genome assembly, were identified. In addition, the improved potential of this high-resolution RH map for BTA6 with respect to comparative mapping was demonstrated. Mapping a large number of genes on BTA6 and cross-referencing them with map locations in corresponding syntenic multi-species chromosome segments (human, mouse, rat, dog, chicken) achieved a refined accurate alignment of conserved segments and evolutionary breakpoints across the species included.

Conclusion

The gene-anchored high-resolution RH map (1 locus/300 kb) for the targeted region of BTA6 presented here will provide a valuable platform to guide high-quality assembling and annotation of the currently existing bovine genome sequence draft to establish the final architecture of BTA6. Hence, a sequence-based map will provide a key resource to facilitate prospective continued efforts for the selection and validation of relevant positional and functional candidates underlying QTL for milk production and growth-related traits mapped on BTA6 and on similar chromosomal regions from evolutionary closely related species like sheep and goat. Furthermore, the high-resolution sequence-referenced BTA6 map will enable precise identification of multi-species conserved chromosome segments and evolutionary breakpoints in mammalian phylogenetic studies.