Manipulating replisome dynamics to enhance lambda Red-mediated multiplex genome engineering.

Disrupting the interaction between primase and helicase in Escherichia coli increases Okazaki fragment (OF) length due to less frequent primer synthesis. We exploited this feature to increase the amount of ssDNA at the lagging strand of the replication fork that is available for λ Red-mediated Multiplex Automatable Genome Engineering (MAGE). Supporting this concept, we demonstrate that MAGE enhancements correlate with OF length. Compared with a standard recombineering strain (EcNR2), the strain with the longest OFs displays on average 62% more alleles converted per clone, 239% more clones with 5 or more allele conversions and 38% fewer clones with 0 allele conversions in 1 cycle of co-selection MAGE (CoS-MAGE) with 10 synthetic oligonucleotides. Additionally, we demonstrate that both synthetic oligonucleotides and accessible ssDNA targets on the lagging strand of the replication fork are limiting factors for MAGE. Given this new insight, we generated a strain with reduced oligonucleotide degradation and increased genomic ssDNA availability, which displayed 111% more alleles converted per clone, 527% more clones with 5 or more allele conversions and 71% fewer clones with 0 allele conversions in 1 cycle of 10-plex CoS-MAGE. These improvements will facilitate ambitious genome engineering projects by minimizing dependence on time-consuming clonal isolation and screening.

Large scale bacterial gene discovery by similarity search.

DNA sequencing efforts frequently uncover genes other than the targeted ones. We have used rapid database scanning methods to search for undescribed eubacterial and archean protein coding frames in regions flanking known genes. By searching all prokaryotic DNA sequences not marked as coding for proteins or stable RNAs against the protein databases, we have identified more than 450 new examples of bacterial proteins, as well as a smaller number of possible revisions to known proteins, at a surprisingly high rate of one new protein or revision for every 24 initial DNA sequences or 8,300 nucleotides examined. Seven proteins are members of families which have not been described in prokaryotic sequences. We also describe 49 re-interpretations of existing sequence data of particular biological significance.

Identifying regulatory networks by combinatorial analysis of promoter elements.

Several computational methods based on microarray data are currently used to study genome-wide transcriptional regulation. Few studies, however, address the combinatorial nature of transcription, a well-established phenomenon in eukaryotes. Here we describe a new approach using microarray data to uncover novel functional motif combinations in the promoters of Saccharomyces cerevisiae. In addition to identifying novel motif combinations that affect expression patterns during the cell cycle, sporulation and various stress responses, we observed regulatory cross-talk among several of these processes. We have also generated motif-association maps that provide a global view of transcription networks. The maps are highly connected, suggesting that a small number of transcription factors are responsible for a complex set of expression patterns in diverse conditions. This approach may be useful for modeling transcriptional regulatory networks in more complex eukaryotes.

Correlation between transcriptome and interactome mapping data from Saccharomyces cerevisiae.

Genomic and proteomic approaches can provide hypotheses concerning function for the large number of genes predicted from genome sequences. Because of the artificial nature of the assays, however, the information from these high-throughput approaches should be considered with caution. Although it is possible that more meaningful hypotheses could be formulated by integrating the data from various functional genomic and proteomic projects, it has yet to be seen to what extent the data can be correlated and how such integration can be achieved. We developed a 'transcriptome-interactome correlation mapping' strategy to compare the interactions between proteins encoded by genes that belong to common expression-profiling clusters with those between proteins encoded by genes that belong to different clusters. Using this strategy with currently available data sets for Saccharomyces cerevisiae, we provide the first global evidence that genes with similar expression profiles are more likely to encode interacting proteins. We show how this correlation between transcriptome and interactome data can be used to improve the quality of hypotheses based on the information from both approaches. The strategy described here may help to integrate other functional genomic and proteomic data, both in yeast and in higher organisms.

A computational analysis of whole-genome expression data reveals chromosomal domains of gene expression.

