Pyrosequencing: history, biochemistry and future

Molecular tools for deciphering the microbial community structure and diversity in rumen ecosystem

Rumen microbial community comprising of bacteria, archaea, fungi, and protozoa is characterized not only by the high population density but also by the remarkable diversity and the most complex microecological interactions existing in the biological world. This unprecedented biodiversity is quite far from full elucidation as only about 15-20 % of the rumen microbes are identified and characterized till date using conventional culturing and microscopy. However, the last two decades have witnessed a paradigm shift from cumbersome and time-consuming classical methods to nucleic acid-based molecular approaches for deciphering the rumen microbial community. These techniques are rapid, reproducible and allow both the qualitative and quantitative assessment of microbial diversity. This review describes the different molecular methods and their applications in elucidating the rumen microbial community.

Experimental analysis of sources of error in evolutionary studies based on Roche/454 pyrosequencing of viral genomes

Factors affecting the reliability of Roche/454 pyrosequencing for analyzing sequence polymorphism in within-host viral populations were assessed by two experiments: 1) sequencing four clonal simian immunodeficiency virus (SIV) stocks and 2) sequencing mixtures in different proportions of two SIV strains with known fixed nucleotide differences. Observed nucleotide diversity and frequency of undetermined nucleotides were increased at sites in homopolymer runs of four or more identical nucleotides, particularly at AT sites. However, in the mixed-strain experiments, the effects on estimated nucleotide diversity of such errors were small in comparison to known strain differences. The results suggest that biologically meaningful variants present at a frequency of around 10% and possibly much lower are easily distinguished from artifacts of the sequencing process. Analysis of the clonal stocks revealed numerous rare variants that showed the signature of purifying selection and that elimination of variants at frequencies of less than 1% reduced estimates of nucleotide diversity by about an order of magnitude. Thus, using a 1% frequency cutoff for accepting a variant as real represents a conservative standard, which may be useful in studies that are focused on the discovery of specific mutations (such as those conferring immune escape or drug resistance). On the other hand, if the goal is to estimate nucleotide diversity, an optimal strategy might be to include all observed variants (even those at less than 1% frequency), while masking out homopolymer runs of four or more nucleotides.

A virtual pyrogram generator to resolve complex pyrosequencing results

We report a freely available software program, Pyromaker, which generates simulated traces for pyrosequencing results based on user inputs. Simulated pyrograms can aid in the analysis of complex pyrosequencing results in which various hypothesized mutations can be tested, and the resultant pyrograms can be matched with the actual pyrogram. We validated the software using the actual pyrograms for common KRAS gene mutations as well as several mutations in the BRAF, GNAS, and p53 genes. We demonstrate that all 18 possible single-base mutations within codons 12 and 13 of KRAS generate unique pyrosequencing traces and highlight the distinctions between them. We further show that all reported codon 12 and 13 complex mutations produce unique pyrograms. However, some complex mutations are indistinguishable from single-base mutations. For complicated pyrograms, Pyromaker was used in two modes, one in which hypothesis-based simulated pyrograms were pattern-matched with the actual pyrograms. In a second strategy with only the pyrogram, Pyromaker was used to identify the underlying mutation by iteratively reconstructing the mutant pyrogram. Either strategy was able to successfully identify the complex mutations, which were confirmed by cloning and sequencing. Using two examples of KRAS codon 12 mutations (specifically GGT→TTT, G12F and GGT→GAG, G12E), we report which combinations of five approaches permit unambiguous mutation identification. The most efficient approach was found to be pyrosequencing with Pyromaker.

Molecular markers for detection and differentiation of Plasmodium falciparum and Plasmodium vivax in human blood samples by pyrosequencing

PCR amplification coupled with pyrosequencing was used to measure molecular markers that could be used to detect and differentiate Plasmodium falciparum and Plasmodium vivax in human blood samples. The detection rates were in agreement with the results of Giemsa-stained film microscopy, which is the current gold standard for detection. This method provides an exciting alternative for malaria diagnosis.

Characterization of the gastrointestinal microbiota in health and inflammatory bowel disease

The enteric bacterial flora play a key role in maintaining health. Inflammatory bowel disease is associated with quantitative and qualitative alterations in the microbiota. Early characterization of the microbiota involved culture-dependent techniques. The advent of metagenomic techniques, however, allows for structural and functional characterization using culture-independent methods. Changes in diversity, together with quantitative alterations in specific bacterial species, have been identified. The functional significance of these changes, and their pathogenic role, remain to be elucidated.

Detection of pathogens in water: from phylochips to qPCR to pyrosequencing

Waterborne pathogens pose a significant threat to human health and a proper assessment of microbial water quality is important for decision making regarding water infrastructure and treatment investments and eventually to provide early warning of disease, particularly given increasing global disasters associated with severe public health risks. Microbial water quality monitoring has undergone tremendous transition in recent years, with novel molecular tools beginning to offer rapid, high-throughput, sensitive and specific detection of a wide spectrum of microbial pathogens that challenge traditional culture-based techniques. High-density microarrays, quantitative real-time PCR (qPCR) and pyrosequencing which are considered to be breakthrough technologies borne out of the ‘molecular revolution’ are at present emerging rapidly as tools of pathogen detection and discovery. Future challenges lie in integrating these molecular tools with concentration techniques and bioinformatics platforms for unbiased guide of pathogen surveillance in water and developing standardized protocols.

GenoType MTBDRsl for molecular detection of second-line-drug and ethambutol resistance in Mycobacterium tuberculosis strains and clinical samples

