Detection of G-Quadruplex DNA Using Primer Extension as a Tool

DNA sequence and structure play a key role in imparting fragility to different regions of the genome. Recent studies have shown that non-B DNA structures play a key role in causing genomic instability, apart from their physiological roles at telomeres and promoters. Structures such as G-quadruplexes, cruciforms, and triplexes have been implicated in making DNA susceptible to breakage, resulting in genomic rearrangements. Hence, techniques that aid in the easy identification of such non-B DNA motifs will prove to be very useful in determining factors responsible for genomic instability. In this study, we provide evidence for the use of primer extension as a sensitive and specific tool to detect such altered DNA structures. We have used the G-quadruplex motif, recently characterized at the BCL2 major breakpoint region as a proof of principle to demonstrate the advantages of the technique. Our results show that pause sites corresponding to the non-B DNA are specific, since they are absent when the G-quadruplex motif is mutated and their positions change in tandem with that of the primers. The efficiency of primer extension pause sites varied according to the concentration of monovalant cations tested, which support G-quadruplex formation. Overall, our results demonstrate that primer extension is a strong in vitro tool to detect non-B DNA structures such as G-quadruplex on a plasmid DNA, which can be further adapted to identify non-B DNA structures, even at the genomic level.

Binding site size limitations of imidazole-pyrrole polyamides for recognition in the minor groove of DNA

The discovery that the three ring polyamide Im-Py-Py-Dp containing imidazole (Im) and pyrrole (Py) carboxamides binds the DNA sequence 5'-(A,T)G(A,T)C(A,T)-3' as an antiparallel dimer offers a new model for the design of ligands for specific recognition of sequences in the minor groove containing both G,C and A,T base pairs. In Chapter 2, experiments are described in which the sequential addition of five N- methylpyrrolecarboxamides to the imidazole-pyrrole polyamide Im-Py-Py-Dp affords a series of six homologous polyamides, Im-(Py)_2-7-Dp, that differ in the size of their binding site, apparent first order binding affinity, and sequence specificity. These results demonstrate that DNA sequences up to nine base pairs in length can be specifically recognized by imidazole-pyrrole polyamides containing three to seven rings by 2:1 polyamide-DNA complex formation in the minor groove. Recognition of a nine base pair site defines the new lower limit of the binding site size that can be recognized by polyamides containing exclusively imidazole and pyrrolecarboxamides. The results of this study should provide useful guidelines for the design of new polyamides that bind longer DNA sites with enhanced affinity and specificity.

In Chapter 3 the design and synthesis of the hairpin polyamide Im-Py-Im-Py-γ-Im- Py-Im-Py-Dp is described. Quantitative DNase I footprint titration experiments reveal that Im-Py-Im-Py-γ-Im-Py-Im-Py-Dp binds six base pair 5'-(A,T)GCGC(A,T)-3' sequences with 30-fold higher affinity than the unlinked polyamide Im-Py-Im-Py-Dp. The hairpin polyamide does not discriminate between A•T and T•A at the first and sixth positions of the binding site as three sites 5'-TGCGCT-3', 5'-TGCGCA-3', and 5 'AGCGCT- 3' are bound with similar affinity. However, Im-Py-Im-Py-γ-Im-Py-Im-PyDp is specific for and discriminates between G•C and C•G base pairs in the 5'-GCGC-3' core as evidenced by lower affinities for the mismatched sites 5'-AACGCA-3', 5'- TGCGTT-3', 5'-TGCGGT-3', and 5'-ACCGCT-3'.

In Chapter 4, experiments are described in which a kinetically stable hexa-aza Schiff base La³⁺ complex is covalently attached to a Tat(49-72) peptide which has been shown to bind the HIV-1 TAR RNA sequence. Although these metallo-peptides cleave TAR site-specifically in the hexanucleotide loop to afford products consistent with hydrolysis, a series of control experiments suggests that the observed cleavage is not caused by a sequence-specifically bound Tat(49-72)-La(L)³⁺ peptide.

Cleavage of DNA by polyamide-seco-CBI conjugates

