Primer Design for Multiplex PCR Using a Genetic Algorithm

[[abstract]]Multiplex Polymerase chain reaction (PCR) is the term used when more than one pair of primers is used in a polymerase chain reaction. The goal of multiplex PCR is to amplify several segments of target DNA simultaneously and thereby to conserve template DNA, save time, and minimize expense. The success of the experiment is dependent on primer design. However, this can be a dreary task as there are many constrains such as melting temperatures, primer length, GC content and complementarity that need to be optimized to obtain a good PCR product. In our investigations, we found few primer design tools for multiplex PCR and there was no suitable tool for our partners who want to use a multiplex PCR genotypic assay. The tool draws on a genetic algorithm where stochastic approaches based on the concept of biological evolution, biological genetics and genetic operations on chromosomes are used to find an optimal solution for multiplex PCR. The presented experimental results indicate that the proposed algorithm is able to find a set of primer pairs that not only obey the design properties but also work in the same tube

Multiplex PCR primer design for gene family using genetic algorithm

The multiplex PCR experiment is to amplify multiple regions of a DNA sequence at the same time by using different primer pairs. Designing feasible primer pairs for multiplex PCR is a tedious task since there are too many constraints to be satisfied. In this paper, a new method for multiplex PCR primer design strategy using genetic algorithm is proposed. The proposed algorithm is able to find a set of suitable primer pairs more efficient and uses a MAP model to speed up the examination of the specificity constraint that is important for gene family sequences. The dry-dock experiment shows that the proposed algorithm finds several sets of primer pairs of gene family sequences for multiplex PCR that not only obey the design properties, but also have specificity

Diversity of the parB and repA genes of the Burkholderia cepacia complex and their utility for rapid identification of Burkholderia cenocepacia

Background: Burkholderia cenocepacia is the most prominent species of the B. cepacia complex (Bcc), a group of nine closely related and difficult to identify bacteria that cause serious infections in patients with cystic fibrosis. Despite its clinical relevance, identification of B. cenocepacia as a single species is unavailable, as it splits by a widely used recA gene-based PCR identification method into discrete phylogenetic subgroups IIIA, IIIB, IIIC and IIID. With the aim of identifying gene targets suitable for unified detection of B. cenocepacia strains, we examined sequence polymorphisms in the repA and parB genes. These essential genes are involved in the replication and partitioning of bacterial replicons, hence we also had the opportunity for the first time to investigate the evolution of the multireplicon (three chromosome) structure of Bcc genomes. Results: Alignment of the repA and parB genes from publicly available Bcc genome sequences enabled the design of primers for their amplification and sequence analysis. Multilocus sequencing typing, a highly discriminatory method for Bcc species and strain discrimination, was used to select strains of unique sequence types (STs) that spanned the known Bcc genetic diversity. Sequence datasets of repA (83 isolates, 67 STs) and parB (120 isolates, 95 STs) genes from the second chromosome were aligned and examined phylogenetically to identify polymorphisms suitable for identification of B. cenocepacia. In contrast to parB, the Bcc repA sequences demonstrated distinct clustering of B. cenocepacia from other species, which enabled the design a species-specific multiplex PCR. The novel single-reaction B. cenocepacia detection method was tested on a panel of 142 different Bcc strains (142 STs) and distinguished recA groups IIIA, IIIB and IIID, from all other Bcc members with 100% sensitivity and 93% specificity. Conclusion: The repA-based multiplex PCR is a useful aid to the rapid identification of the most clinically relevant B. cenocepacia recA subgroups IIIA, IIIB and IIID. Phylogenetic analysis of repA and parB genes demonstrated that acquisition of the second and third replicons of Bcc genomes occurred prior to their differentiation into discrete species and that the sharing of replicons across species had not occurred

Use of Whole Genome Phylogeny and Comparisons in the Development of a Multiplex-PCR Assay to Identify Sequence Type 36 Vibrio parahaemolyticus

