Butyrate Induced IGF2 Activation Correlated with Distinct Chromatin Signatures Due to Histone Modification

Histone modification has emerged as a very important mechanism regulating the transcriptional status of the genome. Insulin-like growth factor 2 (IGF2) is a peptide hormone controlling various cellular processes, including proliferation and apoptosis. H19 gene is closely linked to IGF2 gene, and IGF2 and H19 are reciprocally regulated imprinted genes. The epigenetic signature of H19 promoter (hypermethylation) on the paternal allele plays a vital role in allowing the expression of the paternal allele of IGF2.46 Our previous studies demonstrate that butyrate regulates the expression of IGF2 as well as genes encoding IGF Binding proteins. To obtain further understanding of histone modification and its regulatory potentials in controlling IGF2/H19 gene expression, we investigated the histone modification status of some key histones associated with the expression of IGF2/H19 genes in bovine cells using RNA-seq in combination with Chip-seq technology. A high-resolution map of the major chromatin modification at the IGF2/H19 locus induced by butyrate was constructed to illustrate the fundamental association of the chromatin modification landscape that may play a role in the activation of the IGF2 gene. High-definition epigenomic landscape mapping revealed that IGF2 and H19 have distinct chromatin modification patterns at their coding and promoter regions, such as TSSs and TTSs. Moreover, the correlation between the differentially methylated regions (DMRs) of IGF2/H19 locus and histone modification (acetylation and methylation) indicated that epigenetic signatures/markers of DNA methylation, histone methylation and histone acetylation were differentially distributed on the expressed IGF2 and silenced H19 genes. Our evidence also suggests that butyrate-induced regional changes of histone acetylation statusin the upstream regulation domain of H19 may be related to the reduced expression of H19 and strong activation of IGF2. Our results provided insights into the mechanism of butyrate-induced loss of imprinting (LOI) of IGF2 and regulation of gene expression by histone modification.

Fine mapping for Weaver syndrome in Brown Swiss cattle and the identification of 41 concordant mutations across NRCAM, PNPLA8 and CTTNBP2

Bovine Progressive Degenerative Myeloencephalopathy (Weaver Syndrome) is a recessive neurological disease that has been observed in the Brown Swiss cattle breed since the 1970's in North America and Europe. Bilateral hind leg weakness and ataxia appear in afflicted animals at 6 to 18 months of age, and slowly progresses to total loss of hind limb control by 3 to 4 years of age. While Weaver has previously been mapped to Bos taurus autosome (BTA) 4∶46-56 Mb and a diagnostic test based on the 6 microsatellite (MS) markers is commercially available, neither the causative gene nor mutation has been identified; therefore misdiagnosis can occur due to recombination between the diagnostic MS markers and the causative mutation. Analysis of 34,980 BTA 4 SNPs genotypes derived from the Illumina BovineHD assay for 20 Brown Swiss Weaver carriers and 49 homozygous normal bulls refined the Weaver locus to 48-53 Mb. Genotyping of 153 SNPs, identified from whole genome sequencing of 10 normal and 10 carrier animals, across a validation set of 841 animals resulted in the identification of 41 diagnostic SNPs that were concordant with the disease. Except for one intergenic SNP all are associated with genes expressed in nervous tissues: 37 distal to NRCAM, one non-synonymous (serine to asparagine) in PNPLA8, one synonymous and one non-synonymous (lysine to glutamic acid) in CTTNBP2. Haplotype and imputation analyses of 7,458 Brown Swiss animals with Illumina BovineSNP50 data and the 41 diagnostic SNPs resulted in the identification of only one haplotype concordant with the Weaver phenotype. Use of this haplotype and the diagnostic SNPs more accurately identifies Weaver carriers in both Brown Swiss purebred and influenced herds.

Identification of candidate transcription factor binding sites in the cattle genome

