Edesign: Primer and Enhanced Internal Probe Design Tool for Quantitative PCR Experiments and Genotyping Assays

Analytical PCR experiments preferably use internal probes for monitoring the amplification reaction and specific detection of the amplicon. Such internal probes have to be designed in close context with the amplification primers, and may require additional considerations for the detection of genetic variations. Here we describe Edesign, a new online and stand-alone tool for designing sets of PCR primers together with an internal probe for conducting quantitative real-time PCR (qPCR) and genotypic experiments. Edesign can be used for selecting standard DNA oligonucleotides like for instance TaqMan probes, but has been further extended with new functions and enhanced design features for Eprobes. Eprobes, with their single thiazole orange-labelled nucleotide, allow for highly sensitive genotypic assays because of their higher DNA binding affinity as compared to standard DNA oligonucleotides. Using new thermodynamic parameters, Edesign considers unique features of Eprobes during primer and probe design for establishing qPCR experiments and genotyping by melting curve analysis. Additional functions in Edesign allow probe design for effective discrimination between wild-type sequences and genetic variations either using standard DNA oligonucleotides or Eprobes. Edesign can be freely accessed online at http://www.dnaform.com/edesign2/, and the source code is available for download.     

Nucleotide sequences of the 5' termini of phi X174 mRNAs synthesized in vitro

When using X174 RFI DNA as a template, in vitro, E. coli RNA polymerase synthesizes four major purine triphosphate-containing 5′ end sequences. RNase A digests of α³²P labeled RNA were further digested with spleen exonuclease to remove the bulk of the oligonucleotides with 5′ hydroxyls and then chromatographed on DEAE cellulose to resolve the remaining 5′ terminal oligonucleotides. By application of standard separation and sequence techniques, the major 5′ end sequences were shown to be: pppApUp(Cp), pppApApApUp(Cp), pppApApApApUp(Cp), and pppGpApUp(Gp).     

Automated Analyses of Innate Olfactory Behaviors in Rodents

Olfaction based behavioral experiments are important for the investigation of sensory coding, perception, decision making and memory formation. The predominant experimental paradigms employ forced choice operant assays, which require associative learning and reinforced training. Animal performance in these assays not only reflects odor perception but also the confidence in decision making and memory. In this study, we describe a versatile and automated setup, “Poking-Registered Olfactory Behavior Evaluation System” (PROBES), which can be adapted to perform multiple olfactory assays. In addition to forced choice assays, we employ this system to examine animal’s innate ability for odor detection, discrimination and preference without elaborate training procedures. These assays provide quantitative measurements of odor discrimination and robust readouts of odor preference. Using PROBES, we find odor detection thresholds are at lower concentrations in naïve animals than those determined by forced choice assays. PROBES-based automated assays provide an efficient way of analyzing innate odor-triggered behaviors.     

Structure/Function Analysis of PARP-1 in Oxidative and Nitrosative Stress-Induced Monomeric ADPR Formation

Poly adenosine diphosphate-ribose polymerase-1 (PARP-1) is a multifunctional enzyme that is involved in two major cellular responses to oxidative and nitrosative (O/N) stress: detection and response to DNA damage via formation of protein-bound poly adenosine diphosphate-ribose (PAR), and formation of the soluble 2^nd messenger monomeric adenosine diphosphate-ribose (mADPR). Previous studies have delineated specific roles for several of PARP-1′s structural domains in the context of its involvement in a DNA damage response. However, little is known about the relationship between the mechanisms through which PARP-1 participates in DNA damage detection/response and those involved in the generation of monomeric ADPR. To better understand the relationship between these events, we undertook a structure/function analysis of PARP-1 via reconstitution of PARP-1 deficient DT40 cells with PARP-1 variants deficient in catalysis, DNA binding, auto-PARylation, and PARP-1′s BRCT protein interaction domain. Analysis of responses of the respective reconstituted cells to a model O/N stressor indicated that PARP-1 catalytic activity, DNA binding, and auto-PARylation are required for PARP-dependent mADPR formation, but that BRCT-mediated interactions are dispensable. As the BRCT domain is required for PARP-dependent recruitment of XRCC1 to sites of DNA damage, these results suggest that DNA repair and monomeric ADPR 2^nd messenger generation are parallel mechanisms through which PARP-1 modulates cellular responses to O/N stress.     

Thermal denaturation of B. subtilis DNA in H2O and D2O observed by electron microscopy