The purpose of this study was to evaluate the GenoType MTBDRsl assay (Hain Lifescience GmbH, Nehren, Germany) for its ability to detect resistance to fluoroquinolones (FLQ), injectable second-line antibiotics [kanamycin (KM) and capreomycin (CM)], and ethambutol (EMB) in Mycobacterium tuberculosis clinical strains and directly in clinical samples. A total of 34 clinical strains were characterized with the Bactec 460 TB system. Fifty-four clinical samples from 16 patients (5 were smear negative and 49 were smear positive) were also tested directly. The corresponding isolates of the clinical specimens were also analyzed with the Bactec 460TB. When there was a discrepancy between assays, pyrosequencing was performed. The overall rates of concordance of the MTBDRsl and the Bactec 460TB for the detection of FLQ, KM/CM, and EMB susceptibility in clinical strains were 72.4% (21/29), 88.8% (24/27), and 67.6% (23/34), whereas for clinical samples, rates were 86.5% (45/52), 92.3% (48/52), and 56% (28/50), respectively. In conclusion, the GenoType MTBDRsl assay may be a useful tool for making early decisions regarding KM/CM susceptibility and to a lesser extent regarding FLQ and EMB susceptibility. The test is able to detect mutations in both clinical strains and samples with a short turnaround time. However, for correct management of patients with extensively drug-resistant tuberculosis, results must be confirmed by a phenotypical method.

Identification and subtyping of Francisella by pyrosequencing and signature matching of 16S rDNA fragments

AIMS

To analyse the V1 region of the 16S rDNA gene by a universal pyrosequencing protocol to identify and subtype Francisella in 31 strains from a repository collection and 96 patient isolates.

METHODS AND RESULTS

Pyrosequencing was used to determine the nucleotide sequence of PCR amplification products of the variable region (V1) of the 16S rDNA from 31 repository strains and 96 isolates from Swedish patients with ulceroglandular tularaemia. Pyrosequencing resulted in a 37 nucleotide sequence, specific for Francisella sp., for all repository strains and patient samples analysed. In addition, the isolates could be divided into two groups based on the analysis of a single nucleotide polymorphism in the sequence: one group included Francisella tularensis ssp. tularensis, ssp. holarctica and ssp. mediasiatica, whereas the other group included Francisella tularensis ssp. novicida and other species of Francisella. The analysis of samples taken from patients suffering from ulceroglandular tularaemia revealed that all isolates belonged to the first group comprising subspecies of F. tularensis virulent for humans.

CONCLUSIONS

The pyrosequencing analysis of the 16S rDNA V1 is a useful molecular tool for the rapid identification of suspected isolates of Francisella sp. in clinical or environmental samples.

SIGNIFICANCE AND IMPACT OF THE STUDY

Virulent F. tularensis ssp. causing ulceroglandular tularaemia, or those with a potential to be used in a bioterrorism event, could rapidly be discriminated from subspecies less virulent for humans.

The molecular epidemiology of parasite infections: tools and applications

Molecular epidemiology, broadly defined, is the application of molecular genetic techniques to the dynamics of disease in a population. In this review, we briefly describe molecular and analytical tools available for molecular epidemiological studies and then provide an overview of how they can be applied to better understand parasitic disease. A range of new molecular tools have been developed in recent years, allowing for the direct examination of parasites from clinical or environmental samples, and providing access to relatively cheap, rapid, high throughput molecular assays. At the same time, new analytical approaches, in particular those derived from coalescent theory, have been developed to provide more robust estimates of evolutionary processes and demographic parameters from multilocus, genotypic data. To date, the primary application of molecular epidemiology has been to provide specific and sensitive identification of parasites and to resolve taxonomic issues, particularly at the species level and below. Population genetic studies have also been used to determine the extent of genetic diversity among populations of parasites and the degree to which this diversity is associated with different host cycles or epidemiologically important phenotypes. Many of these studies have also shed new light on transmission cycles of parasites, particularly the extent to which zoonotic transmission occurs, and on the prevalence and importance of mixed infections with different parasite species or intraspecific variants (polyparasitism). A major challenge, and one which is now being addressed by an increasing number of studies, is to find and utilize genetic markers for complex traits of epidemiological significance, such as drug resistance, zoonotic potential and virulence.

Genomics in pediatric endocrinology-genetic disorders and new techniques

In the last few years, there have been remarkable advances in the development of new and more sophisticated genetic techniques. These have allowed a better understanding of the molecular mechanisms of genetically determined pediatric endocrine disorders and are paving the way for a radical change in diagnosis and treatment. This article introduces some of these concepts and some of the genetic techniques being used.

Whole cancer genome sequencing by next-generation methods

Traditional approaches to sequence analysis are widely used to guide therapy for patients with lung and colorectal cancer and for patients with melanoma, sarcomas (eg, gastrointestinal stromal tumor), and subtypes of leukemia and lymphoma. The next-generation sequencing (NGS) approach holds a number of potential advantages over traditional methods, including the ability to fully sequence large numbers of genes (hundreds to thousands) in a single test and simultaneously detect deletions, insertions, copy number alterations, translocations, and exome-wide base substitutions (including known "hot-spot mutations") in all known cancer-related genes. Adoption of clinical NGS testing will place significant demands on laboratory infrastructure and will require extensive computational expertise and a deep knowledge of cancer medicine and biology to generate truly useful "clinically actionable" reports. It is anticipated that continuing advances in NGS technology will lower the overall cost, speed the turnaround time, increase the breadth of genome sequencing, detect epigenetic markers and other important genomic parameters, and become applicable to smaller and smaller specimens, including circulating tumor cells and circulating free DNA in plasma.

Comprehensive next-generation cancer genome sequencing in the era of targeted therapy and personalized oncology