Small molecules that bind to any predetermined DNA sequence in the human genome are potentially useful tools for molecular biology and human medicine. Polyamides containing N-methylimidazole (Im) N-methylpyrrole (Py) are cell permeable small molecules that bind DNA according to a set of "pairing rules" with affinities and specificities similar to many naturally occurring DNA binding proteins. Py-Im polyamides offer a general approach to the chemical regulation of gene expression. We demonstrate here that polyamide containing a DNA alkylating moiety seco-CBI can specifically direct sequence specific DNA alkylation. We can also control the strand of DNA that is alkylated, depending on the enantiomer of seco-CBI used and the orientation of the polyamide relative to the alkylation site (Chapter 2). This class of molecules has been applied to a gene repair system in collaboration with the Baltimore group at Caltech (Chapter 3). Also reported are additional seco-CBI polyamide conjugates synthesized to study other systems (HIV-1 and COX-2) (Appendix 1).

Sequence-function relationships in E. coli transcriptional regulation

Understanding how transcriptional regulatory sequence maps to regulatory function remains a difficult problem in regulatory biology. Given a particular DNA sequence for a bacterial promoter region, we would like to be able to say which transcription factors bind there, how strongly they bind, and whether they interact with each other and/or RNA polymerase, with the ultimate objective of integrating knowledge of these parameters into a prediction of gene expression levels. The theoretical framework of statistical thermodynamics provides a useful framework for doing so, enabling us to predict how gene expression levels depend on transcription factor binding energies and concentrations. We used thermodynamic models, coupled with models of the sequence-dependent binding energies of transcription factors and RNAP, to construct a genotype to phenotype map for the level of repression exhibited by the lac promoter, and tested it experimentally using a set of promoter variants from E. coli strains isolated from different natural environments. For this work, we sought to ``reverse engineer'' naturally occurring promoter sequences to understand how variations in promoter sequence affects gene expression. The natural inverse of this approach is to ``forward engineer'' promoter sequences to obtain targeted levels of gene expression. We used a high precision model of RNAP-DNA sequence dependent binding energy, coupled with a thermodynamic model relating binding energy to gene expression, to predictively design and verify a suite of synthetic E. coli promoters whose expression varied over nearly three orders of magnitude.

However, although thermodynamic models enable predictions of mean levels of gene expression, it has become evident that cell-to-cell variability or ``noise'' in gene expression can also play a biologically important role. In order to address this aspect of gene regulation, we developed models based on the chemical master equation framework and used them to explore the noise properties of a number of common E. coli regulatory motifs; these properties included the dependence of the noise on parameters such as transcription factor binding strength and copy number. We then performed experiments in which these parameters were systematically varied and measured the level of variability using mRNA FISH. The results showed a clear dependence of the noise on these parameters, in accord with model predictions.

Finally, one shortcoming of the preceding modeling frameworks is that their applicability is largely limited to systems that are already well-characterized, such as the lac promoter. Motivated by this fact, we used a high throughput promoter mutagenesis assay called Sort-Seq to explore the completely uncharacterized transcriptional regulatory DNA of the E. coli mechanosensitive channel of large conductance (MscL). We identified several candidate transcription factor binding sites, and work is continuing to identify the associated proteins.

DNA-mediated oxidation of transcription factor p53

Transcription factor p53 is the most commonly altered gene in human cancer. As a redox-active protein in direct contact with DNA, p53 can directly sense oxidative stress through DNA-mediated charge transport. Electron hole transport occurs with a shallow distance dependence over long distances through the π-stacked DNA bases, leading to the oxidation and dissociation of DNA-bound p53. The extent of p53 dissociation depends upon the redox potential of the response element DNA in direct contact with each p53 monomer. The DNA sequence dependence of p53 oxidative dissociation was examined by electrophoretic mobility shift assays using radiolabeled oligonucleotides containing both synthetic and human p53 response elements with an appended anthraquinone photooxidant. Greater p53 dissociation is observed from DNA sequences containing low redox potential purine regions, particularly guanine triplets, within the p53 response element. Using denaturing polyacrylamide gel electrophoresis of irradiated anthraquinone-modified DNA, the DNA damage sites, which correspond to locations of preferred electron hole localization, were determined. The resulting DNA damage preferentially localizes to guanine doublets and triplets within the response element. Oxidative DNA damage is inhibited in the presence of p53, however, only at DNA sites within the response element, and therefore in direct contact with p53. From these data, predictions about the sensitivity of human p53-binding sites to oxidative stress, as well as possible biological implications, have been made. On the basis of our data, the guanine pattern within the purine region of each p53-binding site determines the response of p53 to DNA-mediated oxidation, yielding for some sequences the oxidative dissociation of p53 from a distance and thereby providing another potential role for DNA charge transport chemistry within the cell.

