B. Murdoch

Genome-Wide Survey of SNP Variation Uncovers the Genetic Structure of Cattle Breeds

A survey of genetic diversity of cattle suggests two domestication events in Asia and selection by husbandry. Not Just Dinner on Legs Several thousand years ago, human beings realized the virtues of domesticating wild animals as easy meat. Soon other possibilities became apparent, and as revealed in a series of papers in this issue, early pastoralists became selective about breeding for wool, leather, milk, and muscle power. In two papers, Gibbs et al. report on the bovine genome sequence (p. 522; see the cover, the Perspective by Lewin, and the Policy Forum by Roberts) and trace the diversity and genetic history of cattle (p. 528), while Chessa et al. (p. 532) survey the occurrence of endogenous retroviruses in sheep and map their distribution to historical waves of human selection and dispersal across Europe. Finally, Ludwig et al. (p. 485) note the origins of variation in the coat-color of horses and suggest that it is most likely to have been selected for by humans in need of good-looking transport. The imprints of domestication and breed development on the genomes of livestock likely differ from those of companion animals. A deep draft sequence assembly of shotgun reads from a single Hereford female and comparative sequences sampled from six additional breeds were used to develop probes to interrogate 37,470 single-nucleotide polymorphisms (SNPs) in 497 cattle from 19 geographically and biologically diverse breeds. These data show that cattle have undergone a rapid recent decrease in effective population size from a very large ancestral population, possibly due to bottlenecks associated with domestication, selection, and breed formation. Domestication and artificial selection appear to have left detectable signatures of selection within the cattle genome, yet the current levels of diversity within breeds are at least as great as exists within humans.

Whole genome linkage disequilibrium maps in cattle

BackgroundBovine 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.ResultsLinkage disequilibrium was assessed using r2 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 r2 values. We found that the extent of linkage disequilibrium is no more than 0.5 Mb in these eight breeds of cattle.ConclusionLinkage 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.

Bisphenol A alters early oogenesis and follicle formation in the fetal ovary of the rhesus monkey

Widespread use of the endocrine disrupting chemical bisphenol A (BPA) in consumer products has resulted in nearly continuous human exposure. In rodents, low-dose exposures have been reported to adversely affect two distinct stages of oogenesis in the developing ovary: the events of prophase at the onset of meiosis in the fetal ovary and the formation of follicles in the perinatal ovary. Because these effects could influence the reproductive longevity and success of the exposed individual, we conducted studies in the rhesus monkey to determine whether BPA induces similar disturbances in the developing primate ovary. The routes and levels of human exposure are unclear; hence, two different exposure protocols were used: single daily oral doses and continuous exposure via subdermal implant. Our analyses of second trimester fetuses exposed at the time of meiotic onset suggest that, as in mice, BPA induces subtle disturbances in the prophase events that set the stage for chromosome segregation at the first meiotic division. Our analyses of third-trimester fetuses exposed to single daily oral doses during the time of follicle formation revealed an increase in multioocyte follicles analogous to that reported in rodents. However, two unique phenotypes were evident in continuously exposed animals: persistent unenclosed oocytes in the medullary region and small, nongrowing oocytes in secondary and antral follicles. Because effects on both stages of oogenesis were elicited using doses that yield circulating levels of BPA analogous to those reported in humans, these findings raise concerns for human reproductive health.

Assessment of positional candidate genes myf5 and igf1 for growth on bovine chromosome 5 in commercial lines of Bos taurus.