DNA sequence analysis has become a significant laboratory test in oncology, permitting treatment to become increasingly personalized for both solid tumors and hematologic malignancies. Traditional approaches to sequence analysis, including Sanger sequencing, pyrosequencing and allele-specific PCR, are now widely used to guide therapy for patients diagnosed with lung and colorectal cancer as well as for melanoma, sarcomas (e.g., gastrointestinal stromal tumors) and subtypes of leukemia and lymphoma. Traditional sequence analysis has been limited in bandwidth and throughput and as a result, has been focused exclusively on testing the most common aberrations in key genes or fully sequencing single genes. The massively parallel or next-generation sequencing (NGS) approach to DNA analysis holds a number of potential advantages over the traditional methods, including the ability to fully sequence large numbers of genes (hundreds to thousands) in a single test. Furthermore, NGS can simultaneously detect deletions, insertions, copy number alterations, translocations and exome-wide base substitutions (including known hot-spot mutations) in all known cancer-related genes. However, significant challenges, particularly with respect to demands on expertise and infrastructure, will have to be overcome to translate NGS to the bedside of the cancer patient. Extensive computational expertise is required to bring NGS into clinical context, and a deep knowledge of cancer medicine and cancer biology will be required to generate truly useful, so-called 'clinically actionable' reports for clinicians. While NGS is on the cusp of being launched as a clinical test, it may be expected that the near future will continue to bring major advances in the technology that will lower the overall cost, speed up the turnaround time, increase the breadth of genome sequencing, and detect epigenetic markers and other important genomic parameters, while becoming applicable to smaller and smaller specimens, including circulating tumor cells and circulating free DNA in plasma.

The role of molecular genetic analysis within the diagnostic haemato-oncology laboratory

The identification of the molecular genetic basis to many haematological malignancies along with the increased use of molecularly targeted therapy has heralded an increasing role for molecular genetic-based techniques. Demonstration of acquired changes such as the JAK2 V617F mutation within myeloproliferative neoplasms has quickly moved from a research setting to the diagnostic laboratory. Disease-specific genetic markers, such as the BCR-ABL1 fusion gene in chronic myeloid leukaemia, enable sensitive molecular genetic methods to be applied for the detection and quantification of low-level residual disease, allowing early identification of relapse. Consequently, molecular genetics now plays a crucial role in diagnosis, the identification of prognostic markers and monitoring of haematological malignancies. The development of high-throughput whole-genome approaches offers the potential to rapidly screen newly diagnosed patients for all disease-associated molecular genetic changes.

Ontogeny-driven rDNA rearrangement, methylation, and transcription, and paternal influence

Gene rearrangement occurs during development in some cell types and this genome dynamics is modulated by intrinsic and extrinsic factors, including growth stimulants and nutrients. This raises a possibility that such structural change in the genome and its subsequent epigenetic modifications may also take place during mammalian ontogeny, a process undergoing finely orchestrated cell division and differentiation. We tested this hypothesis by comparing single nucleotide polymorphism-defined haplotype frequencies and DNA methylation of the rDNA multicopy gene between two mouse ontogenic stages and among three adult tissues of individual mice. Possible influences to the genetic and epigenetic dynamics by paternal exposures were also examined for Cr(III) and acid saline extrinsic factors. Variables derived from litters, individuals, and duplicate assays in large mouse populations were examined using linear mixed-effects model. We report here that active rDNA rearrangement, represented by changes of haplotype frequencies, arises during ontogenic progression from day 8 embryos to 6-week adult mice as well as in different tissue lineages and is modifiable by paternal exposures. The rDNA methylation levels were also altered in concordance with this ontogenic progression and were associated with rDNA haplotypes. Sperm showed highest level of methylation, followed by lungs and livers, and preferentially selected haplotypes that are positively associated with methylation. Livers, maintaining lower levels of rDNA methylation compared with lungs, expressed more rRNA transcript. In vitro transcription demonstrated haplotype-dependent rRNA expression. Thus, the genome is also dynamic during mammalian ontogeny and its rearrangement may trigger epigenetic changes and subsequent transcriptional controls, that are further influenced by paternal exposures.

The effect of TP53 codon 72 and RNASEL codon 462 polymorphisms on the development of cervical cancer in Argentine women

Epidemiological evidence suggests that genetic factors, such as variants in cancer suppressor genes, may play an important role in the etiology of cervical carcinoma. TP53 is an outstanding cell cycle regulator, mutated in most human cancers, and RNASEL is thought to be involved in antiviral and apoptotic responses. To determine whether TP53 Arg72Pro and RNASEL Arg462Gln polymorphisms are associated with susceptibility to cervical cancer, a case-control study of 98 cancer patients and 123 healthy controls was conducted. Cervical samples were genotyped for both polymorphisms by pyrosequencing technology. The association between cervical cancer risk and the studied SNPs was evaluated by logistic regression, and potential gene-gene interactions were studied by Multifactor Dimensionality Reduction analysis. In the single-locus analysis, only the heterozygous TP53 Arg72Pro genotype was significantly associated with the risk of developing a cervical carcinoma, while the RNASEL polymorphism showed no association after age adjustment. In addition, the combination of both polymorphisms gives near-null information gain. Consequently, the effect provided by each single nucleotide polymorphism individually is considered higher than the effect resulting from the interaction between these two genes in cervical cancer risk. These results suggest that a heterozygous TP53 Arg72Pro genotype may contribute to cervical cancer susceptibility.

Detecting mutations in PfCRT and PfMDR1 genes among Plasmodium falciparum isolates from Saudi Arabia by pyrosequencing