To determine whether the change in p53 response element occupancy observed in vitro also correlates in cellulo, chromatin immunoprecipition (ChIP) and quantitative PCR (qPCR) were used to directly quantify p53 binding to certain response elements in HCT116N cells. The HCT116N cells containing a wild type p53 were treated with the photooxidant [Rh(phi)2bpy]³⁺, Nutlin-3 to upregulate p53, and subsequently irradiated to induce oxidative genomic stress. To covalently tether p53 interacting with DNA, the cells were fixed with disuccinimidyl glutarate and formaldehyde. The nuclei of the harvested cells were isolated, sonicated, and immunoprecipitated using magnetic beads conjugated with a monoclonal p53 antibody. The purified immounoprecipiated DNA was then quantified via qPCR and genomic sequencing. Overall, the ChIP results were significantly varied over ten experimental trials, but one trend is observed overall: greater variation of p53 occupancy is observed in response elements from which oxidative dissociation would be expected, while significantly less change in p53 occupancy occurs for response elements from which oxidative dissociation would not be anticipated.

The chemical oxidation of transcription factor p53 via DNA CT was also investigated with respect to the protein at the amino acid level. Transcription factor p53 plays a critical role in the cellular response to stress stimuli, which may be modulated through the redox modulation of conserved cysteine residues within the DNA-binding domain. Residues within p53 that enable oxidative dissociation are herein investigated. Of the 8 mutants studied by electrophoretic mobility shift assay (EMSA), only the C275S mutation significantly decreased the protein affinity (KD) for the Gadd45 response element. EMSA assays of p53 oxidative dissociation promoted by photoexcitation of anthraquinone-tethered Gadd45 oligonucleotides were used to determine the influence of p53 mutations on oxidative dissociation; mutation to C275S severely attenuates oxidative dissociation while C277S substantially attenuates dissociation. Differential thiol labeling was used to determine the oxidation states of cysteine residues within p53 after DNA-mediated oxidation. Reduced cysteines were iodoacetamide labeled, while oxidized cysteines participating in disulfide bonds were ¹³C₂D₂-iodoacetamide labeled. Intensities of respective iodoacetamide-modified peptide fragments were analyzed using a QTRAP 6500 LC-MS/MS system, quantified with Skyline, and directly compared. A distinct shift in peptide labeling toward ¹³C₂D₂-iodoacetamide labeled cysteines is observed in oxidized samples as compared to the respective controls. All of the observable cysteine residues trend toward the heavy label under conditions of DNA CT, indicating the formation of multiple disulfide bonds potentially among the C124, C135, C141, C182, C275, and C277. Based on these data it is proposed that disulfide formation involving C275 is critical for inducing oxidative dissociation of p53 from DNA.

Investigations of Pyrrole-Imidazole Polyamide Effects on DNA Replication

Pyrrole–Imidazole polyamides are programmable, cell-permeable small molecules that bind in the minor groove of double-stranded DNA sequence-specifically. Polyamide binding has been shown to alter the local helical structure of DNA, disrupt protein-DNA interactions, and modulate endogenous gene expression. Py–Im polyamides targeted to the androgen receptor-DNA interface have been observed to decrease expression of androgen-regulated genes, upregulate p53, and induce apoptosis in a hormone-sensitive prostate cancer cell line. Here we report that androgen response element (ARE)-targeted polyamides induced DNA replication stress in a hormone-insensitive prostate cancer cell line. The ATR checkpoint kinase was activated in response to this stress, causing phosphorylation of MCM2, and FANCD2 was monoubiquitinated. Surprisingly, little single-stranded DNA was exhibited, and the ATR targets RPA2 and Chk1 were not phosphorylated. We conclude that polyamide induces relatively low level replication stress, and suggest inhibition of the replicative helicase as a putative mechanism based on in vitro assays. We also demonstrate polyamide-induced inhibition of DNA replication in cell free extracts from X. laevis oocytes. In this system, inhibition of chromatin decondensation is observed, preventing DNA replication initiation. Finally, we show that Py-Im polyamides targeted to the ARE and ETS binding sequence downregulate AR- and ERG-driven signaling in a prostate cancer cell line harboring the TMPRSS2-ERG fusion. In a mouse xenograft model, ARE-targeted polyamide treatment reduced growth of the tumor.

