|
|
||||||||
Technical Briefs |
1 Department of Medicine, Thomas Jefferson University, Jefferson Medical College, Philadelphia, PA 19104
2 Unitadi Ricerca in Gastroenterologia, IRCCS Casa Sollievo della Sofferenza, I-71013 San Giovanni Rotondo, Italy
3 Pyrosequencing, Inc., Westborough, MA 01581
4 IRCCS C.S.S., I-71013 San Giovanni Rotondo, Italy
5 C.S.S.-Mendel, 00161 Rome, Italy
6 Dipartimento di Medicina Sperimentale e Patologia, Universitá "La Sapienza", 00161 Rome, Italy
7 Istituto e Servizio di Gastroenterologia, IRCCS Ospedale Maggiore e Universitá di Milano, 20122 Milan, Italy
8 Department of Biomedical Research, Nemours Childrens Clinic, Wilmington, DE 19880
9 Dipartimento di Oncologia, Biologia e Genetica, Universitá di Genova, 16146 Genoa, Italy
aaddress correspondence to this author at: 406 Medical Office Bldg., 1100 Walnut St., Philadelphia, PA 19107; fax 215-955-6905, e-mail paolo.fortina{at}jefferson.edu
We describe here a method for detection of three common single-nucleotide polymorphisms (SNPs) in the inflammatory bowel disease (IBD)-related CARD15 gene. We used a two-step approach with multiplex PCR and pyrosequencing. Although the three loci are spaced far from one another, the built-in multiplexing capability of pyrosequencing reactions to extend for numerous bases was used with appropriate primers and the design of efficient deoxynucleoside triphosphate (dNTP) dispensation orders to screen for the variants simultaneously.
IBDs are chronic relapsing disorders affecting the gastrointestinal tract with a population prevalence in Western countries of 1:1000 (1). The complex trait is subdivided on the basis of clinical and histologic features into two main subtypes: Crohn disease (CD) and ulcerative colitis (UC) (2)(3)(4)(5). Whereas UC is limited to the rectal and colonic mucosal layers, CD can involve any part of the intestine, most frequently the terminal ileum and colon (6)(7)(8)(9). Several candidate loci for IBD have been identified on the basis of whole-genome scans in different populations (10)(11)(12)(13)(14)(15)(16)(17)(18). Of these, the one on chromosome 16 named IBD1 has been confirmed in several independent studies (11)(12)(15)(16)(17)(18). Three SNPs in the CARD15 gene, located in the IBD1 region, have been demonstrated to confer susceptibility to CD, but not to UC, in several populations and to correlate with an ileal site of the disease (6)(7)(8)(9)(10).
In this report, using a previously described approach for SNP/mutation analysis (19), we developed and validated a protocol to characterize simultaneously the R702W (SNP8), G908R (SNP12), and L1007fsinsC (SNP13) variants, which are associated with increased risk of CD (20)(21)(22). Although additional CARD15 variants have been found to be associated to CD, these are the three most common and are found in almost all populations tested (7).
A total of 141, 163, and 61 individuals from 47 CD, 49 UC, and 18 mixed (containing affected members with both CD and UC) families, respectively, and 108 unrelated healthy controls were included in the study after informed consent. Genomic DNA samples were purified from peripheral blood leukocytes according to standard protocols or with the Puregene DNA isolation reagent set (Gentra System).
Multiplex PCR reactions were performed on samples to be genotyped for the three variants (R702W, G908R, and L1007fsinsC) in the CARD15 gene. PCR and Pyrosequencing sequencing primers were designed according to Pyrosequencing guidelines (http://www.pyrosequencing.com). Briefly, three fragments of 106, 117, and 104 bp, respectively, encompassing the three SNPs were PCR-amplified using forward and the reverse primer pairs R702W-F (5'-TTC CTG GCA GGG CTG TT-3') and R702W-R (5'-biotin-CCA GAC ACC AGC GGG C-3'), G908R-F (5'-TAG AGG GAG GAG GAC TGT TAG TT-3') and G908R-R (5'-biotin-CCC CCT CGT CAC CCA CT-3'), and L1007fsinsC-F (5'-TCC TTA CCA GAC TTC CAG GAT G-3') and L1007fsinsC-R (5'-biotin-TCT TCT TTT CCA GGT TGT CCA-3') with one of the two PCR primers (Integrated DNA Technologies) for each SNP containing biotin at the 5'-terminus to facilitate isolation of a single-stranded template for the Pyrosequencing assay. All primers were HPLC-purified.
Because a template hairpin was predicted by SNP Primer Design software Ver. 1.0 (Pyrosequencing AB), a nonmatching T (underlined and bolded above) was added to the 5' end of the nonbiotinylated PCR primers for the G908R and L1007fsinsC SNPs to prevent extension from these templates. The triplex PCR reactions were done in a final volume of 50 µL containing GeneAmp 1x Buffer II (10 mM Tris-HCl, pH 8.3; 50 mM KCl), 2.5 mM MgCl2, 3 U of AmpliTaq Gold DNA polymerase (Applied Biosystems), 200 nM each primer, 1 mM dNTPs (Amersham Pharmacia Biotech), and 50 ng of the genomic DNA. Reactions were performed with initial denaturation at 95 °C for 10 min, followed by 35 cycles at 95 °C for 1 min, annealing at 62 °C for 1 min, and extension at 72 °C for 1 min. Final extension was 72 °C for 10 min.
Bound biotinylated single-stranded DNA was generated using protocols recommended by the manufacturer (PSQTM HS 96A System Sample Preparation Kit; Pyrosequencing AB). The biotinylated triplex PCR products (5 µL) were bound to streptavidin-Sepharose HP beads (2 µL; Amersham Biosciences) and incubated with binding buffer pH 7.6 (38 µL; 10 mmol/L Tris-HCl, 2 mol/L NaCl, 1 mmol/L EDTA, 1 mL/L Tween 20) and diethylpyrocarbonate-treated water (30 µL; Ambion) in a 96-well microtiter plate for 15 min at 40 °C with mixing at 1400 rpm in a Thermomixer R (Eppendorf). The 96-well plate was subjected to vacuum filtration using the Vacuum Prep Tool (Pyrosequencing AB), which consists of a hand-grip with 96-filter probes. The bead-bound PCR products were transferred to the filter at the end of each 96-filter probe, and the remaining liquid was removed by vacuum filtration. The Vacuum Prep Tool was transferred to a trough containing 700 mL/L ethanol for 5 s. Nonbiotinylated strands were removed after treatment with denaturation solution (0.2 mol/L NaOH) for 5 s. Bead-bound biotinylated strands were retained on the filter and rinsed in washing buffer (10 mmol/L Tris-acetate, pH 7.6) for 7 s. The vacuum was removed to release beads from the filter. Bead-bound single-stranded DNA templates for each SNP were resuspended in annealing buffer (12 µL; 20 mM Tris-acetate, pH 7.6; 2 mM magnesium acetate) containing 0.5 µM (6 pmol) of the HPLC-purified Pyrosequencing primers R702W-F (5'-A CAT TGG AGA AGG CCC-3'), G908R-F (5'-CTT TTA ACC TTT TCA G-3'), and L1007fsinsC-F (5'-TGT CAT TCC TTT CAA GGG-3') in wells of a PSQ 96 plate (Pyrosequencing AB). To prevent extension from sequencing primers self-annealing, TG and AA mismatches (underlined and bolded above) were added to the R702W and G908R Pyrosequencing sequencing primers, respectively. Pyrosequencing assay design flexibility and chemistry made it possible to introduce these specific mismatches without affecting sequencing results. The mixture was heated at 94 °C for 2 min and then cooled to room temperature to facilitate annealing of the Pyrosequencing primers to templates. The plate was placed into the PSQ HS 96A System for analysis (Pyrosequencing AB).
An automated Pyrosequencing instrument, PSQ HS 96A (Pyrosequencing AB), was used to perform multiplex genotyping in a single Pyrosequencing reaction at 28 °C. Each well contained Sepharose beads with the bound single-stranded DNA template and its annealed Pyrosequencing primer. The Pyrosequencing reaction began with the addition of enzymes (DNA polymerase exonuclease-deficient, apyrase, luciferase, ATP sulfurylase) and substrates (adenosine 5'-phosphosulfate and luciferin) to each well. Individual dNTPs [TTP, dGTP, dCTP, or deoxyadenosine
-thio-triphosphate (dATP
S)] were sequentially added to the reactions after each elongation reaction step. Excess nucleoside triphosphates from the previous step were degraded by apyrase. Only one dNTP was added at each step according to a dispensation order (CGC TAT CGT GCG CTC TGA GA) entered in the Pyrosequencing SNP Multiplex Entry software (Pyrosequencing AB), which was selected to allow for simultaneous detection of all three SNPs. The sequence of incorporated nucleoside at each elongation step was inferred by measuring light emission as an indicator of incorporation with the particular dispensed dNTP. The resulting sequences were analyzed automatically by the SNP evaluation software, Ver. 1.1 (Pyrosequencing AB). Each SNP in all samples was confirmed by DNA sequence analysis with the ABI BigDye Terminator, Ver. 3.1, Cycle Sequencing Ready Reaction on a ABI 3100 DNA sequencer (Applied Biosystems) according to the manufacturers recommendations.
Three different primers pairs were used for multiplex PCR, and three different Pyrosequencing primers were used for simultaneous detection of the three SNPs in the CARD15 gene in a single Pyrosequencing reaction. For triplex CARD15 SNP genotyping, a dispensation order was determined that allowed each SNP variant to appear as individual peaks in the pyrogram. Pyrograms (left) and the theoretical Pyrosequencing results (right) for the three wild-type regions and six different genotypes are shown in panels A and BD of Fig. 1
, respectively. The contribution to the theoretical sequence readout at each nucleotide position by each of three different SNP regions is indicated in different colors on the right (blue for L1007fsinsC, red for G908R, and brown for R702W). Incorporation of each nucleotide for each allele produces a 0.5 peak height.
|
The blue boxes in panels A and B show results for the L1007fsinsC SNP extension reaction. L1007fsinsC is produced by a one-base insertion at nucleotide 3020 in exon 11 of the gene, which causes a frameshift mutation (20)(21). The wild-type genotype pyrogram for the L1007fsinsC is shown in Fig. 1A
, with the nucleotide dispensation order indicated above; the primer extension produces a double C peak height (e.g., 0.5 x 4 Cs; 2 for each wild-type allele) followed by no signals for G or C at the next two positions. The results for a heterozygote for G insertion at SNP L1007fsinsC are shown in the top portion of Fig. 1B
, with the primer extension producing a single C peak (0.5 x 2 Cs from the wild-type allele), a half-magnitude G peak (0.5 x 1 G from the variant SNP), and a single C peak (0.5 x 2 Cs from the variant SNP). The results for homozygosity for the G insertion are shown in the bottom portion of Fig. 1B
, with the primer extension producing no C peak, followed by a single G peak (corresponding to 0.5 x 2 Gs for the G insertion in both alleles) and a double C peak (0.5 x 4 Cs; 2 from each insertion G allele).
The red boxes in panels A and C of Fig. 1
show results for the G908R SNP extension reaction. G908R is caused by a G-to-C base substitution at nucleotide 2722 in exon 8 of the gene, producing a missense mutation (20)(21). Results for the wild-type genotype are shown in Fig. 1A
, with a quadruplicate G peak (0.5 x 4 Gs in each allele) and a triple C peak (0.5 x 2 Cs from the two wild-type alleles at this locus plus 0.5 x 2 Cs each from the two wild-type alleles at the L1007fsinsC and R702W loci), followed by no G and no C peaks (indicating no mutation). The bottom portion of Fig. 1C
shows the results for homozygosity for the G908R SNP, which produces a double G peak followed by a triple C (as indicated above), a single G, and a single C peak. The top portion in Fig. 1C
shows the results for heterozygosity at G908R with a triple G (0.5 x 4 Gs from the wild-type allele and 0.5 x 2 Gs from the variant) and C peaks (as indicated above), and half G and C peaks attributable to the mutation.
The peaks in the green boxes are results for extension reactions for wild-type, heterozygous, and mutant/mutant genotypes for R702W, a C-to-T base substitution at nucleotide 2104 in exon 4 of the gene, which produces a missense mutation (20)(21). Results are completely concordant with the theoretical expectation (Fig. 1, A and D
).
Allele and genotype frequencies were calculated from the families probands, and differences between patients and controls were tested by the Fisher exact and
2 tests, respectively. The frequency of L1007fsinsC was significantly higher in probands from the CD and mixed families than in controls, but not in probands from the UC families (P = 0.0002, 0.00002, and 0.38, respectively). In addition, the frequency of G908R was significantly higher in probands of CD families (P = 0.01). These frequencies are similar to those observed in other European populations (7). From the same data, we calculated that the relative risk of CD was 4.05 (95% confidence interval, 1.928.53) for heterozygotes and homozygotes for any of the CARD15 SNP variants (Table 1
). No Mendelian inheritance error was found. All 365 samples were scored correctly compared to automated DNA sequence analyses. Our experience demonstrates several advantages of Pyrosequencing (23)(24)(25), including a relatively easy and low-cost set up and a multiplexing capability that facilitates simultaneous detection of DNA variations that map to different sites in the same template. The use of nonproprietary supplies from the hardware manufacturers decreases the cost of Pyrosequencing, and template preparation can be standardized and robotically automated, allowing 96 samples to be sequenced in 25 min.
|
Acknowledgments
We thank the patients and their families who provided blood samples for this study and the Cooperative Investigators of the Italian Group for Inflammatory Bowel Disease (IG-IBD): G. Iaquinto and I. Sorrentini (U.O. di Gastroenterologia, Avellino); M. Campieri and P. Giochetti (Clinica Medica, Policlinico S. Orsola, Bologna); D. Valpiani and E. Ricci (Divisione di Medicina, Forlí); G. Frieri and M.T. Pimpo (Cattedra di Gastroenterologia, Universitá LAquila); S. Giaccari and L. Grasso (Divisione di Gastroenterologia, Galatina); M Saibene and R. De Franchis (Gastroenterologia, Ospedale Maggiore IRCCS, Milano); G. Bianchi Porro (Gastroenterologia, Ospedale Sacco, Milano); G. Riegler and F. Morace (Cattedra di Gastroenterologia, IIa Universitá di Napoli); S. Cucchiara and O. Borrelli (Dipartimento di Pediatria, IIa Universitá di Napoli); F. Castiglione and G. Mazzacca (Cattedra di Gastroenterologia, II Universitá Napoli); C. Prantera and A. Andreoli (Divisione di Gastroenterologia, Ospedale N. Regina Margherita, Roma); C. Papi and L. Capurso (Gastroenterologia, S. Filippo Neri, Roma); R. DIncá and G. Sturniolo (Cattedra di Gastroenterologia, Universitá di Padova); M. Oliva (Gastroenterologia, Ospedale "Cervello", Palermo); G. Lombardi, S. Fiorella, G. Corritore, and T. Latiano (Gastroenterologia, Ospedale CSS-IRCSS, San Giovanni Rotondo); M. Rizzetto, M. Astegiano, and F. Bresso (Dipartimento di Gastroenterologia, Ospedale Molinette, Torino); and A. Pera and M.T. Fiorentini (Divisione di Gastroenterologia, Ospedale Mauriziano, Torino). This work was supported by grants from the Center for Translational Medicine at Thomas Jefferson University (to P.F.), the Gastroenterology Research Unit at IRCCS Casa Sollievo della Sofferenza (to V.A. and A.L.), and the Italian Ministry of Health (RC2003, to B.D.).
References
The following articles in journals at HighWire Press have cited this article:
![]() |
C. A. A. van der Graaf, M. G. Netea, B. Franke, S. E. Girardin, J. W. M. van der Meer, and B. J. Kullberg Nucleotide Oligomerization Domain 2 (Nod2) Is Not Involved in the Pattern Recognition of Candida albicans. Clin. Vaccine Immunol., March 1, 2006; 13(3): 423 - 425. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. L. Roberts, R. B. Gearry, M. D.E. Allington, H. R. Morrin, B. A. Robinson, and F. A. Frizelle Caspase recruitment domain-containing protein 15 mutations in patients with colorectal cancer. Cancer Res., March 1, 2006; 66(5): 2532 - 2535. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Verri, F. Focher, G. Tettamanti, and V. Grazioli Two-Step Genetic Screening of Thrombophilia by Pyrosequencing Clin. Chem., July 1, 2005; 51(7): 1282 - 1284. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |