Determining SNP allele frequencies in DNA pools

Single nucleotide polymorphisms (SNPs) are among the most common types of polymorphism used for genetic association studies. A method to allow the accurate quantitation of their allele frequencies from DNA pools would both increase throughput and decrease costs for large-scale genotyping. However, to date, most DNA pooling studies have concentrated on the use of microsatellite polymorphisms. In the case of SNPs that are restriction fragment length polymorphisms (RFLPs), studies have tended to use methods for the quantitation of allele frequency from pools that rely on densitometric evaluation of bands on an autoradiograph. Radiation-based methods have well-known drawbacks, and we present two alternative methods for the determination of SNP allele frequencies. For RFLPs, we used agarose gel electrophoresis of digested PCR products with ethidium bromide staining combined with densitometric analysis of gel images on a PC. For all types of SNP, we used allele-specific fluorescent probes in the Taqman assay to determine the relative frequencies of two different alleles. Both methods gave accurate and reproducible results, suggesting they are suitable for use in DNA pooling experiments.     

An Optimized Pentaplex PCR for Detecting DNA Mismatch Repair-Deficient Colorectal Cancers

Purpose

Microsatellite instability (MSI) is used to screen colorectal cancers (CRC) for Lynch Syndrome, and to predict outcome and response to treatment. The current technique for measuring MSI requires DNA from normal and neoplastic tissues, and fails to identify tumors with specific DNA mismatch repair (MMR) defects. We tested a panel of five quasi-monomorphic mononucleotide repeat markers amplified in a single multiplex PCR reaction (pentaplex PCR) to detect MSI.

Experimental Design

We investigated a cohort of 213 CRC patients, comprised of 114 MMR-deficient and 99 MMR-proficient tumors. Immunohistochemical (IHC) analysis evaluated the expression of MLH1, MSH2, PMS2 and MSH6. MSI status was defined by differences in the quasi-monomorphic variation range (QMVR) from a pool of normal DNA samples, and measuring differences in allele lengths in tumor DNA.

Results

Amplification of 426 normal alleles allowed optimization of the QMVR at each marker, and eliminated the requirement for matched reference DNA to define MSI in each sample. Using ≥2/5 unstable markers as the criteria for MSI resulted in a sensitivity of 95.6% (95% CI = 90.1–98.1%) and a positive predictive value of 100% (95% CI = 96.6%–100%). Detection of MSH6-deficiency was limited using all techniques. Data analysis with a three-marker panel (BAT26, NR21 and NR27) was comparable in sensitivity (97.4%) and positive predictive value (96.5%) to the five marker panel. Both approaches were superior to the standard approach to measuring MSI.

Conclusions

An optimized pentaplex (or triplex) PCR offers a facile, robust, very inexpensive, highly sensitive, and specific assay for the identification of MSI in CRC.     

Single nucleotide polymorphism (SNP) allele frequency estimation in DNA pools using Pyrosequencing

Identifying the genetic variation underlying complex disease requires analysis of many single nucleotide polymorphisms (SNPs) in a large number of samples. Several high-throughput SNP genotyping techniques are available; however, their cost promotes the use of association screening with pooled DNA. This protocol describes the estimation of SNP allele frequencies in pools of DNA using the quantitative sequencing method Pyrosequencing (PSQ). PSQ is a relatively recently described high-throughput method for genotyping, allele frequency estimation and DNA methylation analysis based on the detection of real-time pyrophosphate release during synthesis of the complementary strand to a PCR product. The protocol involves the following steps: (i) quantity and quality assessment of individual DNA samples; (ii) DNA pooling, which may be undertaken at the pre- or post-PCR stage; (iii) PCR amplification of PSQ template containing the variable sequence region of interest; and (iv) PSQ to determine the frequency of alleles at a particular SNP site. Once the quantity and quality of individual DNA samples has been assessed, the protocol usually requires a few days for setting up pre-PCR pools, depending on sample number. After PCR amplification, preparation and analysis of PCR amplicon by PSQ takes 1 h per plate.     

Design of a novel quantitative PCR (QPCR)-based protocol for genotyping mice carrying the neuroprotective Wallerian degeneration slow (Wld s ) gene

Background

The strain of MeCP2^tm1.1Bird mice is a broadly used model for Rett syndrome. Because males carrying the invalidated MeCP2 locus are sterile, this strain has to be maintained in a heterozygous state. All animals therefore have to be genotyped at every generation to discriminate those carrying the invalidated allele (+/- females and y/- males) from those that do not. This is conveniently carried out by PCR on tail genomic DNA but because the primer pairs described initially for this purpose yield very similar size DNA bands on the WT and the KO alleles, this requires to carry out two independent PCR reactions on tail DNA preparations from all animals.

Results

After cloning and sequencing the PCR fragment amplified on the KO allele, we tested several sets of primers that were designed to yield PCR fragments of different sizes on the KO and WT alleles.

Conclusion

We have thus identified a set of three primers that allows for efficient genotyping of the animals by a single PCR reaction. Furthermore, using of this set of primers also resolves a recurrent problem related to the tendency of one of the initial primers to give rise to a non specific band because of its capacity to anneal at both ends of a repeated genomic element which we have identified as MurvyLTR.     

Estimation of allele frequency in pooled DNA by using PCR–RFLP combined with microchip electrophoresis

Biallelic marker, most commonly single nucleotide polymorphism (SNP), is widely utilized in genetic association analysis, which can be speeded up by estimating allele frequency in pooled DNA instead of individual genotyping. Several methods have shown high accuracy and precision for allele frequency estimation in pools. Here, we explored PCR restriction fragment length polymorphism (PCR–RFLP) combined with microchip electrophoresis as a possible strategy for allele frequency estimation in DNA pools. We have used the commercial available Agilent 2100 microchip electrophoresis analysis system for quantifying the enzymatically digested DNA fragments and the fluorescence intensities to estimate the allele frequencies in the DNA pools. In this study, we have estimated the allele frequencies of five SNPs in a DNA pool composed of 141 previously genotyped health controls and a DNA pool composed of 96 previously genotyped gastric cancer patients with a frequency representation of 10–90% for the variant allele. Our studies show that accurate, quantitative data on allele frequencies, suitable for investigating the association of SNPs with complex disorders, can be estimated from pooled DNA samples by using this assay. This approach, being independent of the number of samples, promises to drastically reduce the labor and cost of genotyping in the initial association analysis.     

Different real time PCR approaches for the fine quantification of SNP's alleles in DNA pools: assays development, characterization and pre-validation

Single nucleotide polymorphisms (SNPs) are becoming the most common type of markers used in genetic analysis. In the present report a SNP has been chosen to test the applicability of Real Time PCR to discriminate and quantify SNPs alleles on DNA pools. Amplification Refractory Mutation System (ARMS) and Mismatch Amplification Mutation Assay (MAMA) has been applied. Each assay has been pre-validated testing specificity and performances (linearity, PCR efficiency, interference limit, limit of detection, limit of quantification, precision and accuracy). Both the approaches achieve a precise and accurate estimation of the allele frequencies on pooled DNA samples in the range from 5 % to 95 % and don't require standard curves or calibrators. The lowest measurement that could be significantly distinguished from the background noise has been determined around the 1 % for both the approaches, allowing to extend the range of quantifications from 1 % to 99 %. Furthermore applicability of Real Time PCR assays for general diagnostic purposes is discussed.     

高频度微卫星不稳定性大肠癌"修饰型"微卫星变化与MLH1及KRAS基因突变密切相关

本研究通过方法学的改良和观察方式的创新试图阐明这种现象的原因.微卫星非传统的检测方法仅能实现微卫星定性检测,我所在的研究组开发了自动片段分析双荧光标识技术,提高了微卫星检测的感度和重复性,并实现了微卫星片段变化长度的定量.小于6碱基的微卫星变化被定义为修饰型微卫星不稳定,大于8碱基的变化被定义为跳跃型微卫星不稳定,它们的电泳谱截然不同.前者表现为在非肿瘤来源微卫星位点基础上的增加或减少,后者表现为距离非肿瘤微卫星片段远隔部位的新波形的出现.通过研究我们发现,在DNA错配修复缺陷细胞系及基因敲除大鼠自发肿瘤样本,仅有修饰型微卫星不稳定性检出;在人类DNA错配修复缺陷细胞系连续80次传代也没有检出跳跃型变化.跳跃型变化不能通过简单重复序列不稳定基础上的增加或减少的累加而获得.在76例散发大肠癌,我们检测了微卫星不稳定性,KRAS基因突变,并对高频度微卫星不稳定性病例的两个主要DNA错配修复基因MSH2和MLHl进行了全长测序.我们发现,在大肠癌,按频度的传统分类与按波形变化的分类有高度的一致性,高频度微卫星不稳定性病例均检测到跳跃型表现,低频度微卫星不稳定性都表现为修饰型变化.在12例高频度微卫星不稳定病例,有三例检出了跳跃型和修饰型同时存在微卫星不稳定的特殊表型,这3例均检出KRAS的突变,更有趣的是该3例病例也同时检出了DNA错配修复基因MLH1的变异.而在其他9例高频度微卫星不稳定病例,KRAS突变及MLH1、MSH2交变未检出.通过对突变谱的分析我们还发现,修饰型微卫星不稳定与KTAS基因12号密码子的转换型突变高度相关,而微卫星稳定的病例检出的KRAS基因12号密码子突变多为颠换型突变.修饰型微卫星不稳定表型检出的高频度转换突变可由DNA错配修复缺陷的分子背景解释.通过本研究,我们认为以波形为基础的微卫星不稳定新分型可能是解决目前微卫星研究领域矛盾的一个选项.一直公认为高频度微卫星不稳定性是"真正"的DNA错配修复缺陷表型,我们的研究提示实际上高频度微卫星的可能是多元的.修饰型微卫星不稳定与DNA错配修复缺陷直接关联,而跳跃型微卫星不稳定的原因尚未阐明.在高频度为微型不稳定中,携带修饰型变化的病例可以通过DNA错配修复系统缺陷来解释其病因.     

成骨细胞特异性表达Cre重组酶转基因小鼠的建立

构建了含有骨钙素基因启动子、Cre重组酶基因和人生长激素基因polyA的转基因载体pOC-Cre, 以显微注射的方法将4.6 kb的转基因片段OC-Cre导入小鼠受精卵。16只子代小鼠中经PCR和Southern杂交鉴定, 有2只小鼠携带外源基因, 整合率为12.5%。为了检测OC-Cre在转基因小鼠中表达的组织特异性, 将转基因首建者小鼠与基因组上携带有LoxP位点的条件性Smad4基因敲除小鼠交配, PCR结果显示, 仅在子代纯合型小鼠骨组织基因组中扩增出了Cre介导重组后的片段。将OC-Cre转基因小鼠与ROSA26报告小鼠交配, 利用LacZ染色对双转基因阳性子代小鼠进行检测, 结果显示Cre重组酶在成骨细胞中特异性表达并介导ROSA基因座LoxP位点间的重组。所有这些结果说明：所建立的OC-Cre转基因小鼠在成骨细胞中特异性表达Cre重组酶, 并能在体内介导成骨细胞基因组上LoxP位点间的重组, 是一种理想的研制成骨细胞特异性基因敲除小鼠的工具小鼠。     

Ribonucleoside and deoxyribonucleoside triphosphate pools during 2-aminopurine mutagenesis in T4 mutator-, wild type-, and antimutator-infected Escherichia coli

Ribonucleoside and deoxyribonucleoside triphosphate pools have been measured in Escherichia coli infected with bacteriophage T4 DNA polymerase mutator, wild type, and antimutator alleles during mutagenesis by the base analogue 2-aminopurine. ATP and GTP pools expand significantly during mutagenesis, while CTP and UTP pools contract slightly. The DNA polymerase (gene 43) alleles and an rII lesion perturb normal dNTP pools more than does the presence of 2-aminopurine. We find no evidence that 2-aminopurine induces mutations indirectly by causing an imbalance in normal dNTP pools. Rather, it seems likely that, by forming base mispairs with thymine and with cytosine, 2-aminopurine is involved directly in causing bidirectional A.T in equilibrium G.C transitions. The ratios for 2-aminopurine deoxyribonucleoside triphosphate/dATP pools are 5-8% for tsL56 mutator and 1-5% for tsL141 antimutator and 43+ alleles. We conclude that the significant differences observed in the frequencies of induced transition mutations in the three alleles can be attributed primarily to the properties of the DNA polymerases with their associated 3'-exonuclease activities in controlling the frequency of 2-aminopurine.cystosine base mispairs.     

MethyLight droplet digital PCR for detection and absolute quantification of infrequently methylated alleles

Aberrant DNA methylation is a common epigenetic alteration found in colorectal adenomas and cancers and plays a role in cancer initiation and progression. Aberrantly methylated DNA loci can also be found infrequently present in normal colon tissue, where they seem to have potential to be used as colorectal cancer (CRC) risk biomarkers. However, detection and precise quantification of the infrequent methylation events seen in normal colon is likely beyond the capability of commonly used PCR technologies. To determine the potential for methylated DNA loci as CRC risk biomarkers, we developed MethyLight droplet digital PCR (ddPCR) assays and compared their performance to the widely used conventional MethyLight PCR. Our analyses demonstrated the capacity of MethyLight ddPCR to detect a single methylated NTRK3 allele from among more than 3125 unmethylated alleles, 25-fold more sensitive than conventional MethyLight PCR. The MethyLight ddPCR assay detected as little as 19 and 38 haploid genome equivalents of methylated EVL and methylated NTRK3, respectively, which far exceeded conventional MethyLight PCR (379 haploid genome equivalents for both genes). When assessing methylated EVL levels in CRC tissue samples, MethyLight ddPCR reduced coefficients of variation (CV) to 6–65% of CVs seen with conventional MethyLight PCR. Importantly, we showed the ability of MethyLight ddPCR to detect infrequently methylated EVL alleles in normal colon mucosa samples that could not be detected by conventional MethyLight PCR. This study suggests that the sensitivity and precision of methylation detection by MethyLight ddPCR enhances the potential of methylated alleles for use as CRC risk biomarkers.     

Evaluation of DNA pooling for the estimation of microsatellite allele frequencies: a case study using striped bass (Morone saxatilis)

Using striped bass (Morone saxatilis) and six multiplexed microsatellite markers, we evaluated procedures for estimating allele frequencies by pooling DNA from multiple individuals, a method suggested as cost-effective relative to individual genotyping. Using moment-based estimators, we estimated allele frequencies in experimental DNA pools and found that the three primary laboratory steps, DNA quantitation and pooling, PCR amplification, and electrophoresis, accounted for 23, 48, and 29%, respectively, of the technical variance of estimates in pools containing DNA from 2-24 individuals. Exact allele-frequency estimates could be made for pools of sizes 2-8, depending on the locus, by using an integer-valued estimator. Larger pools of size 12 and 24 tended to yield biased estimates; however, replicates of these estimates detected allele frequency differences among pools with different allelic compositions. We also derive an unbiased estimator of Hardy-Weinberg disequilibrium coefficients that uses multiple DNA pools and analyze the cost-efficiency of DNA pooling. DNA pooling yields the most potential cost savings when a large number of loci are employed using a large number of individuals, a situation becoming increasingly common as microsatellite loci are developed in increasing numbers of taxa.     

高频度微卫星不稳定性大肠癌“修饰型”微卫星变化与MLH1及KRAS基因突变密切相关

本研究通过方法学的改良和观察方式的创新试图阐明这种现象的原因。微卫星非传统的检测方法仅能实现微卫星定性检测,我所在的研究组开发了自动片段分析双荧光标识技术,提高了微卫星检测的感度和重复性。并实现了微卫星片段变化长度的定量。小于6碱基的微卫星变化被定义为修饰型微卫星不稳定,大于8碱基的变化被定义为跳跃型微卫星不稳定．它们的电泳谱截然不同。前者表现为在非肿瘤来源微卫星位点基础上的增加或减少,后者表现为距离非肿瘤微卫星片段远隔部位的新波形的出现。通过研究我们发现,在DNA错配修复缺陷细胞系及基因敲除大鼠自发肿瘤样本,仅有修饰型微卫星不稳定性检出;在人类DNA错配修复缺陷细胞系连续80次传代也没有检出跳跃型变化。跳跃型变化不能通过简单重复序列不稳定基础上的增加或减少的累加而获得。在76例散发大肠癌,我们检测了微卫星不稳定性,KRAS基因突变,并对高频度微卫星不稳定性病例的两个主要DNA错配修复基因MSH2和MLH1进行了全长测序。我们发现,在大肠癌,按频度的传统分类与按波形变化的分类有高度的一致性,高频度微卫星不稳定性病例均检测到跳跃型表现,低频度微卫星不稳定性都表现为修饰型变化。在12例高频度微卫星不稳定病例,有三例检出了跳跃型和修饰型同时存在微卫星不稳定的特殊表型,这3例均检出KRAS的突变,更有趣的是该3例病例也同时检出了DNA错配修复基因MLH1的变异。而在其他9例高频度微卫星不稳定病例,KRAS突变及MLH1、MSH2突变未检出。通过对突变谱的分析我们还发现,修饰型微卫星不稳定与KRAS基因12号密码子的转换型突变高度相关,而微卫星稳定的病例检出的KRAS基因12号密码子突变多为颠换型突变。修饰型微卫星不稳定表型检出的高频度转换?     

MiRNA genes constitute new targets for microsatellite instability in colorectal cancer

Mismatch repair-deficient colorectal cancers (CRC) display widespread instability at DNA microsatellite sequences (MSI). Although MSI has been reported to commonly occur at coding repeats, leading to alterations in the function of a number of genes encoding cancer-related proteins, nothing is known about the putative impact of this process on non-coding microRNAs. In miRbase V15, we identified very few human microRNA genes with mono- or di-nucleotide repeats (n = 27). A mutational analysis of these sequences in a large series of MSI CRC cell lines and primary tumors underscored instability in 15 of the 24 microRNA genes successfully studied at variable frequencies ranging from 2.5% to 100%. Following a maximum likelihood statistical method, microRNA genes were separated into two groups that differed significantly in their mutation frequencies and in their tendency to represent mutations that may or may not be under selective pressures during MSI tumoral progression. The first group included 21 genes that displayed no or few mutations in CRC. The second group contained three genes, i.e., hsa-mir-1273c, hsa-mir-1303 and hsa-mir-567, with frequent (≥ 80%) and sometimes bi-allelic mutations in MSI tumors. For the only one expressed in colonic tissues, hsa-mir-1303, no direct link was found between the presence or not of mono- or bi-allelic alterations and the levels of mature miR expression in MSI cell lines, as determined by sequencing and quantitative PCR respectively. Overall, our results provide evidence that DNA repeats contained in human miRNA genes are relatively rare and preserved from mutations due to MSI in MMR-deficient cancer cells. Functional studies are now required to conclude whether mutated miRNAs, and especially the miR-1303, might have a role in MSI tumorigenesis.     

SNP genotyping on pooled DNAs: comparison of genotyping technologies and a semi automated method for data storage and analysis

