Structural analysis of DNA sequence: evidence for lateral gene transfer in Thermotoga maritima

The recently published complete DNA sequence of the bacterium Thermotoga maritima provides evidence, based on protein sequence conservation, for lateral gene transfer between Archaea and Bacteria. We introduce a new method of periodicity analysis of DNA sequences, based on structural parameters, which brings independent evidence for the lateral gene transfer in the genome of T.maritima. The structural analysis relates the Archaea-like DNA sequences to the genome of Pyrococcus horikoshii. Analysis of 24 complete genomic DNA sequences shows different periodicity patterns for organisms of different origin. The typical genomic periodicity for Bacteria is 11 bp whilst it is 10 bp for Archaea. Eukaryotes have more complex spectra but the dominant period in the yeast Saccharomyces cerevisiae is 10.2 bp. These periodicities are most likely reflective of differences in chromatin structure.     

DNA-based watermarks using the DNA-Crypt algorithm

Background  

The aim of this paper is to demonstrate the application of watermarks based on DNA sequences to identify the unauthorized use of genetically modified organisms (GMOs) protected by patents. Predicted mutations in the genome can be corrected by the DNA-Crypt program leaving the encrypted information intact. Existing DNA cryptographic and steganographic algorithms use synthetic DNA sequences to store binary information however, although these sequences can be used for authentication, they may change the target DNA sequence when introduced into living organisms.     

Identification of the Maize Gravitropism Gene lazy plant1 by a Transposon-Tagging Genome Resequencing Strategy

Since their initial discovery, transposons have been widely used as mutagens for forward and reverse genetic screens in a range of organisms. The problems of high copy number and sequence divergence among related transposons have often limited the efficiency at which tagged genes can be identified. A method was developed to identity the locations of Mutator (Mu) transposons in the Zea mays genome using a simple enrichment method combined with genome resequencing to identify transposon junction fragments. The sequencing library was prepared from genomic DNA by digesting with a restriction enzyme that cuts within a perfectly conserved motif of the Mu terminal inverted repeats (TIR). Paired-end reads containing Mu TIR sequences were computationally identified and chromosomal sequences flanking the transposon were mapped to the maize reference genome. This method has been used to identify Mu insertions in a number of alleles and to isolate the previously unidentified lazy plant1 (la1) gene. The la1 gene is required for the negatively gravitropic response of shoots and mutant plants lack the ability to sense gravity. Using bioinformatic and fluorescence microscopy approaches, we show that the la1 gene encodes a cell membrane and nuclear localized protein. Our Mu-Taq method is readily adaptable to identify the genomic locations of any insertion of a known sequence in any organism using any sequencing platform.     

Quartz crystal microbalance (QCM) affinity biosensor for genetically modified organisms (GMOs) detection

A DNA piezoelectric sensor has been developed for the detection of genetically modified organisms (GMOs). Single stranded DNA (ssDNA) probes were immobilised on the sensor surface of a quartz crystal microbalance (QCM) device and the hybridisation between the immobilised probe and the target complementary sequence in solution was monitored. The probe sequences were internal to the sequence of the 35S promoter (P) and Nos terminator (T), which are inserted sequences in the genome of GMOs regulating the transgene expression. Two different probe immobilisation procedures were applied: (a) a thiol-dextran procedure and (b) a thiol-derivatised probe and blocking thiol procedure. The system has been optimised using synthetic oligonucleotides, which were then applied to samples of plasmidic and genomic DNA isolated from the pBI121 plasmid, certified reference materials (CRM), and real samples amplified by the polymerase chain reaction (PCR). The analytical parameters of the sensor have been investigated (sensitivity, reproducibility, lifetime etc.). The results obtained showed that both immobilisation procedures enabled sensitive and specific detection of GMOs, providing a useful tool for screening analysis in food samples.     

QCM-based DNA biosensor for detection of genetically modified organisms (GMOs)

