High-throughput discovery of rare human nucleotide polymorphisms by Ecotilling

Human individuals differ from one another at only ~0.1% of nucleotide positions, but these single nucleotide differences account for most heritable phenotypic variation. Large-scale efforts to discover and genotype human variation have been limited to common polymorphisms. However, these efforts overlook rare nucleotide changes that may contribute to phenotypic diversity and genetic disorders, including cancer. Thus, there is an increasing need for high-throughput methods to robustly detect rare nucleotide differences. Toward this end, we have adapted the mismatch discovery method known as Ecotilling for the discovery of human single nucleotide polymorphisms. To increase throughput and reduce costs, we developed a universal primer strategy and implemented algorithms for automated band detection. Ecotilling was validated by screening 90 human DNA samples for nucleotide changes in 5 gene targets and by comparing results to public resequencing data. To increase throughput for discovery of rare alleles, we pooled samples 8-fold and found Ecotilling to be efficient relative to resequencing, with a false negative rate of 5% and a false discovery rate of 4%. We identified 28 new rare alleles, including some that are predicted to damage protein function. The detection of rare damaging mutations has implications for models of human disease.     

Discovery of single nucleotide polymorphisms and mutations by pyrosequencing

Comparative genomics, analyzing variation among individual genomes, is an area of intense investigation. DNA sequencing is usually employed to look for polymorphisms and mutations. Pyrosequencing, a real-time DNA sequencing method, is emerging as a popular platform for comparative genomics. Here we review the use of this technology for mutation scanning, polymorphism discovery and chemical haplotyping. We describe the methodology and accuracy of this technique and discuss how to reduce the cost for large-scale analysis.     

Efficient discovery of DNA polymorphisms in natural populations by Ecotilling

We have adapted the mutation detection technology used in Targeting Induced Local Lesions in Genomes (TILLING) to the discovery of polymorphisms in natural populations. The genomic DNA of a queried individual is mixed with a reference DNA and used to amplify a target 1-kbp region of DNA with asymmetrically labeled fluorescent primers. After heating and annealing, heteroduplexes are nicked at mismatched sites by the endonuclease CEL I and cut strands are visualized using Li-cor gel analyzers. Putative polymorphisms detected in one fluorescence channel can be verified by appearance of the opposite cut strand in the other channel. We demonstrated the efficiency of this technology, called Ecotilling, by the discovery in 150+ individuals of 55 haplotypes in five genes, ranging from sequences differing by a single nucleotide polymorphism to those representing complex haplotypes. The discovered polymorphisms were confirmed by sequencing and included base-pair changes, small insertions and deletions, and variation in microsatellite repeat number. Ecotilling allows the rapid detection of variation in many individuals and is cost effective because only one individual for each haplotype needs to be sequenced. The technology is applicable to any organism including those that are heterozygous and polyploid.     

Single nucleotide polymorphisms in the human E-cadherin gene

We report four DNA variants in the gene coding for the cell adhesion molecule E-cadherin. The polymorphisms affect codons 115, 133, 582 and the 3-noncoding region.     

Single nucleotide polymorphisms in the equine transferrin gene

Single nucleotide polymorphisms (SNPs) in exons 13, 15 and 16 of equine transferrin for common, rare and mutant variants were investigated. Compared with previous work a further 13 SNPs have been identified, allowing for the two previously identified clades to be subdivided into 11 groups. A combination of one or more of eight SNPs can be used to classify the equine variants into these 11 groups, since most are co-inherited. Putative sites of glycosylation in exons 13 and 16 showed no polymorphism, suggesting that presence or absence of sugar moieties does not lead to electrophoretic variation between the variants. Using the 26 SNPs currently identified in transferrin it is still not possible to differentiate variants F1 from F2, or D from H2, which represent 75% of the variants occurring in Thoroughbred equine population. This suggests that further SNPs exist in equine transferrin. The significance of the high level of variation in exon 15 is discussed.     

Single nucleotide polymorphisms in caprine calpastatin gene