Chromosome correlation maps display correlations between the expression patterns of genes on the same chromosome. Using these maps, we show here that adjacent pairs of genes, as well as nearby non-adjacent pairs of genes, show correlated expression independent of their orientation. We present specific examples of adjacent pairs with highly correlated expression patterns, in which the promoter of only one of the two genes contains an upstream activating sequence (UAS) known to be associated with that expression pattern. Finally, we show that genes with similar functions tend to occur in adjacent positions along the chromosomes. Our results suggest that, in certain chromosomal expression domains, an UAS can affect the transcription of genes that are not immediately downstream from it.

Systematic determination of genetic network architecture.

Technologies to measure whole-genome mRNA abundances and methods to organize and display such data are emerging as valuable tools for systems-level exploration of transcriptional regulatory networks. For instance, it has been shown that mRNA data from 118 genes, measured at several time points in the developing hindbrain of mice, can be hierarchically clustered into various patterns (or 'waves') whose members tend to participate in common processes. We have previously shown that hierarchical clustering can group together genes whose cis-regulatory elements are bound by the same proteins in vivo. Hierarchical clustering has also been used to organize genes into hierarchical dendograms on the basis of their expression across multiple growth conditions. The application of Fourier analysis to synchronized yeast mRNA expression data has identified cell-cycle periodic genes, many of which have expected cis-regulatory elements. Here we apply a systematic set of statistical algorithms, based on whole-genome mRNA data, partitional clustering and motif discovery, to identify transcriptional regulatory sub-networks in yeast-without any a priori knowledge of their structure or any assumptions about their dynamics. This approach uncovered new regulons (sets of co-regulated genes) and their putative cis-regulatory elements. We used statistical characterization of known regulons and motifs to derive criteria by which we infer the biological significance of newly discovered regulons and motifs. Our approach holds promise for the rapid elucidation of genetic network architecture in sequenced organisms in which little biology is known.

Multiplex DNA sequencing.

The increasing demand for DNA sequences can be met by replacement of each DNA sample in a device with a mixture of N samples so that the normal throughput is increased by a factor of N. Such a method is described. In order to separate the sequence information at the end of the processing, the DNA molecules of interest are ligated to a set of oligonucleotide "tags" at the beginning. The tagged DNA molecules are pooled, amplified, and chemically fragmented in 96-well plates. The resulting reaction products are fractionated by size on sequencing gels and transferred to nylon membranes. These membranes are then probed as many times as there are types of tags in the original pools, producing, in each cycle of probing, autoradiographs similar to those from standard DNA sequencing methods. Thus, each reaction and gel yields a quantity of data equivalent to that obtained from conventional reactions and gels multiplied by the number of probes used. To date, even after 50 successive probings, the original signal strength and the image quality are retained, an indication that the upper limit for the number of reprobings may be considerably higher.

B lineage--specific interactions of an immunoglobulin enhancer with cellular factors in vivo.

The mouse heavy chain immunoglobulin gene contains a tissue-specific enhancer. The enhancer and flanking sequences were studied in vivo by carrying out dimethyl sulfate protection experiments on living cells, in combination with genomic sequencing. Relative to reactions on naked DNA, there are changes (protections and enhancements) in the reactivity of guanine residues to dimethyl sulfate within the enhancer sequence in myeloma, B, and early B cells, whereas virtually no alterations appear in cells of non-B lineage. Most of the affected residues are in four clusters, in sequences homologous to the octamer 5'CAGGTGGC 3' (C, cytosine; A, adenine; G. guanine; T, thymine). The alterations in the pattern of G reactivity are consistent with the tissue-specific binding of molecules to the mouse immunoglobulin heavy chain enhancer.

Microarray analysis of the transcriptional network controlled by the photoreceptor homeobox gene Crx.