The emergence of chloroquine resistance in Plasmodium falciparum is a significant public health problem where malaria is endemic. We aimed to evaluate the efficacy of pyrosequencing to assess chloroquine resistance among P. falciparum isolates from the southwestern region of Saudi Arabia by analyzing the K76T and N86Y mutations in the P. falciparum chloroquine resistance transporter (PfCRT) and P. falciparum multidrug resistance 1 (PfMDR1) genes, respectively. Blood samples (n = 121) from microscopically positive P. falciparum cases were collected. DNA was extracted, and fragments from each of the genes were amplified by PCR using new sets of primers. The amplicons were sequenced using a pyrosequencer. All of the 121 samples were amplified for assessment of the PfCRT K76T and PfMDR1 N86Y mutations. All of the samples amplified for the PfCRT 76T mutation harbored the ACA codon (121/121; 100%), indicating the presence of the 76T mutation. For the PfMDR1 N86Y mutation, 72/121 samples (59.5%) had the sequence AAT at that position, indicating the presence of the wild-type allele (86N). However, 49/121 samples (40.5%) had a TAT codon, indicating the mutant allele (Y) at position 86. This study shows that pyrosequencing could be useful as a high throughput, rapid, and sensitive assay for the detection of specific single nucleotide polymorphisms in drug-resistant P. falciparum strains. This will help health authorities in malaria-endemic regions to adopt new malaria control strategies that will be applicable for diagnostic and drug resistance assays for malaria and other life-threatening pathogens that are endemic in their respective countries.

Analysing DNA methylation using bisulphite pyrosequencing

Bisulphite pyrosequencing is a quantitative methodology for the investigation of DNA methylation of sequences up to 100-bp in length. Biotin-labelled, single-stranded PCR products generated from bisulphite-treated DNA are used as a template with an internal primer to perform the pyrosequencing reaction. Nucleotides are added in a predetermined order in each pyrosequencing cycle and the amount of incorporated nucleotide results in a proportional emission of light. DNA methylation ratios are calculated from the levels of light emitted from each nucleotide incorporated at individual CpG positions in a strand-dependent manner. The methylation detection limit at individual CpG sites is approximately 5% and the results are displayed as an average methylation level for each CpG position assayed across all amplification products generated during a PCR reaction. As a consequence, bisulphite pyrosequencing allows the identification of heterogeneous DNA methylation patterns but does not provide information at a single allele resolution. This methodology is suited to analyse short DNA sequences such as those typically extracted from formalin-fixed paraffin-embedded specimens. Nevertheless, longer PCR products can be sequenced by serial bisulphite pyrosequencing, which utilises tandem assays along the amplicon. The general information provided is applicable for all formats of current pyrosequencing instruments, however, a specific protocol for the PyroMark Q24 instrument is provided.

Picoliter DNA sequencing chemistry on an electrowetting-based digital microfluidic platform

The results of investigations into performing DNA sequencing chemistry on a picoliter-scale electrowetting digital microfluidic platform are reported. Pyrosequencing utilizes pyrophosphate produced during nucleotide base addition to initiate a process ending with detection through a chemiluminescence reaction using firefly luciferase. The intensity of light produced during the reaction can be quantified to determine the number of bases added to the DNA strand. The logic-based control and discrete fluid droplets of a digital microfluidic device lend themselves well to the pyrosequencing process. Bead-bound DNA is magnetically held in a single location, and wash or reagent droplets added or split from it to circumvent product dilution. Here we discuss the dispensing, control, and magnetic manipulation of the paramagnetic beads used to hold target DNA. We also demonstrate and characterize the picoliter-scale reaction of luciferase with adenosine triphosphate to represent the detection steps of pyrosequencing and all necessary alterations for working on this scale.

Common charge-shift mutation Glu65Lys in K+ channel β₁-Subunit KCNMB1: pleiotropic consequences for glomerular filtration rate and progressive renal disease

BACKGROUND

Glomerular filtration rate (GFR) is a heritable trait, and hyperfiltration (GFR increment in remnant nephrons) may accelerate renal functional decline in chronic kidney disease (CKD). Mesangial and vascular smooth myocytes control GFR by contraction, dependent on voltage-gated Ca(2+) influx, which is controlled by the regulatory β₁-subunit (KCNMB1) of large-conductance heteromeric K+ ('BK') channels. KCNMB1 gain-of-function variant Glu65Lys results in generalized vasorelaxation and thus protection against systemic hypertension. Here we asked whether the Glu65Lys variant influences GFR, in the basal state or during progressive renal decline.

METHODS

We explored Glu65Lys effects on GFR in three populations spanning two ethnicities and two diseases (hypertension and nephrosclerosis). GFR was either estimated (eGFR from serum creatinine) or directly measured (iothalamate clearance).

RESULTS

The 65Lys variant was relatively common, occurring on ∼5-10% of chromosomes in different biogeographic ancestry groups, and 65Lys carriers exhibited higher eGFR in two primary care populations: extreme BP values in Kaiser clinics (p = 0.029, accounting for ∼0.2% of trait variance), or treated hypertensives in VA clinics (p = 0.017, accounting for ∼0.9% of trait variance). In blacks with progressive renal disease (NIDDK AASK), 65Lys carriers displayed a steeper slope in GFR chronic decline (p = 0.030, accounting for ∼0.4% of trait variance), and Glu65Lys genotype also predicted time of onset of renal failure (log rank p = 0.019).

CONCLUSIONS

Common KCNMB1 gain-of-function variant Glu65Lys influences GFR, and 65Lys carriers exhibit not only elevated baseline GFR, but also more rapid GFR decline (and consequent development of renal failure) in CKD. The results suggest that profiling patients at Glu65Lys can assist in gauging renal prognosis as well as selection of rational therapy in hypertension with progressive renal disease.