Development of a mass sensitive quartz crystal microbalance (QCM)-based DNA biosensor for the detection of the hybridization of CaMV 35S promoter sequence (P35S) was investigated for the screening of genetically modified organisms (GMOs). Attention was focused on the choice of the coating chemistry that could be used for the immobilization of probe sequences on the gold surface of the quartz crystal. Two immobilization procedures were tested and compared considering the amount of the immobilized P35S probe and the extent of the hybridization reaction with the target oligonucleotide. In wet chemistry procedure, the interaction between the thiol and gold for the immobilization of a thiolated probe was employed. Direct surface functionalization of piezoelectric quartz crystals were achieved in 13.56 MHz plasma polymerization reactor utilising ethylenediamine (EDA) precursors for the immobilization of amined probes. Results indicated that immobilization of a thiolated probe provides better immobilization characteristics and higher sensitivity for the detection of the hybridization reaction. The thiolated probe was used for the detection of P35S sequence in PCR-amplified DNAs and in real samples of pflp (ferrodoxin like protein)-gene inserted tobacco plants. Fragmentation of the genomic DNAs were achieved by digestion with restriction endonucleases and ultrasonication. The results obtained from the fragmented genomic DNAs demonstrated that it is possible to detect the target sequence directly in non-amplified genomic DNAs by using the developed QCM-based DNA biosensor system. The developed QCM-based DNA biosensor represented promising results for a real-time, label-free, direct detection of DNA samples for the screening of GMOs.     

Capturing genomic signatures of DNA sequence variation using a standard anonymous microarray platform

Comparative genomics, using the model organism approach, has provided powerful insights into the structure and evolution of whole genomes. Unfortunately, only a small fraction of Earth's biodiversity will have its genome sequenced in the foreseeable future. Most wild organisms have radically different life histories and evolutionary genomics than current model systems. A novel technique is needed to expand comparative genomics to a wider range of organisms. Here, we describe a novel approach using an anonymous DNA microarray platform that gathers genomic samples of sequence variation from any organism. Oligonucleotide probe sequences placed on a custom 44 K array were 25 bp long and designed using a simple set of criteria to maximize their complexity and dispersion in sequence probability space. Using whole genomic samples from three known genomes (mouse, rat and human) and one unknown (Gonystylus bancanus), we demonstrate and validate its power, reliability, transitivity and sensitivity. Using two separate statistical analyses, a large numbers of genomic ‘indicator’ probes were discovered. The construction of a genomic signature database based upon this technique would allow virtual comparisons and simple queries could generate optimal subsets of markers to be used in large-scale assays, using simple downstream techniques. Biologists from a wide range of fields, studying almost any organism, could efficiently perform genomic comparisons, at potentially any phylogenetic level after performing a small number of standardized DNA microarray hybridizations. Possibilities for refining and expanding the approach are discussed.     

metaSHARK: software for automated metabolic network prediction from DNA sequence and its application to the genomes of Plasmodium falciparum and Eimeria tenella

The metabolic SearcH And Reconstruction Kit (metaSHARK) is a new fully automated software package for the detection of enzyme-encoding genes within unannotated genome data and their visualization in the context of the surrounding metabolic network. The gene detection package (SHARKhunt) runs on a Linux system and requires only a set of raw DNA sequences (genomic, expressed sequence tag and/or genome survey sequence) as input. Its output may be uploaded to our web-based visualization tool (SHARKview) for exploring and comparing data from different organisms. We first demonstrate the utility of the software by comparing its results for the raw Plasmodium falciparum genome with the manual annotations available at the PlasmoDB and PlasmoCyc websites. We then apply SHARKhunt to the unannotated genome sequences of the coccidian parasite Eimeria tenella and observe that, at an E-value cut-off of 10⁻²⁰, our software makes 142 additional assertions of enzymatic function compared with a recent annotation package working with translated open reading frame sequences. The ability of the software to cope with low levels of sequence coverage is investigated by analyzing assemblies of the E.tenella genome at estimated coverages from 0.5× to 7.5×. Lastly, as an example of how metaSHARK can be used to evaluate the genomic evidence for specific metabolic pathways, we present a study of coenzyme A biosynthesis in P.falciparum and E.tenella.     

Long-range oscillation in a periodic DNA sequence motif may influence nucleosome array formation