Quantitative trait loci for growth traits in beef cattle have been previously reported and fine-mapped in three chromosomal regions of 0 to 30 cM, 55 to 70 cM, and 70 to 80 cM of bovine chromosome 5. In this study, we further examined the association between gene-specific single nucleotide polymorphisms (SNP) of two positional candidate genes, bovine myogenic factor 5 (myf5) and insulin-like growth factor-1 (igf1), in the QTL regions and the birth weight (BWT), preweaning average daily gain (PWADG), and average daily gain on feed (ADGF) in commercial lines of Bos taurus. The QTL regions for the growth traits identified using a haplotype association analysis, which included the gene-specific SNP markers for both genes in this study, were in agreement with previous studies. The gene-specific SNP marker association analysis indicated that the SNP in myf5 had a significant additive effect on PWADG in the M1 line of Beefbooster Inc. (P < 0.10), and a significant additive effect (P < 0.05) and a significant dominance effect (P < 0.10) on ADGF in the M3 line of Beefbooster Inc. When the data from the two commercial lines were pooled, the SNP in myf5 showed a significant association with PWADG (P < 0.10) and with ADGF (P < 0.05). The association between the SNP and BWT, however, did not reach a significance level in the M1 line, the M3 line, or across the lines. For igf1, no significant association between the SNP and the growth traits was detected in either the M1 line or the M3 line, whereas there was only a significant dominance effect (P < 0.10) on BWT detected for the SNP in igfl when the data from the two commercial lines were pooled. These results suggest that myf5 is a strong candidate gene that influences PWADG and ADGF in beef cattle. The SNP of igf1 may not be a causative or close to the causative mutation that affects the three growth traits in the populations of beef cattle examined in this study. Other SNP of igf1 and myf5 or other genes in their respective chromosomal regions, however, should also be studied.

An assessment of population structure in eight breeds of cattle using a whole genome SNP panel

BackgroundAnalyses of population structure and breed diversity have provided insight into the origin and evolution of cattle. Previously, these studies have used a low density of microsatellite markers, however, with the large number of single nucleotide polymorphism markers that are now available, it is possible to perform genome wide population genetic analyses in cattle. In this study, we used a high-density panel of SNP markers to examine population structure and diversity among eight cattle breeds sampled from Bos indicus and Bos taurus.ResultsTwo thousand six hundred and forty one single nucleotide polymorphisms (SNPs) spanning all of the bovine autosomal genome were genotyped in Angus, Brahman, Charolais, Dutch Black and White Dairy, Holstein, Japanese Black, Limousin and Nelore cattle. Population structure was examined using the linkage model in the program STRUCTURE and Fst estimates were used to construct a neighbor-joining tree to represent the phylogenetic relationship among these breeds.ConclusionThe whole-genome SNP panel identified several levels of population substructure in the set of examined cattle breeds. The greatest level of genetic differentiation was detected between the Bos taurus and Bos indicus breeds. When the Bos indicus breeds were excluded from the analysis, genetic differences among beef versus dairy and European versus Asian breeds were detected among the Bos taurus breeds. Exploration of the number of SNP loci required to differentiate between breeds showed that for 100 SNP loci, individuals could only be correctly clustered into breeds 50% of the time, thus a large number of SNP markers are required to replace the 30 microsatellite markers that are currently commonly used in genetic diversity studies.

A physical map of the bovine genome

BackgroundCattle 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.ResultsA 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.ConclusionFurther 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.

Fine mapping quantitative trait loci for feed intake and feed efficiency in beef cattle