The calpains and calpastatin (CAST) make up a major cytosolic proteolytic system, the calpain-calpastatin system, found in mammalian tissues. The relative levels of the components of the calpain-calpastatin system determine the extent of meat tenderization during postmortem storage. Calpastatin (CAST) is a protein inhibitor of the ubiquitous calcium-dependent proteases, μ-calpain, and m-calpain. Polymorphisms in the bovine, ovine and pig CAST gene have been associated with meat tenderness but little is known about how caprine CAST gene may affect goat meat quality traits. In this study we selected different parts of the CAST gene: (1) that have been previously reported to be polymorphic, intron 5 and 12 and 3’UTR; (2) first time explored (exon 3, 7 and 8 and part of intron 7 and 8) to investigate polymorphic status of caprine CAST gene. Using comparative sequencing ten novel SNPs located in exon 3 and intron 5, 7 and 8 were identified. Previously reported SNPs in intron 5, 3’UTR and intron 12 were absent. Sequence analysis revealed a non synonymous amino acid variation in exon 3, which would result in Lys/Arg substitution in the corresponding protein sequence. Considerable variation was detected in intronic regions. Twenty-four InDel were also recognized in intronic regions (15) and 3’UTR (9). All the sequences shared high homology with published bovine and ovine sequences. Three PCR-RFLP loci have been established for further analyzing genetic polymorphism in indigenous goats.     

Y-chromosome polymorphisms and the origins of the European gene pool

Gradients of allele frequencies have long been considered the main genetic characteristic of the European population, but mitochondrial DNA diversity seems to be distributed differently. One Alu insertion (YAP), five tetranucleotide (DYS19, DYS389B, DYS390, DYS391 and DYS393) and one trinucleotide (DYS392) microsatellite loci of the Y chromosome were analysed for geographical patterns in 59 European populations. Spatial correlograms showed clines for most markers, which paralleled the gradients previously observed for two RFLP polymorphisms. Effective separation times between populations were estimated from genetic distances at microsatellite loci. Even after correcting for the possible effects of continuous local gene flow, the most distant Indo-European-speaking populations seem to have separated no more than 7000 years ago. The clinal patterns and the estimated, recent separation times between populations jointly suggest that Y-chromosome diversity in Europe largely reflects a directional demic expansion, which is unlikely to have occurred before the Neolithic period.     

Discovery of single nucleotide polymorphisms in Lycopersicon esculentum by computer aided analysis of expressed sequence tags

Single nucleotide polymorphisms (SNPs) are useful for characterizing allelic variation, for genome-wide mapping, and as a tool for marker-assisted selection. Discovery of SNPs through de novo sequencing is inefficient within cultivated tomato (Lycopersicon esculentum Mill.) because the polymorphism rate is more than ten-fold lower than the sequencing error rate. The availability of expressed sequence tag (EST) data has made it feasible to discover putative SNPs in silico prior to experimental verification. By exploiting redundancy among EST data available for different varieties among 148,373 tomato ESTs, we have identified candidate SNPs for use within cultivated germplasm pools. 1,245 contigs having three EST sequences of Rio Grande and three EST sequences of TA496 were used for SNP discovery. We detected 1 SNP for every 8,500 bases analyzed, with 101 candidate SNPs in 44 genes identified. Sixty-six SNPs could be recognized by restriction enzymes, and subsequent experimental verification using restriction digestion or CEL I digestion confirmed 83% of the putative polymorphisms tested. SNPs between TA496 and Rio Grande have a high probability (53%) of detecting polymorphisms between other L. esculentum varieties. Twenty-six SNPs in 18 unigenes were mapped to specific chromosomes. Two SNPs, LEOH23 and LEOH37, were shown to be linked to quantitative trait loci contributing to fruit color within elite breeding populations. These results suggest that the growing databases of DNA sequence will yield information that facilitates improvement within the germplasm pools that have contributed to productive modern varieties.     

Discovery,validation and characterization of 1039 cattle single nucleotide polymorphisms

