Genomic markers on synthetic genomes

Genome synthesis endows scientists the ability of de novo creating genomes absent in nature, by thorough redesigning DNA sequences and introducing numerous custom features. However, the genome synthesis is a labor‐ and time‐consuming work, and thus it is a challenge to verify and quantify the synthetic genome rapidly and precisely. Thus, specific DNA sequences different from native genomic sequences are designed into synthetic genomes during synthesis, namely genomic markers. Genomic markers can be easily detected by PCR reaction, whole‐genome sequencing (WGS) and a variety of methods to identify the synthetic genome from native one. Here, we review types and applications of genomic markers utilized in synthetic genomes, with the hope of providing a guidance for future works.     

Selective Whole Genome Amplification for Resequencing Target Microbial Species from Complex Natural Samples

Population genomic analyses have demonstrated power to address major questions in evolutionary and molecular microbiology. Collecting populations of genomes is hindered in many microbial species by the absence of a cost effective and practical method to collect ample quantities of sufficiently pure genomic DNA for next-generation sequencing. Here we present a simple method to amplify genomes of a target microbial species present in a complex, natural sample. The selective whole genome amplification (SWGA) technique amplifies target genomes using nucleotide sequence motifs that are common in the target microbe genome, but rare in the background genomes, to prime the highly processive phi29 polymerase. SWGA thus selectively amplifies the target genome from samples in which it originally represented a minor fraction of the total DNA. The post-SWGA samples are enriched in target genomic DNA, which are ideal for population resequencing. We demonstrate the efficacy of SWGA using both laboratory-prepared mixtures of cultured microbes as well as a natural host–microbe association. Targeted amplification of Borrelia burgdorferi mixed with Escherichia coli at genome ratios of 1:2000 resulted in >10⁵-fold amplification of the target genomes with <6.7-fold amplification of the background. SWGA-treated genomic extracts from Wolbachia pipientis-infected Drosophila melanogaster resulted in up to 70% of high-throughput resequencing reads mapping to the W. pipientis genome. By contrast, 2–9% of sequencing reads were derived from W. pipientis without prior amplification. The SWGA technique results in high sequencing coverage at a fraction of the sequencing effort, thus allowing population genomic studies at affordable costs.     

Improving specificity of DNA hybridization-based methods

Methods based on DNA reassociation in solution with the subsequent PCR amplification of certain hybrid molecules, such as coincidence cloning and subtractive hybridization, all suffer from a common imperfection: cross-hybridization between various types of paralogous repetitive DNA fragments. Although the situation can be slightly improved by the addition of repeat-specific competitor DNA into the hybridization mixture, the cross-hybridization outcome is a significant number of background chimeric clones in resulting DNA libraries. In order to overcome this challenge, we developed a technique called mispaired DNA rejection (MDR), which utilizes a treatment of resulting reassociated DNA with mismatch-specific nucleases. We examined the MDR efficiency using cross-hybridization of complex, whole genomic mixtures derived from human and chimpanzee genomes, digested with frequent-cutter restriction enzyme. We show here that both single-stranded DNA-specific and mismatched double-stranded DNA-specific nucleases can be used for MDR separately or in combination, reducing the background level from 60 to 4% or lower. The technique presented here is of universal usefulness and can be applied to both cDNA and genomic DNA subtractions of very complex DNA mixtures. MDR is also useful for the genome-wide recovery of highly conserved DNA sequences, as we demonstrate by comparing human and pygmy marmoset genomes.     

Massively parallel characterization of restriction endonucleases

Restriction endonucleases are highly specific in recognizing the particular DNA sequence they act on. However, their activity is affected by sequence context, enzyme concentration and buffer composition. Changes in these factors may lead to either ineffective cleavage at the cognate restriction site or relaxed specificity allowing cleavage of degenerate ‘star’ sites. Additionally, uncharacterized restriction endonucleases and engineered variants present novel activities. Traditionally, restriction endonuclease activity is assayed on simple substrates such as plasmids and synthesized oligonucleotides. We present and use high-throughput Illumina sequencing-based strategies to assay the sequence specificity and flanking sequence preference of restriction endonucleases. The techniques use fragmented DNA from sequenced genomes to quantify restriction endonuclease cleavage on a complex genomic DNA substrate in a single reaction. By mapping millions of restriction site–flanking reads back to the Escherichia coli and Drosophila melanogaster genomes we were able to quantitatively characterize the cognate and star site activity of EcoRI and MfeI and demonstrate genome-wide decreases in star activity with engineered high-fidelity variants EcoRI-HF and MfeI-HF, as well as quantify the influence on MfeI cleavage conferred by flanking nucleotides. The methods presented are readily applicable to all type II restriction endonucleases that cleave both strands of double-stranded DNA.     

Analyzing the forces binding a restriction endonuclease to DNA using a synthetic nanopore

Restriction endonucleases are used prevalently in recombinant DNA technology because they bind so stably to a specific target sequence and, in the presence of cofactors, cleave double-helical DNA specifically at a target sequence at a high rate. Using synthetic nanopores along with molecular dynamics (MD), we have analyzed with atomic resolution how a prototypical restriction endonuclease, EcoRI, binds to the DNA target sequence—GAATTC—in the absence of a Mg²⁺ ion cofactor. We have previously shown that there is a voltage threshold for permeation of DNA bound to restriction enzymes through a nanopore that is associated with a nanonewton force required to rupture the complex. By introducing mutations in the DNA, we now show that this threshold depends on the recognition sequence and scales linearly with the dissociation energy, independent of the pore geometry. To predict the effect of mutation in a base pair on the free energy of dissociation, MD is used to qualitatively rank the stability of bonds in the EcoRI–DNA complex. We find that the second base in the target sequence exhibits the strongest binding to the protein, followed by the third and first bases, with even the flanking sequence affecting the binding, corroborating our experiments.     

Methylglyoxal, an endogenous aldehyde, crosslinks DNA polymerase and the substrate DNA

Methylglyoxal, a known endogenous and environmental mutagen, is a reactive α-ketoaldehyde that can modify both DNA and proteins. To investigate the possibility that methylglyoxal induces a crosslink between DNA and DNA polymerase, we treated a ‘primed template’ DNA and the exonuclease-deficient Klenow fragment (KF^exo–) of DNA polymerase I with methylglyoxal in vitro. When the reaction mixtures were analyzed by SDS–PAGE, we found that methylglyoxal induced a DNA–KF^exo– crosslink. The specific binding complex of KF^exo– and ‘primed template’ DNA was necessary for formation of the DNA–KF^exo– crosslink. Methylglyoxal reacted with guanine residues in the single-stranded portion of the template DNA. When 2′-deoxyguanosine was incubated with Nα-acetyllysine or N-acetylcysteine in the presence of methylglyoxal, a crosslinked product was formed. No other amino acid derivatives tested could generate a crosslinked product. These results suggest that methylglyoxal crosslinks a guanine residue of the substrate DNA and lysine and cysteine residues near the binding site of the DNA polymerase during DNA synthesis and that DNA replication is severely inhibited by the methylglyoxal-induced DNA–DNA polymerase crosslink.     

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.     

A unique family of Mrr-like modification-dependent restriction endonucleases

Mrr superfamily of homologous genes in microbial genomes restricts modified DNA in vivo. However, their biochemical properties in vitro have remained obscure. Here, we report the experimental characterization of MspJI, a remote homolog of Escherichia coli’s Mrr and show it is a DNA modification-dependent restriction endonuclease. Our results suggest MspJI recognizes ^mCNNR (R = G/A) sites and cleaves DNA at fixed distances (N₁₂/N₁₆) away from the modified cytosine at the 3′ side (or N₉/N₁₃ from R). Besides 5-methylcytosine, MspJI also recognizes 5-hydroxymethylcytosine but is blocked by 5-glucosylhydroxymethylcytosine. Several other close homologs of MspJI show similar modification-dependent endonuclease activity and display substrate preferences different from MspJI. A unique feature of these modification-dependent enzymes is that they are able to extract small DNA fragments containing modified sites on genomic DNA, for example ∼32 bp around symmetrically methylated CG sites and ∼31 bp around methylated CNG sites. The digested fragments can be directly selected for high-throughput sequencing to map the location of the modification on the genomic DNA. The MspJI enzyme family, with their different recognition specificities and cleavage properties, provides a basis on which many future methods can build to decode the epigenomes of different organisms.     

Crystal structure and mechanism of action of the N6-methyladenine-dependent type IIM restriction endonuclease R.DpnI

DNA methylation-dependent restriction enzymes have many applications in genetic engineering and in the analysis of the epigenetic state of eukaryotic genomes. Nevertheless, high-resolution structures have not yet been reported, and therefore mechanisms of DNA methylation-dependent cleavage are not understood. Here, we present a biochemical analysis and high-resolution DNA co-crystal structure of the N(6)-methyladenine (m6A)-dependent restriction enzyme R.DpnI. Our data show that R.DpnI consists of an N-terminal catalytic PD-(D/E)XK domain and a C-terminal winged helix (wH) domain. Surprisingly, both domains bind DNA in a sequence- and methylation-sensitive manner. The crystal contains R.DpnI with fully methylated target DNA bound to the wH domain, but distant from the catalytic domain. Independent readout of DNA sequence and methylation by the two domains might contribute to R.DpnI specificity or could help the monomeric enzyme to cut the second strand after introducing a nick.     

Functional analysis of the homeodomain protein SIX5

SIX5 (previously known as myotonic dystrophy associated homeodomain protein-DMAHP) is a member of the SIX [sine oculis homeobox (Drosophila) homologue] gene family which encodes proteins containing a SIX domain adjacent to a homeodomain. To investigate the DNA binding specificities of these two domains in SIX5, they were expressed as GST fusion proteins, both separately and together. Affinity purified recombinant proteins and cell lysates from bacteria expressing the recombinant proteins were used in gel retardation assays with double stranded oligonucleotides representing putative DNA binding sites. The putative sites included two in the promoter region of DMPK (dystrophia myotonica protein kinase) and the previously characterised murine Six4 DNA binding site in the Na⁺/K⁺ ATPase α1 subunit gene (ATP1A1) regulatory element (ARE). None of the recombinant proteins showed any affinity for the two putative sites in DMPK. However, the two recombinant proteins containing the homeodomain both formed at least one specific complex with the ARE. The recombinant protein containing both domains formed a second specific complex with the ARE, assumed to be a dimer complex. Finally, a whole genome PCR-based screen was used to identify genomic DNA sequences to which SIX5 binds, as an initial stage in the identification of genes regulated by SIX5.     

Trapping DNA-protein binding reactions with neutral osmolytes for the analysis by gel mobility shift and self-cleavage assays

We take advantage of our previous observation that neutral osmolytes can strongly slow down the rate of DNA–protein complex dissociation to develop a method that uses osmotic stress to ‘freeze’ mixtures of DNA–protein complexes and prevent further reaction enabling analysis of the products. We apply this approach to the gel mobility shift assay and use it to modify a self-cleavage assay that uses the nuclease activity of the restriction endonucleases to measure sensitively their specific binding to DNA. At sufficiently high concentrations of neutral osmolytes the cleavage reaction can be triggered at only those DNA fragments with initially bound enzyme. The self-cleavage assay allows measurement of binding equilibrium and kinetics directly in solution avoiding the intrinsic problems of gel mobility shift and filter binding assays while providing the same sensitivity level. Here we compare the self-cleavage and gel mobility shift assays applied to the DNA binding of EcoRI and BamHI restriction endonucleases. Initial results indicate that BamHI dissociation from its specific DNA sequence is strongly linked to water activity with the half-life time of the specific complex increasing ~20-fold from 0 to 1 osmolal betaine.     

A Method for Selectively Enriching Microbial DNA from Contaminating Vertebrate Host DNA

DNA samples derived from vertebrate skin, bodily cavities and body fluids contain both host and microbial DNA; the latter often present as a minor component. Consequently, DNA sequencing of a microbiome sample frequently yields reads originating from the microbe(s) of interest, but with a vast excess of host genome-derived reads. In this study, we used a methyl-CpG binding domain (MBD) to separate methylated host DNA from microbial DNA based on differences in CpG methylation density. MBD fused to the Fc region of a human antibody (MBD-Fc) binds strongly to protein A paramagnetic beads, forming an effective one-step enrichment complex that was used to remove human or fish host DNA from bacterial and protistan DNA for subsequent sequencing and analysis. We report enrichment of DNA samples from human saliva, human blood, a mock malaria-infected blood sample and a black molly fish. When reads were mapped to reference genomes, sequence reads aligning to host genomes decreased 50-fold, while bacterial and Plasmodium DNA sequences reads increased 8–11.5-fold. The Shannon-Wiener diversity index was calculated for 149 bacterial species in saliva before and after enrichment. Unenriched saliva had an index of 4.72, while the enriched sample had an index of 4.80. The similarity of these indices demonstrates that bacterial species diversity and relative phylotype abundance remain conserved in enriched samples. Enrichment using the MBD-Fc method holds promise for targeted microbiome sequence analysis across a broad range of sample types.     

