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.     

A protocol for TILLING and Ecotilling in plants and animals

We describe Targeting-Induced Local Lesions IN Genomes (TILLING), a reverse-genetic strategy for the discovery and mapping of induced mutations. TILLING is suitable for essentially any organism that can be mutagenized. The TILLING procedure has also been adapted for the discovery and cataloguing of natural polymorphisms, a method called Ecotilling. To discover nucleotide changes within a particular gene, PCR is performed with gene-specific primers that are end-labeled with fluorescent molecules. After PCR, samples are denatured and annealed to form heteroduplexes between polymorphic DNA strands. Mismatched base pairs in these heteroduplexes are cleaved by digestion with a single-strand specific nuclease. The resulting products are size-fractionated using denaturing polyacrylamide gel electrophoresis and visualized by fluorescence detection. The migration of cleaved products indicates the approximate location of nucleotide polymorphisms. Throughput is increased and costs are reduced by sample pooling, multi-well liquid handling and automated gel band mapping. Once genomic DNA samples have been obtained, pooled and arrayed, thousands of samples can be screened daily.     

Discovery of nucleotide polymorphisms in the Musa gene pool by Ecotilling

Musa (banana and plantain) is an important genus for the global export market and in local markets where it provides staple food for approximately 400 million people. Hybridization and polyploidization of several (sub)species, combined with vegetative propagation and human selection have produced a complex genetic history. We describe the application of the Ecotilling method for the discovery and characterization of nucleotide polymorphisms in diploid and polyploid accessions of Musa. We discovered over 800 novel alleles in 80 accessions. Sequencing and band evaluation shows Ecotilling to be a robust and accurate platform for the discovery of polymorphisms in homologous and homeologous gene targets. In the process of validating the method, we identified two single nucleotide polymorphisms that may be deleterious for the function of a gene putatively important for phototropism. Evaluation of heterozygous polymorphism and haplotype blocks revealed a high level of nucleotide diversity in Musa accessions. We further applied a strategy for the simultaneous discovery of heterozygous and homozygous polymorphisms in diploid accessions to rapidly evaluate nucleotide diversity in accessions of the same genome type. This strategy can be used to develop hypotheses for inheritance patterns of nucleotide polymorphisms within and between genome types. We conclude that Ecotilling is suitable for diversity studies in Musa, that it can be considered for functional genomics studies and as tool in selecting germplasm for traditional and mutation breeding approaches.     

Use of Ecotilling as an efficient SNP discovery tool to survey genetic variation in wild populations of Populus trichocarpa

Abstract Ecotilling was used as a simple nucleotide polymorphism (SNP) discovery tool to examine DNA variation in natural populations of the western black cottonwood, Populus trichocarpa, and was found to be more efficient than sequencing for large-scale studies of genetic variation in this tree. A publicly available, live reference collection of P. trichocarpa from the University of British Columbia Botanical Garden was used in this study to survey variation in nine different genes among individuals from 41 different populations. A large amount of genetic variation was detected, but the level of variation appears to be less than in the related species, Populus tremula, based on reported statistics for that tree. Genes examined varied considerably in their level of variation, from PoptrTB1 which had a single SNP, to PoptrLFY which had more than 23 in the 1000-bp region examined. Overall nucleotide diversity, measured as (Total), was relatively low at 0.00184. Linkage disequilibrium, on the other hand, was higher than reported for some woody plant species, with mean r2 equal to 0.34. This study reveals the potential of Ecotilling as a rapid genotype discovery method to explore and utilize the large pool of genetic variation in tree species.     

Genetic variation in the odorant receptor OR2J3 is associated with the ability to detect the "grassy" smelling odor, cis-3-hexen-1-ol

The ability to detect many odors varies among individuals; however, the contribution of genotype to this variation has been assessed for relatively few compounds. We have identified a genetic basis for the ability to detect the flavor compound cis-3-hexen-1-ol. This compound is typically described as "green grassy" or the smell of "cut grass," with variation in the ability to detect it linked to single nucleotide polymorphisms (SNPs) in a region on human chromosome 6 containing 25 odorant receptor genes. We have sequenced the coding regions of all 25 receptors across an ethnically mixed population of 52 individuals and identified 147 sequence variants. We tested these for association with cis-3-hexen-1-ol detection thresholds and found 3 strongly associated SNPs, including one found in a functional odorant receptor (rs28757581 in OR2J3). In vitro assays of 13 odorant receptors from the region identified 3 receptors that could respond to cis-3-hexen-1-ol, including OR2J3. This gene contained 5 predicted haplotypes across the 52 individuals. We tested all 5 haplotypes in vitro and several amino acid substitutions on their own, such as rs28757581 (T113A). Two amino acid substitutions, T113A and R226Q, impaired the ability of OR2J3 to respond to cis-3-hexen-1-ol, and together these two substitutions effectively abolished the response to the compound. The haplotype of OR2J3 containing both T113A and R226Q explains 26.4% of the variation in cis-3-hexen-1-ol detection in our study cohort. Further research is required to examine whether OR2J3 haplotypes explain variation in perceived flavor experience and the consumption of foods containing cis-3-hexen-1-ol.     