Feed intake and feed efficiency are economically important traits in beef cattle because feed is the greatest variable cost in production. Feed efficiency can be measured as feed conversion ratio (FCR, intake per unit gain) or residual feed intake (RFI, measured as DMI corrected for BW and growth rate, and sometimes a measure of body composition, usually carcass fatness, RFI(bf)). The goal of this study was to fine map QTL for these traits in beef cattle using 2,194 markers on 24 autosomes. The animals used were from 20 half-sib families originating from Angus, Charolais, and University of Alberta Hybrid bulls. A mixed model with random sire and fixed QTL effect nested within sire was used to test each location (cM) along the chromosomes. Threshold levels were determined at the chromosome and genome levels using 20,000 permutations. In total, 4 QTL exceeded the genome-wise threshold of P < 0.001, 3 exceeded at P < 0.01, 17 at P < 0.05, and 30 achieved significance at the chromosome-wise threshold level (at least P < 0.05). No QTL were detected on BTA 8, 16, and 27 above the 5% chromosome-wise significance threshold for any of the traits. Nineteen chromosomes contained RFI QTL significant at the chromosome-wise level. The RFI(bf) QTL results were generally similar to those of RFI, the positions being similar, but occasionally differing in the level of significance. Compared with RFI, fewer QTL were detected for both FCR and DMI, 12 and 4 QTL, respectively, at the genome-wise thresholds. Some chromosomes contained FCR QTL, but not RFI QTL, but all DMI QTL were on chromosomes where RFI QTL were detected. The most significant QTL for RFI was located on BTA 3 at 82 cM (P = 7.60 x 10(-5)), for FCR on BTA 24 at 59 cM (P = 0.0002), and for DMI on BTA 7 at 54 cM (P = 1.38 x 10(-5)). The RFI QTL that showed the most consistent results with previous RFI QTL mapping studies were on BTA 1, 7, 18, and 19. The identification of these QTL provides a starting point to identify genes affecting feed intake and efficiency for use in marker-assisted selection and management.

Gene Expression in the Fetal Mouse Ovary Is Altered by Exposure to Low Doses of Bisphenol A

Evidence from experimental studies suggests that fetal exposure to the endocrine-disrupting chemical bisphenol A (BPA) has adverse reproductive effects in both males and females. Studies from our laboratory suggest that exposure to the developing female fetus produces a unique, multigenerational effect. Specifically, maternal exposure affects the earliest stages of oogenesis in the developing fetal ovary, and the resulting subtle meiotic defects increase the likelihood that embryos produced by the exposed female in adulthood (i.e., the grandchildren) will be chromosomally abnormal. To understand the impact of BPA on the developing ovary, we conducted expression studies to characterize gene expression changes in the fetal ovary that result from BPA exposure. We first tested the validity of the approach, asking whether we could reliably detect temporal changes in expression levels of meiotic genes in controls. As anticipated, we were able to identify appropriate increases in expression in meiotic, but in few other, genes. Intriguingly, this analysis provided data on a small set of genes for which timing and expression changes suggest that they may have important and heretofore unrecognized meiotic roles. After verifying the utility of our approach, we focused our analysis on BPA-exposed animals. We found modest, but significant, changes in gene expression in the fetal ovaries from exposed fetuses. The first changes were evident within 24 h of exposure, and the most extensive changes correlated with the onset of meiosis. Furthermore, gene ontology analysis suggested that BPA acts to down-regulate mitotic cell-cycle genes, raising the possibility that fetal BPA exposure may act to limit expansion of the primordial germ cell population.

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 polymorphisms influencing feed intake and efficiency in beef cattle

Feed efficiency is an economically important trait in beef cattle. Net feed efficiency, measured as residual feed intake (RFI), is the difference between actual feed intake and the predicted feed intake required for maintenance and gain of the animal. SNPs that show associations with RFI may be useful quantitative trait nucleotides for marker-assisted selection. This study identified associations between SNPs underlying five RFI QTL on five bovine chromosomes (BTA2, 5, 10, 20 and 29) with measures of dry matter intake (DMI), RFI and feed conversion ratio (FCR) in beef cattle. Six SNPs were found to have effects on RFI (P < 0.05). The largest single SNP allele substitution effect for RFI was -0.25 kg/day located on BTA2. The combined effects of the SNPs found significant in this experiment explained 6.9% of the phenotypic variation of RFI. Not all the RFI SNPs showed associations with DMI and FCR even though these traits are highly correlated with RFI (r = 0.77 and r = 0.62 respectively). This shows that these SNPs may be affecting the underlying biological mechanisms of feed efficiency beyond feed intake control and weight gain efficiency. These SNPs can be used in marker-assisted selection but first it will be important to verify these effects in independent populations of cattle.