Haplotype-Phased Synthetic Long Reads from Short-Read Sequencing

Next-generation DNA sequencing has revolutionized the study of biology. However, the short read lengths of the dominant instruments complicate assembly of complex genomes and haplotype phasing of mixtures of similar sequences. Here we demonstrate a method to reconstruct the sequences of individual nucleic acid molecules up to 11.6 kilobases in length from short (150-bp) reads. We show that our method can construct 99.97%-accurate synthetic reads from bacterial, plant, and animal genomic samples, full-length mRNA sequences from human cancer cell lines, and individual HIV env gene variants from a mixture. The preparation of multiple samples can be multiplexed into a single tube, further reducing effort and cost relative to competing approaches. Our approach generates sequencing libraries in three days from less than one microgram of DNA in a single-tube format without custom equipment or specialized expertise.     

Direct cloning of microsatellite loci from Cottus gobio through a simplified enrichment procedure

Applying a simplified enrichment procedure, we have isolated and characterized 177 microsatellite markers for Cottus gobio L., 1758 (Cottidae, Scorpaeniformes, Teleostei). In contrast to using specific probes for the enrichment, we use genomic DNAs of unrelated organisms for cross‐hybridization. This takes advantage of the fact that simple sequences are the only repetitive elements that are abundantly found in all eukaryotic genomes and that any genome usually contains all permutations of microsatellite motifs. This cross‐hybridization principle was employed to enrich genomic libraries of Cottus DNA to obtain a large number of nonredundant microsatellite markers without further screening procedures.     

Recognition of native DNA methylation by the PvuII restriction endonuclease

Recognizing the methylation status of specific DNA sequences is central to the function of many systems in eukaryotes and prokaryotes. Restriction–modification systems have to distinguish between ‘self’ and ‘non-self’ DNA and depend on the inability of restriction endonucleases to cleave their DNA substrates when the DNA is appropriately methylated. These endonucleases thus provide a model system for studying the recognition of DNA methylation by proteins. We have characterized the interaction of R·PvuII with DNA containing the physiologically relevant N4-methylcytosine modification. R·PvuII binds ^N4mC-modified DNA and cleaves it very slowly. Methylated strands in hemimethylated duplexes were cleaved at a higher rate than in fully methylated duplexes, in parallel with a higher binding affinity for hemimethylated DNA. The co-crystal structures of R·PvuII–DNA, together with a mutagenesis study, have implicated specific amino acids in recognition of the methylatable base; one of these is His84. We report that replacing His84 with Ala reduced the rate of cleavage of unmodified DNA but, in contrast, slightly increased the cleavage of ^N4mC-modified DNA.     

Illumina TruSeq Synthetic Long-Reads Empower De Novo Assembly and Resolve Complex,Highly-Repetitive Transposable Elements

High-throughput DNA sequencing technologies have revolutionized genomic analysis, including the de novo assembly of whole genomes. Nevertheless, assembly of complex genomes remains challenging, in part due to the presence of dispersed repeats which introduce ambiguity during genome reconstruction. Transposable elements (TEs) can be particularly problematic, especially for TE families exhibiting high sequence identity, high copy number, or complex genomic arrangements. While TEs strongly affect genome function and evolution, most current de novo assembly approaches cannot resolve long, identical, and abundant families of TEs. Here, we applied a novel Illumina technology called TruSeq synthetic long-reads, which are generated through highly-parallel library preparation and local assembly of short read data and which achieve lengths of 1.5–18.5 Kbp with an extremely low error rate (0.03% per base). To test the utility of this technology, we sequenced and assembled the genome of the model organism Drosophila melanogaster (reference genome strain y; cn, bw, sp) achieving an N50 contig size of 69.7 Kbp and covering 96.9% of the euchromatic chromosome arms of the current reference genome. TruSeq synthetic long-read technology enables placement of individual TE copies in their proper genomic locations as well as accurate reconstruction of TE sequences. We entirely recovered and accurately placed 4,229 (77.8%) of the 5,434 annotated transposable elements with perfect identity to the current reference genome. As TEs are ubiquitous features of genomes of many species, TruSeq synthetic long-reads, and likely other methods that generate long-reads, offer a powerful approach to improve de novo assemblies of whole genomes.     