Validation of a urine test and characterization of the putative genetic mutation for hyperuricosuria in Bulldogs and Black Russian Terriers

OBJECTIVE

To determine whether hyperuricosuria was a predisposing factor for urate urolithiasis in Bulldogs and Black Russian Terriers (BRTs) and to estimate the allele frequency of the Cys181Phe genetic mutation in urate transporter SLC2A9 in these breeds.

ANIMALS

192 Bulldogs, 101 BRTs, 10 Dalmatians, and 9 dogs of other breeds.

PROCEDURES

Uric acid (UA) and creatinine (Cr) concentrations were quantified in urine samples collected from all dogs via midstream catch during natural voiding. Buccal swab or blood samples were also obtained, and DNA was extracted and used to genotype SLC2A9 sequence variants by use of pyrosequencing assays. A urine test for hyperuricosuria was validated in adult dogs by comparing urinary UA:Cr ratios between known hyperuricosuric and nonhyperuricosuric dogs.

RESULTS

Significantly higher UA:Cr ratios were found in some Bulldogs and BRTs, compared with ratios in other dogs from these breeds. These dogs were also homozygous for the SLC2A9 Cys181Phe mutation. The allele frequency of the Cys181Phe mutation was 0.16 in Bulldogs and 0.51 in BRTs. On the basis of these allele frequencies, 3% of the Bulldog population and 27% of the BRT population were estimated to be hyperuricosuric.

CONCLUSIONS AND CLINICAL RELEVANCE

Results suggested the genetic mutation associated with hyperuricosuria, first identified in Dalmatians, also appears to cause hyperuricosuria in Bulldogs and BRTs, indicating that similar management strategies for urate urolithiasis can be used in these breeds. The allele frequency of the mutation was high in both breeds, and DNA testing can be used to select against the mutation.

Environmental proteomics: analysis of structure and function of microbial communities

Prokaryotic and eukaryotic microorganisms make a vital contribution to biogeochemical cycles by decomposing virtually all natural compounds and thereby exert a lasting effect on biosphere and climate. The rapidly growing number of metagenomic sequences together with revolutionary advances in bioinformatics and protein analyses have opened completely new horizons to investigate the molecular basis of such complex processes. Proteomics has contributed substantially to our understanding of individual organisms at the cellular level as it offers excellent possibilities to probe many protein functions and responses simultaneously. However, it has not yet been widely applied in microbial ecology, although most proteins have an intrinsic metabolic function which can be used to relate microbial activities to the identity of defined organisms in multispecies communities. Albeit still in its infancy, environmental proteomics enables simple protein cataloging, comparative and semi-quantitative proteomics, analyses of protein localization, discovery of post-translational modifications, and even determination of amino-acid sequences and genotypes by strain-resolved proteogenomics. This review traces the historical development of environmental proteomics and summarizes milestone publications in the field. In conclusion, we briefly discuss current limitations of microbial community proteomics but also the potential of emerging technologies to shape the future of metaproteome analyses.

Template tailoring: Accurate determination of heterozygous alleles using peptide nucleic acid and dideoxyNTP

Measurement of the length of DNA fragments plays a pivotal role in genetic mapping, disease diagnostics, human identification and forensic applications. PCR followed by electrophoresis is used for DNA length measurement of STRs, a process that requires labeled primers and allelic ladders as standards to avoid machine error. Sequencing-based approaches can be used for STR analysis to eliminate the requirement of labeled primers and allelic ladder. However, the limiting factor with this approach is unsynchronized polymerization in heterozygous sample analysis, in which alleles with different lengths can lead to imbalanced heterozygote peak height ratios. We have developed a rapid DNA length measurement method using peptide nucleic acid and dideoxy dNTPs to "tailor" DNA templates for accurate sequencing to overcome this hurdle. We also devised an accelerated "dyad" pyrosequencing strategy, such that the combined approach can be used as a faster, more accurate alternative to de novo sequencing. Dyad sequencing interrogates two bases at a time by allowing the polymerase to incorporate two nucleotides to DNA template, cutting the analysis time in half. In addition, for the first time, we show the effect of peptide nucleic acid as a blocking probe to stop polymerization, which is essential to analyze the heterozygous samples by sequencing. This approach provides a new platform for rapid and cost-effective DNA length measurement for STRs and resequencing of small DNA fragments.

Association between the catechol-O-methyltransferase Val158Met polymorphism and panic disorder: a replication

The association between the catechol-O-methyltransferase Val158Met polymorphism and panic disorder was studied in a Swedish sample of 211 patients and 452 controls. We found a significant excess of the Val allele in both male and female patients, the latter but not the former finding being in line with previous studies.

Pyrosequencing on nicked dsDNA generated by nicking endonucleases

Although the pyrosequencing method is simple and fast, the step of ssDNA preparation increases the cost, labor, and cross-contamination risk. In this paper, we proposed a method enabling pyrosequencing directly on dsDNA digested by nicking endonucleases (NEases). Recognition sequence of NEases was introduced using artificially mismatched bases in a PCR primer (in the case of genotyping) or a reverse-transcription primer (in the case of gene expression analysis). PCR products were treated to remove excess amounts of primers, nucleotides, and pyrophosphate (PPi) prior to sequencing. After the nicking reaction, pyrosequencing starts at the nicked 3' end, and extension reaction occurs when the added dNTP is complementary to the non-nicked strand. Although the activity of strand displacement by Klenow is limited, approximately 10 bases are accurately sequenced; this length is long enough for genotyping and SRPP-based differential gene expression analysis. It was observed that the signals of two allele-specific bases in a pyrogram from nicked dsDNA are highly quantitative, enabling quantitative determination of allele-specific templates; thus, Down's Syndrome diagnosis as well as differential gene expression analysis was successfully executed. The results indicate that pyrosequencing using nicked dsDNA as templates is a simple, inexpensive, and reliable way in either quantitative genotyping or gene expression analysis.