In the present work we studied the initial part of thermal denaturation curves of B. subtilis DNA in normal and heavy water, observing, by electron microscopy, the “opening” of DNA as a function of temperature. The results support the hypothesis that D₂O plays a stabilizing role on strands separation during thermal denaturation of DNA.     

HEK293T Cells Are Heterozygous for CCR5 Delta 32 Mutation

C-C chemokine receptor 5 (CCR5) is a receptor for chemokines and a co-receptor for HIV-1 entry into the target CD4+ cells. CCR5 delta 32 deletion is a loss-of-function mutation, resistant to HIV-1 infection. We tried to induce the CCR5 delta 32 mutation harnessing the genome editing technique, CRISPR-Cas9 (Clustered Regularly Interspaced Short Palindromic Repeats, CRISPR and CRISPR associated protein 9, Cas9) in the commonly used cell line human embryonic kidney HEK 293T cells. Surprisingly, we found that HEK293T cells are heterozygous for CCR5 delta 32 mutation, in contrast to the wild type CCR5 cells, human acute T cell leukemia cell line Jurkat and human breast adenocarcinoma cell line MDA-MB-231 cells. This finding indicates that at least one human cell line is heterozygous for the CCR5 delta 32 mutation. We also found that in PCR amplification, wild type CCR5 DNA and mutant delta 32 DNA can form mismatched heteroduplex and move slowly in gel electrophoresis.     

DISSECT Method Using PNA-LNA Clamp Improves Detection of EGFR T790m Mutation

Non-small cell lung cancer (NSCLC) patients treated with small molecule EGFR inhibitors, such as gefitinib, frequently develop drug resistance due to the presence of secondary mutations like the T790M mutation on EGFR exon 20. These mutations may originate from small subclonal populations in the primary tumor that become dominant later on during treatment. In order to detect these low-level DNA variations in the primary tumor or to monitor their progress in plasma, it is important to apply reliable and sensitive mutation detection methods. Here, we combine two recently developed methodologies, Differential Strand Separation at Critical Temperature (DISSECT), with peptide nucleic acid-locked nucleic acid (PNA-LNA) polymerase chain reaction (PCR) for the detection of T790M EGFR mutation. DISSECT pre-enriches low-abundance T790M EGFR mutations from target DNA prior to implementing PNA-LNA PCR, a method that can detect 1 mutant allele in a background of 100–1000 wild type alleles. The combination of DISSECT and PNA-LNA PCR enables the detection of 1 mutant allele in a background of 10,000 wild type alleles. The combined DISSECT-PNA-LNA PCR methodology is amenable to adaptation for the sensitive detection of additional emerging resistance mutations in cancer.     

Competitive action between Brassinosteroid and tracheary element differentiation inhibitory factor in controlling xylem cell differentiation

For permanent secondary growth in plants, cell proliferation and differentiation should be strictly controlled in the vascular meristem consisting of (pro)cambial cells. A peptide hormone tracheary element differentiation inhibitory factor (TDIF) functions to inhibit xylem differentiation, while a plant hormone brassinosteroid (BR) promotes xylem differentiation in (pro)cambial cells. However, it remains unclear how TDIF and BR cooperate to regulate xylem differentiation for the proper maintenance of the vascular meristem. In this study, I developed an easy evaluation method for xylem differentiation frequency in a vascular induction system Vascular cell Induction culture System Using Arabidopsis Leaves (VISUAL) by utilizing a xylem-specific luciferase reporter line. In this quantitative system, TDIF suppressed and BR promoted xylem differentiation in a dose-dependent manner, respectively. Moreover, simultaneous treatment of TDIF and BR with (pro)cambial cells revealed that they can cancel their each other’s effect on xylem differentiation, suggesting a competitive relationship between TDIF and BR. Thus, mutual inhibition of “ON” and “OFF” signal enables the fine-tuned regulation of xylem differentiation in the vascular meristem.     

A pre-ribosomal RNA interaction network involving snoRNAs and the Rok1 helicase

