- Technical Note
- Open Access
Three multiplex snapshot assays for SNP genotyping in candidate innate immune genes
BMC Research Notesvolume 6, Article number: 54 (2013)
Innate immune system is the first line of research when studying immune response to diverse infections and autoimmune/inflammatory diseases. This immune response has been reported to be genetically diverse, due to polymorphisms coded by different genes. For this reason, our purpose was to develop a multiplex assay that allows the genotyping of candidate single nucleotide polymorphisms (SNPs) in innate immune genes.
We developed three multiplex PCR panels coupled with the minisequencing (SNaPshot) technique (multiplex PCR, multiplex primer extension, and capillary electrophoresis). The panels were tested in a sample set composed of 100 anonymous DNAs from healthy blood donors living in São Miguel Island (Azores, Portugal). Sixteen relevant SNPs among nine genes of the innate immune system – IL1α, IL1β, IL6, IL10, IL12RB1, TLR2, TLR4, TLR9 and CD14 – were genotyped and validated by direct sequencing, with the exception of one that was undetected by minisequencing. We suggest that these panels can be used in future studies for detection of risk gene variants in several populations and/or diseases.
In summary, we propose a multiplex assay that is able to identify the most frequent candidate SNPs in innate immune genes, using a medium scale genotyping platform. The assays can be used to evaluate the risk gene variants in populations of various geographic origins.
The innate immune system is involved in autoimmune/ inflammatory diseases, and is the first line of host defence against invading organisms. In part, most of these infective agents are successful in surviving, disseminating and infecting hosts, because of their ability to evade the innate immune system. Several studies suggested that, in human susceptibility to infection, host genetic makeup has a major role in innate immune response. This susceptibility has been reported among infectious diseases such as tuberculosis, bacterial meningitis, leprosy and hepatitis[5, 6]. The genetic variation of innate immune response to pathogens usually involves a variety of polymorphisms among different genes, including those that encode toll-like receptors (TLRs) and interleukins (ILs).
The methods used to investigate genetic variation in TLR s and IL s are mostly singleplex PCR-based techniques, such as direct sequencing and real-time PCR. For small and medium-size laboratories, singleplex techniques became expensive and time-consuming when a large number of samples are evaluated. The alternatives to these methodologies, such as minisequencing assay – SNaPshot multiplex PCR –, have been successfully applied to support studies regarding oncogenic related genes[7–9], bacterial serotyping and infectious diseases. Moreover, the SNaPshot is a relatively simple and affordable method that allows more efficient and quick results when it is compared with singleplex techniques. In addition, it only requires the use of a thermal cycler and a capillary electrophoresis, which are equipments commonly available in molecular biology laboratories. In order to apply this technique to candidate genes studies associated with autoimmune and infectious diseases, we developed a multiplex assay designed to genotype 16 single nucleotide polymorphisms (SNPs) among nine genes of the innate immune system: IL1α, IL1β, IL6, IL10, IL12RB1, TLR2, TLR4, TLR9 and CD14.
Selection of candidate innate immune system genes
Candidate genes and SNPs were selected considering their association to human autoimmune and infectious diseases. For this purpose, we used an updated and publically downloadable database – the Human Genome Epidemiology (HuGE) Literature Finder (http://www.hugenavigator.net) – that provides access to human genome epidemiology, including information on population prevalence of genetic variants, gene-disease associations, gene-gene and gene-environment interactions, and evaluation of genetic tests.
Population samples and sequence analysis
The sample set was composed of DNA from 100 healthy blood donors living in São Miguel Island, obtained from the anonymized Azorean DNA bank located at the Hospital of Divino Espirito Santo of Ponta Delgada, EPE, of the same island. This DNA bank was established after approval by the Ethics Committee for Health of the same Hospital, and follows the international ethical guidelines, which include informed consent, confidentiality, anonymity of personal data and abandonment option in case of expressed will.
From this sample set, 25 DNAs that evidence high SNP heterozygosity were selected to perform direct sequence analysis. Each 20 μl PCR reaction contained 100 ng of genomic DNA, 10 μM primers, 100 μM dNTPs (Promega), 25 nM MgCl2 (Qiagen), 1× Q-Solution (Qiagen), 1x buffer (Qiagen), 5 U of HotStart Taq (Qiagen) and sterile water. The PCR started with an enzyme activation step at 94°C for 15 min, then 40 cycles at 94°C for 1 min, 50°C for 1 min and 72°C for 1 min, followed by a final extension step at 72°C for 7 min. Amplification products were purified using the QIAquick PCR Purification Kit (Qiagen), according to the manufacturer’s instructions. Purified products were sequenced, using the same primers of the PCR amplification, with the BigDye Terminator v1.1 cycle sequencing kit (Applied Biosystems) under the following conditions: 1 μl ready reaction mix, 5 μl BigDye sequencing buffer, 3.2 pmol forward or reverse primer, 7 ng DNA and sterile water to a final reaction volume of 20 μl. Cycle sequencing was performed using initial denaturation step at 96°C for 1 min followed by 25 cycles at 96°C for 10 s, 50°C for 5 s and 60°C for 4 min in a GeneAmp® PCR System 2700 (Applied Biosystems). The sequencing products were purified with BigDye XTerminator® Purification Kit and separated by capillary electrophoresis in an automated sequencer (ABI 3130 Genetic Analizer, Applied Biosystems) with a 36 cm length capillary and POP-7™ polymer, according to the manufacturer’s instructions. Data were analyzed with Sequencing Analysis Software version 5.3.1 (Applied Biosystems). The alignment and edition of sequences were carried out using the Bioedit™ software version 7.0.0.
Multiplex PCR primers and SNaPshot probes
The first step of this study was to design a multiplex PCR assay that includes all 16 SNPs among the nine genes of the innate immune system (Table 1). The key point in this step is the combination of primer pairs in a master mix. We used the PrimerPlex biosoft software. This software is automatically linked to the SNP database (dbSNP, http://www.ncbi.nlm.nih.gov/projects/SNP/) and only requests the insert of SNP rs number. Using this information, the program picks the DNA sequences from the dbSNP and mixture them. The output of the program gives the best multiplex primers and probes pairs with the less secondary structures possible (Tables 2 and3). For this reason, it was not possible to use less PCR fragments to capture the SNPs: bigger fragments comprising more SNPs would allow the formation of more secondary structures. In the multiplex primer extension panels, poly-T tails of different lengths were added to each probe at their 5’ ends, to separate extension products by size.
Multiplex PCR amplification
Multiplex PCRs were performed in 10 μl total volume with Type-it® Mutation detect PCR kit (Qiagen) protocol. Thermal cycler conditions were: initial denaturation step at 95°C for 5 min, then 35 cycles at 95°C for 30 s, 50°C (panel 1) or 51°C (panels 2 and 3) for 90 s and 72°C for 30 s, followed by final extension at 68°C for 10 min. PCR products were treated with a mix of 5 U Exonuclease I (Exo I, Fermentas) and 1 U of FastAP™ Thermosensitive Alkaline Phosphatase (SAP, Fermentas) at 37°C for 15 min, and enzyme inactivation at 85°C for 15 min. PCR quality was evaluated after electrophoresis in a 3% agarose gel.
SNaPshot analysis was performed using a SNaPshot Multiplex Kit protocol (Applied Biosystems). Products were treated with 1 U of SAP at 37°C for 60 min and 75°C for 15 min, followed by a denaturation step at 95°C for 5 min. Detection was carried out using 0.5 μl of SNaPshot products mixed with 9 μl of HiDi™ formamide and 0.5 μl of GeneScan-120LIZ size standard (Applied Biosystems). Data were generated after capillary electrophoresis on an automated sequencer (ABI 3130 Genetic Analyzer, Applied Biosystems) with a 36 cm length capillary and POP-7™ polymer and analyzed with GeneMapper analysis software version 3.7 (Applied Biosystems).
Results and discussion
The strategy of the present work included the design of a multiplex assay, in order to detect 16 SNPs among nine genes of the innate immune system. For all panels, the specificity of each amplicon was evaluated first in singleplex and after in multiplex PCR, followed by electrophoresis in an agarose gel. According to the amplicons quality, primer concentrations and annealing temperatures were adjusted in order to obtain bands with equal intensities (Table 2). SNPs detection also implicated the optimisation of the three multiplex primer extension (SNaPshot) panels. In order to achieve more comparable peak heights, we adjusted probe concentrations (Table 3).
Figure1 represents the electropherograms obtained for each panel. Polymorphisms were identified on the basis of peak size and colour. As expected, minor shifts in the electrophoretic mobility were observed due to the incorporated base at the end of each probe and to the POP-7™ polymer. However, these shifts did not interfered with the analysis because poly-T tails increased probes spacing. Some samples evidenced one PCR product peak with more than one colour due to pull-ups, as illustrated in Figure1A, IL1-β polymorphism -511A>G. The peak in blue produced a secondary peak in black or green. This problem is probably due to a bleed-through when an adjacent lane has very high peak heights. Nevertheless, it was resolved by reducing the concentration of the respective probe.
The panels, here described, were evaluated in 100 DNA samples. From the 16 SNPs, we expected to identify 32 alleles and 48 different genotypes. However, we did not identify the G allele of the IL12 c.705 A>G (the last SNP in panel 1) and its inherent genotypes (AG and GG). This missing allele is a failure of the minisequencing reaction, since the Sanger DNA sequencing (validation test) revealed the presence of the two genotypes in our sample (Figures 2A and2B). Indeed, in a study carried out in the Brazilian population, the multiplex protocol, which was optimized for simultaneous detection of 12 autosomal SNPs, revealed that large amplicons (413 and 576 bp) had lower sensitivity to allele detection in multiplex minisequencing. Moreover, two homozygous genotypes were not observed in TLR2, namely: AA for c.2259A>G, and TT for c.2029C>T (panel 2 and 3, respectively). According to the SNP database (dbSNP), these genotypes were not found in the 11 HapMap populations (phase 3), which includes the Europeans (CEU). In fact, the A allele (c.2259A>G) is rare – 0.052 in CEU – and was only detected in homozygosity in tuberculosis patients; the T allele (c.2029C>T) is observed in Asian and African populations, but it seems to be absent in white populations[17, 18]. Nevertheless, in the present work, the T allele was identified in heterozygosity in one individual, as illustrated in Figure1C (c.2029C>T) and was confirmed by sequencing (Figure2C). This finding is not surprising since our previous studies show African influence in São Miguel population[19, 20].
In order to validate our assays, we selected 25 DNA samples for Sanger DNA sequencing. The results confirmed the 44 genotypes previously obtained by multiplex primer extension assay, as well as the two corresponding genotypes (AG and GG) of the missing allele IL12 c.705G.
Considering that these nine candidate innate immune genes have been investigated in autoimmune and infectious diseases through association studies, the three panels here described may allow the detection of at least 15 risk gene variants in all populations. Panel 1, for example, contains 4 SNPs in interleukin genes, which are often related to autoinflammatory diseases. On the other hand, panel 2 and 3 could be used for polymorphism genotyping (5 and 6 SNPs, respectively) in studies of infectious diseases associated with microbial pathogens, whose recognition by the innate immune system depends on Toll-like receptors.
With this study, we were able to develop combined multiplex PCR amplifications and SNaPshot reactions that genotypes 15 SNPs among nine innate immune response genes. These multiplex panels can be useful to investigate candidate immune response genes, using a medium scale genotyping platform, rather than singleplex PCR-based techniques.
Hill AV: Evolution, revolution and heresy in the genetics of infectious disease susceptibility. Philos Trans R Soc Lond B Biol Sci. 2012, 367: 840-849. 10.1098/rstb.2011.0275.
Tobin DM, Vary JC, Ray JP, Walsh GS, Dunstan SJ, Bang ND, Hagge DA, Khadge S, King MC, Hawn TR: The lta4h locus modulates susceptibility to mycobacterial infection in zebrafish and humans. Cell. 2010, 140: 717-730. 10.1016/j.cell.2010.02.013.
van Well GTJ, Sanders MS, Ouburg S, van Furth AM, Morre SA: Polymorphisms in Toll-Like Receptors 2, 4, and 9 are highly associated with hearing loss in survivors of bacterial meningitis. PLoS One. 2012, 7: e35837-10.1371/journal.pone.0035837.
Mira MT, Alcaïs A, Nguyen VT, Moraes MO, Di Flumeri C, Vu HT, Mai CP, Nguyen TH, Nguyen NB, Pham XK: Susceptibility to leprosy is associated with PARK2 and PACRG. Nature. 2004, 427: 636-640. 10.1038/nature02326.
Mbarek H, Ochi H, Urabe Y, Kumar V, Kubo M, Hosono N, Takahashi A, Kamatani Y, Miki D, Abe H: A genome-wide association study of chronic hepatitis B identified novel risk locus in a Japanese population. Hum Mol Genet. 2011, 20: 3884-3892. 10.1093/hmg/ddr301.
Thomas DL, Thio CL, Martin MP, Qi Y, Ge D, O'Huigin C, Kidd J, Kidd K, Khakoo SI, Alexander G: Genetic variation in IL28B and spontaneous clearance of hepatitis C virus. Nature. 2009, 461: U752-U798.
Hurst CD, Zuiverloon TC, Hafner C, Zwarthoff EC, Knowles MA: A SNaPshot assay for the rapid and simple detection of four common hotspot codon mutations in the PIK3CA gene. BMC Res Notes. 2009, 2: 66-10.1186/1756-0500-2-66.
Lurkin I, Stoehr R, Hurst CD, van Tilborg AAG, Knowles MA: Two Multiplex Assays that simultaneously identify 22 possible mutation sites in the KRAS, BRAF, NRAS and PIK3CA genes. PLoS One. 2010, 5: e8802-10.1371/journal.pone.0008802.
Filippini S, Blanco A, Fernandez-Marmiesse A, Alvarez-Iglesias V, Ruiz-Ponte C, Carracedo A, Vega A: Multiplex SNaPshot for detection of BRCA1/2 common mutations in Spanish and Spanish related breast/ovarian cancer families. BMC Med Genet. 2007, 8: 40-
Ben-Darif E, Jury F, De Pinna E, Threlfall EJ, Bolton FJ, Fox AJ, Upton M: Development of a multiplex primer extension assay for rapid detection of Salmonella isolates of diverse serotypes. J Clin Microbiol. 2010, 48: 1055-1060. 10.1128/JCM.01566-09.
Bouakaze C, Keyser C, de Martino SJ, Sougakoff W, Veziris N, Dabernat H, Ludes B: Identification and genotyping of Mycobacterium tuberculosis complex species by use of a SNaPshot minisequencing-based assay. J Clin Microbiol. 2010, 48: 1758-1766. 10.1128/JCM.02255-09.
Yu W, Gwinn M, Clyne M, Yesupriya A, Khoury MJ: A Navigator for Human Genome Epidemiology. Nat Genet. 2008, 40 (2): 124-5. 10.1038/ng0208-124.
Mota-Vieira L, Pacheco PR, Almeida ML, Cabral R, Carvalho J, Branco CC, de Fez L, Peixoto BR, Araújo AL, Mendonça P: Human DNA bank in São Miguel island (Azores), a resource for genetic diversity studies. Progress in Forensic Genetics, Proceedings of the 21st International Congress of Forensic Genetics, 13 - 17 September Ponta Delgada, Volume 1288. Edited by: Amorim A, Côrte-Real F, Morling N. 2005, 388-390.
PREMIER Biosoft software.http://www.premierbiosoft.com/primerplex/index.html,
Carvalho CM, Pena SD: Optimization of a multiplex minisequencing protocol for population studies and medical genetics. Genet Mol Res. 2005, 4: 115-125.
Ogus AC, Yoldas B, Ozdemir T, Uguz A, Olcen S, Keser I, Coskun M, Cilli A, Yegin O: The Arg753Gln polymorphism of the human toll-like receptor 2 gene in tuberculosis disease. Eur Respir J. 2004, 23: 219-223. 10.1183/09031936.03.00061703.
Ioana M, Ferwerda B, Plantinga TS, Stappers M, Oosting M, McCall M, Cimpoeru A, Burada F, Panduru N, Sauerwein R: Different patterns of Toll-like receptor 2 polymorphisms in populations of various ethnic and geographic origins. Infect Immun. 2012, 80: 1917-1922. 10.1128/IAI.00121-12.
Thurow HS, Sarturi CR, Fallavena PR, Paludo FJ, Picanço JB, Fraga LR, Graebin P, de Souza VC, Dias FS, Nóbrega Ode T: Very low frequencies of Toll-like receptor 2 supposed-2029T and 2258A (RS5743708) mutant alleles in southern Brazilian critically ill patients: would it be a lack of worldwide-accepted clinical applications of Toll-like receptor 2 variants?. Genet Test Mol Biomarkers. 2010, 14: 405-419. 10.1089/gtmb.2009.0169.
Pacheco PR, Branco CC, Cabral R, Costa S, Araújo AL, Peixoto BR, Mendonça P, Mota-Vieira L: The Y-chromosomal heritage of the Azores Islands population. Ann Hum Genet. 2005, 69: 145-156. 10.1046/j.1469-1809.2004.00147.x.
Pacheco PR, Branco CC, Gomes CT, Cabral R, Mota-Vieira L: HLA Class I and II profiles in São Miguel Island (Azores): genetic diversity and linkage disequilibrium. BMC Res Notes. 2010, 3: 134-10.1186/1756-0500-3-134.
Zenewicz LA, Abraham C, Flavell RA, Cho JH: Unravelling the genetics of autoimmunity. Cell. 2010, 140: 791-797. 10.1016/j.cell.2010.03.003.
Fang J, Fang D, Silver PB, Wen F, Li B, Ren X, Lin Q, Caspi RR, Su SB: The Role of TLR2, TRL3, TRL4, and TRL9 signalling in the pathogenesis of autoimmune disease in a retinal autoimmunity model. Invest Ophthalmol Vis Sci. 2010, 51: 3092-3099. 10.1167/iovs.09-4754.
This work was supported the “Fundação Calouste Gulbenkian” (grant reference number P-99888: SDH49), and “Direcção Regional da Ciência, Tecnologia e Comunicações (DRCTC)” of Azores Government. The authors wish to thank Drs Claudia Castelo Branco and Renato Pires for useful comments on the article.
The authors declare that they have no competing interests.
LME developed and performed the multiplex experiments and drafted the manuscript. SMB and MJB did the validation assays by Sanger DNA sequencing. LMV provided scientific orientation and revised the manuscript. All authors read and approved the final manuscript.
Sara M Bulhões, Maria J Brilhante contributed equally to this work.