Use of bovine EST data and human genomic sequences to map 100 gene-specific bovine markers

A system to use bovine EST data in conjunction with human genomic sequence to improve the bovine linkage map over the entire genome or on specific chromosomes was evaluated. Bovine EST sequence was used to provide primer sequences corresponding to bovine genes, while human genomic sequence directed primer design to flank introns and produce amplicons of appropriate size for efficient direct sequencing. The sequence tagged sites (STS) produced in this way from the four sires of the MARC reference families were examined for single nucleotide polymorphisms (SNPs) that could be used to map the corresponding genes. With this approach, along with a primer/extension mass spectrometry SNP genotyping assay, 100 ESTs were placed on the bovine genetic linkage map. The first 70 were chosen at random from bovine EST–human genomic comparisons. An additional 30 ESTs were successfully mapped to bovine Chromosome 19 (BTA19), and comparison of the resulting BTA19 map to the position of the corresponding human orthologs on the HSA17 draft sequences revealed differences in the spacing and order of genes. Over 80% of successful amplicons contained SNPs, indicating that this is an efficient approach to generating EST-associated genetic markers. We have demonstrated the feasibility of constructing a linkage map based on SNPs associated with ESTs and the plausibility of utilizing EST, comparative mapping information, and human sequence data to target regions of the bovine genome for SNP marker development.     

Construction of a dense SNP map for bovine chromosome 6 to assist the assembly of the bovine genome sequence

A linkage map was constructed for bovine chromosome 6 (BTA6), using 399 single nucleotide polymorphisms (SNPs) detected primarily from PCR-resequencing. The efficiency of SNP detection was highly dependent on the source of sequence information chosen for primer design (BAC-end sequences, introns or promoters). The SNPs were used to build a linkage map comprising 104 cM on BTA6. The SNP order in the linkage map corresponded very well with radiation hybrid (RH) maps available for BTA6 as well as with expected positions in the human comparative map, but diverged significantly from the current assembly of the bovine genome (Btau_3.1). When performing linkage analysis with the marker order suggested from the Btau_3.1 we observed an expansion of the genetic map from 104 cM to 137 cM, strongly suggesting a reordering of scaffolds in the current version of the bovine genome assembly. The extent of LD on BTA6 was evaluated by calculating the average r ² for SNP pairs separated by given distances. The decline of LD was rapid with distance, such that r ² was 0.1 at 100 kb. Our results indicate that linkage mapping will be a valuable source of information for correcting errors in the current bovine assembly. These errors were sufficiently frequent to be of concern for the accuracy of mapping QTL with panels of SNPs whose positions are based on the current assembly.     

A Global View of 54,001 Single Nucleotide Polymorphisms (SNPs) on the Illumina BovineSNP50 BeadChip and Their Transferability to Water Buffalo

The Illumina BovineSNP50 BeadChip features 54,001 informative single nucleotide polymorphisms (SNPs) that uniformly span the entire bovine genome. Among them, 52,255 SNPs have locations assigned in the current genome assembly (Btau_4.0), including 19,294 (37%) intragenic SNPs (i.e., located within genes) and 32,961 (63%) intergenic SNPs (i.e., located between genes). While the SNPs represented on the Illumina Bovine50K BeadChip are evenly distributed along each bovine chromosome, there are over 14,000 genes that have no SNPs placed on the current BeadChip. Kernel density estimation, a non-parametric method, was used in the present study to identify SNP-poor and SNP-rich regions on each bovine chromosome. With bandwidth = 0.05 Mb, we observed that most regions have SNP densities within 2 standard deviations of the chromosome SNP density mean. The SNP density on chromosome X was the most dynamic, with more than 30 SNP-rich regions and at least 20 regions with no SNPs. Genotyping ten water buffalo using the Illumina BovineSNP50 BeadChip revealed that 41,870 of the 54,001 SNPs are fully scored on all ten water buffalo, but 6,771 SNPs are partially scored on one to nine animals. Both fully scored and partially/no scored SNPs are clearly clustered with various sizes on each chromosome. However, among 43,687 bovine SNPs that were successfully genotyped on nine and ten water buffalo, only 1,159 were polymorphic in the species. These results indicate that the SNPs sites, but not the polymorphisms, are conserved between two species. Overall, our present study provides a solid foundation to further characterize the SNP evolutionary process, thus improving understanding of within- and between-species biodiversity, phylogenetics and adaption to environmental changes.     

A global view of 54,001 single nucleotide polymorphisms (SNPs) on the Illumina BovineSNP50 BeadChip and their transferability to water buffalo

The Illumina BovineSNP50 BeadChip features 54,001 informative single nucleotide polymorphisms (SNPs) that uniformly span the entire bovine genome. Among them, 52,255 SNPs have locations assigned in the current genome assembly (Btau_4.0), including 19,294 (37%) intragenic SNPs (i.e., located within genes) and 32,961 (63%) intergenic SNPs (i.e., located between genes). While the SNPs represented on the Illumina Bovine50K BeadChip are evenly distributed along each bovine chromosome, there are over 14,000 genes that have no SNPs placed on the current BeadChip. Kernel density estimation, a non-parametric method, was used in the present study to identify SNP-poor and SNP-rich regions on each bovine chromosome. With bandwidth = 0.05 Mb, we observed that most regions have SNP densities within 2 standard deviations of the chromosome SNP density mean. The SNP density on chromosome X was the most dynamic, with more than 30 SNP-rich regions and at least 20 regions with no SNPs. Genotyping ten water buffalo using the Illumina BovineSNP50 BeadChip revealed that 41,870 of the 54,001 SNPs are fully scored on all ten water buffalo, but 6,771 SNPs are partially scored on one to nine animals. Both fully scored and partially/no scored SNPs are clearly clustered with various sizes on each chromosome. However, among 43,687 bovine SNPs that were successfully genotyped on nine and ten water buffalo, only 1,159 were polymorphic in the species. These results indicate that the SNPs sites, but not the polymorphisms, are conserved between two species. Overall, our present study provides a solid foundation to further characterize the SNP evolutionary process, thus improving understanding of within- and between-species biodiversity, phylogenetics and adaption to environmental changes.     

Single nucleotide polymorphism (SNP) discovery and linkage mapping of bovine cytokine genes

Polymorphic markers at bovine gene loci facilitate the integration of cattle genetic maps with those of humans and mice. To this end, 31 single nucleotide polymorphism (SNP) markers were developed for seven bovine chemokine genes. Loci were amplified from bovine genomic DNA by the polymerase chain reaction, and candidate amplicons were sequenced to determine their identity. Amplified loci from 24 founding parents and select progeny from a beef cattle reference population were sequenced and analyzed for SNPs. SNP haplotype alleles were determined by examining segregation patterns and used to establish the locus position on the bovine linkage map. Loci for growth-related proteins (GRO3, GRO1, and GROX) were clustered with the related CXC chemokine genes, interleukin (IL) 8, and epithelial cell inflammatory protein 1, at 84 cM from the centromeric end of the bovine chromosome (BTA) 6 linkage group. Bovine loci for a cluster of IL8 receptors, a stromal cell-derived factor 1, interferon-γ, and tumor necrosis factor-α were mapped at 90, 55, 59, and 34 cM, respectively, from the centromeric ends of the BTA 2, 28, 5, and 23 linkage groups. The positions of these bovine loci were compared with those of orthologous loci on the human map to refine the boundaries of conserved synteny. These seven loci provide examples of SNP development in which the efficiency was largely dependent on the availability of bovine genomic or cDNA sequence. The polymorphic nature of these SNP haplotype markers suggests that they will be useful for mapping complex traits in cattle, such as resistance to infectious disease. Received: 30 April 1999 / Accepted: 12 July 1999     

BcSNPdb: bovine coding region single nucleotide polymorphisms located proximal to quantitative trait loci