We have compared the accuracy, efficiency and robustness of three methods of genotyping single nucleotide polymorphisms on pooled DNAs. We conclude that (i) the frequencies of the two alleles in pools should be corrected with a factor for unequal allelic amplification, which should be estimated from the mean ratio of a set of heterozygotes (k); (ii) the repeatability of an assay is more important than pinpoint accuracy when estimating allele frequencies, and assays should therefore be optimised to increase the repeatability; and (iii) the size of a pool has a relatively small effect on the accuracy of allele frequency estimation. We therefore recommend that large pools are genotyped and replicated a minimum of four times. In addition, we describe statistical approaches to allow rigorous comparison of DNA pool results. Finally, we describe an extension to our ACeDB database that facilitates management and analysis of the data generated by association studies.     

Contributions of different mosquito species to the transmission of lymphatic filariasis in central Nigeria: implications for monitoring infection by PCR in mosquito pools

Background

Monitoring and evaluation are essential to the successful implementation of mass drug administration programmes for LF elimination. Monitoring transmission when it is low requires both large numbers of mosquito vectors and sensitive methods for detecting Wuchereria bancrofti infections in them. PCR-based methods are preferred over classical dissections but the best protocol so far achieved detection of one L₃ Wuchereria bancrofti larva in a pool of 35–50 Anopheles mosquitoes. It also lacks consistency and remains still a costly tool. Hence we decided to improve upon this to achieve detection in a pool of 100 or more by enhancing the quality of the template DNA. Prior to this we also evaluated three vector sampling methods in the context of numbers for monitoring.

Methods

Human landing, pyrethrium spray and light traps catches were conducted concurrently at sites in an LF endemic district in Ghana and the numbers obtained compared. Two DNA extraction methods; Bender buffer and phenol/chloroform purification, and DNAeasy Tissue kit (Quaigen Inc) were used on pools of 25, 50, 75 100 and 150 mosquitoes each seeded with one L₃ or its quivalent amount of DNA. Then another set of extracted DNA by the two methods was subjected to Dynal bead purification method (using capture oligonucleotide primers). These were used as template DNA in PCR to amplify W. bancrofti sequences. The best PCR result was then evaluated in the field at five sites by comparing its results (infections per 1000 mosquitoes) with that of dissection of roughly equal samples sizes.

Results

The largest numbers of mosquitoes were obtained with the human landing catches at all the sites sampled. Although PCR detection of one L₃ in pools of 25, 50 and 75 mosquitoes was consistent irrespective of the extraction method, that of one L₃ in 100 was only achieved with the kit-extracted DNA/Dynal bead purification method. Infections were found at only two sites by both dissection and pool-screening being 14.3 and 19 versus 13.4 and 20.1 per 1000 Anopheles mosquitoes respectively, which were not statistically significant

Discussion and conclusion

HLC still remains the best option for sampling for the large numbers of mosquitoes required for monitoring transmission during MDA programmes, when vector population densities are high and classical indices of transmission are required. One – in – 100 detection is an improvement on previous PCR pool-screening methods, which in our opinion was a result of the introduction of the extra step of parasite DNA capture using Dynal/beads. As pool sizes increase the insects DNA will swamp parasite DNA making the latter less available for an efficient PCR, therefore we propose either additional steps of parasite DNA capture or real-time PCR to improve further the pool screening method. The study also attests also to the applicability of Katholi et al's algorithm developed for determining onchocerciasis prevalence in LF studies.     

Patterns of nuclear DNA degeneration over time--a case study in historic teeth samples

The amount of nuclear DNA extracted from teeth of 279 individual red fox Vulpes vulpes collected over a period spanning the last three decades was determined by quantitative polymerase chain reaction (PCR). Although teeth were autoclaved during initial collection, 73.8% of extracts contained sufficient DNA concentration (> 5 pg/ micro L) suitable for reliable microsatellite genotyping but the quantity of nuclear DNA decayed significantly over time in a nonlinear pattern. The success of PCR amplification across four examined canine microsatellites over time was dependent on fragment size. By including data from two different tests for human contamination and from frequencies of allelic dropout and false alleles, the methodological constraints of population genetic studies using microsatellite loci amplified from historic DNA are discussed.     

Rapid SNP allele frequency determination in genomic DNA pools by pyrosequencing