We have experimentally examined the characteristics of nucleosome array formation in different regions of mouse liver chromatin, and have computationally analyzed the corresponding genomic DNA sequences. We have shown that the mouse adenosine deaminase (MADA) gene locus is packaged into an exceptionally regular nucleosome array with a shortened repeat, consistent with our computational prediction based on the DNA sequence. A survey of the mouse genome indicates that <10% of 70 kb windows possess a nucleosome-ordering signal, consisting of regular long-range oscillations in the period-10 triplet motif non-T, A/T, G (VWG), which is as strong as the signal in the MADA locus. A strong signal in the center of this locus, confirmed by in vitro chromatin assembly experiments, appears to cooperate with weaker, in-phase signals throughout the locus. In contrast, the mouse odorant receptor (MOR) locus, which lacks locus-wide signals, was representative of ~40% of the mouse genomic DNA surveyed. Within this locus, nucleosome arrays were similar to those of bulk chromatin. Genomic DNA sequences which were computationally similar to MADA or MOR resulted in MADA- or MOR-like nucleosome ladders experimentally. Overall, we provide evidence that computationally predictable information in the DNA sequence may affect nucleosome array formation in animal tissue.     

Transgene constructs in coho salmon (Oncorhynchus kisutch) are repeated in a head-to-tail fashion and can be integrated adjacent to horizontally-transmitted parasite DNA

Currently, little information is available regarding the molecular organization of integrated transgenes in genetically-engineered fish. We performed a detailed structural analysis of an inserted transgene in one strain (M77) of transgenic coho salmon (Oncorhynchus kisutch) containing a salmon growth hormone gene construct (OnMTGH1). Microinjected DNA was found to have inserted into a single site in the coho salmon genome, and was organized with four complete internal copies and two partial terminal copies of the OnMTGH1 construct. All construct copies were organized in a direct-tandem (head-to-tail) repeat fashion in strain M77 and five additional strains (one also possessed a second recombinant junction fragment). For strain M77, the junctions between the transgene insert and the insertion point within the wild-type genome were cloned from strain-specific cosmid libraries and sequenced, revealing that the transgene insertion was accompanied by a deletion of 587 bp of wild-type DNA as well as a small insertion (19 bp) of unknown DNA upstream and a 14 bp direct- tandem duplication of sequence downstream. Upstream and downstream wild-type DNA sequence contained several repetitive sequence elements based on Southern blot analysis and homology to repetitive sequences in GenBank. In the downstream flank, a pseudogene sequence was also identified which has high homology to the CA membrane protein gene from Schistosoma japonicum, a parasite closely related to Sanguinicola sp. parasites which infect salmonids. Whether the presence of an inserted transgene and the presence of potentially horizontally-transmitted DNA are indicative of a genomic region with a predisposition for insertion of foreign DNA requires further study. The information derived from this transgene structure provides information useful for comparison to other transgenic organisms and for determination of the mechanism of transgene integration in lower vertebrates.     

Characterization of Insertion Sites in Rainbow Papaya, the First Commercialized Transgenic Fruit Crop

Inserts and insert sites in transgenic, papaya ringspot virus (PRSV)-resistant commercial papaya Rainbow and SunUp, were characterized as part of a petition to Japan to allow import of fresh fruit of these cultivars from the U.S. and to provide data for a larger study aimed at understanding the global impact of DNA transformation on whole genome structure. The number and types of inserts were determined by Southern analysis using probes spanning the entire transformation plasmid and their sequences determined from corresponding clones or sequence reads from the whole-genome shotgun (WGS) sequence of SunUp papaya. All the functional transgenes, coding for the PRSV coat protein (CP), neophosphotransferase (nptII) and β-glucuronidase (uidA) were found in a single 9,789 basepair (bp) insert. Only two other inserts, one consisting of a 290 bp nonfunctional fragment of the nptII gene and a 1,533 bp plasmid-derived fragment containing a nonfunctional 222 bp segment of the tetA gene were detected in Rainbow and SunUp. Detection of the same three inserts in samples representing transgenic generations five to eight (R5 to R8) suggests that the three inserts are stably inherited. Five out of the six genomic DNA segments flanking the three inserts were nuclear plastid sequences (nupts). From the biosafety standpoint, no changes to endogenous gene function based on sequence structure of the transformation plasmid DNA insertion sites could be determined and no allergenic or toxic proteins were predicted from analysis of the insertion site and flanking genomic DNA.     

Efficient transfer of base changes from a vector to the rice genome by homologous recombination: involvement of heteroduplex formation and mismatch correction

Gene targeting refers to the alteration of a specific DNA sequence in an endogenous gene at its original locus in the genome by homologous recombination. Through a gene-targeting procedure with positive–negative selection, we previously reported the generation of fertile transgenic rice plants with a positive marker inserted into the Adh2 gene by using an Agrobacterium-mediated transformation vector containing the positive marker flanked by two 6-kb homologous segments for recombination. We describe here that base changes within the homologous segments in the vector could be efficiently transferred into the corresponding genomic sequences of rice recombinants. Interestingly, a few sequences from the host genome were flanked by the changed sequences derived from the vector in most of the recombinants. Because a single-stranded T-DNA molecule in Agrobacterium-mediated transformation is imported into the plant nucleus and becomes double-stranded, both single-stranded and double-stranded T-DNA intermediates can serve in gene-targeting processes. Several alternative models, including the occurrence of the mismatch correction of heteroduplex molecules formed between the genomic DNA and either a single-stranded or double-stranded T-DNA intermediate, are compared to explain the observation, and implications for the modification of endogenous genes for functional genomic analysis by gene targeting are discussed.     

Precise and Heritable Genome Editing in Evolutionarily Diverse Nematodes Using TALENs and CRISPR/Cas9 to Engineer Insertions and Deletions

Exploitation of custom-designed nucleases to induce DNA double-strand breaks (DSBs) at genomic locations of choice has transformed our ability to edit genomes, regardless of their complexity. DSBs can trigger either error-prone repair pathways that induce random mutations at the break sites or precise homology-directed repair pathways that generate specific insertions or deletions guided by exogenously supplied DNA. Prior editing strategies using site-specific nucleases to modify the Caenorhabditis elegans genome achieved only the heritable disruption of endogenous loci through random mutagenesis by error-prone repair. Here we report highly effective strategies using TALE nucleases and RNA-guided CRISPR/Cas9 nucleases to induce error-prone repair and homology-directed repair to create heritable, precise insertion, deletion, or substitution of specific DNA sequences at targeted endogenous loci. Our robust strategies are effective across nematode species diverged by 300 million years, including necromenic nematodes (Pristionchus pacificus), male/female species (Caenorhabditis species 9), and hermaphroditic species (C. elegans). Thus, genome-editing tools now exist to transform nonmodel nematode species into genetically tractable model organisms. We demonstrate the utility of our broadly applicable genome-editing strategies by creating reagents generally useful to the nematode community and reagents specifically designed to explore the mechanism and evolution of X chromosome dosage compensation. By developing an efficient pipeline involving germline injection of nuclease mRNAs and single-stranded DNA templates, we engineered precise, heritable nucleotide changes both close to and far from DSBs to gain or lose genetic function, to tag proteins made from endogenous genes, and to excise entire loci through targeted FLP-FRT recombination.     

Non-referenced genome assembly from epigenomic short-read data

Current computational methods used to analyze changes in DNA methylation and chromatin modification rely on sequenced genomes. Here we describe a pipeline for the detection of these changes from short-read sequence data that does not require a reference genome. Open source software packages were used for sequence assembly, alignment, and measurement of differential enrichment. The method was evaluated by comparing results with reference-based results showing a strong correlation between chromatin modification and gene expression. We then used our de novo sequence assembly to build the DNA methylation profile for the non-referenced Psammomys obesus genome. The pipeline described uses open source software for fast annotation and visualization of unreferenced genomic regions from short-read data.     

A new concept in (adenoviral) oncogenesis: integration of foreign DNA and its consequences

A new concept for viral oncogenesis is presented which is based on experimental work on the chromosomal integration of adenovirus DNA into mammalian genomes. The mechanism of adenovirus DNA integration is akin to non-sequence-specific insertional recombination in which patch homologies between the recombination partners are frequently observed. This reaction has been imitated in a cell-free system by using nuclear extracts from hamster cells and partly purified fractions derived from them. As a consequence of foreign DNA insertion into the mammalian genome, the foreign DNA is extensively de novo methylated in specific patterns, presumably as part of a mammalian host cell defense mechanism against inserted foreign DNA which can be permanently silenced in this way. A further corollary of foreign (adenovirus or bacteriophage λ) DNA integration is seen in extensive changes in cellular DNA methylation patterns at sites far remote from the locus of insertional recombination. Repetitive cellular, retrotransposon-like sequences are particularly, but not exclusively, prone to these increases in DNA methylation. It is conceivable that these changes in DNA methylation are a reflection of a profound overall reorganization process in the affected genomes. Could these alterations significantly contribute to the transformation events during viral or other types of oncogenesis? These sequelae of foreign DNA integration into established mammalian genomes will have to be critically considered when interpreting results obtained with transgenic, knock-out, and knock-in animals and when devising schemes for human somatic gene therapy.The interpretation of de novo methylation as a cellular defense mechanism has prompted investigations on the fate of food-ingested foreign DNA. The gastrointestinal (GI) tract provides a large surface for the entry of foreign DNA into any organism. As a tracer molecule, bacteriophage M13 DNA has been fed to mice. Fragments of this DNA can be found in small amounts (about 1 % of the administered DNA) in all parts of the intestinal tract and in the feces. Furthermore, M13 DNA can be traced in the columnar epithelia of the intestine, in Peyer's plaque leukocytes, in peripheral white blood cells, in spleen, and liver. Authentic M13 DNA has been recloned from total spleen DNA. If integrated, this DNA might elicit some of the described consequences of foreign DNA insertion into the mammalian genome. Food-ingested DNA will likely infiltrate the organism more frequently than viral DNA.     

Repertoires of the Nucleosome-Positioning Dinucleotides

It is generally accepted that the organization of eukaryotic DNA into chromatin is strongly governed by a code inherent in the genomic DNA sequence. This code, as well as other codes, is superposed on the triplets coding for amino acids. The history of the chromatin code started three decades ago with the discovery of the periodic appearance of certain dinucleotides, with AA/TT and RR/YY giving the strongest signals, all with a period of 10.4 bases. Every base-pair stack in the DNA duplex has specific deformation properties, thus favoring DNA bending in a specific direction. The appearance of the corresponding dinucleotide at the distance 10.4 xn bases will facilitate DNA bending in that direction, which corresponds to the minimum energy of DNA folding in the nucleosome. We have analyzed the periodic appearances of all 16 dinucleotides in the genomes of thirteen different eukaryotic organisms. Our data show that a large variety of dinucleotides (if not all) are, apparently, contributing to the nucleosome positioning code. The choice of the periodical dinucleotides differs considerably from one organism to another. Among other 10.4 base periodicities, a strong and very regular 10.4 base signal was observed for CG dinucleotides in the genome of the honey bee A. mellifera. Also, the dinucleotide CG appears as the only periodical component in the human genome. This observation seems especially relevant since CpG methylation is well known to modulate chromatin packing and regularity. Thus, the selection of the dinucleotides contributing to the chromatin code is species specific, and may differ from region to region, depending on the sequence context.     

Evolutionary force of AT‐rich repeats to trap genomic and episomal DNAs into the rice genome: lessons from endogenous pararetrovirus

In plant genomes, the incorporation of DNA segments is not a common method of artificial gene transfer. Nevertheless, various segments of pararetroviruses have been found in plant genomes in recent decades. The rice genome contains a number of segments of endogenous rice tungro bacilliform virus‐like sequences (ERTBVs), many of which are present between AT dinucleotide repeats (ATrs). Comparison of genomic sequences between two closely related rice subspecies, japonica and indica, allowed us to verify the preferential insertion of ERTBVs into ATrs. In addition to ERTBVs, the comparative analyses showed that ATrs occasionally incorporate repeat sequences including transposable elements, and a wide range of other sequences. Besides the known genomic sequences, the insertion sequences also represented DNAs of unclear origins together with ERTBVs, suggesting that ATrs have integrated episomal DNAs that would have been suspended in the nucleus. Such insertion DNAs might be trapped by ATrs in the genome in a host‐dependent manner. Conversely, other simple mono‐ and dinucleotide sequence repeats (SSR) were less frequently involved in insertion events relative to ATrs. Therefore, ATrs could be regarded as hot spots of double‐strand breaks that induce non‐homologous end joining. The insertions within ATrs occasionally generated new gene‐related sequences or involved structural modifications of existing genes. Likewise, in a comparison between Arabidopsis thaliana and Arabidopsis lyrata, the insertions preferred ATrs to other SSRs. Therefore ATrs in plant genomes could be considered as genomic dumping sites that have trapped various DNA molecules and may have exerted a powerful evolutionary force.     