We identified ～13 000 putative single nucleotide polymorphisms (SNPs) by comparison of repeat‐masked BAC‐end sequences from the cattle RPCI‐42 BAC library with whole‐genome shotgun contigs of cattle genome assembly Btau 1.0. Genotyping of a subset of these SNPs was performed on a panel containing 186 DNA samples from 18 cattle breeds including 43 trios. Of 1039 SNPs confirmed as polymorphic in the panel, 998 had minor allele frequency ≥0.25 among unrelated individuals of at least one breed. When Btau 4.0 became available, 974 of these validated SNPs were assigned in silico to known cattle chromosomes, while 41 SNPs were mapped to unassigned sequence scaffolds, yielding one SNP every ～3 Mbp on average. Twenty‐four SNPs identified in Btau 1.0 were not mapped to Btau 4.0. Of the 1015 SNPs mapped to Btau 4.0, 959 SNPs had nucleotide bases identical in Btau 4.0 and Btau 1.0 contigs, whereas 56 bases were changed, resulting in the loss of the in silico SNP in Btau 4.0. Because these 1039 SNPs were all directly confirmed by genotyping on the multi‐breed panel, it is likely that the original polymorphisms were correctly identified. The 1039 validated SNPs identified in this study represent a new and useful resource for genome‐wide association studies and applications in animal breeding.     

Single nucleotide polymorphisms and disease gene mapping

Single nucleotide polymorphisms are the most important and basic form of variation in the genome, and they are responsible for genetic effects that produce susceptibility to most autoimmune diseases. The rapid development of databases containing very large numbers of single nucleotide polymorphisms, and the characterization of haplotypes and patterns of linkage disequilibrium throughout the genome, provide a unique opportunity to advance association strategies in common disease rapidly over the next few years. Only the careful use of these strategies and a clear understanding of their statistical limits will allow novel genetic determinants for many of the common autoimmune diseases to be determined.     

Identification of single nucleotide polymorphisms in the ovine acetyl-CoA carboxylase-alpha gene

Acetyl-CoA carboxylase (ACACA) is the rate-limiting enzyme in the biosynthesis of palmitic acid and long-chain fatty acids. The dietary intake of palmitic acid, which represents approximately 22% of sheep milk fatty acids, increases low-density lipoprotein (LDL) levels and the risk of developing human cardiovascular diseases. Following the candidate gene approach for improving sheep milk composition, and as a first step in assessing the possible influence of the ovine ACACA gene on milk fatty acid composition and its potential use as an animal genetic model of human atherosclerosis disease, we present here an investigation into the genetic variability of the ovine ACACA gene. We sequenced approximately 6.6 kb of ovine ACACA cDNA, including most of the coding sequence of the protein (except 348 bp), in Spanish Churra sheep. A total of 22 synonymous single nucleotide polymorphisms (SNPs) were identified in the analysed sequence, which were genotyped in a set of eight sheep breeds with different productive aptitudes (dairy, meat and double aptitudes). Two of the SNPs identified, SNP03 (c.1450T>C) and SNP15 (c.5134T>C), which appeared to be breed-specific variations, were situated in the gene sequence coding for the biotin-carboxylase (BC) and acetyl-CoA carboxyltransferase (ACCT) domains of the protein, respectively. Particularly interesting is SNP12 (c.4579G>A), which displayed higher frequencies in the dairy-specialised breeds relative to the meat-producing breeds. Moreover, in the dairy breeds studied, the frequency of this SNP showed a positive correlation with the degree of dairy specialisation. A previously described alternative splicing site (Ser-1200) affecting an important regulatory region of the enzyme was observed in one of the Churra animals. Despite the high genetic variability observed in this gene, none of the identified SNPs caused an amino acid change. However, these polymorphisms could be in linkage disequilibrium with other mutations showing a functional effect on the ACACA enzyme. Hence, the characterisations of the allelic variants reported herein lay the groundwork for evaluation of the potential use of these SNPs as genetic markers of fat content and fatty acid composition in sheep dairy products.     