A resource that provides candidate transcription factor binding sites (TFBSs) does not currently exist for cattle. Such data is necessary, as predicted sites may serve as excellent starting locations for future omics studies to develop transcriptional regulation hypotheses. In order to generate this resource, we employed a phylogenetic footprinting approach—using sequence conservation across cattle, human and dog—and position-specific scoring matrices to identify 379,333 putative TFBSs upstream of nearly 8000 Mammalian Gene Collection (MGC) annotated genes within the cattle genome. Comparisons of our predictions to known binding site loci within the PCK1, ACTA1 and G6PC promoter regions revealed 75% sensitivity for our method of discovery. Additionally, we intersected our predictions with known cattle SNP variants in dbSNP and on the Illumina BovineHD 770k and Bos 1 SNP chips, finding 7534, 444 and 346 overlaps, respectively. Due to our stringent filtering criteria, these results represent high quality predictions of putative TFBSs within the cattle genome. All binding site predictions are freely available at http://bfgl.anri.barc.usda.gov/BovineTFBS/ or http://199.133.54.77/BovineTFBS.

Analysis of copy number variations in Holstein cows identify potential mechanisms contributing to differences in residual feed intake

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. In this study, we performed an initial analysis of copy number variations (CNVs) using BovineHD SNP genotyping data from 147 Holstein cows identified as having high or low feed efficiency as estimated by residual feed intake (RFI). We detected 443 candidate CNV regions (CNVRs) that represent 18.4 Mb (0.6 %) of the genome. To investigate the functional impacts of CNVs, we created two groups of 30 individual animals with extremely low or high estimated breeding values (EBVs) for RFI, and referred to these groups as low intake (LI; more efficient) or high intake (HI; less efficient), respectively. We identified 240 (~9.0 Mb) and 274 (~10.2 Mb) CNVRs from LI and HI groups, respectively. Approximately 30-40 % of the CNVRs were specific to the LI group or HI group of animals. The 240 LI CNVRs overlapped with 137 Ensembl genes. Network analyses indicated that the LI-specific genes were predominantly enriched for those functioning in the inflammatory response and immunity. By contrast, the 274 HI CNVRs contained 177 Ensembl genes. Network analyses indicated that the HI-specific genes were particularly involved in the cell cycle, and organ and bone development. These results relate CNVs to two key variables, namely immune response and organ and bone development. The data indicate that greater feed efficiency relates more closely to immune response, whereas cattle with reduced feed efficiency may have a greater capacity for organ and bone development.

Fine mapping of copy number variations on two cattle genome assemblies using high density SNP array

Background

Btau_4.0 and UMD3.1 are two distinct cattle reference genome assemblies. In our previous study using the low density BovineSNP50 array, we reported a copy number variation (CNV) analysis on Btau_4.0 with 521 animals of 21 cattle breeds, yielding 682 CNV regions with a total length of 139.8 megabases.

Results

In this study using the high density BovineHD SNP array, we performed high resolution CNV analyses on both Btau_4.0 and UMD3.1 with 674 animals of 27 cattle breeds. We first compared CNV results derived from these two different SNP array platforms on Btau_4.0. With two thirds of the animals shared between studies, on Btau_4.0 we identified 3,346 candidate CNV regions representing 142.7 megabases (~4.70%) of the genome. With a similar total length but 5 times more event counts, the average CNVR length of current Btau_4.0 dataset is significantly shorter than the previous one (42.7 kb vs. 205 kb). Although subsets of these two results overlapped, 64% (91.6 megabases) of current dataset was not present in the previous study. We also performed similar analyses on UMD3.1 using these BovineHD SNP array results. Approximately 50% more and 20% longer CNVs were called on UMD3.1 as compared to those on Btau_4.0. However, a comparable result of CNVRs (3,438 regions with a total length 146.9 megabases) was obtained. We suspect that these results are due to the UMD3.1 assembly's efforts of placing unplaced contigs and removing unmerged alleles. Selected CNVs were further experimentally validated, achieving a 73% PCR validation rate, which is considerably higher than the previous validation rate. About 20-45% of CNV regions overlapped with cattle RefSeq genes and Ensembl genes. Panther and IPA analyses indicated that these genes provide a wide spectrum of biological processes involving immune system, lipid metabolism, cell, organism and system development.

Conclusion

We present a comprehensive result of cattle CNVs at a higher resolution and sensitivity. We identified over 3,000 candidate CNV regions on both Btau_4.0 and UMD3.1, further compared current datasets with previous results, and examined the impacts of genome assemblies on CNV calling.

Recent applications of DNA sequencing technologies in food, nutrition and agriculture