BACKGROUND: Terminal differentiation of many cell types is controlled and maintained by tissue- or cell-specific transcription factors. Little is known, however, of the transcriptional networks controlled by such factors and how they regulate differentiation. The paired-type homeobox transcription factor, Crx, has a pivotal role in the terminal differentiation of vertebrate photoreceptors. Mutations in the human CRX gene result in either congenital blindness or photoreceptor degeneration and targeted mutation of the mouse Crx results in failure of development of the light-detecting outer segment of photoreceptors. RESULTS: We have characterized the transcriptional network controlled by Crx by microarray analysis of gene expression in developing retinal tissue from Crx(+/+) and Crx(-/-) mice. These data were combined with analyses of gene expression in developing and adult retina, as well as adult brain. The most abundant elements of this network are ten photoreceptor-specific or -enriched genes, including six phototransduction genes. All of the available 5' regulatory regions of the putative Crx targets contain a novel motif that is composed of a head-to-tail arrangement of two Crx-binding-element-like sequences. Analysis of the 5' regions of a set of mouse and human genes suggests that this motif is specific to Crx targets. CONCLUSIONS: This study demonstrates that cDNA microarrays can be successfully used to define the transcriptional networks controlled by transcription factors in vertebrate tissue in vivo.

Regulatory networks revealed by transcriptional profiling of damaged Saccharomyces cerevisiae cells: Rpn4 links base excision repair with proteasomes.

Exposure to carcinogenic alkylating agents, oxidizing agents, and ionizing radiation modulates transcript levels for over one third of Saccharomyces cerevisiae's 6,200 genes. Computational analysis delineates groups of coregulated genes whose upstream regions bear known and novel regulatory sequence motifs. One group of coregulated genes contain a number of DNA excision repair genes (including the MAG1 3-methyladenine DNA glycosylase gene) and a large selection of protein degradation genes. Moreover, transcription of these genes is modulated by the proteasome-associated protein Rpn4, most likely via its binding to MAG1 upstream repressor sequence 2-like elements, that turn out to be almost identical to the recently identified proteasome-associated control element (G. Mannhaupt, R. Schnall, V. Karpov, I. Vetter, and H. Feldmann, FEBS Lett. 450:27-34, 1999). We have identified a large number of genes whose transcription is influenced by Rpn4p.

Quantitative whole-genome analysis of DNA-protein interactions by in vivo methylase protection in E. coli.

A global methylation-based technique was used to identify, display, and quantitate the in vivo occupancy of numerous protein-binding sites within the Escherichia coli genome. The protein occupancy profiles of these sites showed variation across different growth conditions and genetic backgrounds. Of the 25 sites identified in this study, 24 occurred within 5' noncoding regions. Protein occupancy at 13 of these sites was supported by independent biochemical and genetic evidence. Most of the remaining 12 sites fell upstream of genes with no previously known function. A multivariate statistical analysis was utilized to group such uncharacterized genes with well-characterized ones, providing insights into their function based on a common pattern of transcriptional regulation.

Quantifying DNA-protein interactions by double-stranded DNA arrays.

We have created double-stranded oligonucleotide arrays to perform highly parallel investigations of DNA-protein interactions. Arrays of single-stranded DNA oligonucleotides, synthesized by a combination of photolithography and solid-state chemistry, have been used for a variety of applications, including large-scale mRNA expression monitoring, genotyping, and sequence-variation analysis. We converted a single-stranded to a double-stranded array by synthesizing a constant sequence at every position on an array and then annealing and enzymatically extending a complementary primer. The efficiency of second-strand synthesis was demonstrated by incorporation of fluorescently labeled dNTPs (2'-deoxyribonucleoside 5'-triphosphates) and by terminal transferase addition of a fluorescently labeled ddNTP. The accuracy of second-strand synthesis was demonstrated by digestion of the arrayed double-stranded DNA (dsDNA) on the array with sequence-specific restriction enzymes. We showed dam methylation of dsDNA arrays by digestion with DpnI, which cleaves when its recognition site is methylated. This digestion demonstrated that the dsDNA arrays can be further biochemically modified and that the DNA is accessible for interaction with DNA-binding proteins. This dsDNA array approach could be extended to explore the spectrum of sequence-specific protein binding sites in genomes.