Single nucleotide polymorphisms in an intron of the ovine calpastatin gene

Calpastatin is the specific inhibitor of the ubiquitous calcium-dependent proteases mu-calpain and m-calpain. Enzyme assay data from sheep and cattle inversely correlates post-mortem muscle calpastatin levels with ultimate meat tenderness. Genetic markers of meat quality may therefore be found linked to the calpastatin gene (CAST). A three-allele system detected by polymerase chain reaction-single strand conformational polymorphism (PCR-SSCP) has been observed in the ovine CAST. The three allele amplimers have been fully nucleotide sequenced and their differences in terms of single nucleotide polymorphism (SNPs) in the intron region of the amplimer are reported and compared to a consensus sequence of the orthologous region of the cattle CAST. A PCR-RFLP for more rapid CAST genotyping of all three ovine alleles was also developed.     

Four new nucleotide sequence polymorphisms in the LDL receptor gene detected by SSCP analysis

Seven nucleotide sequence polymorphisms were detected within exons of the low-density lipoprotein (LDL) receptor gene using single-strand conformation polymorphism (SSCP) analysis followed by direct sequence analysis on amplified DNA. Four nucleotide changes at nucleotide positions 1617, 1725, 2232, and 2635 were new nucleotide sequence polymorphisms not previously described. The remaining three nucleotide changes were identical with restriction fragment length polymorphisms and a previously reported nucleotide sequence polymorphism. These nucleotide sequence polymorphisms, detectable by SSCP analysis, are useful genetic markers for linkage analysis of the LDL receptor gene and familial hypercholesterolemia.     

Discovery and characterization of single nucleotide polymorphisms in coho salmon,Oncorhynchus kisutch

Molecular population genetic analyses have become an integral part of ecological investigation and population monitoring for conservation and management. Microsatellites have been the molecular marker of choice for such applications over the last several decades, but single nucleotide polymorphism (SNP) markers are rapidly expanding beyond model organisms. Coho salmon (Oncorhynchus kisutch) is native to the north Pacific Ocean and its tributaries, where it is the focus of intensive fishery and conservation activities. As it is an anadromous species, coho salmon typically migrate across multiple jurisdictional boundaries, complicating management and requiring shared data collection methods. Here, we describe the discovery and validation of a suite of novel SNPs and associated genotyping assays which can be used in the genetic analyses of this species. These assays include 91 that are polymorphic in the species and one that discriminates it from a sister species, Chinook salmon. We demonstrate the utility of these SNPs for population assignment and phylogeographic analyses, and map them against the draft trout genome. The markers constitute a large majority of all SNP markers described for coho salmon and will enable both population‐ and pedigree‐based analyses across the southern part of the species native range.     

Discovery and characterization of single nucleotide polymorphisms in Chinook salmon, Oncorhynchus tshawytscha

Molecular population genetics of non-model organisms has been dominated by the use of microsatellite loci over the last two decades. The availability of extensive genomic resources for many species is contributing to a transition to the use of single nucleotide polymorphisms (SNPs) for the study of many natural populations. Here we describe the discovery of a large number of SNPs in Chinook salmon, one of the world's most important fishery species, through large-scale Sanger sequencing of expressed sequence tag (EST) regions. More than 3 Mb of sequence was collected in a survey of variation in almost 132 kb of unique genic regions, from 225 separate ESTs, in a diverse ascertainment panel of 24 salmon. This survey yielded 117 TaqMan (5' nuclease) assays, almost all from separate ESTs, which were validated in population samples from five major stocks of salmon from the three largest basins on the Pacific coast of the contiguous United States: the Sacramento, Klamath and Columbia Rivers. The proportion of these loci that was variable in each of these stocks ranged from 86.3% to 90.6% and the mean minor allele frequency ranged from 0.194 to 0.236. There was substantial differentiation between populations with these markers, with a mean F(ST) estimate of 0.107, and values for individual loci ranging from 0 to 0.592. This substantial polymorphism and population-specific differentiation indicates that these markers will be broadly useful, including for both pedigree reconstruction and genetic stock identification applications.     