Differenzierung von Kammmuscheln durch DNA-Analyse

Abstract In the last years scallops have reached a considerable popularity and the import of scallops into the EU has increased about 20 % over the last fi ve years from some 50.000 t to nearly 63.000 t in the year 2010. Scallops are fi shed or farmed, and traded as fresh or deep frozen product. Recently investigation of scallop products of various origins by determining the species using molecular biological techniques showed that the species had been mislabelled in a considerable proportion of samples. Determination of the species was performed by PCR-based DNA-analysis of mitochondrial DNA followed by (i) sequencing the PCR product and (ii) comparison of the DNA sequence with entries in GenBank using BLAST. The deduced sequences of the analysed samples were considerably different from each other allowing the unambiguous assignment of samples to a certain species. Kurzfassung Die Nachfrage von Kammmuscheln in der EU hat in den letzten fünf Jahren erheblich zugenommen. Der Import stieg von knapp 53.000 t im Jahr 2005 um 20% auf annähernd 63.000 t im Jahr 2010. Gehandelt werden Kammmuscheln sowohl als frische als auch als Tiefkühlware aus Wildfängen und Aquakultur. Untersuchungen von Kammmuschel-Proben aus verschiedenen Ursprungsländern und Bestimmung der Spezies auf molekularbiologischer Basis zeigten, dass ein erheblicher Anteil der Proben falsch deklariert war. Die Bestimmung der Spezies erfolgte durch Vervielfältigung eines Abschnitts des 16S rRNA Gens durch Polymerase- Kettenreaktion (PCR), anschließender Sequenzanalyse der PCR-Produkte und Vergleich der DNA Sequenzen untereinander und mit Dateneintragungen in GenBank. Die DNA-Sequenzen der ermittelten Abschnitte der 16S rRNA der Proben unterschieden sich erheblich voneinander und erlaubten eine eindeutige Zuordnung zu jeweils einer Spezies.

Hardware paralelo reconfigurável para identificação de alinhamentos de sequências de DNA.

Amostras de DNA são encontradas em fragmentos, obtidos em vestígios de uma cena de crime, ou coletados de amostras de cabelo ou sangue, para testes genéticos ou de paternidade. Para identificar se esse fragmento pertence ou não a uma sequência de DNA, é necessário compará-los com uma sequência determinada, que pode estar armazenada em um banco de dados para, por exemplo, apontar um suspeito. Para tal, é preciso uma ferramenta eficiente para realizar o alinhamento da sequência de DNA encontrada com a armazenada no banco de dados. O alinhamento de sequências de DNA, em inglês DNA matching, é o campo da bioinformática que tenta entender a relação entre as sequências genéticas e suas relações funcionais e parentais. Essa tarefa é frequentemente realizada através de softwares que varrem clusters de base de dados, demandando alto poder computacional, o que encarece o custo de um projeto de alinhamento de sequências de DNA. Esta dissertação apresenta uma arquitetura de hardware paralela, para o algoritmo BLAST, que permite o alinhamento de um par de sequências de DNA. O algoritmo BLAST é um método heurístico e atualmente é o mais rápido. A estratégia do BLAST é dividir as sequências originais em subsequências menores de tamanho w. Após realizar as comparações nessas pequenas subsequências, as etapas do BLAST analisam apenas as subsequências que forem idênticas. Com isso, o algoritmo diminui o número de testes e combinações necessárias para realizar o alinhamento. Para cada sequência idêntica há três etapas, a serem realizadas pelo algoritmo: semeadura, extensão e avaliação. A solução proposta se inspira nas características do algoritmo para implementar um hardware totalmente paralelo e com pipeline entre as etapas básicas do BLAST. A arquitetura de hardware proposta foi implementada em FPGA e os resultados obtidos mostram a comparação entre área ocupada, número de ciclos e máxima frequência de operação permitida, em função dos parâmetros de alinhamento. O resultado é uma arquitetura de hardware em lógica reconfigurável, escalável, eficiente e de baixo custo, capaz de alinhar pares de sequências utilizando o algoritmo BLAST.