Organization of shared and unshared sequences in the genomes of chicken endogenous and sarcoma viruses.

The genome of the genetically transmitted endogenous C type virus of chickens, RAV-O, is closely related to that of Rous sarcoma virus (RSV). Nevertheless, these viruses differ widely in oncogenicity and regulation by the host cell. Competitive hybridization analysis of ¹²⁵I-labeled genomic RNA demonstrated that the genome of RAV-O lacks about 35% of the sequences of nondefective RSV which formed hybrids with proviral DNA from RSV-infected cells, and that the genome of transformation-defective deletion mutants of RSV (td RSV) lacks about 15% of these sequences. Conversely, about 12% of the RAV-O sequences forming hybrids with normal chicken cell DNA were not detected in the sarcoma virus. A technique was developed to map the location of these unshared sequences by competitive hybridization. The deletion in the genome of td RSV was seen to begin at about 0.2 and to end at about 0.05 of the genome length from the 3′ end of sarcoma virus RNA, confirming the results of other laboratories using the method of mapping RNAase TI resistance of oligonucleotides. The 35% of RSV sequences missing and/or diverged in the genome of RAV-O were concentrated within 40% of the sarcoma virus genome from the 3′ end, and most of this large section did not appear to form hybrids with chicken DNA under the conditions of these experiments. A low level of hybrid formation was, however, detected between uninfected chicken cellular DNA and a small fraction of the nucleotides in the region of the td deletion. Analysis of RAV-O 3′ end fragments demonstrated that the genomic sequences of RAV-O missing in RSV were concentrated at the 3′ end of the endogenous viral genome. We conclude that the sequence differences between endogenous and sarcoma viruses are largely concentrated in specific regions of the viral genome.     

Transfection of plant mitochondria and in organello gene integration

Investigation and manipulation of mitochondrial genetics in animal and plant cells remains restricted by the lack of an efficient in vivo transformation methodology. Mitochondrial transfection in whole cells and maintenance of the transfected DNA are main issues on this track. We showed earlier that isolated mitochondria from different organisms can import DNA. Exploiting this mechanism, we assessed the possibility to maintain exogenous DNA in plant organelles. Whereas homologous recombination is scarce in the higher plant nuclear compartment, recombination between large repeats generates the multipartite structure of the plant mitochondrial genome. These processes are under strict surveillance to avoid extensive genomic rearrangements. Nevertheless, following transfection of isolated organelles with constructs composed of a partial gfp gene flanked by fragments of mitochondrial DNA, we demonstrated in organello homologous recombination of the imported DNA with the resident DNA and integration of the reporter gene. Recombination yielded insertion of a continuous exogenous DNA fragment including the gfp sequence and at least 0.5 kb of flanking sequence on each side. According to our observations, transfection constructs carrying multiple sequences homologous to the mitochondrial DNA should be suitable and targeting of most regions in the organelle genome should be feasible, making the approach of general interest.     

Comparative genomics: the key to understanding the Human Genome Project

The sequencing of the human genome is well underway. Technology has advanced, such that the total genomic sequence is possible, along with an extensive catalogue of genes via comprehensive cDNA libraries. With the recent completion of the Saccharomyces cerevisiae sequencing project and the imminent completion of that of Caenorhabditis elegans, the most frequently asked question is how much can sequence data alone tell us? The answer is that that a DNA sequence taken in isolation from a single organism reveals very little. The vast majority of DNA in most organisms is noncoding. Protein coding sequences or genes cannot function as isolated units without interaction with noncoding DNA and neighboring genes. This genomic environment is specific to each organism. In order to understand this we need to look at similar genes in different organisms, to determine how function and position has changed over the course of evolution. By understanding evolutionary processes we can gain a greater insight into what makes a gene and the wider processes of genetics and inheritance. Comparative genomics (with model organisms), once the poor relation of the human genome project, is starting to provide the key to unlock the DNA code.