Structural basis of the methylation specificity of R.DpnI

R.DpnI consists of N-terminal catalytic and C-terminal winged helix domains that are separately specific for the Gm6ATC sequences in Dam-methylated DNA. Here we present a crystal structure of R.DpnI with oligoduplexes bound to the catalytic and winged helix domains and identify the catalytic domain residues that are involved in interactions with the substrate methyl groups. We show that these methyl groups in the Gm6ATC target sequence are positioned very close to each other. We further show that the presence of the two methyl groups requires a deviation from B-DNA conformation to avoid steric conflict. The methylation compatible DNA conformation is complementary with binding sites of both R.DpnI domains. This indirect readout of methylation adds to the specificity mediated by direct favorable interactions with the methyl groups and solvation/desolvation effects. We also present hydrogen/deuterium exchange data that support ‘crosstalk’ between the two domains in the identification of methylated DNA, which should further enhance R.DpnI methylation specificity.     

Common patterns in type II restriction enzyme binding sites

Restriction enzymes are among the best studied examples of DNA binding proteins. In order to find general patterns in DNA recognition sites, which may reflect important properties of protein–DNA interaction, we analyse the binding sites of all known type II restriction endonucleases. We find a significantly enhanced GC content and discuss three explanations for this phenomenon. Moreover, we study patterns of nucleotide order in recognition sites. Our analysis reveals a striking accumulation of adjacent purines (R) or pyrimidines (Y). We discuss three possible reasons: RR/YY dinucleotides are characterized by (i) stronger H-bond donor and acceptor clusters, (ii) specific geometrical properties and (iii) a low stacking energy. These features make RR/YY steps particularly accessible for specific protein–DNA interactions. Finally, we show that the recognition sites of type II restriction enzymes are underrepresented in host genomes and in phage genomes.     

Shewanella spp. Genomic Evolution for a Cold Marine Lifestyle and In-Situ Explosive Biodegradation

Shewanella halifaxensis and Shewanella sediminis were among a few aquatic γ-proteobacteria that were psychrophiles and the first anaerobic bacteria that degraded hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX). Although many mesophilic or psychrophilic strains of Shewanella and γ-proteobacteria were sequenced for their genomes, the genomic evolution pathways for temperature adaptation were poorly understood. On the other hand, the genes responsible for anaerobic RDX mineralization pathways remain unknown. To determine the unique genomic properties of bacteria responsible for both cold-adaptation and RDX degradation, the genomes of S. halifaxensis and S. sediminis were sequenced and compared with 108 other γ-proteobacteria including Shewanella that differ in temperature and Na⁺ requirements, as well as RDX degradation capability. Results showed that for coping with marine environments their genomes had extensively exchanged with deep sea bacterial genomes. Many genes for Na⁺-dependent nutrient transporters were recruited to use the high Na⁺ content as an energy source. For coping with low temperatures, these two strains as well as other psychrophilic strains of Shewanella and γ-proteobacteria were found to decrease their genome G+C content and proteome alanine, proline and arginine content (p-value <0.01) to increase protein structural flexibility. Compared to poorer RDX-degrading strains, S. halifaxensis and S. sediminis have more number of genes for cytochromes and other enzymes related to RDX metabolic pathways. Experimentally, one cytochrome was found induced in S. halifaxensis by RDX when the chemical was the sole terminal electron acceptor. The isolated protein degraded RDX by mono-denitration and was identified as a multiheme 52 kDa cytochrome using a proteomic approach. The present analyses provided the first insight into divergent genomic evolution of bacterial strains for adaptation to the specific cold marine conditions and to the degradation of the pollutant RDX. The present study also provided the first evidence for the involvement of a specific c-type cytochrome in anaerobic RDX metabolism.