A label-free assay of exonuclease activity using a pyrosequencing technique

Enzymes with 3'-5' exonuclease activities are important in promoting the accuracy of DNA replication and DNA repair by proofreading. The alteration of the function of these enzymes by endogenous or exogenous effectors could, therefore, have a considerable impact on DNA replication and ultimately on genome integrity. We have developed a label-free high-throughput screening method for quantifying the effects of different reagents on exonuclease activity. The assay is based on a hairpin-forming biotinylated oligonucleotide substrate that contains one or more exonuclease-resistant phosphorothioate nucleotides. The activity and specificity of the selected 3'-5' exonuclease is determined indirectly using a sensitive pyrosequencing reaction after cleanup of the samples. In this pyrosequencing step, the amount of nucleotides filled into each position of the exonucleolytically degraded 3' end of the substrate can be recorded quantitatively and equals the amount of the nucleotides removed by the exonuclease. This system allows the estimation of both processivity and efficiency of the exonuclease activity. We have employed compounds reported in the literature to inhibit the exonuclease activities of either exonuclease III or the large fragment of polymerase I (Klenow fragment) to evaluate the assay.

Progression of chronic kidney disease: Adrenergic genetic influence on glomerular filtration rate decline in hypertensive nephrosclerosis

BACKGROUND

African-Americans are likely to develop hypertension and hypertensive nephrosclerosis. This grave prognosis, coupled with familial aggregation of end-stage renal disease (ESRD) in Blacks, prompts a search for genetic risk factors for ESRD. Recent evidence implicates a crucial role for the sympathetic nervous system in progressive renal disease.

METHODS

We used the African-American Study of Kidney Disease to probe whether beta2-adrenergic receptor (ADRB2) predicts glomerular filtration rate (GFR) decline rate. A total of 580 participants were included. Baseline GFR was 51.2 +/- 0.5 ml/min/1.73 m2. Subjects were randomized in a 2 x 3 block design: to intensively lowered (MAP < or = 92 mm Hg) versus 'usual' (MAP = 102-107 mm Hg) blood pressure goal groups, and also divided by three randomized antihypertensive drugs (ramipril, metoprolol, or amlodipine). We scored 4 SNPs at the ADRB2 locus.

RESULTS

Haplotypes at ADRB2 predicted chronic GFR decline rate, GFR declined more slowly in individuals with haplotype-1 (-804G-->173T-->16Gly-->27GIn), and faster in those who carried haplotype-3 (-804G-->173T-->16Arg-->27Gln). ADRB2 genotype interacted with antihypertensive drug class to influence GFR slope (p = 0.001-0.037). We extended our findings to an independent case/control sample of Black hypertensive ESRD, in which we found that variant Gly16Arg that tagged the GFR slope-determining ADRB2 haplotype also conferred risk for the ESRD trait in Blacks.

CONCLUSIONS

The GFR decline/progression rate in hypertensive renal disease is controlled in part by genetic variation within the adrenergic pathway.

A strategy for single nucleotide polymorphism analysis of chimerism for somatic cell therapy

BACKGROUND AIMS

Chimerism is an important outcome measure in hematopoietic cell transplantation as well as somatic cell therapy. Commonly used methods to estimate chimerism are restricted by either gender or inefficient sensitivity. In principle, real-time polymerase chain reaction (PCR)-based assays can be used to assess single nucleotide polymorphisms (SNP), which are a vast resource of molecular markers, and such assays demonstrate a substantially higher sensitivity (0.001%), but the specificity is unclear because of a low-level signal from mismatched sequences.

METHODS

In this study, we cloned 14 pairs of SNP selected from the SNP HapMap database and examined the specificity and sensitivity of their detection by real-time PCR using two primer/fluorescent probe pairs to allow genotyping of the two possible variant alleles. Clinical donor-recipient pairs from 18 families were used to explore the efficacy of using SNP assays to measure chimerism.

RESULTS

We found that the polymorphic nucleotide influences the ability to distinguish the signal generated by the target and mismatched sequences. Moreover, the specific fluorescent reporter probe can affect the difference in signal intensity between the target and mismatched sequences. Real-time PCR SNP assays can attain a sensitivity of 0.1-0.5% with 100% specificity. When comparing possible clinical donor-recipient pairs, we found an average 3.3 out of 14 SNP were informative.

CONCLUSIONS

By optimal selection of the polymorphic sequences and fluorescent reporter, the real-time PCR SNP assay is superior to the short-tandem repeat chimerism assay and broadly applicable. This strategy may be applied in future clinical trials of bone marrow cell therapy.

Detection of BRAF mutations in thyroid nodules by allele-specific PCR using a dual priming oligonucleotide system

With various methods, BRAF mutations have been detected in 73.4% to 87.1% of papillary thyroid carcinomas (PTCs) in Korea. We assessed the ability of allele-specific polymerase chain reaction using a dual priming oligonucleotide system, compared with direct sequencing and pyrosequencing, to detect BRAF mutations in fine-needle aspiration specimens from 85 patients undergoing thyroidectomy. Final pathologic diagnoses were 55 malignant lesions and 30 benign lesions. We detected BRAF mutations in 47 (90%) of 52 PTCs by at least 1 method. The sensitivity of the dual priming oligonucleotide system (88.5%) was slightly higher than that of direct sequencing (82.7%) and pyrosequencing (86.5%). The specificity and positive predictive value of all 3 methods were 100%. The dual priming oligonucleotide system is a simple, rapid, and reliable method for detecting BRAF mutations. This method may be a useful adjunct tool to improve the preoperative diagnostic accuracy of PTC in fine-needle aspiration biopsy specimens.