Individual genotyping of single nucleotide polymorphisms (SNPs) remains expensive, especially for linkage disequilibrium mapping strategies involving high-throughput SNP genotyping. On one hand, current methods may suit scientific and laboratory needs in regard to accuracy, reproducibility/robustness, and large-scale application. On the other hand, a cheaper and less time-consuming alternative to individual genotyping is the use of SNP allelefrequencies determined in DNA pools. We have developed an accurate and reproducible protocol for allele frequency determination using Pyrosequencing technology in large genomic DNA pools (374 individuals). The measured correlation (R2) in large DNA pools was 0.980. In the context of disease-associated SNPs studies, we compared the allele frequencies between the disease (e.g., type 2 diabetes and obesity) and control groups detected by either individual genotyping or Pyrosequencing of DNA pools. In large pools, the variation between the two methods was 1.5 +/- 0.9%. It may be concluded that the allele frequency determination protocol could reliably detect over 4% differences between populations. The method is economical in regard to amounts of DNA, PCR, and primer extension reagents required. Furthermore, it allows the rapid determination of allelefrequency differences in case/control groups for association studies and susceptibility gene discovery in complex diseases.     

Demonstration of sequence variations in the promoter region of the human cystatin C gene.

Four point mutations in the promoter region of the human cystatin C gene have been detected by direct sequencing of polymerase chain reaction (PCR) amplified DNA. The four base changes are all localized within a short segment of 85 base pairs. Three cystatin C gene alleles could be defined with respect to these promoter mutations; one with the sequence previously published, one carrying three of the mutations and one with all four base substitutions. Two of the observed mutations are involved in a novel Sst II polymorphism and another generates a new Dde I restriction site. A PCR-based assay for analysis of these Sst II and Dde I sites was designed and used to demonstrate Mendelian inheritance of the polymorphisms as well as to determine the frequencies of the cystatin C gene alleles in the population.     

Cost-efficient multiplex PCR for routine genotyping of up to nine classical HLA loci in a single analytical run of multiple samples by next generation sequencing

Background

HLA genotyping by next generation sequencing (NGS) requires three basic steps, PCR, NGS, and allele assignment. Compared to the conventional methods, such as PCR-sequence specific oligonucleotide primers (SSOP) and -sequence based typing (SBT), PCR-NGS is extremely labor intensive and time consuming. In order to simplify and accelerate the NGS-based HLA genotyping method for multiple DNA samples, we developed and evaluated four multiplex PCR methods for genotyping up to nine classical HLA loci including HLA-A, HLA-B, HLA-C, HLA-DRB1/3/4/5, HLA-DQB1, and HLA-DPB1.

Results

We developed multiplex PCR methods using newly and previously designed middle ranged PCR primer sets for genotyping different combinations of HLA loci, (1) HLA-DRB1/3/4/5, (2) HLA-DQB1 (3.8 kb to 5.3 kb), (3) HLA-A, HLA-B, HLA-C, and (4) HLA-DPB1 (4.6 kb to 7.2 kb). The primer sets were designed to genotype polymorphic exons to the field 3 level or 6-digit typing. When we evaluated the PCR method for genotyping all nine HLA loci (9LOCI) using 46 Japanese reference subjects who represented a distribution of more than 99.5% of the HLA alleles at each of the nine HLA loci, all of the 276 alleles genotyped, except for HLA-DRB3/4/5 alleles, were consistent with known alleles assigned by the conventional methods together with relevant locus balance and no excessive allelic imbalance. One multiplex PCR method (9LOCI) was able to provide precise genotyping data even when only 1 ng of genomic DNA was used for the PCR as a sample template.

Conclusions

In this study, we have demonstrated that the multiplex PCR approach for NGS-based HLA genotyping could serve as an alternative routine HLA genotyping method, possibly replacing the conventional methods by providing an accelerated yet robust amplification step. The method also could provide significant merits for clinical applications with its ability to amplify lower quantity of samples and the cost-saving factors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1514-4) contains supplementary material, which is available to authorized users.