Vibrio parahaemolyticus sequence type (ST) 36 strains that are native to the Pacific Ocean have recently caused multi-state outbreaks of gastroenteritis linked to shellfish harvested from the Atlantic Ocean. Whole genome comparisons of 295 genomes of V. parahaemolyticus, including several traced to northeastern US sources, were used to identify diagnostic loci: one putatively encoding an endonuclease (prp), and two others potentially conferring O-antigenic properties (cps and flp). The combination of all three loci was present only in one clade of closely-related strains, of ST36, ST59 and one additional unknown sequence type. However, each locus was also identified outside this clade, with prp and flp occurring in only two non-clade isolates, and cps in four. Based on the distribution of these loci in sequenced genomes, prp could identify clade strains with \u3e99% accuracy, but the addition of one more locus would increase accuracy to 100%. Oligonucleotide primers targeting prp and cps were combined in a multiplex PCR method that defines species using the tlh locus, and determines presence of both the tdh and trh hemolysin-encoding genes which are also present in ST36. Application of the method in vitro to a collection of 94 clinical isolates collected over a four year period in three Northeastern US, and 87 environmental isolates, revealed the prp and cps amplicons were only detected in clinical isolates identified as belonging to the ST36-clade, and in no environmental isolates from the region. The assay should improve detection and surveillance, thereby reducing infections

A new multi locus variable number of tandem repeat analysis scheme for epidemiological surveillance of Xanthomonas vasicola pv. musacearum, the plant pathogen causing bacterial wilt on banana and enset

Xanthomonas vasicola pv. musacearum (Xvm) which causes Xanthomonas wilt (XW) on banana (Musa accuminata x balbisiana) and enset (Ensete ventricosum), is closely related to the species Xanthomonas vasicola that contains the pathovars vasculorum (Xvv) and holcicola (Xvh), respectively pathogenic to sugarcane and sorghum. Xvm is considered a monomorphic bacterium whose intra-pathovar diversity remains poorly understood. With the sudden emergence of Xvm within east and central Africa coupled with the unknown origin of one of the two sublineages suggested for Xvm, attention has shifted to adapting technologies that focus on identifying the origin and distribution of the genetic diversity within this pathogen. Although microbiological and conventional molecular diagnostics have been useful in pathogen identification. Recent advances have ushered in an era of genomic epidemiology that aids in characterizing monomorphic pathogens. To unravel the origin and pathways of the recent emergence of XW in Eastern and Central Africa, there was a need for a genotyping tool adapted for molecular epidemiology. Multi-Locus Variable Number of Tandem Repeat Analysis (MLVA) is able to resolve the evolutionary patterns and invasion routes of a pathogen. In this study, we identified microsatellite loci from nine published Xvm genome sequences. Of the 36 detected microsatellite loci, 21 were selected for primer design and 19 determined to be highly typeable, specific, reproducible and polymorphic with two- to four- alleles per locus on a sub-collection. The 19 markers were multiplexed and applied to genotype 335 Xvm strains isolated from seven countries over several years. The microsatellite markers grouped the Xvm collection into three clusters; with two similar to the SNP-based sublineages 1 and 2 and a new cluster 3, revealing an unknown diversity in Ethiopia. Five of the 19 markers had alleles present in both Xvm and Xanthomonas vasicola pathovars holcicola and vasculorum, supporting the phylogenetic closeliness of these three pathovars. Thank to the public availability of the haplotypes on the MLVABank database, this highly reliable and polymorphic genotyping tool can be further used in a transnational surveillance network to monitor the spread and evolution of XW throughout Africa.. It will inform and guide management of Xvm both in banana-based and enset-based cropping systems. Due to the suitability of MLVA-19 markers for population genetic analyses, this genotyping tool will also be used in future microevolution studies

Two new rapid SNP-typing methods for classifying Mycobacterium tuberculosis complex into the main phylogenetic lineages

There is increasing evidence that strain variation in Mycobacterium tuberculosis complex (MTBC) might influence the outcome of tuberculosis infection and disease. To assess genotype-phenotype associations, phylogenetically robust molecular markers and appropriate genotyping tools are required. Most current genotyping methods for MTBC are based on mobile or repetitive DNA elements. Because these elements are prone to convergent evolution, the corresponding genotyping techniques are suboptimal for phylogenetic studies and strain classification. By contrast, single nucleotide polymorphisms (SNP) are ideal markers for classifying MTBC into phylogenetic lineages, as they exhibit very low degrees of homoplasy. In this study, we developed two complementary SNP-based genotyping methods to classify strains into the six main human-associated lineages of MTBC, the 'Beijing' sublineage, and the clade comprising Mycobacterium bovis and Mycobacterium caprae. Phylogenetically informative SNPs were obtained from 22 MTBC whole-genome sequences. The first assay, referred to as MOL-PCR, is a ligation-dependent PCR with signal detection by fluorescent microspheres and a Luminex flow cytometer, which simultaneously interrogates eight SNPs. The second assay is based on six individual TaqMan real-time PCR assays for singleplex SNP-typing. We compared MOL-PCR and TaqMan results in two panels of clinical MTBC isolates. Both methods agreed fully when assigning 36 well-characterized strains into the main phylogenetic lineages. The sensitivity in allele-calling was 98.6% and 98.8% for MOL-PCR and TaqMan, respectively. Typing of an additional panel of 78 unknown clinical isolates revealed 99.2% and 100% sensitivity in allele-calling, respectively, and 100% agreement in lineage assignment between both methods. While MOL-PCR and TaqMan are both highly sensitive and specific, MOL-PCR is ideal for classification of isolates with no previous information, whereas TaqMan is faster for confirmation. Furthermore, both methods are rapid, flexible and comparably inexpensive

