James Ireland

Analysis of molecular inversion probe performance for allele copy number determination

We have developed a new protocol for using molecular inversion probes to accurately and specifically measure allele copy number. The new protocol provides for significant improvements, including the reduction of input DNA (from 2 μg) by more than 25-fold (to 75 ng total genomic DNA), higher overall precision resulting in one order of magnitude lower false positive rate, and greater dynamic range with accurate absolute copy number up to 60 copies.

Optimal genotype determination in highly multiplexed SNP data

High-throughput genotyping technologies that enable large association studies are already available. Tools for genotype determination starting from raw signal intensities need to be automated, robust, and flexible to provide optimal genotype determination given the specific requirements of a study. The key metrics describing the performance of a custom genotyping study are assay conversion, call rate, and genotype accuracy. These three metrics can be traded off against each other. Using the highly multiplexed Molecular Inversion Probe technology as an example, we describe a methodology for identifying the optimal trade-off. The methodology comprises: a robust clustering algorithm and assessment of a large number of data filter sets. The clustering algorithm allows for automatic genotype determination. Many different sets of filters are then applied to the clustered data, and performance metrics resulting from each filter set are calculated. These performance metrics relate to the power of a study and provide a framework to choose the most suitable filter set to the particular study.

High-throughput method for analyzing methylation of CpGs in targeted genomic regions.

A unique microarray-based method for determining the extent of DNA methylation has been developed. It relies on a selective enrichment of the regions to be assayed by target amplification by capture and ligation (mTACL). The assay is quantitatively accurate, relatively precise, and lends itself to high-throughput determination using nanogram amounts of DNA. The measurements using mTACLs are highly reproducible and in excellent agreement with those obtained by sequencing (r = 0.94). In the present work, the methylation status of >145,000 CpGs from 5,472 promoters in 221 samples was measured. The methylation levels of nearby CpGs are correlated, but the correlation falls off dramatically over several hundred base pairs. In some instances, nearby CpGs have very different levels of methylation. Comparison of normal and tumor samples indicates that in tumors, the promoter regions of genes involved in differentiation and signaling are preferentially hypermethylated, whereas those of housekeeping genes remain hypomethylated. mTACL is a platform for profiling the state of methylation of a large number of CpG in many samples in a cost-effective fashion, and is capable of scaling to much larger numbers of CpGs than those collected here.

High-throughput, high-accuracy array-based resequencing

Although genomewide association studies have successfully identified associations of many common single-nucleotide polymorphisms (SNPs) with common diseases, the SNPs implicated so far account for only a small proportion of the genetic variability of tested diseases. It has been suggested that common diseases may often be caused by rare alleles missed by genomewide association studies. To identify these rare alleles we need high-throughput, high-accuracy resequencing technologies. Although array-based genotyping has allowed genomewide association studies of common SNPs in tens of thousands of samples, array-based resequencing has been limited for 2 main reasons: the lack of a fully multiplexed pipeline for high-throughput sample processing, and failure to achieve sufficient performance. We have recently solved both of these problems and created a fully multiplexed high-throughput pipeline that results in high-quality data. The pipeline consists of target amplification from genomic DNA, followed by allele enrichment to generate pools of purified variant (or nonvariant) DNA and ends with interrogation of purified DNA on resequencing arrays. We have used this pipeline to resequence ≈5 Mb of DNA (on 3 arrays) corresponding to the exons of 1,500 genes in >473 samples; in total >2,350 Mb were sequenced. In the context of this large-scale study we obtained a false positive rate of ≈1 in 500,000 bp and a false negative rate of ≈10%.

Functional single nucleotide polymorphism-based association studies

Association studies hold great promise for the elucidation of the genetic basis of diseases. Studies based on functional single nucleotide polymorphisms (SNPs) or on linkage disequilibrium (LD) represent two main types of designs. LD-based association studies can be comprehensive for common causative variants, but they perform poorly for rare alleles. Conversely, functional SNP-based studies are efficient because they focus on the SNPs with the highest a priori chance of being associated. Our poor ability to predict the functional effect of SNPs, however, hampers attempts to make these studies comprehensive. Recent progress in comparative genomics, and evidence that functional elements tend to lie in conserved regions, promises to change the landscape, permitting functional SNP association studies to be carried out that comprehensively assess common and rare alleles. SNP genotyping technologies are already sufficient for such studies, but studies will require continued genomic sequencing of multiple species, research on the functional role of conserved sequences and additional SNP discovery and validation efforts (including targeted SNP discovery to identify the rare alleles in functional regions). With these resources, we expect that comprehensive functional SNP association studies will soon be possible.

Rapid identification of somatic mutations in colorectal and breast cancer tissues using mismatch repair detection (MRD).

Mismatch repair detection (MRD) was used to screen 93 matched tumor‐normal sample pairs and 22 cell lines for somatic mutations in 30 cancer relevant genes. Using a starting amount of only 150 ng of genomic DNA, we screened 102 kb of sequence for somatic mutations in colon and breast cancer. A total of 152 somatic mutations were discovered, encompassing previously reported mutations, such as BRAF V600E and KRAS G12S, G12V, and G13D, as well as novel mutations, including some in genes in which somatic mutations have not previously been reported, such as MAP2K1 and MAP2K2. The distribution of mutations ranged widely within and across tumor types. The functional significance of many of these mutations is not understood, with patterns of selection only evident in KRAS and BRAF in colon cancer. These results present a novel approach to high‐throughput mutation screening using small amounts of starting material and reveal a mutation spectrum across 30 genes in a large cohort of breast and colorectal cancers. Hum Mutat 29(3), 441–450, 2008.

Large-scale characterization of public database SNPs causing non-synonymous changes in three ethnic groups

Single nucleotide polymorphisms (SNPs) that lead to non-synonymous changes in proteins may have functional effects and be subject to selection. Hence they are of particular interest in the study of genetic diseases. We have genotyped approximately 28,000 such SNPs in three ethnic populations (the HapMap plates) and ten primate species and analyzed these data for evidence of selection. We find SNPs predicted by PolyPhen to be damaging, have lower allele frequencies, and are particularly likely to be population-specific. We have also grouped SNPs by molecular function or biological process of the associated genes and find evidence that selection may be acting in concert on classes of genes.

Discovery of T Cell Receptor β Motifs Specific to HLA–B27–Positive Ankylosing Spondylitis by Deep Repertoire Sequence Analysis

Ankylosing spondylitis (AS), a chronic inflammatory disorder, has a notable association with HLA–B27. One hypothesis suggests that a common antigen that binds to HLA–B27 is important for AS disease pathogenesis. This study was undertaken to determine sequences and motifs that are shared among HLA–B27–positive AS patients, using T cell repertoire next‐generation sequencing.

Discovery of T-Cell Receptor Beta Motifs Specific to HLA-B27(+) Ankylosing Spondylitis by Deep Repertoire Sequence Analysis.

Ankylosing spondylitis (AS), a chronic inflammatory disorder, has a notable association with HLA–B27. One hypothesis suggests that a common antigen that binds to HLA–B27 is important for AS disease pathogenesis. This study was undertaken to determine sequences and motifs that are shared among HLA–B27–positive AS patients, using T cell repertoire next‐generation sequencing.