Next-generation DNA sequencing technologies are able to produce millions of short sequence reads in a highthroughput, cost-effective fashion. The emergence of these technologies has not only facilitated genome sequencing but also changed the landscape of life sciences. This review surveys their recent applications in food, nutrition and agriculture ranging from whole-genome sequencing and resequencing, RNA-seq and ChIP-seq, structural, functional and comparative genomics to metagenomics and epigenetics. We already began to witness broad impacts of these DNA sequencing technologies for solving the complex biological problems in food, nutrition and agriculture. In this article, recent patent-based information is also included.

Copy number variation of individual cattle genomes using next-generation sequencing

Copy number variations (CNVs) affect a wide range of phenotypic traits; however, CNVs in or near segmental duplication regions are often intractable. Using a read depth approach based on next-generation sequencing, we examined genome-wide copy number differences among five taurine (three Angus, one Holstein, and one Hereford) and one indicine (Nelore) cattle. Within mapped chromosomal sequence, we identified 1265 CNV regions comprising ~55.6-Mbp sequence--476 of which (~38%) have not previously been reported. We validated this sequence-based CNV call set with array comparative genomic hybridization (aCGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH), achieving a validation rate of 82% and a false positive rate of 8%. We further estimated absolute copy numbers for genomic segments and annotated genes in each individual. Surveys of the top 25 most variable genes revealed that the Nelore individual had the lowest copy numbers in 13 cases (~52%, χ(2) test; P-value <0.05). In contrast, genes related to pathogen- and parasite-resistance, such as CATHL4 and ULBP17, were highly duplicated in the Nelore individual relative to the taurine cattle, while genes involved in lipid transport and metabolism, including APOL3 and FABP2, were highly duplicated in the beef breeds. These CNV regions also harbor genes like BPIFA2A (BSP30A) and WC1, suggesting that some CNVs may be associated with breed-specific differences in adaptation, health, and production traits. By providing the first individualized cattle CNV and segmental duplication maps and genome-wide gene copy number estimates, we enable future CNV studies into highly duplicated regions in the cattle genome.

Refinement of Bos taurus sequence assembly based on BAC-FISH experiments

Background

The sequencing of the cow genome was recently published (Btau_4.0 assembly). A second, alternate cow genome assembly (UMD2), based on the same raw sequence data, was also published. The two assemblies have been subsequently updated to Btau_4.2 and UMD3.1, respectively.

Results

We compared the Btau_4.2 and UMD3.1 alternate assemblies. Inconsistencies were grouped into three main categories: (i) DNA segments showing almost coincidental chromosomal mapping but discordant orientation (inversions); (ii) DNA segments showing a discordant map position along the same chromosome; and (iii) sequences present in one chromosomal assembly but absent in the corresponding chromosome of the other assembly. The latter category mainly consisted of large amounts of scaffolds that were unassigned in Btau_4.2 but successfully mapped in UMD3.1. We sampled 70 inconsistencies and identified appropriate cow BACs for each of them. These clones were then utilized in FISH experiments on cow metaphase or interphase nuclei in order to disambiguate the discrepancies. In almost all instances the FISH results agreed with the UMD3.1 assembly. Occasionally, however, the mapping data of both assemblies were discordant with the FISH results.

Conclusions

Our work demonstrates how FISH, which is assembly independent, can be efficiently used to solve assembly problems frequently encountered using the shotgun approach.

Assessment of genome integrity with array CGH in cattle transgenic cell lines produced by homologous recombination and somatic cell cloning

Background

Transgenic cattle carrying multiple genomic modifications have been produced by serial rounds of somatic cell chromatin transfer (cloning) of sequentially genetically targeted somatic cells. However, cloning efficiency tends to decline with the increase of rounds of cloning. It is possible that multiple rounds of cloning compromise the genome integrity or/and introduce epigenetic errors in the resulting cell lines, rendering a decline in cloning. To test these possibilities, we performed 9 high density array Comparative Genomic Hybridization (CGH) experiments to test the genome integrity in 3 independent bovine transgenic cell lineages generated from genetic modification and cloning. Our plan included the control hybridizations (self to self) of the 3 founder cell lines and 6 comparative hybridizations between these founders and their derived cell lines with either high or low cloning efficiencies.

Results