Multiple Locus Variable number of tandem repeat Analysis : a molecular genotyping tool for Paenibacillus larvae

American Foulbrood, caused by Paenibacillus larvae, is the most severe bacterial disease of honey bees (Apis mellifera). To perform genotyping of P.larvae in an epidemiological context, there is a need of a fast and cheap method with a high resolution. Here, we propose Multiple Locus Variable number of tandem repeat Analysis (MLVA). MLVA has been used for typing a collection of 209 P.larvae strains from which 23 different MLVA types could be identified. Moreover, the developed methodology not only permits the identification of the four Enterobacterial Repetitive Intergenic Consensus (ERIC) genotypes, but allows also a discriminatory subdivision of the most dominant ERIC type I and ERIC type II genotypes. A biogeographical study has been conducted showing a significant correlation between MLVA genotype and the geographical region where it was isolated

Plasmodium falciparum parasite population structure and gene flow associated to anti-malarial drugs resistance in Cambodia

Background: Western Cambodia is recognized as the epicentre of emergence of Plasmodium falciparum multi-drug resistance. The emergence of artemisinin resistance has been observed in this area since 2008–2009 and molecular signatures associated to artemisinin resistance have been characterized in k13 gene. At present, one of the major threats faced, is the possible spread of Asian artemisinin resistant parasites over the world threatening millions of people and jeopardizing malaria elimination programme efforts. To anticipate the diffusion of artemisinin resistance, the identification of the P. falciparum population structure and the gene flow among the parasite population in Cambodia are essential. Methods: To this end, a mid-throughput PCR-LDR-FMA approach based on LUMINEX technology was developed to screen for genetic barcode in 533 blood samples collected in 2010–2011 from 16 health centres in malaria endemics areas in Cambodia. Results: Based on successful typing of 282 samples, subpopulations were characterized along the borders of the country. Each 11-loci barcode provides evidence supporting allele distribution gradient related to subpopulations and gene flow. The 11-loci barcode successfully identifies recently emerging parasite subpopulations in western Cambodia that are associated with the C580Y dominant allele for artemisinin resistance in k13 gene. A subpopulation was identified in northern Cambodia that was associated to artemisinin (R539T resistant allele of k13 gene) and mefloquine resistance. Conclusions: The gene flow between these subpopulations might have driven the spread of artemisinin resistance over Cambodia

Dynamic telomerase gene suppression via network effects of GSK3 inhibition

<b>Background</b>: Telomerase controls telomere homeostasis and cell immortality and is a promising anti-cancer target, but few small molecule telomerase inhibitors have been developed. Reactivated transcription of the catalytic subunit hTERT in cancer cells controls telomerase expression. Better understanding of upstream pathways is critical for effective anti-telomerase therapeutics and may reveal new targets to inhibit hTERT expression. <b>Methodology/Principal Findings</b>: In a focused promoter screen, several GSK3 inhibitors suppressed hTERT reporter activity. GSK3 inhibition using 6-bromoindirubin-3′-oxime suppressed hTERT expression, telomerase activity and telomere length in several cancer cell lines and growth and hTERT expression in ovarian cancer xenografts. Microarray analysis, network modelling and oligonucleotide binding assays suggested that multiple transcription factors were affected. Extensive remodelling involving Sp1, STAT3, c-Myc, NFκB, and p53 occurred at the endogenous hTERT promoter. RNAi screening of the hTERT promoter revealed multiple kinase genes which affect the hTERT promoter, potentially acting through these factors. Prolonged inhibitor treatments caused dynamic expression both of hTERT and of c-Jun, p53, STAT3, AR and c-Myc. <b>Conclusions/Significance</b>: Our results indicate that GSK3 activates hTERT expression in cancer cells and contributes to telomere length homeostasis. GSK3 inhibition is a clinical strategy for several chronic diseases. These results imply that it may also be useful in cancer therapy. However, the complex network effects we show here have implications for either setting