Ribosome biogenesis in yeast requires 75 small nucleolar RNAs (snoRNAs) and a myriad of cofactors for processing, modification, and folding of the ribosomal RNAs (rRNAs). For the 19 RNA helicases implicated in ribosome synthesis, their sites of action and molecular functions have largely remained unknown. Here, we have used UV cross-linking and analysis of cDNA (CRAC) to reveal the pre-rRNA binding sites of the RNA helicase Rok1, which is involved in early small subunit biogenesis. Several contact sites were identified in the 18S rRNA sequence, which interestingly all cluster in the “foot” region of the small ribosomal subunit. These include a major binding site in the eukaryotic expansion segment ES6, where Rok1 is required for release of the snR30 snoRNA. Rok1 directly contacts snR30 and other snoRNAs required for pre-rRNA processing. Using cross-linking, ligation and sequencing of hybrids (CLASH) we identified several novel pre-rRNA base-pairing sites for the snoRNAs snR30, snR10, U3, and U14, which cluster in the expansion segments of the 18S rRNA. Our data suggest that these snoRNAs bridge interactions between the expansion segments, thereby forming an extensive interaction network that likely promotes pre-rRNA maturation and folding in early pre-ribosomal complexes and establishes long-range rRNA interactions during ribosome synthesis.     

The molecular complexity of mu and pi symbiont DNA of Paramecium aurelia

The molecular size of mu and pi symbionts of Parameciumaurelia has been calculated from renaturation kinetic data. Observed values were 0.78 × 10⁹ daltons for mu particle DNA and 0.81 × 10⁹ daltons for pi particle DNA. Estimates of analytical complexity were 4.45 × 10⁹ and 5.05 × 10⁹ daltons respectively. Based on these data, mu and pi symbionts appear to possess multiple genomes and contain a minimum of 5 or 6 copies of each DNA sequence.     

Quantification of microtubule stutters: dynamic instability behaviors that are strongly associated with catastrophe

Microtubules (MTs) are cytoskeletal fibers that undergo dynamic instability (DI), a remarkable process involving phases of growth and shortening separated by stochastic transitions called catastrophe and rescue. Dissecting DI mechanism(s) requires first characterizing and quantifying these dynamics, a subjective process that often ignores complexity in MT behavior. We present a Statistical Tool for Automated Dynamic Instability Analysis (STADIA) that identifies and quantifies not only growth and shortening, but also a category of intermediate behaviors that we term “stutters.” During stutters, the rate of MT length change tends to be smaller in magnitude than during typical growth or shortening phases. Quantifying stutters and other behaviors with STADIA demonstrates that stutters precede most catastrophes in our in vitro experiments and dimer-scale MT simulations, suggesting that stutters are mechanistically involved in catastrophes. Related to this idea, we show that the anticatastrophe factor CLASP2γ works by promoting the return of stuttering MTs to growth. STADIA enables more comprehensive and data-driven analysis of MT dynamics compared with previous methods. The treatment of stutters as distinct and quantifiable DI behaviors provides new opportunities for analyzing mechanisms of MT dynamics and their regulation by binding proteins.     

Mutation Detection by Real-Time PCR: A Simple,Robust and Highly Selective Method

Background

Molecular tests for diagnosis of disease, particularly cancer, are gaining increased acceptance by physicians and their patients for disease prognosis and selection of treatment options. Gene expression profiles and genetic mutations are key parameters used for the molecular characterization of tumors. A variety of methods exist for mutation analysis but the development of assays with high selectivity tends to require a process of trial and error, and few are compatible with real-time PCR. We sought to develop a real-time PCR-based mutation assay methodology that successfully addresses these issues.

Methodology/Principal Findings

The method we describe is based on the widely used TaqMan® real-time PCR technology, and combines Allele-Specific PCR with a Blocking reagent (ASB-PCR) to suppress amplification of the wildype allele. ASB-PCR can be used for detection of germ line or somatic mutations in either DNA or RNA extracted from any type of tissue, including formalin-fixed paraffin-embedded tumor specimens. A set of reagent design rules was developed enabling sensitive and selective detection of single point substitutions, insertions, or deletions against a background of wild-type allele in thousand-fold or greater excess.

Conclusions/Significance

ASB-PCR is a simple and robust method for assaying single nucleotide mutations and polymorphisms within the widely used TaqMan® protocol for real time RT-PCR. The ASB-PCR design rules consistently produce highly selective mutation assays while obviating the need for redesign and optimization of the assay reagents. The method is compatible with formalin-fixed tissue and simultaneous analysis of gene expression by RT-PCR on the same plate. No proprietary reagents other than those for TaqMan chemistry are required, so the method can be performed in any research laboratory with real-time PCR capability.     

The Evolutionary Divergence of Shiga Toxin-Producing Escherichia coli Is Reflected in Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR) Spacer Composition