We detected similar amounts of differences between the control hybridizations (8, 13 and 39 differences) and the comparative analyses of both "high" and "low" cell lines (ranging from 7 to 57 with a mean of ~20). Almost 75% of the large differences (>10 kb) and about 45% of all differences shared the same type (loss or gain) and were located in nearby genomic regions across hybridizations. Therefore, it is likely that they were not true differences but caused by systematic factors associated with local genomic features (e.g. GC contents).

Conclusions

Our findings reveal that large copy number variations are less likely to arise during genetic targeting and serial rounds of cloning, fortifying the notion that epigenetic errors introduced from serial cloning may be responsible for the cloning efficiency decline.

Applications and case studies of the next-generation sequencing technologies in food, nutrition and agriculture

The next-generation sequencing technologies are able to produce millions of short sequence reads in a high-throughput, cost-effective fashion. The emergence of these technologies has not only facilitated genome sequencing but also started to change the landscape of life sciences. Here, I survey their major applications ranging from whole-genome sequencing and resequencing, single nucleotide polymorphism (SNP) and structural variation discovery, to mRNA and noncoding RNA profiling and protein-nucleic acid interaction assay. These case studies in structural, functional and comparative genomics, metagenomics, and epigenomics are providing a more complete picture of the genome structures and functions. In the near future, we will witness broad impacts of these next-generation sequencing technologies for solving the complex biological problems in food, nutrition and agriculture. In this article, recent patents based information is also included.

Analysis of copy number variations among diverse cattle breeds

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. Here, we describe the first systematic and genome-wide analysis of copy number variations (CNVs) in modern domesticated cattle using array comparative genomic hybridization (array CGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH). The array CGH panel included 90 animals from 11 Bos taurus, three Bos indicus, and three composite breeds for beef, dairy, or dual purpose. We identified over 200 candidate CNV regions (CNVRs) in total and 177 within known chromosomes, which harbor or are adjacent to gains or losses. These 177 high-confidence CNVRs cover 28.1 megabases or approximately 1.07% of the genome. Over 50% of the CNVRs (89/177) were found in multiple animals or breeds and analysis revealed breed-specific frequency differences and reflected aspects of the known ancestry of these cattle breeds. Selected CNVs were further validated by independent methods using qPCR and FISH. Approximately 67% of the CNVRs (119/177) completely or partially span cattle genes and 61% of the CNVRs (108/177) directly overlap with segmental duplications. The CNVRs span about 400 annotated cattle genes that are significantly enriched for specific biological functions, such as immunity, lactation, reproduction, and rumination. Multiple gene families, including ULBP, have gone through ruminant lineage-specific gene amplification. We detected and confirmed marked differences in their CNV frequencies across diverse breeds, indicating that some cattle CNVs are likely to arise independently in breeds and contribute to breed differences. Our results provide a valuable resource beyond microsatellites and single nucleotide polymorphisms to explore the full dimension of genetic variability for future cattle genomic research.

Analysis of recent segmental duplications in the bovine genome

BACKGROUND

Duplicated sequences are an important source of gene innovation and structural variation within mammalian genomes. We performed the first systematic and genome-wide analysis of segmental duplications in the modern domesticated cattle (Bos taurus). Using two distinct computational analyses, we estimated that 3.1% (94.4 Mb) of the bovine genome consists of recently duplicated sequences (>or= 1 kb in length, >or= 90% sequence identity). Similar to other mammalian draft assemblies, almost half (47% of 94.4 Mb) of these sequences have not been assigned to cattle chromosomes.

RESULTS

In this study, we provide the first experimental validation large duplications and briefly compared their distribution on two independent bovine genome assemblies using fluorescent in situ hybridization (FISH). Our analyses suggest that the (75-90%) of segmental duplications are organized into local tandem duplication clusters. Along with rodents and carnivores, these results now confidently establish tandem duplications as the most likely mammalian archetypical organization, in contrast to humans and great ape species which show a preponderance of interspersed duplications. A cross-species survey of duplicated genes and gene families indicated that duplication, positive selection and gene conversion have shaped primates, rodents, carnivores and ruminants to different degrees for their speciation and adaptation. We identified that bovine segmental duplications corresponding to genes are significantly enriched for specific biological functions such as immunity, digestion, lactation and reproduction.

CONCLUSION