Mitochondrial DNA variation, effective female population size and population history of the endangered Chinese sturgeon, Acipenser sinensis

The anadromous Chinese sturgeon (Acipenser sinensis), mainly endemic to the Yangtze River in China, is an endangered fish species. The natural population has declined since the Gezhouba Dam blocked its migratory route to the spawning grounds in 1981. In the near future, the completion of the Three Gorges Dam, the world's largest hydroelectric project, may further impact this species by altering the water flow of the Yangtze River. Little is currently known about the population genetic structure of the Chinese sturgeon. In this study, DNA sequence data were determined from the control region (D-loop) of the mitochondrial genome of adult sturgeons (n = 106) that were collected between 1995-2000. The molecular data were used to investigate genetic variation, effective female population size and population history of the Chinese sturgeon in the Yangtze River. Our results indicate that the reduction in abundance did not change genetic variation of the Chinese sturgeon, and that the population underwent an expansion in the past. AMOVA analysis indicated that 98.7% of the genetic variability occurred within each year's spawning populations, the year of collection had little influence on the diversity of annual temporary samples. The relative large effective female population size (N-ef) indicates that good potential exists for the recovery of this species in the future. Strikingly, the ratio of N-ef to the census female population size (N-f) is unusually high (0.77-0.93). This may be the result of a current bottleneck in the population of the Chinese sturgeon that is likely caused by human intervention.

BarraCUDA - a fast short read sequence aligner using graphics processing units.

BACKGROUND: With the maturation of next-generation DNA sequencing (NGS) technologies, the throughput of DNA sequencing reads has soared to over 600 gigabases from a single instrument run. General purpose computing on graphics processing units (GPGPU), extracts the computing power from hundreds of parallel stream processors within graphics processing cores and provides a cost-effective and energy efficient alternative to traditional high-performance computing (HPC) clusters. In this article, we describe the implementation of BarraCUDA, a GPGPU sequence alignment software that is based on BWA, to accelerate the alignment of sequencing reads generated by these instruments to a reference DNA sequence. FINDINGS: Using the NVIDIA Compute Unified Device Architecture (CUDA) software development environment, we ported the most computational-intensive alignment component of BWA to GPU to take advantage of the massive parallelism. As a result, BarraCUDA offers a magnitude of performance boost in alignment throughput when compared to a CPU core while delivering the same level of alignment fidelity. The software is also capable of supporting multiple CUDA devices in parallel to further accelerate the alignment throughput. CONCLUSIONS: BarraCUDA is designed to take advantage of the parallelism of GPU to accelerate the alignment of millions of sequencing reads generated by NGS instruments. By doing this, we could, at least in part streamline the current bioinformatics pipeline such that the wider scientific community could benefit from the sequencing technology.BarraCUDA is currently available from http://seqbarracuda.sf.net.

Cloning and analysis of 16 Rab genes from macronuclear DNA of Euplotes octocarinatus

Rab proteins belong to the largest family of the Ras superfamily of small GTPase that play an important role in intracellular vesicular traffic. So far, almost 60 members of Rab family have been identified in mammalian cells. To further study the diversity and function of Rab protein in evolution, unicellular protozoa ciliates, Euplotes octocarinatus, were used in this study, Rab genes were screened by PCR method from macronuclear DNA of E. octocarinatus. Sixteen Rab genes were obtained. They share 87.6 - 99.5% identities. Highly conserved GTP-binding domains were found. There are some hot regions that diverse sharply in these genes as well.

DNA diagnosis by capillary electrophoresis and microfabricated electrophoretic devices