The Shiga toxin-producing Escherichia coli (STEC) strains, including those of O157:H7 and the “big six” serogroups (i.e., serogroups O26, O45, O103, O111, O121, and O145), are a group of pathogens designated food adulterants in the United States. The relatively conserved nature of clustered regularly interspaced short palindromic repeats (CRISPRs) in phylogenetically related E. coli strains makes them potential subtyping markers for STEC detection, and a quantitative PCR (qPCR)-based assay was previously developed for O26:H11, O45:H2, O103:H2, O111:H8, O121:H19, O145:H28, and O157:H7 isolates. To better evaluate the sensitivity and specificity of this qPCR method, the CRISPR loci of 252 O157 and big-six STEC isolates were sequenced and analyzed along with 563 CRISPR1 and 624 CRISPR2 sequences available in GenBank. General conservation of spacer content and order was observed within each O157 and big-six serogroup, validating the qPCR method. Meanwhile, it was found that spacer deletion, the presence of an insertion sequence, and distinct alleles within a serogroup are sources of false-negative reactions. Conservation of CRISPR arrays among isolates expressing the same flagellar antigen, specifically, H7, H2, and H11, suggested that these isolates share an ancestor and provided an explanation for the false positives previously observed in the qPCR results. An analysis of spacer distribution across E. coli strains provided limited evidence for temporal spacer acquisition. Conversely, comparison of CRISPR sequences between strains along the stepwise evolution of O157:H7 from its O55:H7 ancestor revealed that, over this ∼7,000-year span, spacer deletion was the primary force generating CRISPR diversity.     

Pgrmc1 (Progesterone Receptor Membrane Component 1) Associates with Epidermal Growth Factor Receptor and Regulates Erlotinib Sensitivity

Tumorigenesis requires the concerted action of multiple pathways, including pathways that stimulate proliferation and metabolism. Epidermal growth factor receptor (EGFR) is a transmembrane receptor-tyrosine kinase that is associated with cancer progression, and the EGFR inhibitors erlotinib/tarceva and tyrphostin/AG-1478 are potent anti-cancer therapeutics. Pgrmc1 (progesterone receptor membrane component 1) is a cytochrome b₅-related protein that is up-regulated in tumors and promotes cancer growth. Pgrmc1 and its homologues have been implicated in cell signaling, and we show here that Pgrmc1 increases susceptibility to AG-1478 and erlotinib, increases plasma membrane EGFR levels, and co-precipitates with EGFR. Pgrmc1 co-localizes with EGFR in cytoplasmic vesicles and co-fractionates with EGFR in high density microsomes. The findings have therapeutic potential because a Pgrmc1 small molecule ligand, which inhibits growth in a variety of cancer cell types, de-stabilized EGFR in multiple tumor cell lines. EGFR is one of the most potent receptor-tyrosine kinases driving tumorigenesis, and our data support a role for Pgrmc1 in promoting several cancer phenotypes at least in part by binding EGFR and stabilizing plasma membrane pools of the receptor.     

A Disintegrin and Metalloprotease 17 Dynamic Interaction Sequence,the Sweet Tooth for the Human Interleukin 6 Receptor

A disintegrin and metalloprotease 17 (ADAM17) is a major sheddase involved in the regulation of a wide range of biological processes. Key substrates of ADAM17 are the IL-6 receptor (IL-6R) and TNF-α. The extracellular region of ADAM17 consists of a prodomain, a catalytic domain, a disintegrin domain, and a membrane-proximal domain as well as a small stalk region. This study demonstrates that this juxtamembrane segment is highly conserved, α-helical, and involved in IL-6R binding. This process is regulated by the structure of the preceding membrane-proximal domain, which acts as molecular switch of ADAM17 activity operated by a protein-disulfide isomerase. Hence, we have termed the conserved stalk region “Conserved ADAM seventeen dynamic interaction sequence” (CANDIS). Finally, we identified the region in IL-6R that binds to CANDIS. In contrast to the type I transmembrane proteins, the IL-6R, and IL-1RII, CANDIS does not bind the type II transmembrane protein TNF-α, demonstrating fundamental differences in the respective shedding by ADAM17.     

Structural basis for the sequence-dependent effects of platinum–DNA adducts