Our results suggest that in most mammalian lineages segmental duplications are organized in a tandem configuration. Segmental duplications remain problematic for genome and assembly and we highlight genic regions that require higher quality sequence characterization. This study provides insights into mammalian genome evolution and generates a valuable resource for cattle genomics research.

Comparative analysis of Alu repeats in primate genomes

Using bacteria artificial chromosome (BAC) end sequences (16.9 Mb) and high-quality alignments of genomic sequences (17.4 Mb), we performed a global assessment of the divergence distributions, phylogenies, and consensus sequences for Alu elements in primates including lemur, marmoset, macaque, baboon, and chimpanzee as compared to human. We found that in lemurs, Alu elements show a broader and more symmetric sequence divergence distribution, suggesting a steady rate of Alu retrotransposition activity among prosimians. In contrast, Alu elements in anthropoids show a skewed distribution shifted toward more ancient elements with continual declining rates in recent Alu activity along the hominoid lineage of evolution. Using an integrated approach combining mutation profile and insertion/deletion analyses, we identified nine novel lineage-specific Alu subfamilies in lemur (seven), marmoset (one), and baboon/macaque (one) containing multiple diagnostic mutations distinct from their human counterparts-Alu J, S, and Y subfamilies, respectively. Among these primates, we show that that the lemur has the lowest density of Alu repeats (55 repeats/Mb), while marmoset has the greatest abundance (188 repeats/Mb). We estimate that approximately 70% of lemur and 16% of marmoset Alu elements belong to lineage-specific subfamilies. Our analysis has provided an evolutionary framework for further classification and refinement of the Alu repeat phylogeny. The differences in the distribution and rates of Alu activity have played an important role in subtly reshaping the structure of primate genomes. The functional consequences of these changes among the diverse primate lineages over such short periods of evolutionary time are an important area of future investigation.

Activation of SIRT1 by resveratrol represses transcription of the gene for the cytosolic form of phosphoenolpyruvate carboxykinase (GTP) by deacetylating hepatic nuclear factor 4alpha

The SIRT1 activators isonicotinamide (IsoNAM), resveratrol, fisetin, and butein repressed transcription of the gene for the cytosolic form of phosphoenolpyruvate carboxykinase (GTP) (PEPCK-C). An evolutionarily conserved binding site for hepatic nuclear factor (HNF) 4alpha (-272/-252) was identified, which was required for transcriptional repression of the PEPCK-C gene promoter caused by these compounds. This site contains an overlapping AP-1 binding site and is adjacent to the C/EBP binding element (-248/-234); the latter is necessary for hepatic transcription of PEPCK-C. AP-1 competed with HNF4alpha for binding to this site and also decreased HNF4alpha stimulation of transcription from the PEPCK-C gene promoter. Chromatin immunoprecipitation experiments demonstrated that HNF4alpha and AP-1, but not C/EBPbeta, reciprocally bound to this site prior to and after treating HepG2 cells with IsoNAM. IsoNAM treatment resulted in deacetylation of HNF4alpha, which decreased its binding affinity to the PEPCK-C gene promoter. In HNF4alpha-null Chinese hamster ovary cells, IsoNAM and resveratrol failed to repress transcription from the PEPCK-C gene promoter; overexpression of HNF4alpha in Chinese hamster ovary cells re-established transcriptional inhibition. Exogenous SIRT1 expression repressed transcription, whereas knockdown of SIRT1 by RNA interference reversed this effect. IsoNAM decreased the level of mRNA for PEPCK-C but had no effect on mRNA for glucose-6-phosphatase in AML12 mouse hepatocytes. We conclude that SIRT1 activation inhibited transcription of the gene for PEPCK-C in part by deacetylation of HNF4alpha. However, SIRT1 deacetylation of other key regulatory proteins that control PEPCK-C gene transcription also likely contributed to the inhibitory effect.

Identification of conserved regulatory elements in mammalian promoter regions: a case study using the PCK1 promoter