Haplotypes at ATM identify coding-sequence variation and indicate a region of extensive linkage disequilibrium

Genetic variation in the human population may lead to functional variants of genes that contribute to risk for common chronic diseases such as cancer. In an effort to detect such possible predisposing variants, we constructed haplotypes for a candidate gene and tested their efficacy in association studies. We developed haplotypes consisting of 14 biallelic neutral-sequence variants that span 142 kb of the ATM locus. ATM is the gene responsible for the autosomal recessive disease ataxia-telangiectasia (AT). These ATM noncoding single-nucleotide polymorphisms (SNPs) were genotyped in nine CEPH families (89 individuals) and in 260 DNA samples from four different ethnic origins. Analysis of these data with an expectation-maximization algorithm revealed 22 haplotypes at this locus, with three major haplotypes having frequencies > or = .10. Tests for recombination and linkage disequilibrium (LD) show reduced recombination and extensive LD at the ATM locus, in all four ethnic groups studied. The most striking example was found in the study population of European ancestry, in which no evidence for recombination could be discerned. The potential of ATM haplotypes for detection of genetic variants through association studies was tested by analysis of 84 individuals carrying one of three ATM coding SNPs. Each coding SNP was detected by association with an ATM haplotype. We demonstrate that association studies with haplotypes for candidate genes have significant potential for the detection of genetic backgrounds that contribute to disease.     

Large-scale detection of rare variants via pooled multiplexed next-generation sequencing: towards next-generation Ecotilling

Common variants, such as those identified by genome-wide association scans, explain only a small proportion of trait variation. Growing evidence suggests that rare functional variants, which are usually missed by genome-wide association scans, play an important role in determining the phenotype. We used pooled multiplexed next-generation sequencing and a customized analysis workflow to detect mutations in five candidate genes for lignin biosynthesis in 768 pooled Populus nigra accessions. We identified a total of 36 non-synonymous single nucleotide polymorphisms, one of which causes a premature stop codon. The most common variant was estimated to be present in 672 of the 1536 tested chromosomes, while the rarest was estimated to occur only once in 1536 chromosomes. Comparison with individual Sanger sequencing in a selected sub-sample confirmed that variants are identified with high sensitivity and specificity, and that the variant frequency was estimated accurately. This proposed method for identification of rare polymorphisms allows accurate detection of variation in many individuals, and is cost-effective compared to individual sequencing.     

Analysis of Human Mitochondrial DNA Polymorphisms in the Japanese Population

The highly polymorphic nature and high amplification efficiency of mitochondrial DNA (mtDNA) is valuable for the analysis of biological evidence in forensic casework, such as the identification of individuals and assignment of race/ethnicity. To be useful, a mtDNA polymorphism database for the Japanese population requires an understanding of the range of haplotype variation and phylogenies of mtDNA sequences. To extend current knowledge on the haplotypes in the Japanese population, this study defines new lineages and provides more detail about some of those previously described. We compared the hypervariable regions (HVRs) of 270 healthy, unrelated Japanese individuals and demonstrated 192 haplotypes. Combining HVR1 and HVR2, the genetic diversity was 0.9935, thus providing a high level of identification capability. Haplogroup status was defined for 160 individuals using HVR1, HVR2, and particular coding region polymorphisms; these individuals belonged to 94 haplotypes, four of which were new lineages. The complete mtDNA sequence was also determined from seven individuals.     

Human pepsinogen A (PGA): an informative gene complex located at 11q13

Summary Human pepsinogen (PGA) exhibits extensive polymorphism that can be detected both at the protein and the DNA level. We describe here two restriction fragment length polymorphisms, EcoRI and BglII, which provide for the detection of three of the most common PGA haplotypes (A, B, and C) in the United States population. The relationship of these polymorphisms to each PGA haplotype was determined by analysis of DNA from individuals exhibiting the corresponding protein phenotypes and by analysis of a series of human × mouse somatic cell hybrids containing the individual chromosome 11 homologous from heterozygous individuals exhibiting the AB and AC protein phenotypes. The use of the BglII polymorphism in combination with previously described EcoRI polymorphism provides a very informative marker of 11q13.     

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.