Identification of novel single nucleotide polymorphisms in the DGAT1 gene of buffaloes by PCR-SSCP

Diacylglycerol O-acyltransferase 1 (DGAT1) is a microsomal enzyme that catalyzes the final step of triglyceride synthesis. The DGAT1 gene is a strong functional candidate for determining milk fat content in cattle. In this work, we used PCR-SSCP (polymerase chain reaction-single-strand conformation polymorphism) and DNA sequencing to examine polymorphism in the region spanning exon 7 to exon 9 of the DGAT1 gene in Murrah and Pandharpuri buffaloes. Three alleles (A, B and C) and four novel single-nucleotide polymorphisms were identified in the buffalo DGAT1 gene. The frequencies of the alleles differed between the two buffalo breeds, with allele C being present in Murrah but not in Pandharpuri buffalo. The allele variation detected in this work may influence DGAT1 expression and function. The results described here could be useful in examining the association between the DGAT1 gene and milk traits in buffalo.     

Single nucleotide polymorphisms of myostatin gene in Chinese domestic horses

The myostatin gene (MSTN) is a genetic determinant of skeletal muscle growth. Single nucleotide polymorphisms (SNP) in MSTN are of importance due to their strong associations with horse racing performances. In this study, we screened the SNPs in MSTN gene in 514 horses from 15 Chinese horse breeds. Six SNPs (g.26 T > C, g.156 T > C, g.587A > G, g.598C > T, g.1485C > T, g.2115A > G) in MSTN gene were detected by sequencing and genotyped using PCR-RFLP method. The g.587A > G and g.598C > T residing in the 5′UTR region were novel SNPs identified by this study. The g.2115A > G which have previously been associated with racing performances were present in Chinese horse breeds, providing valuable genetic information for evaluating the potential racing performances in Chinese domestic breeds. The six SNPs together defined thirteen haplotypes, demonstrating abundant haplotype diversities in Chinese horses. Most of the haplotypes were shared among different breeds with no haplotype restricted to a specific region or a single horse breed. AMOVA analysis indicated that most of the genetic variance was attributable to differences among individuals without any significant contribution by the four geographical groups. This study will provide fundamental and instrumental genetic information for evaluating the potential racing performances of Chinese horse breeds.     

Interferon-gamma and interleukin-4 single nucleotide gene polymorphisms in Paracoccidioidomycosis

The gene polymorphisms interferon-gamma (IFN-γ) +874 T/A and interleukin (IL)-4 −590 C/T have been associated with the altered production of cytokines. Therefore, they might be indicative of the occurrence of Paracoccidioidomycosis (PCM) caused by Paracoccidioides brasiliensis. The analysis of single nucleotide polymorphism (SNP) at position +874 IFN-γ showed an increase occurrence of A/T genotype in both PCM patients and healthy individuals as control (HIC) (56% and 45%, respectively), while the allelic distribution showed 82% of A allele in the patients and 80% in the controls. The SNP of −590 IL-4 showed that C/T genotype was significantly (p < 0.05) more prevalent (39%) in PCM group compared to the HIC group (19%), while IL-4 C/C genotype was significantly less frequent (59%) in the patient group compared to the control group (81%). Otherwise, 41% of PCM patients and 19% of HIC individuals carried the IL-4 T allele. Stimulation of peripheral blood mononuclear cells (PBMC) from PCM patients with cell extract antigenic preparations (PbAg) as well as secreted and surface antigens (MEXO) of P. brasiliensis evidenced that there is no difference in the IFN-γ production related to A and T alleles between PCM and HIC individuals. However, with IL-4 production, PCM patients classified as C phenotype showed two times more IL-4 production than PCM patients classified as T phenotype and HIC controls. In conclusion, our results suggest that functional genetic variants in the IL-4 promoter could influence the production of IL-4 in PCM.