A systematic phylogenetic footprinting approach was performed to identify conserved transcription factor binding sites (TFBSs) in mammalian promoter regions using human, mouse and rat sequence alignments. We found that the score distributions of most binding site models did not follow the Gaussian distribution required by many statistical methods. Therefore, we performed an empirical test to establish the optimal threshold for each model. We gauged our computational predictions by comparing with previously known TFBSs in the PCK1 gene promoter of the cytosolic isoform of phosphoenolpyruvate carboxykinase, and achieved a sensitivity of 75% and a specificity of approximately 32%. Almost all known sites overlapped with predicted sites, and several new putative TFBSs were also identified. We validated a predicted SP1 binding site in the control of PCK1 transcription using gel shift and reporter assays. Finally, we applied our computational approach to the prediction of putative TFBSs within the promoter regions of all available RefSeq genes. Our full set of TFBS predictions is freely available at http://bfgl.anri.barc.usda.gov/tfbsConsSites.

Calibration of mutation rates reveals diverse subfamily structure of galliform CR1 repeats

Chicken Repeat 1 (CR1) repeats are the most abundant family of repeats in the chicken genome, with more than 200,000 copies accounting for approximately 80% of the chicken interspersed repeats. CR1 repeats are believed to have arisen from the retrotransposition of a small number of master elements, which gave rise to the 22 CR1 subfamilies as previously reported in Repbase. We performed a global assessment of the divergence distributions, phylogenies, and consensus sequences of CR1 repeats in the chicken genome. We identified and validated 57 chicken CR1 subfamilies and further analyzed the correlation between these subfamilies and their regional GC contents. We also discovered one novel lineage-specific CR1 subfamilies in turkeys when compared with chickens. We built an evolutionary tree of these subfamilies and concluded that CR1 repeats may play an important role in reshaping the structure of bird genomes.

The genome sequence of taurine cattle: a window to ruminant biology and evolution

To understand the biology and evolution of ruminants, the cattle genome was sequenced to about sevenfold coverage. The cattle genome contains a minimum of 22,000 genes, with a core set of 14,345 orthologs shared among seven mammalian species of which 1217 are absent or undetected in noneutherian (marsupial or monotreme) genomes. Cattle-specific evolutionary breakpoint regions in chromosomes have a higher density of segmental duplications, enrichment of repetitive elements, and species-specific variations in genes associated with lactation and immune responsiveness. Genes involved in metabolism are generally highly conserved, although five metabolic genes are deleted or extensively diverged from their human orthologs. The cattle genome sequence thus provides a resource for understanding mammalian evolution and accelerating livestock genetic improvement for milk and meat production.

Genomic expression analysis by single-cell mRNA differential display of quiescent CD8 T cells from tumour-infiltrating lymphocytes obtained from in vivo liver tumours

We performed a genomic study combining single-cell mRNA differential display and RNA subtractive hybridization to elucidate CD8 T-cell quiescence/ignorance. By comparing actively maintained quiescent CD8 T cells from liver tumour tumour-infiltrating lymphocytes (TILs) with quiescent T cells at the single-cell level, we identified differentially expressed candidate genes by high-throughput screening and comparative analysis of expressed sequence tags (ESTs). While genes for the T-cell receptor, tumour necrosis factor (TNF) receptor, TNF-related apoptosis inducing ligand (TRAIL) and perforin were down-regulated, key genes such as Tob, transforming growth factor (TGF)-beta, lung Krüpple-like factor (LKLF), Sno-A, Ski, Myc, Ets-2 repressor factor (ERF) and RE1-silencing transcription factor (REST/NRSF) complex were highly expressed in the quiescent TIL CD8 cells. Real-time polymerase chain reaction (PCR) further confirmed these results. A regulation model is proposed for actively maintained quiescence in CD8 T cells, including three components: up-regulation of the TGF-beta pathway, a shift in the MYC web and inhibition of the cell cycle.

Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds

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.

Detection of germline and somatic copy number variations in cattle

As a complement to the Bovine HapMap Consortium project, we initiated a systematic study of the copy numbervariation (CNV) within the same cattle population using array comparative genomic hybridization (array CGH). Oligonucleotide CGH arrays were designed and fabricated to cover all chromosomes with an average interval of 6 kb using the latest bovine genome assembly. In the initial screening, three Holstein bulls were selected to represent major paternal lineages of the Holstein breed with some maternal linkages between these lines. Dual-label hybridizations were performed using either Hereford L1 Dominette 01449 or L1 Domino 99375 as reference. The CNVs were represented by gains and losses of normalized fluorescence intensities relative to the reference. The data presented here, for the first time, demonstrated that significant amounts of germline and fewer somatic CNVs exist in cattle, that many CNVs are common both across diverse cattle breeds and among individuals within a breed, and that array CGH is an effective tool to systematically detect bovine CNV. Selected CNVs have been confirmed by independent methods using real-time (RT) PCR. The strategy used in this study, based on genome higher-orderarchitecture variation, is a powerful approach to generating resources for the identification of novel genomic variation and candidate genes for economically important traits.