SNP-RFLPing 2: an updated and integrated PCR-RFLP tool for SNP genotyping

Background

PCR-restriction fragment length polymorphism (RFLP) assay is a cost-effective method for SNP genotyping and mutation detection, but the manual mining for restriction enzyme sites is challenging and cumbersome. Three years after we constructed SNP-RFLPing, a freely accessible database and analysis tool for restriction enzyme mining of SNPs, significant improvements over the 2006 version have been made and incorporated into the latest version, SNP-RFLPing 2.

Results

The primary aim of SNP-RFLPing 2 is to provide comprehensive PCR-RFLP information with multiple functionality about SNPs, such as SNP retrieval to multiple species, different polymorphism types (bi-allelic, tri-allelic, tetra-allelic or indels), gene-centric searching, HapMap tagSNPs, gene ontology-based searching, miRNAs, and SNP500Cancer. The RFLP restriction enzymes and the corresponding PCR primers for the natural and mutagenic types of each SNP are simultaneously analyzed. All the RFLP restriction enzyme prices are also provided to aid selection. Furthermore, the previously encountered updating problems for most SNP related databases are resolved by an on-line retrieval system.

Conclusions

The user interfaces for functional SNP analyses have been substantially improved and integrated. SNP-RFLPing 2 offers a new and user-friendly interface for RFLP genotyping that can be used in association studies and is freely available at http://bio.kuas.edu.tw/snp-rflping2.

Emerging molecular assays for detection and characterization of respiratory viruses

This article describes several emerging molecular assays that have potential applications in the diagnosis and monitoring of respiratory viral infections. These techniques include direct nucleic acid detection by quantum dots, loop-mediated isothermal amplification, multiplex ligation-dependent probe amplification, amplification using arbitrary primers, target-enriched multiplexing amplification, pyrosequencing, padlock probes, solid and suspension microarrays, and mass spectrometry. Several of these systems already are commercially available to provide multiplex amplification and high-throughput detection and identification of a panel of respiratory viral pathogens. Further validation and implementation of such emerging molecular assays in routine clinical virology services will enhance the rapid diagnosis of respiratory viral infections and improve patient care.

Influenza genome analysis using pyrosequencing method: current applications for a moving target

Pyrosequencing is a high-throughput non-gel-based DNA sequencing method that was introduced in the late 1990s. It employs a DNA sequencing-by-synthesis approach based on real-time measurement of pyrophosphate released from incorporation of dNTPs. A cascade of enzymatic reactions proportionally converts the pyrophosphate to a light signal recorded in a form of peaks, known as pyrograms. Routinely, a 45-60-nucleotide sequence is obtained per reaction. Recent improvements introduced in the assay chemistry have extended the read to approximately 100 nucleotides. Since its advent, pyrosequencing has been applied in the fields of microbiology, molecular biology and pharmacogenomics. The pyrosequencing approach was first applied to analysis of influenza genome in 2005, when it played a critical role in the timely detection of an unprecedented rise in resistance to the adamantane class of anti-influenza drugs. More recently, pyrosequencing was successfully applied for monitoring the emergence and spread of influenza A (H1N1) virus resistance to oseltamivir, a newer anti-influenza drug. The present report summarizes known applications of the pyrosequencing approach for influenza genome analysis with an emphasis on drug-resistance detection.

Application of pyrosequencing to the identification of sequence variations in the cystic fibrosis transmembrane conductance regulator gene

BACKGROUND

A high number of mutations associated with cystic fibrosis have been identified in the cystic fibrosis transmembrane conductance regulator (CFTR) gene, most of which are rare and, therefore, hamper extensive molecular diagnosis. In couples undergoing prenatal diagnosis where no mutation is found in one or both partners, additional analysis of intragenic polymorphisms may allow for the identification of fetal alleles associated with cystic fibrosis.

METHODS

We developed novel, rapid and accurate assays for CFTR genotype determination using pyrosequencing technology; a simple, automated and reliable technique with low cost.

RESULTS

Assays were optimized for the identification of the seven most frequent CFTR mutations (p.DeltaF508, p.N1303K, p.G542X, c.2183AA>G, c.1717-1G>A, p.W1282X, p.R1162X) in the Italian population and two common intragenic polymorphisms (rs213950 and rs1800136). Blind validation on 15 known control samples, typed for each sequence variation, allowed correct identification of all 135 genotypes.

CONCLUSIONS

We demonstrated that this procedure is highly specific for the identification of individual CFTR sequence variations associated with cystic fibrosis, allowing both population screening and prenatal diagnosis.

Genomics in pediatric endocrinology--genetic disorders and new techniques

In the last few years, there have been remarkable advances in the development of new and more sophisticated genetic techniques. These have allowed a better understanding of the molecular mechanisms of genetically determined pediatric endocrine disorders and are paving the way for a radical change in diagnosis and treatment. This article introduces some of these concepts and some of the genetic techniques being used.

Rapid high-throughput human leukocyte antigen typing by massively parallel pyrosequencing for high-resolution allele identification