Finding DNA regulatory motifs within unaligned noncoding sequences clustered by whole-genome mRNA quantitation.

Whole-genome mRNA quantitation can be used to identify the genes that are most responsive to environmental or genotypic change. By searching for mutually similar DNA elements among the upstream non-coding DNA sequences of these genes, we can identify candidate regulatory motifs and corresponding candidate sets of coregulated genes. We have tested this strategy by applying it to three extensively studied regulatory systems in the yeast Saccharomyces cerevisiae: galactose response, heat shock, and mating type. Galactose-response data yielded the known binding site of Gal4, and six of nine genes known to be induced by galactose. Heat shock data yielded the cell-cycle activation motif, which is known to mediate cell-cycle dependent activation, and a set of genes coding for all four nucleosomal proteins. Mating type alpha and a data yielded all of the four relevant DNA motifs and most of the known a- and alpha-specific genes.

Selection analyses of insertional mutants using subgenic-resolution arrays.

We describe a method of genome-wide analysis of quantitative growth phenotypes using insertional mutagenesis and DNA microarrays. We applied the method to assess the fitness contributions of Escherichia coli gene domains under specific growth conditions. A transposon library was subjected to competitive growth selection in Luria-Bertani (LB) and in glucose minimal media. Transposon-containing genomic DNA fragments from the selected libraries were compared with the initial unselected transposon insertion library on DNA microarrays to identify insertions that affect fitness. Genes involved in the biosynthesis of nutrients not provided in the growth medium were found to be significantly enriched in the set of genes containing negatively selected insertions. The data also identify fitness contributions of several uncharacterized genes, including putative transcriptional regulators and enzymes. The applicability of this high-resolution array selection in other species is discussed.

RNA expression analysis using a 30 base pair resolution Escherichia coli genome array.

We have developed a high-resolution "genome array" for the study of gene expression and regulation in Escherichia coli. This array contains on average one 25-mer oligonucleotide probe per 30 base pairs over the entire genome, with one every 6 bases for the intergenic regions and every 60 bases for the 4,290 open reading frames (ORFs). Twofold concentration differences can be detected at levels as low as 0.2 messenger RNA (mRNA) copies per cell, and differences can be seen over a dynamic range of three orders of magnitude. In rich medium we detected transcripts for 97% and 87% of the ORFs in stationary and log phases, respectively. We found that 1, 529 transcripts were differentially expressed under these conditions. As expected, genes involved in translation were expressed at higher levels in log phase, whereas many genes known to be involved in the starvation response were expressed at higher levels in stationary phase. Many previously unrecognized growth phase-regulated genes were identified, such as a putative receptor (b0836) and a 30S ribosomal protein subunit (S22), both of which are highly upregulated in stationary phase. Transcription of between 3,000 and 4,000 predicted ORFs was observed from the antisense strand, indicating that most of the genome is transcribed at a detectable level. Examples are also presented for high-resolution array analysis of transcript start and stop sites and RNA secondary structure.

Automated electrotransformation of Escherichia coli on a digital microfluidic platform using bioactivated magnetic beads.