DNA diagnosis is experiencing an impressive progression towards the development of novel technology to identity various clinically relevant categories of genetic changes and to meet the exponential growth of genomics. The introduction of capillary electrophoresis has dramatically accelerated the completion of the first draft of the human DNA sequence in the Human Genome Project, and thus, has become the method of choice for analysis of various genetic variants. The recent development of microfabricated electrophoretic devices has led to the possibility of integrating multiple sample handling with the actual measurement steps required for automation of molecular diagnostics. This review highlights the most recent progress in capillary electrophoresis and electrophoretic microdevices for DNA-based diagnostics, including the important areas of genotyping for point mutation, single nucleotide polymorphisms, short tandem repeats and organism identification. The application of these techniques for infectious and genetic disease diagnosis, as well as forensic identification purpose, are covered. The promising development and the challenges for techinical problems are also discussed.

attB site disruption in marine Actinomyces sp M048 via DNA transformation of a site-specific integration vector

An efficient conjugation method has been developed for the marine Actinomyces sp. isolate M048 to facilitate the genetic manipulation of the chandrananimycin biosynthesis gene cluster. A phi C31-derived integration vector pIJ8600 containing oriT and attP fragments was introduced into strain M048 by bi-parental conjugation from Escherichia coli ET12567 to strain M048. Transformation efficiency was (6.38 +/- 0.41) x 10(-5) exconjugants per recipient spore. Analysis of eight exconjugants showed that the plasmid pIJ8600 was stably integrated at a single chromosomal site (attB) of the Actinomyces genome. The DNA sequence of the attB was cloned and shown to be conserved. The results of antimicrobial activity analysis indicated that the insertion of plasmid pIJ8600 seemed to affect the biosynthesis of antibiotics that could strongly inhibit the growth of E. coli and Mucor miehei (Tu284). HPLC-MS analysis of the extracts indicated that disruption of the attB site resulted in the complete abolition of chandrananimycin A-C production, proving the identity of the gene cluster. Instead of chandrananimycins, two bafilomycins were produced through disruption of the attB site from the chromosomal DNA of marine Actinomyces sp. M048.

Phylogeography of Pinus tabulaeformis Carr. (Pinaceae), a dominant species of coniferous forest in northern China

How coniferous trees in northern China changed their distribution ranges in response to Quaternary climatic oscillations remains largely unknown. Here we report a study of the phylogeography of Pinus tabulaeformis, an endemic and dominant species of coniferous forest in northern China. We examined sequence variation of maternally inherited, seed-dispersed mitochondrial DNA (mtDNA) (nad5 intron 1 and nad4/3-4) and paternally inherited, pollen- and seed-dispersed chloroplast DNA (cpDNA) (rpl16 and trnS-trnG) within and among 30 natural populations across the entire range of the species. Six mitotypes and five chlorotypes were recovered among 291 trees surveyed. Population divergence was high for mtDNA variation (G(ST) = 0.738, N-ST = 0.771) indicating low levels of seed-based gene flow and significant phylogeographical structure (N-ST > G(ST), P < 0.05). The spatial distribution of mitotypes suggests that five distinct population groups exist in the species: one in the west comprising seven populations, a second with a north-central distribution comprising 15 populations, a third with a southern and easterly distribution comprising five populations, a fourth comprising one central and one western population, and a fifth comprising a single population located in the north-central part of the species' range. Each group apart from the fourth group is characterized by a distinct mitotype, with other mitotypes, if present, occurring at low frequency. It is suggested, therefore, that most members of each group apart from Group 4 are derived from ancestors that occupied different isolated refugia in a previous period of range fragmentation of the species, possibly at the time of the Last Glacial Maximum. Possible locations for these refugia are suggested. A comparison of mitotype diversity between northern and southern subgroups within the north-central group of populations (Group 2) showed much greater uniformity in the northern part of the range both within and between populations. This could indicate a northward migration of the species from a southern refugium in this region during the postglacial period, although alternative explanations cannot be ruled out. Two chlorotypes were distributed across the geographical range of the species, resulting in lower levels of among-population chlorotype variation. The geographical pattern of variation for all five chlorotypes provided some indication of the species surviving past glaciations in more than one refugium, although differentiation was much less marked, presumably due to the greater dispersal of cpDNA via pollen.