Bovine coding region single nucleotide polymorphisms located proximal to quantitative trait loci were identified to facilitate bovine QTL fine mapping research. A total of 692,763 bovine SNPs was extracted from 39,432 UniGene clusters, and 53,446 candidate SNPs were found to be a depth >3. In order to validate the in silico SNPs experimentally, 186 animals representing 14 breeds and 100 mixed breeds were analyzed. Genotyping of 40 randomly selected candidate SNPs revealed that 43% of these SNPs ranged in frequency from 0.009 to 0.498. To identify non-synonymous SNPs and to correct for possible frameshift errors in the ESTs at the predicted SNP positions, we designed a program that determines coding regions by protein-sequence referencing, and identified 17,735 nsSNPs. The SNPs and bovine quantitative traits loci informations were integrated into a bovine SNP data: BcSNPdb (http://snugenome.snu.ac.kr/BtcSNP/). Currently there are 43 different kinds of quantitative traits available. Thus, these SNPs would serve as valuable resources for exploiting genomic variation that influence economically and agriculturally important traits in cows.     

Genome‐wide association mapping and pathway analysis of leukosis incidence in a US Holstein cattle population

Bovine leukosis virus is an oncogenic virus that infects B cells, causing bovine leukosis disease. This disease is known to have a negative impact on dairy cattle production and, because no treatment or vaccine is available, finding a possible genetic solution is important. Our objective was to perform a comprehensive genetic analysis of leukosis incidence in dairy cattle. Data on leukosis occurrence, pedigree and molecular information were combined into multitrait GBLUP models with milk yield (MY) and somatic cell score (SCS) to estimate genetic parameters and to perform whole‐genome scans and pathway analysis. Leukosis data were available for 11 554 Holsteins daughters of 3002 sires from 112 herds in 16 US states. Genotypes from a 60K SNP panel were available for 961 of those bulls as well as for 2039 additional bulls. Heritability for leukosis incidence was estimated at about 8%, and the genetic correlations of leukosis disease incidence with MY and SCS were moderate at 0.18 and 0.20 respectively. The genome‐wide scan indicated that leukosis is a complex trait, possibly modulated by many genes. The gene set analysis identified many functional terms that showed significant enrichment of genes associated with leukosis. Many of these terms, such as G‐Protein Coupled Receptor Signaling Pathway, Regulation of Nucleotide Metabolic Process and different calcium‐related processes, are known to be related to retrovirus infection. Overall, our findings contribute to a better understanding of the genetic architecture of this complex disease. The functional categories associated with leukosis may be useful in future studies on fine mapping of genes and development of dairy cattle breeding strategies.     

SNPchip: R classes and methods for SNP array data

High-density single nucleotide polymorphism microarrays (SNP chips) provide information on a subject's genome, such as copy number and genotype (heterozygosity/homozygosity) at a SNP. While fluorescence in situ hybridization and karyotyping reveal many abnormalities, SNP chips provide a higher resolution map of the human genome that can be used to detect, e.g., aneuploidies, microdeletions, microduplications and loss of heterozygosity (LOH). As a variety of diseases are linked to such chromosomal abnormalities, SNP chips promise new insights for these diseases by aiding in the discovery of such regions, and may suggest targets for intervention. The R package SNPchip contains classes and methods useful for storing, visualizing and analyzing high density SNP data. Originally developed from the SNPscan web-tool, SNPchip utilizes S4 classes and extends other open source R tools available at Bioconductor. This has numerous advantages, including the ability to build statistical models for SNP-level data that operate on instances of the class, and to communicate with other R packages that add additional functionality. AVAILABILITY: The package is available from the Bioconductor web page at www.bioconductor.org. SUPPLEMENTARY INFORMATION: The supplementary material as described in this article (case studies, installation guidelines and R code) is available from http://biostat.jhsph.edu/~iruczins/publications/sm/     

HDAPD: a web tool for searching the disease-associated protein structures

Background

PCR-restriction fragment length polymorphism (RFLP) assay is a cost-effective method for SNP genotyping and mutation detection, but the manual mining for restriction enzyme sites is challenging and cumbersome. Three years after we constructed SNP-RFLPing, a freely accessible database and analysis tool for restriction enzyme mining of SNPs, significant improvements over the 2006 version have been made and incorporated into the latest version, SNP-RFLPing 2.

Results

The primary aim of SNP-RFLPing 2 is to provide comprehensive PCR-RFLP information with multiple functionality about SNPs, such as SNP retrieval to multiple species, different polymorphism types (bi-allelic, tri-allelic, tetra-allelic or indels), gene-centric searching, HapMap tagSNPs, gene ontology-based searching, miRNAs, and SNP500Cancer. The RFLP restriction enzymes and the corresponding PCR primers for the natural and mutagenic types of each SNP are simultaneously analyzed. All the RFLP restriction enzyme prices are also provided to aid selection. Furthermore, the previously encountered updating problems for most SNP related databases are resolved by an on-line retrieval system.

Conclusions

The user interfaces for functional SNP analyses have been substantially improved and integrated. SNP-RFLPing 2 offers a new and user-friendly interface for RFLP genotyping that can be used in association studies and is freely available at http://bio.kuas.edu.tw/snp-rflping2.     

The locus for bovine dilated cardiomyopathy maps to chromosome 18

Bovine dilated cardiomyopathy (BDCMP) is a severe and terminal disease of the heart muscle observed in Holstein-Friesian cattle over the last 30 years. There is strong evidence for an autosomal recessive mode of inheritance for BDCMP. The objective of this study was to genetically map BDCMP, with the ultimate goal of identifying the causative mutation. A whole-genome scan using 199 microsatellite markers and one SNP revealed an assignment of BDCMP to BTA18. Fine-mapping on BTA18 refined the candidate region to the MSBDCMP06-BMS2785 interval. The interval containing the BDCMP locus was confirmed by multipoint linkage analysis using the software loki. The interval is about 6.7 Mb on the bovine genome sequence (Btau 3.1). The corresponding region of HSA19 is very gene-rich and contains roughly 200 genes. Although telomeric of the marker interval, TNNI3 is a possible positional and a functional candidate for BDCMP given its involvement in a human form of dilated cardiomyopathy. Sequence analysis of TNNI3 in cattle revealed no mutation in the coding sequence, but there was a G-to-A transition in intron 6 (AJ842179:c.378+315G>A). The analysis of this SNP using the study's BDCMP pedigree did not conclusively exclude TNNI3 as a candidate gene for BDCMP. Considering the high density of genes on the homologous region of HSA19, further refinement of the interval on BTA18 containing the BDCMP locus is needed.     

A genome‐wide scan for copy number variations using high‐density single nucleotide polymorphism array in Simmental cattle

Copy number variations (CNVs) have recently been identified as promising sources of genetic variation, complementary to single nucleotide polymorphisms (SNPs). As a result, detection of CNVs has attracted a great deal of attention. In this study, we performed genome‐wide CNV detection using Illumina Bovine HD BeadChip (770k) data on 792 Simmental cattle. A total of 263 CNV regions (CNVRs) were identified, which included 137 losses, 102 gains and 24 regions classified as both loss and gain, covering 35.48 Mb (1.41%) of the bovine genome. The length of these CNVRs ranged from 10.18 kb to 1.76 Mb, with an average length of 134.78 kb and a median length of 61.95 kb. In 136 of these regions, a total of 313 genes were identified related to biological functions such as transmembrane activity and olfactory transduction activity. To validate the results, we performed quantitative PCR to detect nine randomly selected CNVRs and successfully confirmed seven (77.6%) of them. Our results present a map of cattle CNVs derived from high‐density SNP data, which expands the current CNV map of the cattle genome and provides useful information for investigation of genomic structural variation in cattle.     

Software solutions for the livestock genomics SNP array revolution

Since the beginning of the genomic era, the number of available single nucleotide polymorphism (SNP) arrays has grown considerably. In the bovine species alone, 11 SNP chips not completely covered by intellectual property are currently available, and the number is growing. Genomic/genotype data are not standardized, and this hampers its exchange and integration. In addition, software used for the analyses of these data usually requires not standard (i.e. case specific) input files which, considering the large amount of data to be handled, require at least some programming skills in their production. In this work, we describe a software toolkit for SNP array data management, imputation, genome‐wide association studies, population genetics and genomic selection. However, this toolkit does not solve the critical need for standardization of the genotypic data and software input files. It only highlights the chaotic situation each researcher has to face on a daily basis and gives some helpful advice on the currently available tools in order to navigate the SNP array data complexity.     

An interactive bovine <Emphasis Type="Italic">in silico</Emphasis> SNP database (IBISS)

An interactive bovine in silico SNP (IBISS) database has been created through the clustering and aligning of bovine EST and mRNA sequences. Approximately 324,000 EST and mRNA sequences were clustered to produce 29,965 clusters (producing 48,679 consensus sequences) and 48,565 singletons. A SNP screening regime was placed on variations detected in the multiple sequence alignment files to determine which SNPs are more likely to be real rather than sequencing errors. A small subset of predicted SNPs was validated on a diverse set of bovine DNA samples using PCR amplification and sequencing. Fifty percent of the predicted SNPs in the putative >1 category were polymorphic in the population sampled. The IBISS database represents more than just a SNP database; it is also a genomic database containing uniformly annotated predicted gene mRNA and protein sequences, gene structure, and genomic organization information.     

Genome-Wide Identification of Copy Number Variations in Chinese Holstein

Recent studies of mammalian genomes have uncovered the vast extent of copy number variations (CNVs) that contribute to phenotypic diversity. Compared to SNP, a CNV can cover a wider chromosome region, which may potentially incur substantial sequence changes and induce more significant effects on phenotypes. CNV has been becoming an alternative promising genetic marker in the field of genetic analyses. Here we firstly report an account of CNV regions in the cattle genome in Chinese Holstein population. The Illumina Bovine SNP50K Beadchips were used for screening 2047 Holstein individuals. Three different programes (PennCNV, cnvPartition and GADA) were implemented to detect potential CNVs. After a strict CNV calling pipeline, a total of 99 CNV regions were identified in cattle genome. These CNV regions cover 23.24 Mb in total with an average size of 151.69 Kb. 52 out of these CNV regions have frequencies of above 1%. 51 out of these CNV regions completely or partially overlap with 138 cattle genes, which are significantly enriched for specific biological functions, such as signaling pathway, sensory perception response and cellular processes. The results provide valuable information for constructing a more comprehensive CNV map in the cattle genome and offer an important resource for investigation of genome structure and genomic variation underlying traits of interest in cattle.     

An integrated database to enhance the identification of SNP markers for rice

The National Academy of Agricultural Science (NAAS) has developed a web-based marker database to provide information about SNP markers in rice. The database consists of three major functional categories: map viewing, marker searching and gene annotation. It provides 12,829 SNP markers information including gene location information on 12 chromosomes in rice. The annotation of SNP marker provides information such as marker name, EST number, gene definition and general marker information. Users are assisted in tracing any new structures of the chromosomes and gene positional functions using specific SNP markers. AVAILABILITY: The database is available for free at http://nabic.niab.go.kr/SNP/     

Evaluation of genotyping concordance for commercial bovine SNP arrays using quality‐assurance samples

SNP arrays are widely used in genetic research and agricultural genomics applications, and the quality of SNP genotyping data is of paramount importance. In the present study, SNP genotyping concordance and discordance were evaluated for commercial bovine SNP arrays based on two types of quality assurance (QA) samples provided by Neogen GeneSeek. The genotyping discordance rates (GDRs) between chips were on average between 0.06% and 0.37% based on the QA type I data and between 0.05% and 0.15% based on the QA type II data. The average genotyping error rate (GER) pertaining to single SNP chips, based on the QA type II data, varied between 0.02% and 0.08% per SNP and between 0.01% and 0.06% per sample. These results indicate that genotyping concordance rate was high (i.e. from 99.63% to 99.99%). Nevertheless, mitochondrial and Y chromosome SNPs had considerably elevated GDRs and GERs compared to the SNPs on the 29 autosomes and X chromosome. The majority of genotyping errors resulted from single allotyping errors, which also included the opposite instances for allele ‘dropout’ (i.e. from AB to AA or BB). Simultaneous allotyping errors on both alleles (e.g. mistaking AA for BB or vice versa) were relatively rare. Finally, a list of SNPs with a GER greater than 1% is provided. Interpretation of association effects of these SNPs, for example in genome‐wide association studies, needs to be taken with caution. The genotyping concordance information needs to be considered in the optimal design of future bovine SNP arrays.     

Genetic mapping of bovine T-cell receptor complex loci

Thymic-derived lymphocytes (T cells) recognize and respond to antigens through the mediation of the T-cell receptor complex. Here we report polymorphism at each of the five loci that encode the different components of the T-cell receptor complex in cattle. These genes were mapped on the bovine genome by genetic linkage analysis in the International Bovine Reference Panel (IBRP). These mapping data provide additional type I markers for linking the bovine genetic map with the human and mouse maps and also permit investigation of the effect of T-cell receptor polymorphism on immune responsiveness and disease susceptibility.     

RTPrimerDB: the real-time PCR primer and probe database

The real-time polymerase chain reaction (PCR) methodology has become increasingly popular for nucleic acids detection and/or quantification. As primer/probe design and experimental evaluation is time-consuming, we developed a public database application for the storage and retrieval of validated real-time PCR primer and probe sequence records. The integrity and accuracy of the data are maintained by linking to and querying other reference databases. RTPrimerDB provides free public access through the Web to perform queries and submit user based information. Primer/probe records can be searched for by official gene symbol, nucleotide sequence, type of application, detection chemistry, LocusLink or Single Nucleotide Polymorphism (SNP) identifier, and submitter's name. Each record is directly linked to LocusLink, dbSNP and/or PubMed to retrieve additional information on the gene/SNP for which the primers/probes are designed. Currently, the database contains primer/probe records for human, mouse, rat, fruit fly and zebrafish, and all current detection chemistries such as intercalating dyes (SYBR Green I), hydrolysis probes (Taqman), adjacent hybridizations probes and molecular beacons. Real-time PCR primer/probe records are available at http://www.realtimeprimerdatabase.ht.st.     

The bovine and ovine genomes contain multiple sequences homologous to the alpha-lactalbumin-encoding gene

Bovine and ovine (pseudo)genes homologous to the alpha-lactalbumin-encoding gene are described. In both cases, sequence analysis reveals homology extending downstream from exon 2. Southern analysis indicates the presence of a family of alpha-lactalbumin-related sequences in the bovine genome.