The differences in efficacy and molecular mechanisms of platinum based anti-cancer drugs cisplatin (CP) and oxaliplatin (OX) have been hypothesized to be in part due to the differential binding affinity of cellular and damage recognition proteins to CP and OX adducts formed on adjacent guanines in genomic DNA. HMGB1a in particular exhibits higher binding affinity to CP-GG adducts, and the extent of discrimination between CP- and OX-GG adducts is dependent on the bases flanking the adducts. However, the structural basis for this differential binding is not known. Here, we show that the conformational dynamics of CP- and OX-GG adducts are distinct and depend on the sequence context of the adduct. Molecular dynamics simulations of the Pt-GG adducts in the TGGA sequence context revealed that even though the major conformations of CP- and OX-GG adducts were similar, the minor conformations were distinct. Using the pattern of hydrogen bond formation between the Pt–ammines and the adjacent DNA bases, we identified the major and minor conformations sampled by Pt–DNA. We found that the minor conformations sampled exclusively by the CP-GG adduct exhibit structural properties that favor binding by HMGB1a, which may explain its higher binding affinity to CP-GG adducts, while these conformations are not sampled by OX-GG adducts because of the constraints imposed by its cyclohexane ring, which may explain the negligible binding affinity of HMGB1a for OX-GG adducts in the TGGA sequence context. Based on these results, we postulate that the constraints imposed by the cyclohexane ring of OX affect the DNA conformations explored by OX-GG adduct compared to those of CP-GG adduct, which may influence the binding affinities of HMG-domain proteins for Pt-GG adducts, and that these conformations are further influenced by the DNA sequence context of the Pt-GG adduct.     

The Construction and Use of Log-Odds Substitution Scores for Multiple Sequence Alignment

Most pairwise and multiple sequence alignment programs seek alignments with optimal scores. Central to defining such scores is selecting a set of substitution scores for aligned amino acids or nucleotides. For local pairwise alignment, substitution scores are implicitly of log-odds form. We now extend the log-odds formalism to multiple alignments, using Bayesian methods to construct “BILD” (“Bayesian Integral Log-odds”) substitution scores from prior distributions describing columns of related letters. This approach has been used previously only to define scores for aligning individual sequences to sequence profiles, but it has much broader applicability. We describe how to calculate BILD scores efficiently, and illustrate their uses in Gibbs sampling optimization procedures, gapped alignment, and the construction of hidden Markov model profiles. BILD scores enable automated selection of optimal motif and domain model widths, and can inform the decision of whether to include a sequence in a multiple alignment, and the selection of insertion and deletion locations. Other applications include the classification of related sequences into subfamilies, and the definition of profile-profile alignment scores. Although a fully realized multiple alignment program must rely upon more than substitution scores, many existing multiple alignment programs can be modified to employ BILD scores. We illustrate how simple BILD score based strategies can enhance the recognition of DNA binding domains, including the Api-AP2 domain in Toxoplasma gondii and Plasmodium falciparum.     

Distinct Patterns of DNA Damage Response and Apoptosis Correlate with Jak/Stat and PI3Kinase Response Profiles in Human Acute Myelogenous Leukemia

Background

Single cell network profiling (SCNP) utilizing flow cytometry measures alterations in intracellular signaling responses. Here SCNP was used to characterize Acute Myeloid Leukemia (AML) disease subtypes based on survival, DNA damage response and apoptosis pathways.

Methodology and Principal Findings

Thirty four diagnostic non-M3 AML samples from patients with known clinical outcome were treated with a panel of myeloid growth factors and cytokines, as well as with apoptosis-inducing agents. Analysis of induced Jak/Stat and PI3K pathway responses in blasts from individual patient samples identified subgroups with distinct signaling profiles that were not seen in the absence of a modulator. In vitro exposure of patient samples to etoposide, a DNA damaging agent, revealed three distinct “DNA damage response (DDR)/apoptosis” profiles: 1) AML blasts with a defective DDR and failure to undergo apoptosis; 2) AML blasts with proficient DDR and failure to undergo apoptosis; 3) AML blasts with proficiency in both DDR and apoptosis pathways. Notably, AML samples from clinical responders fell within the “DDR/apoptosis” proficient profile and, as well, had low PI3K and Jak/Stat signaling responses. In contrast, samples from clinical non responders had variable signaling profiles often with in vitro apoptotic failure and elevated PI3K pathway activity. Individual patient samples often harbored multiple, distinct, leukemia-associated cell populations identifiable by their surface marker expression, functional performance of signaling pathway in the face of cytokine or growth factor stimulation, as well as their response to apoptosis-inducing agents.

Conclusions and Significance

Characterizing and tracking changes in intracellular pathway profiles in cell subpopulations both at baseline and under therapeutic pressure will likely have important clinical applications, potentially informing the selection of beneficial targeted agents, used either alone or in combination with chemotherapy.