Transplantation and, notably, hematopoietic stem cell transplantation require high-resolution human leukocyte antigen (HLA) typing and, because of the heterozygous genomic DNA samples, are dependent on clonal analytical methods. High-resolution HLA typing is a necessity for accomplishing the best possible histocompatibility match between donor and recipient, because mismatches strongly increase the risk of severe acute graft-versus-host disease. We describe the development and first application in a clinical setting of a novel, HLA sequence-based typing method by exploring the next-generation sequencing technology as provided by the Genome Sequencer FLX system (Roche/454 Life Sciences, Branford, CT). The developed system allows for ambiguity-free, high-throughput, high-resolution HLA-A and -B typing with the potential for automation. Primers and Genome Sequencer FLX specific adapters were lengthened with donor-identifying barcode sequences to identify each of eight Caucasian reference donors within one single multiplex sequencing run. Compared with normal SBT HLA typing, results indicate that every patient was identified correctly with an average of 1000 reads per amplicon. Furthermore, current investments for increased read lengths and fully automated molecular diagnostic software tools, using original GS-FLX data file formats, will enhance this novel HLA typing strategy in the near future.

Application of massive parallel sequencing to whole genome SNP discovery in the porcine genome

BACKGROUND

Although the Illumina 1 G Genome Analyzer generates billions of base pairs of sequence data, challenges arise in sequence selection due to the varying sequence quality. Therefore, in the framework of the International Porcine SNP Chip Consortium, this pilot study aimed to evaluate the impact of the quality level of the sequenced bases on mapping quality and identification of true SNPs on a large scale.

RESULTS

DNA pooled from five animals from a commercial boar line was digested with DraI; 150-250-bp fragments were isolated and end-sequenced using the Illumina 1 G Genome Analyzer, yielding 70,348,064 sequences 36-bp long. Rules were developed to select sequences, which were then aligned to unique positions in a reference genome. Sequences were selected based on quality, and three thresholds of sequence quality (SQ) were compared. The highest threshold of SQ allowed identification of a larger number of SNPs (17,489), distributed widely across the pig genome. In total, 3,142 SNPs were validated with a success rate of 96%. The correlation between estimated minor allele frequency (MAF) and genotyped MAF was moderate, and SNPs were highly polymorphic in other pig breeds. Lowering the SQ threshold and maintaining the same criteria for SNP identification resulted in the discovery of fewer SNPs (16,768), of which 259 were not identified using higher SQ levels. Validation of SNPs found exclusively in the lower SQ threshold had a success rate of 94% and a low correlation between estimated MAF and genotyped MAF. Base change analysis suggested that the rate of transitions in the pig genome is likely to be similar to that observed in humans. Chromosome X showed reduced nucleotide diversity relative to autosomes, as observed for other species.

CONCLUSION

Large numbers of SNPs can be identified reliably by creating strict rules for sequence selection, which simultaneously decreases sequence ambiguity. Selection of sequences using a higher SQ threshold leads to more reliable identification of SNPs. Lower SQ thresholds can be used to guarantee sufficient sequence coverage, resulting in high success rate but less reliable MAF estimation. Nucleotide diversity varies between porcine chromosomes, with the X chromosome showing less variation as observed in other species.

Analysis of splicing patterns by pyrosequencing

Several different mRNAs can be produced from a given pre-mRNA by regulated alternative splicing, or as the result of deregulations that may lead to pathological states. Analysing splicing patterns is therefore of importance to describe and understand developmental programs, cellular responses to internal or external cues, or human diseases. We describe here a method, Pyrosequencing Analysis of Splicing Patterns (PASP), that combines RT–PCR and pyrosequencing of PCR products. We demonstrated that: (i) Ratios of two pure RNAs mixed in various proportions were accurately measured by PASP; (ii) PASP can be adapted to virtually any splicing event, including mutually exclusive exons, complex patterns of exon skipping or inclusion, and alternative 3′ terminal exons; (iii) In extracts from different organs, the proportions of RNA isoforms measured by PASP reflected those measured by other methods. The PASP method is therefore reliable for analysing splicing patterns. All steps are done in 96-wells microplates, without gel electrophoresis, opening the way to high-throughput comparisons of RNA from several sources.

Prim-SNPing: a primer designer for cost-effective SNP genotyping

Many kinds of primer design (PD) software tools have been developed, but most of them lack a single nucleotide polymorphism (SNP) genotyping service. Here, we introduce the web-based freeware "Prim-SNPing," which, in addition to general PD, provides three kinds of primer design functions for cost-effective SNP genotyping: natural PD, mutagenic PD, and confronting two-pair primers (CTPP) PD. The natural PD and mutagenic PD provide primers and restriction enzyme mining for polymerase chain reaction-restriction fragment of length polymorphism (PCR-RFLP), while CTPP PD provides primers for restriction enzyme-free SNP genotyping. The PCR specificity and efficiency of the designed primers are improved by BLAST searching and evaluating secondary structure (such as GC clamps, dimers, and hairpins), respectively. The length pattern of PCR-RFLP using natural PD is user-adjustable, and the restriction sites of the RFLP enzymes provided by Prim-SNPing are confirmed to be absent within the generated PCR product. In CTPP PD, the need for a separate digestion step in RFLP is eliminated, thus making it faster and cheaper. The output of Prim-SNPing includes the primer list, melting temperature (Tm) value, GC percentage, and amplicon size with enzyme digestion information. The reference SNP (refSNP, or rs) clusters from the Single Nucleotide Polymorphism database (dbSNP) at the National Center for Biotechnology Information (NCBI), and multiple other formats of human, mouse, and rat SNP sequences are acceptable input. In summary, Prim-SNPing provides interactive, user-friendly and cost-effective primer design for SNP genotyping. It is freely available at http://bio.kuas.edu.tw/prim-snping.