This paper reports on the use of a digital microfluidic platform to perform multiplex automated genetic engineering (MAGE) cycles on droplets containing Escherichia coli cells. Bioactivated magnetic beads were employed for cell binding, washing, and media exchange in the preparation of electrocompetent cells in the electrowetting-on-dieletric (EWoD) platform. On-cartridge electroporation was used to deliver oligonucleotides into the cells. In addition to the optimization of a magnetic bead-based benchtop protocol for generating and transforming electrocompetent E. coli cells, we report on the implementation of this protocol in a fully automated digital microfluidic platform. Bead-based media exchange and electroporation pulse conditions were optimized on benchtop for transformation frequency to provide initial parameters for microfluidic device trials. Benchtop experiments comparing electrotransformation of free and bead-bound cells are presented. Our results suggest that dielectric shielding intrinsic to bead-bound cells significantly reduces electroporation field exposure efficiency. However, high transformation frequency can be maintained in the presence of magnetic beads through the application of more intense electroporation pulses. As a proof of concept, MAGE cycles were successfully performed on a commercial EWoD cartridge using variations of the optimal magnetic bead-based preparation procedure and pulse conditions determined by the benchtop results. Transformation frequencies up to 22% were achieved on benchtop; this frequency was matched within 1% (21%) by MAGE cycles on the microfluidic device. However, typical frequencies on the device remain lower, averaging 9% with a standard deviation of 9%. The presented results demonstrate the potential of digital microfluidics to perform complex and automated genetic engineering protocols.

Predicting regulons and their cis-regulatory motifs by comparative genomics.

We have combined and compared three techniques for predicting functional interactions based on comparative genomics (methods based on conserved operons, protein fusions and correlated evolution) and optimized these methods to predict coregulated sets of genes in 24 complete genomes, including Saccharomyces cerevisiae, Caenorhabditis elegans and 22 prokaryotes. The method based on conserved operons was the most useful for this purpose. Upstream regions of the genes comprising these predicted regulons were then used to search for regulatory motifs in 22 prokaryotic genomes using the motif-discovery program AlignACE. Many significant upstream motifs, including five known Escherichia coli regulatory motifs, were identified in this manner. The presence of a significant regulatory motif was used to refine the members of the predicted regulons to generate a final set of predicted regulons that share significant regulatory elements.

In situ localized amplification and contact replication of many individual DNA molecules.

We describe a method to clone and amplify DNA by performing the polymerase chain reaction (PCR) in a thin polyacrylamide film poured on a glass microscope slide. The polyacrylamide matrix retards the diffusion of the linear DNA molecules so that the amplification products remain localized near their respective templates. At the end of the reaction, a number of PCR colonies, or 'polonies', have formed, each one grown from a single template molecule. As many as 5 million clones can be amplified in parallel on a single slide. If an Acrydite modification is included at the 5' end of one of the primers, the amplified DNA will be covalently attached to the polyacrylamide matrix, allowing further enzymatic manipulations to be performed on all clones simultaneously. We describe techniques to make replicas of these polony slides, and high throughput sequencing protocols for this technology. Other applications are also discussed.

Overcoming Challenges in Engineering the Genetic Code.

Withstanding 3.5 billion years of genetic drift, the canonical genetic code remains such a fundamental foundation for the complexity of life that it is highly conserved across all three phylogenetic domains. Genome engineering technologies are now making it possible to rationally change the genetic code, offering resistance to viruses, genetic isolation from horizontal gene transfer, and prevention of environmental escape by genetically modified organisms. We discuss the biochemical, genetic, and technological challenges that must be overcome in order to engineer the genetic code.

Complementary recognition in condensed DNA: accelerated DNA renaturation.

The functional consequences of DNA condensation are investigated. The recognition of complementary strands is profoundly modified by this critical phenomenon. (1) Condensation of denatured DNA greatly accelerates the kinetics of DNA renaturation. We propose a unifying explanation for the effects of several accelerating solvents studied here including polymers, di- and multivalent cations, as well as effects seen with the phenol emulsions and single-stranded nucleic acid binding proteins. Optimal conditions for renaturation at or above the calculated three dimensional diffusion limit are theoretically consistent with a limited search space in the condensed phases. (2) In addition to these effects on association of two single strands, similar condensation acceleration effects can be seen in strand exchange experiments with double stranded DNA without proteins. These may model a mechanism of recombinational protein function.