Characterization of a novel microdeletion polymorphism on BTA5 in cattle

We present a detailed breakpoint mapping and population frequency analysis of a 214-kb microdeletion that removes multiple olfactory receptor genes. Using progressive rounds of PCR assays, we mapped the upstream and downstream breakpoints of this microdeletion event to approximately 1 and 12 kb genomic regions, respectively. We developed PCR-based genotyping assays, characterized a dairy cattle panel of 96 samples and found that the frequency of the deletion allele was over 51%. Our results indicated that this microdeletion is an ancient event occurring in one of the earlier founders, and that it has been stably inherited across generations in the North American dairy cattle population.

Genomic divergences among cattle, dog and human estimated from large-scale alignments of genomic sequences

BACKGROUND

Approximately 11 Mb of finished high quality genomic sequences were sampled from cattle, dog and human to estimate genomic divergences and their regional variation among these lineages.

RESULTS

Optimal three-way multi-species global sequence alignments for 84 cattle clones or loci (each >50 kb of genomic sequence) were constructed using the human and dog genome assemblies as references. Genomic divergences and substitution rates were examined for each clone and for various sequence classes under different functional constraints. Analysis of these alignments revealed that the overall genomic divergences are relatively constant (0.32-0.37 change/site) for pairwise comparisons among cattle, dog and human; however substitution rates vary across genomic regions and among different sequence classes. A neutral mutation rate (2.0-2.2 x 10(-9) change/site/year) was derived from ancestral repetitive sequences, whereas the substitution rate in coding sequences (1.1 x 10(-9) change/site/year) was approximately half of the overall rate (1.9-2.0 x 10(-9) change/site/year). Relative rate tests also indicated that cattle have a significantly faster rate of substitution as compared to dog and that this difference is about 6%.

CONCLUSION

This analysis provides a large-scale and unbiased assessment of genomic divergences and regional variation of substitution rates among cattle, dog and human. It is expected that these data will serve as a baseline for future mammalian molecular evolution studies.

Genome resources and comparative analysis tools for cardiovascular research

Disorders of the cardiovascular system are often caused by the interaction of genetic and environmental factors that jointly contribute to individual susceptibility. Genomic data and bioinformatics tools generated from genome projects, coupled with functional verification, offer novel approaches to study both rare single-gene and complex multigenic cardiovascular diseases. These approaches include gene mapping using genome variation, especially single-nucleotide polymorphisms and comparative genomics within and between species. This chapter illustrates the major genome resources, associated bioinformatics tools, and their potential application in cardiovascular research.

Neutralization and buffering capacity of leaves of sugar maple, largetooth aspen, paper birch and balsam fir

We compared the acidity, the external acid neutralizing capacity and the buffering capacity of leaves of four commercially important tree species, largetooth aspen (Populus grandidentata Michx.), sugar maple (Acer saccharum Marsh.), paper birch (Betula papyrifera Marsh.) and balsam fir (Abies balsamea (L.) Mill), at two sites of contrasting soil fertility in southern Quebec. External acid neutralizing capacity (ENC) of leaves was determined by measuring the change in pH induced by soaking fresh leaves in an acidic solution (pH 4.0) for two hours. The ENC was highest for largetooth aspen (14.3 micro equiv H(+) g(-1)), and lowest for sugar maple and balsam fir (< 5 micro equiv H(+) g(-1)). The buffering capacity index (BCI) was determined by measuring the amount of acid necessary to produce a change of 5 micro equiv H(+) in the leaf homogenate. The BCI ranged from 883 micro equiv H(+) g(-1) for largetooth aspen to less than 105 micro equiv H(+) g(-1) for sugar maple and balsam fir. Leaves of sugar maple and balsam fir had a lower internal pH and a higher percentage of ENC over BCI than paper birch and largetooth aspen. Overall, ENC was correlated with the concentration of all leaf nutrients except Ca, and BCI was correlated with Mg, N and Ca. The site effect was relatively unimportant for all variables.